一种CT图像交互分割方法及装置

摘要

本发明提供了一种CT图像交互分割方法及装置，其中方法包括：调整CT强度阈值范围，获取CT强度值属于强度阈值范围内的待显示部位，并对待显示部位进行标记着色，利用体绘制技术模拟等值面显示，显示原始数据体绘制视图；通过标记工具确定原始数据体绘制视图中的待分割区域，使用填充工具对待分割区域进行计算，生成初步分割结果，并使用腐蚀膨胀工具调整初步分割结果；在初步分割结果的基础上，使用标记工具选中待分割区域中未分割出的区域，生成分割结果，显示分割体绘制视图。由此可以以直观的交互手段简化分割过程，节省大量的分割时间，更利于用户使用。

权利要求

1.一种CT图像交互分割方法，其特征在于，包括：

调整CT强度阈值范围，获取所述CT强度值属于所述强度阈值范围内的待显示部位，并对所述待显示部位进行标记着色，利用体绘制技术模拟等值面显示，显示原始数据体绘制视图；

通过标记工具确定所述原始数据体绘制视图中的待分割区域，使用填充工具对所述待分割区域进行计算，生成初步分割结果，并使用腐蚀膨胀工具调整所述初步分割结果；

在所述初步分割结果的基础上，使用标记工具选中所述待分割区域中未分割出的区域，生成分割结果，显示分割体绘制视图；

其中：所述通过标记工具确定所述原始数据体绘制视图中的待分割区域包括：

在所述原始数据体绘制视图中，通过拖动鼠标，标记拖动轨迹附近的区域作为所述待分割区域，所述拖动轨迹终止于深度突变前的位置；其中，所述拖动轨迹为:在拖动的某一时刻，视线通过所述鼠标指针的延长线与当前所显示的等值面产生若干交点，取其中最先相交的两点的中点，将所有时刻的中点连接成的折线。

2.根据权利要求1所述的方法，其特征在于，

所述显示原始数据体绘制视图的同时，还包括：

显示所述原始数据体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图，且在所述原始数据体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或所述显示分割体绘制视图的同时，还包括：

显示所述分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图，且在所述分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示。

3.根据权利要求1所述的方法，其特征在于，所述使用填充工具对所述待分割区域进行计算，生成初步分割结果包括：利用鼠标在所述原始数据体绘制视图中设置种子点，选中与所述种子点相关联的区域，并利用洪水填充算法或者区域生长算法，计算所述初步分割结果。

4.根据权利要求3所述的方法，其特征在于，所述使用腐蚀膨胀工具调整所述初步分割结果包括：调整与所述种子点相关联的区域的边缘，对与所述种子点相关联的区域执行三维的形态学腐蚀/膨胀算法，调整所述初步分割结果。

5.一种CT图像交互分割装置，其特征在于，包括：

显示模块，用于调整CT强度阈值范围，获取所述CT强度值属于所述强度阈值范围内的待显示部位，并对所述待显示部位进行标记着色，利用体绘制技术模拟等值面显示，显示原始数据体绘制视图；

初步分割模块，用于通过标记工具确定所述原始数据体绘制视图中的待分割区域，使用填充工具对所述待分割区域进行计算，生成初步分割结果，并使用腐蚀膨胀工具调整所述初步分割结果；

分割模块，用于在所述初步分割结果的基础上，使用标记工具选中所述待分割区域中未分割出的区域，生成分割结果，通知所述显示模块显示分割体绘制视图；

其中：所述初步分割模块具体用于通过如下方式通过标记工具确定所述原始数据体绘制视图中的待分割区域：初步分割模块，还用于在所述原始数据体绘制视图中，通过拖动鼠标，标记拖动轨迹附近的区域作为所述待分割区域，所述拖动轨迹终止于深度突变前的位置；其中，所述拖动轨迹为:在拖动的某一时刻，视线通过所述鼠标指针的延长线与当前所显示的等值面产生若干交点，取其中最先相交的两点的中点，将所有时刻的中点连接成的折线。

6.根据权利要求5所述的装置，其特征在于，

