CN100466981C - 图像重建方法和x射线计算机断层摄影设备 - Google Patents
图像重建方法和x射线计算机断层摄影设备 Download PDFInfo
- Publication number
- CN100466981C CN100466981C CNB2005101373058A CN200510137305A CN100466981C CN 100466981 C CN100466981 C CN 100466981C CN B2005101373058 A CNB2005101373058 A CN B2005101373058A CN 200510137305 A CN200510137305 A CN 200510137305A CN 100466981 C CN100466981 C CN 100466981C
- Authority
- CN
- China
- Prior art keywords
- data
- projection
- filtering
- weighting
- filter
- 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.)
- Expired - Fee Related
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
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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/027—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Abstract
一种在扫描对象的计算机断层摄影(CT)重建点确定图像数据值的图像重建方法,包括:获取(900)主体的投影数据,用一维斜坡滤波器滤波(901)该投影数据,以生成斜坡滤波数据,将距离倒数加权的反向投影算子作用(902)于该斜坡滤波数据,以在CT图像的重建点生成图像数据值。
Description
技术领域
本发明主要涉及医学图像的重建。更具体地,本发明涉及新的图像重建方法和X射线计算机断层摄影设备,以利用混合型滤波提高重建图像质量和重建效率。
背景技术
本发明包括各种技术的使用,在以下参考列表记载的文档中参考并描述,在整个说明书中由括号中相应的参考标号表示:
参考文献
[1]C.R.Crawford和K.F.King,Computerized tomography scanning withsimultaneous patient translation,Med.Phys.,17,967-982,1990年。
[2]P.E.Danielsson,P.Edholm,J.Eriksson和M.Seger,Towards exact3D-econstruction for helical cone-beam scanning of long objects。A new detectorarrangement and a new completeness condition,Proc.on Fully 3D ImageReconstruction in Radiology and Nuclear Med.,141-144,1997年。
[3]L.A.Feldkamp,L.C.Davis和J.W.Kress。Practical cone-beamalgorithm,J.Opt.Soc.Am,vol.1,pp.612-619,1984年。
[4]G.T.Herman和A.Naparstek,Fast image reconstruction based on aRadon inversion formula appropriate for rapidly collected data,SIAM J.Appl.Math.,vol.33,pp.511-533,1977年。
[5]A.Katsevich,Theoretically Exact FBP-Type Inversion Algorithm forSpiral CT,第六次无线电与核医学中的完整三维图像重建国际会议,PacificGrove,CA,USA,2001年10月-11月。
[6]A.Katsevich。An Improved Exact Filtered Backprojection Algorithm forSpiral Computed Tomography,Advances in Applied Mathematics,V.32-4,pp.625-825,2004年5月。
[7]A.Katsevich。Analysis of an Exact Inversion Algorithm for SpiralCone-Beam CT,Physics in Medicine and Biology,2002年,vol.47,pp.2583-2598。
[8]H.Kudo,F.Noo,M.Defrise和R.Clackdoyle,New super-short-scanreconstruction algorithm for fan-beam and cone-beam tomography,IEEENSS-MIC 2002年,M5-3。
[9]H.Kudo,F.Noo,M.Defrise和T.Rodet,New Approximate FilteredBackprojection Algorithm for Helical Cone-Beam CT with Redundant Data,IEEE NSS-MIC 2003年,M14-330。
[10]H.Kudo,T.Rodet,F.Noo和M.Defrise,Exact and approximatealgorithm for helical cone-beam ct with redundant data,Phys.Med.Biol.49(2004年),2525-2546。
[11]A.C.Kak和M.Slaney,Principles of computerized tomographicimaging,IEEE Press,纽约,1988年。
[12]A.Katsevich和A.Zamyatin,Analysis of a family of exact inversionformulas for cone beam CT,已提交。
[13]A.V.Lakshminarayanan,Reconstruction from divergent ray data,Tech.Rep.92,Dept.of Computer Science,布法罗,SUNY,1975年。
[14]F.Natterer,The mathematics of computerized tomography,纽约:Wiley,1986年。
[15]F.Noo,M.Defrise,R.Clackdoyle,H.Kudo,Image reconstruction fromfan-beam projections on less then a short scan,Phys.Med.Biol.,47,2525-2546,2002年。
[16]D.L.Parker,Optimal short scan convolution reconstruction forfan-beam CT,Med.Phys.9(2),1982年,254-257。
[17]J.D.Pack,F.Noo和H.Kudo,Investigation of saddle trajectories forcardiac CT imaging in cone-beam geometry,Phys.Med.Biol.,49,No11,2317-2336,2004年。
[18]M.D.Silver,A method for including redundant data in computedtomography,Med.Phys.,27,773-774,2000年。
[19]M.D.Silver,K.Taguchi和K.S.Han,Field-of-view dependent helicalpitch in multi-slice CT Proc.of SPIEMed.Imag.Conf.4320,839-850,2001年。
[20]K.Taguchi,Temporal resolution and the evaluation of candidatealgorithms for four-dimensional CT,Med.Phys.,30(4),640-650,2003年。
[21]K.Taguchi,B.S.Chiang和M.D.Silver,A new weighting scheme forcone-beam helical CT to reduce the image noise,Phys.Med.Biol.49,2351-2364,2004年。
[22]K.C.Tam,S.Samarasekera和F.Sauer,Exact cone-beam CT with aspiral scan,Phys.Med.Biol.,43,1015-1024,1998年。
[23]G.Wang,C.R.Crawford,W.A.Kalender,Multi-row-detector andcone-beam spiral/helical CT.IEEE Trans.Med.Imaging 19:817-821,2000年。
[24]G.Wang,T.H.Lin,P.Cheng和D.M.Shinozaki,A general cone-beamreconstruction algorithm,IEEE Trans.Med.Imaging,MI-12,486-96,1993年。
[25]A.Zamyatin,Analysis of Cone Beam Reconstruction in ComputerTomography,PhD Dissertation,中佛罗里达大学,2003年。
由计算机断层摄影(CT)装置创建的重建图像的质量和效率对该CT装置的总体性能来说是很重要的。重建图像时使用的算法影响质量和效率。
图1表示CT装置的一种可能的源轨迹的实例,其中可得到重建函数f(x),在二维空间中x=(x1,x2),在三维空间中x=(x1,x2,x3)。图1的x3轴用z轴表示。图1所示的螺旋轨迹由方程a(t)=(R cos t,R sin t,tH/2π)给出,其中H称为螺距(helical pitch)。对于一个图像切片的重建,仅需要有限的投影范围。通常这种投影范围,表示为Λ,集中在由t0=x3·2π/H限定的图像切片上。可根据s改写该轨迹方程(如a(s)=(R cos s,R sin s,sH/2π),s∈Λ)其中s=t-t0。有时投影角β(或λ)代替参数s使用。注意β=s mod 2π。图1所示的L(s,x)表示从重建点x到焦点a(s)的水平距离。
图1的虚线圆圈表示视野(FOV)。FOV为包含重建对象或病人的圆形区域。通常FOV由其半径rFOV限定。FOV直径的典型值是500mm。感兴趣区(theregion of interest)(ROI)可以是整个FOV或FOV的一部分。
检测器几何形状和数据结构
用g(s,Θ)表示在Θ方向上沿着来自源位置a(s)的射线的总衰减(或线积分)。该函数表示发生重建前通过CT扫描获得的数据。如图2所示,γ表示扇形角,α表示圆锥角。在扇形射束的情况下,Θ=γ。在圆锥射束的情况下,Θ=(γ,α)。也可参照该检测器表示Θ角。参照该检测器,ΘD=(u,v)(或在扇形射束的情况下ΘD=u),因此ΘD取决于检测器几何形状。实际上,圆锥射束数据是检测器坐标g(s,ΘD)=g(s,u,v)的函数,但一些理论公式将其作为扇形角和圆锥角g(s,Θ)=g(s,γ,α),即无检测器表达式(detector-free notation)的函数。这里使用无检测器表达式。
有许多具有不同几何形状的检测器类型:等角的、等距的、不等距的、平的、圆柱的、球的、倾斜的、旋转的、PI掩蔽的(PI-masked)等等。图2表示等角的圆柱形检测器,图3表示平的等距共线检测器。
卷积
为清楚起见,下面使用算子符号。每个算子、卷积或反向投影大体上是独立于其自变量(argument)而定义的。算子用粗斜体表示,以区别于函数。相关算子定义如下。
傅立叶变换表示为
G(s,ω)=FT[g(s,γ)],
其中傅立叶变换应用于第二变量γ。
希耳伯特滤波算子(H)表示为
H[g(s,γ)]=g(s,γ)*h(γ)
H[g(s,γ)]=FT-1[G(s,ω)H(ω)],其中
h(γ)=-1/(πγ),和H(ω)=i sign ω。
斜坡滤波算子(Q)表示为
Q[g(s,γ)]=g(s,γ)*q(γ)
Q[g(s,γ)]=FT-1[G(s,ω)Q(ω)],其中
q(γ)=FT-1[Q(ω)],和Q(ω)=|ω|。
改进的希耳伯特滤波算子(Hm)(“改进的”表示核函数是h(sinγ)而不是h(γ))表示为
Hm[g(s,γ)]=g(s,γ)*h(sinγ)。
改进的斜坡滤波算子(Qm)表示为
Qm[g(s,γ)]=g(s,γ)*q(sinγ),
其中 以及Qm(ω)=FT[q(sinγ)]。
最后,伴有直流漂移(Qm0)的改进的斜坡滤波算子表示为
Qm0[g(s,γ)]=FT-1[G(s,ω)Qm(ω)+G(s,0)/4π]。
斜坡滤波一般用于滤波的反向投影(FBP)算法。FBP是计算机断层摄影(CT)中主要的重建方法。与希耳伯特滤波相比,具有如下优点:(1)对于同一噪声级,具有更高的分辨率(更清晰的图象),(2)为了在重建中平衡噪声,对分辨率进行更多的控制。第二个优点对于临床应用是非常重要的。在本发明中,主要贡献是从斜坡滤波的数据重建图像。斜坡滤波保持图像的所有高频分量。希耳伯特滤波为低频校正,使圆锥射束赝像(artifact)减少。
反向投影
一般地,反向投影算子作用于数据g(s,Θ),其中该算子将该数据映射到图像空间。由如下方程表示:
这里,首先在Θ=Θ(s,x)计算g(s,Θ),其中Θ(s,x)是来自横跨x的a(s)的射线。然后该方程沿着投影范围Λ上的源轨迹对所有这种射线进行积分。为了找到Θ(s,x)使用如下方程:
其中za(s)是a(s)的z坐标。反向投影通常应用于滤波的数据。
有两种方法加权该反向投影算子。该反向投影算子可由L的倒数或L2的倒数加权。这两种方法在如下两个方程中表示。
与1/L相比,1/L2的反向投影加权导致噪声恶化和图像上的点展开函数(BSF)的均匀性。
该反向投影范围由Λ表示。用BPJ算子符号中的上角标表示反向投影范围如下:
图4A-4E表示五个不同的重建范围。图4A表示全扫描,图4B表示过扫描,图4C表示短扫描,图4D表示超短扫描,图4E表示一种弹隍范围。
冗余补偿加权
可使用加权因子w(s,γ)修正数据中的冗余。对扇形射束数据不必在整个区域上扫描。对于扇形射束几何结构,g(s,γ)=g(s+π+2γ,-γ)。对于扇形射束数据,整个扫描会导致对毒千射线计数两次。因此,该数据需要由函数w(s,γ)=1/2加权。当跨过FOV的任何线与重建段相交至少一次时,满足二维数据充分条件[14]。当重建范围Λ≥π+2arcsin(rFOV/R)时,满足二维数据充分杂件。对扇形射柬数据来说,覆盖Λ=π+2arcsin(rFOV/R)的投影集称为“最小完整数据集”[16]。如下术语在现有技术中已经使用并在这里使用:
短扫描,或半扫描重建:Λ=π+2Г,
arcsin(rFOV/R)≤Г<π/2;
超短扫描重建:Λ<π+2arcsin(rFOV/R);
全扫描重建:Λ=2π;和
过扫描重建:Λ>2π。
短扫描的Parker加权
从关系g(s,γ)=g(s+π+2γ,-γ),其中-γm≤γ≤γm可以看到,对于精确扇形射束重建,仅π+2γm的重建范围已经足够,且对于区域π-2γ<s<π+2γm中的数据,数据g(s,γ),0<s<2γm,-2γ是冗余的。这里γm=arcsin(rFOV/R)为该检测器允许的最大扇形角。参考[16]指出,在Parker加权中,应该以这样的方式加权最小完整数据区的数据,即尽可能对间隔均匀分配。参考[16]提出了如下加权函数:
wP(s,γ)+wP(s+π+2γ,-γ)=1
γm=arcsin(rFOV/R)为该检测器允许的最大扇形角。如果没有定义,则该加权为零,例如,如果s<0或s>π+2γm,则wP(s,γ)=0。在BPJ算子符号中用上角标P表示反向投影范围π+2γm:
通用Parker加权(MHS加权)
Parker加权是MHS加权的一种特殊情况。在MHS加权中,γm应至少为arcsin(rFOV/R)。通过实质上增大的γm,可获得更大的重建范围,从而获得更佳的噪声特性[18,19]。在wP(β,γ)的方程中,通过用虚拟最大扇形角Γ代替物理最大扇形角γm,得到如下方程:
在BPJ算子符号中用上角标MHS表示反向投影范围π+2Γ:
过扫描加权
由于g(s,γ)=g(s+2π,γ),过扫描加权用于反向投影范围Λ=2πn+Δ,其中n=1,2,......,0<Δ<2π。加权函数如下。
在过扫描加权中,可使用MHS加权和Parker加权代替三角函数sin2或cos2、多项式3x2-2x3[1]或某些其它光滑函数。在BPJ算子符号中,用上角标OS表示反向投影范围2πn+Δ:
Noo的加权
Noo的加权的优点在于可使用任意重建范围Λ(s0,s1),其中s0和s1为该重建段的起点和终点。这种加权具有小于半扫描的重建范围,可用于ROI重建。这种重建范围称为短扫描。该加权函数如下:
例如,在简单的短扫描情况下,
在过扫描情况下:
函数c(β)如下:
其中Δs为平滑滤波间隔。Noo建议取Δs=10°[15]。注意具有大Δs(≈50°)的Noo的加权等价于Parker加权。Noo的加权可使用任意反向投影范围Λ,如图4A-4E所示[15]。
准圆锥射束加权
扇形射束加权可扩展至圆锥射束数据[21]。一旦计算出扇形射束加权w(β,γ),作为圆锥角的函数,对其加权,并对其标准化,以得到准圆锥射束加权函数wQ3D(β,γ,α)。必须基于该数据的有效性(精确性)定义这种加权,有效性加权由
wVal(α)=1-α(3x2-2x3)表示,其中
且α1=tan-1[t1D/(2R)];α2=tan-1[t2D/(2R)]。
两个圆锥角(α1和α2)定义有效曲线的拐点。为了补偿冗余采样,有效加权wVal(α)与扇形射束加权函数w(β,γ)相结合,然后标准化。有效加权wVal(α)可以是任意值;但是,如果选择参数t1和t2,则wVal(α)是有意义的,从而wVal(α)将全加权分配给有效(测量的)射线和,将较少的加权分配给无效(未测量的)射线和,并对该转换滤波。因此,准圆锥射束加权函数为:
其中在所有补偿位置求和,以使
Tam窗口加权
Tam窗口是由来自焦点a(s)的螺旋线的上部和下部投影限定的检测器的一部分。图5表示仅包含非冗余数据的完整的πTam窗口[22,2]。图5表示rH=1.0的情况。虚线表示当rH=1.25时物理检测器的边界。这里rH=H/D,其中D为检测器宽度。π Tam窗口的加权函数如下:
3π Tam窗口的加权函数如下:
在BPJ算子符号中,用上角标π仅表示π Tam窗口的反向投影,根据螺距,在BPJ算子符号中,用上角标TW表示π Tam窗口或3π Tam窗口的反向投影:
注意Parker和Noo的加权函数最初都是为扇形射束数据引入的和圆锥角的独立。这些加权函数的使用产生圆锥射束的赝像。Tam窗口加权为移位常量或呈投影角独立提供三维加权,即同一加权用于每个投影。
Tam窗口加权的优点是:(1)真实的圆锥射束加权,(2)移位常量,即加权函数wT(γ,v)对所有投影来说是相同的,独立于位置z,(3)执行简单。Tam窗口加权的缺点是:(1)未使用冗余数据(未使用测量的数据的某部分),(2)Tam窗口是固定的(因此对应于π和3π,仅两个螺距是最佳的),(3)利用不同的反向投影范围,重建不同的图像象素,产生较少的空间均匀图像。
在图6中,虚线表示当该源沿该螺旋线移动时,两个不同图像点沿该检测器的路径。粗线表示在反向投影中使用的数据,以重建A点和B点。这样,如从图6所见,A点的重建比B点的重建使用更多的数据。
扩展的Tam窗口加权
利用平滑滤波函数,Tam窗口可在z方向上扩展。该扩展的Tam窗口在图7中表示。扩展的Tam窗口加权可失去理论精确度,但具有更高的数据使用率。在BPJ算子符号中,用上角标ETW表示扩展的Tam窗口加权:
滤波线的方向
图8A-8F表示不同的滤波线。起初,沿检测器行或段(图8A)滤波。对于螺旋轨迹,沿着平行于源位置点的螺旋线的切线的方向(切线)滤波,使圆锥射束赝像减少。由HTan和QTan表示的算子沿切线(图8B)滤波。在旋转滤波中(图8C),中间滤波线与该螺旋线相切,最顶部和最底部的滤波线是水平的(平的)。它们之间的其它滤波线逐步旋转,以形成平滑滤波线族。由HRot和QRot表示的算子沿旋转线滤波。精确的算法[12]使用Katsevich的滤波线族。可在图8D-8F中看到Katsevich滤波线,沿Katsevich的曲线族滤波表示为HKat和QKat。
在实际操作中,如果检测器行不平行于滤波方向,则滤波需要从检测器点阵变换(rebinning)到滤波器点阵。每个变换包括内插,产生更光滑的图像并可能失去细节。HTan和QTan(及其它)算子的使用使圆锥射束赝像减少,其代价是重建时间增加且分辨率减小。
加权和卷积的顺序
在操作中,执行加权和卷积的顺序是非常重要的。冗余加权w(s(x3),γ)取决于重建切片z位置x3。如果卷积前必须加权,即滤波应用于加权数据子集,则对于每个切片,重建范围中的所有投影需要重加权和重卷积。另一方面,如果加权在卷积之后,则对于所有图像切片,该数据可以仅卷积一次,且每个切片仅需重加权。因此,在后一种情况下,重建所需的卷积数减少了很多。
重建算法
斜坡滤波一般用于重建医学图像的FBP算法。扇形射束数据的原始FBP算法具有如下形式:
FBP算法首先在等距共线检测器的情况下发展起来[13],后来在等角射线[11]的情况下得以扩展[4]。它使用1/L2加权的斜坡滤波。可以看到1/L2的反向投影加权与1/L相比,导致噪声恶化和整个图像上的PSF均匀性。该原始算法起初是为全扫描轨迹设计的,后来扩展至短扫描轨迹[16]:
再后来[3]是为圆形圆锥射束的几何结构设计的:
(cosΘ=cosγcosα)。
卷积前应用短扫描轨迹的加权函数[FBP-P],即将滤波应用于加权的数据子集。由此,在该重建范围中所有投影的每个切片需要重加权并重卷积。卷积后进行加权的算法是更有效的。这是因为对所有数据图像切片来说,该数据仅卷积一次,然后重加权。因此,如果卷积后加权,重建所需的卷积数会减少很多。因此,卷积后应用加权函数的算法可提供显著的计算优点。
Feldkamp算法扩展至螺旋轨迹[24],后来扩展至任意扫描段[18]如下:
该算法的流程图如图13A所示。
上述讨论的算法限于圆形或螺旋源轨迹。不限于圆形或螺旋轨迹的算法是优选的。由于它可应用于其它轨迹如鞍形轨迹,因此可提供更通用的算法。
Katsevich[5,7]引入了FBP类型的精确圆锥射束算法。
它使用该圆锥射束数据的偏微分的改进的希耳伯特变换,而不使用一般的斜坡滤波。该算法的流程图如图13C所示。后来Katsevich将其方程概括为一个,不带项,并也发展了3π情况下的算法。在上述算法中,在滤波线的特殊族中必须进行滤波[12],在卷积前后,需要附加的非一般的(non-trivial)变换步骤。图8D-8F表示Katsevich滤波线族。然而出于实际目的,通用的Feldkamp算法相当适用,且具有相对较少检测器行的现有的扫描仪不需要这种复杂性。但是基于希耳伯特滤波的重建提供了一些Feldkamp算法所不具备的非常好的特性[3]。Katsevich的算法与Feldkamp算法相比,显示出基于希耳伯特滤波的重建提供了一些非常好的特性。
后来引入基于希耳伯特变换重建的扇形射束重建算法[15]:
参考[15]指出对于精确重建,仅需覆盖Λ=π的投影。这使得小于短扫描的重建成为可能。对于短扫描重建,需要整个FOV上的投影集(π+2γm)。另外,该数据充分性条件放松。[NDCK-FB]算法的其它优点是卷积后执行冗余扇形射束数据的加权。由于每个切片的数据不必重卷积,这使得[NDCK-FB]算法与任何Feldkamp类型的算法相比更加有效。
在发明人估算噪声流体圆柱形模型的过程中,发现了Noo的算法的另一个优点。可以证明,与Feldkamp算法相比,噪声变化在整个图像上更加均匀。对于Noo的算法,PSF也具有较少的空间变化。这可以解释为反向投影加权是距离的倒数,而不是距离平方的倒数,且减小了所谓放大作用。
基于希耳伯特变换的算法与利用斜坡滤波[8]的算法相比,由于附加的数值微分步骤,引入了更强的平滑滤波。Kudo[8]提出的扇形和圆锥射束数据的算法包括斜坡和希耳伯特过滤两个滤波:
[KNDC]
在参考[10]中,将[NDCK-FB]和[KNDC-FB]两算法概括为具有平的检测器的圆锥射束的情况如下:
[KRND-2]
Kudo的算法有几个缺点。首先,其平方倒数加权和仅对特殊加权函数起作用的缺点类似于Feldkamp的算法。其次,Kudo的算法包括对加权函数取偏微分,这使得该算法缺乏实用性。除了Katsevich的算法,上述算法对圆形轨迹上的扇形射束扫描是精确的,对圆锥射束扫描是近似的。
公知的重建算法有几个缺点。上述讨论算法中没有一个具有所有下列有益的方面:(1)1/L反向投影加权,(2)卷积后应用冗余补偿加权函数,(3)独立于所使用的加权函数的类型的算法,(4)独立于源轨迹的算法,以及(5)混合型过滤滤波的结合(斜坡型滤波和希耳伯特型滤波)。
发明内容
因此,为了克服相关技术中重建算法的问题,本发明寻求提供一种方法、系统和计算机程序产品,以在扫描对象的计算机断层摄影重建点确定图像数据值。
因此提供了一种方法、系统和计算机程序产品,以在扫描对象的计算机断层摄影重建点确定图像数据值,包括:(1)获取扫描对象的投影数据;(2)用一维斜坡型滤波器滤波获取的投影数据,以生成斜坡滤波数据;(3)将具有距离倒数加权的反向投影算子作用于该斜坡滤波数据,以在CT图像的重建点生成图像数据值。
另外,按照本发明的一个实施例,上述方法还包括:将投影减法作用于获取的投影数据,以生成减法数据;对该减法数据应用希耳伯特型滤波器,以生成希耳伯特滤波数据;将投影加法作用于希耳伯特滤波数据和斜坡滤波数据,以生成滤波的数据;将冗余加权作用于滤波的数据,以生成加权数据,其中应用该反向投影算子的步骤作用于该加权数据,以在CT图像的重建点生成图像数据值。
按照本发明的另一方面,提供了一种X射线计算机断层摄影(CT)系统,以在重建点确定图像数据值,包括:(1)CT扫描单元,用于生成扫描对象的投影数据,该扫描单元包括生成X射线的X射线源和具有检测器元件的检测器,该检测器用于产生该投影数据;以及(2)处理器,包括:滤波单元,用于对该投影数据应用斜坡型滤波器,以生成斜坡滤波数据;以及反向投影单元,用于对该斜坡滤波数据应用具有距离倒数加权的反向投影算子,以在重建点生成图像数据值。
另外,按照本发明的一个实施例,上述系统还包括:投影减法单元,用于将投影减法作用于该投影数据,以生成减法数据;希耳伯特滤波单元,用于对该减法数据应用希耳伯特型滤波器,以生成希耳伯特滤波数据;投影加法单元,用于将投影加法作用于希耳伯特滤波数据和斜坡滤波数据,以生成滤波的数据;加权单元,用于将冗余加权作用于滤波的数据,以生成加权数据,其中该反向投影单元用于将具有距离倒数加权的反向投影算子作用于由该加权单元产生的加权数据。
通过下列附图和优选实施例的具体描述,本发明的其它方法、系统和计算机程序产品对本领域普通技术人员来说是显而易见的。
附图说明
容易获得本发明更完整的说明及其许多附带的优点,结合附图参照以下具体描述,更容易理解这些内容,其中:
图1表示具有螺旋几何结构的源轨迹;
图2表示等角圆柱形检测器;
图3表示共线检测器;
图4A-4E分别表示全扫描、过扫描、短扫描、超短扫描和弹性重建范围;
图5表示π Tam窗口;
图6表示Tam窗口重建范围;
图7表示扩展的Tam窗口;
图8A-8F表示滤波线;
图9A表示按照本发明,在全扫描的重建点确定图像数据值的方法;
图9B是执行图9A所示方法的系统;
图10A表示按照本发明,在弹性投影范围、超短投影范围、短扫描或过扫描的重建点确定图像数据值的方法;
图10B是执行图10A所示方法的系统;
图11是描述可在本发明中使用的不同源轨迹、投影范围、加权函数以及滤波方向的表格;
图12表示CT设备;以及
图13A-13D表示四种重建算法。
具体实施方式
本发明针对扇形射束数据的精确算法和圆锥射束数据的准精确算法。将这些算法用于由CT装置产生的图像重建中。下面给出本发明的实施例。
扇形射束数据的2π投影范围(全扫描)方程:
扇形射束数据的弹性投影范围,超短扫描,短扫描或过扫描方程:
圆锥射束数据的2π投影范围(全扫描)方程:
圆锥射束数据的弹性投影范围,超短扫描,短扫描或过扫描方程:
本发明的上述实施例针对等角检测器给出。本发明的算法不仅限于等角检测器,还对其它检测器的几何结构起作用如共线检测器、等距检测器、非等距检测器、平面检测器、圆柱形检测器、球形检测器、倾斜的和PI掩蔽检测器。
另外,本发明的算法独立于源轨迹。源轨迹不需限于圆形或螺旋轨迹。例如可使用鞍形轨迹。
图9A表示全扫描的本发明的重建算法的一个实施例。在步骤900,通过公知方法获取CT投影数据。函数g(s,γ)或g(s,Θ)分别为图9A中步骤900表示的扇形射束数据或圆锥射束数据的投影数据。
在步骤901,该数据用斜坡型滤波器滤波,生成滤波的数据。该斜坡型滤波器可以是斜坡滤波器、改进的斜坡滤波器或者具有DC偏置改进的斜坡滤波器。
在步骤902,1/L加权的反向投影应用于该斜坡滤波的数据,以在CT图像中产生图像数据值。
在步骤903,确定在该CT图像中是否有需要重建的其它图像数据值(CT数)。如果有其它图像数据值,重复步骤902,直至在该CT图像中没有需要重建的其它图像数据值。
在步骤904,通过按照重建点设置该图像数据值,从步骤902输出的图像值用于生成CT图像。
图9B表示在CT图像中重建图像数据值的系统。CT扫描单元951产生投影数据。该CT扫描单元可向处理器960或向存储单元952传输该投影数据。处理器960用于直接从CT扫描单元951或从存储单元952存取接收该投影数据。处理器960包括滤波单元961,反向投影单元962和输出单元963。滤波单元961用于对该投影数据应用斜坡型滤波器,以生成斜坡滤波的数据。反向投影单元962用于使用反向投影算子对斜坡滤波的数据应用1/L加权,以在重建点上生成图像数据值。然后输出单元963向显示器970、存储单元971或两者输出图像数据值。斜坡滤波单元961和反向投影单元962均可存储并恢复来自存储单元952的信息。
图10A表示本发明的重建算法的另一个实施例。对弹性范围、超短扫描、短扫描或过扫描的扇形射束和圆锥射束方程使用g(s,γ)或g(s,Θ)分别表示图10A中步骤1000表示的扇形射束数据或圆锥射束数据的投影数据。
在步骤1001,投影减法应用于投影数据,生成减法数据。投影减法是在本发明的上述算法中对偏微分项的应用。
在步骤1002,希耳伯特型滤波器应用于该减法数据,以生成希耳伯特滤波数据。在步骤1002,可使用希耳伯特滤波器或改进的希耳伯特滤波器。
在步骤1003,斜坡型滤波器应用于投影数据,生成斜坡滤波数据。
在步骤1004,希耳伯特滤波数据和斜坡滤波数据相结合,生成滤波数据。
在步骤1005,冗余加权函数应用于该滤波数据。冗余加权函数w(s,γ)在所述重建算法中不是特定的,可在Parker加权、通用Parker加权(MHS)、Noo的加权、过扫描加权、准圆锥射束加权以及Tam窗口加权{wP(s,γ),wMHS(s,γ),wN(s,γ),wOS(s,γ),w3D(s,γ),wT(γ,v)}中自由选择。
图11的表格包括使用哪种加权函数,用哪种类型的射束数据。图11表示本发明的算法对于扇形射束数据和圆锥射束数据,如何作用于不同的源轨迹、不同的检测器几何结构,可使用什么投影范围,可使用什么加权函数以及可使用什么滤波方向。投影范围Λ决定加权函数的选择。对扇形射束数据,超短扫描使用Noo加权;短扫描使用Parker加权、MHS加权或Noo加权;过扫描使用OS加权或Noo加权;弹性扫描使用Noo加权。除此之外当该投影范围是短扫描或过扫描时,可使用Q3D加权,圆锥射束数据与扇形射束数据相同。Tam窗口加权也可与圆锥射束数据一起使用。
在步骤1006,该加权数据通过距离的倒数加权作用于反向投影算子,以生成图像数据值。
在步骤1007,确定在该CT图像中是否有需要重建的其它图像数据值。如果有其它图像数据值,重复步骤1005-1007,直至在该CT图像中没有其它图像数据值。
在步骤1008,输出该图像数据值,通过按照重建点设置该图像数据值,生成CT图像。
图10B表示在CT图像中处理数据图像值的系统。CT扫描单元951获取扇形射束数据或圆锥射束数据。该CT扫描单元可向处理器960或向存储单元952传输该射束数据。处理器960直接从CT扫描单元951或从存储单元952存取接收该射束数据。在该实施例中,处理器960包括斜坡滤波装置961,投影减法单元964,希耳伯特滤波单元965,投影附加单元966,加权单元967,反向投影单元962以及输出单元963。投影减法单元964对来自CT扫描单元951或存储单元952的射束数据应用投影减法,以生成减法数据。然后希耳伯特滤波单元965对该减法数据应用希耳伯特型滤波器,以生成希耳伯特滤波数据。斜坡滤波单元961对来自CT扫描单元951或存储单元952的射束数据应用斜坡型滤波器,以生成斜坡滤波数据。投影附加单元966将该斜坡滤波数据与该希耳伯特滤波数据结合,以生成滤波数据。加权单元967对该滤波数据应用加权函数,以生成加权数据。反向投影单元962对该加权数据应用1/L加权的反向投影算子,以按照重建点生成数据图像值。输出单元963向显示器970、存储单元971或两者输出该数据图像值。希耳伯特滤波单元965,投影附加单元966,加权单元967,反向投影单元962以及输出单元963均可存储并恢复来自存储单元952的信息。
可以在扫描模式下应用这里描述的实施例,其中X射线源在圆形、螺旋形或鞍形轨迹上移动。该算法独立于该检测器的几何结构,即可使用等角的、等距的、非等距的、平的、圆柱的、球的、倾斜的、旋转的、PI掩蔽的等等。该方程独立于所使用的滤波线的类型;水平的,切向的、旋转的、Katsevich或任意其它合理的滤波线族。可使用超短扫描、短扫描、过扫描以及满足该充分条件的任意轨迹。
本发明与其它CT图像重建算法的比较
重建算法可应用于不同的源轨迹,不同的检测器几何结构使用不同的滤波方向和不同的加权函数。但是,重建流和主要步骤对每个算法是唯一的。为了更清楚地看到该算法的主要特征,可使用算子表达式。下面对具有冗余加权螺旋圆锥射束几何结构的方程进行比较,该加权表示最实用的关注区。对等角检测器几何结构重写所有的方程。表示以下算法的流程图在图13A-D中表示,按照本发明的算法在图10A中表示。
[KRND]
[本发明的一个实施例]
空间均匀性取决于该反向投影加权。1/L反向投影加权的算法具有更佳的空间均匀性。因此,NDCK、Katsevich和本发明的算法产生比GFDK和KRND具有更佳的空间均匀性的图像。另一方面,在Katsevich的算法中,不同的图像象素具有不同的投影范围,导致更差的空间均匀性。图10和图13A-D的方框表示每个算法的反向投影加权。
由于本段的讨论不具有圆锥射束赝像,因此与[Katsevich]无关。所有其它算法仅对具有圆形轨迹的扇形射束几何结构是精确的(即在二维空间中)。通过扩展到三维空间,它们还可为各个相关小圆锥角和螺距提供相当良好的结果。但是,该算法在三维空间中不能同样良好地执行。圆锥射束赝像遮盖在物体的一个侧面上,并在另一个侧面上闪烁。[NDCK]和[KRND]伴有切向滤波,可减少但不会消除该圆锥射束赝像。由于切向滤波(以及旋转滤波)还可应用于[GFDK]和本发明的算法,因此不能视为仅[NDCK]和[KRND]的优点。由于[NDCK]和本发明的算法包含项,即相邻投影之间的微分,因此具有较少的圆锥射束赝像。这就补偿了螺旋数据的不一致性。注意圆锥射束赝像在Katsevich的算法中不出现,该算法也具有项。
圆锥射束赝像的主要原因之一是对螺旋圆锥射束数据应用扇形射束冗余加权。使用Tam窗口冗余加权函数使圆锥射束赝像显著减少。表1表示不同算法可使用哪些加权。
OS | MHS | Noo | Tam窗口 | |
GFDK | * | * | - | - |
NDCK | * | * | * | * |
KRND | * | * | * | - |
Katsevich | - | - | - | * |
本发明 | * | * | * | * |
表1.不同算法的冗余加权
注意仅[NDCK]和本发明使用所有加权。
当对CT重建算法相互比较时,不管是否有弹性重建范围和软件实施的简化,因子诸如体积重建的速度是很重要的。
体积重建速度主要由在切片重建环中执行多少操作确定。在图13A-D中,[GFDK]的切片重建环包含滤波,这意味着在每个图像切片中对相同的投影重卷积多次。所有其它算法更加有效:每个投影仅卷积一次(对于[NKDC]和Katsevich)或两次(对于[KRND]和图10A中本发明的一个实施例);不需要重卷积。但是,[KRND]切片重建环比[NKDC]、[Katsevich]和本发明更复杂。注意反向投影是反向投影中最需要计算的部分,每个切片周期的仅一个反向投影是商用CT重建设备所强烈期望的。
弹性重建范围意味着源轨迹的任意子集可用于精确重建,该轨迹在xy平面上的投影满足二维数据充分条件。在结构上[NKDC]、[KRND]和本发明具有弹性重建范围。[Katsevich]由于使用Tam窗口加权,因此没有弹性重建范围。弹性重建范围还意味着超短扫描的可能性。具有这种可能性的算法为:[NKDC]、[KRND]和本发明。
实施的简化对商用CT重建算法是关键的。一个重要的准则是,每个步骤(滤波、加权、反向投影)可表示为单独的模块。另外,[NKDC]中的角度微分取决于它如何实施,并会变得相当复杂。在本发明中,仅有的微分为简单的投影减法。
数字数据微分导致分辨率降低。在[GFDK]、[KRND]和本发明中,由于图像的主要部分从斜坡滤波数据重建,因此没有分辨率损失,且不经过微分步骤。
表2表示在此情况下所述算法的主要特征:
表2 性能比较
图12表示X射线计算机断层摄影装置,可用于获得由本发明的方法所处理的数据。由台架1组成的投影数据测量系统包括X射线源3,用于生成近似圆锥形的X射线流的圆锥射束,还包括二维阵列型X射线检测器5,包括以二维方式排列的多个检测器元件,即多个以一维多行堆叠方式排列的元件。X射线源3和二维阵列型X射线检测器5安装在旋转环2上,位于床6的滑动板上的主体(subject)的相对面上。二维阵列型X射线检测器5安装在旋转环2上。每个检测器元件与一个通道相对应。来自X射线源3的X射线通过X射线滤波器4指向主体。通过二维阵列型X射线检测器5,对已经通过该主体的X射线作为电信号进行检测。
X射线控制器8向高电压发生器7提供触发信号。高电压发生器7按照接收触发信号的时序向X射线源3典型地为X射线管提供高电压。这导致X射线从X射线源3发出。台架/床控制器9同步地控制台架1的旋转环2的旋转和床6的滑动板的滑动。系统控制器10组成整个系统的控制中心,控制X射线控制器8和台架/床控制器9,如从该主体所见,使X射线源3执行所谓螺旋扫描,其中它沿着螺旋路线移动。特别地,旋转环2以固定角速度连续旋转,同时该滑动板以固定速度移动,X射线以固定间隔角连续或间歇地从X射线源3发出。
二维阵列型X射线检测器5的输出信号由每个通道的数据收集单元11放大并转化为数字信号,以生成投影数据。从数据收集单元11输出的投影数据提供给重建处理单元12。重建处理单元12利用该投影数据找到反向投影数据,该反向投影数据在每个体素中反射该X射线的吸收。如第一实施例所述,在使用X射线的圆锥射束的螺旋扫描系统中,成像区域(有效视野)呈绕旋转轴中心半径为ω圆柱形。重建处理单元12在该成像区域中定义了多个体素(三维象素),并为每个体素找到反向投影数据。利用该反向投影数据编码的三维图像数据或断层摄影数据发送到显示装置14,在那里在视觉上显示为三维图像或断层摄影图像。该扫描装置包括单元1-11。
为了这种描述的目的,我们应该定义一种图像,以表示物理场景,其中该图像已经由某种成像技术生成。成像技术的实例可包括电视机或CCD照相机或X射线,声纳或超声波成像装置。记录图像的初始媒介可以是电子固态装置,光学记录胶片或某种其它装置,例如可光激发的荧光。然后记录的图像可通过电子结合(如在CCD信号的情况下)或机械/光手段(如在对光学记录胶片进行数字化或对来自可光激发的荧光的数据进行数字化的情况下)转化为数字形式。
利用一般的通用计算机或按照本发明的教导的可编程的微处理器可方便地实施本发明的所有实施例,这对计算机领域的普通技术人员是显而易见的。合适的软件可以由基于本发明教导的普通编程人员容易地制定,这对软件领域的普通技术人员是显而易见的。特别地,计算机外壳可容纳主板,包含CPU、存储器(如DRAM、ROM、EPROM、EEPROM、SRAM、SDRAM以及FlashRAM)及其它可选的特定目的的逻辑器件(如ASICS)或可配置逻辑器件(如GAL和可再编程FPGA)。该计算机还包括多个输入装置,(如键盘和鼠标),以及控制显示器的显示卡。另外,该计算机可包括软盘驱动器;其它可移动的媒介装置(如小型盘、磁带和可移动的磁光媒介);还包括硬盘或其它固定的高密度的媒介驱动器,并利用合适的装置总线连接(如SCSI总线、增强IDE总线或超DMA总线)。该计算机还可包括小型盘阅读器,小型盘读/写单元或小型盘光盘机,它们可连接于同一装置总线或连接于另一个装置总线。例如,典型地,在计算机系统上执行的计算机程序包括获得该主体投影数据的装置,生成斜坡滤波数据的装置,用于通过一维斜坡滤波获得的投影数据,还包括对所述斜坡滤波数据应用反向投影处理的装置,其中用对应于从焦点到重建点的距离的倒数的加权处理对该反向投影进行处理,以在CT图像中的重建点生成图像数据值。
与本发明相关的计算机可读媒介包括小型盘、硬盘、软盘、磁带、磁光盘、PROMs(如EPROM、EEPROM、Flash EPROM)、DRAM、SRAM、SDRAM等。存储在这些计算机可读媒介的任何一个或结合体上,本发明包括控制计算机硬件和使计算机与用户交互的软件。这种软件可包括但不限于装置驱动器、操作系统和用户应用程序,如开发工具。本发明的计算机程序产品包括存储计算机程序指令的任何计算机可读媒介(如计算机代码装置),当由计算机执行时,引起该计算机执行本发明的方法。本发明的计算机代码装置可以是任意可注释或可执行的代码机构,包括但不限于描述语言、注释、动态链接库、Java类以及完整的可执行程序。此外,为了更好的性能、稳定性和/或更低的成本,可对本发明的处理部分进行分配(如(1)多重CPUs或(2)至少一个CPU和至少一个可配置逻辑器件)。例如,可在第一计算机上选择轮廓或图像,并发送至第二计算机进行远程诊断。
本发明还可通过应用特定集成电路的准备工作或通过相互连接合适的常规元件电路网络得以实施,这对本领域普通技术人员来说是显而易见的。
按照本发明的图像数据源可以是任何合适的图像获取装置,如X射线机器或CT设备。另外,如果获取的数据未转换成数字形式,则可将其数字化。可选地,正在获取和处理的图像数据源可以是存储由图像获取装置产生的数据的存储器,该存储器可以是本地或远程的,其中可使用数据通信网络,如PACS(图片存档计算机系统),以存取按照本发明处理的图像数据。
根据上述教导,本发明的许多修改和变化都是可能的。因此可以理解,在所附的权利要求书的范围内,可以对本发明在特定描述之外实施。
Claims (20)
1、一种在扫描对象的计算机断层摄影图像中的重建点确定图像数据值的重建方法,其特征在于,包括:
获取扫描对象的投影数据;
用一维斜坡型滤波器滤波获取的投影数据,以生成斜坡滤波数据;以及
将具有对应从焦点至重建点距离的倒数的加权处理的反向投影处理作用于所述斜坡滤波数据,以在计算机断层摄影图像中的重建点生成图像数据值。
2、如权利要求1所述的重建方法,其特征在于,进一步包括:
在医学图像中对多个重建点重复所述的作用步骤,以获得多个图像数据值;以及
通过按照所述多个重建点设置所述多个图像数据值,生成计算机断层摄影图像。
3、如权利要求1所述的重建方法,其特征在于,进一步包括:
将投影减法作用于获取的投影数据,以生成减法数据;
将希耳伯特滤波器作用于该减法数据,以生成希耳伯特滤波数据;
将投影加法作用于希耳伯特滤波数据和斜坡滤波数据,以生成滤波的数据;以及
将冗余加权作用于所述滤波的数据,以生成加权数据,
其中,应用反向投影处理的步骤作用于所述加权数据,以在计算机断层摄影图像的重建点生成图像数据值。
4、如权利要求3所述的重建方法,其特征在于,进一步包括:
对计算机断层摄影图像中多个重建点重复应用冗余加权的步骤和应用反向投影处理的步骤;
通过按照所述多个重建点设置所述多个图像数据值,生成计算机断层摄影图像。
5、如权利要求1所述的重建方法,其特征在于,该滤波步骤包括:
用改进的斜坡滤波器滤波所述投影数据,在改进的斜坡滤波器中用核函数q(sinγ)代替q(γ)。
6、如权利要求1所述的重建方法,其特征在于,该滤波步骤包括:
用具有DC偏置的改进的斜坡滤波器滤波所述投影数据,在改进的斜坡滤波器中用核函数q(sinγ)代替q(γ)。
7、如权利要求3所述的重建方法,其特征在于,应用希耳伯特滤波器的步骤包括:
对所述减法数据应用改进的希耳伯特滤波器,在改进的希耳伯特滤波器中用核函数h(sinγ)代替h(γ)。
8、如权利要求1所述的重建方法,其特征在于,该获取步骤包括:
利用具有圆形源轨迹、螺旋源轨迹、和鞍形源轨迹之一的计算机断层摄影系统获取该投影数据。
9、如权利要求3所述的重建方法,其特征在于,应用冗余加权的步骤包括:
对所述滤波数据应用(1)Parker加权,(2)MHS加权,(3)OS加权,(4)Noo加权,(5)Q3D加权,(6)Tam窗口加权中的一种。
10、如权利要求1所述的重建方法,其特征在于,该斜波滤波的步骤包括:
利用(1)水平滤波,(2)切向滤波,(3)旋转滤波,(4)Katsevich滤波中的一种对获取的投影数据进行滤波。
11、如权利要求3所述的重建方法,其特征在于,应用希耳伯特滤波器的步骤包括:
对减法数据应用(1)水平滤波,(2)切向滤波,(3)旋转滤波,(4)Katsevich滤波中的一种。
12、如权利要求1所述的重建方法,其特征在于,该获取步骤包括:
利用满足Λ≥π+2arcsin(rFOV/R)的投影范围获取该投影数据,其中Λ是投影范围,rFOV是视野(FOV)的径向量度,R是源轨迹的径向量度。
13、一种X射线计算机断层摄影设备,用于在扫描对象的计算机断层摄影图像中的重建点确定图像数据值,其特征在于,包括:
计算机断层摄影扫描单元,用于生成扫描对象的投影数据,该计算机断层摄影扫描单元包括生成X射线的X射线源和具有检测器元件的检测器,该检测器用于产生扫描对象的投影数据;以及
处理器,包括:
滤波单元,用于对由计算机断层摄影扫描单元生成的投影数据应用斜坡滤波器,以生成斜坡滤波数据;以及
反向投影单元,用于对所述斜坡滤波数据应用具有对应从焦点至重建点的距离的倒数的加权处理的反向投影处理,以在计算机断层摄影图像中的重建点生成图像数据值。
14、如权利要求13的X射线计算机断层摄影设备,其特征在于,该处理器进一步包括:
投影减法单元,用于将投影减法作用于由计算机断层摄影扫描单元生成的投影数据,以生成减法数据;
希耳伯特滤波单元,用于对由投影减法单元产生的减法数据应用希耳伯特滤波器,以生成希耳伯特滤波数据;
投影加法单元,用于将投影加法作用于(1)由希耳伯特滤波单元产生的希耳伯特滤波数据,以及(2)由滤波单元产生的斜坡滤波数据,以生成滤波的数据;以及
加权单元,用于将冗余加权作用于由投影加法单元产生的滤波的数据,以生成加权数据,
其中,该反向投影单元用于将具有距离的倒数的加权的反向投影处理作用于由该加权单元产生的加权数据。
15、如权利要求13的X射线计算机断层摄影设备,其特征在于,该滤波单元用于对由计算机断层摄影扫描单元产生的投影数据应用改进的斜坡滤波器,在改进的斜坡滤波器中用核函数q(sinγ)代替q(γ)。
16、如权利要求13的X射线计算机断层摄影设备,其特征在于,该滤波单元用于对由计算机断层摄影扫描单元产生的投影数据应用具有DC偏置的改进的斜坡滤波器,在改进的斜坡滤波器中用核函数q(sinγ)代替q(γ)。
17、如权利要求14的X射线计算机断层摄影设备,其特征在于,该希耳伯特滤波单元用于对由投影减法单元产生的减法数据应用改进的希耳伯特滤波器,在改进的希耳伯特滤波器中用核函数h(sinγ)代替h(γ)。
18、如权利要求13的X射线计算机断层摄影设备,其特征在于,该计算机断层摄影扫描单元利用圆形源轨迹、螺旋源轨迹和鞍形轨迹之一生成投影数据。
19、如权利要求14的X射线计算机断层摄影设备,其特征在于,该加权单元对由投影加法单元产生的滤波的数据应用(1)Parker加权,(2)MHS加权,(3)OS加权,(4)Noo加权,(5)Q3D,以及(6)Tam窗口加权中的一种。
20、如权利要求13的X射线计算机断层摄影设备,其特征在于,该扫描单元利用满足Λ≥π+2arcsin(rFOV/R)的投影范围生成该投影数据,其中Λ是投影范围,rFOV是视野(FOV)的径向量度,R是源轨迹的径向量度。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/951,650 US7424088B2 (en) | 2004-09-29 | 2004-09-29 | Image reconstruction method using Hilbert transform |
US10/951,650 | 2004-09-29 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1788684A CN1788684A (zh) | 2006-06-21 |
CN100466981C true CN100466981C (zh) | 2009-03-11 |
Family
ID=35431883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005101373058A Expired - Fee Related CN100466981C (zh) | 2004-09-29 | 2005-09-29 | 图像重建方法和x射线计算机断层摄影设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7424088B2 (zh) |
EP (1) | EP1643446B1 (zh) |
JP (1) | JP5019193B2 (zh) |
CN (1) | CN100466981C (zh) |
DE (1) | DE602005001310T2 (zh) |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006073584A2 (en) * | 2004-11-24 | 2006-07-13 | Wisconsin Alumni Research Foundation | Cone-beam filtered backprojection image reconstruction method for short trajectories |
WO2006058124A1 (en) * | 2004-11-24 | 2006-06-01 | Wisconsin Alumni Research Foundation | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data |
US20060198491A1 (en) * | 2005-03-04 | 2006-09-07 | Kabushiki Kaisha Toshiba | Volumetric computed tomography system for imaging |
JP2008539930A (ja) * | 2005-05-12 | 2008-11-20 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 超短スキャン及び最新データのより強い重み付けを実行する連続コンピュータ断層撮影 |
US7653230B2 (en) * | 2006-02-21 | 2010-01-26 | General Electric Company | Methods and systems for image reconstruction using low noise kernel |
US7409033B2 (en) * | 2006-05-31 | 2008-08-05 | The Board Of Trustees Of The Leland Stanford Junior University | Tomographic reconstruction for x-ray cone-beam scan data |
WO2008064367A2 (en) * | 2006-11-24 | 2008-05-29 | Kabushiki Kaisha Toshiba | Method and system for tomographic reconstruction in medical imaging using the circle and line trajectory |
JP5426075B2 (ja) * | 2007-01-31 | 2014-02-26 | 株式会社東芝 | X線ct装置 |
JP2009022450A (ja) * | 2007-07-18 | 2009-02-05 | Ge Medical Systems Global Technology Co Llc | X線ct装置および画像作成方法 |
US8494111B2 (en) * | 2007-11-02 | 2013-07-23 | Kabushiki Kaisha Toshiba | System and method for image reconstruction for helical cone beam computed tomography with factorized redundancy weighting |
DE102007054371B4 (de) * | 2007-11-14 | 2014-09-04 | Siemens Aktiengesellschaft | Verfahren zur Bilddarstellung von Projektionsdaten |
US7792238B2 (en) * | 2008-02-18 | 2010-09-07 | General Electric Company | Method and system for reconstructing cone-beam projection data with reduced artifacts |
US8284892B2 (en) * | 2008-12-22 | 2012-10-09 | General Electric Company | System and method for image reconstruction |
WO2010109352A1 (en) * | 2009-03-26 | 2010-09-30 | Koninklijke Philips Electronics N.V. | Method and apparatus for computed tomography image reconstruction |
WO2010121043A2 (en) * | 2009-04-15 | 2010-10-21 | Virginia Tech Intellectual Properties, Inc. | Exact local computed tomography based on compressive sampling |
US9665953B2 (en) | 2009-06-30 | 2017-05-30 | Analogic Corporation | Efficient quasi-exact 3D image reconstruction algorithm for CT scanners |
JP5433334B2 (ja) * | 2009-07-27 | 2014-03-05 | 株式会社東芝 | X線ct装置 |
EP2315178A1 (en) * | 2009-10-12 | 2011-04-27 | Siemens Aktiengesellschaft | Reconstruction of 3D image datasets from X-ray cone-beam data |
US8824760B2 (en) * | 2009-10-20 | 2014-09-02 | University Of Utah Research Foundation | Modification and elimination of back projection weight during the CT image reconstruction |
US8731266B2 (en) * | 2009-12-17 | 2014-05-20 | General Electric Company | Method and system for correcting artifacts in image reconstruction |
US8204172B1 (en) * | 2010-03-17 | 2012-06-19 | General Electric Company | System and method of prior image constrained image reconstruction using short scan image data and objective function minimization |
DE102010020770A1 (de) * | 2010-05-17 | 2011-11-17 | Siemens Aktiengesellschaft | Verfahren zur Reduzierung von Bildartefakten, insbesondere von Metallartefakten, in CT-Bilddaten |
DE102011076348B4 (de) * | 2011-05-24 | 2014-02-13 | Siemens Aktiengesellschaft | Verfahren und Computertomographiesystem zur Erzeugung tomographischer Bilddatensätze |
DE102011076547A1 (de) * | 2011-05-26 | 2012-11-29 | Siemens Aktiengesellschaft | Verfahren zum Gewinnen eines 3D-Bilddatensatzes zu einem Bildobjekt |
US8805037B2 (en) | 2011-05-31 | 2014-08-12 | General Electric Company | Method and system for reconstruction of tomographic images |
US9076255B2 (en) * | 2011-05-31 | 2015-07-07 | General Electric Company | Method and system for reconstruction of tomographic images |
US8837797B2 (en) * | 2012-01-10 | 2014-09-16 | Kabushiki Kaisha Toshiba | Spatial resolution improvement in computer tomography (CT) using iterative reconstruction |
CN102590243B (zh) * | 2012-02-17 | 2013-09-04 | 重庆大学 | 一种铁路铸件全身ct扫描成像方法 |
US8798350B2 (en) * | 2012-03-22 | 2014-08-05 | Kabushiki Kaisha Toshiba | Method and system for reconstruction algorithm in cone beam CT with differentiation in one direction on detector |
US8948337B2 (en) | 2013-03-11 | 2015-02-03 | General Electric Company | Computed tomography image reconstruction |
KR20150022176A (ko) | 2013-08-22 | 2015-03-04 | 삼성전자주식회사 | 영상 처리 유닛, 영상 재구성 방법 및 방사선 영상 장치 |
CN103714560B (zh) * | 2013-12-27 | 2016-10-05 | 哈尔滨工业大学深圳研究生院 | 一种基于Katsevich算法的图像重建方法和系统 |
CN110251047B (zh) * | 2014-03-28 | 2022-01-18 | 直观外科手术操作公司 | 手术植入物的定量三维成像和打印 |
JP6854237B2 (ja) | 2014-03-28 | 2021-04-07 | インテュイティブ サージカル オペレーションズ, インコーポレイテッド | 視野内の器具の定量的三次元視覚化 |
US10555788B2 (en) | 2014-03-28 | 2020-02-11 | Intuitive Surgical Operations, Inc. | Surgical system with haptic feedback based upon quantitative three-dimensional imaging |
US9510787B2 (en) * | 2014-12-11 | 2016-12-06 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for reconstructing sampled signals |
CN106097411B (zh) * | 2016-06-29 | 2019-07-05 | 沈阳开普医疗影像技术有限公司 | Ct机图像重建方法及高分辨ct扫描机 |
EP3559911B1 (en) * | 2016-12-21 | 2020-05-20 | Koninklijke Philips N.V. | Redundancy weighting for short scan off-center detector x-ray tomography |
WO2019033390A1 (en) * | 2017-08-18 | 2019-02-21 | Shenzhen United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR IMAGE RECONSTRUCTION |
DE112020002334T5 (de) * | 2019-05-10 | 2022-02-17 | Koninklijke Philips N.V. | Kegelstrahlartefaktkorrektur für gesteuerte bildgebung |
CN111489408B (zh) * | 2020-04-21 | 2022-11-11 | 赛诺威盛科技(北京)股份有限公司 | 基于重叠3d半扫描数据的ct图像重建方法及装置、终端 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6018561A (en) * | 1998-07-27 | 2000-01-25 | Siemens Corporate Research, Inc. | Mask boundary correction in a cone beam imaging system using simplified filtered backprojection image reconstruction |
US6452996B1 (en) * | 2001-03-16 | 2002-09-17 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus utilizing generalized helical interpolation algorithm |
CN1418599A (zh) * | 2001-11-13 | 2003-05-21 | 皇家菲利浦电子有限公司 | 荧光检查计算机断层方法 |
CN1460452A (zh) * | 2002-05-22 | 2003-12-10 | Ge医疗系统环球技术有限公司 | 三维反投影方法和一种x射线计算机层析成像装置 |
CN1483384A (zh) * | 2002-08-13 | 2004-03-24 | GEҽ��ϵͳ����������˾ | 多行检测器x-射线计算断层成像设备 |
DE10248766A1 (de) * | 2002-10-18 | 2004-07-15 | Siemens Ag | Verfahren zur Bilderstellung für einen mehrzelligen Spiral-Computertomographen und Computer-Tomographie-Gerät zur Durchführung dieses Verfahrens |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5926521A (en) * | 1998-03-31 | 1999-07-20 | Siemens Corporate Research, Inc. | Exact region of interest cone beam imaging using 3D backprojection |
US6574299B1 (en) * | 2001-08-16 | 2003-06-03 | University Of Central Florida | Exact filtered back projection (FBP) algorithm for spiral computer tomography |
DE10304662A1 (de) * | 2003-02-05 | 2004-08-19 | Siemens Ag | Verfahren zur Erzeugung von Bildern in der Computertomographie mit einem 3D-Bildrekonstruktionsverfahren |
CN100359537C (zh) * | 2003-02-19 | 2008-01-02 | 皇家飞利浦电子股份有限公司 | 利用螺旋相对运动和锥形光束的计算机化x线断层摄影方法 |
US6990167B2 (en) * | 2003-08-29 | 2006-01-24 | Wisconsin Alumni Research Foundation | Image reconstruction method for divergent beam scanner |
-
2004
- 2004-09-29 US US10/951,650 patent/US7424088B2/en not_active Expired - Fee Related
-
2005
- 2005-09-14 JP JP2005267253A patent/JP5019193B2/ja active Active
- 2005-09-19 DE DE602005001310T patent/DE602005001310T2/de active Active
- 2005-09-19 EP EP05020346A patent/EP1643446B1/en not_active Expired - Fee Related
- 2005-09-29 CN CNB2005101373058A patent/CN100466981C/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6018561A (en) * | 1998-07-27 | 2000-01-25 | Siemens Corporate Research, Inc. | Mask boundary correction in a cone beam imaging system using simplified filtered backprojection image reconstruction |
US6452996B1 (en) * | 2001-03-16 | 2002-09-17 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus utilizing generalized helical interpolation algorithm |
CN1418599A (zh) * | 2001-11-13 | 2003-05-21 | 皇家菲利浦电子有限公司 | 荧光检查计算机断层方法 |
CN1460452A (zh) * | 2002-05-22 | 2003-12-10 | Ge医疗系统环球技术有限公司 | 三维反投影方法和一种x射线计算机层析成像装置 |
CN1483384A (zh) * | 2002-08-13 | 2004-03-24 | GEҽ��ϵͳ����������˾ | 多行检测器x-射线计算断层成像设备 |
DE10248766A1 (de) * | 2002-10-18 | 2004-07-15 | Siemens Ag | Verfahren zur Bilderstellung für einen mehrzelligen Spiral-Computertomographen und Computer-Tomographie-Gerät zur Durchführung dieses Verfahrens |
Non-Patent Citations (1)
Title |
---|
A new framework of image reconstruction from fan beamprojections. Guang-Hong Chen.Medical Physics,Vol.30 No.6. 2003 * |
Also Published As
Publication number | Publication date |
---|---|
JP5019193B2 (ja) | 2012-09-05 |
US20060067457A1 (en) | 2006-03-30 |
CN1788684A (zh) | 2006-06-21 |
EP1643446B1 (en) | 2007-06-06 |
JP2006095297A (ja) | 2006-04-13 |
EP1643446A1 (en) | 2006-04-05 |
DE602005001310D1 (de) | 2007-07-19 |
US7424088B2 (en) | 2008-09-09 |
DE602005001310T2 (de) | 2008-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100466981C (zh) | 图像重建方法和x射线计算机断层摄影设备 | |
US7251307B2 (en) | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data | |
Clackdoyle et al. | Quantitative reconstruction from truncated projections in classical tomography | |
Noo et al. | Exact helical reconstruction using native cone-beam geometries | |
De Man et al. | Distance-driven projection and backprojection | |
US7203272B2 (en) | Cone-beam filtered backprojection image reconstruction method for short trajectories | |
US6292525B1 (en) | Use of Hilbert transforms to simplify image reconstruction in a spiral scan cone beam CT imaging system | |
CN100469317C (zh) | 半径内图像依存检测器列滤波方法、x射线ct装置和程序 | |
CN1572253B (zh) | 锥面射束重建设备以及计算机断层摄影设备 | |
EP0502187A1 (en) | Parallel processing method and apparatus based on the algebra reconstruction technique for reconstructing a three-dimensional computerized tomography | |
EP1644897B1 (en) | A fourier tomographic image reconstruction method for fan-beam data | |
US7305061B2 (en) | Efficient image reconstruction algorithm for the circle and arc cone beam computer tomography | |
Bruder et al. | Single-slice rebinning reconstruction in spiral cone-beam computed tomography | |
US7197105B2 (en) | Efficient image reconstruction algorithm for the circle and line cone beam computed tomography | |
Weng et al. | A reconstruction algorithm for helical cone-beam SPECT | |
JP3006722B2 (ja) | コンピュータ断層撮影装置 | |
Grangeat et al. | Indirect cone-beam three-dimensional image reconstruction | |
Ramamurthi et al. | Exact 3d cone-beam reconstruction from projections obtained over a wobble trajectory on a c-arm | |
Zeng et al. | Review of recent developments in cone-beam ct reconstruction algorithms for long-object problem | |
JP3204651B2 (ja) | コンピュータ断層撮影装置 | |
JP2772289B2 (ja) | コンピュータ断層撮影装置 | |
Hu et al. | A simple derivation and analysis of a helical cone beam tomographic algorithm for long object imaging via a novel definition of region of interest | |
Mazin et al. | A Fourier rebinning algorithm for cone beam CT | |
Yu et al. | A rebinning-type backprojection-filtration algorithm for image reconstruction in helical cone-beam CT | |
Rodet et al. | Parallel algorithm based on a frequential decomposition for dynamic 3D computed tomography |
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 | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160720 Address after: Japan Tochigi Patentee after: TOSHIBA MEDICAL SYSTEMS Corp. Address before: Tokyo, Japan Patentee before: Toshiba Corp. Patentee before: TOSHIBA MEDICAL SYSTEMS Corp. |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090311 |
|
CF01 | Termination of patent right due to non-payment of annual fee |