CN104471615B - 利用正则化的迭代图像重建 - Google Patents
利用正则化的迭代图像重建 Download PDFInfo
- Publication number
- CN104471615B CN104471615B CN201380017943.9A CN201380017943A CN104471615B CN 104471615 B CN104471615 B CN 104471615B CN 201380017943 A CN201380017943 A CN 201380017943A CN 104471615 B CN104471615 B CN 104471615B
- Authority
- CN
- China
- Prior art keywords
- scaling
- regularization term
- projection data
- image
- regularization
- 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
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 description 8
- 230000005855 radiation Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 230000004075 alteration Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011480 coordinate descent method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002574 poison Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
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
-
- 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/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)
- Image Processing (AREA)
- Image Generation (AREA)
Abstract
一种方法,包括利用缩放值缩放利用正则化的迭代重建算法的更新算法的正则化项。所述缩放值在至少一个维度中是可变化的,从而在至少一个维度中使所述迭代重建的所述正则化变化。所述方法还包括至少基于更新算法、变化的经缩放的正则化项以及投影数据来迭代重建图像。
Description
技术领域
以下总体涉及利用正则化的迭代图像重建,并且具体应用于计算机断层摄影(CT)。以下也适用于其他成像模态,例如混合式CT/PET(正电子发射断层摄影)、混合式CT/MR(磁共振)、数字X射线和/或其他成像模态。
背景技术
CT扫描器包括X射线管,其被安装在可旋转的机架上,所述机架围绕检查区域关于z轴旋转。探测器阵列在在检查区域对面与X射线管相对,对向一角度弧。X射线管发射贯穿检查区域的辐射。探测器阵列探测贯穿检查区域的辐射,并且生成指示所述辐射的投影数据。重建器使用迭代或非迭代的重建算法处理投影数据,所述重建算法生成指示检查区域的图像。
利用正则化的迭代图像重建算法一直基于代价函数,其包括数据比较项和图像噪声惩罚项。这样的代价函数的一般公式化为:Ψ(x)=-L(Ax/y)+β·R(x),其中,Ψ(x)表示代价函数,L(Ax/y)表示可能性项,所述可能性项将前向投影图像(Ax,其中,A是前向投影算子,并且x是图像)与测量数据(y)进行比较,R(x)表示粗糙度惩罚项,其在重建图像(x)中惩罚噪声(或“粗糙度”),并且β表示正则化项,其控制正则化的强度。在没有惩罚项的情况下,算法将收敛到非常噪杂的图像以匹配在数据中存在的噪声。
一般,正则化项β是自由参数,并且手动确定其值。能够选择β的值以允许在最终的重建图像中更多或更少的噪声。遗憾的是,β的值趋向于数据集依赖,并且因此,对于每个患者应当被确定,以便针对不同尺寸的患者和剂量水平实现一致的噪声水平。即,如果在对于不同尺寸的解剖结构,和/或具有不同管电流的不同尺寸的患者采集的数据集上使用相同值的β,则将不可能实现在产生的迭代重建图像上的一致的图像质量。
本文描述的各方面解决了以上提到的问题和其他问题。
发明内容
在一个方面,一种方法包括利用缩放值来缩放利用正则化的迭代重建算法的更新算法的正则化项。缩放值在至少一个维度中是可变的,从而在至少一个维度中使迭代重建的正则化变化。所述方法还包括迭代重建至少基于更新算法、变化的缩放正则化项以及投影数据的图像。
在另一方面,一种重建器包括缩放因子算法库和缩放因子生成器,所述缩放因子算法库具有至少一个缩放因子算法,并且所述缩放因子生成器基于至少一个缩放因子算法生成正则化项缩放因子。更新算法库包括采用正则化项的至少一个迭代重建更新算法。重建器处理扫描的投影数据,并且基于至少一个迭代重建算法和正则化项缩放因子来重建图像,其中,缩放因子在至少一个维度中使迭代重建算法的正则化项变化。
在另一方面,一种计算机可读存储介质被编码有计算机可读指令,所述计算机可读指令在由处理器执行时,令所述处理器:重建投影数据,从而使用利用正则化的迭代重建算法生成图像,在所述利用正则化的迭代重建算法中,基于在投影数据中的探测到的光子的数目,迭代重建算法的更新算法的正则化项在至少一个维度中变化。
在另一方面,所述方法包括缩放利用正则化的迭代重建算法的更新算法的正则化项,使得在探测到的光子的总数目或平均数目更高的图像区域中的正则化的缩放抵消更新算法的更强更新值。
本发明可以采取各种部件和部件的布置,以及各种步骤和步骤的布置的形式。附图仅用于图示优选实施例的目的,并不应被解读为对本发明的限制。
附图说明
图1示意性地图示了包括重建器的范例性成像系统,所述重建器采用利用变化的正则化的迭代重建算法以重建图像。
图2示意性地图示了使用利用变化的正则化的迭代重建算法重建图像的范例性重建器。
图3图示了使用利用变化的正则化的迭代重建算法重建图像的范例性方法。
具体实施方式
以下总体涉及利用正则化的迭代图像重建。如以下更加详细地描述,利用本文描述的方法,利用随函数变化的值缩放正则化项β,所述函数重建投影数据的光子数目的反向投影。这样,经缩放的正则化考虑患者尺寸、解剖结构尺寸和/或管电流,其允许重建具有一致的图像质量的图像,所述一致的图像质量不取决于患者尺寸、解剖结构尺寸和/或管电流。
图1图示了示范性成像系统100,例如计算机断层摄影(CT)扫描器。
所述成像系统100包括固定机架102和旋转机架104,所述旋转机架104由固定机架102可旋转地支撑,并且围绕检查区域106关于z轴旋转。
诸如X射线管的辐射源108由旋转机架104可旋转地支撑,并且随旋转机架104旋转,并且发射贯穿检查区域106的辐射。
一维或二维辐射敏感探测器阵列110在检查区域106对面与辐射源108相对,对向一角度弧。探测器阵列110探测贯穿检查区域106的辐射,并且生成指示所述辐射的投影数据。
重建器112处理投影数据,并且重建指示所述投影数据的体积图像数据。在该范例中,重建器112使用诸如利用正则化的迭代重建算法116的来自重建算法库114的一个或多个重建算法来处理投影数据。如以下更详细地描述,在一个实例中,重建器112确定针对利用正则化的迭代重建算法116的正则化项(β)的缩放因子(s),并且在迭代重建期间利用缩放因子来缩放正则化项,其基于投影数据在至少一个维度(例如,一个、两个、三个等维度)使正则化变化。
结果,正则化对患者和/或扫描参数具有特异性,考虑患者尺寸、解剖结构尺寸和/或管电流,并且提供不论患者尺寸、解剖结构尺寸和/或管电流的穿过患者的一致的图像质量(例如,关于图像噪声)。应当理解,用于给定扫描的缩放因子能够从多个预定的缩放因子中选择。例如,在一个实例中,算子可以能够在低、中或高缩放之间选择,取决于感兴趣的正则化水平,其可以是诸如优选算子、患者的尺寸、利用的扫描协议、最初侦察扫描和/或其他信息的信息的函数。
诸如卧榻的对象支撑物118支撑检查区域106中的物体或对象(诸如人或动物对象)。
通用计算系统或计算机充当操作者控制台120。所述控制台120包括诸如监测器的人类可读输出设备和诸如键盘、鼠标等的输入设备。驻留在控制台120上的软件允许操作者经由图形用户界面(GUI)或以其他方式与成像系统100交互和/或操作成像系统100。例如,控制台120允许操作者识别感兴趣的重建算法,例如利用正则化的迭代重建算法116和/或其他算法,以从多个可获得的缩放因子中选择正则化缩放因子等。
应当理解,经由一个或多个处理器(中心处理单元(CPU)、微处理器(μ-CPU)等)能够实施重建器112,所述一个或多个处理器执行嵌入或编码在诸如物理存储器或其他非暂态介质的计算机可读存储介质上的一个或多个计算机可执行指令。额外地或备选地,计算机可执行指令中的至少一个能够由载波、信号或其他瞬态介质运载。
此外,重建器112能够是控制台120、本地计算设备和/或远程(远离系统100,例如,位于不同房间、工厂、州等)计算设备的部分。
图2图示了重建器112的范例。
缩放因子生成器202接收来自探测器阵列110(图1)的投影数据,并且基于接收到的投影数据使用缩放因子算法库206的缩放因子算法204来生成缩放因子s。合适的缩放因子算法包括,当与正则化项β组合时,其使正则化项β变化,使得正则化在至少一个维度变化的缩放因子算法。
在方程1中示出了这样的缩放因子算法的非限制性范例:
方程1:
s=(BP[yi])n,
其中,yi表示投影数据,BP表示反向投影算子,并且n表示正实数。一般,yi是与针对在当前扫描之内的每个整合期中的每个数据点的探测到的光子数目(或信号强度)成比例的数,并且i是在正弦图中的每个数据点上的指数。另一非限制性范例将是多项式函数,
s=a*(BP[yi])k+b*(BP[yi])k-1+…g,其中,k是正整数。本文也预期其他缩放因子算法。
基于缩放因子s、投影数据、初始图像(例如,一切归零、初始反向投影图像、去噪图像等)、停止准则(例如,时限、连续重建之间的不同值等)、用于正则化项β的预定值(例如,默认的、用户定义的、经验确定的等)以及图像更新算法库212的图像更新算法210,图像更新器208迭代重建图像。
在方程2中示出了合适的图像更新算法210的非限制性范例:
方程2:
其中,表示当前图像,表示先前图像,M表示投影子集的数目,BP表示反向投影算子,表示前向投影图像,yi表示投影数据,ai表示沿着线积分i的前向投影算子A的权重的总和,表示惩罚函数(例如,胡贝尔(Huber)惩罚或其他惩罚),表示图像体素梯度(或体素之间的差异),wk表示体素之间的权重(例如,中心到中心的距离或其他距离等),k表示相邻像素(其所有或其子集),并且β表示正则化项。
如在方程3中所示出的,图像更新器208利用方程1缩放方程2中的正则化项β:
方程3:
这意味着,在探测到的光子的总(或平均)数目为高的图像区域中,正则化强度将更强,其抵消在方程3的分子的左侧中给出的更强更新值,并且在探测到的光子的数目为低的图像区域中,将减小正则化强度,在这些区域中再次平衡来自减小的更新的改变。
方程2的更新算法是基于Poison噪声惩罚。对于高斯噪声惩罚,在方程4中示出了合适的图像更新算法210的非限制性实例:
方程4:
其中,fi表示与数据测量结果yi相关联的线积分。如在方程5中所示出的,图像更新器208利用方程1缩放方程4中的正则化项β。
方程5:
利用方程3和方程5,图像更新器208每次迭代同时更新整个图像。图像更新器208每次迭代时能够备选地执行由体素到体素的更新或由体素组到体素组的更新,和/或整个图像的子集的更新。一般,更新算法210能够基于解决一般代价函数Ψ(x)=-L(Ax/y)+β·R(x)的任何迭代重建更新算法,其中,利用缩放因子缩放正则化项β。
在J-B.Thibault、K.D.Sauer、C.A.Bouman、J.Hsieh的“A three-dimensionalstatistical approach to improved image quality for multislice helical CT”(Med.Phys,34(11),2007年11月)中讨论了替代地在体素水平工作的算法的非限制性范例。其中公开的算法被称为迭代坐标下降法或ICD。在J.A.Fessler的“Statistical imagereconstruction methods for transmission tomography”(Handbook of MedicalImaging,卷2:Medical Image Processing and Analysis,第1.70页,SPIE,贝灵汉,2000年)中讨论了另一非限制性范例。
同样地,图像更新器208将利用由缩放因子生成器202生成的缩放因子缩放在这些算法中的正则化项。
图3图示了根据本文描述的实施例的范例性方法。
应当理解,本文描述的方法中的动作的顺序并非是限制性的。这样,本文预期其他顺序。另外,可以省略一个或多个动作和/或可以包括一个或多个额外的动作。
在302处,获得采用正则化的迭代重建更新算法。结合以上方程讨论这样的算法的非限制性范例。
在304处,获得来自扫描的投影数据。
在306处,获得初始图像。如本文所讨论的,初始图像能够是一切归零、滤波反向投影图像、去噪图像和/或其他图像。
在308处,获得正则化项。如本文所讨论的,正则化项能够是默认项、用户指定项等。
在310处,基于所获得的投影数据生成正则化项缩放因子。如本文所描述的,正则化项缩放因子是基于投影数据的光子(或噪声)的数目的反向投影,并且一个非限制性范例包括计算投影数据的反向投影的平方根。
在312处,正则化项缩放因子被并入更新算法。
在314处,至少基于所获得的迭代重建更新算法、所获得的投影数据、所获得的初始图像、所获得的正则化项以及所生成的正则化项缩放因子来重建图像。
以上可以以计算机可读指令的方式实施,所述计算机可读指令被编码或嵌入在计算机可读存储介质上,当(多个)计算机处理器执行所述计算机可读指令时,所述计算机可读指令令(多个)处理器执行所描述的动作。额外地或备选地,计算机可读指令中的至少一个由信号、载波或其他暂态介质运载。
如本文所描述的,缩放使正则化变化,使得在探测到的光子的总数目或平均数目更高的图像区域中,正则化强度更强,并且在探测到的光子的数目更低的图像区域中,正则化强度更不强。
这样,在探测到的光子的总数目或平均数目更高的图像区域中的缩放抵消更新算法的更强更新值,并且在探测到的光子的总数目或平均数目更低的图像区域中的缩放抵消更新算法的更不强的更新值。
这允许重建不取决于患者尺寸、解剖结构尺寸和/或管电流具有一致的图像质量的图像。
已经参考优选实施例描述了本发明。他人在阅读和理解以上具体实施方式的情况下可能想到修改或替代。本文旨在将本发明解释为包括所有这种修改和替代,只要它们落入权利要求书及其等价方案的范围之内。
Claims (19)
1.一种用于利用正则化的迭代图像重建的方法,包括:
根据从CT扫描器接收到的投影数据计算缩放值;
利用所述缩放值缩放利用正则化的迭代重建算法的更新算法的正则化项,
其中,所述缩放值在至少一个维度中是能变化的,从而在所述至少一个维度中使所述迭代重建的正则化项变化;
至少基于所述更新算法、所述变化的经缩放的正则化项以及所述投影数据来迭代重建图像。
2.根据权利要求1所述的方法,还包括:
通过计算所述投影数据的反向投影的平方根来计算所述缩放值。
3.根据权利要求1至2中的任一项所述的方法,其中,在所述投影数据中的噪声与所述投影数据的强度成比例。
4.根据权利要求1至2中的任一项所述的方法,其中,所述更新算法包括常数正则化项,并且还包括在每个更新时将所述常数正则化项乘以缩放值。
5.根据权利要求1至2中的任一项所述的方法,其中,所述缩放值使所述正则化项变化,使得在探测到的光子的总数目或平均数目更高的所述图像的区域中,所述正则化项的正则化强度更强,并且在所述探测到的光子的数目更低的所述图像的区域中,所述正则化项的正则化强度更不强。
6.根据权利要求5所述的方法,其中,在探测到的光子的所述总数目或所述平均数目更高的所述图像的所述区域中的所述正则化项的所述缩放抵消所述更新算法的所述投影数据的更强更新值。
7.根据权利要求5所述的方法,其中,在探测到的光子的所述总数目或所述平均数目更低的所述图像的所述区域中的所述正则化项的所述缩放抵消所述更新算法的所述投影数据的更不强的更新值。
8.根据权利要求1至2中的任一项所述的方法,其中,所述缩放在至少两个维度中使所述正则化项变化。
9.根据权利要求1至2中的任一项所述的方法,其中,所述缩放在至少三个维度中使所述正则化项变化。
10.一种重建器(112),包括:
缩放因子算法库(206),其具有至少一个缩放因子算法(204);
缩放因子生成器(202),其基于所述至少一个缩放因子算法生成正则化项缩放因子并根据从CT扫描器接收到的投影数据来生成所述正则化项缩放因子;
更新算法库(212),其具有采用正则化项的至少一个迭代重建更新算法(210);以及
图像更新器(208),其处理从所述CT扫描器接收到的所述投影数据,并且基于所述至少一个迭代重建更新算法和所述正则化项缩放因子来重建图像,
其中,所述缩放因子在至少一个维度中使所述迭代重建更新算法的正则化项变化。
11.根据权利要求10所述的重建器,其中,所述图像更新器在处理所述投影数据时将所述正则化项乘以所述缩放因子。
12.根据权利要求10-11中的任一项所述的重建器,其中,所述缩放因子生成器将所述正则化项缩放因子生成为所述投影数据的反向投影的平方根。
13.根据权利要求10至11中的任一项所述的重建器,其中,所述投影数据的噪声与所述投影数据的强度成比例。
14.根据权利要求10至11中的任一项所述的重建器,其中,所述缩放因子使所述正则化项变化,使得在探测到的光子的总数目或平均数目更高的所述图像的区域中,所述正则化项的正则化强度更强,并且在所述探测到的光子的数目更低的所述图像的区域中,所述正则化项的正则化强度更不强。
15.根据权利要求14所述的重建器,其中,在探测到的光子的所述总数目或所述平均数目更高的所述图像的所述区域中的正则化项的缩放抵消所述更新算法的所述投影数据的更强更新值,并且在探测到的光子的所述总数目或所述平均数目更低的所述图像的所述区域中的正则化项的缩放抵消所述更新算法的所述投影数据的更不强的更新值。
16.根据权利要求10至11中的任一项所述的重建器,还包括:
控制台(120),其可视地呈现多个缩放因子算法,并且接收指示所述多个缩放因子算法中的用户选择的一个的信号,其中,所述缩放因子生成器基于所述多个缩放因子算法中的所述用户选择的一个来生成所述正则化项缩放因子。
17.一种用于利用正则化的迭代图像重建的方法,包括:
根据从CT扫描器接收到的投影数据计算缩放值;
处理所述投影数据,从而使用利用正则化的迭代重建算法来重建图像,在所述利用正则化的迭代重建算法中,基于所述投影数据中的探测到的光子的数目,所述迭代重建算法的更新算法的正则化项根据所计算的缩放值在至少一个维度中变化。
18.根据权利要求17所述的方法,还包括:
计算所述投影数据的反向投影的平方根;并且
采用所计算的值以在所述至少一个维度中缩放所述正则化项并使其变化。
19.一种用于利用正则化的迭代图像重建的方法,包括:
缩放利用正则化的迭代重建算法的更新算法的正则化项,使得在探测到的光子的总数目或平均数目更高的所述图像的区域中的所述正则化项的缩放抵消所述更新算法的投影数据的更强更新值,其中,所述缩放是根据从CT扫描器接收到的投影数据计算的值的;以及
至少基于所述更新算法、经缩放的正则化项以及所述投影数据来迭代重建图像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261617094P | 2012-03-29 | 2012-03-29 | |
US61/617,094 | 2012-03-29 | ||
PCT/IB2013/052298 WO2013144804A1 (en) | 2012-03-29 | 2013-03-22 | Iterative image reconstruction with regularization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104471615A CN104471615A (zh) | 2015-03-25 |
CN104471615B true CN104471615B (zh) | 2018-05-01 |
Family
ID=48521371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380017943.9A Active CN104471615B (zh) | 2012-03-29 | 2013-03-22 | 利用正则化的迭代图像重建 |
Country Status (6)
Country | Link |
---|---|
US (1) | US9449404B2 (zh) |
EP (1) | EP2831844B1 (zh) |
JP (1) | JP6158910B2 (zh) |
CN (1) | CN104471615B (zh) |
RU (1) | RU2014143489A (zh) |
WO (1) | WO2013144804A1 (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3262612B1 (en) * | 2015-02-25 | 2019-01-30 | Koninklijke Philips N.V. | A method for the reconstruction of quantitative iodine maps using energy resolved tomography |
US10740877B2 (en) * | 2016-04-05 | 2020-08-11 | Koninklijke Philips N.V. | Medical image processing |
JP7159167B2 (ja) * | 2017-01-06 | 2022-10-24 | コーニンクレッカ フィリップス エヌ ヴェ | 陽電子放出断層撮影(pet)撮像における結果のロバスト性を向上させるための標準取込値(suv)誘導再構成制御 |
CN107638189B (zh) * | 2017-10-25 | 2021-06-01 | 东软医疗系统股份有限公司 | Ct成像方法和装置 |
EP3574476A1 (en) | 2018-01-31 | 2019-12-04 | Mitos GmbH | Method for image reconstruction of an object, in particular based on computed-tomography image reconstruction, and apparatus, system and computer program product for the same |
JP6950661B2 (ja) | 2018-10-29 | 2021-10-13 | 株式会社デンソー | 電池監視装置 |
US10685461B1 (en) | 2018-12-20 | 2020-06-16 | Canon Medical Systems Corporation | Apparatus and method for context-oriented iterative reconstruction for computed tomography (CT) |
WO2023114923A1 (en) * | 2021-12-17 | 2023-06-22 | The Regents Of The University Of California | Fusion of deep-learning based image reconstruction with noisy image measurements |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1180876A (zh) * | 1996-10-11 | 1998-05-06 | 通用电气公司 | 计算机x-射线断层透视系统中用于螺旋图像重构的方法和设备 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7272205B2 (en) * | 2004-11-17 | 2007-09-18 | Purdue Research Foundation | Methods, apparatus, and software to facilitate computing the elements of a forward projection matrix |
US7970195B2 (en) * | 2005-03-16 | 2011-06-28 | Koninklijke Philips Electronics N.V. | Method and device for the iterative reconstruction of tomographic images |
WO2006097871A1 (en) | 2005-03-17 | 2006-09-21 | Philips Intellectual Property & Standards Gmbh | Method and device for the iterative reconstruction of cardiac images |
US8538099B2 (en) | 2005-03-23 | 2013-09-17 | General Electric Company | Method and system for controlling image reconstruction |
DE102009014723B4 (de) * | 2009-03-25 | 2012-10-25 | Siemens Aktiengesellschaft | Kontrastabhängige Regularisierungsstärke bei der iterativen Rekonstruktion von CT-Bildern |
US8111893B2 (en) | 2009-06-09 | 2012-02-07 | Wisconsin Alumni Research Foundation | Method for dynamic prior image constrained image reconstruction |
US8705831B2 (en) | 2009-09-24 | 2014-04-22 | Koninklijke Philips N.V. | System and method for generating an image of a region of interest |
US8498465B2 (en) | 2009-09-29 | 2013-07-30 | The Board Of Trustees Of The Leland Stanford Junior University | Accurate determination of the shape and localization of metallic object(s) in X-ray CT imaging |
US8483463B2 (en) * | 2010-05-19 | 2013-07-09 | Wisconsin Alumni Research Foundation | Method for radiation dose reduction using prior image constrained image reconstruction |
DE102010022306A1 (de) | 2010-06-01 | 2011-12-01 | Siemens Aktiengesellschaft | Iterative CT-Bildrekonstruktion in Kombination mit einem vierdimensionalen Rauschfilter |
WO2012017345A1 (en) * | 2010-08-04 | 2012-02-09 | Koninklijke Philips Electronics N.V. | Method and system for iterative image reconstruction |
US8731269B2 (en) * | 2011-10-19 | 2014-05-20 | Kabushiki Kaisha Toshiba | Method and system for substantially reducing artifacts in circular cone beam computer tomography (CT) |
JP6223990B2 (ja) * | 2011-11-23 | 2017-11-01 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | イメージ領域のノイズ除去 |
-
2013
- 2013-03-22 EP EP13725202.9A patent/EP2831844B1/en active Active
- 2013-03-22 RU RU2014143489A patent/RU2014143489A/ru not_active Application Discontinuation
- 2013-03-22 US US14/385,504 patent/US9449404B2/en active Active
- 2013-03-22 CN CN201380017943.9A patent/CN104471615B/zh active Active
- 2013-03-22 JP JP2015502507A patent/JP6158910B2/ja active Active
- 2013-03-22 WO PCT/IB2013/052298 patent/WO2013144804A1/en active Application Filing
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1180876A (zh) * | 1996-10-11 | 1998-05-06 | 通用电气公司 | 计算机x-射线断层透视系统中用于螺旋图像重构的方法和设备 |
Non-Patent Citations (1)
Title |
---|
自适应正则MAP的CT图像重建方法;何玲君等;《计算机工程与应用》;20111231;第47卷(第28期);198-200 * |
Also Published As
Publication number | Publication date |
---|---|
CN104471615A (zh) | 2015-03-25 |
WO2013144804A1 (en) | 2013-10-03 |
EP2831844B1 (en) | 2017-08-23 |
EP2831844A1 (en) | 2015-02-04 |
JP6158910B2 (ja) | 2017-07-05 |
JP2015516831A (ja) | 2015-06-18 |
RU2014143489A (ru) | 2016-05-27 |
US20150049930A1 (en) | 2015-02-19 |
US9449404B2 (en) | 2016-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104471615B (zh) | 利用正则化的迭代图像重建 | |
EP3195265B1 (en) | Iterative image reconstruction with a sharpness driven regularization parameter | |
US9159122B2 (en) | Image domain de-noising | |
EP2992504B1 (en) | De-noised reconstructed image data edge improvement | |
US9600924B2 (en) | Iterative reconstruction of image data in CT | |
EP2691932B1 (en) | Contrast-dependent resolution image | |
CN103907132B (zh) | 图像数据处理 | |
US9025838B2 (en) | Apparatus and method for hybrid reconstruction of an object from projection data | |
EP3270788B1 (en) | System and method for computed tomography | |
JP2013543784A (ja) | 順投影装置 | |
US9299172B2 (en) | Image generation apparatus | |
EP2646985A1 (en) | Iterative reconstruction algorithm with a constant variance based weighting factor | |
Li et al. | Multienergy cone-beam computed tomography reconstruction with a spatial spectral nonlocal means algorithm | |
EP3552181B1 (en) | Image noise estimation using alternating negation | |
WO2014170780A1 (en) | Iterative image reconstruction with tissue dependent regularization | |
Anthoine et al. | On the efficiency of proximal methods in CBCT and PET |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |