牟希農(nóng)

(甘肅中醫(yī)藥大學(xué)定西校區(qū) 醫(yī)學(xué)教學(xué)部,甘肅 定西 743000)

超聲成像是現(xiàn)有醫(yī)學(xué)成像技術(shù)中常用的診斷方法，而超聲成像過程中由信號的差異與強(qiáng)弱產(chǎn)生的噪聲會降低超聲圖像的質(zhì)量，進(jìn)而會影響臨床醫(yī)生的判斷和治療[1-2]，因此研究超聲圖像去噪技術(shù)具有重要意義。近年來一些研究者將數(shù)字圖像濾波技術(shù)應(yīng)用到醫(yī)學(xué)影像圖像的去噪中，并取得了良好效果[3-4]。例如：肖磊等[5]提出了一種沿邊緣切線擴(kuò)散去噪的方法，該方法很好地去除了醫(yī)學(xué)超聲圖像的邊緣噪聲，較好地保持了圖像細(xì)節(jié)信息，克服了多次迭代運(yùn)算的煩瑣，但對其中心區(qū)域的噪聲抑制不理想。袁樂民[6]運(yùn)用中值濾波算法對含有噪聲的醫(yī)學(xué)圖像去噪方法進(jìn)行了探討，取得了好的濾波效果，但該方法對高密度噪聲抑制效果不好。梅玲等[7]采用加權(quán)中值濾波算法對MRI中的高密度椒鹽噪聲和高斯噪聲進(jìn)行了去除研究，該方法能夠較好地抑制高密度噪聲，有效地保護(hù)了圖像細(xì)節(jié)，但運(yùn)算量大。張雪峰等[8]提出了一種基于分?jǐn)?shù)階的中值濾波算法，該方法可以有效抑制圖像中的椒鹽噪聲，達(dá)到增強(qiáng)圖像功效，能較好地保留圖像邊緣突出的細(xì)節(jié)，但該算法對高斯噪聲不敏感。李瑞琦等[9]將小波變換和中值濾波相結(jié)合，能夠很好地消除高頻噪聲和椒鹽噪聲，但存在無法抑制大噪聲的問題。方斌等[10]采用小波閾值方法能夠有效抑制脈沖噪聲，該方法僅對椒鹽噪聲具有較好的峰值信噪比和邊緣保持指數(shù)，但對高斯噪聲不理想。陳軍[11]采用小波分析對醫(yī)學(xué)影像圖像的高斯噪聲抑制問題進(jìn)行了實驗研究，在小噪聲環(huán)境下效果較好，但噪聲密度較高時存在條紋效應(yīng)，使圖像邊緣產(chǎn)生間斷現(xiàn)象。何文[12]通過小波分析和模糊理論相結(jié)合的算法對超聲圖像的高頻噪聲進(jìn)行有效地抑制,使圖像質(zhì)量得到改善，但該方法存在閾值選擇的限制，使圖像清晰度受影響。傅博等[13]提出了邊緣導(dǎo)向的非局部圖像均值濾波算法，該方法能有效抑制高強(qiáng)度噪聲干擾環(huán)境下圖像的邊緣修復(fù)能力，但對圖像細(xì)節(jié)的保持性不高。劉傳義等[14]采用余弦相似度非局部均值濾波方法研究了圖像去噪和邊緣結(jié)構(gòu)信息保持問題，能夠反映圖像中的鄰域相似性，克服光照不均的影響。但算法僅在相對小噪聲干擾情況下去噪效果好，對大噪聲圖像效果不佳?；葆摰萚15]提出了自適應(yīng)去噪方法,該方法采用比率距離計算超聲圖像像素間的相似度距離，達(dá)到對濾波窗口的調(diào)整而有效地過濾了超聲圖像的斑點(diǎn)噪聲，提高了圖像的邊緣細(xì)節(jié)，但該算法對高斯噪聲的抑制不理想。基于上述研究，為進(jìn)一步提升去噪性能、減少運(yùn)算量，提出了一種基于閾值鄰域均值算法，并采用仿真軟件對含有高密度的椒鹽和高斯噪聲的胰腺超聲圖像去噪進(jìn)行探討，驗證了該算法的實效性。

1 基于閾值鄰域均值的算法分析

醫(yī)學(xué)影像圖像中多數(shù)像素的灰度與鄰近像素的灰度存在很大的灰度相關(guān)性。但噪聲具有隨機(jī)性，它產(chǎn)生的灰度與鄰點(diǎn)灰度無關(guān)，噪聲通常在其高頻部分呈現(xiàn)。所以，減弱高頻成分增強(qiáng)低頻成分是抑制噪聲的有效辦法。

