CN101930595A

CN101930595A - 图像处理方法和图像处理设备

Info

Publication number: CN101930595A
Application number: CN2010102123012A
Authority: CN
Inventors: 近江裕行
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2009-06-26
Filing date: 2010-06-28
Publication date: 2010-12-29
Also published as: KR20110000537A; EP2267655A2; JP2011005050A; KR101264182B1; US20100329533A1; EP2267655A3; CN103544684A

Abstract

本发明涉及一种图像处理方法和图像处理设备。该图像处理设备用于对由X射线摄像系统拍摄的图像进行色调转换，该图像处理设备提取用作图像基准的基准区域，并且如果存在图像变化，则进行色调转换，从而在所提取的基准区域中抑制由图像变化所引起的对比度变化，并且在其它区域中反映该对比度变化。

Description

图像处理方法和图像处理设备

技术领域

本发明涉及一种用于将X射线图像变换为具有最佳色调的输出图像的X射线图像处理方法。更具体地，本发明涉及用于创建色调转换曲线以限定输出图像的对比度的图像处理方法(即，色调映射(tone-mapping)方法)和图像处理设备，其中，所创建的色调转换曲线至少部分包含前一图像的色调转换曲线。

背景技术

对于X射线图像，已经提出了各种类型的色调转换，以提高输入X射线图像的输出版本，从而力图提高医生的诊断能力。进行用于最优化X射线图像的观察区域的色调转换，以提高医生的诊断能力。为了进行该最优化，已经提出了从X射线图像中提取相对于噪声和运动具有更高鲁棒性的观察区域。已经使用S形函数作为色调曲线以最优化该区域的输出图像。当拍摄X射线图像作为(由帧序列组成的)运动图像时，还需要考虑这些帧之间的对比度变化并使其稳定。由于X射线变化以及被摄体变化，因而引起帧之间的对比度变化。X射线变化是指即使在恒定摄像条件下也产生的X射线的量的变化、以及由于X射线操作或控制引起的摄像条件的变化。例如，被摄体变化可以是指肺野由于呼吸而移入和移出摄像区域，或者是指注入造影剂。

迄今为止，已通过使X射线更稳定来应对由X射线的量的差异所引起的帧之间的对比度变化。另一方面，已通过分析被摄体的方法或通过控制色调转换曲线来应对被摄体变化。

分析被摄体的示例方法包括涉及根据输入图像的像素值创建直方图、提取最小值和最大值并在时间方向上对这些值进行滤波的方法(例如，参见日本专利3334321)。在时间方向上对最小值和最大值进行滤波使得能够抑制对对比度变化的灵敏度并使其稳定。

用于控制色调转换的示例方法包括涉及通过分析输入图像检测场景变化、并基于场景变化所需的时间将新创建的色调转换曲线与过去的色调转换曲线进行合成的方法(例如，参见日本专利4050547)。根据该方法，可以通过在时间方向上对色调转换曲线进行滤波来抑制对对比度变化的灵敏度并使其稳定。

然而，上述传统技术中存在有以下问题。即使通过对直方图的最小值和最大值(表示特别暗和特别亮的区域的最小值和最大值)进行滤波来抑制图像整体的对比度变化，也不能抑制正被注视观察的区域中的变化，并且(诸如正在使用输出图像进行诊断的医生等的)用户敏锐地感受到该状况。响应于此而抑制观察区域中的变化的尝试导致了饱和和限幅，这是因为不能反映运动图像中的对比度变化并且不能创建最适于输入帧的输出图像。

期望一种在反映图像整体的对比度变化的同时抑制正被注视观察的区域中的对比度变化的方法。此外，存在如下多种情况：在图像序列内基于场景变化所需的时间来控制色调转换曲线的情况下，期望与该场景变化所需的时间无关地控制对比度变化。

发明内容

