王 強，張合新，孟 飛，熊 鵬

?

基于模糊系統(tǒng)的散斑噪聲濾波器

王 強，張合新，孟 飛，熊 鵬

（第二炮兵工程大學(xué)控制工程系，陜西 西安， 710025）

為了能有效濾除散斑噪聲，提出了一種以模糊系統(tǒng)為基礎(chǔ)的散斑噪聲濾波器。這種濾波器包括了一個模糊推理系統(tǒng)，一個邊緣檢測和膨脹模塊和一個圖像合成器。模糊推理系統(tǒng)包括5個輸入和1個輸出，負(fù)責(zé)對散斑噪聲圖像進(jìn)行濾波處理，其參數(shù)通過克隆選擇優(yōu)化算法訓(xùn)練優(yōu)化得到。邊緣檢測和膨脹模塊用于區(qū)分圖像的邊緣區(qū)域和光滑區(qū)域，圖像合成器則是根據(jù)邊緣檢測和膨脹模塊獲取的信息分區(qū)域?qū)V波輸出圖像進(jìn)行合成。文中將該方法與其他幾種常用的散斑噪聲濾除方法進(jìn)行了比較，實驗結(jié)果表明，與其他方法相比該方法可以顯著減少圖像的散斑噪聲，同時保留邊緣、細(xì)節(jié)等有價值的信息。

散斑噪聲濾波器；模糊系統(tǒng)；圖像處理

0 引言

成像系統(tǒng)在獲取圖像過程中會受各種干擾因素的影響，噪聲則是最常見的干擾因素。在信號和圖像處理領(lǐng)域有很多噪聲模型以及不同噪聲對于圖像質(zhì)量的影響分析。其中有一種是散斑噪聲模型，這是一種乘性噪聲。在這種模型中，圖像質(zhì)量被降低是因為圖像中得到的像素值被乘以了噪聲分量。合成孔徑雷達(dá)，超聲波、激光等成像系統(tǒng)都會受到散斑噪聲的影響。在這些系統(tǒng)中，散斑噪聲的產(chǎn)生是物體表面反射信號干涉產(chǎn)生的結(jié)果。在最后得到的圖像中幾乎所有的像素都會受到該噪聲的影響，最終導(dǎo)致的結(jié)果便是圖像中產(chǎn)生整塊的斑紋以致于降低了圖像的信噪比和分辨率。一些常用的散斑噪聲濾波器的目的在于濾除噪聲的同時不丟失有用的信息，保留邊緣等的有價值的細(xì)節(jié)。

在此之前，已經(jīng)有很多用于散斑噪聲濾除的方法已經(jīng)被提出。其中基于圖像的局部統(tǒng)計濾波的Lee方法，基于同態(tài)映射的同態(tài)濾波，小波軟閥值算法和基于數(shù)學(xué)形態(tài)的多方向形態(tài)濾波及其它們的改進(jìn)型是最常見的方法[1]。在噪聲圖像中，細(xì)節(jié)或邊緣與噪聲之間的區(qū)別并不明顯。這種信息的不確定性和不完整性給濾波操作帶來了困難，而模糊理論正好能適應(yīng)這些非確定性，因此在圖像濾波過程中應(yīng)用模糊理論可以取得較好的噪聲去除效果[2]。近年來一些以模糊邏輯為基礎(chǔ)的方法被提出用于噪聲濾除。例如：F. Russo等提出了FF方法，將模糊系統(tǒng)理論應(yīng)用于圖像濾波。羅海馳等提出了一種基于自適應(yīng)神經(jīng)模糊推理系統(tǒng)[3]用于濾除椒鹽脈沖噪聲。而本文提出了一種以模糊推理系統(tǒng)為基礎(chǔ)的散斑噪聲濾波器。實驗結(jié)果表明，本文所提出的濾波器能有效濾除散斑噪聲，而且能夠很好保護(hù)邊緣、細(xì)節(jié)等信息。

1 散斑噪聲基本性質(zhì)及常用濾除方法

1.1 散斑噪聲基本性質(zhì)