显示模块，还用于在显示原始数据体绘制视图的同时，显示所述原始数据体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图，且在所述原始数据体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或显示分割体绘制视图的同时，显示所述分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图，且在所述分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示。

7.根据权利要求5所述的装置，其特征在于，所述初步分割模块具体用于通过如下方式使用填充工具对所述待分割区域进行计算，生成初步分割结果：所述初步分割模块，还用于利用鼠标在所述原始数据体绘制视图中设置种子点，选中与所述种子点相关联的区域，并利用洪水填充算法或者区域生长算法，计算所述初步分割结果。

8.根据权利要求7所述的装置，其特征在于，所述初步分割模块具体用于通过如下方式使用腐蚀膨胀工具调整所述初步分割结果：所述初步分割模块，还用于调整与所述种子点相关联的区域的边缘，对与所述种子点相关联的区域执行三维的形态学腐蚀/膨胀算法，调整所述初步分割结果。

说明书

一种CT图像交互分割方法及装置

技术领域

[0001] 本发明涉及数字化医疗领域，尤其涉及一种CT(Computed Tomography，计算机断层扫 描)图像交互分割方法及装置。

背景技术

[0002] 在CT图像分割领域，目前国际上比较流行的软件有美国犹他大学的Seg3D和德国癌症 研究中心的MITK工作台等。这些软件提供滤波、阈值等工具，可生成对目标对象的初步分 割。接下来需要人工地对每张切片涂抹、擦除标记，并查看其它视图确认分割效果，从而得 到最终的分割。实际分割数据时，仍需大量的人工交互操作，例如需要使用画笔工具逐一涂 抹图像切片等，而一个CT图像的切片层数常多达数百层，使用此方式完成数据分割非常耗 时，完成整个CT中的目标对象分割通常需要1～2小时，且难以分割血管等细小结构。而近 年的研究工作则根据分割对象的不同提出了若干种自动分割方法，但这些方法适用性较窄， 精度也较低，难以直接应用于实际操作中。

发明内容

[0003] 本发明旨在提供一种克服上述问题或者至少部分地解决上述问题的一种CT图像交互分 割方法及装置，能够通过直观、简单的交互得到准确的分割结果。

[0004] 为达到上述目的，本发明的技术方案具体是这样实现的：

[0005] 本发明的一个方面提供了一种CT图像交互分割方法，包括：调整CT强度阈值范围，获 取CT强度值属于强度阈值范围内的待显示部位，并对待显示部位进行标记着色，利用体绘 制技术模拟等值面显示，显示原始数据体绘制视图；通过标记工具确定原始数据体绘制视图 中的待分割区域，使用填充工具对待分割区域进行计算，生成初步分割结果，并使用腐蚀膨 胀工具调整初步分割结果；在初步分割结果的基础上，使用标记工具选中待分割区域中未分 割出的区域，生成分割结果，显示分割体绘制视图。

[0006] 另外，显示原始数据体绘制视图的同时，还包括：显示原始数据体绘制视图的矢状面切 片、冠状面切片和横截面切片的二维视图，且在原始数据体绘制视图的矢状面切片、冠状面 切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或显示分割体绘 制视图的同时，还包括：显示分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二 维视图，且在分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当 前全局阈值的区域进行着色显示。

[0007] 另外，通过标记工具确定原始数据体绘制视图中的待分割区域包括：在原始数据体绘制 视图中，通过拖动鼠标，标记拖动轨迹附近的区域作为待分割区域，拖动轨迹终止于深度突 变前的位置；其中，拖动轨迹为：在拖动的某一时刻，视线通过鼠标指针的延长线与当前所 显示的等值面产生若干交点，取其中最先相交的两点的中点，将所有时刻的中点连接成的折 线。

[0008] 另外，使用填充工具对待分割区域进行计算，生成初步分割结果包括：利用鼠标在原始 数据体绘制视图中设置种子点，选中与种子点相关联的区域，并利用洪水填充算法或者区域 生长算法，计算初步分割结果。

[0009] 另外，使用腐蚀膨胀工具调整初步分割结果包括：调整与种子点相关联的区域的边缘， 对与种子点相关联的区域执行三维的形态学腐蚀/膨胀算法，调整初步分割结果。

