CN111265231B - 分布式光源ct图像重建方法与系统 - Google Patents
分布式光源ct图像重建方法与系统 Download PDFInfo
- Publication number
- CN111265231B CN111265231B CN201910301697.9A CN201910301697A CN111265231B CN 111265231 B CN111265231 B CN 111265231B CN 201910301697 A CN201910301697 A CN 201910301697A CN 111265231 B CN111265231 B CN 111265231B
- Authority
- CN
- China
- Prior art keywords
- projection data
- light source
- scanning
- linear
- sections
- 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.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4007—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
- A61B6/4014—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
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段扫描投影数据pI(l,t)~pN(l,t),其中pN(l,t)表示进行第N段扫描过程中,直线光源阵列中坐标为l的光源点闪烁时,探测器上坐标为t处的单元采集到的投影数据值。
在一些实施例中,所述根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,包括:
对所述N段扫描投影数据pI(l,t)~pN(l,t)进行平移变换,得到N段平移后的投影数据qI(l,t)~qN(l,t);
对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据QI(l,t′I)~QN(l,t′N),t′N表示进行第N段扫描过程中,在直线光源阵列中坐标为l的光源点闪烁时,经过点(x,y)的射线在等效探测器上的采样位置;
对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像fI(x,y)~fN(x,y),(x,y)表示以物体中心为原点所建立的坐标系坐标;
对所述N段扫描投影数据的重建图像进行加和,得到最终重建图像f(x,y)。
在一些实施例中,对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据,包括:
分别将所述N段平移后的投影数据qI(l,t)~qN(l,t)乘上权重因子,用卷积核h沿着直线分布式探测器所在的直线方向做一维卷积得到N段加权滤波后的数据QI(l,t′I)~QN(l,t′N)。
在一些实施例中,对所述N段平移后的投影数据进行加权处理,以消除多段投影数据冗余重叠所带来的重建伪影。
在一些实施例中,对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像fI(x,y)~fN(x,y),包括:
沿直线分布式光源所在的直线方向反投影QI(l,t′I)~QN(l,t′N),得到所述N段扫描投影数据的重建图像fI(x,y)~fN(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.基于直线分布式光源和探测器的重建算法——滤波反投影(flteredbackprojection)
此处先介绍无限长直线分布式光源和探测器的情况以利于理解。成像系统几何结构如图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)=fI(x,y)+fII(x,y) (6)
其中,
其中qI(l,t)与qII(l,t)分别是原始采集的投影pI(l,t)和pII(l,t)的平移变换,而pI(l,t)和pII(l,t)分别是第一段扫描和第二段扫描采集得到的原始投影数据,WI(l,t)和WII(l,t)是为了解决两段投影数据冗余重叠问题而提出的权重因子,分别对应于第一段投影数据和第二段投影数据。
在此基础上,如图3所示,本实施例重建算法包括:
a,分别对于原始投影pI(l,t)和pII(l,t)进行平移变换,得到qI(l,t)与qII(l,t);
b,分别将投影qI(l,t)与qII(l,t)乘上和用卷积核h沿着t方向(第一段直线分布式探测器的分布方向和第二段直线分布式探测器的分布方向,也即各段直线分布式探测器所在的直线方向)做一维卷积得到QI(l,t′I)和QII(l,t′II);
c,沿l方向(第一段直线分布式光源的分布方向和第二段直线分布式光源的分布方向,也即各段直线分布式光源所在的直线方向)反投影QI(l,t′I)和QII(l,t′II),分别对应几何修正参数为和得到第一段扫描投影数据的重建图像fI(x,y)和第二段扫描投影数据的重建图像fII(x,y)
d,将两段重建结果fI(x,y)和fII(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)的推倒过程类似,可以得到
其中,lm和tm分别为直线分布式光源和探测器的半长,对于第二段有限长度直线分布式光源和探测器,如图2中(b)所示,其扫描几何与圆轨道平行束扫描集合的关系与第一段相比会有一些差异,为
最后得到两段有限长度直线分布式光源和探测器的静态CT重建公示:
f(x,y)=fI(x,y)+fII(x,y) (24)
注意,WI(l,t)和WII(l,t)是本发明为了解决两段投影数据冗余重叠问题而提出的权重因子,分别对应于第一段投影数据和第二段投影数据,具体细节将于下一部分介绍。
3.权重因子
由于直线分布式光源和探测器都是有限长度的,若仅用一段的投影数据进行物体重建会带来有限角度的问题,故本发明中采用两段直线分布式光源和探测器扫描得到的投影数据进行物体的重建,这两段投影会存在投影的重叠情况如图4所示,多余的投影数据(即超过180度范围的投影数据)区域是H,K,R,F,分别与区域H’,K’,R’,F’重叠,为了让待重建物体中的每一点在180度的某一角度内有且只有一条射线照射,这里提出了对投影进行加权的方法,如公式中(9)和(10)中的WI和WII分别对应第一段扫描和第二段扫描投影的权重因子。
和区域H重叠的H’应该满足如下条件:θ1=θ2+π和s1=-s2,利用公式(15)和公式(22)条件转化为得到:
进一步化简可以得到区域H’的表达式为:
区域F’和区域F的重叠方式一致,同理,可以得到区域F’的表达式为:
和区域K重叠的K’满足的条件是:θ1=θ2和s1=s2,利用公式(15)和公式(22)代入(l1,t1)和(l2,t2),满足的条件转化为:
进一步化简可以得到区域K’的表达式为:
区域R’和区域R的重叠方式与K’和K的重叠方式一致,同理,可以得到区域R’的表达式为:
为了保证待重建物体的每一个点都有180度的角度覆盖范围,并且在180度的覆盖范围内的每个角度有且只有一条射线照射,权重因子WI和WII应该满足如下条件:
同时为了防止加权滤波后反投影重建图像出现明显的条纹伪影WI和WII在边界处还应该保持连续。
在第二段扫描中,如图4中(b)所示,区域分割线有5条,Line1~Line5,设这五条线的方程分别为:l2=knt2+bn,n=1,2,3,4,5。由于第二段扫描提供范围的投影数据,并且根据公式(22)可知,k1=1,b1=D。Linel与直线t2=-tm的交点为P1,Linel与Line2的交点为P2,Line2与Line3的交点为P3,Line3与Line1的交点为P4,Line4与直线l2=-lm的交点为P5,Line4与Line5的交点为P6,Line5与直线t2=tm的交点为P7。设点P1~P7的坐标为(t2-n,l2-n),n=1,2,3,4,5,6,7。根据区域之间的对应关系可以计算得到这些点的坐标,如下表所示:
两点确定一条直线,根据点P1~P7的坐标,可以计算得到Line2~Line4的方程。
满足公式(35)和公式(36)以及连续性条件的WI和WII权重因子的形式有很多种表达,这里提供一种比较合适的加权方法:
其中函数fs(x)的表达式如下:
可由公式(36),计算得到WI对应的表达式。
采用本发明直线分布式光源和探测器的静态CT可以实现扫描成像,重建断层图像的一典型实例如图5所示。本实例中,光源阵列和探测器阵列中都是等间距排列的,间距为2.5mm,光源到探测器的距离为2000mm,Shepp-Logan头模型由系统的设置根据解析方法获得前向投影,然后用提出的算法进行重建,重建像素大小为1.5mm2。图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 (6)
1.一种分布式光源CT图像重建方法,其特征在于,包括:
获取N段扫描投影数据;以及
根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像,其中,N≥2,
所述获取N段扫描投影数据,包括:对物体进行N段扫描,获取N段扫描投影数据pI(l,t)~pN(l,t),其中pN(l,t)表示进行第N段扫描过程中,直线光源阵列中坐标为l的光源点闪烁时,探测器上坐标为t处的单元采集到的投影数据值,
对所述N段扫描投影数据pI(l,t)~pN(l,t)进行平移变换,得到N段平移后的投影数据qI(l,t)~qN(l,t),
对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据QI(l,t′I)~QN(l,t′N),t′N表示进行第N段扫描过程中,在直线光源阵列中坐标为l的光源点闪烁时,经过点(x,y)的射线在等效探测器上的采样位置,
对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像fI(x,y)~fN(x,y),(x,y)表示以物体中心为原点所建立的坐标系坐标,
对所述N段扫描投影数据的重建图像进行加和,得到最终重建图像f(x,y),
在N为2时,所述根据获取N段扫描投影数据,利用加权滤波反投影重建算法,获取重建图像的重建公式为:
f(x,y)=fI(x,y)+fII(x,y)
其中,
2.根据权利要求1所述的方法,其特征在于,所述获取N段扫描投影数据,包括:利用N段直线分布式光源及N段直线分布式探测器采用静态扫描模式获取N段扫描投影数据。
3.根据权利要求2所述的方法,其特征在于,在获取N段扫描投影数据之前,还包括:对直线分布式光源之间的间距及夹角,直线分布式光源与直线分布式探测器之间的间距,直线分布式探测器之间的夹角进行标定。
4.根据权利要求1所述的方法,其特征在于,对所述N段平移后的投影数据进行加权滤波处理,得到N段加权滤波后的数据,包括:
分别将所述N段平移后的投影数据qI(l,t)~qN(l,t)乘上权重因子,用卷积核h沿着直线分布式探测器所在的直线方向做一维卷积得到N段加权滤波后的数据QI(l,t′I)~QN(l,t′N)。
5.根据权利要求4所述的方法,其特征在于,对所述N段平移后的投影数据进行加权处理,以消除多段投影数据冗余重叠所带来的重建伪影。
6.根据权利要求4所述的方法,其特征在于,对所述N段加权滤波后的数据进行反投影,得到所述N段扫描投影数据的重建图像fI(x,y)~fN(x,y),包括:
沿直线分布式光源所在的直线方向反投影QI(l,t′I)~QN(l,t′N),得到所述N段扫描投影数据的重建图像fI(x,y)~fN(x,y)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910301697.9A CN111265231B (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910301697.9A CN111265231B (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111265231A CN111265231A (zh) | 2020-06-12 |
CN111265231B true CN111265231B (zh) | 2021-08-31 |
Family
ID=70991036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910301697.9A Active CN111265231B (zh) | 2019-04-15 | 2019-04-15 | 分布式光源ct图像重建方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111265231B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111839568B (zh) * | 2020-07-22 | 2023-12-12 | 重庆大学 | 一种新型大视场直线扫描ct系统及图像重建方法 |
CN111982939B (zh) * | 2020-07-27 | 2024-03-22 | 清华大学 | 可移动式多段直线光源ct成像系统及方法 |
CN113223112B (zh) * | 2021-06-22 | 2022-07-15 | 清华大学 | 广义等角探测器ct的图像解析重建方法 |
CN113533392B (zh) * | 2021-07-12 | 2022-08-26 | 重庆大学 | 一种组合扫描cl成像方法 |
CN114324425A (zh) * | 2022-01-04 | 2022-04-12 | 北京航星机器制造有限公司 | 一种多视角静态ct系统 |
WO2024046424A1 (zh) * | 2022-08-31 | 2024-03-07 | 同方威视技术股份有限公司 | 用于静态ct设备的成像方法、静态ct设备、电子设备和介质 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
CN103901488A (zh) * | 2012-12-27 | 2014-07-02 | 同方威视技术股份有限公司 | 固定式ct装置 |
CN103948395A (zh) * | 2007-07-19 | 2014-07-30 | 北卡罗来纳大学查珀尔希尔分校 | 固定 x 射线数字化断层合成或断层摄影系统和相关方法 |
CN104483711A (zh) * | 2014-12-17 | 2015-04-01 | 同方威视技术股份有限公司 | 基于分布式光源的辐射成像系统 |
CN204422777U (zh) * | 2014-12-17 | 2015-06-24 | 同方威视技术股份有限公司 | 基于分布式光源的辐射成像系统 |
CN104749197A (zh) * | 2013-12-26 | 2015-07-01 | 清华大学 | Ct系统及其方法 |
CN104768467A (zh) * | 2012-08-17 | 2015-07-08 | 北卡罗来纳大学教堂山分校 | 具有分布式x射线源阵列的固定台架计算机断层扫描系统和方法 |
CN105806858A (zh) * | 2014-12-31 | 2016-07-27 | 北京固鸿科技有限公司 | Ct检测方法和ct设备 |
CN106153648A (zh) * | 2016-06-22 | 2016-11-23 | 深圳先进技术研究院 | 一种静态高分辨显微ct成像系统及应用其的成像方法 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
CN106651982A (zh) * | 2016-12-16 | 2017-05-10 | 西安交通大学 | 一种基于阵列x射线源和探测器的ct图像重建方法 |
CN106896120A (zh) * | 2015-12-18 | 2017-06-27 | 清华大学 | 多模态检测系统和方法 |
CN108201447A (zh) * | 2017-11-21 | 2018-06-26 | 深圳先进技术研究院 | 一种静态ct系统几何参数校正方法 |
CN108577876A (zh) * | 2018-02-28 | 2018-09-28 | 西安交通大学 | 一种多边形静止ct及其工作方法 |
CN110428478A (zh) * | 2019-07-15 | 2019-11-08 | 清华大学 | 交替光源扇束x射线ct采样方法及装置 |
-
2019
- 2019-04-15 CN CN201910301697.9A patent/CN111265231B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103948395A (zh) * | 2007-07-19 | 2014-07-30 | 北卡罗来纳大学查珀尔希尔分校 | 固定 x 射线数字化断层合成或断层摄影系统和相关方法 |
CN102456227A (zh) * | 2010-10-28 | 2012-05-16 | 清华大学 | Ct图像重建方法及装置 |
CN104768467A (zh) * | 2012-08-17 | 2015-07-08 | 北卡罗来纳大学教堂山分校 | 具有分布式x射线源阵列的固定台架计算机断层扫描系统和方法 |
CN103901488A (zh) * | 2012-12-27 | 2014-07-02 | 同方威视技术股份有限公司 | 固定式ct装置 |
CN104749197A (zh) * | 2013-12-26 | 2015-07-01 | 清华大学 | Ct系统及其方法 |
CN104483711A (zh) * | 2014-12-17 | 2015-04-01 | 同方威视技术股份有限公司 | 基于分布式光源的辐射成像系统 |
CN204422777U (zh) * | 2014-12-17 | 2015-06-24 | 同方威视技术股份有限公司 | 基于分布式光源的辐射成像系统 |
CN105806858A (zh) * | 2014-12-31 | 2016-07-27 | 北京固鸿科技有限公司 | Ct检测方法和ct设备 |
CN106896120A (zh) * | 2015-12-18 | 2017-06-27 | 清华大学 | 多模态检测系统和方法 |
CN106153648A (zh) * | 2016-06-22 | 2016-11-23 | 深圳先进技术研究院 | 一种静态高分辨显微ct成像系统及应用其的成像方法 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
CN106651982A (zh) * | 2016-12-16 | 2017-05-10 | 西安交通大学 | 一种基于阵列x射线源和探测器的ct图像重建方法 |
CN108201447A (zh) * | 2017-11-21 | 2018-06-26 | 深圳先进技术研究院 | 一种静态ct系统几何参数校正方法 |
CN108577876A (zh) * | 2018-02-28 | 2018-09-28 | 西安交通大学 | 一种多边形静止ct及其工作方法 |
CN110428478A (zh) * | 2019-07-15 | 2019-11-08 | 清华大学 | 交替光源扇束x射线ct采样方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111265231A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111265231B (zh) | 分布式光源ct图像重建方法与系统 | |
US7215731B1 (en) | Fast backprojection/reprojection with hexagonal segmentation of image | |
US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
CN105596022B (zh) | 锥束ct动态成像的方法和设备 | |
US6256370B1 (en) | Method and apparatus for performing tomosynthesis | |
US8625870B2 (en) | Method and system for supplementing detail image in successive multi-scale reconstruction | |
US20120155728A1 (en) | Tomographic iterative reconstruction | |
JP2001515378A (ja) | らせん走査ctスキャナにおけるオンライン画像再構成 | |
CA2425323A1 (en) | System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit | |
JP2009529394A (ja) | ラドンデータから画像関数を再構成する方法 | |
EP1644897B1 (en) | A fourier tomographic image reconstruction method for fan-beam data | |
CN107822652B (zh) | 用于重建光谱结果图像数据的方法 | |
CN102456227A (zh) | Ct图像重建方法及装置 | |
US7602879B2 (en) | Method for increasing the resolution of a CT image during image reconstruction | |
US8948337B2 (en) | Computed tomography image reconstruction | |
JPH08263638A (ja) | 物体の断層写真像を作成するシステム及び方法 | |
US20160210762A1 (en) | Method for tomographic reconstruction | |
CN102062740B (zh) | 锥束ct扫描成像方法及系统 | |
CN111839568B (zh) | 一种新型大视场直线扫描ct系统及图像重建方法 | |
US8385620B2 (en) | Method and system for multi-detector-row computed tomography utilizing projection data up-sampling with shift | |
CN112070704A (zh) | 一种基于紧小波框架的双正则化有限角ct图像重建方法 | |
Costin et al. | A 2D multiresolution image reconstruction method in X-ray computed tomography | |
CN111583389A (zh) | 一种基于cad模型的不完全扫描ct图像重建方法 | |
JP6243219B2 (ja) | 画像生成装置および放射線断層撮影装置並びにプログラム | |
US4769757A (en) | Improved CT scanner having distortion-free image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |