专利名称：边缘保护的噪声滤波的制作方法
技术领域：
本发明涉及一种用于处理图像的图像处理装置、图像处理方法和图像处理计算机程序。
背景技术：
C. K. Chu 等人的文章 “Edge-Preserving Smoothers for ImageProcessing，，, Journal of The American Statistical Association,第 93 卷，第 442 期，第526-541页，1998年6月公开了一种滤波方法，该滤波方法适于执行对图像的滤波，从而在保护边缘的同时降低不期望的噪声。然而，如果在图像中CNR (对比噪声比)过小，那么边缘可能看不到，而且可能无法保护，导致图像的模糊，并且因此导致降低的图像质量
发明内容
因此，本发明的目的在于提供用于处理图像的图像处理装置、图像处理方法和图像处理计算机程序，其中，能够改进处理后的图像的图像质量。在本发明的第一方面中，提出了一种用于处理图像的图像处理装置，其中该图像处理装置包括-用于提供对象的第一图像以及同一对象的第二图像的图像提供单元，-用于对第一图像进行滤波的滤波单元，其中，滤波单元适于根据第一图像的图像值之间的第一相似度以及根据第二图像的图像值之间的第二相似度对第一图像进行滤波。根据图像值之间的相似度的滤波允许根据图像值属于对象的同一部分(例如，如果对象是人或动物的话，属于对象的组织或骨骼物质)的可能性(Iikeliness)对第一图像进行滤波。例如，如果第一图像的或者第二图像的图像值之间的相似度大，那么对应的图像值就很可能属于对象的同一部分。相反，如果第一图像的和第二图像的图像值之间的相似度较小，那么它们更可能属于对象的不同部分，并且它们例如被对象的边缘隔开。即使由于噪声或其他干扰，第一图像的图像值比起在没有噪声或其他干扰时它们所应当的那样更加相似，在滤波期间也可以考虑图像值示出对象的同一部分的可能性，这是因为滤波单元适于不仅根据第一图像的图像值之间的第一相似度还根据第二图像的图像值之间的相似度来对第一图像进行滤波，该第二图像可能未被与第一图像相同的程度地干扰或者可能根本未被干扰。这允许提高第一图像的图像质量。图像提供单元可以是已经存储了第一图像和第二图像的存储单元，或者图像提供单元可以是诸如计算机断层扫描系统、磁共振成像系统、核成像系统、超声成像系统等的成像模态。如果图像提供单元是成像模态，那么也可以将图像处理装置看作是生成对象的图像并且对这些图像进行处理的成像系统。滤波单元优先地还适于根据第一图像的图像值之间的第一相似度并且根据第二图像的图像值之间的第二相似度对第二图像进行滤波。第一图像和第二图像可以是二维图像、三维图像或四维图像。它们示出同一对象，即，整个对象或者仅对象的一部分，例如，仅对象的切片，如果对象在运动，则优选的是在同一时间和/或在同一运动状态下的。优先基于采集的对象的同一数据(尤其是投影数据)来重建第一图像和第二图像。滤波单元优先适于根据所述第一图像的图像值之间的差异来确定第一相似度，并且根据第二图像的图像值之间的差异来确定第二相似度，其中随着差异的增大，尤其是随着绝对差异的增大，相似度降低。而且，可以采用其他的相似性度量来用于确定图像值之间的相似度，例如，图像值的比值，其中，随着该比值与I的偏差的缩小，相似度提高。滤波单元优先地适于执行滤波以致第一图像和/或第二图像被降噪。优选地，滤波单元适于对第一图像的像元的图像值进行加权平均，其中，为了确定第一 图像的第一像元的平均图像值，对包括第一像元的第一区域内的第一区域像元的图像值进行加权并且平均，其中，采用根据下述内容的权重对第一区域像元的图像值进行加权-作为a)第一区域像元的图像值与b)第一图像的第一像元的图像值之间的相似度的第一相似度-作为a)第二图像中的与第一区域对应的第二区域内的第二区域像元的图像值与b)第二图像中的与第一像元对应的第二像元的图像值之间的相似度的第二相似度。这允许通过简单的数学运算并且因而以低计算工作量对第一图像进行滤波。进一步优选的是，调整滤波单元使得如果a)第一区域像元与第一像元之间的第一相似度和b)第二区域像元与第二像元之间的第二相似度中的一个减小，并且如果a)第一区域像元与第一像元之间的第一相似度和b)第二区域像元与第二像元之间的第二相似度中的另一个保持不变或者减小，则第一区域像元的图像值的权重减小。还可以调整滤波单元使得如果a)第一区域像元与第一像元之间的第一相似度和b)第二区域像元与第二像元之间的第二相似度中的一个减小，并且如果a)第一区域像元与第一像元之间的第一相似度和b)第二区域像元与第二像元之间的第二相似度中的另一个增大，则第一区域像元的图像值的权重减小。因此，如果第一相似度和第二相似度中的至少一个减小，则优选地，第一区域元素的图像值的权重减小。这允许以简单的方式考虑比其他图像值强的在图像的至少一幅中的具有较大的相似度的图像值。因此，即使由于噪声或其他干扰的原因，相似度在图像之一内比没有噪声或其他干扰的情况下所应当的相似度更大，相应的图像值也能够接受小的权重。进一步优选的是，第一相似度是a)第一区域像元的图像值和围绕相应的第一区域像元的像元的图像值与b)第一像元的图像值和围绕第一图像的第一像元的像元的图像值之间的相似度，并且其中，第二相似度是a)第二图像中的与第一区域对应的第二区域内的第二区域像元的图像值和围绕相应的第二区域像元的像元的图像值与b)第二图像中的与第一像元对应的第二像元的图像值和围绕第二像元的像元的图像值之间的相似度。因而，为了确定相应的相似度，不仅考虑图像值的对，还考虑图像值的群组，其中，在这些不同的图像值群组之间确定相似度。这允许例如，在确定相应的相似度的同时考虑不同的图像值群组中所示的结构。进一步优选的是，如果第一图像将被降噪，则第一相似度是第一图像的图像值相似的概率，并且其中，如果第二图像将被降噪，则第二相似度是第二图像的图像值相似的概率。进一步优选的是，滤波单元适于基于图像值之间的差异来确定图像值相似的概率。这允许如果相应的图像将被降噪，则直接根据图像值将相似的概率来对该图像值进行滤波。进一步优选的是，滤波单元适于根据第一图像和第二图像中的至少一个中的像元之间的距离对第一图像进行滤波。因为如果不同的像元的距离增大，那么它们属于对象的不同部分的可能性也增大，因而通过根据第一图像和第二图像中的至少一个中的像元之间的距离对第一图像进行滤波，能够进一步改进滤波，例如使得边缘得到保护，由此进一步提高了经处理的图像的质量。进一步优选的是，图像提供单元适于提供第一图像和第二图像使得它们示出对象的不同特性。进一步优选的是，图像提供单元适于提供谱计算机断层摄影图像，其中，第一图像指示对象的第一物理特性，而第二图像指示对象的第二物理特性，第一物理特性和第二物理特性互不相同。进一步优选的是，第一图像和第二图像之一是康普顿(Compton)图像，而第一图像和第二图像中的另一个是光电图像。进一步优选的是，第一图像和第二图像之一是康普顿图像或光电图像，而第一图像和第二图像中的另一个是康普顿图像和光电图像的组合。图像提供单元还可以适于提供微分相位衬度计算机断层摄影图像，其中，第一图像和第二图像之一指示对象的吸收，而第一图像和第二图像中的另一个指示由对象诱发的 横贯对象的辐射的相移。第一图像能够显示对象的总衰减系数，而第二图像能够显示对象的折射率(refractive index)的实部。如果所述第一图像和第二图像表现对象的不同特性，则很可能的是，如果在第一图像中的某一区域内存在噪声，在第二图像的对应区域内可能不存在相同程度的噪声或者可能根本不存在噪声。因而，如果例如，由于噪声的原因在第一图像的某一区域内看不到边缘，那么在第二图像的对应区域内可能看得到该边缘，并且能够在滤波过程期间通过根据第一图像的图像值之间的第一相似度和第二图像的图像值之间的第二相似度对第一图像和/或第二图像进行滤波来考虑该边缘。就本发明的另一方面而言，提出了一种用于处理图像的图像处理方法，其中，该图像处理方法包括-提供对象的第一图像和同一对象的第二图像，-根据第一图像的图像值之间的第一相似度以及根据第二图像的图像值之间的第二相似度对第一图像进行滤波。在本发明的另外的方面中，提出了一种用于处理图像的图像处理计算机程序，其中，该计算机程序包括程序代码模块，当在控制如权利要求I所述的图像处理装置的计算机上运行该计算机程序时，该程序代码模块使该图像处理装置执行如权利要求13所述的图像处理方法的步骤。应当理解的是，权利要求I的图像处理装置、权利要求13的图像处理方法和权利要求14的图像处理计算机程序具有类似的和/或等同的优选实施例，该优选实施例尤其是从属权利要求中定义的。应当理解的是，本发明的优选实施例也可以是从属权利要求与相应的独立权利要求的任意组合。本发明的这些和其他方面将根据下文描述的实施例而变得清楚，并将参考下文描述的实施例而得到阐述。
在下面的图中图I示意性并且示范性地示出了用于处理图像的图像处理装置的实施例，图2示意性并且示范性地示出了图像提供单元的实施例，图3示意性并且示范性地示出了多色辐射源的发射谱，图4示意性并且示范性地示出了光电谱、康普顿谱和具有k边缘的谱，图5示意性和示范性地示出了微分相位衬度计算机断层摄影系统的实施例的要素，以及图6示意性地示出了说明用于处理图像的图像处理方法的实施例的流程图。
具体实施例方式图I示意性并且示范性地示出了用于处理图像的图像处理装置I的实施例。图像·处理装置I包括用于提供对象的第一图像以及同一对象的第二图像的图像提供单元2。图像处理装置I还包括用于对第一图像进行滤波的滤波单元3，其中该滤波单元3适于根据第一图像的图像值之间的第一相似度和第二图像的图像值之间的第二相似度对第一图像进行滤波。滤波单元3优先适于还根据第一图像的图像值之间的第一相似度以及第二图像的图像值之间的第二相似度对第二图像进行滤波。经滤波的图像能够在显示器4上显示出来。第一图像和第二图像可以是显示同一对象的二维图像、三维图像或四维图像，该对象可以是整个对象或仅对象的一部分，例如，仅对象的切片。如果对象是如人的心脏那样的运动对象，那么第一图像和第二图像显示优选处于同一运动状态的对象。在该实施例中，基于投影数据，特别是基于计算机断层摄影投影数据，来重建第一图像和第二图像，其中，根据同一投影数据来重建第一图像和第二图像。图像提供单元2适于提供第一图像和第二图像，使得它们显示对象的不同特性，特别是显示对象的不同物理特性。在实施例中，图像提供单元2适于提供谱计算机断层摄影图像，其中第一图像指示对象的第一物理特性，而第二图像指示对象的第二物理特性，其中第一物理特性和第二物理特性彼此不同。例如，第一图像和第二图像之一可以是康普顿图像，而第一图像和第二图像中的另一个可以是光电图像，或者第一图像和第二图像之一可以是康普顿图像或光电图像，而第一图像和第二图像中的另一个可以是康普顿图像和光电图像的组合。图像提供单元2可以是存储了第一图像和第二图像的存储单元，或者图像提供单元2也可以是允许通过数据连接接收第一图像和/或第二图像的图像接收单元。图像提供单元2也可以是生成第一图像和第二图像的成像模态(如谱计算机断层摄影系统)。在下文中将参考图2示意性并且示范性地描述这样的谱计算机断层摄影系统。谱计算机断层摄影系统包括能够围绕旋转轴R旋转的台架(gantry) 14,旋转轴R平行于z方向延伸。将诸如X射线管之类的辐射源15安装到台架14上。在高实施例中，辐射源15发射多色辐射。辐射源15设置有准直器设备16，准直器设备16根据辐射源15发射的辐射形成锥形辐射束17。在其他实施例中，准直器设备16可以适于形成具有另一形状，例如具有扇形的辐射束。辐射在圆柱检查区5中的兴趣区域中横贯对象(图2中未示出)，该对象例如是患者或技术对象。在横贯兴趣区域之后，辐射束17入射到探测设备6上，在这一实施例中，探测设备6具有二维探测表面，其中探测设备6安装在台架14上。在另一实施例中，探测设备6可以包括一维探测表面。探测设备6优先包括能量分辨探测器元件。能量分辨探测器元件优先是光子计数探测器元件，该光子计数探测器元件例如按照对入射光子进行计数的原理工作，并且输出表示不同能量窗内的光子的数量的信号。例如，在Llopart X.等人的“First testmeasurements of a 64k pixel readout chip working in a single photon countingmode”，Nucl. Inst, and Meth. A, 509 (1-3) : 157-163，2003 中，以及在 Llopart, X.等人的“Medipix2:A 64_k pixel readout chip with 55 μ m square elements working in asingle photon counting mode”，IEEE Trans. Nucl. Sci. 49 (5) : 2279-2283，2002 中描述了这样的能量分辨探测设备。优选地，调整能量分辨探测器元件，使得每一探测器元件提供对应至少两个能量窗的至少两个能量分辨探测信号。然而，具有甚至更高的能量分辨率，以便增强成像系统的灵敏度和噪声鲁棒性是有利的。 由电机7以优选恒定但是可调节的角速度驱动台架14。提供另外的电机8，来用于使对象(例如安置在检查区5内的患者台上的患者)平行于旋转轴或者z轴的方向发生位移。例如，这些电机7、8由控制单元9控制，使得辐射源15和检查区5，特别是兴趣区域沿螺旋轨迹相对于彼此移动。对象或者检查区5 (特别是兴趣区域)不发生移动，而是辐射源15旋转，即辐射源15相对于兴趣区域沿圆形轨迹行进，也是可以的。在该实施例中，由探测设备6获取的数据是取决于能量的探测数据，将该探测数据提供给重建单元10，来用于根据所提供的取决于能量的探测数据重建兴趣区域的至少第一图像和第二图像。重建单元10也优先由控制单元9控制。将重建的图像提供给滤波单元3，以用于对至少第一图像进行滤波。重建单元10包括用于根据探测数据来确定至少两个衰减分量的计算单元12。探测数据的衰减分量是仅由一种物理效应或者由几种物理效应引起的，例如，一个衰减分量可以是由康普顿效应引起的分量，而另一衰减分量可以是由光电效应引起的分量。另外的衰减分量可以是由兴趣区域内存在的k边缘引起的分量。可替代地或附加地，衰减分量可以是由例如兴趣区域内的某种材料的吸收引起的分量。例如，一衰减分量可以是由某种材料的吸收引起的分量，而另一衰减分量可以是由另一种材料的吸收引起的分量。重建单元10还包括反投影单元13，反投影单元13用于对计算出的探测数据的衰减分量进行反投影，从而重建出对应的衰减分量图像。例如，可以将康普顿图像重建成第一图像和第二图像之一，并且可以将光电图像重建成第一图像和第二图像中的另一个。到达计算单元12的输入是多个能量仓(在该实施例中最少是三个能量仓)的能量分辨探测数据屯。这些探测数据Cli示出第i个能量仓匕的谱灵敏度Di (E)。此外，多色辐射源15的发射谱T(E)通常是已知的或者能够测量的。在图3中示意性并且示范性地示出了多色辐射源的这样的发射谱T(E)的示例。在计算单元12中，将探测数据Cli的生成模型化为具有谱P (E)的光电效应、具有谱C(E)的康普顿效应和具有谱K (E)的k边缘的线性组
口 ο在图4中示范性并且示意性地示出了谱P(E)、C(E)和K(E)。可以通过下述方程系来将探测数据的生成模型化
Cli= f ClET(E)Di (E) exp [_( P 光子 P (E)+ P 康普顿 C (E) + P k 边缘 K (E))]⑴其中，P光子、P *普顿和P k边缘分别是光电分量、康普顿分量和k边缘分量的密度长度乘积(density length product)。由于对于至少三个能量仓匕，b2, b3可获得至少三个探测信号Cl1, d2, d3，因而形成了至少三个方程的方程系，该方程系具有三个未知数，该三个未知数是三个密度长度乘积，因而能够在计算单元12中使用已知数值法解出该三个未知数。如果可获得多于三个的能量仓，那么优选的是，使用考虑了测量结果的噪声统计的最大似然法。通常，三个能量仓足够了。然而，为了提高灵敏度和噪声鲁棒性，优选的是，具有更多能量仓的更多探测数据。对于每一衰减分量而言，即对于每一密度长度乘积而言，可以通过对相应的密度长度乘积进行反投影来重建对象的衰减分量图像。可替代地或附加地，可以将不同密度长度乘积的图像组合，或者在反投影之前将不同密度长度乘积组合，以生成组合图像，该组合图像是由组合的衰减分量生成的对象图像。在另一实施例中，图像提供单元2能够适于提供微分相位衬度计算机断层摄影图像，其中第一图像和第二图像之一指示对象的吸收，而第一图像和第二图像中的另一个指·示由对象诱发的横贯对象的辐射的相移。特别是，第一图像优先显示对象的总衰减系数，而第二图像优先显示对象的折射率的实部。而且，可以将这些第一图像和第二图像存储在可以为图像提供单元2的存储单元内，或者可以通过也可以为图像提供单元2的图像接收单元接收这些第一图像和第二图像。然而，图像提供单元2也可以是微分相位衬度计算机断层摄影系统，其生成第一图像和第二图像，并且在下文中将参考图5对其进行示范性和示意性地描述。图5示出了一种微分相位衬度计算机断层摄影系统的示范性实施例的三维表示。与源光栅34相邻布置相当大的X射线源32。既然由于X射线源32与所发射的辐射的波长相对的尺寸所致，可以认为X射线源32是不相干的，那么应用源光栅34来用于提供多个单一的相干X射线源。X辐射35从X射线源32沿光轴37的方向放射，其有可能构成x射线的扇束或锥束。在图5中未描绘X射线射束的相应形状。X辐射35抵达对象36，横贯对象36，并且随后抵达分束器光栅38。分束器光栅38的沟槽或缝隙相对于分束器光栅的实区(阻挡区域)改变了经过的电磁辐射的相位。相应地，实现了相移Φ，特别是，实现了相移η。将分析器光栅(analyzer grating) 40布置在分束器光栅38和x射线探测器42之间。源自分束器光栅38朝X射线探测器的方向的多个波抵达分析器光栅40，并且随后在X射线探测器42的表面上产生强度调制图案。通过将分束器光栅38与分析器光栅40相对移动，从而使光栅相对于彼此发生位移，尤其是分束器光栅38和分析器光栅40的光栅周期的一部分，由此通过图像探测器42可获得多个透射图像。因而，在透射图像的采集期间，执行所谓的相位步进，其中，在一个周期内例如以四到十步改变光栅38、40的相对位置，并且针对每一相对位置测量透射图像。根据该系列的透射图像，通过针对每一探测器像素对所测量的强度进行简单地求均值来获得常规的吸收图像。与常规的X射线透射成像中一样，这一测量涉及根据比尔定律沿X射线路径在总的衰减系数上的线性积分。通过在相位步进循环上对测量的强度的动态特性进行分析而得到有关折射率的实部的信息。该强度优先随分析器光栅40的光栅周期而周期性变化。典型地，可以假定正弦变化。正弦的相位与碰撞分束器光栅38的相位波前的梯度有关。根据物理学，该相位波前的梯度与对象的折射率的实部的投影的一阶导数线性相关。因而，能够基于正弦的相位来确定对象的折射率的实部的投影。通过对所确定的对象折射率的实部的投影进行反投影，能够重建出显示对象的折射率实部的图像。因此，微分相位衬度计算机断层摄影系统能够重建显示对象的总衰减系数的第一图像以及显示对象的折射率实部的第二图像。例如在Timm Weitkamp等人的文章“X-ray phase imaging with agrating interferometer”，OPTICS EXPRESS,第 13 卷，第 16 期，第 6296-6304 页(2005年8月)中公开了对已知微分相位衬度计算机断层摄影技术的更详细描述，通过引用将该文章并入本文。再次参考图1，滤波单元3优先适于执行滤波，使得通过使用下列的滤波方程来对第一图像进行降噪
Σ Li (xi)-u (xj))L2w ixi)~u (xj))K(xi-xj)
Xj& 其中，表示在第一图像的位置足处的第一像元的平均图像值，Ai表示包括位置足处的第一像元的第一区域，A表示处于第一区域Ai内的第一区域像元的位置，
表示a)第一区域像元 的图像值和可选的围绕位置七处的相应第一区
域元素的像元的图像值与b)位置足处的第一像元的图像值和可选的围绕第一图像的第一像元的像元的图像值之间的第一相似度。函数/^7(2)汉)-#( ))表示a)第二图像中与第一区域对应的第二区域内的i/立
置处的第二区域像元的图像值和可选的围绕位置士处的第二区域元素的像元的图像值与b)第二图像中的于处于第一图像内的相同位置上的第一像元对应的位置名处的第二像元的图像值和可选的围绕第二像元的像元的图像值之间的第二相似度。函数Ι( , )是取决于位置足和七处的像元之间的距离的距离函数，而表示第一图像内的位置士处的未经滤波的图像值。可以将第一相似度L1、第二相似度L2和距离函数K的组合除以方程(2)中所示的分母看作是根据方程(2)执行加权平均的权重。因而，根据方程(2)，滤波单元3优先适于对第一图像的像元的图像值进行加权平均，其中，为了确定第一图像的第一像元的平均图像值，对第一区域(包括第一像元)内的第一区域像元的图像值进行加权并求平均，其中利用根据第一相似度以及根据第二相似度的权重对第一区域像元的图像值进行加权，该第一相似度是a)第一区域像元的图像值和围绕相应的第一像元的可选的像元的图像值与b)第一像元的图像值和围绕第一图像的第一像元的像元的图像值之间的相似度，该第二相似度是a)第二图像中的与第一区域对应的第二区域内的第二区域像元的图像值和可选的围绕相应的第二区域像元的像元的图像值与b)第二图像中的与第一像元对应的第二像元的图像值和可选的围绕第二像元的像元的图像值之间的相似度。可以将矢量 Hi.,,)昏作是以无,/为中心的图像块(image patch)。例如，在二维情况下，如果考虑可选的围绕相应的第一像元的像元的图像值和可选的围绕相应的第二像元的像元的图像值，那么该图像块可以是以位置处的像素为中心的9X9子图像。该块也可以具有其他尺寸。例如，该补片可以具有5X5或者7X7的尺寸。此外，也可以采用三维块。如果不考虑可选的图像值，即，如果只考虑位置名和七处的单一图像值之间的差，那么相应的矢量将分别只包括位置？或士处的相应的单个图像值。那么该块尺寸为1X1。如果在两幅以上的图像的基础上执行滤波，那么滤波单元3优先适于根据下列滤波方程对第一图像进行滤波
权利要求
1.一种用于处理图像的图像处理装置，所述图像处理装置(I)包括 -图像提供单元(2)，其用于提供对象的第一图像以及同一对象的第二图像， -用于对所述第一图像进行滤波的滤波单元(3)，其中，所述滤波单元(3)适于根据所述第一图像的图像值之间的第一相似度以及根据所述第二图像的图像值之间的第二相似度对所述第一图像进行滤波。
2.根据权利要求I所述的图像处理装置，其中，所述滤波单元(3)适于对所述第一图像的像元的图像值进行加权平均，其中，为了确定所述第一图像的第一像元的平均图像值，对包括所述第一像元的第一区域内的第一区域像元的图像值进行加权并且求平均，其中，采用取决于下述内容的权重对所述第一区域像元的图像值进行加权 -作为a)所述第一区域像元的图像值与b)所述第一图像的所述第一像元的图像值之间的相似度的第一相似度； -作为a)所述第二图像中的与所述第一区域对应的第二区域内的第二区域像元的图像值与b)所述第二图像中的与所述第一像元对应的第二像元的图像值之间的相似度的第二相似度。
3.根据权利要求2所述的图像处理装置，其中，调整所述滤波单元(3)，使得如果a)所述第一区域像元与所述第一像元之间的所述第一相似度和b)所述第二区域像元与所述第二像元之间的所述第二相似度中的一个降低，以及如果a)所述第一区域像元与所述第一像元之间的所述第一相似度和b)所述第二区域像元与所述第二像元之间的所述第二相似度中的另一个保持不变或者降低，则第一区域像元的图像值的权重降低。
4.根据权利要求2所述的图像处理装置，其中，所述第一相似度是a)所述第一区域像元的图像值和围绕相应的第一区域像元的像元的图像值与b)所述第一像元的图像值和围绕所述第一图像的所述第一像元的像元的图像值之间的相似度，并且其中，所述第二相似度是a)所述第二图像中的与所述第一区域对应的第二区域内的第二区域像元的图像值和围绕相应的第二区域像元的像元的图像值与b)所述第二图像中的与所述第一像元对应的第二像元的图像值和围绕所述第二像元的像元的图像值之间的相似度。
5.根据权利要求I所述的图像处理装置，其中，如果所述第一图像将被降噪，则所述第一相似度是所述第一图像的图像值相似的概率，并且其中，如果所述第二图像将被降噪，则所述第二相似度是所述第二图像的图像值相似的概率。
6.根据权利要求5所述的图像处理装置，其中，所述滤波单元(3)适于基于图像值之间的差异来确定所述图像值相似的概率。
7.根据权利要求I所述的图像处理装置，其中，所述滤波单元(3)适于根据所述第一图像和所述第二图像中的至少一个中的像元之间的距离来对所述第一图像进行滤波。
8.根据权利要求I所述的图像处理装置，其中，所述图像提供单元(2)适于提供所述第一图像和所述第二图像，使得它们显示所述对象的不同特性。
9.根据权利要求8所述的图像处理装置，其中，所述图像提供单元(2)适于提供谱计算机断层摄影图像，其中，所述第一图像指示所述对象的第一物理特性，而所述第二图像指示所述对象的第二物理特性，所述第一物理特性和所述第二物理特性互不相同。
10.根据权利要求9所述的图像处理装置，其中，所述第一图像和所述第二图像中的一个是康普顿图像，而所述第一图像和所述第二图像中的另一个是光电图像。
11.根据权利要求9所述的图像处理装置，其中，所述第一图像和所述第二图像中的一个是康普顿图像或光电图像，而所述第一图像和所述第二图像中的另一个是康普顿图像和光电图像的组合。
12.根据权利要求I所述的图像处理装置，其中，所述图像提供单元(2)适于提供微分相位衬度计算机断层摄影图像，其中，所述第一图像和所述第二图像中的一个指示所述对象的吸收，而所述第一图像和所述第二图像中的另一个指示由所述对象诱发的贯穿所述对象的辐射的相移。
13.一种用于处理图像的图像处理方法，所述图像处理方法包括 -提供对象的第一图像和同一对象的第二图像， -根据所述第一图像的图像值之间的第一相似度以及根据所述第二图像的图像值之间的第二相似度对所述第一图像进行滤波。
14.一种用于处理图像的图像处理计算机程序，所述计算机程序包括程序代码模块，当在控制根据权利要求I所述的图像处理装置的计算机上运行所述计算机程序时，所述程序代码模块使所述图像处理装置执行根据权利要求13所述的图像处理方法的步骤。
全文摘要
本发明涉及一种图像处理装置(1)，其包括用于提供对象的第一图像和同一对象的第二图像的图像提供单元(2)，还包括根据第一图像的图像值之间的第一相似度以及根据第二图像的图像值之间的第二相似度对第一图像进行滤波的滤波单元(3)。即使由于噪声的原因对第一图像和第二图像之一的图像值发生了干扰，这仍然允许根据所述图像值属于对象的同一部分(例如，如果对象是人或者动物的话，属于对象的组织或骨骼物质)的可能性对第一图像进行滤波，从而提高了例如滤波过程的边缘保护特性。
文档编号G06T5/00GK102906783SQ201180025245
公开日2013年1月30日 申请日期2011年5月16日 优先权日2010年5月21日
发明者T·克勒, B·J·布伦德尔 申请人:皇家飞利浦电子股份有限公司