本发明涉及医学图像处理方法，尤其涉及一种基于多方向轮廓的医学图像分割方法。
背景技术：
现有技术中，在临床外科医生的合作过程中，了解到临床医生在疾病诊断、治疗决策以及手术规划过程中，需要掌握组织以及病灶的解剖结构信息，例如肝脏肿瘤的位置、体积、形状等。目前多是利用手工标注的方式，结合阈值连通、区域增长等图像处理方法，对ct或者mri图像进行逐层分割病灶区域。随着技术发展进步，图像分辨率越来越高，层间距越来越小，这就给手工分割带来了挑战。1)逐层手工标注是一个非常耗时耗力的工作，分割一副3d图像通常耗费数小时甚至更多的时间；而且，分割结果受标注人的主观影响较大，不同标注人对同一个对象的分割结果、甚至同一标注人对同一对象的多次分割都会存在较大的差异；2)逐层标注方式未考虑三维信息，使得分割体数据的边界非常粗糙，影响建模精度和可视化效果；3)对于肿瘤等病灶组织，经常存在模糊边缘，这使得标注该区域时，难以确定边界；而肿瘤组织内部灰度分布具有异质性，使得阈值连通、区域增长等自动分割方法都不能完整分割该区域，需要借助大量手工后处理，增加工作量。
虽然基于轮廓约束的智能重建方法能够较好的解决以上问题，但是仍然存在如下缺陷：1)操作较为繁琐。为了获取精确的重建结果，需要分别从三维图像的横断面、冠状面、矢状面、甚至是斜切面等多个方向切面上进行手工分割以提取轮廓；2)从多个方向切面上提取轮廓时，未考虑轮廓之间的空间对应，使得重建的三维模型存在凹陷，影响重建精度。
技术实现要素：
本发明要解决的技术问题在于，针对现有技术的不足，提供一种基于多方向轮廓的医学图像分割方法，该方法利用一次操作即时提取待分割对象在多个不同方向切面上的轮廓，从而建立形状约束，并通过插值重建技术进行曲面重建，有效提高建模效率和精度。
为解决上述技术问题，本发明采用如下技术方案。
一种基于多方向轮廓的医学图像分割方法，包括：步骤s1，数据准备，获取序列薄层的ct或mri医学图像；步骤s2，导入序列医学图像，对于图像进行mpr重建，并在视窗中显示医学图像，可以在显示视图中显示图像横断面、冠状面和矢状面，并且可以根据需要来添加和调节任意方向切面进行图像显示；步骤s3，对医学图像中的分割对象选取多角度切面，利用鼠标在任意切面上待分割对象内部设置一个种子点，在多角度切面上使用魔术棒分割方法分割阈值近似区域，并依据区域边缘提取各个切面上的2d轮廓，得到多方向轮廓；步骤s4，利用生成的多方向轮廓作为形状约束，基于插值重建算法实现目标区域的表面网格重建，获得网格重模型，然后根据空间位置，将这个模型映射至图像空间，得到二值分割mask图像；步骤s5，在多角度图像切面视窗中，逐层判断二值分割mask图像是否准确，若是，则执行步骤s7，若否，则执行步骤s6；步骤s6，编辑切面轮廓，选择二值分割mask图像中分割不精确的切面，根据二值分割mask图像中的边缘信息生成轮廓线和关键点，利用鼠标对轮廓线上的关键点进行局部调节，调节结束后，重新执行步骤s4；步骤s7，输出目标区域的二值分割mask图像和表面网格，分割算法结束。
优选地，所述步骤s1、步骤s2中，对输入ct、mri等医学影像序列图像进行mpr重建后，支持调节任意方向进行多角度切面图像显示浏览。
优选地，所述步骤s3中，基于1个种子点，通过魔术棒工具在多角度切面上分割提取2d轮廓集合。
优选地，所述步骤s4中，以多个不同方向的2d轮廓作为形状约束，通过插值重建算法实现目标区域的表面网格重建，并实现网格模型映射为二值分割mask图像。
优选地，所述步骤s6中，利用轮廓编辑工具，从二值分割mask图像的切面分割图像中提取轮廓线，并支持基于关键点的轮廓线局部调节。
本发明公开的基于多方向轮廓的医学图像分割方法中，通过一次操作，同时产生多个方向上待分割对象的轮廓，并以这些轮廓建立形状约束，通过插值重建算法实现精确建模，对于重建结果，提供工具进行迭代式编辑和重建，进行局部优化，提高分割精度。该方法对于图像灰度依赖性小，可以适用于不同类型图像的分割建模。相比现有技术而言，本发明利用一次操作即时提取待分割对象在多个不同方向切面上的轮廓，从而建立形状约束，并通过插值重建技术进行曲面重建，有效提高了建模效率和精度。
附图说明
图1为本发明医学图像分割方法的流程图。
图2为优选实施例中横断面子宫包后块区域增长分割图像。
图3为矢状面上子宫包后块区域增长分割图像。
图4为冠状面上子宫包后块区域增长分割图像。
图5为三个切面的轮廓线。
图6为初始子宫包后块的表面网络模型。
图7为在冠状面新添加一个轮廓线的图示。
图8为优化后子宫包后块的表面网络模型。
具体实施方式
下面结合附图和实施例对本发明作更加详细的描述。
本发明公开了一种基于多方向轮廓的医学图像分割方法，请参照图1，其包括有如下步骤：
步骤s1，数据准备，获取序列薄层的ct或mri医学图像；
步骤s2，导入序列医学图像，对于图像进行mpr重建，并在视窗中显示医学图像，可以在显示视图中显示图像横断面、冠状面和矢状面，并且可以根据需要来添加和调节任意方向切面进行图像显示；
步骤s3，对医学图像中的分割对象选取多角度切面，利用鼠标在任意切面上待分割对象内部设置一个种子点，在多角度切面上使用魔术棒分割方法分割阈值近似区域，并依据区域边缘提取各个切面上的2d轮廓，得到多方向轮廓；
步骤s4，利用生成的多方向轮廓作为形状约束，基于插值重建算法实现目标区域的表面网格重建，获得网格重模型，然后根据空间位置，将这个模型映射至图像空间，得到二值分割mask图像；
步骤s5，在多角度图像切面视窗中，逐层判断二值分割mask图像是否准确，若是，则执行步骤s7，若否，则执行步骤s6；
步骤s6，编辑切面轮廓，选择二值分割mask图像中分割不精确的切面，根据二值分割mask图像中的边缘信息生成轮廓线和关键点，利用鼠标对轮廓线上的关键点进行局部调节，调节结束后，重新执行步骤s4；
步骤s7，输出目标区域的二值分割mask图像和表面网格，分割算法结束。
上述医学图像分割方法中，通过一次操作，同时产生多个方向上待分割对象的轮廓，并以这些轮廓建立形状约束，通过插值重建算法实现精确建模，对于重建结果，提供工具进行迭代式编辑和重建，进行局部优化，提高分割精度。该方法对于图像灰度依赖性小，可以适用于不同类型图像的分割建模。相比现有技术而言，本发明利用一次操作即时提取待分割对象在多个不同方向切面上的轮廓，从而建立形状约束，并通过插值重建技术进行曲面重建，有效提高了建模效率和精度。
本实施例中，所述步骤s1中，所述医学图像是由医学影像检测设备上产生的个体化序列的图像。进一步地，所述医学影像检测设备为ct、mri或pet设备。
本实施例的所述步骤s1、步骤s2中，对输入ct、mri等医学影像序列图像进行mpr重建后，支持调节任意方向进行多角度切面图像显示浏览。
作为一种优选方式，所述步骤s3中，基于1个种子点，通过魔术棒工具在多角度切面上分割提取2d轮廓集合。
进一步地，所述步骤s4中，以多个不同方向的2d轮廓作为形状约束，通过插值重建算法实现目标区域的表面网格重建，并实现网格模型映射为二值分割mask图像。
本实施例的所述步骤s6中，利用轮廓编辑工具，从二值分割mask图像的切面分割图像中提取轮廓线，并支持基于关键点的轮廓线局部调节。
本发明基于用户设置的种子点，同时在横断面、冠状面和矢状面方向上进行魔术棒分割，并提取相应的轮廓作为插值重建算法输入。其次，利用提取的多方向轮廓信息作为形状约束，通过插值重建算法实现目标区域的3d表面网格重建，并生成二值分割mask图像。此外，利用轮廓编辑工具对插值重建分割结果中不精确或者有误的层面进行局部调整后再次进行插值重建。并通过迭代式策略进行不断优化，直至得到满意的分割结果。
本发明公开的基于多方向轮廓的医学图像分割方法，其实际应用过程中，可参考如下步骤：
1、数据准备的过程主要的工作是，获取序列薄层的ct或mri医学图像；
2、对导入图像进行mpr重建，并在视窗中显示。可以在显示视图中显示图像横断面、冠状面和矢状面，并且可以根据需要来添加和调节任意方向切面进行图像显示；
3、调节图像使得分割对象在选取的多角度切面上清晰，利用鼠标在任意切面上待分割对象内部设置一个种子点，并在多角度切面上进行魔术棒分割，并提取2d轮廓集合；
4、利用生成的多方向轮廓集合作为形状约束，通过插值重建算法实现目标区域的表面网格重建，并映射至图像空间，得到二值分割mask图像；
5、逐层观察二值分割mask是否准确。如果无需编辑2d切面轮廓，则输出目标区域的二值分割和表面网格，分割算法结束，如果分割存在误差或者错误，则需对二值分割切面进行调整；
6、选择二值分割mask中分割不精确的切面，利用轮廓编辑工具提取二值分割图像切面的轮廓线和关键点，利用鼠标对轮廓线关键点并进行局部调节。调节结束后，重新执行插值重建，更新网格模型和二值分割mask。
本发明公开的基于多方向轮廓的医学图像分割方法，以子宫包后块为例，其处理过程包括：
首先，调整ct图像，使得子宫包后块在横断面，矢状面和冠状面方向上可见。并选择横断面，在子宫包后块内部设置种子点，并在三个方向切面图像中进行魔术棒分割，结果如图2、图3、图4。
请参照图5，分别从三个方向切面图像的魔术棒分割结果中提取三个轮廓线，作为插值重建算法的输入。
请参照图6，基于当前三个方向上的轮廓线，通过插值重建算法，生成子宫包后块的表面网格模型以及相应的二值分割图像。
请参照图7，观察当前分割，对于分割不精确的层面，利用轮廓编辑工具对分割层面的轮廓进行局部微调，更新形状轮廓约束。
请参照图8，编辑完成后，重新进行插值重建，即可得到优化后的子宫包后块的表面网格模型以及相应的二值分割图像。
反复执行局部微调、形状约束和插值重建，直至生成子宫包后块的表面网络模型以及相应的二值分割图像满足要求后停止。
本发明公开的基于多方向轮廓的医学图像分割方法，其相比现有技术而言的有益效果在于，本发明操作简单、高效，通过一次设置种子点，自动得到多个方向上的轮廓，并快速重建得到表面网格和二值分割，同时，本发明重建分割精度高，重建算法的输入轮廓自动生成，保持轮廓之间的空间对应，避免轮廓二义性导致的重建错误，能够生成更为平滑的网格。
以上所述只是本发明较佳的实施例，并不用于限制本发明，凡在本发明的技术范围内所做的修改、等同替换或者改进等，均应包含在本发明所保护的范围内。