為了能夠更好地去除散斑噪聲，首先對散斑噪聲的性質(zhì)進(jìn)行分析研究。散斑噪聲可以近似為一種乘性噪聲，具有負(fù)指數(shù)分布特性[4]?？梢杂萌缦碌脑肼暷Ｐ蛠斫票硎荆?/p>

＝0×sn(1)

式中：為接收到的信號強度，通過圖像中的像素灰度值來體現(xiàn)；0是理想無噪聲信號強度，它反映的是目標(biāo)的真實信息；sn是散斑作用系數(shù)；á?表示求平均值。

1.2 常用的散斑噪聲去除方法

1.2.1 同態(tài)濾波

在前文中提到，散斑噪聲是一種乘性的非線性噪聲，乘性噪聲存在的一個問題就是它與信號的相倚性比較強，而同態(tài)變換可以使被乘性噪聲污染的信號實現(xiàn)信噪分離，使乘性噪聲轉(zhuǎn)變成為加性噪聲，削弱噪聲與信號的相倚性[5]。同態(tài)濾波方法就是首先對信號進(jìn)行同態(tài)變換，然后用一種合適的線性濾波器抑制噪聲，最后進(jìn)行逆變換，恢復(fù)信號。由于在進(jìn)行濾波的時候采用的是一種線性濾波器，勢必會造成圖像邊緣和細(xì)節(jié)上的模糊和缺失。

1.2.2 局部統(tǒng)計濾波

局部統(tǒng)計濾波算法由Lee J. S.提出所以這算法又叫Lee濾波器[6]。Lee濾波器假定像素的采樣均值和變差等于其鄰域的局部均值和變差，并且將圖像的強度函數(shù)用泰勒展開，取其一階使其線性化，然后根據(jù)最小平方誤差估計在被散斑噪聲污染信號中獲得逼近信號。它的實質(zhì)就是用一個線性函數(shù)來逼近一個非線性函數(shù)，用這種方法處理后肯定會造成圖像邊緣的模糊。

Lee濾波算法可以用下式來表示：

式中：(,)是濾波器逼近的結(jié)果；(,)是含噪圖像信號。

1.2.3 小波閾值算法

小波變換是繼傅里葉分析之后另一種有效信號分析工具，它的實質(zhì)在于減小以至于完全剔除由噪聲產(chǎn)生的系數(shù)，同時最大限度地保留有效信號產(chǎn)生的系數(shù)，再由處理后的小波系數(shù)重構(gòu)原信號，獲得原信號的最佳估計。小波閾值算法主要分為以下3步[7]：

①應(yīng)用Mallat算法對信號進(jìn)行分解：選擇合適的基小波，并確定尺度，然后進(jìn)行小波分解。

②閾值法修正小波系數(shù)：對每一尺度選擇適當(dāng)?shù)拈撝担匆欢ǖ牟呗?，對其尺度Mallat算法對信號的分解的系數(shù)進(jìn)行

③信號重建：通過小波逆變換，利用原信號尺度的近似和經(jīng)閾值處理后的小波系數(shù)重建原信號。

2 模糊推理系統(tǒng)

在本文所提出的方法中使用的模糊推理系統(tǒng)是一個5輸入單輸出的一階Sugeno模糊模型，對于它的輸入選擇了廣義鐘型隸屬函數(shù)，輸出則選用了線性函數(shù)。這個模糊推理系統(tǒng)包含了10條規(guī)則，其模糊規(guī)則集如下：

① if (1?11)&(2?12)&(3?13)&(4?14)

&(5?15)

then1＝111＋122＋133＋144＋155＋16

② if (1?21)&(2?22)&(3?23)&(4?24)

&(5?25)

then1＝211＋222＋233＋244＋255＋26

③ if (1?31)&(2?32)&(3?33)&(4?34)

&(5?35)

then1＝311＋322＋333＋344＋355＋36

……

⑩ if (1?101)&(2?102)&(3?103)&(4?104) &(5?105) then1＝1011＋1022＋1033＋1044＋1055＋106(7)

式中：X是模糊系統(tǒng)的輸入；Q表示系統(tǒng)根據(jù)第條規(guī)則得到的輸出；d,是結(jié)論參數(shù)；M,表示第個輸入的第個隸屬函數(shù)。對于輸入X，定義的廣義鐘型隸屬函數(shù)描述如下：