期望提供在通过抑制运动图像的闪烁获得稳定色调的同时针对(一系列帧或构成运动图像的图像内的)各运动图像帧获得最佳色调的设备和方法。还期望提供在创建没有感觉不自然的观察区域的图像的同时反映(即，考虑)每一帧的对比度变化的设备和方法。期望提供在具有不同亮度或对比度等级的一系列时序帧或同步帧中、或者甚至在具有使观察者不能辨别具有较低对比度等级的其它低亮度区域的高亮度区域的单个图像中获得最佳色调的设备和方法。

根据本发明的一方面，提供了一种图像处理方法，用于对由X射线摄像系统所生成的图像序列的图像进行色调转换，所述图像处理方法包括：提取步骤，用于从所述图像序列的第二图像中提取在所述图像序列的图像中用作基准区域的区域；以及色调转换步骤，用于在相对于所述图像序列中的第一图像在所述第二图像中存在对比度变化的情况下，对所述第二图像进行色调转换，从而在所述基准区域中抑制在所述第一图像和所述第二图像之间产生的对比度变化，并且在所述基准区域以外的区域中反映所述对比度变化。

根据本发明的另一方面，提供了一种图像处理设备，用于对由X射线摄像系统所生成的图像序列的图像进行色调转换，所述图像处理设备包括：提取部件，用于从所述图像序列的第二图像中提取在所述图像序列的图像中用作基准区域的区域；以及色调转换部件，用于在相对于所述图像序列中的第一图像在所述第二图像中存在对比度变化的情况下，对所述第二图像进行色调转换，从而在所述基准区域中抑制所述对比度变化，并且在所述基准区域以外的区域中反映所述对比度变化。

通过以下(参考附图)对典型实施例的说明，本发明的其它特征将变得明显。

附图说明

图1是示出根据本发明的图像处理设备的硬件结构的示例的框图。

图2是示出图1所示的图像处理单元104的详细结构的框图。

图3A和图3B示出两个X射线图像帧中的像素值的直方图。

图4A和图4B示出合成色调转换曲线的特性。

图5是示出根据本发明的反馈系数的图。

图6是示出第一实施例中的图像处理的流程图。

图7A和图7B是示出基准区域中的合成色调转换曲线和反馈系数的图。

图8是示出基于反馈系数对色调转换曲线进行合成的处理的示意图。

图9是示出第二实施例中的图像处理的流程图。

具体实施方式

下文中，将参考附图详细说明用于实现本发明的实施例。在实施例中，将说明用于创建如下图像的方法，在该图像中，通过将正被观察的区域的反馈系数设置为大值而在该正被观察的区域中抑制对对比度变化的灵敏度，并且通过将剩余区域的反馈系数设置为小值而在该剩余区域中反映(即，考虑)对比度变化。“反映”的含义是指在非基准区域中以某种方式识别到对比度变化。例如，如果显示处理完成，则显示对比度变化，因此在该显示中反映已出现的对比度变化。另一方面，如果没有进行完整的显示处理并且仅进行没有引起显示的内部处理，则计算对比度变化，并且在合成色调转换曲线的处理中考虑该对比度变化。

首先，将使用图1说明根据本发明的图像处理设备的硬件结构的示例。图1所示的示例是在个人计算机(PC)上实现图像处理设备的情况下的结构。如图1所示，图像处理设备100包括：中央处理单元，即CPU 101；只读存储器，即ROM 102；随机存取存储器，即RAM 103；图像处理单元104；硬盘驱动器，即HDD 105；输入/输出接口，即输入/输出I/F 106；以及网络接口，即网络I/F 107。组件101～107经由系统总线108相连接。CPU101根据存储在ROM 102和HDD 105等中的计算机程序来控制整体设备。ROM 102是存储启动程序和控制数据等的存储器。RAM 103是当CPU 101执行处理时对程序进行展开的、且定义有各种表和工作区等的存储器。

