CN103163165A - 一种二代ct扫描成像方法 - Google Patents
一种二代ct扫描成像方法 Download PDFInfo
- Publication number
- CN103163165A CN103163165A CN2013100634021A CN201310063402A CN103163165A CN 103163165 A CN103163165 A CN 103163165A CN 2013100634021 A CN2013100634021 A CN 2013100634021A CN 201310063402 A CN201310063402 A CN 201310063402A CN 103163165 A CN103163165 A CN 103163165A
- Authority
- CN
- China
- Prior art keywords
- translation
- field
- scan
- scanning
- ray
- 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
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种二代CT扫描成像方法,它包括有射线源、线阵探测器和扫描转台,具体步骤为:1)初始化CT扫描参数;2)对工件进行平移扫描;3)扫描转台旋转一个角度,继续对工件进行平移扫描,直到扫描转台旋转完一周;4)对扫描得到的数据进行重排;5)采用SART迭代重建算法重建CT图像。本发明在保证CT图像质量不降低的情况下,比传统的二代扫描方式速度快,有效的提高了CT的扫描效率。
Description
技术领域
本发明涉及CT扫描技术领域,特别是一种二代CT的扫描成像方法。
背景技术
计算机断层成像(Computed Tomography,CT)扫描方式有:一代扫描方式、二代扫描方式(TR)、三代扫描方式(RO)和锥束扫描方式等。锥束扫描方式使用的是面阵探测器,锥束扫描方式具有采集数据量大,全身扫描速度快的优点,但是由于面阵探测器结构上的原因,锥束扫描方式有明显的缺点,即无法设计制造出有效的后准直器,不能限制X射线的散射和串扰,且动态范围小,尤其是可接受的射线能量低,一般只接受能量较低的X射线,对大尺寸、高密度的工件不适用。
二代扫描方式和三代扫描方式使用的是线阵探测器,线阵探测器单元一般比面阵探测器单元大,线阵探测器的动态范围大,可接受的射线能量高,不仅适用于X射线球管,而且适用于电子直线加速器;二代和三代扫描方式一般使用后准直器,有利于屏蔽散射射线。三代扫描方式只需要转台和扇束相对旋转,是一种常用的扫描方式,它的优点是扫描射线利用率高,扫描速度快,缺点是很难去掉环状伪影。二代扫描方式需要作平移运动和旋转运动,二代扫描方式的优点是不需要进行探测器一致性校正,由扫描的数据重建的CT图像无环形伪影,因此在工业领域应用广泛;二代扫描方式缺点是扫描时间长,射线利用率低。
目前根据二代CT扫描方式的特点,分析并比较了二代扫描的平移和分度运动之间的关系,提出了提高系统扫描效率的两种方案,其中将单向平移改为双向平移的扫描方式能显著提高扫描的效率。虽然改进的扫描方式使扫描时间下降了,但是仍然采集了许多冗余的数据,因此,需要寻找一种更高效率的二代CT扫描方式,提高工件的CT扫描检测速度,而且与传统的二代扫描方式相比,必须保证新的二代扫描方式不会导致CT图像质量降低。
发明内容
本发明的目的就是提供一种二代CT扫描成像方法,它对线阵CT机的二代扫描方式进行了改进,缩短了CT扫描平移距离,扫描检测速度快。
本发明的目的是通过这样的技术方案实现的,它包括有用于发射射线的射线源、用于接收光源的线阵探测器和用于承载待测工件的扫描转台,具体步骤如下:
1)将工件放置在扫描转台上,初始化CT扫描参数,设置CT扫描的视场大小,视场大小小于扫描转台中心到后准直器的距离;
2)计算有效视场大小,将射线扇束的中心线平移到有效视场的左边并与有效视场相切;
3)平移射线源,直到射线扇束的中心线移到有效视场的右边并与有效视场相切,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
4)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤5);
5)平移射线源,使射线扇束从有效视场的右端平移到左端,当射线扇束中心线平移到与有效视场左边缘相切时,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
6)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤3);
7)对步骤3)和步骤5)采集到的投影数据进行重排;
8)初始化重建参数,对步骤7)重排后的数据采用SART迭代重建算法重建CT图像。
进一步,步骤1中初始化CT扫描参数包括有被检工件的视场直径、平移步距和数据采集时间。
进一步,有效视场大小的计算方法为:设R为视场半径,L为射线扇束的宽度;
当2R>L时,有效视场直径为2R,该视场为大视场;
当2R<L时,有效视场直径为L,该视场为小视场;
当2R=L时,有效视场直径为2R或L,该视场为临界视场。
进一步,步骤7)中投影数据进行重排的方法为:将每个线阵探测器平移一次采集的数据提出来重排成一个平行束,自右至左扫描的投影数据进行左右颠倒;若扇束有N个线阵探测器,一次平移得到N个平行束;对所有平移进行重排,若有M个平移,则重排后有M×N个平行束。
进一步,步骤8)中初始化重建的参数包括有有效线阵探测器的个数、有效扇束张角、平移步距和重建图像的大小。
进一步,步骤8)中所述SART迭代重建算法的步骤为:
8-1)对所有的投影数据取对数操作;
8-2)利用以下公式,对图像进行迭代重建:
式中,分别为第 与 次子迭代过程中的第个像素值,表示第条射线穿过第个像素的长度,表示第条射线穿过第个像素的长度, 表示第 条射线的实测投影值, 为第 条射线的模拟投影值,表示重建图像中像素的个数,表示取遍所有像素, 表示第个投影视角下所有射线索引的集合,为松弛因子。
由于采用了上述技术方案,本发明具有如下的优点:
本发明对现有线阵CT机的二代扫描方式进行了改进,使CT扫描平移距离缩短,扫描效率提高,减少了平移床身的距离,降低了硬件成本,CT机的控制代码只需要更改起止点即可;图像重建算法采用SART算法,能有效的抑制伪影和噪声,重建图像质量好,分辨率高。本发明具有扫描检测速度快、成本低、图像质量好的优点。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
本发明的附图说明如下。
图1为本发明扫描成像的流程框图;
图2为本发明在CT扫描大视场时的状态示意图;
图3为本发明在CT扫描小视场时的状态示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
一种二代CT扫描成像方法,包括有用于发射射线的射线源、用于接收光源的线阵探测器和用于承载待测工件的扫描转台,具体步骤如下:
1)将工件放置在扫描转台上,初始化CT扫描参数,设置CT扫描的视场大小,视场大小小于扫描转台中心到后准直器的距离;
2)计算有效视场大小,将射线扇束的中心线平移到有效视场的左边并与有效视场相切;
3)平移射线源,直到射线扇束的中心线移到有效视场的右边并与有效视场相切,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
4)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤5);
5)平移射线源,使射线扇束从有效视场的右端平移到左端,当射线扇束中心线平移到与有效视场左边缘相切时,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
6)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤3);
7)对步骤3)和步骤5)采集到的投影数据进行重排;
8)初始化重建参数,对步骤7)重排后的数据采用SART迭代重建算法重建CT图像。
本发明对现有线阵CT机的二代扫描方式进行了改进,使CT扫描平移距离缩短,扫描效率提高,减少了平移床身的距离,降低了硬件成本,CT机的控制代码只需要更改起止点即可;图像重建算法采用SART算法,能有效的抑制伪影和噪声,重建图像质量好,分辨率高。本发明具有扫描检测速度快、成本低、图像质量好的优点。
步骤1中初始化CT扫描参数包括有被检工件的视场直径、平移步距和数据采集时间。
有效视场大小的计算方法为:设R为视场半径,L为射线扇束的宽度;
当2R>L时,有效视场直径为2R,该视场为大视场;
当2R<L时,有效视场直径为L,该视场为小视场;
当2R=L时,有效视场直径为2R或L,该视场为临界视场。
步骤7)中投影数据进行重排的方法为:将每个线阵探测器平移一次采集的数据提出来重排成一个平行束,自右至左扫描的投影数据进行左右颠倒;若扇束有N个线阵探测器,一次平移得到N个平行束;对所有平移进行重排,若有M个平移,则重排后有M×N个平行束。
步骤8)中初始化重建的参数包括有有效线阵探测器的个数、有效扇束张角、平移步距和重建图像的大小。
步骤8)中所述SART迭代重建算法的步骤为:
8-1)对所有的投影数据取对数操作;
8-2)利用以下公式,对图像进行迭代重建:
式中,分别为第 与 次子迭代过程中的第个像素值,表示第条射线穿过第个像素的长度,表示第条射线穿过第个像素的长度, 表示第 条射线的实测投影值, 为第 条射线的模拟投影值,表示重建图像中像素的个数,表示取遍所有像素, 表示第个投影视角下所有射线索引的集合,为松弛因子。
迭代的具体步骤如下:
取第一个平行束投影,生成对应的模拟投影数据,计算实际投影与模拟投影的差值,将差值乘上松驰因子后反投影到图像上,获得修正图像。同理,取与第一个平行束正交的第二个平行束投影,在修正的图像上进行模拟投影,将实际投影与模拟投影之差乘以松驰因子,再反投影到图像上,获得修正图像;如果有个平行束投影,那么经过次修正将获得图像。
将估计值的模拟投影与实际投影之间的距离作为判断是否进行下一轮迭代的判据;如果距离大于给定阈值,则需要进行下一轮迭代,将赋为初值,按上述相同的步骤进行迭代。如果计算的距离小于给定的阈值,则停止迭代,将重建结果保存到体数据中。
例如:对于CT扫描中的临界视场,L=2R,传统的扫描方式中,一次平移的距离大于L+2R=2L,而本发明的扫描过程中,一次平移的距离为L,因此一次平衡的时间减少了一半多,效率提高了一倍多,冗余数据也显著的减少。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1. 一种二代CT扫描成像方法,包括有用于发射射线的射线源、用于接收光源的线阵探测器和用于承载待测工件的扫描转台,其特征在于,具体步骤如下:
1)将工件放置在扫描转台上,初始化CT扫描参数,设置CT扫描的视场大小,视场大小小于扫描转台中心到后准直器的距离;
2)计算有效视场大小,将射线扇束的中心线平移到有效视场的左边并与有效视场相切;
3)平移射线源,直到射线扇束的中心线移到有效视场的右边并与有效视场相切,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
4)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤5);
5)平移射线源,使射线扇束从有效视场的右端平移到左端,当射线扇束中心线平移到与有效视场左边缘相切时,停止当前平移扫描,线阵探测器在平移空间上均匀的采集投影数据;
6)旋转射线源,使射线扇束绕有效视场中心旋转一个扇角的角度,并判断射线源是否已经旋转完一周,若是则转向步骤7),若不是则转向步骤3);
7)对步骤3)和步骤5)采集到的投影数据进行重排;
8)初始化重建参数,对步骤7)重排后的数据采用SART迭代重建算法重建CT图像。
2. 如权利要求1所述的一种二代CT扫描成像方法,其特征在于:步骤1中初始化CT扫描参数包括有被检工件的视场直径、平移步距和数据采集时间。
3. 如权利要求1所述的一种二代CT扫描成像方法,其特征在于,有效视场大小的计算方法为:设R为视场半径,L为射线扇束的宽度;
当2R>L时,有效视场直径为2R,该视场为大视场;
当2R<L时,有效视场直径为L,该视场为小视场;
当2R=L时,有效视场直径为2R或L,该视场为临界视场。
4. 如权利要求1所述的一种二代CT扫描成像方法,其特征在于,步骤7)中投影数据进行重排的方法为:将每个线阵探测器平移一次采集的数据提出来重排成一个平行束,自右至左扫描的投影数据进行左右颠倒;若扇束有N个线阵探测器,一次平移得到N个平行束;对所有平移进行重排,若有M个平移,则重排后有M×N个平行束。
5. 如权利要求1所述的一种二代CT扫描成像方法,其特征在于,步骤8)中初始化重建的参数包括有有效线阵探测器的个数、有效扇束张角、平移步距和重建图像的大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310063402.1A CN103163165B (zh) | 2013-02-28 | 2013-02-28 | 一种二代ct扫描成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310063402.1A CN103163165B (zh) | 2013-02-28 | 2013-02-28 | 一种二代ct扫描成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103163165A true CN103163165A (zh) | 2013-06-19 |
CN103163165B CN103163165B (zh) | 2014-11-26 |
Family
ID=48586428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310063402.1A Active CN103163165B (zh) | 2013-02-28 | 2013-02-28 | 一种二代ct扫描成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103163165B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103357588A (zh) * | 2013-08-06 | 2013-10-23 | 重庆大学 | X射线工业ct设备准直片分组方法 |
CN103455989A (zh) * | 2013-09-24 | 2013-12-18 | 南京大学 | 一种结合超声图像提高有限角度ct成像质量的方法 |
CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN104698016A (zh) * | 2015-03-26 | 2015-06-10 | 重庆大学 | 一种多次平移的交错螺旋工业ct扫描成像方法 |
CN105136823A (zh) * | 2015-07-07 | 2015-12-09 | 重庆大学 | 大口径管道壁外部ct局部扫描成像方法 |
CN105809732A (zh) * | 2016-03-21 | 2016-07-27 | 重庆大学 | 一种管道ct图像的体绘制方法 |
CN106691485A (zh) * | 2016-07-18 | 2017-05-24 | 山东省肿瘤防治研究院 | 一种放射影像引导下诊治肺肿瘤装置 |
CN107764846A (zh) * | 2017-10-20 | 2018-03-06 | 重庆大学 | 一种正交直线扫描的cl成像系统及分析方法 |
CN110428478A (zh) * | 2019-07-15 | 2019-11-08 | 清华大学 | 交替光源扇束x射线ct采样方法及装置 |
CN111839568A (zh) * | 2020-07-22 | 2020-10-30 | 重庆大学 | 一种新型大视场直线扫描ct系统及图像重建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070230657A1 (en) * | 2006-04-04 | 2007-10-04 | Walter Garms | Apparatus and method for controlling start and stop operations of a computed tomography imaging system |
CN102331433A (zh) * | 2011-05-30 | 2012-01-25 | 重庆大学 | 大尺寸工业长管道管壁的外部螺旋锥束ct扫描成像方法 |
CN102590243A (zh) * | 2012-02-17 | 2012-07-18 | 重庆大学 | 一种铁路铸件全身ct扫描成像方法 |
-
2013
- 2013-02-28 CN CN201310063402.1A patent/CN103163165B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070230657A1 (en) * | 2006-04-04 | 2007-10-04 | Walter Garms | Apparatus and method for controlling start and stop operations of a computed tomography imaging system |
CN102331433A (zh) * | 2011-05-30 | 2012-01-25 | 重庆大学 | 大尺寸工业长管道管壁的外部螺旋锥束ct扫描成像方法 |
CN102590243A (zh) * | 2012-02-17 | 2012-07-18 | 重庆大学 | 一种铁路铸件全身ct扫描成像方法 |
Non-Patent Citations (3)
Title |
---|
冯鹏等: "工业X-CT二代扫描运动控制系统及其仿真", 《重庆大学学报》 * |
魏东波等: "大尺寸构件工业CT成像方法", 《北京航空航天大学学报》 * |
魏彪等: "X 光工业用CT 成像二代扫描运动控制系统设计", 《光学精密工程》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103357588B (zh) * | 2013-08-06 | 2015-04-29 | 重庆大学 | X射线工业ct设备准直片分组方法 |
CN103357588A (zh) * | 2013-08-06 | 2013-10-23 | 重庆大学 | X射线工业ct设备准直片分组方法 |
CN103455989A (zh) * | 2013-09-24 | 2013-12-18 | 南京大学 | 一种结合超声图像提高有限角度ct成像质量的方法 |
CN103776851B (zh) * | 2014-02-07 | 2016-08-17 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN104698016B (zh) * | 2015-03-26 | 2017-04-05 | 重庆大学 | 一种多次平移的交错螺旋工业ct扫描成像方法 |
CN104698016A (zh) * | 2015-03-26 | 2015-06-10 | 重庆大学 | 一种多次平移的交错螺旋工业ct扫描成像方法 |
CN105136823A (zh) * | 2015-07-07 | 2015-12-09 | 重庆大学 | 大口径管道壁外部ct局部扫描成像方法 |
CN105136823B (zh) * | 2015-07-07 | 2017-12-05 | 重庆大学 | 大口径管道壁外部ct局部扫描成像方法 |
CN105809732A (zh) * | 2016-03-21 | 2016-07-27 | 重庆大学 | 一种管道ct图像的体绘制方法 |
CN105809732B (zh) * | 2016-03-21 | 2018-12-21 | 重庆大学 | 一种管道ct图像的体绘制方法 |
CN106691485A (zh) * | 2016-07-18 | 2017-05-24 | 山东省肿瘤防治研究院 | 一种放射影像引导下诊治肺肿瘤装置 |
CN107764846A (zh) * | 2017-10-20 | 2018-03-06 | 重庆大学 | 一种正交直线扫描的cl成像系统及分析方法 |
CN107764846B (zh) * | 2017-10-20 | 2020-04-14 | 重庆大学 | 一种正交直线扫描的cl成像系统及分析方法 |
CN110428478A (zh) * | 2019-07-15 | 2019-11-08 | 清华大学 | 交替光源扇束x射线ct采样方法及装置 |
CN111839568A (zh) * | 2020-07-22 | 2020-10-30 | 重庆大学 | 一种新型大视场直线扫描ct系统及图像重建方法 |
CN111839568B (zh) * | 2020-07-22 | 2023-12-12 | 重庆大学 | 一种新型大视场直线扫描ct系统及图像重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103163165B (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103163165B (zh) | 一种二代ct扫描成像方法 | |
JP5198443B2 (ja) | 画像の分解能を高めるシステム及び方法 | |
RU2510080C2 (ru) | Устройство для обработки изображения, способ обработки изображения и среда долговременного хранения информации | |
US8284892B2 (en) | System and method for image reconstruction | |
EP1800264B1 (en) | Image reconstruction with voxel dependent interpolation | |
US9824467B2 (en) | Iterative image reconstruction | |
US8805037B2 (en) | Method and system for reconstruction of tomographic images | |
JP4576032B2 (ja) | 2回パス式コーン・ビーム画像再構成の方法及び装置 | |
US20110228898A1 (en) | Inverse geometry volume computed tomography systems | |
US9042514B2 (en) | Dose reduction via dynamic collimation adjustment for targeted field of view and/or digital tilt CT | |
US7978810B2 (en) | Imaging method for variable pitch spiral CT and a CT machine for carrying out the method | |
US8615121B2 (en) | Reconstruction of projection data to generate tomographic images having improved frequency characteristics | |
US8948337B2 (en) | Computed tomography image reconstruction | |
CN103180879A (zh) | 用于从投影数据对对象进行混合重建的设备和方法 | |
CN1742296A (zh) | 锥形光束计算机断层成像的立体像素驱动的螺旋重建方法 | |
CN110057847A (zh) | Tr层析扫描投影重排方法及装置 | |
Park et al. | A fully GPU-based ray-driven backprojector via a ray-culling scheme with voxel-level parallelization for cone-beam CT reconstruction | |
Banjak | X-ray computed tomography reconstruction on non-standard trajectories for robotized inspection | |
CN106097411B (zh) | Ct机图像重建方法及高分辨ct扫描机 | |
Gonzales et al. | Rectangular computed tomography using a stationary array of CNT emitters: initial experimental results | |
US7272205B2 (en) | Methods, apparatus, and software to facilitate computing the elements of a forward projection matrix | |
Varslot et al. | Fast high-resolution micro-CT with exact reconstruction methods | |
US8379948B2 (en) | Methods and systems for fast iterative reconstruction using separable system models | |
KR101914522B1 (ko) | 이진 이동 빔 차단기를 이용하여 산란을 보정하는 콘빔 전산화 단층 촬영 시스템 | |
CN109685867A (zh) | 一种ct系统、图像重建方法、装置、设备及介质 |
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 |