式中：a,，b,，c,為前提參數(shù)，它們和d,一起確定隸屬函數(shù)的形狀，這些參數(shù)在訓(xùn)練過程中不斷優(yōu)化，關(guān)于訓(xùn)練過程這將在后面進(jìn)行說明。

模糊推理系統(tǒng)的輸出等于各Q的加權(quán)平均：

式中：加權(quán)系數(shù)表示第條規(guī)則的激勵強度，的計算方法如下：

3 克隆選擇優(yōu)化算法

克隆選擇優(yōu)化算法（Clonal selection optima- zation algorithm，CSOA）的思想靈感來自于人類免疫系統(tǒng)的克隆選擇原理，它能夠在尋求全局最優(yōu)解的過程中避免陷入局部最優(yōu)的求解，在求解具有多模態(tài)性質(zhì)或時變性質(zhì)的優(yōu)化問題中能夠很好的體現(xiàn)CSOA的這一優(yōu)勢[8]。因此，CSOA正被廣泛的應(yīng)用于工程實踐中。CSOA的原理程序框圖如圖1所示。

初始種群隨機(jī)產(chǎn)生，待求解的問題由抗原表示，可行的解由抗體表示。在算法中根據(jù)抗體和抗原的親和度對抗體進(jìn)行排序，然后根據(jù)親和度不同按不同比例進(jìn)行克隆?？寺?shù)量的計算方法如下：

C＝round(/) (11)

式中：是克隆系數(shù)；是可行解的排序號。從式(1)中可以看出克隆的數(shù)量與親和度有關(guān)，親和度越高復(fù)制數(shù)量越多，反之則越少。于此相反的是：當(dāng)抗體的親和度較高時抗體的變異就小，當(dāng)抗體的親和度較小時，抗體的變異就大。最終這些抗體將被重新評估與抗原的親和度，最合適的抗體將出現(xiàn)在下一代群體中。這個循環(huán)直到達(dá)到終止條件時才會停止。

圖1 CSOA原理程序框圖

4 模糊濾波器

前文中所提到的模糊推理系統(tǒng)通常作為像小波變換這種以模糊邏輯為基礎(chǔ)方法的一部分出現(xiàn)。在這些方法中模糊推理系統(tǒng)通常被用于系數(shù)的調(diào)整與優(yōu)化。但是單獨一個模糊推理系統(tǒng)也可以用于數(shù)字圖像的噪聲濾波只要為它提供訓(xùn)練所需的圖像。在模糊推理系統(tǒng)參數(shù)被準(zhǔn)確找到的前提下，一旦模糊推理系統(tǒng)通過訓(xùn)練，只要一個簡單優(yōu)化算法就可以使系統(tǒng)輸出的圖像與沒有噪聲影響的圖像幾乎沒有什么差別。然而模糊推理系統(tǒng)的參數(shù)確定十分困難，在許多情況中，在有一些微小誤差的情況下也是能得到很不錯的去噪效果的。

4.1 模糊推理系統(tǒng)的訓(xùn)練

模糊推理系統(tǒng)輸入與輸出之間的關(guān)系由用于輸入的隸屬函數(shù)決定。在系統(tǒng)的訓(xùn)練階段，通過所提供的訓(xùn)練數(shù)據(jù)可以訓(xùn)練得到?jīng)Q定模糊規(guī)則的隸屬函數(shù)的系數(shù)。本文模糊濾波器中的模糊系統(tǒng)的參數(shù)則通過CSOA進(jìn)行優(yōu)化。訓(xùn)練的步驟如圖2所示，當(dāng)訓(xùn)練完成時輸出的圖像無限趨近于原始的無噪圖像。

本文所提出模糊濾波器使用的模糊推理系統(tǒng)有5個輸入和1個輸出。模糊規(guī)則庫包含了10條規(guī)則，對于模糊系統(tǒng)規(guī)則庫的數(shù)量可以進(jìn)行試探性的實驗，它可以是大于或小于10。實驗結(jié)果表明，增加規(guī)則的數(shù)量可以略微提高所提出的模糊濾波器的性能，但同時也增加了訓(xùn)練的時間以及該系統(tǒng)的復(fù)雜性。因此，在決定規(guī)則的數(shù)量時應(yīng)盡量平衡系統(tǒng)性能和系統(tǒng)的復(fù)雜性之間的關(guān)系。

圖2 模糊系統(tǒng)訓(xùn)練

當(dāng)前以模糊濾波器為基礎(chǔ)的方法已經(jīng)成功地應(yīng)用于數(shù)字圖像的脈沖噪聲的去除。脈沖噪聲通常只會將整幅圖像中的某些像素的像素值變?yōu)?或255，但是散斑噪聲卻能影響圖像中某個區(qū)域內(nèi)的所有像素。因此，合適訓(xùn)練圖像對于構(gòu)造一個用于散斑噪聲濾除的模糊濾波器具有非常重要的作用。訓(xùn)練圖像必須同時包含含有不同灰度值的光滑的區(qū)域和邊緣區(qū)域。因此，訓(xùn)練圖像可以通過在各種人造的無噪聲圖像上添加散斑噪聲來獲得，當(dāng)模糊系統(tǒng)通過這些圖像訓(xùn)練后，它就可以用于各種散斑圖像的濾波。

圖3(a)和3(b)所示的圖像就是在計算機(jī)中生成的用于訓(xùn)練的圖像。圖3(a)是一幅無噪訓(xùn)練圖，尺寸是120×120像素，其中的每個小方塊大小是8×8像素，每個小方塊中的像素值都相同隨機(jī)均與分布于0和255之間。含噪訓(xùn)練圖3(b)則是在無噪圖像上添加了方差為0.04的散斑噪聲后得到的。噪聲方差值對于模糊系統(tǒng)的訓(xùn)練不是十分關(guān)鍵，但是根據(jù)以往的實驗經(jīng)驗可知，當(dāng)輸入的含噪圖像的噪聲方差接近于訓(xùn)練圖像的噪聲方差時才能得到最優(yōu)的濾波結(jié)果。在通常情況下不可能知道要處理圖像的噪聲方差，所以訓(xùn)練圖像的方差大小的選取全憑經(jīng)驗，不能選太大或太小[9-11]。

圖3 訓(xùn)練圖

在使用訓(xùn)練圖像對系統(tǒng)進(jìn)行訓(xùn)練的時候選擇了一種特定的鄰域。這種特定的像素鄰域在濾波窗口中的形式如圖4(a)中所示，它的應(yīng)用方法則表示在圖4(b)中。其中TN作為模糊系統(tǒng)的輸入從含噪訓(xùn)練圖像中獲取，系統(tǒng)輸出用TO表示。在鄰域的選擇上不止本文中所示的這一種方法，但是實驗結(jié)果表明：使用本文的鄰域選取方法是本文濾波器能達(dá)到最好的效果。

圖4 濾波窗口及其應(yīng)用

4.2 本文濾波器

作為一種噪聲濾波器，在使用的過程中最容易出現(xiàn)的問題就是濾波造成圖像的邊緣模糊。在濾波處理中光滑區(qū)域受的影響往往較小，然而在邊緣和細(xì)節(jié)處有可能會產(chǎn)生信息丟失或者邊緣畸變，最終導(dǎo)致的結(jié)果就是圖像的理解性變差和解釋性變得復(fù)雜[12-13]。如果要想獲得好的濾波效果，那么濾波器就必須對邊緣和細(xì)節(jié)信息十分敏感。

本文所提出的濾波器的結(jié)構(gòu)如圖5所示。在通常情況下，通過訓(xùn)練的模糊濾波器可以滿足對散斑噪聲的濾波需求，但是為了提高濾波器的性能，模糊系統(tǒng)采用了一種級聯(lián)方式。這樣處理的結(jié)果是：光滑區(qū)域更加平滑，但是與此同時細(xì)節(jié)和邊緣區(qū)域也被平滑處理了。為了避免邊緣平滑的發(fā)生，在濾波器中加了一個邊緣檢測和膨脹模塊，該模塊中獲取到的信息將用于最后的圖像復(fù)原。

