CN103356218A - 一种x射线计算机断层扫描方法和系统 - Google Patents

一种x射线计算机断层扫描方法和系统 Download PDF

Info

Publication number
CN103356218A
CN103356218A CN2012101015005A CN201210101500A CN103356218A CN 103356218 A CN103356218 A CN 103356218A CN 2012101015005 A CN2012101015005 A CN 2012101015005A CN 201210101500 A CN201210101500 A CN 201210101500A CN 103356218 A CN103356218 A CN 103356218A
Authority
CN
China
Prior art keywords
interest
tomography
angle
region
roi
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101015005A
Other languages
English (en)
Other versions
CN103356218B (zh
Inventor
张纪庄
陆静
田毅
董建威
王薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Shanghai Medical Equipment Ltd
Original Assignee
Siemens Shanghai Medical Equipment Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Shanghai Medical Equipment Ltd filed Critical Siemens Shanghai Medical Equipment Ltd
Priority to CN201210101500.5A priority Critical patent/CN103356218B/zh
Publication of CN103356218A publication Critical patent/CN103356218A/zh
Application granted granted Critical
Publication of CN103356218B publication Critical patent/CN103356218B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种X射线计算机断层扫描扫描方法和系统。方法包括:确定三维感兴趣区的粗略边界,并计算得到当前三维感兴趣区;对第一个断层采用完全扫描法扫描,从得到的图像中获取该断层中感兴趣区的精确边界信息;利用该精确边界信息对当前三维感兴趣区中相应断层中的感兴趣区的几何信息进行校正;基于校正后的当前三维感兴趣区估算得到下一个断层中感兴趣区的几何信息,并确定该断层中扫描所需的最优放射起始角和最小扫描角度,基于此对该断层中的感兴趣区进行部分扫描。该方案能够在保证三维ROI的重建图像质量的前提下有效减少X射线的总放射剂量。

Description

一种X射线计算机断层扫描方法和系统
技术领域
本发明涉及医学成像领域,尤其是一种X射线计算机断层扫描方法、以及一种X射线计算机断层扫描系统。
背景技术
在进行X射线计算机断层扫描(Computed Tomography,CT)过程中,通常将待检对象需要进行高质量图像重建的器官或者组织称为感兴趣区(region of interest,ROI),该ROI为三维立体区域。基于现有技术,无论希望扫描并重建的ROI位于待检对象(如病人)的哪个部位或其大小如何,都要根据固定的XY准直器(phi-collimator)开放程度、固定的放射量和某种扫描方法决定的固定扫描角度,对待检对象进行逐一断层(slide)的放射扫描。
在传统的完全扫描法中,对每个断层都要进行扫描角度(irradiation duration)γ为360°的完全扫描。在部分扫描法中,采用某固定扫描角度γ(例如240°)来对该断层中的测量域(Field of Measurement,,FOM)进行扫描,FOM决定了直接对待检对象施加的放射剂量和X射线照射区域的大小。而在改进的部分扫描法中(参见上海西门子医疗器械有限公司于2010年12月19日申请的,申请号为201010286350.0,发明名称为“一种X射线计算机断层扫描系统和方法”的中国发明专利申请),在该断层中要进行扫描角度
Figure BDA0000149419260000011
的放射扫描,其中
Figure BDA0000149419260000012
是由待检对象的尺寸大小所确定的一个角度值。
由此可见,根据现有的CT扫描技术,即使三维ROI仅为待检对象很小的一个器官或者组织,总的放射剂量也不能减少。
发明内容
有鉴于此,本发明一方面提出了一种X射线计算机断层扫描方法,另一方面提出了一种X射线计算机断层扫描系统,用以在保证三维ROI的重建图像质量的前提下显著地减小每个断层的扫描角度,从而有效地减少X射线的总放射剂量。
本发明提出的X射线计算机断层扫描方法,包括:
A、根据预扫描出的定位图像,确定待检对象中三维感兴趣区的粗略边界;
B、根据所述粗略边界,计算得到当前三维感兴趣区;
C、确定扫描方向;
D、对所述扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像;
E、将第一个断层作为当前断层,从所述完全扫描图像中获取当前断层中感兴趣区的精确边界信息;
F、利用所述精确边界信息对当前三维感兴趣区中当前断层中的感兴趣区的几何信息进行校正,并基于校正后的当前三维感兴趣区通过插值算法估算得到下一个断层中感兴趣区的几何信息;
G、将下一个断层作为当前断层,根据所述当前断层中感兴趣区的几何信息,确定当前断层中扫描所需的最优放射起始角α0opt和最小扫描角度γmin;
H、根据所述最优放射起始角α0opt和最小扫描角度γmin对所述当前断层中的感兴趣区进行部分扫描,得到当前断层的部分扫描图像;
J、若三维感兴趣区未扫描完成,则从所述部分扫描图像中获取当前断层中感兴趣区的精确边界信息;并返回执行步骤F。
在本发明的一个实施方式中,所述粗略边界由可识别边界和宽泛估计边界构成。
在本发明的一个实施方式中,所述确定扫描方向为:将三维感兴趣区中对应所述可识别边界的一侧作为所述三维感兴趣区的扫描起始侧。
在本发明的一个实施方式中,该方法进一步包括:建立感兴趣区数据库;
所述步骤B之前进一步包括:检索所述感兴趣区数据库,在所述感兴趣区数据库中存在所述待检对象的所述三维感兴趣区的几何信息时,利用所述三维感兴趣区的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。
在本发明的一个实施方式中,所述步骤B之前进一步包括:在所述感兴趣区数据库中不存在所述待检对象的所述三维感兴趣区的几何信息时,利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维感兴趣区的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。
在本发明的一个实施方式中,步骤G中所述根据当前断层中感兴趣区的几何信息,确定当前断层中扫描所需的最优放射起始角α0opt和最小扫描角度γmin包括:
根据当前断层中感兴趣区的几何信息,计算放射起始角α0与放射附加角δ之间的关系曲线;
根据所述关系曲线,确定当前断层中的最优放射起始角α0opt及其对应的最小放射附加角δmin;
根据所述最小放射附加角δmin,计算当前断层中的最小扫描角度γmin。
其中,所述根据当前断层中感兴趣区的几何信息,计算放射起始角α0与放射附加角δ之间的关系曲线包括:
在0°-360°范围内逐一给定一放射起始角;
对每个给定的放射起始角α0,从所述断层中X射线球管运行圆周轨迹上角度为给定放射起始角α0的起始点向XY扫描平面上的所述感兴趣区引两条外切线,并与所述球管运行圆周轨迹相交于两个交点;
将从所述起始点开始在所述球管沿机架旋转方向经过的交点作为结束点;
计算所述起始点到所述结束点所转过的圆心角与180°的差值,即为所述给定放射起始角α0所对应的放射附加角δ的值;
根据各给定放射起始角α0及其对应的放射附加角δ,得到放射起始角α0与放射附加角δ之间的关系曲线;
或者,所述根据当前断层中感兴趣区的几何信息,计算放射起始角α0与放射附加角δ之间的关系曲线包括:
在0°-360°范围内逐一给定一放射起始角;
对每个给定的放射起始角α0,根据如下公式计算放射附加角δ:
δ = 2 arcsin ( ROI _ offset · sin ( ROI _ angle - α 0 ) R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - arcsin ( ROI _ radius R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) )
其中,ROI_adius是当前断层中感兴趣区的半径,ROI_offset是当前断层中感兴趣区相对于扫描中心的偏移量,ROI_angle是当前断层中感兴趣区的位置角,Rfocus是X射线球管运行圆周轨迹的半径;
根据各给定放射起始角α0及其对应的放射附加角δ,得到放射起始角α0与放射附加角δ之间的关系曲线。
本发明提出的X射线计算机断层扫描系统,包括:
一个粗略边界确定单元,用于根据预扫描出的定位图像,确定待检对象中三维感兴趣区的粗略边界;
一个三维感兴趣区计算单元,用于根据所述粗略边界,计算得到当前三维感兴趣区;
一个扫描方向确定单元,用于确定扫描方向;
一个扫描器,用于对所述扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像;对除第一个断层之外的当前断层中的感兴趣区,根据所述当前断层中的最优放射起始角α0opt和最小扫描角度γmin对所述当前断层中的感兴趣区进行部分扫描,得到当前断层的部分扫描图像;
一个精确边界获取单元,用于从所述完全扫描图像中获取第一个断层中感兴趣区的精确边界信息;在三维感兴趣区未扫描完成时,从所述部分扫描图像中获取当前断层中感兴趣区的精确边界信息;
一个三维感兴趣区校正单元,用于利用所述第一个断层中感兴趣区的精确边界信息对当前三维感兴趣区中第一个断层中的感兴趣区的几何信息进行校正;利用所述当前断层中感兴趣区的精确边界信息对当前三维感兴趣区中当前断层中的感兴趣区的几何信息进行校正;
一个断层几何信息确定单元,用于基于校正后的当前三维感兴趣区通过插值算法估算得到下一个断层中感兴趣区的几何信息;
一个断层扫描参数确定单元,用于将下一个断层作为当前断层,根据所述当前断层中感兴趣区的几何信息,确定当前断层中扫描所需的最优放射起始角α0opt和最小扫描角度γmin。
在本发明的一个实施方式中,所述粗略边界确定单元根据预扫描出的定位图像,确定待检对象中三维感兴趣区的由可识别边界和宽泛估计边界构成的粗略边界。
在本发明的一个实施方式中,所述扫描方向确定单元将三维感兴趣区中对应所述可识别边界的一侧作为所述三维感兴趣区的扫描起始侧。
在本发明的一个实施方式中,该系统进一步包括:一个感兴趣区数据库和一个粗略边界修正单元;
所述粗略边界修正单元,用于检索所述感兴趣区数据库,在所述感兴趣区数据库中存在所述待检对象的所述三维感兴趣区的几何信息时,利用所述三维感兴趣区的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界;
所述三维感兴趣区计算单元根据所述粗略边界修正单元修正后的粗略边界,计算得到当前三维感兴趣区。
在本发明的一个实施方式中,所述粗略边界修正单元进一步在所述感兴趣区数据库中不存在所述待检对象的所述三维感兴趣区的几何信息时,利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维感兴趣区的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。
在本发明的一个实施方式中,所述断层扫描参数确定单元包括:
一个关系曲线计算单元,用于根据当前断层中感兴趣区的几何信息,计算放射起始角α0与放射附加角δ之间的关系曲线;
一个起始角和附加角确定单元,用于根据所述关系曲线,确定当前断层中的最优放射起始角α0opt及其对应的最小放射附加角δmin;
一个最小扫描角度计算单元,用于根据所述最小放射附加角δmin,计算当前断层中的最小扫描角度γmin。
在本发明的一个实施方式中,所述关系曲线计算单元以如下方式逐一得到所述关系曲线上给定放射起始角α0所对应的放射附加角δ,继而得出所述关系曲线:
从当前断层中球管运行圆周轨迹上角度为给定放射起始角α0的起始点向XY扫描平面上的所述感兴趣区引两条外切线,并与所述球管运行圆周轨迹相交于两个交点;
将从所述起始点开始在所述球管沿机架旋转方向经过的交点作为结束点;
计算所述起始点到所述结束点所转过的圆心角与180°的差值,即为所述给定放射起始角α0所对应的放射附加角δ的值;
或者,根据当前断层中的感兴趣区和给定的放射起始角α0,所述关系曲线计算单元根据如下公式计算放射附加角δ,继而得出所述关系曲线:
δ = 2 arcsin ( ROI _ offset · sin ( ROI _ angle - α 0 ) R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - arcsin ( ROI _ radius R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) )
其中,ROI_adius是该二维断层中所述感兴趣区的半径,ROI_offset是该二维断层中所述感兴趣区相对于扫描中心的偏移量,ROI_angle是该二维断层中所述感兴趣区的位置角,Rfocus是所述球管运行圆周轨迹的半径。
在本发明的一个实施方式中,该系统进一步包括:一个处理器和一个计算机可读介质;
所述粗略边界确定单元、初始三维感兴趣区计算单元、扫描方向确定单元、精确边界获取单元、三维感兴趣区校正单元、断层几何信息确定单元和断层扫描参数确定单元中的任一个或任意组合为存储在所述计算机可读介质中的程序代码;
所述处理器,用于执行所述计算机可读介质中的所述程序代码。
从上述方案中可以看出,由于本发明实施例中依次利用每个断层的扫描图像校正该断层的ROI区域几何信息,并根据校正后的三维ROI区域通过插值算法估算下一个断层中ROI区域的几何信息,从而可以得到每个断层中较为精确的ROI区域,保证了扫描成像的质量;并且根据不同三维ROI区域,可以调整每个XY扫描平面上的X射线的放射起始角和扫描角度,这样在机架(gantry)的整个扫描过程中,X射线球管只是在每个XY扫描平面中较小的扫描角度范围内处于开启状态。因此,通过本发明的技术方案,可以利用每个断层中较小的扫描角度,即较少的放射剂量实现三维ROI的图像重建,节约了能耗,提高了图像重建的时间分辨率(temporal resolution),而且延长了X射线球管的寿命。
而且根据本发明的特点,当对于偏移扫描中心较多的器官或组织进行成像时,或者对于介入性治疗等需要对特定位置高质量重复成像的特殊医疗处理,可以得到更为显著的效果。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为本发明实施例中的X射线计算机断层扫描方法的示例性流程图。
图2a和图2b为本发明实施例中的CT设备预扫描出的定位图像;图2c中示出了本实施例所用实例中某下一个断层中ROI的几何信息的示意图。
图3a和图3b为本发明实施例中列举的根据定位图像确定三维ROI的粗略边界的一个实例的示意图。
图4是本发明具体实施方式中x-y扫描平面的放射起始角和放射附加角的几何关系示意图。
图5a至图5c为本发明一个具体实施例中x-y扫描平面中ROI断层的几何示意图和放射起始角α0与放射附加角δ的关系曲线图。
图6a至图6c为本发明另一个具体实施例中x-y扫描平面中ROI断层的几何示意图和放射起始角α0与放射附加角δ的关系曲线图。
图7为本发明实施例中X射线计算机断层扫描系统的示例性结构图。
图8为图7所示系统中断层扫描参数确定单元的一个内部结构示意图。
其中,附图标记如下:
101-根据预扫描出的定位图像,确定待检对象中三维ROI的粗略边界
102-判断ROI数据库中是否存在所述待检对象的所述三维ROI的几何信息
103-对所述粗略边界进行配准修正
104-对所述粗略边界进行配准修正
105-根据所述粗略边界,计算得到当前三维ROI
106-确定扫描方向
107-对所述扫描方向上的第一个断层采用完全扫描法进行扫描
108-将第一个断层作为当前断层,获取当前断层中ROI的精确边界信息
109-校正当前断层中的ROI的几何信息,并估算下一个断层中ROI的几何信息
110-将下一个断层作为当前断层,确定所需的最优放射起始角和最小扫描角度
111-对所述当前断层中的ROI进行部分扫描,得到当前断层的部分扫描图像
112-判断整个三维ROI是否扫描完成
113-从所述部分扫描图像中获取当前断层中ROI域的精确边界信息
114-结束
c1-可识别边界                    c2-宽泛估计边界
c3-修正边界                      h-垂直长度
c11隔膜顶端                      c12肋骨
ROI-感兴趣区                     FOM-测量域
701-粗略边界确定单元             702-三维感兴趣区计算单元
703-扫描方向确定单元             704-扫描器
705-精确边界获取单元             706-三维感兴趣区校正单元
707-断层几何信息确定单元         708-断层扫描参数确定单元
709-感兴趣区数据库               710-粗略边界修正单元
801-关系曲线计算单元             802-起始角和附加角确定单元
803-最小扫描角度计算单元
具体实施方式
在本发明中,定义病床的竖直升降方向为y方向,病床进出机架的水平方向为z方向,z方向与y方向垂直,而与y方向和z方向都正交的方向则为x方向。另外,x-y平面就是机架旋转平面,z方向就是垂直于机架旋转平面的方向。
在本发明实施例中,假设只需要针对三维ROI所对应的器官或者组织进行高质量的图像重建。该ROI区域可以包括部分器官或组织、整个器官或组织,或者医疗检查中所关注的几个器官或组织的集合。
本发明实施例中,为了在保证三维ROI的重建图像质量的前提下显著地减小每个断层中x-y平面内的扫描角度,考虑在根据预扫描出的定位图像(topography image)确定三维ROI后,进一步根据每个断层的扫描图像(tomography image)对所述三维ROI进行校正,并基于校正后的三维ROI估算下一个断层中ROI区域的几何信息,从而可以得到每个断层中较为精确的ROI的几何信息,进而可根据该断层中的ROI的几何信息,确定用于该断层扫描的x-y平面内的最优扫描起始角和最小扫描角度,从而可以利用较少的放射剂量实现三维ROI的图像重建。
此外,为了进一步保证三维ROI的重建图像质量,本发明实施例中可建立一个三维ROI数据库,并可首先根据预扫描出的定位图像确定三维ROI的粗略边界,之后利用该数据库中的相应三维ROI的信息对所述粗略边界进行配准修正,以得到修正后的粗略边界,之后再根据该修正后的粗略边界确定三维ROI,进而再如上所述,根据每个断层的扫描图像对所述三维ROI进行校正。
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
图1为本发明实施例中一种X射线计算机断层扫描方法的示例性流程图。如图1所示,该方法包括如下步骤:
步骤101,根据预扫描出的定位图像,确定待检对象中三维ROI的粗略边界。
在CT设备开始正式扫描前,通常需要对待检对象预扫描若干幅定位图像(topographyimage)。图2a和图2b示出了本发明实施例中所使用的定位图像的实例。通常可以对待检对象扫描如图2a所示的正位定位图像和如图2b所示的侧位定位图像等。本步骤中,由于定位图像相对断层扫描图像(tomography image)来说是一个较为粗略的图,基于不同的三维ROI的特点及对比度,通常只能从定位图像中确定待检对象三维ROI大概的位置和大小,从而得到粗略的扫描范围,即三维ROI的粗略边界。该粗略边界有可能都能从定位图像中识别出,即此时该粗略边界完全由可识别边界构成,也有可能只能从定位图像中识别出一部分,而另一部分则只能设置一个较大的力求把所有可能边界包括在内的宽泛估计边界,即此时该粗略边界由可识别边界和宽泛估计边界构成。如图2a所示,在所述正位定位图像中只能得到由可识别边界c1和宽泛估计边界c2构成的粗略边界。
此外,图3a和图3b为本发明实施例中列举的根据定位图像确定三维ROI的粗略边界的一个实例。本实例中以三维ROI为肝脏的情况为例,可以看出,肝脏的一部分边界可以根据图3a和图3b中的结构信息确定,即根据对比信息可以确定隔膜顶端c11和肋骨c12的位置,由于右肝叶的顶端从隔膜的位置开始,左肝叶的顶端靠近心脏,此外肝脏总是位于肋骨内侧,因此可以在图3a中确定出可识别边界c1。肝脏的另一部分边界由于缺少足够的对比信息,因此只能根据一定的垂直长度h得到一个宽泛估计边界c2。其中,该垂直长度h可根据人种以及待检对象健康与否等状态确定。肝脏的正常长度大概是20~30厘米,但患病的肝脏可能会更大更长,因此可根据充足的统计数据确定肝脏的垂直长度安全偏移量,用于确定肝脏的宽泛边界。图3b中的肝脏粗略边界的确定方法与之类似,此处不再赘述。
步骤102,检索一ROI数据库,判断所述ROI数据库中是否存在所述待检对象的所述三维ROI的几何信息,若存在,则执行步骤103;否则执行步骤104。
本实施例中,可建立一ROI数据库,用于存储大量检查对象的各种三维ROI的几何信息,该三维ROI几何信息可以是从断层扫描图像中获取的,也可以是基于任何其他图像处理方法得到的。
在所述ROI数据库中不存在所述待检对象的所述三维ROI的几何信息时,也可不执行步骤104,而是直接执行步骤105。
步骤103,利用所述三维ROI的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。其中,配准方法通常可以采用图像配准的常用方法,此处不再详述。修正方法也可有多种,如可将所述三维ROI的几何信息和所述粗略边界进行加权平均等。
本实施例中,对图2a和图2b所示由可识别边界c1和宽泛估计边界c2构成的粗略边界进行修正后,可得到如图2a中所示的由可识别边界c1和修正边界c3构成的粗略边界。
步骤104,利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维ROI的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。其中,配准方法同样可以采用图像配准常用方法,此处不再详述。修正方法也可有多种,如可将所述三维ROI的几何信息和所述粗略边界进行加权平均等。
步骤105,根据所述粗略边界,计算得到当前三维ROI。
对于某具体医疗检查应用,三维ROI可以根据该粗略边界,运用特定的图像识别算法来自动识别出当前三维ROI区域。其中,相关图像识别算法可以采用常用图像识别算法,此处不再详述。
为了便于分析和理解,我们可以将三维ROI理解为若干个断层在z方向上堆叠而成的三维空间体,其中每个断层构成一个x-y扫描平面,又可称为二维x-y平面ROI或ROI断层。由此,本步骤中也可以根据该粗略边界计算各二维断层中的ROI的几何参数,从而确定出三维ROI。
对于z轴上的每个位置(用变量zi表示),图2(a)中所示的待检对象正位定位图像可用来确定三维ROI在x方向上的大小和位置,该大小和位置用变量Lx(zi)和Cx(zi)表示,其中Lx(zi)和Cx(zi)分别表示正位定位图像中在zi位置处线断层的长度和中心位置;同时,图2(b)中所示的待检图像侧位定位图像可用来确定三维ROI在y方向上的大小和位置,该大小和位置用变量Ly(zi)和Cy(zi)表示,其中Ly(zi)和Cy(zi)分别表示侧位定位图像中在zi位置处线断层的长度和中心位置。
为了方便分析和计算,可将三维ROI区域在每个x-y平面上的断层均指定为一圆形。这样,对于z轴上的每个zi位置,有关ROI某一个x-y平面断层的几个主要几何参数就可以根据定位图像中粗略边界的变量参数,通过如下公式(1)-(4)计算得到:
ROI _ radius | z i = max { L x ( z i ) 2 , L y ( z i ) 2 } - - - ( 1 )
ROI _ offset | z i = C x ( z i ) 2 + C y ( z i ) 2 - - - ( 2 )
cos ( ROI _ angle ) = C x ( z i ) C x ( z i ) 2 + C y ( z i ) 2 - - - ( 3 )
sin ( ROI _ angle ) = C y ( z i ) C x ( z i ) 2 + C y ( z i ) 2 - - - ( 4 )
其中,ROI_radius是ROI断层的圆半径;ROI_offset是ROI断层相对于相应x-y平面的扫描中心(ISO中心)的偏移量;ROI_angle是ROI的位置角,其在0°-360°之间取值。本领域技术人员可知,通过公式(3)和公式(4)的结果,可以确定出ROI_angle的值。
这样,根据ROI在各x-y断层上的几何参数(ROI_radius,ROI_offset和ROI_angle的值),可以构建出当前三维ROI的立体模型。
步骤106,确定扫描方向。
本步骤中,可将三维ROI中对应所述可识别边界的一侧作为所述三维ROI的扫描起始侧。
对应图2a和图2b所示实例以及图3a和图3b所示实例,其扫描方向可以是从头到脚。
步骤107,对所述扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像。
本步骤中的完全扫描方法与现有技术中的常规完全扫描方法一样,此处不再赘述。
步骤108,将第一个断层作为当前断层,从所述完全扫描图像中获取当前断层中ROI的精确边界信息。
步骤109,利用所述精确边界信息对当前三维ROI中当前断层中的ROI的几何信息进行校正,并基于校正后的当前三维ROI估算得到下一个断层中ROI的几何信息。
本步骤中,可根据校正后的当前断层中的ROI的几何信息,利用插值算法,如外插值法,估算得到下一个断层中ROI的几何信息。
图2c中示出了本实施例所用实例中某下一个断层中ROI的几何信息的示意图。
步骤110,将下一个断层作为当前断层,根据所述当前断层中ROI的几何信息,确定当前断层中扫描所需的最优放射起始角α0opt和最小扫描角度γmin。
对于每个断层,本实施例中可以确定使所需扫描角度γ最小的最优放射起始角α0opt以及所对应的最小扫描角度γmin。
经过本发明的发明人研究发现,在只需要针对ROI在x-y扫描平面上的断层进行高质量的图像重建的情况下,从适当的放射起始角α0开始进行一定扫描角度的扫描即可得到对ROI进行数据重建所需的最小数据量,扫描角度γ可按照如下公式(5)来计算:
γ=180°+δ                    ……………………(5)
其中δ在本发明中称为放射附加角,其是由ROI断层和放射起始角α0共同确定的角度值。
通过研究可知,由于放射附加角δ通常小于现有技术中常规所设定的60°扇角或
Figure BDA0000149419260000111
甚至可能是负角度,因此相比起现有技术,本发明可以显著地减小扫描各ROI断层所需的扫描角度γ,从而减小总的放射剂量。
本步骤中,可首先根据当前断层中ROI的几何关系,计算出放射起始角α0与放射附加角δ之间的关系曲线。
放射起始角α0与放射附加角δ之间的关系曲线至少有如下两种实现方式。
方式一:通过几何方式来确定某给定放射起始角α0所对应的放射附加角δ的值。
图4是本发明实施例中x-y扫描平面中放射起始角和放射附加角的几何关系示意图。
如图4所示,当X射线球管从球管运行圆周轨迹(如图4中外圈虚线所示)上的位置1作为放射起始点开始发射X射线并沿圆周轨迹进行顺时针转动时,当设定以从ISO中心4指向3点钟方向的坐标轴作为0°坐标轴位置,则该位置1在x-y平面中所对应的转角即为放射起始角α0。图4中同时示出了在此x-y扫描平面中确定出的、需要进行高质量图像重建的ROI断层的位置。
如图4所示,从位置1引出ROI断层的圆形区域的两条外切线,并与球管运行圆周轨迹相交于位置2和位置3两点。经过研究发现,将X射线球管沿球管沿机架旋转方向相遇的交点位置,将该点作为放射结束点,在图4中即为位置2的点,即可得到对ROI断层进行数据重建所需的最小数据量。由此可见,放射起始点到放射结束点所转过的圆心角与180°的差值,即为给定放射起始角α0所对应的放射附加角δ的值。
由于在图4所示的实例中δ角度为负值,因此图4中位置2与ISO中心4的连线和位置1与ISO中心4连线的延长线之间的夹角大小即为放射附加角δ的绝对值大小。
此外,从图4可以看出,X射线球管从放射起始点(位置1)到放射结束点(位置2)所旋转过的角度即为扫描角度γ,根据公式(5),放射附加角δ与扫描角度γ具有180°的角度差,因而可以通过确定放射附加角δ来确定扫描角度γ。由此关系可知,δ值越小则X射线球管在该x-y扫描平面中的扫描角度γ越小,则对于该ROI断层实施的放射剂量越小。
方式二:通过数学公式来确定某给定放射起始角α0所对应的放射附加角δ的值。
对于某给定的放射起始角α0,本领域技术人员也可以根据在公式(1)-(4)中所确定的ROI断层的数学参数,通过数学公式来计算出给定放射起始角α0所对应的放射附加角δ的值。
对于在x-y平面中确定的ROI断层,可以通过公式(1)-(4)中所确定的ROI_radius、ROI_offset和ROI_angle三个变量的值来确定该ROI的准确位置。本领域技术人员可以理解,根据上述三个变量所确定的ROI位置和给定的放射起始角α0,便可以根据如下公式(6)计算出放射附加角δ的值:
δ = 2 arcsin ( ROI _ offset · sin ( ROI _ angle - α 0 ) R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - arcsin ( ROI _ radius R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - - - ( 6 )
其中,Rfocus是球管运行圆周轨迹的半径,即图4中外圈虚线圆的半径。
如上所述,本领域技术人员可以用图4中的几何作图法或者根据所给定的ROI断层的数学参数来确定给定放射起始角α0所对应的放射附加角δ。当放射起始角在α0在0°-360°范围内逐一取值时,就可以得到在给定ROI断层的条件下,放射起始角α0与放射附加角δ的关系曲线,如图5a和图6a所示。
图5a至图5c和图6a至图6c分别示出了两个不同x-y扫描平面中ROI断层的几何示意图和放射起始角α0与放射附加角δ的关系曲线图。在图5b、图5c、图6b、图6c中,分别示出了对于两个不同ROI断层的x-y平面几何示意图,其中图5a至图5c中的ROI断层较大,且覆盖了x-y扫描平面的ISO中心,而图6a至图6c中的ROI断层较小,且远离x-y扫描平面的ISO中心。图5a和图6a分别针对两个ROI断层,示出了放射起始角α0与放射附加角δ的关系曲线。本领域技术人员可以用图4中的几何作图法或者根据所给定的ROI断层的数学参数,通过公式(6)来得到图5a和图6a所示的关系曲线。
通过公式(5)可知,放射附加角δ越小则X射线球管的扫描角度γ越小,则针对ROI断层所需的放射剂量越小,因而本领域技术人员需要确定出x-y扫描平面中使放射附加角δ达到最小值δmin的最优放射起始角α0opt。
之后,根据计算出的放射起始角α0与放射附加角δ之间的关系曲线,本步骤中可以将关系曲线上放射附加角δ为最小值时的放射起始角α0确定为的最优放射起始角α0opt。同时也可以确定出最小放射附加角δmin的值。进而,可以根据所确定的最小放射附加角δmin,利用公式(5)求得最小扫描角度γmin。
步骤111,根据所述最优放射起始角α0opt和最小扫描角度γmin对所述当前断层中的ROI进行部分扫描,得到当前断层的部分扫描图像。
下面,我们以图5a至图5c和图6a至图6c为例,说明本具体实施例中的X射线计算机断层扫描系统在对x-y扫描平面中的ROI断层进行扫描时,相对于现有技术的显著技术效果。
如图5a至图5c所示,这里所给定的ROI断层较大而且覆盖x-y扫描平面的ISO中心,它的位置参数分别是ROI_radius=150mm、ROI_offset=50mm、ROI_angle=60°。此时,最优放射起始角α0opt所对应的最小放射附加角δmin=19.35°,即此时的最小扫描角度γmin为180°+19.35°=199.35°。对比传统的完全扫描法,扫描角度γ和所对应的放射剂量减少了大约44.6%,而对比部分扫描法,扫描角度γ和所对应的放射剂量减少了大约17%。
如图6a至图6c所示,这里所给定的ROI断层较小而且远离x-y扫描平面的ISO中心,它的位置参数分别是ROI_radius=50mm、ROI_offset=200mm、ROI_angle=230°。此时,最优放射起始角α0opt所对应的最小放射附加角δmin=-29.2°,即此时的最小扫描角度γmin为180°+(-29.2°)=154.7°。对比传统的完全扫描法,扫描角度γ和所对应的放射剂量减少了大约57%,而对比部分扫描法,扫描角度γ和所对应的放射剂量减少了大约37.2%。
由图5a至图5c和图6a至图6c中不同的ROI断层情况可以看出,当ROI断层并不覆盖x-y扫描平面的ISO中心的时候,δmin为负值,此时扫描角度γ可以小于180°,因此x-y扫描平面中用于扫描ROI断层的放射剂量可以得到显著的减小。特别是,当ROI断层的范围较小而且靠近x-y扫描平面的FOM边缘的时候,放射剂量可以得到更大幅度的减小。
步骤112,判断整个三维ROI是否扫描完成,若未完成,则执行步骤113;否则,在步骤114结束当前扫描。
步骤113,从所述部分扫描图像中获取当前断层中ROI的精确边界信息;并返回执行步骤109。
也可不执行上述步骤102至步骤104,则步骤105中可直接根据步骤101中的粗略边界,计算得到当前三维ROI。
以上对本发明实施例中的X射线计算机断层扫描方法进行了详细的描述,下面再对本发明实施例中基于所述X射线计算机断层扫描方法实现的X射线计算机断层扫描系统进行详细描述。
图7为本发明实施例中X射线计算机断层扫描系统的示例性结构图。如图7所示,该系统包括:一个粗略边界确定单元701、一个三维感兴趣区计算单元702、一个扫描方向确定单元703、一个扫描器704、一个精确边界获取单元705、一个三维感兴趣区校正单元706、一个断层几何信息确定单元707和一个断层扫描参数确定单元708。
粗略边界确定单元701,用于根据预扫描出的定位图像,确定待检对象中三维ROI的粗略边界。
三维感兴趣区计算单元702,用于根据所述粗略边界,计算得到当前三维ROI。
扫描方向确定单元703,用于确定扫描方向。
扫描器704,用于对所述扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像;对除第一个断层之外的当前断层中的ROI,根据所述当前断层中的最优放射起始角α0opt和最小扫描角度γmin对所述当前断层中的ROI进行部分扫描,得到当前断层的部分扫描图像。
本实施例中,将X射线计算机断层扫描系统中实现X射线扫描的必要器件一起称为扫描器704。扫描器704可包括X射线计算机断层扫描系统的旋转机架、对应安装在旋转机架上的X射线球管和X射线探测器,还包括控制旋转机架2带动X射线球管3和X射线探测器4旋转并控制X射线球管3和X射线探测器4进行X射线扫描及采集的控制装置等。
精确边界获取单元705,用于从所述完全扫描图像中获取第一个断层中ROI的精确边界信息;在三维ROI未扫描完成时,从所述部分扫描图像中获取当前断层中ROI的精确边界信息。
三维感兴趣区校正单元706,用于利用所述第一个断层中ROI的精确边界信息对当前三维ROI中第一个断层中的ROI的几何信息进行校正;利用所述当前断层中ROI的精确边界信息对当前三维ROI中当前断层中的ROI的几何信息进行校正。
断层几何信息确定单元707,用于基于校正后的当前三维ROI通过插值算法(如外插值法)估算得到下一个断层中ROI的几何信息。
断层扫描参数确定单元708,用于将下一个断层作为当前断层,根据所述当前断层中ROI的几何信息,确定当前断层中扫描所需的最优放射起始角α0opt和最小扫描角度γmin。
与本发明实施例中的X射线计算机断层扫描方法一致,所述粗略边界确定单元701根据预扫描出的定位图像,确定出的待检对象中三维ROI的粗略边界,可以是完全由可识别边界构成的粗略边界,也可以是由可识别边界和宽泛估计边界构成的粗略边界。
同样,在确定扫描方向时,所述扫描方向确定单元703可将三维ROI中对应所述可识别边界的一侧作为所述三维ROI的扫描起始侧。
本发明实施例中的X射线计算机断层扫描系统可进一步包括:一个感兴趣区数据库709和一个粗略边界修正单元710。
所述粗略边界修正单元710,用于检索所述感兴趣区数据库709,在所述感兴趣区数据库709中存在所述待检对象的所述三维ROI的几何信息时,利用所述三维ROI的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。此时,所述三维感兴趣区计算单元702可根据所述粗略边界修正单元710修正后的粗略边界,计算得到当前三维ROI。
同样,所述粗略边界修正单元710可进一步在所述感兴趣区数据库709中不存在所述待检对象的所述三维ROI的几何信息时,可利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维ROI的几何信息对所述粗略边界进行配准修正,得到修正后的粗略边界。
与本发明实施例中的X射线计算机断层扫描方法中的描述类似,所述断层扫描参数确定单元708可如图8所示,包括:一个关系曲线计算单元801、一个起始角和附加角确定单元802和一个最小扫描角度计算单元803。
关系曲线计算单元801,用于根据当前断层中ROI的几何信息,计算放射起始角α0与放射附加角δ之间的关系曲线。
起始角和附加角确定单元802,用于根据所述关系曲线,确定当前断层中的最优放射起始角α0opt及其对应的最小放射附加角δmin。
最小扫描角度计算单元803,用于根据所述最小放射附加角δmin,计算当前断层中的最小扫描角度γmin。
与本发明实施例中的X射线计算机断层扫描方法中的描述相似,所述关系曲线计算单元801可以如下方式逐一得到所述关系曲线上给定放射起始角α0所对应的放射附加角δ,继而得出所述关系曲线:
从当前断层中球管运行圆周轨迹上角度为给定放射起始角α0的起始点向XY扫描平面上的所述ROI引两条外切线,并与所述球管运行圆周轨迹相交于两个交点;
将从所述起始点开始在所述球管沿机架旋转方向经过的交点作为结束点;
计算所述起始点到所述结束点所转过的圆心角与180°的差值,即为所述给定放射起始角α0所对应的放射附加角δ的值;
或者,所述关系曲线计算单元801也可根据当前断层中的ROI和给定的放射起始角α0,所述关系曲线计算单元根据如下公式计算放射附加角δ,继而得出所述关系曲线:
δ = 2 arcsin ( ROI _ offset · sin ( ROI _ angle - α 0 ) R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - arcsin ( ROI _ radius R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) )
其中,ROI_radius是该二维断层中所述ROI的半径,ROI_offset是该二维断层中所述ROI相对于扫描中心的偏移量,ROI_angle是该二维断层中所述ROI的位置角,Rfocus是所述球管运行圆周轨迹的半径。
该系统可进一步包括:一个处理器和一个计算机可读介质(图中未示出)。此时,所述粗略边界确定单元701、初始三维感兴趣区计算单元702、扫描方向确定单元703、精确边界获取单元705、三维感兴趣区校正单元706、断层几何信息确定单元707和断层扫描参数确定单元708中的任一个或任意组合为存储在所述计算机可读介质中的程序代码。
所述处理器,用于执行所述计算机可读介质中的所述程序代码。
本发明公开了一种X射线计算机断层扫描扫描方法和一种X射线计算机断层扫描扫描系统。方法包括:确定三维感兴趣区的粗略边界,并计算得到当前三维感兴趣区;对第一个断层采用完全扫描法扫描,从得到的图像中获取该断层中感兴趣区的精确边界信息;利用该精确边界信息对当前三维感兴趣区中相应断层中的感兴趣区的几何信息进行校正;基于校正后的当前三维感兴趣区估算得到下一个断层中感兴趣区的几何信息,并确定该断层中扫描所需的最优放射起始角和最小扫描角度,基于此对该断层中的感兴趣区进行部分扫描,从得到的图像中获取该断层中感兴趣区的精确边界信息,并返回执行上述对当前三维感兴趣区进行校正的步骤。该方案能够在保证三维ROI的重建图像质量的前提下有效减少X射线的总放射剂量。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (12)

1.一种X射线计算机断层扫描方法,包括:
根据一定位图像,确定待检对象中三维感兴趣区的粗略边界,并根据所述粗略边界计算得到当前三维感兴趣区;
对一扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像,从所述完全扫描图像中获取第一个断层中感兴趣区的精确边界信息,利用所述精确边界信息对第一个断层中感兴趣区的几何信息进行校正,并基于校正后的当前三维感兴趣区估算得到下一个断层中感兴趣区的几何信息;
将下一个断层作为当前断层,根据当前断层中感兴趣区的几何信息,确定当前断层中扫描所需的最优放射起始角和最小扫描角度,根据所述最优放射起始角和最小扫描角度对所述当前断层中的感兴趣区进行部分扫描,得到当前断层的部分扫描图像。
2.根据权利要求1所述的方法,其特征在于,所述粗略边界由可识别边界和宽泛估计边界构成,或者全部由可识别边界构成。
3.根据权利要求2所述的方法,其特征在于,该方法进一步包括:将三维感兴趣区中对应所述可识别边界的一侧作为扫描起始侧。
4.根据权利要求2所述的方法,其特征在于,该方法在计算当前三维感兴趣区之前进一步包括:检索一感兴趣区数据库,在所述感兴趣区数据库中存在所述待检对象的所述三维感兴趣区的几何信息时,利用所述三维感兴趣区的几何信息对所述粗略边界进行配准修正;优选地,在所述感兴趣区数据库中不存在所述待检对象的所述三维感兴趣区的几何信息时,利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维感兴趣区的几何信息对所述粗略边界进行配准修正。
5.根据权利要求1至4中任一项所述的方法,其特征在于,所述确定最优放射起始角和最小扫描角度的步骤包括:
根据当前断层中感兴趣区的几何信息,计算放射起始角与放射附加角之间的关系曲线;
根据所述关系曲线,确定当前断层中的最优放射起始角及其对应的最小放射附加角;
根据所述最小放射附加角,计算当前断层中的最小扫描角度。
6.根据权利要求5所述的方法,其特征在于,所述计算放射起始角与放射附加角之间的关系曲线包括:
在0°-360°范围内逐一给定一放射起始角;对每个给定的放射起始角,从所述断层中X射线球管运行圆周轨迹上为给定放射起始角的起始点向机架旋转平面上的所述感兴趣区引两条外切线,并与所述球管运行圆周轨迹相交于两个交点;将从所述起始点开始在所述球管沿机架旋转方向经过的交点作为结束点;计算所述起始点到所述结束点所转过的圆心角与180°的差值,即为所述给定放射起始角所对应的放射附加角的值;根据各给定放射起始角及其对应的放射附加角,得到放射起始角与放射附加角之间的关系曲线;
或者,所述计算放射起始角与放射附加角之间的关系曲线包括:
在0°-360°范围内逐一给定一放射起始角;对每个给定的放射起始角,根据如下公式计算放射附加角:
δ = 2 arcsin ( ROI _ offset · sin ( ROI _ angle - α 0 ) R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) ) - arcsin ( ROI _ radius R focus 2 + ROI _ offset 2 - 2 R focus · ROI _ offset · cos ( ROI _ angle - α 0 ) )
其中,ROI_radius是当前断层中感兴趣区的半径,ROI_offset是当前断层中感兴趣区相对于扫描中心的偏移量,ROI_angle是当前断层中感兴趣区的位置角,Rfocus是X射线球管运行圆周轨迹的半径,α0为放射起始角,δ为放射附加角;
根据各给定放射起始角及其对应的放射附加角,得到放射起始角与放射附加角之间的关系曲线。
7.根据权利要求5所述的方法,其特征在于,该方法进一步包括:从所述当前断层的部分扫描图像中获取当前断层中感兴趣区的精确边界信息,利用所述精确边界信息对当前断层中感兴趣区的几何信息进行校正,并基于校正后的当前三维感兴趣区通过估算得到下一个断层中感兴趣区的几何信息。
8.一种X射线计算机断层扫描系统,包括:
一个粗略边界确定单元(701),用于根据一定位图像,确定待检对象中三维感兴趣区的粗略边界;
一个三维感兴趣区计算单元(702),用于根据所述粗略边界,计算得到当前三维感兴趣区;
一个扫描器(704),用于对一扫描方向上的第一个断层采用完全扫描法进行扫描,得到第一个断层的完全扫描图像;对除第一个断层之外的当前断层中的感兴趣区,根据所述当前断层中的最优放射起始角和最小扫描角度进行部分扫描,得到当前断层的部分扫描图像;
一个精确边界获取单元(705),用于从所述完全扫描图像中获取第一个断层中感兴趣区的精确边界信息,及从所述部分扫描图像中获取当前断层中感兴趣区的精确边界信息;
一个三维感兴趣区校正单元(706),用于利用所述第一个断层中感兴趣区的精确边界信息对当前三维感兴趣区中第一个断层中的感兴趣区的几何信息进行校正;利用所述当前断层中感兴趣区的精确边界信息对当前三维感兴趣区中当前断层中的感兴趣区的几何信息进行校正;
一个断层几何信息确定单元(707),用于基于校正后的当前三维感兴趣区估算得到下一个断层中感兴趣区的几何信息;
一个断层扫描参数确定单元(708),用于将下一个断层作为当前断层,根据所述当前断层中感兴趣区的几何信息,确定当前断层中扫描所需的最优放射起始角和最小扫描角度。
9.根据权利要求8所述的系统,其特征在于,所述粗略边界确定单元(701)根据所述定位图像,确定待检对象中三维感兴趣区的由可识别边界和宽泛估计边界构成、或者全部由可识别边界构成的粗略边界。
10.根据权利要求9所述的系统,其特征在于,进一步包括一个扫描方向确定单元(703),用于将三维感兴趣区中对应所述可识别边界的一侧作为所述三维感兴趣区的扫描起始侧。
11.根据权利要求9所述的系统,其特征在于,该系统进一步包括:一个感兴趣区数据库(709)和一个粗略边界修正单元(710);
其中,所述粗略边界修正单元(710),用于检索所述感兴趣区数据库(709),在所述感兴趣区数据库(709)中存在所述待检对象的所述三维感兴趣区的几何信息时,利用所述三维感兴趣区的几何信息对所述粗略边界进行配准修正;优选地,所述粗略边界修正单元(710)进一步在所述感兴趣区数据库(709)中不存在所述待检对象的所述三维感兴趣区的几何信息时,利用相同人种、性别中相似体重、年龄和疾病的待检对象的相同三维感兴趣区的几何信息对所述粗略边界进行配准修正;
所述三维感兴趣区计算单元(702)根据所述粗略边界修正单元(710)修正后的粗略边界,计算得到当前三维感兴趣区。
12.根据权利要求8至11所述的系统,其特征在于,所述断层扫描参数确定单元(708)包括:
一个关系曲线计算单元(801),用于根据当前断层中感兴趣区的几何信息,计算放射起始角与放射附加角之间的关系曲线;
一个起始角和附加角确定单元(802),用于根据所述关系曲线,确定当前断层中的最优放射起始角及其对应的最小放射附加角;
一个最小扫描角度计算单元(803),用于根据所述最小放射附加角,计算当前断层中的最小扫描角度。
CN201210101500.5A 2012-03-31 2012-03-31 一种x射线计算机断层扫描方法和系统 Active CN103356218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210101500.5A CN103356218B (zh) 2012-03-31 2012-03-31 一种x射线计算机断层扫描方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210101500.5A CN103356218B (zh) 2012-03-31 2012-03-31 一种x射线计算机断层扫描方法和系统

Publications (2)

Publication Number Publication Date
CN103356218A true CN103356218A (zh) 2013-10-23
CN103356218B CN103356218B (zh) 2015-06-24

Family

ID=49359108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210101500.5A Active CN103356218B (zh) 2012-03-31 2012-03-31 一种x射线计算机断层扫描方法和系统

Country Status (1)

Country Link
CN (1) CN103356218B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104739431A (zh) * 2013-12-31 2015-07-01 上海西门子医疗器械有限公司 确定曝光范围的方法和计算机断层扫描设备
CN107468262A (zh) * 2016-06-07 2017-12-15 西门子保健有限责任公司 通过x光设备的身体部分的三维图像生成
CN107924565A (zh) * 2015-06-24 2018-04-17 光线搜索实验室公司 处理图像数据的系统和方法
CN111281430A (zh) * 2018-12-06 2020-06-16 深圳迈瑞生物医疗电子股份有限公司 超声成像方法、设备及可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1650807A (zh) * 2002-11-05 2005-08-10 西门子公司 对周期运动检查对象进行计算机断层造影的方法及设备
US20050276375A1 (en) * 2004-06-15 2005-12-15 Canon Kabushiki Kaisha X-ray CT apparatus and image-taking method therefor
CN101011258A (zh) * 2005-11-15 2007-08-08 Ge医疗系统环球技术有限公司 X射线ct设备和x射线ct透视设备
US20080118131A1 (en) * 2006-11-22 2008-05-22 General Electric Company Method and system for automatically identifying and displaying vessel plaque views
CN101987020A (zh) * 2009-08-04 2011-03-23 Ge医疗系统环球技术有限公司 倾斜图像扫描方法和重建方法及装置
CN102846331A (zh) * 2011-06-30 2013-01-02 上海西门子医疗器械有限公司 一种x射线计算机断层扫描系统和方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1650807A (zh) * 2002-11-05 2005-08-10 西门子公司 对周期运动检查对象进行计算机断层造影的方法及设备
US20050276375A1 (en) * 2004-06-15 2005-12-15 Canon Kabushiki Kaisha X-ray CT apparatus and image-taking method therefor
CN101011258A (zh) * 2005-11-15 2007-08-08 Ge医疗系统环球技术有限公司 X射线ct设备和x射线ct透视设备
US20080118131A1 (en) * 2006-11-22 2008-05-22 General Electric Company Method and system for automatically identifying and displaying vessel plaque views
CN101987020A (zh) * 2009-08-04 2011-03-23 Ge医疗系统环球技术有限公司 倾斜图像扫描方法和重建方法及装置
CN102846331A (zh) * 2011-06-30 2013-01-02 上海西门子医疗器械有限公司 一种x射线计算机断层扫描系统和方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104739431A (zh) * 2013-12-31 2015-07-01 上海西门子医疗器械有限公司 确定曝光范围的方法和计算机断层扫描设备
CN107924565A (zh) * 2015-06-24 2018-04-17 光线搜索实验室公司 处理图像数据的系统和方法
CN107468262A (zh) * 2016-06-07 2017-12-15 西门子保健有限责任公司 通过x光设备的身体部分的三维图像生成
US10245001B2 (en) 2016-06-07 2019-04-02 Siemens Healthcare Gmbh Generation of a three-dimensional reconstruction of a body part by an X-ray machine
CN107468262B (zh) * 2016-06-07 2019-06-21 西门子保健有限责任公司 通过x光设备的身体部分的三维图像生成和适用于其的x光设备
CN111281430A (zh) * 2018-12-06 2020-06-16 深圳迈瑞生物医疗电子股份有限公司 超声成像方法、设备及可读存储介质
CN111281430B (zh) * 2018-12-06 2024-02-23 深圳迈瑞生物医疗电子股份有限公司 超声成像方法、设备及可读存储介质

Also Published As

Publication number Publication date
CN103356218B (zh) 2015-06-24

Similar Documents

Publication Publication Date Title
US10614634B2 (en) System and method for image composition
US7983457B2 (en) Method and system for automatically determining regions in a scanned object
US9427286B2 (en) Method of image registration in a multi-source/single detector radiographic imaging system, and image acquisition apparatus
US8977026B2 (en) Methods and systems for locating a region of interest in an object
US20080037843A1 (en) Image segmentation for DRR generation and image registration
US20090226066A1 (en) Method and apparatus for aligning a multi-modality imaging system
US20070003020A1 (en) Systems and methods for compensating for table sag
EP1667580A2 (en) 2d/3d image registration in image-guided radiosurgery
Ruchala et al. Limited-data image registration for radiotherapy positioning and verification
US8199879B2 (en) Methods of scatter correction of x-ray projection data 1
CN103356218B (zh) 一种x射线计算机断层扫描方法和系统
US8953739B2 (en) Image processing apparatus and X-ray CT system
CN111166332B (zh) 一种基于磁共振波谱和磁共振影像的放疗靶区勾画方法
US20110110570A1 (en) Apparatus and methods for generating a planar image
CN102397078A (zh) 一种x射线计算机断层扫描系统和方法
Xia et al. Patient‐bounded extrapolation using low‐dose priors for volume‐of‐interest imaging in C‐arm CT
US20220054862A1 (en) Medical image processing device, storage medium, medical device, and treatment system
US20230169668A1 (en) Systems and methods for image registration
CN102846331B (zh) 一种x射线计算机断层扫描系统和方法
US20230083704A1 (en) Systems and methods for determining examination parameters
US9504437B2 (en) Diagnostic imaging apparatus and control method of the same
US20220148180A1 (en) Medical image processing device, storage medium, medical device, and treatment system
JP6821839B1 (ja) 二軸デジタルトモシンセシスシステムに用いられる幾何学的補正方法及びそのシステム
CN102846332B (zh) 一种x射线计算机断层扫描系统和方法
Robertson et al. A block matching‐based registration algorithm for localization of locally advanced lung tumors

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant