CN103299345A - 双能量断层摄影成像系统 - Google Patents
双能量断层摄影成像系统 Download PDFInfo
- Publication number
- CN103299345A CN103299345A CN2011800646244A CN201180064624A CN103299345A CN 103299345 A CN103299345 A CN 103299345A CN 2011800646244 A CN2011800646244 A CN 2011800646244A CN 201180064624 A CN201180064624 A CN 201180064624A CN 103299345 A CN103299345 A CN 103299345A
- Authority
- CN
- China
- Prior art keywords
- projection
- covariance
- inter
- data
- imaging system
- 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
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
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种用于对对象成像的成像系统(30)。投影数据提供单元(31)提供包括至少两个成分的对象的采集的谱投影数据,并且重建单元(10)通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。由于基于所采集的谱投影数据和指示相关噪声的惩罚项更新所述至少两个中间成分图像,因此在迭代重建过程中惩罚相关噪声。最终获得的所述对象的成分图像因此较少地受到相关噪声的破坏,具有改进的图像质量。
Description
技术领域
本发明涉及一种用于对对象成像的成像系统、成像方法和成像计算机程序。
背景技术
Alvarez和Macowski的文章“Energy-selective Reconstructions in X-rayComputerized Tomography”(《Phys.Med.Biol》,第21卷,第5期,第733-744页,1976年)中公开了一种计算机断层摄影系统,该系统包括发射多色辐射的辐射源以及用于生成指示已经穿透对象后的辐射的投影数据的能量分辨探测器。在辐射源和对象相对于彼此移动时采集投影数据,以便采集不同投影方向的投影数据。向所采集的投影数据应用成分分解技术,以便确定分别由光电效应和康普顿效应造成的光电效应衰减值和康普顿效应衰减值。光电效应衰减值被反向投影,以生成对象的光电效应图像,并且康普顿效应衰减值被反向投影,以生成对象的康普顿图像。
光电效应图像和康普顿效应图像被噪声破坏,因此具有降低的图像质量。
发明内容
本发明的目的是提供一种用于对对象成像的成像系统、成像方法和成像计算机程序,其中,所述对象的重建图像的质量能够得以改进。
在本发明的第一方面中,提出了一种用于对对象成像的成像系统,其中,所述成像系统包括:
-投影数据提供单元,其用于提供包括至少两个成分的对象的采集的谱投影数据,
-重建单元,其通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。
由于基于所采集的谱投影数据和指示所述至少两个中间成分图像之间的相关噪声的惩罚项更新所述至少两个中间成分图像,因此所述至少两个中间成分图像之间的相关噪声在迭代重建过程中得以惩罚。最终获得的所述对象的成分图像因此受到相关噪声的更少破坏,具有了改进的图像质量。
所述对象能够包括两个或更多成分。
不同的成分能够是不同的基础材料,例如骨骼、软组织、金属、如碘的造影剂,和/或它们能够涉及不同的物理效应,例如光电效应、康普顿效应、K边缘效应等。
所述投影数据提供单元能够是存储单元,在其中已经存储了采集的投影数据,或者它能够是诸如X射线计算机断层摄影系统、X射线C形臂系统、核成像系统等的采集单元的投影数据采集单元所述核成像系统为例如正电子发射断层摄影系统或单光子发射计算机断层摄影系统。
优选地,所提供的采集的谱投影数据对应于通过如下方式已经采集的投影数据:将发射多色辐射的辐射源和所述对象相对于彼此进行移动,并且在所述辐射源和所述对象相对于彼此移动的同时,以谱方式探测所述投影数据,所述投影数据指示已经穿透所述对象之后的辐射。
所述至少两个中间成分图像之间的相关噪声能够是正相关噪声或负相关噪声。
优选地,所述重建单元适于通过优化a)似然项与b)所述惩罚项之和,来更新所述至少两个中间成分图像,所述似然项依赖于定义了贯穿所述中间成分图像的前向投影产生类似于所采集的谱投影数据的模拟谱投影数据的似然的似然函数。所述惩罚项是这样的项:其能够通过在迭代重建过程中修改所述中间成分图像被优化,从而减小了所述至少两个中间成分图像之间的相关噪声。进一步优选地,所述惩罚项依赖于对不同中间成分图像中对应图像元素的协方差的估计。对所述协方差的估计允许以相对简单的方式,通过所述惩罚项考虑相关噪声。
在优选实施例中,所述重建单元适于a)生成在不考虑所述惩罚项的情况下从所采集的投影数据重建的协方差估计成分图像,并且从所述协方差估计成分图像估计所述协方差,或b)从所采集的投影数据估计所述协方差。具体而言,所述重建单元能够适于通过在不考虑所述惩罚项的情况下,向所采集的投影数据应用成分分解技术以生成所述成分的线积分,并且通过从所述成分的所述线积分重建所述协方差估计成分图像,来生成所述协方差估计成分图像。优选地,所述重建单元适于将所述协方差估计成分图像的图像值求平均,并且适于基于经平均的图像值估计不同中间成分图像中的对应图像元素的协方差。这允许以相对低的计算工作量估计不同中间成分图像中对应图像元素的协方差。
同样优选地,所述重建单元适于通过如下方式估计不同中间成分图像中的对应图像元素的协方差:
-对所采集的投影数据应用基于似然的成分分解,其中,通过调整所述对象的能调整模型来优化依赖于定义了贯穿所述模型的前向投影产生类似于所采集的谱投影数据的模拟谱投影数据的似然的似然函数的似然项,由此生成所述成分的线积分,其中,所述模型表示至少两个成分的分布,
-依赖于经优化的似然项的曲率,确定不同成分的对应线积分的协方差,
-通过反向投影所确定的所述不同成分的对应线积分的协方差,估计不同中间成分图像中的对应图像元素的协方差。优选地,所述似然项是所述似然函数的对数,具体而言,是所述似然函数的负对数,并且优选地依赖于所述似然函数的对数在最优位置处的曲率确定不同成分的对应线积分的协方差。换言之,优选地,将所述似然函数的对数的曲率求得为对应于所述似然函数的优化对数的似然估计。这允许以相对高的准确度估计所述协方差。因此能够相对准确地考虑所述中间成分图像之间的相关噪声,由此进一步提高最终重建的最终成分图像的质量。
进一步优选地,所述惩罚项依赖于相同中间成分图像的图像元素的差。具体而言,所述惩罚项依赖于相同中间成分图像的相邻图像元素的差。因此,能够考虑相同中间成分图像的图像元素的局部差,其中,与较小的局部差相比较大的局部差优选地获得更多的惩罚,即,例如,所述惩罚项随着绝对局部差的增大而增大。由于假设相同中间成分图像的相邻图像元素更可能相对相似,而不是非常不同,因此所述惩罚项对局部差的依赖性能够进一步改进成分图像的质量。优选地,所述惩罚项依赖于相同中间成分图像的相邻图像元素的差的平方和/或不同成分图像的相邻图像元素的差的乘积,即,针对不同成分图像,计算相同成分图像的相邻图像元素的差,并且将不同成分图像的所得差相乘,以生成所述惩罚项可以依赖的乘积。
优选地,所述惩罚项还指示在所述中间成分图像中的每一个自身的噪声。由于不仅考虑了相关噪声,而且考虑了在所述中间成分图像中的每一个自身的噪声,因此能够进一步提高最终重建的最终成分图像的质量。
在本发明的另一方面中,提出了一种用于对对象成像的成像方法,其中,所述成像方法包括:
-提供包括至少两个成分的对象的采集的谱投影数据,
-通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。
在本发明的另一方面中,提出了一种用于对对象成像的成像计算机程序,其中,所述成像计算机程序包括程序代码模块,当在控制根据权利要求1所述的成像装置的计算机上运行所述成像计算机程序时,所述程序代码模块用于令所述成像装置执行根据权利要求11所述的成像方法的各步骤。
应当理解,根据权利要求1所述的成像系统、根据权利要求11所述的成像方法以及根据权利要求12所述的成像计算机程序具有类似和/或等同的优选实施例,如从属权利要求中定义的。
应当理解,本发明的优选实施例也能够是从属权利要求与相应独立权利要求的任意组合。
本发明的这些和其他方面将从下文描述的实施例变得显而易见,并参考下文描述的实施例加以阐述。
附图说明
在以下附图中:
图1示意性且示范性地示出了用于对对象成像的成像系统,以及
图2示出了示范性说明用于对对象成像的成像方法的实施例的流程图。
具体实施方式
图1示意性且示范性地示出了作为计算机断层摄影系统30的用于对对象成像的成像系统。计算机断层摄影系统30包括能够绕平行于z方向延伸的旋转轴R旋转的机架1。在机架1上安装有辐射源2,辐射源2在本实施例中为X射线管。辐射源2配备有准直器3,准直器3在本实施例中形成了来自由辐射源2生成的辐射的锥形辐射束4。所述辐射穿透包括对象的检查区域5,所述对象例如为人、动物或技术对象,其中,检查区域5在本实施例中为圆柱形的。所述对象包括至少两个成分,例如,软组织和骨骼。在已经穿透检查区域5之后,辐射束4入射到探测设备6上,探测设备6包括二维探测表面。探测设备6安装在机架1上。
计算机断层摄影系统30包括两个电动机7、8。由电动机7以优选恒定但能调节的角速度驱动机架1。电动机8被提供用于使布置于检查区域5中患者台上的诸如人的所述对象平行于旋转轴R或z轴的方向发生位移。这些电动机7、8受控制单元9的控制,例如,使得辐射源2和所述对象相对彼此沿着螺旋轨迹移动。然而,也能够沿另一轨迹执行相对移动。例如,也能够不移动所述对象,而只旋转辐射源2,即所述辐射源相对于所述对象沿着圆形轨迹移动。此外,在另一实施例中,准直器3能够适于形成另一射束形状,特别是扇形束,并且探测设备6能够包括这样的探测表面:所述探测表面的形状对应于另一射束形状,特别是扇形束。
辐射源2适于发射多色辐射,并且探测设备6为用于采集谱投影数据的能量分辨探测设备。
能量分辨探测设备的工作,例如,基于对入射光子计数并且输出示出了在特定能量箱中单位能量光子数目的信号的原理。例如,在Llopart X.等人的文章“First test measurements of a 64k pixel readout chip wordkin in asingle photon counting mode”(《Nuclear Instruments and Methods in PhysicsResearch》,A部分,第509卷,第1-3期,第157-163页,2003年)中,以及Llopart X.等人的文章“Medipix2:A 64-k pixel readout chip with 55 mumsquare elements working in a single photon counting mode”(《IEEETransactions on Nuclear Science》,第49卷,第5期,第2279-2283页,2000年)中,描述了这种能量分辨探测设备,这两篇文章以引用方式并入本文。
在辐射源2和所述对象的相对移动过程中,探测设备6依赖于在探测设备6的探测表面上入射的辐射而生成投影数据。因此,辐射源2,用于相对于所述对象移动辐射源2的元件,具体而言,电动机7、8,和机架1,以及探测设备6形成用于提供感兴趣区域的投影数据的探测值提供单元31。
所采集的投影数据被提供给用于迭代地重建所述对象的至少两个最终成分图像的重建单元10,之后所述至少两个最终成分图像被示出在显示器11上。
重建单元10适于通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项,更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。具体而言,重建单元10适于通过将a)似然项-log(L)与b)惩罚项P的和Ψ最小化,来更新所述至少两个中间成分图像,所述似然项依赖于定义了贯穿所述中间成分图像的前向投影产生类似于所采集的谱投影数据的模拟谱投影数据的似然的似然函数L,即,重建单元10适于使以下等式最小化:
ψ=-log(L)+P。 (1)
负对数似然-log(L)例如能够由以下等式定义,其中,噪声模型由泊松分布给出:
其中,表示在由指数j指示的成分图像中由指数i指示的图像元素,NS表示单位几何射线的独立谱测量的数量,例如,就光子计数而言为能量箱的数量,指数m指示相应的独立谱测量,例如,相应的能量箱,NR表示在投影数据的相应采集过程中采集的几何射线的总数量,指数r指示相应的几何射线,λmr表示谱测量的探测值,例如,由m指示的能量箱和由r指示的几何射线的探测值,并且表示探测值λmr的期望值。在这个实施例中,可以将对成分图像的依赖性明确写作:
其中,在等式(3)的第二行中,使用对应于由i指示的图像元素的几何贡献的图像元素权重将成分图像的前向投影包括至由r指示的几何射线的线积分。此外,表示即将穿透待成像对象之前的辐射谱,δm(E)表示第m个能量箱的谱灵敏度,表示沿由r指示的辐射线贯穿第j′个中间成分图像的吸收值的线积分,并且fj′(E)表示由j′指示的成分的依赖能量的吸收。
其中,α是正常数,NB表示成分数量,例如NB可以是2,如果所述对象包括骨骼和软组织,在i上的和是中间成分图像的所有Nv图像元素的和,在k上的和是在具有指数i的相应实际图像元素的邻域中的图像元素的和,其中,在邻域中的图像元素由Ni表示,并且wik是补偿在具有指数i的相应图像元素与在邻域中具有指数k的图像元素之间的不同距离的几何权重。随着距离的增加,几何权重wik优选地减少,例如,几何权重wik能够是逆距离。能够选择正常数α,以获得相对于似然L的负算法的惩罚项P的期望考虑。能够通过(例如)校准重建来获得优化的正常数α,其中,针对不同的常数α重建成分图像,并且选择特定α,这产生包括优化的噪声降低的重建成分图像。
如果在等式(4)中j与j′等同,则所述协方差降低了相应成分图像的相应图像元素中的噪声的方差。所述方差指示中间成分图像中的每一个自身的噪声,即,惩罚项P不但指示在不同中间成分图像之间的相关噪声,而且指示中间成分图像自身中的每一个自身的噪声。
为了估计所述协方差,重建单元10适于生成协方差估计成分图像,所述协方差估计成分图像是在不考虑所述惩罚项的情况下从所采集的投影数据重建的,并且重建单元10适于估计来自所述协方差估计成分图像的协方差。通过在不考虑所述惩罚项的情况下向所采集的投影数据应用成分分解技术以生成成分的线积分,并且通过从所述成分的线积分重建所述协方差估计成分图像,来生成所述协方差估计成分图像。使用依赖于相应成分的依赖能量的衰减的、将探测值描述为贡献组合的探测值模型,通过求解针对作为所采集的谱投影数据的多个依赖能量的检测值的等式组,而向所采集的投影数据应用所述成分分解技术。具体而言,重建单元10能够适于求解以下等式组:
其中,表示沿着由r指示的几何射线或辐射线贯穿由j指示的协方差估计成分图像的吸收值的线积分。
如果能量箱的数量至少等于成分的数量,则能够利用已知的数值方法求解所述等式组,其中,参量δm(E)和fj(E)是已知的,并且求解所述等式组的结果是线积分辐射谱和谱灵敏度δm(E)是成像系统的特性,并且从例如对应的测量获知。成分的谱吸收fj(E),例如,骨骼和软组织的谱吸收,也从测量和/或从文献获知。之后,能够通过例如在所确定的线积分上执行经过滤的反向投影或另一重建算法,来重建所述协方差估计成分图像。
为了估计来自所述协方差估计成分图像的协方差,重建单元10优选地适于将所述协方差估计成分图像的图像值求平均,并且适于基于经平均的图像值估计不同中间成分图像中的对应图像元素的协方差。例如,能够向所述协方差估计成分图像中的每一个应用平均滤波器,这为相应协方差估计成分图像的图像元素的每一个确定了相应图像元素和邻近图像元素的平均值。之后,能够根据以下等式估计两个成分图像之间的协方差:
在另一实施例中,所述重建单元还能够适于从所采集的投影数据估计所述协方差。具体而言,所述重建单元能够适于通过将基于似然的成分分解应用在所采集的投影数据上,来估计不同中间成分图像中的对应图像元素的协方差,其中,通过调整所述对象的能调整模型来优化依赖于定义了贯穿所述模型的前向投影产生类似于所采集的光谱投影数据的模拟光谱投影数据的似然的似然函数的似然项,由此生成由j指示的成分的线积分其中,所述模型表示所述对象的不同成分的分布。在这个实施例中,所述似然项是所述似然函数的负对数,其能够由以下等式定义:
之后,所述重建单元能够依赖于优化的似然项的曲率确定不同成分的对应线积分的协方差,具体而言,根据以下等式:
所述重建单元还能够适于通过反向投影所确定的不同成分的对应线积分的协方差,来估计不同中间成分图像中的对应图像元素的协方差。
在已经估计了所述协方差后,重建单元10通过迭代地重建成分图像,而将等式(1)最小化。
在这个实施例中,重建单元10适于基于坐标下降法,具体而言,基于Newton-Raphson方法,迭代地更新成分图像,其中,局部梯度和局部曲率根据以下等式确定更新步骤:
在等式(9)中,“+”指示在图像数据上加强非负约束的可能性。对于i=k,符号δik等于1,否则,其等于0。要强调的是,以上表达式中出现的协方差和方差优选地不被更新,并且在迭代过程中能被考虑为常数,尽管他们依赖于在图像中的位置。
在下文中,将参照图2所示的流程图示范性地描述用于对对象成像的成像方法的实施例。
在步骤101中,提供包括至少两个成分的对象的采集的谱投影数据。所采集的谱投影数据能够存储在存储单元中,并且能够从所述存储单元检索以提供相同的所采集的谱投影数据。然而,在步骤101中也能够执行实际采集,其中,例如,发射多色辐射的辐射源和所述对象相对于彼此移动,并且在所述辐射源和所述对象相对于彼此移动时,以谱方式探测指示已经穿透所述对象后的辐射的投影数据。在步骤102中,如上文所述,估计不同中间成分图像中对应图像元素的协方差,其中,在估计所述协方差同时,还生成成分的吸收值的线积分。在步骤103中,通过执行若干迭代步骤,重建所述对象的至少两个最终成分图像,其中,在迭代步骤中,基于所采集的谱投影数据和指示相关噪声的惩罚项,更新至少两个中间成分图像。在这个实施例中,所述惩罚项包括在所述至少两个中间成分图像之间的协方差,以便指示相关噪声。具体而言,通过将迭代重建算法应用于从所采集的投影数据确定的不同成分的线积分,来间接地使用所采集的谱投影数据。执行迭代重建以使得上述等式(1)最小化。在步骤104中,在显示器11上示出所得的至少两个最终成分图像。
用于计算机断层摄影成像的非谱迭代重建技术,例如,基于测量数据的似然的技术,能够在保留边缘结构的清晰度的同时,利用在图像中加强特定量平滑度的某种调整。所述调整能够是在图像域中惩罚最近的近邻之间的大图像值差的某种局部算子。在谱计算机断层摄影成像中,能够通过一种方式或另一方式将衰减分解到由像特定基础材料的特定成分造成的衰减中,所述特定基础材料例如为骨骼、软组织和/或造影剂材料。如果与迭代重建组合,在这种意义中增加了调整的可能性:图像域中的惩罚不但能够涉及在各种成分图像中的噪声方差,而且能够涉及在不同成分图像之间的相当干扰性的反相关性。通过将惩罚项加到依赖于似然函数的似然项,上文描述的成像系统和成像方法能够用于降低成分图像中的反相关噪声,其中,所述惩罚项指示不同成分图像之间正相关和/或负相关,即反相关的噪声,并且指示在不同成分图像中的每一个自身的噪声。
通常,像基础材料图像的谱分解的成分图像能够有相当大的噪声,因为不同成分,具体而言,特别是在生物组织中的不同基础材料具有非常相似的谱衰减性质。这种事实因此能够严重影响临床实践中谱分解的有用性。上文描述的成像系统和方法能够使用例如在投影域中获得的成分分解的经过滤的反向投影重建降低噪声,具体而言,相比较于分析重建的图像。
除了惩罚成分图像中的噪声本身之外,所述成像系统和成像方法在迭代过程中优选地惩罚成分图像中正相关噪声和/或负相关的噪声,即反相关噪声。这能够通过使用在迭代过程中针对不同成分图像的噪声方差和噪声协方差的空间分辨的估计来完成。这些估计能够通过,在第一步骤中,生成针对彼此独立的所有测量的似然估计,以及对所采集的投影数据上的方差和协方差的估计,而获得。在第二步骤中,能够使用经过滤的反向投影的分析形式,以获得针对成分图像的噪声方差和噪声协方差的空间分辨的估计。
尽管在上文描述的实施例中,所述成像系统是计算机断层摄影系统,但所述成像系统也能够是使用投影数据以迭代重建对象的图像的另一成像系统。例如,所述成像系统能够是X射线C形臂成像系统,或者像单光子发射计算机断层摄影系统或正电子发射断层摄影系统的核成像系统。
尽管在上文描述的实施例中,所述不同成分例如是像骨骼和软组织的不同基础材料,但所述不同成分也能够涉及像光电效应、康普顿效应、K边缘效应等不同的物理效应。
尽管在上文描述的实施例中,所述探测设备是能量分辨探测设备,但在其他实施例中,也能够使用非能量分辨探测设备,其中,相同的非能量分辨探测设备由具有不同能谱的辐射交替辐照,或其中,不同的非能量分辨探测设备由具有不同能谱的不同辐射辐照。例如,具体而言,如果所述成像系统是计算机断层摄影系统,那么能够使用双能量技术,或像四能量技术的双能量技术的组合。
尽管在上文描述的实施例中,描述了用于优化所述似然项,具体而言,用于优化等式(2)和等式(9)的特定措施,但也能够使用用于优化所述似然项的其他措施,像传输期望最大化算法(EM)(Transmission ExpectationMaximization algorithm),可分离的抛物面型替代函数算法(SPS)(SeparableParaboidal Surrogates algorithm)或加权最小二乘方法(LS)。
通过研究附图、说明书和权利要求书,本领域的技术人员在实施所要求保护的发明时能够理解和实现所公开实施例的其他变型。
在权利要求书中,词语“包括”不排除其他元件或步骤,并且单数词“一”或“一个”不排除多个。
单个单元或设备可以完成权利要求中记载的若干项目的功能。在互不相同的从属权利要求中记载的特定措施并不指示不能有利地使用这些措施的组合。
能够由任何其他数量的单元或设备执行由一个或若干单元或设备执行的,像协方差的估计或重建的计算。例如,能够由单个单元或由任何其他数量的不同单元执行步骤102和103。能够将根据所述成像方法的计算和/或对所述成像系统的控制实施为计算机程序的程序代码模块和/或实施为专用硬件。
可以在适当的介质上存储和/或分配计算机程序,所述适当的介质诸如与其他硬件一起供应或作为其他硬件一部分供应的光学存储介质或固态介质,但所述计算机程序也可以以其他形式分配,诸如经由互联网或其他有线或无线的远程通信系统。
权利要求中的任何附图标记不应被解释为限制范围。
本发明涉及一种用于对对象成像的成像系统。投影数据提供单元提供包括至少两个成分的对象的采集的谱投影数据,并且重建单元通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。由于基于所采集的谱投影数据和指示相关噪声的惩罚项更新所述至少两个中间成分图像,因此在迭代重建过程中惩罚相关噪声。因此,最后获得的所述对象的成分图像较少地受到相关噪声的破坏,具有改进的图像质量。
Claims (12)
1.一种用于对对象成像的成像系统,所述成像系统包括:
-投影数据提供单元(31),其用于提供包括至少两个成分的对象的采集的谱投影数据,
-重建单元(10),其通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。
2.根据权利要求1所述的成像系统,其中,所述重建单元(10)适于通过优化a)似然项与b)所述惩罚项之和,来更新所述至少两个中间成分图像,所述似然项依赖于定义了贯穿所述中间成分图像的前向投影产生类似于所采集的谱投影数据的模拟谱投影数据的似然的似然函数。
3.根据权利要求1所述的成像系统,其中,所述惩罚项依赖于对不同中间成分图像中对应图像元素的协方差的估计。
4.根据权利要求3所述的成像系统,其中,所述重建单元(10)适于a)生成在不考虑所述惩罚项的情况下从所采集的投影数据重建的协方差估计成分图像,并且从所述协方差估计成分图像估计所述协方差,或b)从所采集的投影数据估计所述协方差。
5.根据权利要求4所述的成像系统,其中,所述重建单元(10)适于通过在不考虑所述惩罚项的情况下,向所采集的投影数据应用成分分解技术以生成所述成分的线积分,并且通过从所述成分的所述线积分重建所述协方差估计成分图像,来生成所述协方差估计成分图像。
6.根据权利要求4所述的成像系统,其中,所述重建单元(10)适于将所述协方差估计成分图像的图像值求平均,并且适于基于经平均的图像值估计不同中间成分图像中的对应图像元素的协方差。
7.根据权利要求4所述的成像系统,其中,所述重建单元(10)适于通过如下方式估计不同中间成分图像中的对应图像元素的协方差:
-对所采集的投影数据应用基于似然的成分分解,其中,通过调整所述对象的能调整模型来优化依赖于定义了贯穿所述模型的前向投影产生类似于所采集的谱投影数据的模拟谱投影数据的似然的似然函数的似然项,由此生成所述成分的线积分,其中,所述模型表示至少两个成分的分布,
-依赖于经优化的似然项的曲率,确定不同成分的对应线积分的协方差,
-通过反向投影所确定的所述不同成分的对应线积分的协方差,估计不同中间成分图像中的对应图像元素的协方差。
8.根据权利要求1所述的成像系统,其中,所述惩罚项依赖于相同中间成分图像的图像元素的差。
9.根据权利要求1所述的成像系统,其中,所述惩罚项还指示所述中间成分图像中的每一个自身的噪声。
10.根据权利要求1所述的成像系统,其中,所提供的采集的谱投影数据对应于通过如下方式已经采集的投影数据:将发射多色辐射的辐射源和所述对象相对于彼此进行移动,并且在所述辐射源和所述对象相对于彼此移动的同时,以谱方式探测所述投影数据,所述投影数据指示已经穿透所述对象之后的辐射。
11.一种用于对对象成像的成像方法,所述成像方法包括:
-提供包括至少两个成分的对象的采集的谱投影数据,
-通过执行若干迭代步骤迭代地重建所述对象的至少两个最终成分图像,在所述若干迭代步骤中基于所采集的谱投影数据和惩罚项更新至少两个中间成分图像,所述惩罚项指示所述至少两个中间成分图像之间的相关噪声。
12.一种用于对对象成像的成像计算机程序,所述成像计算机程序包括程序代码模块,当在控制根据权利要求1所述的成像装置的计算机上运行所述成像计算机程序时所述程序代码模块用于令所述成像装置执行根据权利要求11所述的成像方法的各步骤。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11150468.4 | 2011-01-10 | ||
EP11150468 | 2011-01-10 | ||
PCT/IB2011/056004 WO2012095713A1 (en) | 2011-01-10 | 2011-12-29 | Dual -energy tomographic imaging system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103299345A true CN103299345A (zh) | 2013-09-11 |
CN103299345B CN103299345B (zh) | 2016-10-19 |
Family
ID=45531492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180064624.4A Active CN103299345B (zh) | 2011-01-10 | 2011-12-29 | 双能量断层摄影成像系统 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9165384B2 (zh) |
EP (1) | EP2663964A1 (zh) |
CN (1) | CN103299345B (zh) |
WO (1) | WO2012095713A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107077748A (zh) * | 2014-10-20 | 2017-08-18 | 皇家飞利浦有限公司 | 针对谱图像迭代重建的开始图像 |
CN107209946A (zh) * | 2015-01-20 | 2017-09-26 | 皇家飞利浦有限公司 | 图像数据分割和显示 |
CN107810518A (zh) * | 2015-06-26 | 2018-03-16 | 皇家飞利浦有限公司 | 利用相关噪声在图像上的边缘探测 |
CN113039582A (zh) * | 2018-11-16 | 2021-06-25 | 皇家飞利浦有限公司 | 用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9305379B2 (en) * | 2012-01-10 | 2016-04-05 | The Johns Hopkins University | Methods and systems for tomographic reconstruction |
US9978158B2 (en) | 2013-08-30 | 2018-05-22 | Koninklijke Philips N.V. | Spectral projection data de-noising with anti-correlation filter |
CN106462985B (zh) * | 2014-09-15 | 2018-09-21 | 皇家飞利浦有限公司 | 利用锐度驱动的正则化参数的迭代图像重建 |
US9687207B2 (en) * | 2015-04-01 | 2017-06-27 | Toshiba Medical Systems Corporation | Pre-reconstruction calibration, data correction, and material decomposition method and apparatus for photon-counting spectrally-resolving X-ray detectors and X-ray imaging |
CN109074635B (zh) * | 2016-05-03 | 2023-05-30 | 皇家飞利浦有限公司 | 用于对向量值图像去噪的设备和方法 |
US10559100B2 (en) | 2017-05-22 | 2020-02-11 | Prismatic Sensors Ab | Method and devices for image reconstruction |
FR3077908B1 (fr) * | 2018-02-09 | 2021-03-19 | Trixell | Procede d'estimation parametrique d'un objet a estimer dans une image numerique, et procede de suppression de l'objet dans l'image numerique |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003071483A2 (en) * | 2002-02-20 | 2003-08-28 | The Regents Of The University Of Michigan | Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method |
US20040184574A1 (en) * | 2002-07-23 | 2004-09-23 | Xiaoye Wu | Method and apparatus for generating a density map using dual-energy CT |
CN101433464A (zh) * | 2007-11-12 | 2009-05-20 | Ge医疗系统环球技术有限公司 | X射线ct装置和图像处理装置 |
Family Cites Families (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002067201A1 (en) * | 2001-02-15 | 2002-08-29 | The Regents Of The University Of Michigan | Statistically reconstructing an x-ray computed tomography image with beam hardening corrector |
US6507633B1 (en) | 2001-02-15 | 2003-01-14 | The Regents Of The University Of Michigan | Method for statistically reconstructing a polyenergetic X-ray computed tomography image and image reconstructor apparatus utilizing the method |
DE10121441A1 (de) * | 2001-05-02 | 2002-11-07 | Philips Corp Intellectual Pty | Computertomograph |
IL156569A (en) * | 2003-06-22 | 2009-11-18 | Ultraspect Ltd | Improved single-plane plane emission simulation imaging method |
CN100563570C (zh) * | 2004-07-07 | 2009-12-02 | 皇家飞利浦电子股份有限公司 | 心脏锥面光束ct重建中条纹伪影的减少 |
US7684643B2 (en) | 2004-10-26 | 2010-03-23 | Siemens Medical Solutions Usa, Inc. | Mutual information regularized Bayesian framework for multiple image restoration |
JP5091865B2 (ja) * | 2005-09-26 | 2012-12-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像処理システム及び画像処理方法 |
US7983462B2 (en) * | 2005-11-22 | 2011-07-19 | Purdue Research Foundation | Methods and systems for improving quality of an image |
US20070206880A1 (en) | 2005-12-01 | 2007-09-06 | Siemens Corporate Research, Inc. | Coupled Bayesian Framework For Dual Energy Image Registration |
US7760848B2 (en) | 2006-09-08 | 2010-07-20 | General Electric Company | Method and system for generating a multi-spectral image of an object |
US20090116722A1 (en) * | 2007-10-25 | 2009-05-07 | Yunqiang Chen | Method and system for soft tissue image reconstruction in gradient domain |
US8600136B2 (en) * | 2008-09-19 | 2013-12-03 | Koninklijke Philips N.V. | Method for generation of attenuation map in PET-MR |
US8705831B2 (en) * | 2009-09-24 | 2014-04-22 | Koninklijke Philips N.V. | System and method for generating an image of a region of interest |
US8363917B2 (en) * | 2009-10-14 | 2013-01-29 | General Electric Company | System and method of image artifact reduction in fast kVp switching CT |
US8396273B2 (en) * | 2009-12-30 | 2013-03-12 | General Electric Company | Noise reduction method for dual-energy imaging |
DE102010020770A1 (de) * | 2010-05-17 | 2011-11-17 | Siemens Aktiengesellschaft | Verfahren zur Reduzierung von Bildartefakten, insbesondere von Metallartefakten, in CT-Bilddaten |
EP2583241A1 (en) * | 2010-06-21 | 2013-04-24 | Koninklijke Philips Electronics N.V. | Method and system for noise reduction in low dose computed tomography |
CN103229211B (zh) * | 2010-11-25 | 2016-05-11 | 皇家飞利浦电子股份有限公司 | 前向投影设备 |
US8401266B2 (en) * | 2010-11-29 | 2013-03-19 | General Electric Company | Method and system for correlated noise suppression in dual energy imaging |
CN103262124B (zh) * | 2010-12-17 | 2016-11-09 | 皇家飞利浦电子股份有限公司 | 用于对感兴趣区域成像的成像系统 |
US8699812B2 (en) * | 2011-06-09 | 2014-04-15 | General Electric Company | System and method for quality improvement in CT imaging |
CN103620393B (zh) * | 2011-06-21 | 2016-08-24 | 皇家飞利浦有限公司 | 成像装置 |
US8416914B2 (en) * | 2011-07-08 | 2013-04-09 | General Electric Company | System and method of iterative image reconstruction for computed tomography |
CN103717137B (zh) * | 2011-07-28 | 2016-11-23 | 皇家飞利浦有限公司 | 图像生成装置 |
JP6223990B2 (ja) * | 2011-11-23 | 2017-11-01 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | イメージ領域のノイズ除去 |
US9305379B2 (en) * | 2012-01-10 | 2016-04-05 | The Johns Hopkins University | Methods and systems for tomographic reconstruction |
US8923583B2 (en) * | 2012-06-22 | 2014-12-30 | General Electric Company | Methods and systems for performing model-based iterative reconstruction |
KR20140086627A (ko) * | 2012-12-28 | 2014-07-08 | 삼성전자주식회사 | 영상 생성 방법 및 장치 |
US9076237B2 (en) * | 2013-03-12 | 2015-07-07 | Wisconsin Alumni Research Foundation | System and method for estimating a statistical noise map in x-ray imaging applications |
US9269168B2 (en) * | 2013-03-15 | 2016-02-23 | Carestream Health, Inc. | Volume image reconstruction using data from multiple energy spectra |
-
2011
- 2011-12-29 WO PCT/IB2011/056004 patent/WO2012095713A1/en active Application Filing
- 2011-12-29 US US13/978,591 patent/US9165384B2/en active Active
- 2011-12-29 EP EP11813449.3A patent/EP2663964A1/en not_active Withdrawn
- 2011-12-29 CN CN201180064624.4A patent/CN103299345B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003071483A2 (en) * | 2002-02-20 | 2003-08-28 | The Regents Of The University Of Michigan | Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method |
US20040184574A1 (en) * | 2002-07-23 | 2004-09-23 | Xiaoye Wu | Method and apparatus for generating a density map using dual-energy CT |
CN101433464A (zh) * | 2007-11-12 | 2009-05-20 | Ge医疗系统环球技术有限公司 | X射线ct装置和图像处理装置 |
Non-Patent Citations (2)
Title |
---|
JOONKI NOH ET AL.: "Statistical Sinogram Restoration in Dual-Energy CT for PET Attenuation Correction", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》, vol. 28, no. 11, 30 November 2009 (2009-11-30), pages 1688 - 1702, XP011280942, DOI: doi:10.1109/TMI.2009.2018283 * |
顾海峰 等: "双源CT双能量成像及其初步应用初探", 《医疗卫生设备》, vol. 30, no. 6, 30 June 2009 (2009-06-30), pages 63 - 65 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107077748A (zh) * | 2014-10-20 | 2017-08-18 | 皇家飞利浦有限公司 | 针对谱图像迭代重建的开始图像 |
CN107077748B (zh) * | 2014-10-20 | 2020-10-27 | 皇家飞利浦有限公司 | 针对谱图像迭代重建的开始图像 |
CN107209946A (zh) * | 2015-01-20 | 2017-09-26 | 皇家飞利浦有限公司 | 图像数据分割和显示 |
CN107209946B (zh) * | 2015-01-20 | 2020-09-25 | 皇家飞利浦有限公司 | 图像数据分割和显示 |
CN107810518A (zh) * | 2015-06-26 | 2018-03-16 | 皇家飞利浦有限公司 | 利用相关噪声在图像上的边缘探测 |
CN107810518B (zh) * | 2015-06-26 | 2022-06-28 | 皇家飞利浦有限公司 | 图像处理系统和方法 |
CN113039582A (zh) * | 2018-11-16 | 2021-06-25 | 皇家飞利浦有限公司 | 用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 |
Also Published As
Publication number | Publication date |
---|---|
CN103299345B (zh) | 2016-10-19 |
WO2012095713A9 (en) | 2013-05-23 |
WO2012095713A1 (en) | 2012-07-19 |
US9165384B2 (en) | 2015-10-20 |
EP2663964A1 (en) | 2013-11-20 |
US20130287279A1 (en) | 2013-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103299345A (zh) | 双能量断层摄影成像系统 | |
US6754298B2 (en) | Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method | |
CN101681520B (zh) | Pet局部断层摄影 | |
Qi et al. | A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data | |
Lewitt et al. | Overview of methods for image reconstruction from projections in emission computed tomography | |
JP6437652B2 (ja) | 材料特性評価のための断層撮影再構成 | |
EP3050028B1 (en) | Joint reconstruction of electron density images. | |
Frese et al. | Quantitative comparison of FBP, EM, and Bayesian reconstruction algorithms for the IndyPET scanner | |
CN101849247B (zh) | 用于生成衰减图的设备和方法 | |
JP2018140165A (ja) | 医用画像生成装置 | |
US8571291B2 (en) | Combination weight applied to iterative reconstruction in image reconstruction | |
CN107315019A (zh) | 射线透射和荧光ct成像系统和成像方法 | |
CN103180879B (zh) | 用于从投影数据对对象进行混合重建的设备和方法 | |
US9600867B2 (en) | Image processing apparatus and method for filtering an image | |
CN103262124B (zh) | 用于对感兴趣区域成像的成像系统 | |
US20230059132A1 (en) | System and method for deep learning for inverse problems without training data | |
CN101646957B (zh) | 各向同性分辨率的图像重建 | |
CN101622645A (zh) | 用于产生衰减分量的投影系统 | |
CN106133792A (zh) | 图像生成装置 | |
CN108280860A (zh) | 减少多能量ct成像中的光谱通道数目 | |
Elbakri et al. | Statistical x-ray-computed tomography image reconstruction with beam hardening correction | |
Boudjelal et al. | PDEs on graphs for image reconstruction on positron emission tomography | |
CN107251095A (zh) | 图像重建系统、方法和计算机程序 | |
El Khettabi et al. | An inverse problem for three-dimensional x-ray scatter/transmission imaging | |
Chen et al. | A regularized iterative reconstruction algorithm for X-Ray Diffraction 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 |