射线图像诊断装置1的设定的用户，一边目视 窗口 500, 一边在该山峰的右端附近设定第一阈值502即可。
[0086] 这种情况下，本实施方式的X射线图像诊断装置1如图1所示还具备操作输入部 14。并且，通过该操作输入部14等来进行该设定。即，在本实施方式的X射线图像诊断装 置1中，能够在从显示输出部11提示上述直方图的同时，输入来自用户的图像处理部10的 处理的设定变更。
[0087] 此外，窗口 510是输入了运动/变形的影像的情况下的直方图的例子。期望相似 度直方图511也集中在相似度大的区域。由于影像发生运动/变形，因此该山峰表示相似 度变大后的区域。即，如果针对该山峰应用空间方向噪声减低，则能够有效果地去除噪声。 因此，用户一边目视窗口 510,一边在该山峰的左端附近设定第二阈值512即可。
[0088] 如上所述，本实施方式的放射线图像诊断装置1具备：照射放射线的放射线照射 部；检测所述放射线照射部所照射出的放射线而生成图像的检测部；针对由所述检测部生 成的所述图像进行图像处理的图像处理部10。所述图像处理部10具备：移动检测部201， 其使用由所述检测部生成的拍摄时间不同的多个图像中的、时间上靠前的图像和时间上靠 后的图像，进行移动检测，针对一张图像检测一个移动信息；相似度计算部202,其针对每 个局部区域，计算将所述时间上靠前的图像仅变更所述移动信息的量的位置后的图像与所 述时间上靠后的图像的差值，作为相似度；区域判定部，其根据所述相似度进行区域判定， 该区域判定决定在所述图像上进行时间方向的噪声降低处理的区域和进行空间方向的噪 声降低处理的区域；噪声降低处理部，其根据所述区域判定部的所述区域判定结果，针对各 区域进行所述噪声降低处理。
[0089] 并且，所述区域判定部通过针对所述相似度进行阈值处理，决定进行时间方向的 噪声降低处理的区域和进行空间方向的噪声降低处理的区域。
[0090] 进行所述时间方向的噪声降低处理的区域和进行所述空间方向的噪声降低处理 的区域可以一部分或全部重复。
[0091] 根据以上说明的本实施方式的放射线图像诊断装置1以及该图像处理方法，能够 更恰当地选择进行时间方向的噪声降低处理和空间方向的噪声降低处理的每一个或进行 这两者的区域，并能够更恰当地兼顾降低噪声和抑制模糊发生。由此，能够实现更佳的放射 线透视图像的噪声降低。
[0092] 此外，在以上说明的第一实施方式中，在进行噪声降低处理的前段处理中，通过大 范围的移动检测仅检测出一个大范围的移动信息。然而，本实施方式并不妨碍放射线图像 诊断装置1在上述处理的前后进行检测局部移动的处理。此外，本实施方式并不妨碍放射 线图像诊断装置1具有进行使用局部移动的检测的其他噪声降低处理的模式。
[0093] 〈〈第二实施方式》
[0094] 本实施方式的X射线图像诊断装置基本上具有与第一实施方式的X射线图像诊断 装置1同样的结构。但是，图像处理部10的处理不同。以下，针对本实施方式，着眼于与第 一实施方式不同的结构进行说明。
[0095] 图6是表示第二实施方式的图像处理部10的处理框图。在第二实施方式的图像 处理部10中，在第一区域和第二区域的区域判定中，针对相似度计算部202计算出的相似 度，进行使用LUT(L〇〇kUpTable:查找表）的变换处理。LUT可以预先分别存储图像处理 部10的第一区域判定部600或第二区域判定部601，也可以存储在图像处理部10内未图示 的存储部。
[0096] 省略图6的图像处理部10中的、已经说明的图2所示的附带相同符号的结构和具 有相同功能的部分的说明。第一区域判定部600针对相似度计算部202计算出的相似度，进 行使用第一LUT的变换来决定第一区域。考虑图像的噪声水平等来预先设定第一LUT。数 学式（10)是坐标X处的第一区域R、b的计算式，Sa，b表示相似度，L1表示基于第一LUT的 变换。
[0097][式 10]
[0099] 在第一实施方式中，第一区域R^b取2值，然而在本实施方式中取0至1的多个值 (实数）。第一LUT作为整体指定斜率为负的形状。图7是第一LUT的一例。LUT的设定方 法是任意的，可以是用于一般的LUT的设定的基于η点指定的曲线近似。
[0100] 如果更简便地进行设定，则可以使用S形（sigmoid)函数生成LUT形状。数学式 (11)是表示第一LUT即L1的形状的式子，X是LUT的输入值，β郴β2是对LUT的形状进 行控制的参数。
[0101] [式 11]
[0103] 第二区域判定部601针对相似度计算部202计算出的相似度，进行使用第二 LUT(查找表）的变换来决定第二区域。考虑图像的噪声水平等来预先设定第二LUT。数学 式（12)是坐标X处的第二区域R2a,b的计算式，Sa,b表示相似度，L2表示基于第二LUT的变 换。
[0104][式 12]
[0106] 在第一实施方式中，第二区域R\b取2值，然而在本实施方式中取0至1的多个值 (实数）。第二LUT作为整体指定斜率为正的形状。图8是第二LUT的一例。LUT的设定方 法是任意的，可以是用于一般的LUT的设定的基于η点指定的曲线近似。
[0107] 如果更简便地进行设定，则可以使用S形函数生成LUT形状。数学式（13)是表示 第二LUT即L2的形状的式子，X是LUT的输入值，β郴β2是对LUT的形状进行控制的参 数。
[0108][式 13]
[0110] 显示输出部11进行去除噪声后的图像的显示。也可以代替该显示或与该显示包 括该显示，可视化地显示相似度直方图以及第一LUT或第二LUT。
[0111] 图9是将相似度直方图和第一LUT可视化时的显示例。
[0112] 在窗口 900中，配置表现相似度直方图的棒图表901，以及表现第一LUT或第二 LUT的曲线902。
[0113] 如以上那样，通过使用LUT，基于相似度，能够更细化地表现第一区域划分和第二 区域划分。图10(a)-图10(c)中表不本实施方式的第一区域划分和第二区域划分的一例。
[0114] 图10(a)-图10(c)是表不第一区域划分和第二区域划分、以及应用时间方向噪声 去除和空间方向噪声去除的区域的示意图。例如，如图示，检测为第一区域划分1000的中 央附近为"0"，除此以外为"1"，第二区域划分1001的中央附近为"1"，除此以外为"0"。
[0115] 此时，应用各噪声降低处理的范围1002如以下。被检测为第一区域划分1000的 "0"的区域和被检测为第二区域划分1001的"1"的区域发生重叠的区域1011，应用空间方 向噪声降低处理，不应用时间方向噪声降低处理。被检测为第一区域划分1000的"1"的区 域和被检测为第二区域划分1001的"〇"的区域发生重叠的区域1013,应用时间方向噪声降 低处理，不应用空间方向噪声降低处理。被检测为第一区域划分1000的"1"的区域和被检 测为第二区域划分1001的"1"的区域发生重叠的区域1012,应用时间方向噪声降低处理和 空间方向噪声降低处理。此外，被检测为第一区域划分1000的"〇"的区域和被检测为第二 区域划分1001的"〇"的区域发生重叠的区域1014,不进行噪声降低。
[0116] 如图10(a)-图10(c)中所示的例子，本实施方式的处理能够比第一实施方式更 细化地设定第一区域划分和第二区域划分，例如能够设定以下4个区域：（1)进行时间方向 的噪声降低处理，不进行空间方向的噪声降低处理的区域，（2)进行空间方向的噪声降低处 理，不进行空间方向的噪声降低处理的区域，（3)进行时间方向的噪声降低处理和空间方向 的噪声降低处理两者的区域，以及（4)不进行时间方向的噪声降低处理和空间方向的噪声 降低处理两者的区域。
[0117] 如以上所述，与第一实施方式同样地，本实施方式的放射线图像诊断装置1具备： 放射线照射部、检测部、图像处理部10。所述区域判定部针对所述相似度使用预先由所述放 射线图像生成装置1所存储的查找表，由此，决定进行时间方向的噪声降低处理的区域和 进行空间方向的噪声降低处理的区域。
[0118] 进行所述时间方向的噪声降低处理的区域和进行所述空间方向的噪声降低处理 的区域可以一部分或全部重复。
[0119] 通过以上结构，第一区域和第二区域能够以多值表现，并能够通过第一LUT和第 二LUT更细化地调整该值。因此，通过分别独立地设定应用时间方向噪声降低和空间方向 噪声降低的范围和强度，能够更有效地降低X射线透视图像的噪声。
[0120] 此外，如第一实施方式中说明那样，一边观察相似度直方图一边对第一阈值和第 二阈值进行设定是恰当的。用户也通过将第一LUT和第二LUT与相似度直方图重叠，能够 视觉上变更各LUT的形状，并能够针对图像选择恰当的噪声降低处理。
[0121] 最后，说明上述各实施方式的基于图像处理部10的图像处理的流程。图11是上 述实施方式的图像处理的处理流程。这里，以第一实施方式的结构为例进行说明。
[0122] 图像处理部10,针对使用放射线拍摄到的图像，接收拍摄时间不同的多个图像的 输入（步骤S1101)。
[0123] 移动检测部201使用图像处理部10接收到的多个图像中的、时间上靠前的图像和 时间上靠后的图像来进行移动检测，针对一张图像检测一个移动信息（步骤S1102)。
[0124] 相似度计算部202针对每个局部区域，计算将时间上靠前的图像仅变更所述移动 信息的量的位置后的图像与所述时间上靠后的图像的差值，作为相似度（步骤1103)。
[0125] 第一区域判定部203和第二区域判定部205根据所述相似度来决定拍摄的图像的 第一区域和第二区域（步骤