[0010] 另外，传递函数包括：

[0011]

[0012]

[0013] 其中，I为CT强度值，T(A)、T(B)表示分别映射到两个颜色标签的阈值；且CT强度值 满足强度阈值范围的部分被设置为不透明，不满足强度阈值范围的部分被设置为透明。

[0014] 本发明另一方面提供了一种CT图像交互分割装置，包括：显示模块，用于调整CT强度 阈值范围，获取CT强度值属于强度阈值范围内的待显示部位，并对待显示部位进行标记着 色，利用体绘制技术模拟等值面显示，显示原始数据体绘制视图；初步分割模块，用于通过 标记工具确定原始数据体绘制视图中的待分割区域，使用填充工具对待分割区域进行计算， 生成初步分割结果，并使用腐蚀膨胀工具调整初步分割结果；分割模块，用于在初步分割结 果的基础上，使用标记工具选中待分割区域中未分割出的区域，生成分割结果，通知显示模 块显示分割体绘制视图。

[0015] 另外，显示模块，还用于在显示原始数据体绘制视图的同时，显示原始数据体绘制视图 的矢状面切片、冠状面切片和横截面切片的二维视图，且在原始数据体绘制视图的矢状面切 片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或显 示分割体绘制视图的同时，显示分割体绘制视图的矢状面切片、冠状面切片和横截面切片的 二维视图，且在分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于 当前全局阈值的区域进行着色显示。

[0016] 另外，初步分割模块具体用于通过如下方式通过标记工具确定原始数据体绘制视图中的 待分割区域：初步分割模块，还用于在原始数据体绘制视图中，通过拖动鼠标，标记拖动轨 迹附近的区域作为待分割区域，所述拖动轨迹终止于深度突变前的位置；其中，拖动轨迹为： 在拖动的某一时刻，视线通过鼠标指针的延长线与当前所显示的等值面产生若干交点，取其 中最先相交的两点的中点，将所有时刻的中点连接成的折线。

[0017] 另外，初步分割模块具体用于通过如下方式使用填充工具对待分割区域进行计算，生成 初步分割结果：初步分割模块，还用于利用鼠标在原始数据体绘制视图中设置种子点，选中 与种子点相关联的区域，并利用洪水填充算法或者区域生长算法，计算初步分割结果。

[0018] 另外，初步分割模块具体用于通过如下方式使用腐蚀膨胀工具调整初步分割结果：初步 分割模块，还用于调整与种子点相关联的区域的边缘，对与种子点相关联的区域执行三维的 形态学腐蚀/膨胀算法，调整初步分割结果。

[0019] 另外，传递函数包括：

[0020]

[0021]

[0022] 其中，I为CT强度值，T(A)、T(B)表示分别映射到两个颜色标签的阈值；且CT强度值 满足强度阈值范围的部分被设置为不透明，不满足强度阈值范围的部分被设置为透明。

[0023] 由此可见，通过本发明提供的CT图像交互分割方法及装置，可以替代目前Seg3D、MITK 工作台等软件中针对CT图像的分割方法，以直观的交互手段简化分割过程，节省大量的分 割时间，更利于用户使用。

附图说明

[0024] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图 作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的 普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。

[0025] 图1为本发明实施例提供的CT图像交互分割方法的流程图；

[0026] 图2为本发明实施例提供的CT图像交互分割装置的结构示意图；

[0027] 图3为本发明实施例提供的显示多种颜色的传递函数示意图。

具体实施方式

[0028] 下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例 性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。 相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传 达给本领域的技术人员。

[0029] 针对CT图像数据的分割问题，本发明提出了一种CT图像交互分割方法，包括了原始数 据与分割结果的直观显示。

[0030] 图1示出了本发明实施例提供的CT图像交互分割方法的流程图，参见图1，CT图像交 互分割方法包括：

[0031] S101，调整CT强度阈值范围，获取CT强度值属于强度阈值范围内的待显示部位，并对 待显示部位进行标记着色，利用体绘制技术模拟等值面显示，显示原始数据体绘制视图。