4.2.1 邊緣檢測與膨脹模塊

圖5 本文模糊濾波器

在數(shù)字圖像中通過邊緣算子可以很容易地檢測出一幅圖像的邊緣區(qū)域，但是在一幅被噪聲污染的圖像中就不是那么容易了，因為噪聲會在圖像中產(chǎn)生一些偽邊緣[14-15]。在本文濾波器中，在圖像進(jìn)入邊緣檢測與膨脹模塊之前已經(jīng)濾除了散斑噪聲的影響，所以只要一些簡單的邊緣檢測算子就可以很容易地檢測出圖像的邊緣。在這里選用Sobel算子作為邊緣檢測算子。

圖像的膨脹處理是為了檢測出圖像的邊緣附近的區(qū)域，而且在本文方法中要對邊緣圖像要進(jìn)行兩次的膨脹處理，并在邊緣附近標(biāo)記出3×3像素大小的膨脹區(qū)。經(jīng)過這樣處理后圖像整體被分成了3個部分：邊緣區(qū)域、邊緣過渡區(qū)、光滑區(qū)。

4.2.2 圖像合成器

圖像的邊緣及其膨脹處理后的信息將用于圖像的復(fù)原。圖像的邊緣區(qū)域的像素值由第一次濾波輸出圖來確定，而邊緣過渡區(qū)的像素值則用兩次濾波后的輸出圖來確定。在膨脹區(qū)之外的光滑區(qū)像素值則通過計算兩張濾波處理后的輸出圖的均值來確定。從像素值的計算方法上來看，一次濾波輸出的對第一次膨脹區(qū)域貢獻(xiàn)更大，兩次濾波的對輸出對第二次膨脹區(qū)域貢獻(xiàn)更大。實驗結(jié)果表明通過這樣處理模糊濾波器的性能能有很大提升。

5 實驗結(jié)果及分析

為了驗證本文中所提出的模糊濾波器的性能，采用一些常用的、人為加入散斑噪聲的標(biāo)準(zhǔn)數(shù)字圖像作為測試對象。在測試中選用了3張不同的圖像（Lena，Bridge，Boats），每張圖像中又加入了均值為零但是方差不同的3種散斑噪聲。

為了驗證該濾波器的性能將它與現(xiàn)在國內(nèi)外使用效果較好的濾波方法進(jìn)行對比，包括了最小均方誤差Lee濾波器（Minimum mean squared error Lee filter，Lee MMSEF），F(xiàn)rost濾波器（Frost filter，F(xiàn)ROSTF），ADF濾波器（Anisotropic diffusion filter，ADF），SRADF濾波器（Anisotropic diffusion speckle noise removal filter，SRADF），圖像平滑Lee濾波器（Image smoothing Lee filter，Lee ISF），中值濾波器（Median filter，MEDF）和均值濾波器（Mean filter，MF）。而比較的標(biāo)準(zhǔn)主要是噪聲的濾除效果，和濾波器對邊緣及細(xì)節(jié)處的保護(hù)能力。濾波器對噪聲的濾除效果肉眼可以很明顯得看出來。至于邊緣和細(xì)節(jié)保護(hù)能力選用均方誤差（Mean Squared Error，MSE）作為評價參數(shù)[16]。

圖6中(a)～(j)分別是標(biāo)準(zhǔn)圖、噪聲方差為0.04的未經(jīng)處理的圖、均值濾波圖、中值濾波圖、Lee MMSEF濾波圖、FROSTF濾波圖、ADF濾波圖、SRADF濾波圖、Lee ISF濾波圖和本文方法濾波圖。首先我們從散斑噪聲濾除效果上將本文方法與其他方法來進(jìn)行對比分析。從以上幾張圖中可以看到圖(c)、(d)、(f)、(g)中還明顯存在散斑噪聲，說明均值濾波、中值濾波、FROSTF濾波、ADF濾波這幾種方法對散斑噪聲的濾除能力較弱，從功能上來說就不能滿足實際應(yīng)用的需要。在考慮了濾波效果之后要對剩下的幾種濾波器在邊緣和細(xì)節(jié)的保護(hù)能力上來進(jìn)行對比分析。其中，圖(h)的邊緣和細(xì)節(jié)非常模糊，而且在處理完之后整幅圖像顏色偏暗。圖(i)的整體效果還可以，但是與圖(j)相比，在人物的面部細(xì)節(jié)上與帽子的輪廓上它的模糊感覺更加明顯。相比之下圖(e)和圖(j)的效果最好，也就是說在處理方差為0.04的散斑噪聲時Lee MMSEF濾波器和本文提出的濾波器效果最好。