图像处理单元104对输入X射线图像进行(稍后详述的)色调转换等的图像处理。在图1中，将图像处理单元104实现为专用图像处理板，但是还可以实现为软件模块。换句话说，可以根据目的来适当地实现图像处理单元104。HDD 105存储操作系统(OS)和应用程序软件等。输入/输出I/F 106是与显示器等的输出设备和键盘或鼠标等的输入设备的接口。网络I/F 107是与局域网(LAN)等的外部网络的接口。

图像处理单元104连接至X射线摄像系统的网络140。该网络140可以包括控制局域网(control area network，CAN)，或者可以包括光纤。X射线生成设备110、医疗监视器120和X射线传感器(平面检测器)130连接至网络140。此外，影像存档和通信系统(picture archiving and communication system，PACS)和用于存储X射线图像的医学影像设备内置(intra-modality)硬盘设备也可以连接至网络140。可以通过将命令从图像处理单元104发出至X射线摄像系统中的X射线生成设备110或发出至X射线摄像系统中的X射线传感器130来控制被摄体的摄像。

接着，将使用图2来说明图1所示的图像处理单元104的详细结构。图像处理单元104包括图像输入单元201、色调转换曲线计算单元202、基准区域提取单元203、色调转换曲线合成单元204、色调转换曲线存储单元205、色调转换单元206、以及图像输出单元207。图像输入单元201输入要处理的X射线图像，并且进行(稍后所述的)色调转换的处理之前所需的处理。这里，所需的处理例如包括校正X射线传感器特性或校正系统特性。图像输入单元201还根据需要进行用以抑制随机噪声的处理和图像增强。

色调转换曲线计算单元202计算用于对由图像输入单元201处理后的X射线图像进行色调转换的色调转换曲线。基准区域提取单元203提取色调转换曲线合成单元204在对色调转换曲线进行合成时所参考的“要被注视观察的区域”。色调转换曲线合成单元204对前一帧的色调转换曲线和当前帧的色调转换曲线进行合成。色调转换曲线存储单元205保存由色调转换曲线合成单元204合成后的色调转换曲线。该色调转换曲线存储单元205可以位于RAM 103中。

色调转换单元206使用由色调转换曲线合成单元204合成后的色调转换曲线，对输入X射线图像进行色调转换。图像输出单元207对处理后的图像进行所需的处理，并且将该图像输出到医疗监视器120或硬盘设备等。这里，“所需的处理”例如可以包括监视器伽玛转换或几何转换等。

将使用图3A和图3B以及图4A和图4B来说明上述结构中的图像处理单元104的具体处理。图3A是X射线图像序列的第N-1个图像或帧中的像素值的直方图。“序列”的含义是指按时间顺序拍摄的帧序列或者同时拍摄的多个帧。另外可选地，帧序列可以是从单个X射线图像提取出的多个帧，其中该多个帧对于图像中的像素具有不同的亮度或强度等级。例如，如果正被X射线照射的被摄体在曝光至X射线期间发生移动，则可以有效地从图像中提取具有不同强度等级的帧序列。

图3A的点划线表示对于该直方图的最佳或线性色调转换曲线。这里，“最佳”是指直方图范围分布在整个输出范围上。如稍后所述，输出范围是构成作为图像输出单元207的输出的图像的像素值的范围，并且可以缺少输入像素的最大值和最小值。

图3B是X射线图像序列的第N个帧中的像素值的直方图。从图3A和图3B可以看出，由于被摄体变化或造影剂的注入等引起的可应用输入像素值的变化，因而相对于第N-1个帧，沿着x轴(输入像素值)的直方图的范围在第N个帧中改变。因为输入像素值的范围依赖于X射线传感器130所接收到的光的亮度，因此与如上所述的X射线变化相比，输入像素值的范围的变化更可能是由被摄体变化引起的。该直方图的最佳或线性色调转换曲线如由双点划线所示。注意，在图3A和图3B中，附图标记301表示正被注视观察的内脏等的区域。

这里，通过色调转换曲线的倾斜度(gradient)来确定各帧的对比度。当从第N-1个帧改变为第N个帧时，作为考虑了正被X射线照射的被摄体的亮度变化的结果，色调转换曲线的倾斜度降低，并且对比度减小。通过“考虑”第N-1个帧和第N个帧之间的被摄体的变化，理解了一些替代方案。比较图3A和图3B，第N-1个图像的最大和最小输入像素值(即，具有极少出现的低值和高值的输入像素值)不存在。可能因为缺少造影剂等的高对比度被摄体或者因为X射线传感器130的设置，原本就不存在这些值，或者可能在图像处理期间将(亮度或强度)像素值钳位在特定阈值之间。例如，可以对该阈值进行设置，以去除由金属移植引起的特别亮的片或暗片。然而，正是由于拍摄了图像并对其进行处理，因而本实施例的结果是对于第一图像(N-1)，输出像素值相对于输入像素值的倾斜度与对于第二图像(N)的倾斜度相比更陡，其中第二图像没有考虑到所有输入像素值(即，甚至没有考虑具有较高值和较低值的像素)。

由此，在图4A所示的传统技术中，通过对第N-1和第N个帧的色调转换曲线取平均，来创建合成色调转换曲线(实线)。实际进行该取平均的方式是将图3B的直方图中的最大和最小像素值反馈至图像输入单元201，并且针对第N个图像的输出考虑这些值而进行取平均。

然而，即使利用由实线示出的该合成色调转换曲线，正被注视观察的区域401中的图像闪烁也是明显的。

由此，如图4B所示，通过对合成色调转换曲线进行加权以使其更近似接近于正被注视观察的区域401中的第N-1个帧的色调转换曲线，在该区域401中维持第N-1个帧的对比度。因此，对于观察者来说，该区域中的较不明显的对比度更加可见。此外，通过对合成色调转换曲线进行加权以使得随着相对注视观察的区域401的距离增加而接近第N个帧的色调转换曲线，在该区域401之外的区域中显示第N个帧的对比度。这样，在注视观察的区域之外的区域中可以看见诸如由造影剂所引起的对比度的较大的极值。

可以通过以下步骤来实现以上操作：获得第N-1个帧中的多个像素的对比度值；获得第N个帧中的多个像素的对比度值；以及有效地生成如下的第三帧，其中，该第三帧在其基准区域中包含第N-1个帧的多个像素的对比度值，并且在该基准区域以外的区域中包含第N个帧的多个像素的对比度值。然而，为了获得基准区域和基准区域之外的区域的对比度值之间的平滑转变，优选在第三帧中具有如下色调转换曲线，其中，该色调转换曲线没有精确地与基准区域中的第N-1个帧的色调转换曲线相同，而是接近该基准区域中的第N-1个帧的色调转换曲线；并且没有精确地与基准区域之外的第N个帧的色调转换曲线相同，而是接近该基准区域之外的第N个帧的色调转换曲线，在两个对比度值倾斜度之间逐渐弯曲。这可以通过将第N-1和第N个帧的色调转换曲线的(如图4A中的实线所示的)平均值乘以用于对第三帧的色调转换曲线进行期望调整的第三曲线(图5中的点线)来进行。该第三曲线已知为反馈系数α。

更具体地，如图5所示，对反馈系数α进行设置。这里，对反馈系数α进行设置，从而当将该反馈系数α乘以第N-1和第N个帧的色调转换曲线的平均值时，如图4B所示，α的值越大，由此得到的合成色调转换曲线接近第N-1个帧的色调转换曲线，且α的值越小，由此得到的合成色调转换曲线接近第N个帧的色调转换曲线。稍后将进一步说明该反馈系数α。因此，可以抑制正被注视观察的区域401中的图像闪烁，并且反映图像整体的对比度变化。

通过色调转换曲线计算单元202来选择反馈系数α，以使得由此得到期望的色调曲线。优选地该反馈系数α在基准区域(正被注视观察的区域)处具有最大值，并且在该区域之外具有最小值。以下将详细进行说明。