[0032] 具体地，本发明提出了使用体绘制技术模拟等值面显示，并为数据标记着色的方法。

[0033] 设计如下传递函数，可使体绘制只显示大于阈值的部分：

[0034]

[0035] 即，CT强度值大于阈值的部分被设置为不透明，小于阈值的部分被设置为透明。调节阈 值即可模拟显示不同阈值下的等值面，体绘制是由显卡计算的，所以计算非常迅速，可实时 显示。

[0036] 同时为了为有标记数据显示多种颜色，本发明构造了一个特殊的传递函数：

[0037]

[0038]

[0039] 其中，I为CT强度值，T(A)、T(B)表示分别映射到两个颜色标签的阈值；且CT强度值 满足强度阈值范围的部分被设置为不透明，不满足强度阈值范围的部分被设置为透明。

[0040] 具体地，图3为显示多种颜色传递函数示意图：

[0041] 当需要将分割的区域以其他颜色标记时，只需将此区域连同其附近的透明区域的CT值 取相反数即可。将附近透明区域也取相反数的目的是利用体绘制的线性插值，使区域边界显 得更为平滑。

[0042] S102，通过标记工具确定原始数据体绘制视图中的待分割区域，使用填充工具对待分割 区域进行计算，生成初步分割结果，并使用腐蚀膨胀工具调整初步分割结果。

[0043] 具体地，本发明提供了若干方便的交互分割工具，使用这些工具可选中体数据中符合条 件的某些体素，选中后，可将这些体素①标记为分割、②从分割中排除或是③隐藏。

[0044] 作为本发明实施例的一个可选实施方式，分割工具可以包括：标记工具、填充工具、区 域锁定工具、腐蚀膨胀工具等。

[0045] 其中，作为本发明实施例的一个可选实施方式，通过标记工具确定原始数据体绘制视图 中的待分割区域包括：在原始数据体绘制视图中，通过拖动鼠标，标记拖动轨迹附近的区域 作为待分割区域，拖动轨迹终止于深度突变前的位置；其中，拖动轨迹为：在拖动的某一时 刻，视线通过鼠标指针的延长线与当前所显示的等值面产生若干交点，取其中最先相交的两 点的中点，将所有时刻的中点连接成的折线。具体地，标记工具用于选中体数据中鼠标指定 的区域。在原始数据体绘制视图中，使用鼠标沿待分割对象拖动，即可标记拖动轨迹附近的 区域。其中，拖动轨迹定义如下：在拖动的某一时刻，视线通过鼠标指针的延长线与当前所 显示的等值面将产生若干交点，取其中最先相交的两点的中点，所有时刻的这样的一系列中 点连接成的折线即为拖动轨迹。该工具标记出的区域是一个管状区域，管的轴线为上述拖动 轨迹，直径为上述最先相交两点所求出的一系列距离的中位值的1.5倍。特别地，当鼠标拾 取坐标的深度发生突变时，拖动轨迹将终止于深度突变前的位置。例如，开始时鼠标沿血管 标记，在某一时刻移出了此血管区域，那么，拾取点将转移至等值面的其他位置，或是无拾 取点，此时拾取深度发生了突变，拖动轨迹将终止并停止接受鼠标响应。该管状区域半径可 能大于实际欲标记区域的尺寸，从而将一些背景区域包含在内，所以，该工具将仅保留与拖 动轨迹相连通的区域。由此，可以通过标记工具确定待分割区域。