為了更加直接地說明各個濾波器對圖像邊緣及細(xì)節(jié)的保護(hù)能力，同時也為了試驗這些濾波器對圖像噪聲是不是具有選擇性，本文對選用的3張標(biāo)準(zhǔn)圖片分別加入了標(biāo)準(zhǔn)差為0.01，0.04和0.07的散斑噪聲，然后計算每種濾波器對這9張圖處理后的MSE值，最后利用這些數(shù)據(jù)分析各個濾波器的性能。結(jié)果如表1和表2所示。

從表1和表2中可以看出本文提出的濾波器在邊緣和細(xì)節(jié)的保護(hù)上相比于其他幾種濾波器做得更好。在上文中提到當(dāng)散斑噪聲方差為0.04的時Lee MMSEF濾波器和本文提出的濾波器效果都比較好，但是從表中數(shù)據(jù)可以看到當(dāng)方差為0.01和0.07時Lee MMSEF濾波器的性能就明顯不如本文算法了，這說明了本文提出的濾波器對噪聲沒有選擇性。綜上所述，從對噪聲濾除的效果和邊緣、細(xì)節(jié)信息保護(hù)這兩方面來看本文提出的濾波器的性能更加全面、穩(wěn)定。

6 結(jié)束語

本文提出了一種以模糊系統(tǒng)為基礎(chǔ)的散斑噪聲濾波器。該濾波器首先通過模糊推理系統(tǒng)實現(xiàn)對含噪圖像的去噪處理，然后用邊緣檢測和膨脹模塊區(qū)分圖像的邊緣和平滑區(qū)，最后根據(jù)區(qū)分的不同區(qū)域在圖像合成器中采用不同的方法實現(xiàn)圖像合成、復(fù)原。其中的模糊推理系統(tǒng)的參數(shù)通過人工構(gòu)造的計算機(jī)圖像進(jìn)行訓(xùn)練，克隆選擇算法優(yōu)化得到。由實驗結(jié)果可知，與其他濾波器相比該濾波器在有效濾除散斑噪聲的同時可以很好地保護(hù)圖像的邊緣、細(xì)節(jié)等重要敏感信息。

表1 Lena、Bridge和 Boats濾波處理后的MSE值

表2 濾波輸出圖的MSE均值

[1] 蔣立輝, 趙春暉, 王騏. 散斑噪聲抑制算法比較研究[J]. 激光與紅外, 2002, 32(6): 431-434.

JIANG Lihui, ZHAO Chunhui, WANG Qi . Study and comparison of speckle noise suppressing algorithms[J]., 2002, 32(6): 431-434.

[2] 朱立新, 楊擴(kuò), 秦加合. 一種新的基于自適應(yīng)神經(jīng)網(wǎng)絡(luò)模糊推理系統(tǒng)的圖像濾波器[J]. 計算機(jī)科學(xué), 2014, 41(6): 211-214.

ZHU Lixin, YANG Kuo, QIN Jiahe. New noise digital images filter based on adaptive neuro-fuzzy nnference System[J]., 2014, 41(6): 211-214.

[3] 羅海馳, 李岳陽, 孫俊. 一種基于自適應(yīng)神經(jīng)模糊推理系統(tǒng)的圖像濾波方法[J]. 計算機(jī)科學(xué), 2013, 40(7): 302-306.

LUO Haichi, LI Yueyang, SUN Jun. Filtering method for images based on adaptive neuro-fuzzy inference system[J]., 2013, 40(7): 302-306.

