射線圖像診斷裝置1的設(shè)定的用戶,一邊目視 窗口 500, 一邊在該山峰的右端附近設(shè)定第一閾值502即可。
[0086] 這種情況下,本實(shí)施方式的X射線圖像診斷裝置1如圖1所示還具備操作輸入部 14。并且,通過該操作輸入部14等來進(jìn)行該設(shè)定。即,在本實(shí)施方式的X射線圖像診斷裝 置1中,能夠在從顯示輸出部11提示上述直方圖的同時,輸入來自用戶的圖像處理部10的 處理的設(shè)定變更。
[0087] 此外,窗口 510是輸入了運(yùn)動/變形的影像的情況下的直方圖的例子。期望相似 度直方圖511也集中在相似度大的區(qū)域。由于影像發(fā)生運(yùn)動/變形,因此該山峰表示相似 度變大后的區(qū)域。即,如果針對該山峰應(yīng)用空間方向噪聲減低,則能夠有效果地去除噪聲。 因此,用戶一邊目視窗口 510,一邊在該山峰的左端附近設(shè)定第二閾值512即可。
[0088] 如上所述,本實(shí)施方式的放射線圖像診斷裝置1具備:照射放射線的放射線照射 部;檢測所述放射線照射部所照射出的放射線而生成圖像的檢測部;針對由所述檢測部生 成的所述圖像進(jìn)行圖像處理的圖像處理部10。所述圖像處理部10具備:移動檢測部201, 其使用由所述檢測部生成的拍攝時間不同的多個圖像中的、時間上靠前的圖像和時間上靠 后的圖像,進(jìn)行移動檢測,針對一張圖像檢測一個移動信息;相似度計(jì)算部202,其針對每 個局部區(qū)域,計(jì)算將所述時間上靠前的圖像僅變更所述移動信息的量的位置后的圖像與所 述時間上靠后的圖像的差值,作為相似度;區(qū)域判定部,其根據(jù)所述相似度進(jìn)行區(qū)域判定, 該區(qū)域判定決定在所述圖像上進(jìn)行時間方向的噪聲降低處理的區(qū)域和進(jìn)行空間方向的噪 聲降低處理的區(qū)域;噪聲降低處理部,其根據(jù)所述區(qū)域判定部的所述區(qū)域判定結(jié)果,針對各 區(qū)域進(jìn)行所述噪聲降低處理。
[0089] 并且,所述區(qū)域判定部通過針對所述相似度進(jìn)行閾值處理,決定進(jìn)行時間方向的 噪聲降低處理的區(qū)域和進(jìn)行空間方向的噪聲降低處理的區(qū)域。
[0090] 進(jìn)行所述時間方向的噪聲降低處理的區(qū)域和進(jìn)行所述空間方向的噪聲降低處理 的區(qū)域可以一部分或全部重復(fù)。
[0091] 根據(jù)以上說明的本實(shí)施方式的放射線圖像診斷裝置1以及該圖像處理方法,能夠 更恰當(dāng)?shù)剡x擇進(jìn)行時間方向的噪聲降低處理和空間方向的噪聲降低處理的每一個或進(jìn)行 這兩者的區(qū)域,并能夠更恰當(dāng)?shù)丶骖櫧档驮肼暫鸵种颇:l(fā)生。由此,能夠?qū)崿F(xiàn)更佳的放射 線透視圖像的噪聲降低。
[0092] 此外,在以上說明的第一實(shí)施方式中,在進(jìn)行噪聲降低處理的前段處理中,通過大 范圍的移動檢測僅檢測出一個大范圍的移動信息。然而,本實(shí)施方式并不妨礙放射線圖像 診斷裝置1在上述處理的前后進(jìn)行檢測局部移動的處理。此外,本實(shí)施方式并不妨礙放射 線圖像診斷裝置1具有進(jìn)行使用局部移動的檢測的其他噪聲降低處理的模式。
[0093] 〈〈第二實(shí)施方式》
[0094] 本實(shí)施方式的X射線圖像診斷裝置基本上具有與第一實(shí)施方式的X射線圖像診斷 裝置1同樣的結(jié)構(gòu)。但是,圖像處理部10的處理不同。以下,針對本實(shí)施方式,著眼于與第 一實(shí)施方式不同的結(jié)構(gòu)進(jìn)行說明。
[0095] 圖6是表示第二實(shí)施方式的圖像處理部10的處理框圖。在第二實(shí)施方式的圖像 處理部10中,在第一區(qū)域和第二區(qū)域的區(qū)域判定中,針對相似度計(jì)算部202計(jì)算出的相似 度,進(jìn)行使用LUT(L〇〇kUpTable:查找表)的變換處理。LUT可以預(yù)先分別存儲圖像處理 部10的第一區(qū)域判定部600或第二區(qū)域判定部601,也可以存儲在圖像處理部10內(nèi)未圖示 的存儲部。
[0096] 省略圖6的圖像處理部10中的、已經(jīng)說明的圖2所示的附帶相同符號的結(jié)構(gòu)和具 有相同功能的部分的說明。第一區(qū)域判定部600針對相似度計(jì)算部202計(jì)算出的相似度,進(jìn) 行使用第一LUT的變換來決定第一區(qū)域??紤]圖像的噪聲水平等來預(yù)先設(shè)定第一LUT。數(shù) 學(xué)式(10)是坐標(biāo)X處的第一區(qū)域R、b的計(jì)算式,Sa,b表示相似度,L1表示基于第一LUT的 變換。
[0097][式 10]
[0099] 在第一實(shí)施方式中,第一區(qū)域R^b取2值,然而在本實(shí)施方式中取0至1的多個值 (實(shí)數(shù))。第一LUT作為整體指定斜率為負(fù)的形狀。圖7是第一LUT的一例。LUT的設(shè)定方 法是任意的,可以是用于一般的LUT的設(shè)定的基于η點(diǎn)指定的曲線近似。
[0100] 如果更簡便地進(jìn)行設(shè)定,則可以使用S形(sigmoid)函數(shù)生成LUT形狀。數(shù)學(xué)式 (11)是表示第一LUT即L1的形狀的式子,X是LUT的輸入值,β郴β2是對LUT的形狀進(jìn) 行控制的參數(shù)。
[0101] [式 11]
[0103] 第二區(qū)域判定部601針對相似度計(jì)算部202計(jì)算出的相似度,進(jìn)行使用第二 LUT(查找表)的變換來決定第二區(qū)域??紤]圖像的噪聲水平等來預(yù)先設(shè)定第二LUT。數(shù)學(xué) 式(12)是坐標(biāo)X處的第二區(qū)域R2a,b的計(jì)算式,Sa,b表示相似度,L2表示基于第二LUT的變 換。
[0104][式 12]
[0106] 在第一實(shí)施方式中,第二區(qū)域R\b取2值,然而在本實(shí)施方式中取0至1的多個值 (實(shí)數(shù))。第二LUT作為整體指定斜率為正的形狀。圖8是第二LUT的一例。LUT的設(shè)定方 法是任意的,可以是用于一般的LUT的設(shè)定的基于η點(diǎn)指定的曲線近似。
[0107] 如果更簡便地進(jìn)行設(shè)定,則可以使用S形函數(shù)生成LUT形狀。數(shù)學(xué)式(13)是表示 第二LUT即L2的形狀的式子,X是LUT的輸入值,β郴β2是對LUT的形狀進(jìn)行控制的參 數(shù)。
[0108][式 13]
[0110] 顯示輸出部11進(jìn)行去除噪聲后的圖像的顯示。也可以代替該顯示或與該顯示包 括該顯示,可視化地顯示相似度直方圖以及第一LUT或第二LUT。
[0111] 圖9是將相似度直方圖和第一LUT可視化時的顯示例。
[0112] 在窗口 900中,配置表現(xiàn)相似度直方圖的棒圖表901,以及表現(xiàn)第一LUT或第二 LUT的曲線902。
[0113] 如以上那樣,通過使用LUT,基于相似度,能夠更細(xì)化地表現(xiàn)第一區(qū)域劃分和第二 區(qū)域劃分。圖10(a)-圖10(c)中表不本實(shí)施方式的第一區(qū)域劃分和第二區(qū)域劃分的一例。
[0114] 圖10(a)-圖10(c)是表不第一區(qū)域劃分和第二區(qū)域劃分、以及應(yīng)用時間方向噪聲 去除和空間方向噪聲去除的區(qū)域的示意圖。例如,如圖示,檢測為第一區(qū)域劃分1000的中 央附近為"0",除此以外為"1",第二區(qū)域劃分1001的中央附近為"1",除此以外為"0"。
[0115] 此時,應(yīng)用各噪聲降低處理的范圍1002如以下。被檢測為第一區(qū)域劃分1000的 "0"的區(qū)域和被檢測為第二區(qū)域劃分1001的"1"的區(qū)域發(fā)生重疊的區(qū)域1011,應(yīng)用空間方 向噪聲降低處理,不應(yīng)用時間方向噪聲降低處理。被檢測為第一區(qū)域劃分1000的"1"的區(qū) 域和被檢測為第二區(qū)域劃分1001的"〇"的區(qū)域發(fā)生重疊的區(qū)域1013,應(yīng)用時間方向噪聲降 低處理,不應(yīng)用空間方向噪聲降低處理。被檢測為第一區(qū)域劃分1000的"1"的區(qū)域和被檢 測為第二區(qū)域劃分1001的"1"的區(qū)域發(fā)生重疊的區(qū)域1012,應(yīng)用時間方向噪聲降低處理和 空間方向噪聲降低處理。此外,被檢測為第一區(qū)域劃分1000的"〇"的區(qū)域和被檢測為第二 區(qū)域劃分1001的"〇"的區(qū)域發(fā)生重疊的區(qū)域1014,不進(jìn)行噪聲降低。
[0116] 如圖10(a)-圖10(c)中所示的例子,本實(shí)施方式的處理能夠比第一實(shí)施方式更 細(xì)化地設(shè)定第一區(qū)域劃分和第二區(qū)域劃分,例如能夠設(shè)定以下4個區(qū)域:(1)進(jìn)行時間方向 的噪聲降低處理,不進(jìn)行空間方向的噪聲降低處理的區(qū)域,(2)進(jìn)行空間方向的噪聲降低處 理,不進(jìn)行空間方向的噪聲降低處理的區(qū)域,(3)進(jìn)行時間方向的噪聲降低處理和空間方向 的噪聲降低處理兩者的區(qū)域,以及(4)不進(jìn)行時間方向的噪聲降低處理和空間方向的噪聲 降低處理兩者的區(qū)域。
[0117] 如以上所述,與第一實(shí)施方式同樣地,本實(shí)施方式的放射線圖像診斷裝置1具備: 放射線照射部、檢測部、圖像處理部10。所述區(qū)域判定部針對所述相似度使用預(yù)先由所述放 射線圖像生成裝置1所存儲的查找表,由此,決定進(jìn)行時間方向的噪聲降低處理的區(qū)域和 進(jìn)行空間方向的噪聲降低處理的區(qū)域。
[0118] 進(jìn)行所述時間方向的噪聲降低處理的區(qū)域和進(jìn)行所述空間方向的噪聲降低處理 的區(qū)域可以一部分或全部重復(fù)。
[0119] 通過以上結(jié)構(gòu),第一區(qū)域和第二區(qū)域能夠以多值表現(xiàn),并能夠通過第一LUT和第 二LUT更細(xì)化地調(diào)整該值。因此,通過分別獨(dú)立地設(shè)定應(yīng)用時間方向噪聲降低和空間方向 噪聲降低的范圍和強(qiáng)度,能夠更有效地降低X射線透視圖像的噪聲。
[0120] 此外,如第一實(shí)施方式中說明那樣,一邊觀察相似度直方圖一邊對第一閾值和第 二閾值進(jìn)行設(shè)定是恰當(dāng)?shù)?。用戶也通過將第一LUT和第二LUT與相似度直方圖重疊,能夠 視覺上變更各LUT的形狀,并能夠針對圖像選擇恰當(dāng)?shù)脑肼暯档吞幚怼?br>[0121] 最后,說明上述各實(shí)施方式的基于圖像處理部10的圖像處理的流程。圖11是上 述實(shí)施方式的圖像處理的處理流程。這里,以第一實(shí)施方式的結(jié)構(gòu)為例進(jìn)行說明。
[0122] 圖像處理部10,針對使用放射線拍攝到的圖像,接收拍攝時間不同的多個圖像的 輸入(步驟S1101)。
[0123] 移動檢測部201使用圖像處理部10接收到的多個圖像中的、時間上靠前的圖像和 時間上靠后的圖像來進(jìn)行移動檢測,針對一張圖像檢測一個移動信息(步驟S1102)。
[0124] 相似度計(jì)算部202針對每個局部區(qū)域,計(jì)算將時間上靠前的圖像僅變更所述移動 信息的量的位置后的圖像與所述時間上靠后的圖像的差值,作為相似度(步驟1103)。
[0125] 第一區(qū)域判定部203和第二區(qū)域判定部205根據(jù)所述相似度來決定拍攝的圖像的 第一區(qū)域和第二區(qū)域(步驟