[0046] 作为本发明实施例的一个可选实施方式，使用填充工具对待分割区域进行计算，生成初 步分割结果包括：利用鼠标在原始数据体绘制视图中设置种子点，选中与种子点相关联的区 域，并利用洪水填充算法或者区域生长算法，计算初步分割结果。具体地，填充工具可快速 选中与种子点相关联的一个连续区域。种子点可直接在体绘制视图中使用鼠标点击选取，选 取的点是鼠标在等值面上拾取的点。通常，期望的种子点位置都位于等值面包围区域内部而 非表面，因此，本发明可调节种子点的深度，方法是滚动鼠标滚轮，在体绘制视图中会显示 x、y、z轴的参考线以表明种子点位置；二维视图亦会调整显示位置以便观察。例如可以分 别使用红、绿、蓝线表示种子点处平行于x、y、z坐标轴的直线。而判断与种子点“相关联” 的原则是洪水填充(flood fill)算法或区域生长(region growing)算法，根据使用场景不同可选择 不同填充原则。其中，洪水填充算法的对象是高于当前全局阈值的体素，所有与种子点相连 通的体素将被填充；区域生长算法与当前全局阈值无关，从种子点开始不断将与已填充区域 相连且灰度差异不大的体素包含进来，当没有更多满足条件的体素时生长停止。同时，区域 锁定工具通常与填充工具结合使用，以解决填充工具容易选中过多区域的问题。先使用其他 工具选中一个区域作为当前关注的作用范围，当该工具处于启用状态时，填充工具将被限制 为仅在作用范围内执行算法。由此，可以通过填充工具获取初步分割结果。

[0047] 进一步地，作为本发明实施例的一个可选实施方式，使用腐蚀膨胀工具调整初步分割结 果包括：调整与种子点相关联的区域的边缘，对与种子点相关联的区域执行三维的形态学腐 蚀/膨胀算法，调整初步分割结果。具体地，腐蚀膨胀工具可用于调整选中区域的边缘，对选 中区域执行三维的形态学腐蚀/膨胀算法。使用腐蚀算法可使选中区域边缘收缩，使用膨胀算 法可使选中区域边缘扩张。由此，可以通过腐蚀膨胀工具调整初步分割结果，使其更加平滑。

[0048] S103，在初步分割结果的基础上，使用标记工具选中待分割区域中未分割出的区域，生 成分割结果，显示分割体绘制视图。

[0049] 具体地，如有未分割出的区域，还可以继续选择未分割的区域，直至全部待分割区域中 全部选中完毕。

[0050] 同时，用户从模拟等值面仅能观察到体数据中等于某个阈值的体素分布，而难以观察其 他未显示的数据，为方便用户在分割过程中掌握这些信息，作为本发明实施例的一个可选实 施方式，显示原始数据体绘制视图的同时，还包括：显示原始数据体绘制视图的矢状面切片、 冠状面切片和横截面切片的二维视图，且在原始数据体绘制视图的矢状面切片、冠状面切片 和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或显示分割体绘制视 图的同时，还包括：显示分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视 图，且在分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于当前全 局阈值的区域进行着色显示。本发明除体绘制视图外提供了三个二维视图，分别用于显示矢 状面、冠状面、横截面切片的图像。二维视图中，高于当前全局阈值的区域(即当前体绘制视 图中被显示出的数据)例如以绿色显示，其余部分例如以灰度图显示。

[0051] 此外，原始数据体绘制视图和/或分割体绘制视图内也可显示这些切片图像，显示模式可 选择显示整张切片，也可仅显示高于阈值的区域。实际使用时，可根据使用场景灵活选择显 示模式。

[0052] 由此可见，通过本发明提供的CT图像交互分割方法，可以替代目前Seg3D、MITK工作 台等软件中针对CT图像的分割方法，以直观的交互手段简化分割过程，节省大量的分割时 间，更利于用户使用。

[0053] 以下提供一种针对某头颈部CT图像中血管的分割方案，分割过程为：

[0054] (1)调整合适阈值，使血管可以清晰显示；

[0055] (2)在体绘制界面中，单击鼠标拾取一条较粗的血管，使用填充工具的区域生长算法生 成初步分割，在体绘制界面显示切片，并依照切片使用腐蚀膨胀工具调整分割区域；

[0056] (3)找到尚未分割出的细小血管，使用标记工具选中这些血管；在体绘制界面与切片查 看界面可确认当前拾取点的准确性，如不准确可随时调整。

[0057] (4)完成分割，并生成等值面网格。

[0058] 图2示出了本发明实施例提供的CT图像交互分割装置的结构示意图，该CT图像交互分 割装置应用于上述方法，以下仅对CT图像交互分割装置的结构进行简单说明，其他未尽事 宜，请参照上述CT图像交互分割方法中的相关描述，参见图2，CT图像交互分割装置，包 括：

