Ct成像中定位物体的方法以及设备的制作方法
【专利摘要】提出了一种在CT成像中定位多个感兴趣物体的方法和设备。利用三幅近似垂直视角下的投影图像，实现感兴趣物体的三维空间定位。该方法对于扫描视野中存在多个感兴趣物的情况，仍然能够快速的确定CT图像中的感兴趣区物体的位置，并且无需对CT图像进行预重建，算法中也不涉及迭代步骤，方法快速、有效，能够很好地适用于工程应用。
【专利说明】CT成像中定位物体的方法以及设备
【技术领域】
[0001]本发明的实施例主要涉及辐射成像，更具体地，涉及在CT成像中定位物体的方法以及设备。
【背景技术】
[0002]自从1972年Hounsf ield发明了第一台CT机，CT技术给医学诊断和工业无损检测带来了革命性的影响，CT已经成为医疗、生物、航空航天、国防等行业重要的检测手段之一。随着技术的进步，CT扫描模式和成像方法也在不断地改进，三维锥束CT已经成为研究和应用的主流。X射线锥束CT已经在医学临床、安全检查、无损检测等领域得到了广泛的应用，特别是在医学临床诊断中，CT已经成为不可或缺的检查手段之一。
[0003]1989年，螺旋CT开始投入医学临床应用，由于螺旋CT的巨大优势，使得它逐步替代了以前的断层CT，螺旋CT相对于断层CT的优势在于:螺旋CT可以连续不间断地采集投影数据，并通过专门设计的重建算法得到物体的三维体数据，使得CT扫描的时间大大缩短，提供了重建图像的Z轴分辨率，减少了运动伪迹。螺旋CT在临床应用中取得了巨大的成功，很快，螺旋CT使用的探测器从单层发展到双层、四层、8层、16层，并逐渐发展到32层、64层、128层，到2007年Toshiba公司已经率先推出了 320层的螺旋CT。螺旋飞速发展的同时，另外一项技术也在悄然进步:平板探测器技术。不同于多层螺旋CT使用的很多独立探测器模块拼装成的阵列探测器，平板探测器直接使用整块的大面积闪烁体，后面封装大规模光敏单元阵列(例如(XD、CM0S、TFT等)，然后经过A/D转换得到X射线强度数据。平板探测器技术的发展促使一类新的锥束CT(CBCT)得以面世，使用了平板探测器的CBCT系统一般只需旋转一周即可完成较大区域(例如30cm*30cm)的CT扫描，重建出该扫描视野(FOV)范围内的三维CT图像。
[0004]在安全检查和工业无损检查领域，CT技术在近年来也获得了快速的发展，例如基于双能技术安检CT由于具备了很好的区分物质的能力，已经得到了安检领域的认可，正在逐步推广；而针对工业无损检测领域的工业CT也在空间分辨率、密度分辨率方面取得了较大的进步。
[0005]在对CT图像重建的过程中，如何准确地确定例如金属物体之类感兴趣物体的位置，也是一个重要的研究课题。
【发明内容】
[0006]考虑到现有技术的一个或多个问题，提出了 CT成像中定位物体的方法以及设备。
[0007]根据本发明的实施例，一种在CT成像中定位多个感兴趣物体的方法，包括步骤:
[0008]利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像；
[0009]在垂直于所述观察方向的锥束投影中选择两幅投影图像，即第二投影图像和第三投影图像，所述第二投影图像和所述第三投影图像基本上正交；[0010]确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置；
[0011]基于多个感兴趣物体在所述第一投影图像、第二投影图像和所述第三投影图像中的位置计算每个感兴趣物体在三维空间中的位置。
[0012]根据本发明的实施例，计算第一投影平面的步骤包括:
[0013]利用投影数据和雷当变换数值的关系，基于投影数据计算出与第一投影图像对应的雷当数据的偏导数；
[0014]利用滤波反投影算法和所述雷当数据的偏导数计算出第一投影图像。
[0015]根据本发明的实施例，选择两幅投影图像的步骤包括:
[0016]基于所述第一投影图像选择第二投影图像和第三投影图像，以使得第二投影图像和第三投影图像中的多个感兴趣物体之间的重叠区域最小。
[0017]根据本发明的实施例，基于所述第一投影图像选择第二投影图像和第三投影图像的步骤包括:
[0018]对所述第一投影图像进行分割，得到只包含感兴趣区域信息的二值化图像；
[0019]对所述二值化图像进行扇束的前向投影，得到扇束投影的正弦图，此处所用到的扇束的扇角需要等于锥束系统中心层对应光源靶点的张角；
[0020]对于正弦图的每列，通过寻峰算法计算峰的个数；
[0021]对于所有峰值个数等于第一投影图像中感兴趣物体个数的投影角度，选择二者之间角度为90度的投影角度，从而确定第二投影图像和第三投影图像。
[0022]根据本发明的实施例，确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置的步骤包括:
[0023]在第一投影图像中对各个感兴趣物体进行分割，确定每个感兴趣物体区域在第一图像中的重心；
[0024]在第二投影图像和第三投影图像中对各个感兴趣物体区域进行分割，并且确定每个感兴趣物体在第二投影图像和第三投影图像中的重心。
[0025]根据本发明的实施例，一种在CT成像中定位多个感兴趣物体的设备，包括:
[0026]利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像的装置；
[0027]在垂直于所述观察方向的锥束投影中选择两幅投影图像，即第二投影图像和第三投影图像，所述第二投影图像和所述第三投影图像基本上正交的装置；
[0028]确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置的装置；
[0029]基于多个感兴趣物体在所述第一投影图像、第二投影图像和所述第三投影图像中的位置计算每个感兴趣物体在三维空间中的位置的装置。
[0030]通过上述实施例，能够从CT投影数据中确定感兴趣物体的位置。
【专利附图】
【附图说明】
[0031]为了更好的理解本发明，将根据以下附图对本发明进行详细描述:
[0032]图1是根据本发明实施方式的CT设备的结构示意图；[0033]图2示出了如图1所示的计算机数据处理器的结构框图；
[0034]图3示出了根据本发明第一实施方式的控制器的结构框图；
[0035]图4示出了根据本发明实施例的定位物体的方法过程的示意图；
[0036]图5是描述CBCT系统的示意图；
[0037]图6是描述通过投影数据合成俯视角平行束投影图的算法的示意图；
[0038]图7示出了是定位所有备选MAP的示意图；
[0039]图8是描述由俯视图选择两幅正交的水平投影图；
[0040]图9示出了根据本发明另一实施例的消除CT成像中的伪影的方法流程图；
[0041]图10示出了伪影消除的实验结果。
【具体实施方式】
[0042]下面将详细描述本发明的具体实施例，应当注意，这里描述的实施例只用于举例说明，并不用于限制本发明。在以下描述中，为了提供对本发明的透彻理解，阐述了大量特定细节。然而，对于本领域普通技术人员显而易见的是:不必采用这些特定细节来实行本发明。在其他实例中，为了避免混淆本发明，未具体描述公知的电路、材料或方法。
[0043]在整个说明书中，对“ 一个实施例”、“实施例”、“ 一个示例”或“示例”的提及意味着:结合该实施例或示例描述的特定特征、结构或特性被包含在本发明至少一个实施例中。因此，在整个说明书的各个地方出现的短语“在一个实施例中”、“在实施例中”、“一个示例”或“示例”不一定都指同一实施例或示例。此外，可以以任何适当的组合和、或子组合将特定的特征、结构或特性组合在一个或多个实施例或示例中。此外，本领域普通技术人员应当理解，在此提供的示图都是为了说明的目的，并且示图不一定是按比例绘制的。应当理解，当称“元件”“连接到”或“耦接”到另一元件时，它可以是直接连接或耦接到另一元件或者可以存在中间元件。相反，当称元件“直接连接到”或“直接耦接到”另一元件时，不存在中间元件。
[0044]根据本发明的一些实施例，能够在CT成像中定位多个感兴趣物体。利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像。然后在垂直于所述观察方向的锥束投影中选择两幅投影图像，即第二投影图像和第三投影图像，所述第二投影图像和所述第三投影图像基本上正交。确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置。基于多个感兴趣物体在所述第一投影图像、第二投影图像和所述第三投影图像中的位置计算每个感兴趣物体在三维空间中的位置。
[0045]根据本发明的另外实施例，能够消除CT成像中的伪影。从投影数据中计算所述多个特定物体在三个彼此基本上垂直的投影图上的位置从而确定所述多个特定物体在三维空间中的位置。利用所述多个特定物体的重心坐标对所有角度下的投影图像进行物体分害I]。用分割了该多个特定物体后的投影数据对该多个特定物体的投影数据进行修复。利用修复后的投影数据进行CT重建获得不包含该多个特定物体的CT图像。
[0046]图1是根据本发明实施方式的CT设备的结构示意图。如图1所示，根据本实施方式的CT设备包括:发出检查用X射线的射线源10，诸如X光机；承载机构40，其承载被检物品围绕轴Z转动，并且可以升降，使得被检物体进入检测区域，从而由射线源10发出的射线能够透过被检物体；探测和采集装置30，它是具有整体模块结构的探测器及数据采集器，例如平板探测器，用于探测透射被检物品的射线，获得模拟信号，并且将模拟信号转换成数字信号，从而输出被检查物体针对X射线的投影数据；控制器50，它用于控制整个系统的各个部分同步工作；以及计算机数据处理器60，它用来处理由数据采集器采集的数据，对数据进行处理并重建，输出结果。
[0047]如图1所示，射线源10置于可放置被检物体的承载机构40 —侧，探测和采集装置30置于承载机构40的另一侧，包括探测器和数据采集器，用于获取被检物品投射数据和/或多角度投影数据。数据采集器中包括数据放大成形电路，它可工作于(电流)积分方式或脉冲(计数)方式。探测和采集装置30的数据输出电缆与计算机数据处理器60连接，根据触发命令将采集的数据存储在计算机数据处理器60中。
[0048]另外，检查设备还可以包括由金属制成的筒状物体通道20，它设置在承载机构40上，被检物品放置在被检物体通道中。
[0049]图2示出了如图1所示的计算机数据处理器60的结构框图。如图2所示，数据采集器所采集的数据通过接口单元68和总线64存储在存储器61中。只读存储器(ROM) 62中存储有计算机数据处理器的配置信息以及程序。随机存取存储器(RAM)63用于在处理器66工作过程中暂存各种数据。另外，存储器61中还存储有用于进行数据处理的计算机程序。内部总线64连接上述的存储器61、只读存储器62、随机存取存储器63、输入装置65、处理器66、显示装置67和接口单元68。
[0050]在用户通过诸如键盘和鼠标之类的输入装置65输入的操作命令后，计算机程序的指令代码命令处理器66执行预定的数据处理算法，在得到数据处理结果之后，将其显示在诸如LCD显示器之类的显示装置67上，或者直接以诸如打印之类硬拷贝的形式输出处理结果。
[0051]图3示出了根据本发明实施方式的控制器的结构框图。如图3所示，控制器50包括:控制单元51，根据来自计算机60的指令，来控制射线源10、承载机构40和探测和采集装置30 ;触发信号产生单元52，用于在控制单元的控制下产生用来触发射线源10、探测和采集装置30以及承载机构40的动作的触发命令；第一驱动电机55，它在根据触发信号产生单元52在控制单元51的控制下产生的触发命令驱动承载机构40上升或者下降；高度信息获取单元53，它随着承载机构40的运动，向控制单元51反馈承载机构的高度信息；第二驱动电机56，它根据触发信号产生单元52在控制单元51的控制下产生的触发命令驱动承载机构40旋转；角度信息获取单元54，它在承载机构40旋转过程中获取承载机构40的旋转角度，反馈给控制单元51。根据本发明的实施方式，上述的高度信息获取单元53和角度信息获取单元54都是光电码盘，它具备抗干扰的优点。
[0052]下面详细介绍基于三视图的定位CT成像中特定物体(感兴趣物体)的方法。图4示出了根据本发明实施例的定位物体的方法过程的示意图。
[0053]设CBCT的系统如下:重建区域以三维直角坐标系Oxyz描述，其中坐标原点O为系统的旋转中心。S(i3)为锥形束射线源所在的位置，其中β为系统的旋转角度值。平板探测器30置于旋转中心另一侧，并随光源同步转动。探测器上的投影数据以Ρ(β，u，ν)表示，其中U、ν为平板探测器上的直角坐标。CBCT系统示意图如图5所示。
[0054]在步骤S41:利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像。
[0055]通过采集的锥束CT(CBCT)投影数据，计算成像视野内的俯视角平行束投影图。例如利用投影数据和雷当变换数值的关系，基于投影数据计算出与第一投影图像对应的雷当数据的偏导数。例如，可以通过公式(I)计算俯视图(本例中对应于第一投影图像)对应的雷当值(Radon)的偏导数:
【权利要求】
1.一种在CT成像中定位多个感兴趣物体的方法，包括步骤: 利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像； 在垂直于所述观察方向的锥束投影中选择两幅投影图像，即第二投影图像和第三投影图像，所述第二投影图像和所述第三投影图像基本上正交； 确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置； 基于多个感兴趣物体在所述第一投影图像、第二投影图像和所述第三投影图像中的位置计算每个感兴趣物体在三维空间中的位置。
2.如权利要求1所述的方法，其中，计算第一投影平面的步骤包括: 利用投影数据和雷当变换数值的关系，基于投影数据计算出与第一投影图像对应的雷当数据的偏导数； 利用滤波反投影算法和所述雷当数据的偏导数计算出第一投影图像。
3.如权利要求1所述的方法，其中，选择两幅投影图像的步骤包括: 基于所述第一投影图像选择第二投影图像和第三投影图像，以使得第二投影图像和第三投影图像中的多个感兴趣物体之间的重叠区域最小。
4.如权利要求3所述的方法，其中，基于所述第一投影图像选择第二投影图像和第三投影图像的步骤包括: 对所述第一投影图像进行分割，`得到只包含感兴趣区域信息的二值化图像； 对所述二值化图像进行扇束的前向投影，得到扇束投影的正弦图，此处所用到的扇束的扇角需要等于锥束系统中心层对应光源靶点的张角； 对于正弦图的每列，通过寻峰算法计算峰的个数； 对于所有峰值个数等于第一投影图像中感兴趣物体个数的投影角度，选择二者之间角度为90度的投影角度，从而确定第二投影图像和第三投影图像。
5.如权利要求1所述的方法，其中，确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置的步骤包括: 在第一投影图像中对各个感兴趣物体进行分割，确定每个感兴趣物体区域在第一图像中的重心； 在第二投影图像和第三投影图像中对各个感兴趣物体区域进行分割，并且确定每个感兴趣物体在第二投影图像和第三投影图像中的重心。
6.一种在CT成像中定位多个感兴趣物体的设备，包括: 利用投影合成算法从投影数据中计算出其观察方向垂直于CT断层平面的第一投影图像的装置； 在垂直于所述观察方向的锥束投影中选择两幅投影图像的装置，即选择第二投影图像和第三投影图像，所述第二投影图像和所述第三投影图像基本上正交； 确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置的装置； 基于多个感兴趣物体在所述第一投影图像、第二投影图像和所述第三投影图像中的位置计算每个感兴趣物体在三维空间中的位置的装置。
7.如权利要求6所述的设备，其中，计算第一投影平面的装置包括: 利用投影数据和雷当变换数值的关系，基于投影数据计算出与第一投影图像对应的雷当数据的偏导数的装置； 利用滤波反投影算法和所述雷当数据的偏导数计算出第一投影图像的装置。
8.如权利要求6所述的设备，其中，选择两幅投影图像的装置包括: 基于所述第一投影图像选择第二投影图像和第三投影图像，以使得第二投影图像和第三投影图像中的多个感兴趣物体之间的重叠区域最小的装置。
9.如权利要求8所述的设备，其中，基于所述第一投影图像选择第二投影图像和第三投影图像的装置包括: 对所述第一投影图像进行分割，得到只包含感兴趣区域信息的二值化图像的装置； 对所述二值化图像进行扇束的前向投影，得到扇束投影的正弦图的装置，此处所用到的扇束的扇角需要等于锥束系统中心层对应光源靶点的张角； 对于正弦图的每列，通过寻峰算法计算峰的个数的装置； 对于所有峰值个数等于第一投影图像中感兴趣物体个数的投影角度，选择二者之间角度为90度的投影角度，从而确定第二投影图像和第三投影图像的装置。
10.如权利要求6所述的设备，其中，确定多个感兴趣物体在所述第一投影图像、所述第二投影图像和所述第三投影图像中的位置的装置包括: 在第一投影图像中对各个感兴趣物体进行分割，确定每个感兴趣物体区域在第一图像中的重心的装置； 在第二投影图像和第三投影图像中对各个感兴趣物体区域进行分割，并且确定每个感兴趣物体在第二投影图像和第三投影图像中的重心的装置。
【文档编号】G01V5/00GK103713329SQ201210375312
【公开日】2014年4月9日 申请日期:2012年9月29日 优先权日:2012年9月29日
【发明者】李亮, 张丽, 陈志强, 赵自然, 邢宇翔, 肖永顺, 王清礼 申请人:清华大学, 同方威视技术股份有限公司