CN101398397A - 多次扫描模式的ct成像方法 - Google Patents
多次扫描模式的ct成像方法 Download PDFInfo
- Publication number
- CN101398397A CN101398397A CNA2007101755884A CN200710175588A CN101398397A CN 101398397 A CN101398397 A CN 101398397A CN A2007101755884 A CNA2007101755884 A CN A2007101755884A CN 200710175588 A CN200710175588 A CN 200710175588A CN 101398397 A CN101398397 A CN 101398397A
- Authority
- CN
- China
- Prior art keywords
- detector
- image
- data
- theta
- scanning
- 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
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
Abstract
本发明是一种基于多次扫描的CT成像方法,用于CT设备的扫描视野不能覆盖被检测物体时的CT成像(即超视野CT成像)。该发明包括多次扫描采集得到投影数据集和运用扫描重建算法处理投影数据集获得被扫描物体的CT图像,其特征在于对被测物体扫描时,允许每次扫描所覆盖的区域存在相叠部分,处理投影数据集的重建算法不需要对数据重排。其主要优点是被测物体的大小不受CT扫描视野限制,扫描模式灵活,由于重建算法不需要对数据重排,避免了重排所需的插值计算和重排引起的重建图像降质,因此重建速度快,重建图像分辨率高。
Description
技术领域
本发明属于X射线CT技术领域,用于CT设备的扫描视野不能覆盖被检测物体时的CT成像(即超视野CT成像)。
背景技术
在工业无损检测中,常遇到使用短的线阵探测器检测大尺寸物体的问题。针对大尺寸被测物体不能被扫描视野完全覆盖的扫描问题,传统的解决办法是采用二代扫描的方式,即平移加旋转的扫描方式,但它在采集被测物体边缘数据时,要求回转区与视场接近外切,不仅造成视场和射线的浪费,也增加了扫描时间和数据保存传输的成本。
针对二代扫描方式存在的问题,人们提出了旋转加平移的扫描方式。旋转加平移的扫描方式减少了资源浪费,提高了扫描速度,并且与三代扇束扫描方式相比,由于从光源发出的射线不是各向同性的,采用旋转加平移的扫描方法可以减小长的线阵探测器上不同位置的探测器单元探测到的射线强度的差异。
旋转加平移方式的重建方法与二代扫描方式的重建算法相似,都是基于滤波反投影(FBP)思想,为了重建CT图像,这些算法不得不将多组多次扫描数据通过插值重排成平行束或扇形束数据。众所周知,数据重排不仅大大增加了计算量,而且会降低重建图像的分辨率,因此二代扫描方式和旋转加平移扫描方式都存在数据计算量大,图像重建速度慢,重建图像分辨率低的缺点。
发明内容
本发明的目的在于提供一种在CT设备的扫描视野不能覆盖被检测物体时的CT成像(即超视野CT成像)方法,其能够快速重建高分辨率的被测物体断层图像。
为了实现上述发明目的,本发明的技术方案如下:
一种多次扫描模式的CT成像方法,包括多次扫描采集得到投影数据集和运用扫描重建算法处理投影数据集获得被扫描物体的CT图像,其特征在于对被测物体扫描时,允许每次扫描所覆盖的区域存在相叠部分,处理投影数据集的重建算法不需要对数据重排,具体实现步骤如下:
(1)根据被测物体的直径确定CT系统的几何扫描参数,以及转台中心需要平移的距离和次数,转台中心平移的次数即偏置次数M;
(2)将转台中心依次移至不同的位置,在每个位置上,转台旋转时射线源出束并由探测器采集数据,直到数据采集完毕得到M组数据,
扫描时,射线源和探测器固定不动,被扫描物体放置在转台上,转台旋转中心的初始位置O位于中心射线上,转台旋转中心由O沿垂直于中心射线的x1轴方向或反方向分别平移至Om,其中m=1,2,L,M,它们满足|OO1|<lRO/(2RD),且|OmOm+1|<lRO/RD,其中m=1,2,L,M-1,l为探测器长度,RO为光源位置S0到O点的距离,RD为S0到探测器的距离,当转台分别绕着Om旋转一周时,探测器采集穿过被测物体后的射线流量,最终获得旋转中心分别在M个Om位置时采集的投影数据,共采集M组投影数据;
(3)针对射束强度的非一致性、探测器的非一致性和探测器本底,对采集的M组投影数据进行校正;
(4)步骤(2)采集得到的M组投影数据经过步骤(3)校正后,使用下面的重建算法对各组投影数据分别进行图像重建并得到被测物体的M个部分DBP(differentiated backprojection)图像,上述M个部分DBP图像累加,可得到被测物体的整个DBP图像,然后利用有限区域上的Hilbert变换的逆,即可获得整个被测物体的断层图像,
记光源焦点位置为S0,针对第m次扫描,转台相对于初始位置的平移距离为hm,即|OOm|=hm,β是转台逆时针旋转的角度,我们定义第m次扫描时虚拟探测器坐标轴为um,其方向均为β=(cosβ,sinβ),而坐标原点位于O,记u0,m为点x0在第m个虚拟探测器轴上的投影位置,用pm(β,um)分别表示第m组扇束获得的投影数据,pm(β,um)在um∈[-lRO/(2RD),lRO/(2RD)]时是已知的,其中m=1,2,L,M,
那么由第m组数据pm(β,um)重建被测物体的第m个部分DBP图像bq,m(x0)(m=1,2,L,M)表示如下,
其中θ是任意一个角度,x0=(x1,x2),β⊥=(-sinβ,cosβ), (m=1,2,L,M),T0=0, ε是一个小的正数,kε(r)是无穷可微的,且具有性质:(i)当r≥ε时,ke(r)=1,(ii)当-ε<r<ε时,ke(r)是单调递增的,(iii)当r≤-ε时,ke(r)=0,例如,
其中w(r)为磨光函数,由下式定义,
上述M个部分DBP图像累加,即可得到整个DBP图像,
记Hθf(x0)表示图像沿θ=(cosθ,sinθ)方向的Hilbert变换,那么根据DBP图像与Hilbert变换之间的关系bθ(x0)=-2πHθf(x0),通过如下有限区间上的Hilbert逆变换得到断层图像,
其中θ=(cosθ,sinθ)和θ⊥=(-sinθ,cosθ),t=x·θ,Lt,Ut(满足Ut>Lt)和Ct是一个与t有关的常数,Ct的值可以由直线上x·θ=t的区间x·θ⊥∈[Lt,Lt+εt]∪[Ut-εt,Ut]上的Hθf(x)决定。
上述技术方案的优点是:被测物体的大小不受CT扫描视野限制;允许多次扫描数据相叠,扫描模式灵活;重建算法不需要对数据重排,避免了重排所需的插值计算和重排引起的重建图像降质,因此重建速度快,重建图像分辨率高。被扫描物体的直径可超出扫描视野的限制,被扫描物体的直径一般能比扫描视野大3-5倍。重建算法对于线阵探测器所采集的数据是精确重建算法;对于面探测器所采集的数据,对中平面的CT图像是精确重建的,对非中平面的CT图像是近似重建的。
下面对本发明作进一步的说明:如图2(a)所示,设S0为射线源的焦点,粗黑线表示等距线阵探测器,OD为探测器的中心。设直线ODS0垂直于线阵探测器,右手坐标系Ox1x2定义如下:将位于直线ODS0上的O设为坐标原点,x2轴与向量ODS0同向,x1轴垂直于向量ODS0,记RO为S0与O之间的距离,RD为S0与OD之间的距离,l为探测器的长度。
为了方便公式的推导,我们引入等价扫描模式的虚拟扫描模式,如图2(b)所示。显然,“由射线源和探测器构成的扇束固定时,转台沿着x1轴的平移或绕O点的顺时针旋转”分别等价于“转台固定时,相应扇束沿着x1轴的反向的平移或绕着O点的逆时针旋转”。因此,图2(a)所示的扫描模式等价于图2(b)所示的虚拟扫描模式:被测物体固定,S1与S2相应的扇束同时绕O做360度旋转时探测器采集数据,其中位于左侧的扇束的焦点S1在以O为圆心、以为半径的圆周上运动;而位于右侧的扇束的焦点S2在以O为圆心、以为半径的圆周上运动;S1与S2的连线平行于探测器所在的直线。类似地,可以得到转台单侧多次偏置的RT扫描模式等价的虚拟扫描模式。
记第m个光源焦点为Sm,扇束的每次平移距离为hm,即|OOm|=hm;β是直线O1S1与x2轴所形成的逆时针角。我们定义M个虚拟探测器坐标轴,第m个记为um,其方向均为β=(cosβ,sinβ),而坐标原点分别位于Om(Om位于x1轴正方向上)。记u0,m为焦点Sm下点x0在第m个虚拟探测器轴上的投影位置,用pm(β,um)分别表示以Sm为焦点的第m组扇束获取的投影数据,需记住:pm(β,um)在um∈[-lRO/(2RD),lRO/(2RD)]时是已知的,其中m=1,2,L,M。
附图说明
图1为本发明多次扫描模式的CT成像方法的实现流程图;
图2是转台左侧两次偏置的RT扫描模式示意图,图2(a)为两次扫描模式中,转台中心在不同位置时的扫描示意图,O1和O2为转台中心,虚线表示扫描时转台所在的实际位置,图2(b)为两次扫描模式中,转台中心固定不动,射线源和探测器移动的等价模式示意图;
图3为用于测试本发明CT成像方法效果的模型照片;
图4为两次扫描模式中模型的扫描数据图,图4(a)为对应图2(a)中转台位于位置O1初时扫描到的投影数据;图4(b)为对应图2(a)中转台位于位置O2初时扫描到的投影数据;
图5为两次扫描模式中由扫描数据重建的DBP图像,图5(a)为由图4(a)数据重建的部分DBP图像,图5(b)为由图4(b)数据重建的部分DBP图像,图5(c)为两个部分DBP图像(图5(a)和图5(b))累加后得到的整体DBP图像;
图6为本发明重建图像与重排算法重建图像的比较图,图6(a)为重排方法两次扫描重建的图像,图6(b)为图6(a)的局部放大图像,图6(c)为本发明方法两次扫描重建的图像,图6(d)为图6(c)的局部放大图像。
具体实施方式
多次扫描模式的CT成像方法,步骤(1)中转台中心需要平移的距离和次数的具体计算方法是:如图2(a)所示,RO为射线源S0与转台旋转中心O之间的距离,RD为射线源S0与探测器中心OD之间的距离,l为探测器的长度,因此单次扫描的扫描视野半径 假设被测物的最大半径为ObjectR,则转台需要平移的次数为 。本实施例以两次扫描模式为例即M=2。
步骤(2)中转台所移至的位置,参照图2(a),图2(a)中O1,O2的位置可以根据转台平移的距离确定,O1和O2所确定直线与射线源和探测器中心所确定的直线垂直,交点为O,转台分别在O1和O2处旋转采集得到两组数据Idata_1(β,u),Idata_2(β,u)。
步骤(3)中数据非一致性和探测器本底的具体校正方法是:当射线源停止时,探测器采集暗场数据,多次采集平均后得到暗场的平均数据Inoise(u),移动转台将被测物体垂直于主射线方向移开,使得扫描视野中没有被测物,打开射线源和探测器,采集亮场的数据,多次采集平均后得到亮场的平均数据Inormal(u),将亮场的平均数据Inormal(u)减去暗场的平均数据Inoise(u),即对亮场数据进行暗电流校正,得到校正后的数据I0(u),将两次扫描数据Idata_1(β,u)和Idata_2(β,u)减去暗场的平均数据Inoise(u),即对两次扫描数据进行暗电流校正,得到校正后的数据I_1(β,u),I_2(β,u),然后除亮场的校正后数据I0(u),取对数得到两次扫描的投影数据p1(b,u),p2(b,u),在等价模式中投影数据表示成p1(b,u1),p2(b,u2)。
对于步骤(3)得到的两组数据,运用本发明重建算法进行处理,具体实现方法如下:对投影数据p1(b,u1),p2(b,u2),当b固定时数据对u1,u2求数值导数,得到投影数据的导数,然后做加权反投影,得到被测物体的两个部分DBP图像,累加后得到被测物体的整个DBP图像,再利用有限区域上的Hilbert变换的逆,从而完成了被测物体断层图像的重建。
根据公式(1),(2)可得到M=2即两次扫描的DBP重建公式,取θ=0:
其中β⊥=(-sinβ,cosβ),
累加后获得被测物体的整个DBP图像:
b0(x0)=b0,1(x0)+b0,2(x0)
对被测物体的整个DBP图像利用有限区域上的Hilbert变换的逆,即公式(3)获得被测物体的重建断层图像。
参照图6,被扫描物体直径大约是207mm,图6给出了经过本发明的两次扫描模式所重建的结果与传统的重排方法重建结果的比较。系统的扫描参数为有效探测器的长度为l=1172×0.127mm,RO=1660mm,RD=1900mm,单次圆轨道扫描的扫描视野半径Rfov=64.97mm,两次扫描模式的扫描视野半径大约是259.88mm。图6(a)为重排方法重建图像,图6(b)为图6(a)的局部放大图像,图6(c)为本发明方法两次扫描重建图像,图6(d)为图6(c)的局部放大图像。
Claims (4)
1.一种多次扫描模式的CT成像方法,包括多次扫描采集得到投影数据集和运用扫描重建算法处理投影数据集获得被扫描物体的CT图像,其特征在于对被测物体扫描时,允许每次扫描所覆盖的区域存在相叠部分,处理投影数据集的重建算法不需要对数据重排,具体实现步骤如下:
(1)根据被测物体的直径确定CT系统的几何扫描参数,以及转台中心需要平移的距离和次数,转台中心平移的次数即偏置次数M;
(2)将转台中心依次移至不同的位置,在每个位置上,转台旋转时射线源出束并由探测器采集数据,直到数据采集完毕得到M组数据,
扫描时,射线源和探测器固定不动,被扫描物体放置在转台上,转台旋转中心的初始位置O位于中心射线上,转台旋转中心由O沿垂直于中心射线的x1轴方向或反方向分别平移至Om,其中m=1,2,L,M,它们满足|OO1|<lRO/(2RD),且|OmOm+1|<lRO/RD,其中m=1,2,L,M-1,l为探测器长度,RO为光源位置S0到O点的距离,RD为S0到探测器的距离,当转台分别绕着Om旋转一周时,探测器采集穿过被测物体后的射线流量,最终获得旋转中心分别在M个Om位置时采集的投影数据,共采集M组投影数据;
(3)针对射束强度的非一致性、探测器的非一致性和探测器本底,对采集的M组投影数据进行校正;
(4)步骤(2)采集得到的M组投影数据经过步骤(3)校正后,使用下面的重建算法对各组投影数据分别进行图像重建并得到被测物体的M个部分DBP图像,上述M个部分DBP图像累加,可得到被测物体的整个DBP图像,然后利用有限区域上的Hilbert变换的逆,即可获得整个被测物体的断层图像,
记光源焦点位置为S0,针对第m次扫描,转台相对于初始位置的平移距离为hm,即|OOm|=hm,β是转台逆时针旋转的角度,我们定义第m次扫描时虚拟探测器坐标轴为um,其方向均为β=(cosβ,sinβ),而坐标原点位于O,记u0,m为点x0在第m个虚拟探测器轴上的投影位置,用pm(β,um)分别表示第m组扇束获得的投影数据,pm(β,um)在um∈[-lRO/(2RD),lRO/(2RD)]时是已知的,其中m=1,2,L,M,
那么由第m组数据pm(β,um)重建被测物体的第m个部分DBP图像bq,m(x0)(m=1,2,L,M)表示如下,
其中θ是任意一个角度,x0=(x1,x2),β⊥=(-sinβ,cosβ), (m=1,2,L,M),T0=0, ε是一个小的正数,kε(r)是无穷可微的,且具有性质:(i)当r≥ε时,ke(r)=1,(ii)当-ε<r<ε时,ke(r)是单调递增的,(iii)当r≤-ε时,ke(r)=0,如,
其中w(r)为磨光函数,由下式定义,
上述M个部分DBP图像累加,即可得到整个DBP图像,
记Hθf(x0)表示图像沿θ=(cosθ,sinθ)方向的Hilbert变换,那么根据DBP图像与Hilbert变换之间的关系bθ(x0)=-2πHθf(x0),通过如下有限区间上的Hilbert逆变换得到断层图像,
(3)
其中θ=(cosθ,sinθ)和θ⊥=(-sinθ,cosθ),t=x·θ,Lt,Ut(满足Ut>Lt)和Ct是一个与t有关的常数,Ct的值可以由直线上x·θ=t的区间x·θ⊥∈[Lt,Lt+εt]∪[Ut-εt,Ut]上的Hθf(x)决定。
2.根据权利要求1所述的CT成像方法,其特征在于扫描模式所采集的多组数据允许是单侧覆盖的,即多次扫描的视野在每个角度下只覆盖被测物体在旋转中心一侧的部分。
3.根据权利要求1或2所述的CT成像方法,其特征在于步骤(3)中数据校正不但考虑了探测器的非一致性和探测器本底,而且还对射束在不同时刻、不同方向的强度的非一致性进行了校正。
4.根据权利要求3所述的CT成像方法,其特征在于重建算法对于线阵探测器所采集的数据是精确重建算法;对于面探测器所采集的数据,对中平面的CT图像是精确重建的,对非中平面的CT图像是近似重建的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101755884A CN101398397B (zh) | 2007-09-30 | 2007-09-30 | 多次扫描模式的ct成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101755884A CN101398397B (zh) | 2007-09-30 | 2007-09-30 | 多次扫描模式的ct成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101398397A true CN101398397A (zh) | 2009-04-01 |
CN101398397B CN101398397B (zh) | 2012-05-30 |
Family
ID=40517102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101755884A Expired - Fee Related CN101398397B (zh) | 2007-09-30 | 2007-09-30 | 多次扫描模式的ct成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101398397B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062740A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院深圳先进技术研究院 | 锥束ct扫描成像方法及系统 |
CN104274201A (zh) * | 2014-10-10 | 2015-01-14 | 深圳先进技术研究院 | 乳腺层析成像方法和系统及成像设备和图像采集处理方法 |
CN104597061A (zh) * | 2015-01-28 | 2015-05-06 | 中国工程物理研究院应用电子学研究所 | 基于虚拟探测器的大视场锥束ct成像方法 |
CN105266839A (zh) * | 2015-11-23 | 2016-01-27 | 山东科技大学 | 一种三源圆轨迹半径不同的大视野ct成像方法 |
CN105761226A (zh) * | 2016-02-24 | 2016-07-13 | 赛诺威盛科技(北京)有限公司 | 一种超视野ct扫描图像的补偿重建方法 |
US9993215B2 (en) | 2014-11-26 | 2018-06-12 | Shenyang Neusoft Medical Systems Co., Ltd. | CT image correction |
CN110146524A (zh) * | 2019-04-10 | 2019-08-20 | 清华大学 | 基于逆康普顿散射源的ct扫描与重建方法以及成像系统 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN110621985A (zh) * | 2017-07-03 | 2019-12-27 | 株式会社岛津制作所 | X线计算机断层装置 |
CN110693513A (zh) * | 2018-10-27 | 2020-01-17 | 上海联影医疗科技有限公司 | 多模态医学系统的控制方法、系统及存储介质 |
CN111595875A (zh) * | 2020-04-21 | 2020-08-28 | 中国工程物理研究院应用电子学研究所 | 一种大视场ct成像重建方法 |
CN113533392A (zh) * | 2021-07-12 | 2021-10-22 | 重庆大学 | 一种组合扫描cl成像方法 |
CN113827262A (zh) * | 2021-09-29 | 2021-12-24 | 纳米维景(上海)医疗科技有限公司 | 一种新型ct扫描区域重建方法、存储介质及终端设备 |
CN115598152A (zh) * | 2021-07-07 | 2023-01-13 | 同方威视技术股份有限公司(Cn) | 射线成像系统以及射线成像方法 |
US11654304B2 (en) | 2018-12-17 | 2023-05-23 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining a region of interest of a subject |
CN116543071A (zh) * | 2023-07-06 | 2023-08-04 | 有方(合肥)医疗科技有限公司 | 大视野ct成像方法及装置 |
CN117838169A (zh) * | 2024-03-08 | 2024-04-09 | 江苏一影医疗设备有限公司 | 一种基于站立位cbct的成像方法、系统和设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6856666B2 (en) * | 2002-10-04 | 2005-02-15 | Ge Medical Systems Global Technology Company, Llc | Multi modality imaging methods and apparatus |
CN100457041C (zh) * | 2005-11-17 | 2009-02-04 | 中国科学院高能物理研究所 | X射线折射衬度ct数据采集方法和重建方法 |
-
2007
- 2007-09-30 CN CN2007101755884A patent/CN101398397B/zh not_active Expired - Fee Related
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062740A (zh) * | 2010-12-01 | 2011-05-18 | 中国科学院深圳先进技术研究院 | 锥束ct扫描成像方法及系统 |
CN102062740B (zh) * | 2010-12-01 | 2013-02-13 | 中国科学院深圳先进技术研究院 | 锥束ct扫描成像方法及系统 |
CN104274201A (zh) * | 2014-10-10 | 2015-01-14 | 深圳先进技术研究院 | 乳腺层析成像方法和系统及成像设备和图像采集处理方法 |
US9993215B2 (en) | 2014-11-26 | 2018-06-12 | Shenyang Neusoft Medical Systems Co., Ltd. | CT image correction |
CN104597061A (zh) * | 2015-01-28 | 2015-05-06 | 中国工程物理研究院应用电子学研究所 | 基于虚拟探测器的大视场锥束ct成像方法 |
CN104597061B (zh) * | 2015-01-28 | 2017-10-31 | 中国工程物理研究院应用电子学研究所 | 基于虚拟探测器的大视场锥束ct成像方法 |
CN105266839A (zh) * | 2015-11-23 | 2016-01-27 | 山东科技大学 | 一种三源圆轨迹半径不同的大视野ct成像方法 |
CN105266839B (zh) * | 2015-11-23 | 2018-04-13 | 山东科技大学 | 一种三源圆轨迹半径不同的大视野ct成像方法 |
CN105761226A (zh) * | 2016-02-24 | 2016-07-13 | 赛诺威盛科技(北京)有限公司 | 一种超视野ct扫描图像的补偿重建方法 |
CN105761226B (zh) * | 2016-02-24 | 2018-05-18 | 赛诺威盛科技(北京)有限公司 | 一种超视野ct扫描图像的补偿重建方法 |
CN110621985A (zh) * | 2017-07-03 | 2019-12-27 | 株式会社岛津制作所 | X线计算机断层装置 |
CN110621985B (zh) * | 2017-07-03 | 2022-03-11 | 株式会社岛津制作所 | X线计算机断层装置 |
CN110693513A (zh) * | 2018-10-27 | 2020-01-17 | 上海联影医疗科技有限公司 | 多模态医学系统的控制方法、系统及存储介质 |
US11654304B2 (en) | 2018-12-17 | 2023-05-23 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining a region of interest of a subject |
CN110146524A (zh) * | 2019-04-10 | 2019-08-20 | 清华大学 | 基于逆康普顿散射源的ct扫描与重建方法以及成像系统 |
CN110146524B (zh) * | 2019-04-10 | 2021-09-28 | 清华大学 | 基于逆康普顿散射源的ct扫描与重建方法以及成像系统 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN111595875A (zh) * | 2020-04-21 | 2020-08-28 | 中国工程物理研究院应用电子学研究所 | 一种大视场ct成像重建方法 |
CN115598152A (zh) * | 2021-07-07 | 2023-01-13 | 同方威视技术股份有限公司(Cn) | 射线成像系统以及射线成像方法 |
CN115598152B (zh) * | 2021-07-07 | 2024-04-19 | 同方威视技术股份有限公司 | 射线成像系统以及射线成像方法 |
CN113533392A (zh) * | 2021-07-12 | 2021-10-22 | 重庆大学 | 一种组合扫描cl成像方法 |
CN113827262A (zh) * | 2021-09-29 | 2021-12-24 | 纳米维景(上海)医疗科技有限公司 | 一种新型ct扫描区域重建方法、存储介质及终端设备 |
CN116543071A (zh) * | 2023-07-06 | 2023-08-04 | 有方(合肥)医疗科技有限公司 | 大视野ct成像方法及装置 |
CN116543071B (zh) * | 2023-07-06 | 2023-09-19 | 有方(合肥)医疗科技有限公司 | 大视野ct成像方法及装置 |
CN117838169A (zh) * | 2024-03-08 | 2024-04-09 | 江苏一影医疗设备有限公司 | 一种基于站立位cbct的成像方法、系统和设备 |
CN117838169B (zh) * | 2024-03-08 | 2024-05-24 | 江苏一影医疗设备有限公司 | 一种基于站立位cbct的成像方法、系统和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN101398397B (zh) | 2012-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101398397B (zh) | 多次扫描模式的ct成像方法 | |
US7251307B2 (en) | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data | |
US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
EP0492895B1 (en) | Reconstructing 3-D images | |
US7203272B2 (en) | Cone-beam filtered backprojection image reconstruction method for short trajectories | |
US7672427B2 (en) | Imaging system | |
EP0492896A2 (en) | Reconstructing a 3-D image | |
US20120301004A1 (en) | Computed Tomography Imaging Process And System | |
US6292525B1 (en) | Use of Hilbert transforms to simplify image reconstruction in a spiral scan cone beam CT imaging system | |
US8724889B2 (en) | Method and apparatus for CT image reconstruction | |
CN101634638B (zh) | 一种探测器偏置的大视野锥束x射线倾斜扫描三维数字成像方法 | |
EP1644897B1 (en) | A fourier tomographic image reconstruction method for fan-beam data | |
WO2004019279A2 (en) | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system | |
US20100119034A1 (en) | Computed tomography device and method using circular-pixel position-adaptive interpolation | |
CN102973291B (zh) | C型臂半精确滤波反投影断层成像方法 | |
Cai et al. | Edge guided image reconstruction in linear scan CT by weighted alternating direction TV minimization | |
US8948337B2 (en) | Computed tomography image reconstruction | |
CN111265231A (zh) | 分布式光源ct图像重建方法与系统 | |
CN103714578A (zh) | 针对半覆盖螺旋锥束ct的单层重排滤波反投影重建方法 | |
US20060034417A1 (en) | Efficient image reconstruction algorithm for the circle and line cone beam computed tomography | |
CN102062740B (zh) | 锥束ct扫描成像方法及系统 | |
US8861829B2 (en) | Method and system for reconstruction of tomographic images | |
US20090154639A1 (en) | System and method for image reconstruction for helical cone beam computed tomography with factorized redundancy weighting | |
CN116630460A (zh) | 源直线扫描轨迹的探测器行向微分高质量图像重建方法 | |
CN102376096A (zh) | Pi线选取和采样方法和装置以及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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120530 Termination date: 20150930 |
|
EXPY | Termination of patent right or utility model |