第一实施例

将使用图6来说明第一实施例中用以从X射线摄像系统获取X射线图像并对该X射线图像进行色调转换的图像处理。首先，通过图像输入单元201从X射线摄像系统输入要经过色调转换的X射线图像(S601)。接着，进行考虑到X射线传感器130的特性和X射线摄像系统的特性的校正，作为预处理(S602)。对X射线传感器130的特性的校正可以包含进行偏移校正或缺陷校正等。对X射线摄像系统的特性的校正可以包含进行调制传递函数(modulation transfer function，MTF)改善或网格线(grid line)校正等。另外，除了校正X射线传感器130的特性和系统特性以外，还可以根据需要进行用于抑制随机噪声或系统噪声的噪声抑制处理以及用于增强边缘的增强处理等。

这里，预处理后的X射线图像是原始图像。然后，进行场景变化检测(S603)。这里，场景变化是指正被X射线照射的被摄体变化或正被注视观察的观察区域在帧之间变化。在图像的辉度由于X射线操作等导致不稳定的情况下，也检测到场景变化。关于检测方法，如果整个图像的平均辉度超过预定阈值，或者如果X射线管电压或管电流的变化超过预定阈值，则检测到场景变化。这里，如果存在场景变化，则处理直接进入S607。另一方面，如果不存在场景变化，则处理进入S604，并且通过基准区域提取单元203从原始图像提取被摄体区域。

在S604中，首先，从原始图像中检测照射野(treatment field)以外的或不存在被摄体的区域，并且将剩余区域识别为被摄体区域。用于识别照射野的方法包括涉及得到轮廓并计算微分值的方法以及使用神经网络的方法。另一方面，用于检测不存在被摄体的区域的方法可以涉及创建像素值的直方图并基于像素的辉度值进行检测。因此，可以使用这些方法提取被摄体区域。除进行照射野以外的以及不存在被摄体的区域的检测之外，在根据需要去除由于被摄体中的移植金属等而引起的图像中的假象(artefact)之后，可以进行被摄体区域的识别。可以通过示出移植金属或其它反射/高密度材料的区域中的像素的高辉度值来确定这类假象。因此，可以提取直方图中的非常亮的像素值，以去除这些种类的图像假象。因此，基于像素辉度的提取处理可以产生图3A所示的直方图形状。

接着，基于所提取的被摄体区域来提取基准区域(S605)。这里，基准区域是要被注视观察的区域301、401。可以使用诸如内脏的图像形式的表现等的解剖学要素来指定该基准区域。在第一实施例中，使用摄像区域信息(即，与图像中的期望区域有关的信息)来指定解剖学要素。创建表示被摄体区域的像素值的直方图，并且基于摄像区域信息和直方图的形状来确定基准区域。例如，在对腹部区域进行摄像的情况下，可以将该区域大体地划分为肠、肠以外的脏器、骨头、以及其它区域。因此，将自动判别分析应用于直方图，以将该直方图划分为四个区域，并将上述的各解剖学结构分配至区域。将分配给在本示例中作为最受关注的区域的肠的直方图范围确定为基准区域。

注意，当提取基准区域时，除了上述定义的解剖学要素以外，还可以使用摄像技术信息(即，与摄像技术有关的信息)。在这种情况下，创建表示被摄体区域的像素值的直方图，并且基于摄像技术信息和直方图的形状来确定该基准区域。例如，在进行肾血管造影(renal angiography)的摄像技术的情况下，可以将该区域大体地划分为肾血管、肾脏、肾脏以外的脏器、以及其它区域。因此，将自动判别分析应用于该直方图，以将该直方图划分为五个区域，并将各解剖学结构分配至区域。因为该基准区域依赖于所使用的摄像技术(在这种情况下为肾血管造影术)，因此该基准区域被确定为与血管造影术有关的区域。因此，接着将分配给作为在血管造影术中最受关注的区域的肾血管和肾脏的直方图范围确定为基准区域。

可选地，当提取基准区域时，可以使用统计学要素。例如，创建直方图作为统计学要素，并且在该直方图内，可以将累积直方图的40％和60％的位置之间的区域确定为基准区域。可选地，可以将直方图范围自身的40％和60％的位置之间的区域确定为基准区域。

按照如下说明当提取基准区域时所使用的统计学要素的示例。可以使用包含被摄体区域的中心的预定ROI(region ofinterest，关注区域)作为统计学要素。例如，对包含被摄体区域的中心的N*N矩形ROI进行设置，并计算ROI内的像素值的直方图。将ROI内的直方图的累积直方图的40％和60％的位置之间的区域确定为基准区域。可选地，以上述ROI中的平均值作为基准区域的中心像素，并基于基准区域的中心像素，可以将规定像素范围确定为基准区域。

接着，针对所获得的基准区域计算反馈系数(S606)。如图5所示，该反馈系数可以是反馈系数在基准区域内达到其最大值的函数。具体地，可以通过诸如以下等式1等的三次函数来近似该反馈系数，其中，α_min是反馈系数的最小值，x是针对相应点处的反馈系数的当前像素值，x_max是原始图像中的最大像素值，并且x_basis是反馈系数在基准区域内达到其最大值处的原始图像的像素值。k是依赖于与基准区域的距离的加权系数。

x≤x_basis：α＝k₁x³+k₂x²+k₃x+α_min

x＞x_basis：α＝k₁(x_max-x)³+k₂(x_max-x)²+k₃(x_max-x)+α_min (1)

x_basis被确定为基准区域范围中的中间点或者基准区域范围中的累积直方图的50％的位置。等式1的函数可以用于如图7A所示的对比度变化，但是不能应用于如图7B所示的对比度变化。通过基于反馈系数的最小值α_min和最大值α_max利用样条插值、多项式插值、或可选的N维函数进行近似，来计算图7B所示的情况下的反馈系数的函数。根据本实施例，为了确定当前的色调转换曲线，如上所述使用前一色调转换曲线。为了确保在基准区域之外的图像的区域内可识别从前一图像到当前图像的对比度变化而在该基准区域内抑制该对比度变化，最大反馈系数值α_max期望为0.5以上。

接着，通过色调转换曲线计算单元202来计算色调转换曲线(S607)。这里，预先确定直线或S形函数等的用作色调转换曲线的基础的基本形状。计算色调转换曲线，从而将在S604计算出的被摄体区域分配至上述基本形状。

接着，色调转换曲线合成单元204基于在S606计算出的反馈系数，针对原始图像的每一个像素值，对所保存的一帧前的过去的色调转换曲线和在S607计算出的新的色调转换曲线进行合成(S608)，由此有效地创建包含应用于原始图像的每一个像素值的合成色调曲线的第三帧。

图8示出基于反馈系数对色调转换曲线进行合成的处理。将新创建的第N个帧的色调曲线乘以1-α，并且将第N-1个帧的色调曲线乘以α。将这两个乘积相加在一起以生成合成色调转换曲线。通过以下等式2来表示合成色调转换曲线Tc_merge，其中，Tc_new是新的(第N个帧的)色调转换曲线，Tc_old是过去的(第N-1个帧的)色调转换曲线，并且x是原始图像的像素值。注意，在第一个帧不存在过去的色调转换曲线的情况下，以α(x)＝0来计算新的色调转换曲线。如果在S 603检测到场景变化，则也以α(x)＝0来计算新的色调转换曲线。

Tc_merge(x)＝α(x)Tc_old(x)+(1-α(x))Tc_new(x) (2)

