用于谱ct的图像处理的制作方法
【专利摘要】一种方法，包括基于噪声模型针对谱图像的（一个或多个）体素估计结构模型；将结构模型拟合到关于所述（一个或多个）体素的3D邻域；基于所述拟合和预定的模型选择标准来选择针对所述（一个或多个）体素的所述结构模型中的一个；以及基于选择的结构模型来对所述（一个或多个）体素进行降噪，产生一组降噪谱图像。另一种方法，包括基于降噪谱图像、经分解的降噪谱图像、碘图和对比度增强因子，针对对应于不同能量范围的一组谱图像中的每个能量图像生成虚拟对比度增强中间图像；以及通过模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟对比度增强图像。本文中还描述了用于生成虚拟非对比度图像、骨与钙化分割图以及碘图的途径，用于多能量成像研究。
【专利说明】用于谱CT的图像处理
【技术领域】
[0001]下文总体上涉及计算机断层摄影(CT)，并且更具体地，涉及谱CT。
【背景技术】
[0002]CT扫描器一般包括发射电离辐射的X射线管，所述电离辐射穿过检查区域及其中的目标或对象的一部分，并照射在所述检查区域间设置的与所述X射线管相对的探测器阵列。所述探测器产生指示探测到的辐射的投影数据。能够重建所述数据，以生成指示所述目标或对象的所述部分的体积图像数据。对于谱CT，所述投影数据包括被同时采集并且对应于不同光子能量范围的信号。存在若干种途径用于执行谱CT。例如，所述CT扫描器可以包括两个或多个源，至少一个源被配置为在至少两个不同kVp之间切换，和/或具有能量分辨率探测器的探测器阵列。
[0003]对于谱CT，两种采集的信号能够用于确定每个信号的光电贡献和康普顿贡献，并且通过其光电贡献和康普顿贡献的值来识别未知的材料。通常，由于两个基函数的任何两个线性独立的总和跨越整个衰减系数空间，任何材料均能够由两种基础材料的线性组合来表示。这在诸如碘的材料中尤其有效，这样的材料具有接近于诊断能量范围的平均值的k缘能量。此外，额外的谱信息改善了关于被扫描目标及其材料组成能够被确定的定量信息。所述基础材料也允许生成单色图像、材料取消图像、有效原子数图像以及电子密度图像。
[0004]再者，CT扫描器发射电离辐射。遗憾的是，电离辐射可以损伤或杀死细胞和/或增加癌症的风险。文献已指出，来自CT的剂量水平通常超过来自常规放射照相和荧光检查的剂量水平。然而，不能只将针对特定成像过程的辐射剂量降低，这是因为较低的剂量导致增加的图像噪声以及因此更模糊或不清晰的图像。此外，谱CT图像已固有地比常规非谱图像的噪声更大。例如，在双能研究中，每个图像均大致基于对应的非谱常规扫描的辐射剂量的一半。此外，对材料分解的估计是基于在其之间具有窄角度的两个向量之间的投影。这两个因素(即，大噪声与窄角度)的组合，显著放大了经估计的材料分解中的噪声。
[0005]对比度增强CT研究捕获所给予的放射对比度材料通过脉管组织的输送。大体上，针对对比度增强CT，将放射对比度材料团经静脉给予患者，并且扫描包括感兴趣脉管组织的所述患者中的感兴趣区域。随着所述放射对比度材料流过所述脉管组织，所述放射对比度材料引起所述感兴趣脉管组织中的X射线密度暂时增加，得到增强的数据。然而，在给予对比度材料后，一些患者经历特殊效应，并且某些患者会经历严重且有潜在生命威胁的过敏反应。对比度材料也可能诱导肾损伤，并且一些患者已发展为其肾功能的急性恶化。通常，较大的对比度材料体积导致较高对比度噪声比(CNR)的图像，而较低的体积导致较低CNR的图像。遗憾的是，随着所述对比度材料体积增加，与其相关联的风险也增加。
【发明内容】
[0006]本申请的各方面解决上述问题以及其他问题。
[0007]根据一个方面，一种方法，包括针对对应于不同能量范围的一组谱图像中的谱图像的一个或多个体素估计局部噪声值；产生针对所述图像的噪声模型；基于对应的噪声模型，针对所述谱图像的体素估计局部结构模型；将一组所述局部结构模型拟合到关于所述图像中的体素的所述图像中的三维体素邻域；基于所述拟合和预定的模型选择标准，选择针对所述体素的所述局部结构模型中的一个；以及，通过利用基于选择的局部结构模型来估计的值代替所述体素的值，基于选择的局部结构模型，对所述体素进行降噪，其中，在该组谱图像的多个谱图像的多个所述体素被降噪，产生一组降噪谱图像。
[0008]在另一方面中，一种计算装置，包括噪声估计器，其估计对应于不同能量范围的一组谱图像中的谱图像的噪声模式，其中，所述噪声模式用于针对所述谱图像的体素估计局部结构模型；模型拟合器，其将一组所述局部结构模型拟合到关于所述图像中的体素的所述图像中的三维体素邻域；以及，模型选择器，其基于所述拟合和预定的模型选择标准，选择针对所述体素的所述局部结构模型中的一个。
[0009]在另一方面中，一种方法，包括基于降噪谱图像的概率分解，生成钙概率图；通过执行所述钙概率图的全变差函数最小化，增强所述钙概率图；以及，基于增强的钙概率图和预定的阈值，生成表不骨与I丐分割的二值掩模。
[0010]在另一方面中，一种方法，包括基于降噪谱图像的向量分解，生成一幅或多幅碘分布图；以及，基于所述一幅或多幅碘分布图和表示所述骨与钙分割的二值掩模，估计碘图。
[0011]在另一方面中，一种方法，包括基于降噪谱图像、经分解的降噪谱图像、碘图和对比度增强因子，针对对应于不同能量范围的一组谱图像中的每幅能量图像，生成虚拟对比度增强中间图像；以及通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟对比度增强图像。
[0012]在另一方面中，一种方法，包括基于所述降噪谱图像、经分解的降噪谱图像和碘图，针对对应于不同能量范围的一组谱图像中的每个能量图像，生成虚拟非对比度中间图像；以及通过将模拟的部分体积效应与所述中间虚拟对比度增强图像，生成最终虚拟非对比度图像。
[0013]本领域普通技术人员在阅读和理解以下详细描述时，将认识到本发明其他的方面。
【专利附图】
【附图说明】
[0014]本发明可以采取各种部件与部件的布置，以及各个步骤与步骤的安排。附图仅是出于图示优选实施例的目的，并且不应被解读为限制本发明。
[0015]图1示意性地图示了与降噪器和图像处理器连接的成像系统。
[0016]图2示意性地图示了降噪器的范例。
[0017]图3示意性地图示了降噪器的谱噪声去除器的范例。
[0018]图4示意性地图示了图像处理器的范例。
[0019]图5示出了双能研究的能量图/能量散点图以及若干材料响应向量的范例。
[0020]图6示出了能量图中的两个材料响应向量以及从测量点到所述两个向量的较短距离。
[0021]图7示意性地图示了图像处理器的范例分割器。
[0022]图8示意性地图示了图像处理器的范例材料图生成器。[0023]图9示意性地图示了图像处理器的范例虚拟对比度增强图像生成器。
[0024]图10示意性地图示了图像处理器的范例虚拟非对比度图像生成器。
[0025]图11图示了用于对谱图像进行降噪的范例方法。
[0026]图12图示了用于确定针对降噪谱图像的骨与钙分割二值掩模的范例方法。
[0027]图13图示了用于确定针对降噪谱图像的碘图的范例方法。
[0028]图14图示了用于基于降噪谱图像来确定虚拟非对比度图像的范例方法。
[0029]图15图示了用于基于降噪谱图像来确定虚拟对比度增强图像的范例方法。
【具体实施方式】
[0030]图1图示诸如计算机断层摄影(CT)扫描器的成像系统100，其被配置为用于谱CT成像。成像系统100包括固定机架102和由固定机架102可旋转地支撑的旋转机架104。旋转机架104关于纵轴或z轴围绕检查区域106旋转。
[0031]系统100包括至少一个诸如X射线管的辐射源108，其由旋转机架104支撑，并且连同旋转机架104 —起关于检查区域106旋转。至少一个福射源108发射穿过检查区域106的辐射。在具有至少两个辐射源108时，每个源能够被配置为发射具有不同平均发射谱的辐射。额外地或备选地，至少两个辐射源108中的一个或多个能够被配置为在扫描期间可控地在至少两个不同发射电压(kVp)之间切换。多源和/或多kVp切换能够用于谱CT米集。
[0032]辐射敏感探测器阵列110位于跨检查区域106与至少一个辐射源108相对。辐射敏感探测器阵列110包括探测器像素的阵列，所述探测器像素探测穿过检查区域106的辐射，并生成指示所述辐射的投影数据。辐射敏感探测器阵列110能够包括常规探测器和/或能量分辨率探测器，诸如，直接转换探测器和/或基于闪烁体的多谱探测器，所述多谱探测器包括具有不同X射线能量灵敏度的至少两个闪烁体，所述闪烁体分别以光学方式被附接到具有对应的光学灵敏度的至少两个光传感器(例如双层探测器)。所述能量分辨率探测器能够用于谱CT采集。
[0033]重建器112重建所述投影数据，并生成指示检查区域106及其中的目标或对象的所述部分的体积图像数据。在采集谱数据时(例如，在所述投影数据包括经由多个源、kVp切换和/或能量分辨率探测器，同时采集并且对应于不同能量范围的至少两个测量结果时)，重建器112能够针对所述不同能量范围中的每个重建个体谱图像和/或基于对应于所述不同能量范围中的两个或多个的个体谱图像重建组合图像。重建器112也能够采用常规的非谱重建算法。
[0034]诸如卧榻的对象支撑体114支撑检查区域106中的对象(例如人类或动物)或目标，并且能够用于在扫描之前、期间和/或之后，相对于X轴、y轴和/或Z轴和检查106定位所述对象。通用计算系统充当操作者控制台116，并且包括诸如显示器的输出设备和诸如键盘、鼠标和/或其他的输入设备。驻留在控制台116上的软件允许所述操作者控制系统100的运行，例如，允许所述操作者选择谱成像协议、启动扫描等。
[0035]计算装置118包括一个或多个处理器，所述处理器执行被嵌入或编码在诸如物理存储器的计算机可读存储介质上的一个或多个计算机可执行指令。额外地或备选地，所述计算机可执行指令中的一个或多个能够由信号或载波来承载，并由所述一个或多个处理器来执行。在图示的实施例中，所述计算机可执行指令包括用于实施降噪器120和/或图像处理器122的指令。在另一实施例中，经由控制台116和/或其他设备实施降噪器120和/或图像处理器122。
[0036]降噪器120被配置为对来自重建谱图像的谱噪声进行降噪，去除或减少来自所述重建谱图像的谱噪声，并产生降噪重建谱图像。如下文更详细描述的，在一种情况中，降噪器120去除或减少谱噪声，同时保留基本的谱信息和目标结构。在一种情况中，这允许针对给定图像品质的辐射剂量减少。备选地，针对给定剂量，能够增强图像品质。备选地，能够实现剂量减少与图像品质增强的组合。额外地或备选地，降噪器120能够对估计的单色图像进行降噪，使用光电分量与康普顿分量的适当组合，能够模拟所述估计的单色图像，以估计任何keV图像。
[0037]图像处理器122处理所述降噪重建谱图像和/或所述单色图像。如下文更详细描述的，这包括以下中的一个或多个:执行骨和/或钙分割；创建研究中碘的定量分布的碘图；生成虚拟对比度增强(VCE)图像和/或生成虚拟非对比度(VNC)图像。这样的骨与钙化分割能够高度利用额外的定量谱信息，在束硬化校正算法内利用这样的骨与钙化分割，与单色图像重建算法一起利用这样的骨与钙化分割等。所述碘图提供在研究中碘的改善的定量分布。
[0038]虚拟增强对比度允许针对给定图像品质减少给予患者的造影剂的量。备选地，其允许节省在其中已漏掉了来自给药的扫描计时的研究，并且得到的图像具有次优的图像品质，这会导致重复扫描以及更多的造影剂。备选地，其允许临床医师经由鼠标、键盘等手动微调图像处理参数，以实时调查图像并获得期望的可视化结果。所述VNC图像可以消除对非对比度扫描的需要，其能够减少辐射暴露，节约时间，并且延长管寿命。
[0039]计算装置118也包括用户接口 124，其允许用户与计算装置118交互。在一种情况中，这包括允许临床医师选择上面提到的图像处理特征(即骨与钙分割、碘图生成、VNC图像生成和/或VCE图像生成)中的哪个，以用于给定研究。用户接口 124也允许所述临床医师设定和/或改变各个图像处理参数。例如，所述临床医师能够使用用户接口 124来改变针对所述研究降噪的量。这能够利用实时呈现的结果实时动态地完成。
[0040]亦即，在观看结果时，临床医师能够改变参数，其作为响应，调用所述计算装置，以基于改变的参数处理所述降噪重建谱图像，并视觉呈现结果。可以是用户可配置的其他参数包括，但不限于，用于VCE图像处理的对比度增强因子、影响用于VCE和VNC图像生成的模拟部分体积效应的侵蚀性的参数、用于选择用于降噪的去局部结构模型的阈值、针对骨与钙分割的比例因子、用于将局部结构模型拟合到体素的加权因子等。
[0041]数据储存库126能够用于存储所述重建图像、降噪重建图像和/或经处理的重建图像和/或降噪重建图像，并且能够由控制台116、图像处理器122、降噪器120和/或其他设备中的一个或多个访问。数据储存库126能够在系统100本地、远离系统100、为分布式等。数据储存库126可以包括数据库、服务器、图片归档与通信系统(PACS)、放射科信息系统(RIS)、医院信息系统(HIS)、电子病历(EMR)和/或其他电子存储设备或存储器。
[0042]图2示意性地图示了降噪器120的范例。大体上，在该实施例中，降噪器120被配置为针对对应于不同能量范围的一组谱图像中的谱图像确定噪声模式，以及基于所述噪声模式来减少所述谱图像的谱噪声。图示的降噪器120接收一组谱图像作为输入，其可以包括来自扫描器100、储存库126和/或其他位置的一组重建谱图像，和/或一组估计的单色图像。
[0043]噪声估计器202估计谱图像的每个体素的局部噪声值，并且基于所述局部噪声值来生成针对所述谱图像的噪声模型或模式，并且所述噪声模型用于估计所述谱图像中的结构。噪声估计器202能够使用已知的和/或其他途径来估计噪声。合适的途径包括，但不限于，蒙特卡洛估计；诸如在 Wunderlich 和 Noo 的 Phys.Med.Biol.53(2008 年)，2472 至 2493页中讨论的分析途径；基于图像的途径，诸如在2010年10月29日提交的题为“ENHANCEDIMAGE DATA/DOSE REDUCTION” 的 PCT 专利申请序列号 PCT/IB2009/054913 中描述的途径，在此通过引用将其全文并入；和/或其他途径。
[0044]谱噪声去除器204基于估计的噪声模型，从所述谱图像去除谱图像噪声，生成降噪谱图像，同时保留不同能量图像中的基本的谱信息和/或解剖结构，从而改善所述谱图像的信噪比。联系图3描述了这样的范例，在该范例中，谱噪声去除器204包括模型拟合器302，针对所述谱图像中的所述体素中的一个或多个，所述模型拟合器302将基于所述估计噪声模式确定的局部结构模型拟合到关于体素的三维体素区域或邻域。
[0045]谱噪声去除器204也包括模型选择器304，其基于被存储在标准存储器306中的预定的选择标准来选择针对图像中的每个体素的局部结构模型。一旦针对每个体素来选择模型，由谱噪声去除器204利用所述模型来对所述谱图像中的噪声进行降噪或去除，其中，新估计的体素的值是由选择的模型确定的值，并且代替所述体素的原始值。得到的谱图像包括降噪谱图像，或针对所述不同能量范围的图像品质增强的谱图像。
[0046]参考图1、 图2和图3，针对降噪器120的范例噪声去除途径由此得到。对于该范
例，表示体中的体素，其中，通过能量氏获得所述体积。等式I包括能够用于拟
合局部结构模型的最小二乘法:
[0047]等式1:
[0048]PEe =argmin^ J Σ{p)) ir:..r.k.P i'=-n j'~-n k'=—n，
[0049]其中，Mlr k,(p)为针对所述体积中的(i+i'，j+j'，k+k')体素的模型值，并且
Wi, rJ, 为权重因子。所述权重因子能够被视为局部化核心，其为两个权重函数的乘积，如等式2中所示:
[0050]等式2:
[0051]
空间HU
WiIjjc’ = W i\j\k'W i\j\k'，
[0052]其中，表示根据邻域到所述体素的空间距离，对所述邻域的权重；并且wHUrj'x表示根据邻域在亨氏单位(HU)空间中到所述体素的强度距离，对所述邻域的权重。
[0053]基于等式3能够确定所述W2Ivx函数:
[0054]等式3:
[0055]
【权利要求】
1.一种方法，包括: 针对对应于不同能量范围的一组谱图像中的谱图像的一个或多个体素估计局部噪声值，产生针对所述谱图像的噪声模型； 基于对应的噪声模型来估计针对所述谱图像的体素的局部结构模型； 基于预定的模型选择标准来选择针对所述体素的所述局部结构模型中的一个；以及通过用基于选择的局部结构模型估计的值代替所述体素的值，基于选择的局部结构模型对所述体素进行降噪， 其中，对该组谱图像中的多幅谱图像的多个所述体素进行降噪，产生一组降噪谱图像。
2.如权利要求1所述的方法，还包括，在所述的基于预定的模型选择标准来选择针对所述体素的所述局部结构模型中的一个之前: 将一组所述局部结构模型拟合到关于所述图像中的体素在所述图像中的三维体素邻域；以及 基于所述拟合和所述预定的模型选择标准来选择针对所述体素的所述局部结构模型中的所述一个。
3.如权利要求1至2中任一项所述的方法，其中，利用在第一剂量扫描期间采集的数据生成所述谱图像，并且对所述谱图像的所述降噪产生具有与利用在第二剂量扫描期间采集的数据生成的谱图像处于相同的图像噪声水平的图像噪声的量的一组图像，其中，所述第二剂量扫描的剂量高于所述第一剂量扫描的剂量。
4.如权利要求1至3中任一项所述的方法，还包括: 使用最小二乘最小化来拟合所述模型；以及 利用加权因子对所述最小二乘最小化进行加权，其中，所述加权因子包括第一加权分量，其基于相邻的体素与所述体素之间的体素强度距离，对关于所述体素的所述三维体素邻域进行加权。
5.如权利要求4所述的方法，其中，所述第一加权分量为所述体素的所述局部噪声值的函数。
6.如权利要求4至5中任一项所述的方法，其中，所述加权因子权重包括第二加权分量，其基于所述相邻的体素和所述体素之间的空间距离，对关于所述体素的所述三维邻域进行加权。
7.如权利要求1至6中任一项所述的方法，其中，所述局部结构模型包括至少两个噪声模型，其中，所述至少两个局部结构模型至少包括对均质区域进行建模的常数模型和对非均质区域进行建模的二阶多项式。
8.如权利要求6至7中任一项所述的方法，还包括: 基于所述至少两个局部结构模型的局部标准差的比率与预定的阈值之间的预定的关系，选择经拟合的噪声模型的局部结构模型，以对所述体素进行降噪。
9.如权利要求1至8中任一项所述的方法，还包括: 基于对所述降噪谱图像的概率分解，生成钙概率图； 通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及 基于增强的钙概率图和预定的阈值，生成表示骨与钙分割的二值掩模。
10.如权利要求9所述的方法，还包括:基于对所述降噪谱图像的向量分解，生成一幅或多幅碘分布图；以及 基于所述一幅或多幅碘分布图和所述二值掩模来估计碘图。
11.如权利要求10所述的方法，还包括: 基于所述降噪谱图像、经分解的降噪谱图像、所述碘图和对比度增强因子，生成针对每个能量的虚拟对比度增强中间图像；以及 通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟对比度增强图像。
12.如权利要求11所述的方法，其中，所述对比度增强因子具有补偿造影剂的减少的值。
13.如权利要求10所述的方法，还包括: 基于所述降噪谱图像、经分解的降噪谱图像和所述碘图，生成针对每个能量的虚拟非对比度中间图像；以及 通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟非对比度图像。
14.一种计算装置(118)，包括: 噪声估计器(202)，其估计对应于不同能量范围的一组谱图像的谱图像的噪声模式，其中，所述噪声模式用于估计针对所述谱图像的体素的局部结构模型；以及 模型选择器(304)，其基于预定的模型选择标准来选择针对所述体素的所述局部结构模型中的一个。
15.如权利要求14所述的计算装置，还包括: 模型拟合器(302)，其将一组所述局部结构模型拟合到关于所述图像中的体素在所述图像中的三维体素邻域，其中，所述模型选择器基于所述拟合和所述预定的模型选择标准来选择针对所述体素的所述局部结构模型中的所述一个。
16.如权利要求14至15中任一项所述的计算装置,还包括: 谱噪声去除器(204)，其通过利用基于选择的局部结构模型估计的值代替所述体素的值，基于所述选择的局部结构模型来对所述体素进行降噪，其中，对该组谱图像中的多幅谱图像的多个体素进行降噪，产生一组降噪谱图像。
17.如权利要求14至16中任一项所述的计算装置，还包括: 材料分析器(402)，其基于对所述降噪谱图像的概率分解，生成钙概率图； 增强器(704)，其通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及 掩模估计器(706)，其基于增强的钙概率图和预定的阈值，生成表示骨与钙分割的二值掩模。
18.如权利要求17所述的计算装置，其中，所述材料分析器基于对所述降噪谱图像的向量分解，生成一幅或多幅碘分布图，并且，还包括: 碘图估计器(802)，其基于所述一幅或多幅碘分布图和所述二值掩模来估计碘图。
19.如权利要求18所述的计算装置，还包括: 中间虚拟对比度增强图像生成器(902)，其基于所述降噪谱图像、经分解的降噪谱图像、所述碘图和对比度增强因子，生成针对每个能量的虚拟对比度增强中间图像；以及最终虚拟对比度增强图像估计器(904)，其通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟对比度增强图像。
20.如权利要求18所述的计算装置，还包括: 中间虚拟非对比度图像生成器(1002)，其基于所述降噪谱图像、经分解的降噪谱图像和所述碘图，生成针对每个能量的虚拟非对比度中间图像；以及 最终虚拟非对比度图像生成器(1004)，其通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟非对比度图像。
21.—种方法,包括: 基于对降噪谱图像的概率分解，生成钙概率图； 通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及 基于增强的钙概率图和预定的阈值，生成表示骨与钙分割的二值掩模。
22.—种方法,包括: 基于对降噪谱图像的向量分解，生成一幅或多幅碘分布图；以及 基于所述一幅或多幅碘分布图和表示骨与钙分割的二值掩模，估计碘图。
23.如权利要求21所述的方法，还包括: 基于对所述降噪谱图像的概率分解，生成钙概率图； 通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及 基于增强的钙概率图和预定的阈值，生成所述二值掩模。
24.—种方法,包括: 基于降噪谱图像、经分解的降噪谱图像、碘图和对比度增强因子，生成针对对应于不同能量范围的一组谱图像的每个能量图像的虚拟对比度增强中间图像；以及 通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟对比度增强图像。
25.如权利要求23所述的方法，还包括: 基于对所述降噪谱图像的向量分解，生成一幅或多幅碘分布图；以及 基于所述一幅或多幅碘分布图和表示所述骨与钙分割的二值掩模来估计所述碘图。
26.如权利要求24所述的方法，还包括: 基于对所述降噪谱图像的概率分解，生成钙概率图； 通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及 基于增强的钙概率图和预定的阈值，生成所述二值掩模。
27.—种方法,包括: 基于降噪谱图像、经分解的降噪谱图像和碘图，生成对应于不同能量范围的一组谱图像的每个能量图像的虚拟非对比度中间图像；以及 通过将模拟的部分体积效应与所述中间虚拟对比度增强图像进行合并，生成最终虚拟非对比度图像。
28.如权利要求26所述的方法，还包括: 基于对所述降噪谱图像的向量分解，生成一幅或多幅碘分布图；以及 基于所述一幅或多幅碘分布图和表示骨与钙分割的二值掩模来估计所述碘图。
29.如权利要求27所述的方法，还包括:基于对所述降噪谱图像的概率分解，生成钙概率图； 通过执行对所述钙概率图的全变差函数变化最小化，增强所述钙概率图；以及基于增强的钙概率图和预定的阈值，生成所述二值掩模。
【文档编号】G06T5/00GK103649990SQ201280034976
【公开日】2014年3月19日 申请日期:2012年7月10日 优先权日:2011年7月15日
【发明者】L·戈申 申请人:皇家飞利浦有限公司