假定有一幅m×n個像素的二維含噪圖像f(x,y),平滑處理后得到一幅圖像灰度值g(x,y)。則鄰域平均濾波算法可表示為[16]

(1)

式(1)中，x,y=0,1,2,…,N-1,M為S的像素總數(shù)，S為(x,y)空心鄰域中所有像素點(diǎn)坐標(biāo)的集合。平滑后的圖像g(x,y)中的每個像素的灰度值均由包含在(x,y)鄰域內(nèi)的像素灰度值的平均值來決定。h(x,y)為窗口函數(shù)，其灰度值g(x,y)的卷積表達(dá)式為

g(x,y)=f(x,y)×h(x,y)=

(2)

鄰域平均算法去噪具有算法簡便、運(yùn)算速度快的優(yōu)勢，但它會導(dǎo)致圖像細(xì)節(jié)丟失大而產(chǎn)生朦朧感。閾值鄰域均值算法可以克服鄰域平均算法的缺陷，提高圖像的細(xì)節(jié)信息[17]，其閾值鄰域平均法表示如下：

g(x,y)=

(3)

其中，T為閾值，由Bayes理論[18]有

(4)

(5)

(6)

如圖1所示為對含有高密度椒鹽和高斯噪聲的超聲圖像采用閾值鄰域均值去噪的算法主要程序流程。

2 閾值鄰域均值算法在醫(yī)學(xué)超聲圖像中的實現(xiàn)

為驗證閾值鄰域均值算法的有效性，運(yùn)用MATLAB平臺[19-20]將中值、均值、閾值鄰域均值等3種算法的去噪效果進(jìn)行了仿真實驗比較。實驗用圖像來源于醫(yī)學(xué)圖像數(shù)據(jù)集，如圖2為胰腺超聲原始圖，加入噪聲分別為高密度椒鹽和高斯等2種噪聲。對2種噪聲分別利用3種方法進(jìn)行噪聲去除實驗，實驗結(jié)果如圖3(a)—(d)和圖4(a)—(d)所示。

圖3(a)是在圖2中加有均值為0、方差為0.25的高斯噪聲的胰腺超聲圖像；圖3 (b)是對圖3(a)通過中值去噪后的圖像；圖3(c)是對圖3(a)通過均值去噪后的圖像；圖3(d)是對圖3(a)通過本文算法后的圖像。對比實驗結(jié)果，觀察到：中值非線性和均值線性濾波算法在濾除高斯噪聲方面效果不佳，而本研究中提出的算法能有效地抑制其高斯噪聲，能很好地保護(hù)醫(yī)學(xué)影像圖像的細(xì)節(jié)，使圖像更明晰。

圖4(a)是在圖2中加有噪聲密度為0.25椒鹽噪聲的胰腺超聲圖像；圖4(b)是對圖4(a)經(jīng)中值非線性去噪后的圖像；圖4(c)是對圖4(a)經(jīng)均值線性去噪后的圖像；圖4(d)是對圖4(a)經(jīng)本文算法去噪后的圖像。由圖4可以觀察到：中值非線性和均值線性去噪方法都能在一定程度上抑制椒鹽噪聲，中值法比均值法較好地保持了醫(yī)學(xué)影像圖像的原有細(xì)節(jié)，但所得醫(yī)學(xué)影像圖像較模糊、清晰度欠佳。本文提出的方法對椒鹽噪聲的濾除效果最佳，能很好地保護(hù)醫(yī)學(xué)影像圖像的細(xì)節(jié)信息，使圖像清晰度更高。

為了進(jìn)一步驗證本文提出的算法在胰腺超聲圖像中的去噪聲效果，用均值、中值和本文提出的算法等3種方法將圖像進(jìn)行去噪聲運(yùn)算，通過峰值信噪比(peak signal to noise ratio,PSNR)衡量3種方法的去噪功能，各方法的PSNR指標(biāo)見表1。由表1可見，本文算法的峰值信噪比顯著大于其他兩種算法，表明該方法去噪性能優(yōu)于其他兩種方法。由以上可知，本文算法不僅能有效去除胰腺超聲圖像中含有高密度的椒鹽和高斯噪聲，又可有效保持圖像的細(xì)節(jié)。

3 結(jié)語

研究表明，運(yùn)用閾值鄰域均值濾波法既能夠有效抑制胰腺超聲圖像中含有的高密度的椒鹽和高斯噪聲，又可有效保持圖像的細(xì)節(jié)信息。因此，本方法對醫(yī)學(xué)影像圖像的去噪研究有良好的參考價值。在今后的研究中，筆者將對醫(yī)學(xué)影像圖像中的其他類型噪聲的去除進(jìn)行研究。