接着，将由色调转换曲线合成单元204合成后的色调转换曲线保存到色调转换曲线存储单元205(S609)。色调转换单元206使用合成色调转换曲线，对原始图像进行色调转换(S610)。这里，在输出图像之前根据需要进行后处理(S611)。注意，后处理可以包含比特转换、几何转换、或P值转换。当将图像输出至医疗监视器120时，还进行监视器伽玛转换等的处理。

最后，图像输出单元207将已经过了S610的色调转换和S611的后处理的图像输出至医疗监视器120、HDD 105或医学影像设备内置硬盘设备等(S612)。

根据第一实施例，可以使图像在正被注视观察的区域中稳定，并且可以创建反映整体图像或图像序列上的对比度变化的图像。实际上，例如，创建由像素强度或发光度的预定范围所限制的图像。使用该限制范围来示出关注区域之外的区域中的像素强度的大的范围(即，从非常暗到非常亮)，但是在该关注区域内，像素强度的较小的范围(排除强度的极值)在相同的所限制的预定范围上“展开”，以使对比度(即，辉度之间的差)更容易看见。例如，亮的假象在关注区域之外可视为亮的像素，但是在关注区域中，将看不见极值的辉度，并且观察者的眼睛将能够辨认出更精细的特征。结果，可以提高可视性，从而改善医生的诊断精度和手术精度。

第二实施例

接着，将参考附图来说明根据本发明的第二实施例。在第一实施例中，基于反馈系数对过去的色调转换曲线和新的色调转换曲线进行合成，然后保存合成色调转换曲线，但是在第二实施例中，与图6的步骤S608和S609相比，如由图9的步骤S908和S909的颠倒所示，在进行合成之前保存色调转换曲线。

第二实施例中的图像处理设备和X射线摄像系统的结构与图1和图2所示的第一实施例中的结构相同，并且将省略对其的说明。这里，将使用图9来说明第二实施例中的用以从X射线摄像系统中获取X射线图像并对该X摄像图像进行色调转换的图像处理。注意，图9所示的S901～S907和S910～S912的处理与图6所示的S601～S607和S610～S612的处理相同。因此，将说明S908和S909的处理。

色调转换曲线计算单元202将在S907计算出的新的色调转换曲线保存至色调转换曲线存储单元205(S908)。这里，所保存的新的色调转换曲线等同于在进行合成之前的色调转换曲线。接着，基于在S906计算出的反馈系数，对在S907计算出的新的色调转换曲线和已保存的过去的色调转换曲线进行合成(S909)。

这里，通过以下等式3来表示合成色调转换曲线Tc_merge，其中，Tc_new是新的色调转换曲线，Tc_oldmerge是过去的色调转换曲线的组合，并且x是原始图像的像素值。

Tc_merge(x)＝α(x)Tc_oldmerge(x)+(1-α(x))Tc_new(x)

Tc_oldmerge(n)＝kTc_old(n-1)+(1-k)Tc_old(n-2) (3)

Tc_old(n-1)是在S907计算出的第N-1个帧中的色调转换曲线。对过去在S907创建的色调转换曲线进行合成，并且对由此得到的(过去的)色调转换曲线进行合成，作为Tc_oldmerge(n)。

其它实施例

还可以通过系统或设备的计算机(或诸如CPU(中央处理单元)或MPU(微处理器单元)等的装置)或通过以下方法来实现本发明的各方面，其中该计算机读出并执行记录在存储器设备上的程序以进行上述实施例的功能，并且系统或设备的计算机通过例如读出并执行记录在存储器设备上的程序以进行上述实施例的功能来进行该方法的各步骤。为此，例如，经由网络或者从用作存储器设备的各种类型的记录介质(例如，计算机可读介质)向计算机提供程序。

尽管已经参考典型实施例说明了本发明，但是应该理解，本发明不局限于所公开的典型实施例，而是包括所有这类修改、等同结构和功能，并落入权利要求书的范围。

Claims

1.一种图像处理方法，用于对由X射线摄像系统所生成的图像序列的图像进行色调转换，所述图像处理方法包括：