[4] 李曉峰, 徐軍, 羅積軍, 等. 激光主動成像圖像噪聲分析與抑制[J]. 紅外與激光工程, 2011, 40(2): 332-337.

LI Xiaofeng, XU Jun, LUO Jijun, et al. Noise analyzing and denoising of intensity image for laser active imaging system[J]., 2011, 40(2): 332-337.

[5] LEE J S, WEN J H, AINSWORTH T L, et al. Improved sigma ?lter for speckle ?ltering of SAR imagery[J]., 2009, 47(4): 202-213.

[6] PUVANATHASAN P, BIZHEVA K. Interval type-II fuzzy anisotropic diffusion algorithm for speckle noise reduction inoptical coherence tomography[J]., 2009, 17(5): 733-746.

[7] ZHANG Y, CHENG H D, TIAN J, et al. Fractional subpixel diffusion and fuzzy logic approach for ultrasound speckle reduction[J]., 2010, 43(2): 2962-2970.

[8] 唐正, 胡珉. 空間自適應(yīng)免疫克隆選擇優(yōu)化算法[J]. 計算機(jī)應(yīng)用, 2009, 29(2): 561-564.

TANG Zheng, HU Min. Space self-adaptive inmune clonal selection optimization algorithm[J]., 2009, 29(2): 561-564.

[9] YUKSEL M E, BASTURK A.Application of type-2 fuzzy logic ?ltering to reduce noise in color images[J].2012, 56(7): 25-35.

[10] AMIRMAZLAGHANI M, AMINDAVAR H. A novel sparse method for despeckling SAR images[J].2012,50(4):5024-5032.

[11] HAJIABOLI M R. An anisotropic fourth-order diffusion ?lter for image noise removal[J].2011,92(8): 177-191.

[12] 王燦進(jìn), 孫濤, 陳娟. 基于像素點分類的激光主動成像混合濾波[J]. 中國激光, 2014, 41(3): 1-7.

WANG Canjin, SUN Tao, CHEN Juan. Hybrid filter based on pixel classification for laser active imaging[J]., 2014, 41(3): 1-7.

[13] ELAMVAZUTHI I, ZAIN M L B M, BEGAM K M. Despeckling of ultrasound images of bone fracture using multiple?ltering algorithms[J].2013,57(6): 152-168.

[14] HUANG B, BU P, WANG X, et al. Speckle reduction in parallel optical coherence tomography by spatialcompounding[J].2013,45(2): 69-73.

[15] ORTIZ S H C, CHIU T,FOX M D. Ultrasound image enhancement: a review2012, 25(3): 419-428.

[16] 韓紅光. NSCT域紅外圖像改進(jìn)非局部均值濾波算法[J]. 紅外技術(shù), 2015, 37(1) :34-38.

HAN Hong-guang. Improved non-local means infrared image filtering algorithm based on NSCT[J]., 2015, 37(2): 34-38.

A Novel Fuzzy Filter for Speckle Noise Removal

WANG Qiang，ZHANG Hexin，MENG Fei，XIONG Peng

（Department of Automation, The Second Artillery Engineering University, Xi￠an 710025, China）

In this paper, a novel fuzzy system-based filter for speckle noise removal is proposed. The proposed filter consists of a fuzzy inference system, an edge detection and dilation unit, and an image combiner. The fuzzy inference system includes 5 inputs and 1 output, and it is responsible for filtering the speckle noisy image. The edge detection and dilation unit is used for classifying the uniform areas and edge areas. The image combiner unites the output images according to the information coming from the edge detection and dilation unit. The training phase of the fuzzy inference system is implemented by using the clonal selection optimization algorithm with appropriate training data. The performance of the proposed method is compared with popular speckle noise removal filters available in the literature by performing extensive simulations. The experimental results show that the proposed method can significantly reduce the speckle noise from digital images while preserving edges, and valuable details.

speckle noise filtering，fuzzy inference system，image processing

TN751

A

1001-8891(2016)05-0415-07

2015-12-28；

2016-02-29.

王強（1993-），男，碩士研究生，主要從事激光主動成像及目標(biāo)識別方向研究。E-mail：15158001096@163.com。

國家自然科學(xué)基金（61203189）。