WO2020210941A1 - 分布式光源ct图像重建方法与系统 - Google Patents
分布式光源ct图像重建方法与系统 Download PDFInfo
- Publication number
- WO2020210941A1 WO2020210941A1 PCT/CN2019/082672 CN2019082672W WO2020210941A1 WO 2020210941 A1 WO2020210941 A1 WO 2020210941A1 CN 2019082672 W CN2019082672 W CN 2019082672W WO 2020210941 A1 WO2020210941 A1 WO 2020210941A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- light source
- projection data
- segment
- linear
- detector
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
Abstract
一种分布式光源CT图像重建方法与系统,分布式光源CT图像重建系统包括:N段直线分布式光源,用于发出射线;以及N段直线分布式探测器,用于获取N段扫描投影数据,N≥2;其中,分布式光源CT图像重建系统利用N段扫描投影数据获取重建图像。分布式光源CT图像重建方法与系统具有更高的空间分辨率和成像速度,能够获得精确的物体重建图像。
Description
本发明属于辐射成像技术领域,具体的涉及一种分布式光源CT图像重建方法与系统。
传统的CT,在成像时,物体需要旋转,通常需要使用滑环,噪音和振动较大,对X射线机、高压模块、探测器等部件的机械强度要求较高,扫描速度也不能较好的满足使用需求。
此外,现有的CT还存在投影数据角度缺失、投影数据重叠、空间分辨率第、重建效果差等问题。
发明内容
(一)要解决的技术问题
本发明提供了一种分布式光源CT图像重建方法与系统,以至少部分解决以上所存在的技术问题。
(二)技术方案
根据本发明的一个方面,提供了一种分布式光源CT图像重建系统,包括:
N段直线分布式光源,用于发出射线;以及
N段直线分布式探测器,用于获取N段扫描投影数据,N≥2;
其中,所述分布式光源CT图像重建系统利用所述N段扫描投影数据获取重建图像。
在一些实施例中,所述系统包括2段直线分布式光源,分别为第一直线分布式光源和第二直线分布式光源、以及2段直线分布式探测器,分别为第一直线分布式探测器和第二直线分布式探测器;其中,第一直线分布式光源和第一直线分布式探测器分别位于待测物体两侧且相对设置,第二直线分布式光源和第二直线分布式探测器分别位于待测物体两侧且相对设置。
根据本发明的另一个方面,提供了一种分布式光源CT图像重建方法,包括:
获取N段扫描投影数据;以及
根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,其中,N≥2。
在一些实施例中,所述获取N段扫描投影数据,包括:利用N段直线分布式光源及N段直线分布式探测器采用静态扫描模式获取N段扫描投影数据。
在一些实施例中,在获取N段扫描投影数据之前,还包括:对直线分布式光源之间的间距及夹角,直线分布式光源与直线分布式探测器之间的间距,直线分布式探测器之间的夹角进行标定。
在一些实施例中,所述获取N段扫描投影数据,包括:对物体进行N段扫描,获取N段扫描投影数据p
I(l,t)~p
N(l,t),其中p
N(l,t)表示进行第N段扫描过程中,直线光源阵列中坐标为l的光源点闪烁时,探测器上坐标为t处的单元采集到的投影数据值。
在一些实施例中,所述根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,包括:
对所述N段扫描投影数据p
I(l,t)~p
N(l,t)进行平移变换,得到N段平移后的投影数据q
I(l,t)~q
N(l,t);
对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据Q
I(l,t′
I)~Q
N(l,t′
N),t′
N表示进行第N段扫描过程中,在直线光源阵列中坐标为l的光源点闪烁时,经过点(x,y)的射线在等效探测器上的采样位置;
对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像f
I(x,y)~f
N(x,y),(x,y)表示以物体中心为原点所建立的坐标系坐标;
对所述N段扫描投影数据的重建图像进行加和,得到最终重建图像f(x,y)。
在一些实施例中,对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据,包括:
分别将所述N段平移后的投影数据q
I(l,t)~q
N(l,t)乘上权重因子,用卷积核h沿着直线分布式探测器所在的直线方向做一维卷积得到N段加权 滤波后的数据Q
I(l,t′
I)~Q
N(l,t′
N)。
在一些实施例中,对所述N段平移后的投影数据进行加权处理,以消除多段投影数据冗余重叠所带来的重建伪影。
在一些实施例中,对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像f
I(x,y)~f
N(x,y),包括:
沿直线分布式光源所在的直线方向反投影Q
I(l,t′
I)~Q
N(l,t′
N),得到所述N段扫描投影数据的重建图像f
I(x,y)~f
N(x,y)。
(三)有益效果
从上述技术方案可以看出,本发明一种分布式光源CT图像重建方法与系统至少具有以下有益效果其中之一:
(1)本发明不需要将投影数据重排为平行束,具有更高的空间分辨率。对于直线分布式光源和探测器有限长度导致的有限角度问题,本发明利用两段光源和探测器的扫描模式(也可以是多段直线分布式光源和探测器,与之类似),有效解决了有限长度导致的有限角度问题,并且针对多段投影数据的重叠问题,通过投影数据加权,有效解决了多段投影数据冗余重叠问题。在多段直线分布式光源和探测器的静态CT系统中,本发明系统及方法能够获得精确的物体重建图像。
(2)本发明采用直线分布式光源和探测器,并且是静态扫描模式,有着快速的扫描速度,同时针对本发明系统提出的解析重建算法相对迭代重建方法有着很大的速度优势,所以本发明在速度要求较高的领域有着很大的潜能。
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为依据本发明实施例直线分布式光源和探测器的静态CT成像系 统。
图2依据本发明实施例两段直线分布式光源和探测器的扫描几何示意图。
图3依据本发明实施例重建算法流程示意图。
图4依据本发明实施例两段扫描投影的重叠分析示意图。
图5依据本发明实施例直线分布式光源和探测器的静态CT模拟重建成像示意图。
<符号说明>
11,12-直线分布式光源,21,22-直线分布式探测器,3待测物体。
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明提出一种分布式光源CT图像重建系统,包括:N段直线分布式光源,用于发出射线;以及N段直线分布式探测器,用于获取N段扫描投影数据,N≥2;其中,所述分布式光源CT图像重建系统利用所述N段扫描投影数据获取重建图像。
与传统的CT相比,本系统采用直线分布式光源和探测器阵列,并且在成像时,物体不需要旋转,同时直线光源阵列和直线探测器阵列是静止的,避免了滑环的使用,不但减少了噪音和振动,也降低了X射线机、高压模块、探测器等部件的机械强度要求,使整体系统更为安全、可靠,同时具备较快的扫描速度。此种系统具有应用于成像速度要求较高领域的潜力。
本发明还提供了一种分布式光源CT图像重建方法,包括:
获取N段扫描投影数据;以及
根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,其中,N≥2。
其中,所述获取N段扫描投影数据,包括:利用N段直线分布式光源及N段直线分布式探测器采用静态扫描模式获取N段扫描投影数据
进一步的,在获取N段扫描投影数据之前,还包括:对直线分布式光源之间的间距及夹角,直线分布式光源与直线分布式探测器之间的间距, 直线分布式探测器之间夹角进行标定。
本发明提出的基于直线分布式光源和探测器阵列的静态CT成像系统及其对应的滤波反投影重建算法,通过N段直线分布式光源和探测器,有效解决了投影数据角度缺失问题。通过加权方法,很好地解决了投影数据之间的数据冗余重叠问题,能够实现精确的物体重建。
为便于描述,本发明实施例以两段直线分布式光源和探测器为例进行介绍。
在一具体实施例中,为了解决投影数据角度缺失,采用两段直线分布式光源和探测器(多段直线分布式光源和探测器也是可以的,本实施例选择两段),如图1所示,所述两段直线分布式光源是互相垂直的,所述两段直线分布式探测器也是互相垂直的,第二段光源和探测器获得的投影等价于将物体旋转-90度,这样,两段光源和探测器能确保待成像物体的每个点都有至少180度的投影角度覆盖范围。直线分布式光源(X射线源),根据物体尺寸和应用背景而定。直线分布式探测器位于射线源对面。
对于N(N>2)段的情况,各相邻探测器之间的夹角及各相邻光源之间的夹角可以相等,也可以不相等,每段直线分布式光源和直线分布式探测器在实际空间中的角度覆盖范围之和能够无缝覆盖至少180度,以利于更加精确的物体图像重建。
在另一实施例中,本发明提出一种基于上述实施例直线分布式光源和探测器的静态CT系统的重建算法。将待测物体放在直线分布式光源和直线分布式探测器之间,在数据采集过程中,待测物体保持不动,如图1中(a)所示,第一段扫描时,光源阵列上的光源点依次闪烁,探测器阵列持续采集投影数据,当最后一个光源点闪烁结束后,进行第二段扫描,第二段扫描的光源探测器排布如图1中(b)所示,光源点如第一段扫描时一样依次闪烁。将获取的两段投影数据,利用本重建算法,就可以获得精确的重建图像。
以下详细介绍本发明系统结构及方法流程。
1.基于直线分布式光源和探测器的重建算法——滤波反投影(fltered backprojection)
此处先介绍无限长直线分布式光源和探测器的情况以利于理解。成像 系统几何结构如图2所示,当直线分布式光源和探测器无限长时,图像中每一点的重建值可以表示为:
其中,
这里直线分布式光源和探测器都是等距排列的,(x,y)表示以物体中心为原点所建立的坐标系坐标,数据p(l,t)表示光源阵列中为坐标为l的光源点闪烁时,探测器阵列中坐标为t的单元采集到的投影值,值得注意的是,对于p(l,t)来说,l和t分别以光源和探测器阵列的中心为坐标原点。而数据q(l,t)是光源阵列中坐标为l的光源点闪烁时,等效探测器上坐标为t处的数值(等效探测器是将直线分布式探测器平移到物体中心,以物体中心为坐标原点)。D和D′分别表示光源阵列和探测器阵列到物体中心的距离。
由于实际应用中,直线分布式光源和探测器都是有限长度的,所以按照以上的方法进行重建会带来有限角度的问题,即待成像物体中部分点并不能满足180度投影数据的要求。为了解决有限角度的问题,确保让待成像物体的每个点都至少有180度的投影覆盖范围,本发明系统共有两段直线分布式光源和探测器,这两段在空间中是互相垂直的,如图1所示,即第二个扫描段获得的投影数据,相当于在第一扫描段中将待成像物体旋转-90度。这样,两个扫描段的投影数据满足了待成像物体的每一个点都至 少有180度投影数据的要求,但这同时也带来投影数据重叠的问题,即待成像物体中的某些点会被相同的射线照射两次或以上,对应这个问题本发明又提出了相应的投影加权方法。本发明中,对于两段有限长的直线分布式光源和探测器,重建公式为:
f(x,y)=f
I(x,y)+f
II(x,y) (6)
其中,
其中q
I(l,t)与q
II(l,t)分别是原始采集的投影p
I(l,t)和p
II(l,t)的平移变换,而p
I(l,t)和p
II(l,t)分别是第一段扫描和第二段扫描采集得到的原始投影数据,W
I(l,t)和W
II(l,t)是为了解决两段投影数据冗余重叠问题而提出的权重因子,分别对应于第一段投影数据和第二段投影数据。
在此基础上,如图3所示,本实施例重建算法包括:
a,分别对于原始投影p
I(l,t)和p
II(l,t)进行平移变换,得到q
I(l,t)与q
II(l,t);
b,分别将投影q
I(l,t)与q
II(l,t)乘上
和
用卷积核h沿着t方向(第一段直线分布式探测器的分布方向和第二段直线分布式探测器的分布方向,也即各段直线分布式探测器所在的直线方向)做一维卷积得到Q
I(l,t′
I)和Q
II(l,t′
II);
c,沿l方向(第一段直线分布式光源的分布方向和第二段直线分布式光源的分布方向,也即各段直线分布式光源所在的直线方向)反投影Q
I(l,t′
I)和Q
II(l,t′
II),分别对应几何修正参数为
和
得到第一段扫描投影数据的重建图像f
I(x,y)和第二段扫描投影数据的重建图像f
II(x,y)
d,将两段重建结果f
I(x,y)和f
II(x,y)进行加和,得到最后的精确重建图像f(x,y)。
将原始投影数据p(l,t)进行平移变换得到q(l,t)的目的是为了得到等效探测器上的投影数据。用卷积核h滤波的目的与标准FBP重建算法中的滤波相同(算法推导出来需要滤波才能准确重建,同时也是为了消除仅投影产生的星条状伪影)。与重排为平行束的算法相比,本发明能充分利用每一个有效数据,提高了图像分辨率。需要说明的是,重建算法中的滤波的概念与图像处理中滤波的概念不同,它不代表传统意义上的滤除某个频率的信号,而是由于重建算法本身的需要而产生的一种变换。
2.滤波反投影算法的推导
a,投影数据重排为圆轨道平行束扫描
以第一段扫描几何为例,如图2中(a)所示,对于直线分布式光源和探测器的静态CT系统,直线扫描数据重排为圆轨道平行束扫描的重排公式为:
这里g(θ,s)=∫f(x,y)δ(xcosθ+ysinθ-s)dxdy,表示在圆轨道平行束扫描中,扫描角度为θ,距离旋转中心为s的投影数据。q(l,t)表示原始投影p(l,t)经过平移变换后得到的等效探测器上的投影数据,见公式(11),而p(l,t)的物理意义是:直线光源阵列位置l处的光源闪烁时,探测器阵列坐标位置为t的投影值。
利用公式(15)就可以实现扫描投影数据重排为圆轨道平行束扫描下的投影数据。虽然l和t的采样可以是均匀的,但是对应圆轨道平行束扫描下的角度θ和探测器位置s采样都是不均匀的。因此,重排需要在角度方向和探测器方向做插值,造成重建图像分辨率的降低。
b,直接滤波反投影重建
在圆轨道平行束扫描下的滤波反投影重建公式为
利用公式(15),用参数(l,t)替换(θ,s),得到
证明如下:
这样就得到了基于无限长直线分布式光源和无限长直线分布式探测器的静态CT的重建公式:
实际上,光源阵列和探测阵列都是有限长度的,而仅用一段光源和探测器得到的扫描数据会存在有限角度的问题,所以本系统中采用两段直线分布式光源和探测器进行数据获取,并对于投影进行加权以消除两段投影重叠而带来的重建伪影。
对于第一段有限长度直线分布式光源和探测器,其扫描集合如图2中(a)所示,与公式(20)的推倒过程类似,可以得到
其中,l
m和t
m分别为直线分布式光源和探测器的半长,对于第二段有限长度直线分布式光源和探测器,如图2中(b)所示,其扫描几何与圆轨道平行束扫描集合的关系与第一段相比会有一些差异,为
最后得到两段有限长度直线分布式光源和探测器的静态CT重建公示:
f(x,y)=f
I(x,y)+f
II(x,y) (24)
注意,W
I(l,t)和W
II(l,t)是本发明为了解决两段投影数据冗余重叠问题而提出的权重因子,分别对应于第一段投影数据和第二段投影数据,具体细节将于下一部分介绍。
3.权重因子
由于直线分布式光源和探测器都是有限长度的,若仅用一段的投影数据进行物体重建会带来有限角度的问题,故本发明中采用两段直线分布式光源和探测器扫描得到的投影数据进行物体的重建,这两段投影会存在投影的重叠情况如图4所示,多余的投影数据(即超过180度范围的投影数据)区域是H,K,R,F,分别与区域H’,K’,R’,F’重叠,为了让待重建物体中的每一点在180度的某一角度内有且只有一条射线照射,这里提出了对投影进行加权的方法,如公式中(9)和(10)中的W
I和W
II分别对应第一段扫描和第二段扫描投影的权重因子。
和区域H重叠的H’应该满足如下条件:θ
1=θ
2+π和s
1=-s
2,利用公式(15)和公式(22)条件转化为得到:
进一步化简可以得到区域H’的表达式为:
区域F’和区域F的重叠方式一致,同理,可以得到区域F’的表达式为:
和区域K重叠的K’满足的条件是:θ
1=θ
2和s
1=s
2,利用公式(15)和公式(22)代入(l
1,t
1)和(l
2,t
2),满足的条件转化为:
进一步化简可以得到区域K’的表达式为:
区域R’和区域R的重叠方式与K’和K的重叠方式一致,同理,可以得到区域R’的表达式为:
为了保证待重建物体的每一个点都有180度的角度覆盖范围,并且在180度的覆盖范围内的每个角度有且只有一条射线照射,权重因子W
I和W
II应该满足如下条件:
同时为了防止加权滤波后反投影重建图像出现明显的条纹伪影W
I和W
II在边界处还应该保持连续。
在第二段扫描中,如图4中(b)所示,区域分割线有5条,Line1~Line5,设这五条线的方程分别为:l
2=k
nt
2+b
n,n=1,2,3,4,5。由于第二段扫描提供
范围的投影数据,并且根据公式(22)可知,k
1=1,b
1=D。Line1与直线t
2=-t
m的交点为P1,Line1与Line2的交点为P2,Line2与Line3 的交点为P3,Line3与Line1的交点为P4,Line4与直线l
2=-l
m的交点为P5,Line4与Line5的交点为P6,Line5与直线t
2=t
m的交点为P7。设点P1~P7的坐标为(t
2-n,l
2-n),n=1,2,3,4,5,6,7。根据区域之间的对应关系可以计算得到这些点的坐标,如下表所示:
两点确定一条直线,根据点P1~P7的坐标,可以计算得到Line2~Line4的方程。
满足公式(35)和公式(36)以及连续性条件的W
I和W
II权重因子的形式有很多种表达,这里提供一种比较合适的加权方法:
其中函数fs(x)的表达式如下:
可由公式(36),计算得到W
I对应的表达式。
采用本发明直线分布式光源和探测器的静态CT可以实现扫描成像,重建断层图像的一典型实例如图5所示。本实例中,光源阵列和探测器阵列中都是等间距排列的,间距为2.5mm,光源到探测器的距离为2000mm,Shepp-Logan头模型由系统的设置根据解析方法获得前向投影,然后用提出的算法进行重建,重建像素大小为1.5mm
2。图5中(a)Shepp-Logan头模型原始模拟图,灰度窗[1,1.04];图5中(b)只由第一段扫描投影获得的重建图像,灰度窗[0.6,1.0];图5中(c)只有第二段扫描投影获得的重建图像,灰度窗[0.4,0.8];图5中(d)由两段扫描投影获得的重建图像,但是未采用权重因子处理投影重叠问题,灰度窗[1,1.6];图5中(e)由两段扫描投影获得的重建图像,采用了权重因子处理投影的重叠问题,灰度窗[1,1.04];图5中(f)重建图像(a)和(e)中心竖直剖面图。
当然,本发明光源和探测器的段数并不限于两段,多段即可,大于两段的情况与两段类似,此处不再赘述。
为到达精确的图像重建,X射线成像系统的应能够精确测量或标定以下系统参数:直线分布式X射线光源单元之间的间距,光源和探测器之间的间距,两段扫描之间的夹角,探测器物理尺寸,包括单个探测器物理尺寸和探测器阵列的物理尺寸等。
此外,整个成像系统的控制、数据传输、图像重建以及数据处理由计算机工作站完成,扫描控制信息、位置信息、投影数据等通过数据采集系统输入到计算机工作站中,由工作站完成物体的重建工作,最后在显示器上显示出来。
至此,已经结合附图对本实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明有了清楚的认识。
需要说明的是,在附图或说明书正文中,未绘示或描述的实现方式,均为所属技术领域中普通技术人员所知的形式,并未进行详细说明。
还需要说明的是,本文可提供包含特定值的参数的示范,但这些参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应值。实施例中提到的方向用语,仅是参考附图的方向,并非用来限制本发明的保护范围。此外,除非特别描述或必须依序发生的步骤,上述步骤的顺序并无限制于以上所列,且可根据所需设计而变化或重新安排。并且上述实施例可基于设计及可靠度的考虑,彼此混合搭配使用或与其他实施例混合搭配使用,即不同实施例中的技术特征可以自由组合形成更多的实施例。
应注意,贯穿附图,相同的元素由相同或相近的附图标记来表示。在以下描述中,一些具体实施例仅用于描述目的,而不应该理解为对本发明有任何限制,而只是本发明实施例的示例。在可能导致对本发明的理解造成混淆时,将省略常规结构或构造。应注意,图中各部件的形状和尺寸不反映真实大小和比例,而仅示意本发明实施例的内容。
当然,根据实际需要,本发明方法还包含其他的步骤,由于同本发明的创新之处无关,此处不再赘述。
说明书与权利要求中所使用的序数例如“第一”、“第二”、“第三”等的用词,以修饰相应的元件,其本身并不意味着该元件有任何的序数,也不代表某一元件与另一元件的顺序、或是制造方法上的顺序,该些序数的使用仅用来使具有某命名的一元件得以和另一具有相同命名的元件能做出清楚区分。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
- 一种分布式光源CT图像重建系统,其特征在于,包括:N段直线分布式光源,用于发出射线;以及N段直线分布式探测器,用于获取N段扫描投影数据,N≥2;其中,所述分布式光源CT图像重建系统利用所述N段扫描投影数据获取重建图像。
- 根据权利要求1所述的系统,其特征在于,所述系统包括2段直线分布式光源,分别为第一直线分布式光源和第二直线分布式光源、以及2段直线分布式探测器,分别为第一直线分布式探测器和第二直线分布式探测器;其中,第一直线分布式光源和第一直线分布式探测器分别位于待测物体两侧且相对设置,第二直线分布式光源和第二直线分布式探测器分别位于待测物体两侧且相对设置。
- 一种分布式光源CT图像重建方法,其特征在于,包括:获取N段扫描投影数据;以及根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,其中,N≥2。
- 根据权利要求3所述的方法,其特征在于,所述获取N段扫描投影数据,包括:利用N段直线分布式光源及N段直线分布式探测器采用静态扫描模式获取N段扫描投影数据。
- 根据权利要求4所述的方法,其特征在于,在获取N段扫描投影数据之前,还包括:对直线分布式光源之间的间距及夹角,直线分布式光源与直线分布式探测器之间的间距,直线分布式探测器之间的夹角进行标定。
- 根据权利要求4所述的方法,其特征在于,所述获取N段扫描投影数据,包括:对物体进行N段扫描,获取N段扫描投影数据p I(l,t)~p N(l,t),其中p N(l,t)表示进行第N段扫描过程中,直线光源阵列中坐标为l的光源点闪烁时,探测器上坐标为t处的单元采集到的投影数据值。
- 根据权利要求6所述的方法,其特征在于,所述根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,包括:对所述N段扫描投影数据p I(l,t)~p N(l,t)进行平移变换,得到N段平移后的投影数据q I(l,t)~q N(l,t);对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据Q I(l,t′ I)~Q N(l,t′ N),t′ N表示进行第N段扫描过程中,在直线光源阵列中坐标为l的光源点闪烁时,经过点(x,y)的射线在等效探测器上的采样位置;对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像f I(x,y)~f N(x,y),(x,y)表示以物体中心为原点所建立的坐标系坐标;对所述N段扫描投影数据的重建图像进行加和,得到最终重建图像f(x,y)。
- 根据权利要求7所述的方法,其特征在于,对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据,包括:分别将所述N段平移后的投影数据q I(l,t)~q N(l,t)乘上权重因子,用卷积核h沿着直线分布式探测器所在的直线方向做一维卷积得到N段加权滤波后的数据Q I(l,t′ I)~Q N(l,t′ N)。
- 根据权利要求8所述的方法,其特征在于,对所述N段平移后的投影数据进行加权处理,以消除多段投影数据冗余重叠所带来的重建伪影。
- 根据权利要求8所述的方法,其特征在于,对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像f I(x,y)~f N(x,y),包括:沿直线分布式光源所在的直线方向反投影Q I(l,t′ I)~Q N(l,t′ N),得到所述N段扫描投影数据的重建图像f I(x,y)~f N(x,y)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2019/082672 WO2020210941A1 (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2019/082672 WO2020210941A1 (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020210941A1 true WO2020210941A1 (zh) | 2020-10-22 |
Family
ID=72836919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2019/082672 WO2020210941A1 (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2020210941A1 (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
CN103065359A (zh) * | 2013-01-14 | 2013-04-24 | 厦门大学 | 一种光学成像的三维轮廓重建系统及其重建方法 |
CN104749648A (zh) * | 2013-12-27 | 2015-07-01 | 清华大学 | 多能谱静态ct设备 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
US20180235562A1 (en) * | 2017-02-17 | 2018-08-23 | Toshiba Medical Systems Corporation | Combined sinogram- and image-domain material decomposition for spectral computed tomography (ct) |
-
2019
- 2019-04-15 WO PCT/CN2019/082672 patent/WO2020210941A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
CN103065359A (zh) * | 2013-01-14 | 2013-04-24 | 厦门大学 | 一种光学成像的三维轮廓重建系统及其重建方法 |
CN104749648A (zh) * | 2013-12-27 | 2015-07-01 | 清华大学 | 多能谱静态ct设备 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
US20180235562A1 (en) * | 2017-02-17 | 2018-08-23 | Toshiba Medical Systems Corporation | Combined sinogram- and image-domain material decomposition for spectral computed tomography (ct) |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111265231B (zh) | 分布式光源ct图像重建方法与系统 | |
US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
Abella et al. | Software architecture for multi-bed FDK-based reconstruction in X-ray CT scanners | |
Stevens et al. | Alignment of a volumetric tomography system | |
JP4415762B2 (ja) | 断層撮影装置 | |
CN107764846B (zh) | 一种正交直线扫描的cl成像系统及分析方法 | |
CN107328798B (zh) | 一种新型icl系统及实现方法 | |
WO2007076681A1 (fr) | Systeme de scannage x-ct | |
US7209535B2 (en) | Fourier space tomographic image reconstruction method | |
US20120106832A1 (en) | Method and apparatus for ct image reconstruction | |
JP2008520257A (ja) | ボクセルに依存する補間を用いる画像再構成 | |
JP2001095794A (ja) | 斜めのスライスを使用するコンピュータ断層撮影法において画像を再構成するための装置および方法 | |
US7602879B2 (en) | Method for increasing the resolution of a CT image during image reconstruction | |
CN1722178A (zh) | 图像的三维重建方法和系统 | |
US8948337B2 (en) | Computed tomography image reconstruction | |
US20130251224A1 (en) | Method and system for reconstruction algorithm in cone beam ct with differentiation in one direction on detector | |
CN1605322A (zh) | 使扫描图像中的模糊最小化的方法和设备 | |
CN110520052B (zh) | 局部ct图像生成方法 | |
Sun et al. | Reduction of artifacts induced by misaligned geometry in cone-beam CT | |
Ametova et al. | Software-based compensation of instrument misalignments for X-ray computed tomography dimensional metrology | |
CN113533392B (zh) | 一种组合扫描cl成像方法 | |
CN111739113B (zh) | 直线分布式光源和探测器的ct图像重建方法及装置 | |
WO2020210941A1 (zh) | 分布式光源ct图像重建方法与系统 | |
CN112070704A (zh) | 一种基于紧小波框架的双正则化有限角ct图像重建方法 | |
Kong et al. | Analytic reconstruction approach for parallel translational computed tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19925515 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19925515 Country of ref document: EP Kind code of ref document: A1 |