一种ct图像扫描床去除方法及装置的制造方法
【专利摘要】本发明涉及一种CT图像扫描床去除方法及装置。该CT图像扫描床去除方法，包括：步骤a：通过图像处理设备主线程读取三维CT图像，统计CT设备内核数量，并初始化子算法；步骤b：图像处理设备主线程从输入的三维CT图像中提取二维扫描图像，通过共享内存，图像处理设备自动将二维扫描图像分配给内核，实现多线程并行处理对二维扫描图像进行去床操作；步骤c：图像处理设备结束并行操作，输出去床后的三维CT图像。本发明实施例的CT图像扫描床去除方法及装置采用的图像分割算法非常有效和精确，去除扫描床的同时不丢失身体部位的数据信息。
【专利说明】
一种CT图像扫描床去除方法及装置
技术领域
[0001] 本申请涉及图像分割技术领域，特别涉及一种CT图像扫描床去除方法及装置。
【背景技术】
[0002] 随着医学影像技术的发展，CT摄影技术逐渐成熟，被扫描的图像更加清晰，产生的 数据量也不断增大。现今较常用的CT图像的空间分辨率是512*512,仅一张扫描图像 (slice)就有26万多个像素点，若灰度值以8字节进行存储，其数据量达到2兆字节，而对于 一个CT体数据，其扫描的s 1 i ce数一般大于100张，也就是说一个常用的三维CT扫描图像的 数据量要超过200兆字节，待处理图像庞大的数据量和现有医学图像分割算法的局限性等 严重影响临床治疗的效率，图像分割的加速是临床诊断实时性的基础。
[0003] 图像分割中常见的加速方法主要有硬件加速和软件加速。硬件加速是利用高配置 设备的大内存、大容量、多CHJ等特点提升图像分割的速度，但是其缺陷在于，（1)需根据实 际的应用，设计硬件，增加设备成本，而且维修费用较高、难度较大；（2)现有的图像分割算 法局限性的影响，加速效果不明显。软件加速是通过深入研究图像分割算法原理，进行算法 本身的简化或创新，如减少内循环或降采样预处理图像等，但是其缺陷在于，（1)需要研究 算法的本质，从算法的基本原理出发重新编写代码，受算法复杂性、多样性的影响，极难实 现且耗时；（2)提速有限，如图像分割过程中不可避免的图像预处理、灰度统计或多层循环 等。这种方法适用于研究者进行图像分割算法创新，且局限性较大。
[0004] CT扫描床是用来配合扫描装置完成全方位的扫描，扫描床具有上下前后等移动功 能，在使用中根据不同的扫描位置来调整扫描床以达到合理的位置。但是在实践中，拍摄的 CT图像通常含有扫描床的影像，扫描床的影像会对CT图像产生干扰，影响临床诊断的准确 性，因此去除CT扫描床是CT图像处理的基础。目前，CT扫描床去除方法为通过具有内置去床 算法的CT设备实现。具有内置去床算法的CT设备是厂家根据设备中扫描床的型号特点，添 加去床程序，拍摄的结果图像直接为去床后的CT图像。但是采用具有内置去床算法的CT设 备由于内置去床算法不具有通用性，不同厂家之间的设备算法不同，不能通用，另外，去床 算法不可见，研究者和医生不能根据实际需要修改算法;另外，具有内置去床算法的CT设备 通常采用硬件加速或者软件加速，其加速效果不明显，去除扫描床影像耗时长，且去床效果 不好。
【发明内容】
[0005] 本申请提供了一种CT图像扫描床去除方法及装置，以解决现有技术中内置去床算 法不具有通用性、去除扫描床影像耗时长且去床效果不好的技术问题。
[0006] 为了解决上述问题，本发明提供了如下技术方案:一种CT图像扫描床去除方法，包 括：
[0007] 步骤a:通过图像处理设备主线程读取三维CT图像，统计CT设备内核数量，并初始 化子算法；
[0008] 步骤b:图像处理设备主线程从输入的三维CT图像中提取二维扫描图像，通过共享 内存，图像处理设备自动将二维扫描图像分配给内核，实现多线程并行处理对二维扫描图 像进行去床操作；
[0009] 步骤C:图像处理设备结束并行操作，输出去床后的三维CT图像。
[0010] 本发明实施例采取的技术方案还包括:所述步骤b具体包括：
[0011] 步骤bl:从输入的三维CT图像中提取二维扫描图像，读取二维扫描图像，对读取的 二维扫描图像进行分割；
[0012] 步骤b2:提取二维扫描图像中的目标区域图像信息；
[0013] 步骤b3:对提取的目标区域图像信息进行形态学开运算；
[0014] 步骤b4:获取二维扫描图像中的目标区域图像灰度信息；
[0015] 步骤b5:将各线程获取的二维扫描图像中的目标区域图像灰度信息进行结合，去 床扫描床信息。
[0016] 本发明实施例采取的技术方案还包括:在所述步骤bl中，对读取的二维扫描图像 进行分割采用大津阈值分割。
[0017] 本发明实施例采取的技术方案还包括:在所述步骤b2中，提取二维扫描图像中的 目标区域图像信息包括提取二维扫描图像中的身体部分信息。
[0018] 本发明实施例采取的技术方案还包括:在所述步骤b3中，获取二维扫描图像中的 目标区域图像灰度信息为:获取二维扫描图像中身体部分的灰度信息，去除CT图像扫描床 fg息。
[0019] 本发明实施例采取的另一技术方案为:一种CT图像扫描床去除装置，包括图像读 取模块、图像处理模块和图像输出模块，所述图像读取模块读取三维CT图像，统计CT设备内 核数量，并初始化子算法，所述图像处理模块从输入的三维CT图像中提取二维扫描图像，通 过共享内存，自动将二维扫描图像分配给内核，实现多线程并行处理对二维扫描图像进行 去床操作，所述图像输出模块用于结束并行操作，输出去床后的三维CT图像。
[0020] 本发明实施例采取的技术方案还包括:所述图像处理模块包括图像分割模块、图 像提取模块、图像运算模块、信息获取模块和图像结合模块，所述图像分割模块读取二维扫 描图像，对读取的二维扫描图像进行阈值分割，所述图像提取模块提取二维扫描图像中的 目标区域图像信息，所述图像运算模块对提取的目标区域图像信息进行形态学开运算，所 述信息获取模块获取二维扫描图像中的目标区域图像灰度信息，所述图像结合模块将各线 程获取的二维扫描图像中的目标区域图像灰度信息进行结合，去床扫描床信息。
[0021] 本发明实施例采取的技术方案还包括:所述图像分割模块对读取的二维扫描图像 进行阈值分割米用大津阈值分割。
[0022] 本发明实施例采取的技术方案还包括:所述图像提取模块提取二维扫描图像中的 目标区域图像信息具体为:提取二维扫描图像中的目标区域图像信息包括提取二维扫描图 像中的身体部分信息。
[0023] 本发明实施例采取的技术方案还包括:所述信息获取模块获取二维扫描图像中的 目标区域图像灰度信息具体为:获取二维扫描图像中身体部分的灰度信息，去除CT图像扫 描床信息。
[0024] 本发明实施例的CT图像扫描床去除方法及装置采用的图像分割算法非常有效和 精确，去除扫描床的同时不丢失身体部位信息；另外，本发明实施例的CT图描床去除方法及 装置大幅度提高扫描床去除速度，满足实时性要求。
【附图说明】
[0025] 图1是本发明实施例的CT图像扫描床去除方法的流程图；
[0026] 图2是本发明实施例的CT图像扫描床去除方法对二维扫描图像进行去床操作的流 程图；
[0027] 图3是本发明另一实施例的CT图像扫描床去除方法的流程图；
[0028] 图4是本发明实施例的CT图像扫描床去除装置的结构示意图；
[0029] 图5是本发明实施例的CT图像扫描床去除方法的实验结果图；
[0030] 图6是本发明实施例的CT图像扫描床去除方法三维数据分割的准确性度量。
【具体实施方式】
[0031] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对 本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不 用于限定本发明。
[0032] 请参阅图1，图1是本发明实施例的CT图像扫描床去除方法的流程图。本发明实施 例的CT图像扫描床去除方法包括：
[0033]步骤10:通过图像处理设备(主线程)读取三维CT图像，统计CT设备内核数量，并初 始化子算法；
[0034] 在步骤10中，读取三维CT图像的图像处理设备可以设置在CT设备中，也可以设置 在CT设备外或者独立于CT设备。
[0035]步骤20:图像处理设备主线程从输入的三维CT图像中提取二维扫描图像，通过共 享内存，图像处理设备自动将二维扫描图像分配给内核，实现多线程并行处理对二维扫描 图像进行去床操作；
[0036]步骤30:图像处理设备(主线程)结束并行操作，输出去床后的三维CT图像。
[0037]请参阅图2,图2是本发明实施例的CT图像扫描床去除方法对二维扫描图像进行去 床操作的流程图。本发明实施例的CT图像扫描床去除方法对二维扫描图像进行去床操作具 体包括：
[0038]步骤210:从输入的三维CT图像中提取二维扫描图像，读取二维扫描图像，对读取 的二维扫描图像进行大津(0TSU)阈值分割；
[0039] 在步骤210中，根据去床算法原理，对读取的二维扫描图像进行大津(0TSU)阈值分 害J。大津(0TSU)阈值分割是按图像的灰度特性，将图像分成背景和目标两部分。背景和目标 之间的类间方差越大，说明构成图像的两部分的差别越大，当部分目标错分为背景或部分 背景错分为目标都会导致两部分差别变小。因此，使类间方差最大的分割意味着错分概率 最小。
[0040] 步骤220:提取二维扫描图像中的目标区域图像信息；
[0041] 在步骤220中，提取二维扫描图像中的目标区域图像信息包括提取二维扫描图像 中的身体部分信息。
[0042]步骤230:对提取的目标区域图像信息进行形态学开运算；
[0043]在步骤230中，形态学主要是获取物体拓扑和结构信息，通过物体和结构元素相互 作用的某些运算，得到物体更本质的形态。在图像处理中的应用主要是:利用形态学的基本 运算，对图像进行观察和处理，从而达到改善图像质量的目的。图像形态学中的腐蚀和膨胀 能很好的对二值图像进行减噪处理。腐蚀的具体操作是：用一个结构元素(一般是3X3的大 小)扫描图像中的每一个像素，用结构元素中的每一个像素与其覆盖的像素做"与"操作，如 果都为1，则该像素为1，否则为〇。膨胀的具体操作是：用一个结构元素(一般是3 X 3的大小） 扫描图像中的每一个像素，用结构元素中的每一个像素与其覆盖的像素做"与"操作，如果 都为〇,则该像素为〇,否则为1。腐蚀的作用是消除物体边界点，使目标缩小，可以消除小于 结构元素的噪声点;膨胀的作用是将与物体接触的所有背景点合并到物体中，使目标增大， 可添补目标中的空洞。开运算是先腐蚀后膨胀的过程，可以消除图像上细小的噪声，并平滑 物体边界。
[0044]步骤240:获取二维扫描图像中的目标区域图像灰度信息；
[0045]在步骤240中，获取二维扫描图像中的目标区域图像灰度信息为:获取二维扫描图 像中身体部分的灰度信息，去除CT图像扫描床信息；
[0046] 步骤250:将各线程获取的二维扫描图像中的目标区域图像灰度信息进行结合，输 出三维去床CT扫描图像。
[0047] 请参阅图3,图3是本发明另一实施例的CT图像扫描床去除方法的流程图。本发明 实施例的CT图像扫描床去除方法可以针对并行CT扫描床进行去除，也可以针对非并行CT扫 描床进行去除。本发明实施例的CT图像扫描床去除方法如果针对非并行CT扫描床进行去 除，则具体方法包括以下步骤：
[0048]步骤40:通过CT设备读入含有扫描床信息的三维CT图像；
[0049]步骤50:根据去床算法原理，读取扫描图三维CT图像，对三维CT图像依次进行大津 (0TSU)阈值分割、提取前景图像区域(包括CT图像中的身体部位和扫描床）、形态学开运算 等图像分割过程；
[0050]其中，大津(0TSU)阈值分割是按图像的灰度特性，将图像分成背景和目标两部分。 背景和目标之间的类间方差越大，说明构成图像的两部分的差别越大，当部分目标错分为 背景或部分背景错分为目标都会导致两部分差别变小。因此，使类间方差最大的分割意味 着错分概率最小。形态学主要是获取物体拓扑和结构信息，通过物体和结构元素相互作用 的某些运算，得到物体更本质的形态。在图像处理中的应用主要是：利用形态学的基本运 算，对图像进行观察和处理，从而达到改善图像质量的目的。图像形态学中的腐蚀和膨胀能 很好的对二值图像进行减噪处理。腐蚀的具体操作是：用一个结构元素（一般是3X3的大 小)扫描图像中的每一个像素，用结构元素中的每一个像素与其覆盖的像素做"与"操作，如 果都为1，则该像素为1，否则为〇。膨胀的具体操作是：用一个结构元素(一般是3 X 3的大小） 扫描图像中的每一个像素，用结构元素中的每一个像素与其覆盖的像素做"与"操作，如果 都为〇,则该像素为〇,否则为1。腐蚀的作用是消除物体边界点，使目标缩小，可以消除小于 结构元素的噪声点;膨胀的作用是将与物体接触的所有背景点合并到物体中，使目标增大， 可添补目标中的空洞。开运算是先腐蚀后膨胀的过程，可以消除图像上细小的噪声，并平滑 物体边界。
[0051 ]步骤60:获取分割结果图，输出三维去床CT扫描图像。
[0052] 请参阅图4,图4是本发明实施例的CT图像扫描床去除装置的结构示意图。本发明 实施例的CT图像扫描床去除装置包括图像读取模块、图像处理模块和图像输出模块。图像 读取模块读取三维CT图像，统计CT设备内核数量，并初始化子算法。图像处理模块从输入的 三维CT图像中提取二维扫描图像，通过共享内存，自动将二维扫描图像分配给内核，实现多 线程并行处理对二维扫描图像进行去床操作。图像输出模块用于结束并行操作，输出去床 后的三维CT图像。图像处理模块包括图像分割模块、图像提取模块、图像运算模块、信息获 取模块和图像结合模块。图像分割模块读取二维扫描图像，对读取的二维扫描图像进行大 津(0TSU)阈值分割。图像分割模块根据去床算法原理，对读取的二维扫描图像进行大津 (0TSU)阈值分割。大津(0TSU)阈值分割是按图像的灰度特性，将图像分成背景和目标两部 分。背景和目标之间的类间方差越大，说明构成图像的两部分的差别越大，当部分目标错分 为背景或部分背景错分为目标都会导致两部分差别变小。因此，使类间方差最大的分割意 味着错分概率最小。
[0053] 图像提取模块提取二维扫描图像中的目标区域图像信息，提取二维扫描图像中的 目标区域图像信息包括提取二维扫描图像中的身体部分信息。
[0054]图像运算模块对提取的目标区域图像信息进行形态学开运算。形态学主要是获取 物体拓扑和结构信息，通过物体和结构元素相互作用的某些运算，得到物体更本质的形态。 在图像处理中的应用主要是:利用形态学的基本运算，对图像进行观察和处理，从而达到改 善图像质量的目的。图像形态学中的腐蚀和膨胀能很好的对二值图像进行减噪处理。腐蚀 的具体操作是：用一个结构元素(一般是3X3的大小）扫描图像中的每一个像素，用结构元 素中的每一个像素与其覆盖的像素做"与"操作，如果都为1，则该像素为1，否则为0。膨胀的 具体操作是：用一个结构元素(一般是3X3的大小）扫描图像中的每一个像素，用结构元素 中的每一个像素与其覆盖的像素做"与"操作，如果都为〇，则该像素为〇，否则为1。腐蚀的作 用是消除物体边界点，使目标缩小，可以消除小于结构元素的噪声点;膨胀的作用是将与物 体接触的所有背景点合并到物体中，使目标增大，可添补目标中的空洞。开运算是先腐蚀后 膨胀的过程，可以消除图像上细小的噪声，并平滑物体边界。
[0055] 信息获取模块获取二维扫描图像中的目标区域图像灰度信息。信息获取模块获取 二维扫描图像中身体部分的灰度信息，去除CT图像扫描床信息。
[0056] 图像结合模块将各线程获取的二维扫描图像中的目标区域图像灰度信息进行结 合，去床扫描床彳目息。
[0057]对本发明实施例的CT图像扫描床去除方法进行临床实验验证，可以理解，临床实 验验证用于进一步说明本发明的有益效果，对发明的实施方式和保护范围并不构成限制。
[0058] 请参阅图5,图5是本发明实施例的CT图像扫描床去除方法的实验结果图。图5第一 行给出原始的带有扫描床的CT图像，其中（A)为三维图像，可以清晰的看到身体左侧的薄 板，（B)为轴面切向图像，扫描床近似为两条弯曲的曲线，（C)为矢状面图像，扫描床近似为 一条垂线与身体近乎平行；图5第二行给出经本发明算法去床后的结果图，从视觉上，本发 明提出的扫描床去除程序能够很好地去除扫描床，且几乎没有误侵蚀现象。
[0059] 每个断层图像的平均消耗时间公式如下：
其中tcl为第i个断层图像分割所需要的时间。
[0061 ]图像分割准确率参数计算公式如下：
[0063]图像分割误差率参数计算公式如下：
[0064]假阳性(false positive，FP):表示本发明提出的算法未能成功去除扫描床的误 差率。
[0066] 假阴性(false negative，FN):表示本发明提出的算法对身体部位(mask)的误侵 蚀率。
[0068]其中| ? |用来统计三维数据内的点个数，G为手动分割的金标准，而S为分割结果。
[0069] 请参阅图6,图6是本发明实施例的CT图像扫描床去除方法三维数据分割的准确性 度量。整体上来看，分割的平均精度达到99%，本发明实施例的CT图像扫描床去除方法分割 算法非常有效和精确;假阳性和假阴性的平均值分别为0.4%和1.63%，说明本发明的去床 算法能精确的去除扫描床，而且几乎不损伤身体部位。
[0070] 本发明实施例的CT图像扫描床去除方法软件实现为Visual Studio 2010、?，采 用(^1^进行加速实现。实验机器为8核1加€1_(^^(了1〇，主频3.76抱，内存166。
[0072] 上述表格比较了手动分割时间，没有引入本加速策略的分割运行时间，以及引入 本加速策略的时间。通过分析发现，本发明实施例的CT图像扫描床去除方法能在0.29秒内 对单张分辨率为[512,512]的图像进行去床操作，其去床的速度是未加速的2.72倍，是手动 分割的400倍以上，本发明实施例的CT图像扫描床去除方法大幅度提高扫描床去除速度，满 足实时性要求。
[0073] 本发明实施例的CT图像扫描床去除方法及装置采用的图像分割算法非常有效和 精确，去除扫描床的同时不丢失身体部位信息；另外，本发明实施例的CT图描床去除方法及 装置大幅度提高扫描床去除速度，满足实时性要求。
[0074]对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。 对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的 一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明 将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一 致的最宽的范围。
【主权项】
1. 一种CT图像扫描床去除方法，包括： 步骤a:通过图像处理设备主线程读取三维CT图像，统计CT设备内核数量，并初始化子 算法； 步骤b:图像处理设备主线程从输入的三维CT图像中提取二维扫描图像，通过共享内 存，图像处理设备自动将二维扫描图像分配给内核，实现多线程并行处理对二维扫描图像 进行去床操作； 步骤C:图像处理设备结束并行操作，输出去床后的三维CT图像。2. 根据权利要求1所述的CT图像扫描床去除方法，其特征在于，所述步骤b具体包括： 步骤b 1:从输入的三维CT图像中提取二维扫描图像，读取二维扫描图像，对读取的二维 扫描图像进行分割； 步骤b2:提取二维扫描图像中的目标区域图像信息； 步骤b3:对提取的目标区域图像信息进行形态学开运算； 步骤b4:获取二维扫描图像中的目标区域图像灰度信息； 步骤b5:将各线程获取的二维扫描图像中的目标区域图像灰度信息进行结合，去床扫 描床信息。3. 根据权利要求2所述的CT图像扫描床去除方法，其特征在于，在所述步骤bl中，对读 取的二维扫描图像进行分割采用大津阈值分割。4. 根据权利要求2所述的CT图像扫描床去除方法，其特征在于，在所述步骤b2中，提取 二维扫描图像中的目标区域图像信息包括提取二维扫描图像中的身体部分信息。5. 根据权利要求1所述的CT图像扫描床去除方法，其特征在于，在所述步骤b3中，获取 二维扫描图像中的目标区域图像灰度信息为:获取二维扫描图像中身体部分的灰度信息， 去除CT图像扫描床彳目息。6. -种CT图像扫描床去除装置，其特征在于:包括图像读取模块、图像处理模块和图像 输出模块，所述图像读取模块读取三维CT图像，统计CT设备内核数量，并初始化子算法，所 述图像处理模块从输入的三维CT图像中提取二维扫描图像，通过共享内存，自动将二维扫 描图像分配给内核，实现多线程并行处理对二维扫描图像进行去床操作，所述图像输出模 块用于结束并行操作，输出去床后的三维CT图像。7. 根据权利要求6所述的CT图像扫描床去除装置，其特征在于，所述图像处理模块包括 图像分割模块、图像提取模块、图像运算模块、信息获取模块和图像结合模块，所述图像分 割模块读取二维扫描图像，对读取的二维扫描图像进行阈值分割，所述图像提取模块提取 二维扫描图像中的目标区域图像信息，所述图像运算模块对提取的目标区域图像信息进行 形态学开运算，所述信息获取模块获取二维扫描图像中的目标区域图像灰度信息，所述图 像结合模块将各线程获取的二维扫描图像中的目标区域图像灰度信息进行结合，去床扫描 床信息。8. 根据权利要求7所述的CT图像扫描床去除装置，其特征在于，所述图像分割模块对读 取的二维扫描图像进行阈值分割采用大津阈值分割。9. 根据权利要求7或8所述的CT图像扫描床去除装置，其特征在于，所述图像提取模块 提取二维扫描图像中的目标区域图像信息具体为:提取二维扫描图像中的目标区域图像信 息包括提取二维扫描图像中的身体部分信息。10.根据权利要求7或8所述的CT图像扫描床去除装置，其特征在于，所述信息获取模块 获取二维扫描图像中的目标区域图像灰度信息具体为:获取二维扫描图像中身体部分的灰 度ig息，去除CT图像扫描床彳g息D
【文档编号】G06T7/00GK105931251SQ201610319007
【公开日】2016年9月7日
【申请日】2016年5月12日
【发明人】余绍德, 陈璐明, 姬治华, 江帆, 伍世宾, 谢耀钦, 王磊
【申请人】中国科学院深圳先进技术研究院