像13和Ib間的移動(dòng)方向。數(shù) 學(xué)式⑵是以輸入圖像Ia為基準(zhǔn)的輸入圖像I3和Ib間的移動(dòng)信息Vaib的計(jì)算式，Ψ表示 搜索范圍的坐標(biāo)的集合。搜索范圍典型地指定水平20像素、垂直20像素程度的矩形區(qū)域的 坐標(biāo)。但是，也可以指定考慮了輸入圖像中設(shè)想的被檢體的最大移動(dòng)速度、計(jì)算數(shù)學(xué)式（1) 和數(shù)學(xué)式（2)的速度而決定的區(qū)域。
[0039] [式 2]
[0041] 相似度計(jì)算部202以幀存儲(chǔ)器200中存儲(chǔ)的2張輸入圖像和由移動(dòng)檢測部201計(jì) 算出的移動(dòng)信息作為輸入，在將拍攝時(shí)間不同的2張輸入圖像中時(shí)間上靠前的輸入圖像僅 偏移該移動(dòng)信息的量的位置后，針對每個(gè)局部區(qū)域計(jì)算并輸出圖像間的相似度。
[0042] 更具體地，在各坐標(biāo)計(jì)算局部區(qū)域的SAD并將其作為相似度。因此，相似度的值越 小，圖像越相似。數(shù)學(xué)式⑶是輸入圖像13和I別的坐標(biāo)X處的相似度sa，b的計(jì)算式，ΩX 是坐標(biāo)X處的局部區(qū)域的坐標(biāo)的集合，Va，b是由移動(dòng)檢測部201計(jì)算出的圖像13和I,間間 的移動(dòng)信息。針對Ωχ典型地指定以坐標(biāo)X為中心的水平7像素、垂直7像素的坐標(biāo)的集 合。但是，也可以考慮圖像中的物體的圖案的細(xì)節(jié)、圖像的噪聲水平、計(jì)算速度等，變更其范 圍和形狀。
[0043] [式 3]
[0045] 第一區(qū)域判定部203和第二區(qū)域判定部205,是根據(jù)相似度計(jì)算部202計(jì)算出的相 似度進(jìn)行區(qū)域判定的區(qū)域判定部，該區(qū)域判定決定在輸入圖像上進(jìn)行時(shí)間方向的噪聲降低 處理的區(qū)域和進(jìn)行空間方向的噪聲降低處理的區(qū)域。這里，將輸入圖像劃分為（1)進(jìn)行時(shí) 間方向的噪聲降低處理，不進(jìn)行空間方向的噪聲降低處理的區(qū)域，（2)進(jìn)行空間方向的噪聲 降低處理，不進(jìn)行時(shí)間方向的噪聲降低處理的區(qū)域，以及（3)進(jìn)行時(shí)間方向的噪聲降低處 理和空間方向的噪聲降低處理兩者的區(qū)域這3個(gè)區(qū)域中的至少2個(gè)區(qū)域。
[0046] 第一區(qū)域判定部203針對相似度計(jì)算部202計(jì)算出的相似度Sa』，進(jìn)行使用第一閾 值的閾值處理，決定第一區(qū)域劃分。考慮圖像的噪聲水平等預(yù)先設(shè)定第一閾值。數(shù)學(xué)式（4) 是坐標(biāo)X處的第一區(qū)域劃分R\b的計(jì)算式，Saib表示相似度，T1表示第一閾值。
[0047] [式 4]
[0049] 第一區(qū)域劃分為"1"的坐標(biāo)的集合表示相似度?。▓D像相似）的區(qū)域。也就 是說，是當(dāng)使用移動(dòng)檢測部201檢測出的移動(dòng)信息Va，b對輸入圖像13和I,進(jìn)行了移動(dòng)補(bǔ)償 時(shí)，接近靜止的狀態(tài)的區(qū)域。即，以提取接近靜止?fàn)顟B(tài)的區(qū)域的方式，設(shè)定第一閾值T1。
[0050] 第二區(qū)域判定部205針對相似度計(jì)算部202計(jì)算出的相似度Saib，進(jìn)行使用第一閾 值的閾值處理，決定第二區(qū)域劃分?？紤]圖像的噪聲水平等預(yù)先設(shè)定第二閾值。數(shù)學(xué)式（5) 是坐標(biāo)X處的第二區(qū)域劃分R\b的計(jì)算式，Saib表示相似度，T2表示第二閾值。
[0051][式 5]
[0053] 第二區(qū)域劃分R2a,b為"1"的坐標(biāo)的集合表示相似度大（圖像不相似）的區(qū)域。也 就是說，是當(dāng)使用移動(dòng)檢測部201檢測出的移動(dòng)信息對輸入圖像1 3和Ib進(jìn)行了移動(dòng)補(bǔ) 償時(shí)，進(jìn)行運(yùn)動(dòng)、變形的區(qū)域。即，以提取進(jìn)行運(yùn)動(dòng)、變形的區(qū)域的方式，設(shè)定第二閾值T2。
[0054] 時(shí)間方向噪聲降低部204和空間方向噪聲降低部206,是根據(jù)第一區(qū)域判定部203 和第二區(qū)域判定部205的判定結(jié)果，針對各區(qū)域進(jìn)行噪聲降低處理的噪聲降低處理部。
[0055] 時(shí)間方向噪聲降低部204輸入幀存儲(chǔ)器中存儲(chǔ)的多張輸入圖像、由移動(dòng)檢測部 201計(jì)算出的移動(dòng)信息和由第一區(qū)域判定部203決定的第一區(qū)域劃分。使用移動(dòng)信息對輸 入的多張圖像進(jìn)行移動(dòng)補(bǔ)償，同時(shí)針對第一區(qū)域劃分為1的區(qū)域計(jì)算輸入圖像的幀之間的 加權(quán)平均，由此進(jìn)行噪聲降低，并輸出降低噪聲后的圖像。
[0056] 數(shù)學(xué)式（6)是進(jìn)行了坐標(biāo)X處的時(shí)間方向噪聲降低的圖像ITIME的計(jì)算式， Ιη(1 <η<A)表示輸入圖像，cn(l<η<A)表示每個(gè)幀的加權(quán)系數(shù)，\"表示最新的幀與 各幀間的移動(dòng)，R\in表示最新的幀與各幀間的相似度。
[0057][式 6]
[0060] 這里，如果輸入圖像^至IA是時(shí)間上連續(xù)的幀，則當(dāng)利用移動(dòng)檢測部201計(jì)算鄰 接的幀之間的移動(dòng)信息\2、v2,3、…、v(A1)，A時(shí)，利用這些的累計(jì)加法運(yùn)算求取Vlin。數(shù)學(xué)式 (7)是輸入圖像^至In間的移動(dòng)信息V的計(jì)算式。
[0061][式7]
[0063] 其中，νΜ= 0(無移動(dòng)）
[0064] 此外，以總和為1的方式預(yù)先設(shè)置加權(quán)系數(shù)cn。數(shù)學(xué)式（8)是加權(quán)系數(shù)cn的限制 式。
[0065][式 8]
[0067] 如以上，本實(shí)施方式的時(shí)間方向噪聲降低部204使用由移動(dòng)檢測部201檢測出的 大范圍的移動(dòng)信息針對輸入的多張圖像進(jìn)行移動(dòng)補(bǔ)償，同時(shí)通過輸入圖像的幀間的加權(quán)平 均實(shí)現(xiàn)噪聲的降低。因此，相比于未使用移動(dòng)信息地進(jìn)行時(shí)間方向的加權(quán)平均的情況，能夠 更加恰當(dāng)?shù)亟档驮肼暋?br>[0068] 空間方向噪聲降低部206以由時(shí)間方向噪聲降低部204針對第一區(qū)域劃分為1的 區(qū)域進(jìn)行時(shí)間方向噪聲去除后的圖像、以及由第二區(qū)域判定部205決定的第二區(qū)域劃分作 為輸入，針對輸入的圖像中的第二區(qū)域劃分為1的區(qū)域進(jìn)行空間方向的噪聲降低，并輸出 噪聲降低后的圖像。
[0069] 作為空間方向的噪聲降低法，采用難以發(fā)生在時(shí)間方向的噪聲降低處理中可能發(fā) 生的物體運(yùn)動(dòng)所導(dǎo)致的模糊的方法。即，當(dāng)使用僅用輸入的1張圖像來除去噪聲的方法時(shí)， 原理上不會(huì)發(fā)生該模糊。例如，高斯濾波器、雙邊濾波器、非局部均值濾波器、小波收縮等方 法是合適的。但是，如果是針對由物體的運(yùn)動(dòng)導(dǎo)致的模糊很強(qiáng)的噪聲降低方法，則也可以使 用利用多個(gè)幀等其他方法。
[0070] 以下，作為一例，說明使用高斯濾波器進(jìn)行空間方向的噪聲降低的例子。空間方向 的噪聲降低應(yīng)用于進(jìn)行了時(shí)間方向的噪聲降低處理后的圖像中的第二區(qū)域劃分為1的區(qū) ±或。數(shù)學(xué)式（9)是生成坐標(biāo)X處的進(jìn)行空間方向噪聲降低后的圖像ISPAeE的計(jì)算式，1?， 示進(jìn)行時(shí)間方向的噪聲去除后的圖像，R\2表示輸入圖像I:和I2間的第二區(qū)域劃分，α表 示高斯核（gaussiankernel)的大小（對噪聲除去的強(qiáng)度進(jìn)行控制的參數(shù)），Σ表示應(yīng)用 高斯核的范圍的坐標(biāo)的集合。考慮圖像的噪聲水平、數(shù)學(xué)式（9)的計(jì)算速度，決定高斯核的 大小α和高斯核的應(yīng)用范圍Σ。
[0071][式 9]
[0075] 相似度直方圖計(jì)算部207計(jì)算由相似度計(jì)算部202計(jì)算出的相似度的頻數(shù)分布， 并將其輸出為相似度直方圖。直方圖的小區(qū)間（bin)的數(shù)目為輸入圖像的灰階數(shù)程度。但 是，可以考慮后述的直方圖窗口的尺寸（像素?cái)?shù)）、第一（第二）閾值的設(shè)定精度來進(jìn)行改 變。
[0076] 顯示輸出部11除了除去噪聲后的圖像之外，還可視化地顯示相似度直方圖、以及 第一閾值或第二閾值。
[0077] 圖3是將相似度直方圖可視化時(shí)的顯示例。
[0078] 窗口 300中，在表現(xiàn)了相似度直方圖的棒圖表301上，重疊地配置表現(xiàn)第一閾值或 第二閾值的縱線302。
[0079] 棒圖表301以橫軸為相似度，以縱軸為頻數(shù)，表現(xiàn)相似度直方圖?？v線302與棒圖 表301的橫軸（相似度）的位置相符合地表現(xiàn)第一（第二）閾值。
[0080] 通過以上結(jié)構(gòu)，本實(shí)施方式的圖像處理部10能夠通過第一閾值和第二閾值的調(diào) 整，在由檢測部生成的圖像中設(shè)定以下3個(gè)區(qū)域：⑴進(jìn)行時(shí)間方向的噪聲降低處理，不進(jìn) 行空間方向的噪聲降低處理的區(qū)域，（2)進(jìn)行空間方向的噪聲降低處理，不進(jìn)行時(shí)間方向的 噪聲降低處理的區(qū)域，以及（3)進(jìn)行時(shí)間方向的噪聲降低處理和空間方向的噪聲降低處理 兩者的區(qū)域。
[0081] 圖4(a) -圖4(c)是表示第一區(qū)域劃分和第二區(qū)域劃分，以及應(yīng)用時(shí)間方向噪聲去 除和空間方向噪聲去除的區(qū)域的示意圖。例如，如圖示，作為第一區(qū)域劃分400,檢測為中 央附近為"〇"，除此以外為"1"，作為第二區(qū)域劃分401，檢測為中央附近為"1"，除此以外為 "0"。
[0082] 此時(shí)，應(yīng)用各噪聲降低處理的范圍402如以下。被檢測為第一區(qū)域劃分400的"0" 的區(qū)域和被檢測為第二區(qū)域劃分401的" 1"的區(qū)域發(fā)生重疊的區(qū)域411，應(yīng)用空間方向噪聲 降低處理，不應(yīng)用時(shí)間方向噪聲降低處理。被檢測為第一區(qū)域劃分400的"1"的區(qū)域和被 檢測為第二區(qū)域劃分401的"0"的區(qū)域發(fā)生重疊的區(qū)域413,應(yīng)用時(shí)間方向噪聲降低處理， 不應(yīng)用空間方向噪聲降低處理。被檢測為第一區(qū)域劃分400的"1"的區(qū)域和被檢測為第二 區(qū)域劃分401的"1"的區(qū)域發(fā)生重疊的區(qū)域412,應(yīng)用空間方向噪聲降低處理和時(shí)間方向噪 聲降低處理。
[0083]S卩，在本實(shí)施方式的處理中，通過第一閾值和第二閾值選擇，能夠恰當(dāng)?shù)卮_定時(shí)間 方向的噪聲降低有效果的區(qū)域（靜止的區(qū)域等）、空間方向的噪聲降低有效果的區(qū)域（進(jìn) 行運(yùn)動(dòng)或變形的區(qū)域等）、應(yīng)用時(shí)間、空間兩者的噪聲降低有效果的區(qū)域（略微變形的區(qū)域 等），并能夠針對各個(gè)區(qū)域應(yīng)用恰當(dāng)?shù)脑肼暯档吞幚怼?br>[0084] 此外，還能夠以視覺上更易懂的方式來進(jìn)行第一閾值和第二閾值的設(shè)定。圖5(a) 和圖5(b)圖示了這種情況的設(shè)定情形，是向X射線圖像診斷裝置1的顯示輸出部11輸出 的設(shè)定畫面的一部分（窗口）。該設(shè)定畫面可以與X射線攝像圖像的顯示重疊，也可以并列 顯示。此外，也可以不顯示X射線攝像圖像地單獨(dú)顯示。
[0085] 例如，窗口 500是輸入了靜止的影像的情況下的直方圖的例子。此時(shí)，期望相似度 直方圖501成為在相似度小的區(qū)域中集中的類型。該山峰表示影像中包含的噪聲成分。如 果是靜止的影像且未承載噪聲，則相似度全部為零。即，如果針對該山峰應(yīng)用時(shí)間方向噪聲 減低，則能夠有效地去除噪聲。因此，進(jìn)行X