[0059] 显示模块201，用于调整CT强度阈值范围，获取CT强度值属于强度阈值范围内的待显 示部位，并对待显示部位进行标记着色，利用体绘制技术模拟等值面显示，显示原始数据体 绘制视图；

[0060] 初步分割模块202，用于通过标记工具确定原始数据体绘制视图中的待分割区域，使用 填充工具对待分割区域进行计算，生成初步分割结果，并使用腐蚀膨胀工具调整初步分割结 果；

[0061] 分割模块203，用于在初步分割结果的基础上，使用标记工具选中待分割区域中未分割 出的区域，生成分割结果，通知显示模块201显示分割体绘制视图。

[0062] 由此可见，通过本发明提供的CT图像交互分割装置，可以替代目前Seg3D、MITK工作 台等软件中针对CT图像的分割方法，以直观的交互手段简化分割过程，节省大量的分割时 间，更利于用户使用。

[0063] 由于用户从模拟等值面仅能观察到体数据中等于某个阈值的体素分布，而难以观察其他 未显示的数据，为方便用户在分割过程中掌握这些信息，作为本发明实施例的一个可选实施 方式，显示模块201，还用于在显示原始数据体绘制视图的同时，显示原始数据体绘制视图 的矢状面切片、冠状面切片和横截面切片的二维视图，且在原始数据体绘制视图的矢状面切 片、冠状面切片和横截面切片的二维视图中高于当前全局阈值的区域进行着色显示；和/或显 示分割体绘制视图的同时，显示分割体绘制视图的矢状面切片、冠状面切片和横截面切片的 二维视图，且在分割体绘制视图的矢状面切片、冠状面切片和横截面切片的二维视图中高于 当前全局阈值的区域进行着色显示。

[0064] 作为本发明实施例的一个可选实施方式，初步分割模块202具体用于通过如下方式通过 标记工具确定原始数据体绘制视图中的待分割区域：初步分割模块202，还用于在原始数据 体绘制视图中，通过拖动鼠标，标记拖动轨迹附近的区域作为待分割区域，拖动轨迹终止于 深度突变前的位置；其中，拖动轨迹为：在拖动的某一时刻，视线通过鼠标指针的延长线与 当前所显示的等值面产生若干交点，取其中最先相交的两点的中点，将所有时刻的中点连接 成的折线。由此，可以通过标记工具确定待分割区域。

[0065] 作为本发明实施例的一个可选实施方式，初步分割模块202具体用于通过如下方式使用 填充工具对待分割区域进行计算，生成初步分割结果：初步分割模块202，还用于利用鼠标 在原始数据体绘制视图中设置种子点，选中与种子点相关联的区域，并利用洪水填充算法或 者区域生长算法，计算初步分割结果。由此，可以通过填充工具获取初步分割结果。

[0066] 作为本发明实施例的一个可选实施方式，初步分割模块202具体用于通过如下方式使用 腐蚀膨胀工具调整初步分割结果：初步分割模块202，还用于调整与种子点相关联的区域的 边缘，对与种子点相关联的区域执行三维的形态学腐蚀/膨胀算法，调整初步分割结果。由此， 可以通过腐蚀膨胀工具调整初步分割结果，使其更加平滑。

[0067] 作为本发明实施例的一个可选实施方式，传递函数包括：

[0068]

[0069]

[0070] 其中，I为CT强度值，T(A)、T(B)表示分别映射到两个颜色标签的阈值；且CT强度值 满足强度阈值范围的部分被设置为不透明，不满足强度阈值范围的部分被设置为透明。

[0071] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。 因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的 形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储 介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形 式。

[0072] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和 /或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和 /或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指 令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一 个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流 程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

[0073] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工 作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制 造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指 定的功能。

[0074] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或 其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编 程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多 个方框中指定的功能的步骤。

[0075] 在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接 口和内存。

[0076] 存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易 失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。

[0077] 计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术 来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机 的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态 随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可 擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器 (CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他 磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中 的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和 载波。

[0078] 以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申 请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进 等，均应包含在本申请的权利要求范围之内。