提取步骤，用于从所述图像序列的第二图像中提取在所述图像序列的图像中用作基准区域的区域；以及

色调转换步骤，用于在相对于所述图像序列中的第一图像在所述第二图像中存在对比度变化的情况下，对所述第二图像进行色调转换，从而在所述基准区域中抑制在所述第一图像和所述第二图像之间产生的对比度变化，并且在所述基准区域以外的区域中反映所述对比度变化。

2.根据权利要求1所述的图像处理方法，其特征在于，还包括：

第一获得步骤，用于获得所述第一图像中的多个像素的对比度值；以及

第二获得步骤，用于获得所述第二图像中的多个像素的对比度值，

其中，所述色调转换步骤包括生成如下第三图像，其中，所述第三图像在所述基准区域中的对比度值接近所述第一图像中的多个像素的对比度值，并且所述第三图像在所述基准区域以外的区域中的对比度值接近所述第二图像中的多个像素的对比度值。

3.根据权利要求2所述的图像处理方法，其特征在于，生成所述第三图像的步骤包括：

获得所述第一图像的色调转换曲线的步骤；

计算步骤，用于计算所述第二图像的色调转换曲线；

合成步骤，用于基于反馈系数(α)，对所述第一图像的色调转换曲线和所述第二图像的色调转换曲线进行合成，以获得合成色调转换曲线，从而生成所述第三图像。

4.根据权利要求3所述的图像处理方法，其特征在于，所述反馈系数(α)包括依赖于与所述基准区域的距离进行加权的乘数。

5.根据权利要求3所述的图像处理方法，其特征在于，根据以下等式计算所述反馈系数(α)：

x≤x_basis：α＝k₁x³+k₂x²+k₃x+α_min

x＞x_basis：α＝k₁(x_max-x)³+k₂(x_max-x)²+k₃(x_max-x)+α_min，

其中，x是输入像素值，x_basis是所述基准区域中的中间像素，x_max是原始图像中的最大像素值，α_min是最小反馈系数值，以及k是根据x与x_basis之间的距离进行加权的加权系数。

6.根据权利要求1所述的图像处理方法，其特征在于，还包括：

合成步骤，用于基于反馈系数(α)，对所述图像序列的所述第一图像的色调转换曲线和所述第二图像的色调转换曲线进行合成，其中，与针对所述基准区域以外的区域的反馈系数相比，针对所述基准区域的反馈系数被设置为较大的值，

其中，在所述色调转换步骤中，使用在所述合成步骤中合成得到的合成色调转换曲线来进行色调转换。

7.根据权利要求1所述的图像处理方法，其特征在于，所述基准区域是包含所关注的解剖学区域的图像的图像区域。

8.根据权利要求1所述的图像处理方法，其特征在于，通过摄像区域信息或摄像技术信息来确定所述基准区域的位置。

9.根据权利要求1所述的图像处理方法，其特征在于，所述基准区域是使用统计学要素确定出的区域。

10.一种图像处理设备，用于对由X射线摄像系统所生成的图像序列的图像进行色调转换，所述图像处理设备包括：

提取部件，用于从所述图像序列的第二图像中提取在所述图像序列的图像中用作基准区域的区域；以及

色调转换部件，用于在相对于所述图像序列中的第一图像在所述第二图像中存在对比度变化的情况下，对所述第二图像进行色调转换，从而在所述基准区域中抑制所述对比度变化，并且在所述基准区域以外的区域中反映所述对比度变化。

11.根据权利要求10所述的图像处理设备，其特征在于，还包括：

获得部件，用于获得所述第一图像中的多个像素的对比度值，并且获得所述第二图像中的多个像素的对比度值，

其中，所述色调转换部件用于生成如下第三图像，其中，所述第三图像在所述基准区域中的对比度值接近所述第一图像中的多个像素的对比度值，并且所述第三图像在所述基准区域以外的区域中的对比度值接近所述第二图像中的多个像素的对比度值。