張魯松，張壽明

（昆明理工大學(xué) 信息工程與自動(dòng)化學(xué)院，云南 昆明 650500）

0 引言

我國(guó)是電力生產(chǎn)大國(guó)，在熱電廠鍋爐運(yùn)行中，尤其是在鍋爐的點(diǎn)火與滅火環(huán)節(jié)，極易發(fā)生由于火焰檢測(cè)不當(dāng)、煤粉與油供給量不平衡而導(dǎo)致的滅火與“放炮”等惡性事故，輕則導(dǎo)致機(jī)組停機(jī)影響發(fā)電質(zhì)量，重則會(huì)導(dǎo)致?tīng)t膛爆炸等危及人身安全的事故[1]。為了保證發(fā)電環(huán)節(jié)能安全高效地進(jìn)行，大型熱電廠均配置爐膛安全系統(tǒng)（Furnace Safeguard Supervisory System，簡(jiǎn)稱(chēng)FSSS）。其替代了原始的窺火孔人工觀察火焰的模式，做到可以實(shí)時(shí)清晰地觀察爐膛內(nèi)煤粉與燃料的燃燒情況。

基于火焰電視的爐膛火焰檢測(cè)系統(tǒng)，使?fàn)t膛監(jiān)視工作更為安全，減少了人工定期巡查作業(yè)，減少了安全隱患，并且相比較傳統(tǒng)人工監(jiān)視方式，通過(guò)火焰電視監(jiān)視的方式更加直觀全面。但爐膛周?chē)墓ぷ鳝h(huán)境以及爐膛內(nèi)煤粉與燃油的復(fù)雜燃燒會(huì)使基于 CCD工業(yè)相機(jī)的爐膛火焰電視系統(tǒng)在采集圖像與視頻的同時(shí)產(chǎn)生大量的圖像噪聲，震動(dòng)與煤粉顆粒以及燃燒結(jié)焦產(chǎn)生的圖像缺失、圖像模糊以及煤粉燃燒時(shí)閃爍產(chǎn)生的脈沖噪聲會(huì)大大影響爐膛火焰的監(jiān)視工作[2-3]。

爐膛火焰圖像相比較傳統(tǒng)火焰圖像有高溫區(qū)與低溫區(qū)邊界不明顯、高溫區(qū)“火云”形狀變化快且隨機(jī)性強(qiáng)的特點(diǎn)，為使計(jì)算機(jī)可以對(duì)圖片細(xì)節(jié)完整、圖片內(nèi)容豐富的火焰圖像進(jìn)行分析，因此需要對(duì)火焰電視采集的原始火焰圖像進(jìn)行噪聲濾除處理，為使計(jì)算機(jī)可以從噪聲圖像中獲得更多的圖片信息，國(guó)內(nèi)外對(duì)噪聲濾波[4-6]與噪聲圖像的邊緣提取進(jìn)行了廣泛研究[7-10]。本文通過(guò)matlab仿真實(shí)現(xiàn)ABDND算法并對(duì)典型煤粉燃燒灰度圖像進(jìn)行濾波實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果表明濾波后圖像細(xì)節(jié)保留完整，噪聲濾波效果良好。

1 經(jīng)典噪聲濾波算法

1.1 多幀平均法

設(shè)由攝像頭采集的圖像為I（x,y），I（x,y）是由原圖像 f （ x, y）與噪聲 n （ x, y ）疊加形成，假設(shè) n （ x, y ）內(nèi)元素都不相關(guān)且均值為零，則采集圖像 I（ x, y ）可表示為：

在Δt時(shí)間內(nèi)采集M幀圖像，對(duì)M幀圖像求平均值，可得到Δt時(shí)間內(nèi)的平均圖像 I（ x, y ）可表示為：

其中 Ij（x,y） 為第j幀采集的圖像， WJ表示權(quán)值參數(shù)。

1.2 鄰域平均法

鄰域平均法是一種線(xiàn)性濾波算法，鄰域平均法又被稱(chēng)為均值濾波。算法實(shí)現(xiàn)方法是使數(shù)字圖像內(nèi)像素點(diǎn)的灰度值等于該點(diǎn)鄰域內(nèi)像素點(diǎn)的平均灰度值。設(shè) f（x, y）是一幅m×n的數(shù)字圖像的灰度圖像，經(jīng)鄰域平均法處理后的圖像為g（x, y ），則g（x, y ）可表示為：

其中，S為像素點(diǎn)（x, y ）為中心鄰域內(nèi)像素點(diǎn)的集合，M表示集合S內(nèi)像素點(diǎn)的個(gè)數(shù)。

1.3 中值濾波法

中值濾波是基于排序統(tǒng)計(jì)理論的一種可以有效抑制噪聲的非線(xiàn)性濾波算法，以灰度圖像為例，其基本原理是將數(shù)字圖像或數(shù)字信號(hào)中一點(diǎn)的灰度值用該點(diǎn)領(lǐng)域中各點(diǎn)灰度值的中值代替，使周?chē)袼氐幕叶戎到咏鎸?shí)值，從而消除孤立的噪聲點(diǎn)。方法是用某種結(jié)構(gòu)的二維滑動(dòng)模板（通常為ww×的方形區(qū)域），將模板內(nèi)的像素灰度值大小進(jìn)行升序或降序排列，生成二維單調(diào)上升（或下降）的二維數(shù)據(jù)序列，設(shè)（）g.xy為輸出后圖像，則可表示為：

其中f （ x, y）為原圖像，w為二維模板的邊界值。

1.4 自適應(yīng)中值濾波法

自適應(yīng)中值濾波算法在中值濾波法的基礎(chǔ)上增加了濾波窗口大小調(diào)節(jié)算法，使其可以先檢測(cè)噪聲像素再對(duì)噪聲進(jìn)行濾波。設(shè) WT為當(dāng)前中值濾波窗口的邊緣值， fxy為像素點(diǎn)（x, y ）處的灰度值， fmax為濾波窗口 WT內(nèi)所有像素點(diǎn)的灰度值的最大值， fmin為濾波窗口 WT內(nèi)所有像素點(diǎn)的灰度值的最小值，fmed為濾波窗口 WT內(nèi)所有像素點(diǎn)的灰度值的平均值， Wmax為濾波窗口 WT的邊緣最大值， Wmin為濾波窗口 WT的邊緣最小值，具體實(shí)現(xiàn)算法如下：

當(dāng)M1＞0、M2＜0時(shí)，算法步驟由步驟（1）轉(zhuǎn)為步驟（2），否則，濾波窗口 WT的邊緣值增大，若滿(mǎn)足 Wmin≤ WT≤ Wmax，則重復(fù)執(zhí)行步驟（1），否則窗口輸出灰度值為fxy，當(dāng)N1＞0、N2＜0時(shí)，則輸出窗口灰度值為 fxy，負(fù)責(zé)輸出灰度值 fmed[11-12]。

綜上所述經(jīng)典濾波算法中，多幀平均法在實(shí)際應(yīng)用中會(huì)因?yàn)閹瑪?shù)增多而影響處理速度，并且煤粉在爐膛內(nèi)燃燒時(shí)煤粉燃燒火焰的變化頻率高、隨機(jī)性高等特點(diǎn)使多幀平均法并不適用于爐膛火焰脈沖噪聲濾除。鄰域平均法雖然算法簡(jiǎn)單且處理速度快，但其不足是經(jīng)鄰域平均法處理后的圖像邊緣會(huì)不完整且變得模糊，圖像細(xì)節(jié)不能很好地保留。中值濾波法是應(yīng)用較為普遍的脈沖濾波算法，自適應(yīng)中值濾波算法在中值濾波算法的基礎(chǔ)上做出改進(jìn)，濾波效果更佳，本文以自適應(yīng)濾波算法與ABDND算法進(jìn)行濾波仿真實(shí)驗(yàn)并對(duì)濾波結(jié)果進(jìn)行對(duì)比。

2 ABDND開(kāi)關(guān)中值濾波算法

ABDND（advanced boundary discriminative noise detection algorithm）是由A.K. Tripathi等人提出的一種基于自適應(yīng)開(kāi)關(guān)中值濾波算法（ASMF）[13]與邊緣鑒別噪聲檢測(cè)算法（BAND）[14]的改進(jìn)算法，其不僅可以準(zhǔn)確判斷噪聲像素，還可以準(zhǔn)確檢測(cè)隨機(jī)脈沖噪聲。ABDND算法主要由三個(gè)噪聲處理階段組成，階段一為噪聲像素檢測(cè)階段；階段二為噪聲像素核實(shí)階段；階段三為噪聲像素濾除階段[15]。其具體實(shí)現(xiàn)過(guò)程如下：

（1）噪聲像素檢測(cè)階段

1. 獲取圖像所有像素點(diǎn)的灰度直方圖，并對(duì)每個(gè)像素點(diǎn)相鄰像素點(diǎn)的灰度值直方圖數(shù)值做減法處理處理，得到差值序列ΔiD :

2. 根據(jù)差值序列ΔiD，得到差值中的正最大值和負(fù)最大值，根據(jù)公式（8）確定判斷噪聲像素的兩個(gè)邊界1T和2T：

3. 采用邊界尺寸大小為w×w的濾波窗口對(duì)圖像像素點(diǎn)進(jìn)行遍歷（本文中設(shè)定W=11），濾波窗口的中心點(diǎn)覆蓋于待檢測(cè)像素點(diǎn)，并統(tǒng)計(jì)該窗口內(nèi)像素點(diǎn)的最大灰度值 Smax和最小灰度值 Smin，根據(jù)式（9）得到濾波圖像的噪聲標(biāo)識(shí)矩陣 r（i, j），初步判斷該像素點(diǎn)是否是噪聲像素。

其中， si,j是圖像中點(diǎn)（i, j）的像素灰度值，r（i,j）=1表示該點(diǎn)被初步判定為噪聲像素點(diǎn)，r（i, j ） =0表示該點(diǎn)不是噪聲像素點(diǎn)。

（2）噪聲像素核實(shí)階段

對(duì)于臨時(shí)噪聲像素作以下算法處理：

1. 以臨時(shí)噪聲像素點(diǎn)為中心，取邊界尺寸大小為w×w的窗口，設(shè)w遠(yuǎn)小于W（本文中設(shè)w=3），隨后統(tǒng)計(jì)窗口中非噪聲像素點(diǎn)的數(shù)量 df，若df＞w，計(jì)算窗口內(nèi)所有非噪聲像素點(diǎn)與中心臨時(shí)噪聲像素點(diǎn)的像素差值，其和值記為 S，隨后執(zhí)行步驟3）；如果df≤w ，執(zhí)行步驟2）。

2. 令 w=w+2，執(zhí)行步驟 1），濾波窗口的長(zhǎng)度逐步遞增 2，但需保證濾波窗口的長(zhǎng)度最大值小于等于W。

3. 結(jié)合圖像的最大像素灰度值MAX、最小像素灰度值MIN以及點(diǎn)（i, j）的灰度值 si,j，確定閾值T3。若S小于 T3，則認(rèn)為該點(diǎn)為非噪聲像素點(diǎn)，同時(shí)在噪聲標(biāo)示矩陣 r（i, j）中修改為 r（i, j ） =0；反之，則該點(diǎn)是噪聲像素點(diǎn)，保持矩陣 r（i, j） 中該點(diǎn)的灰度值不變。

4. 將w×w濾波窗口移至下一個(gè)臨時(shí)噪聲像素點(diǎn)，按照步驟 1-3判斷下一個(gè)噪聲像素點(diǎn)是否是噪聲像素，對(duì)所有臨時(shí)噪聲像素點(diǎn)逐一判斷，排除非噪聲像素點(diǎn)到最終的噪聲像素點(diǎn)。

（3）噪聲像素濾除階段

經(jīng)過(guò)階段一和階段二的檢測(cè)與核實(shí)后，整個(gè)圖像的像素點(diǎn)被分為噪聲像素點(diǎn)和非噪聲像素點(diǎn)（信號(hào)像素）兩類(lèi)。對(duì)于非噪聲像素點(diǎn)，無(wú)需濾波，保持其灰度值不變即可；對(duì)于噪聲像素點(diǎn)，需利用濾波窗口中的非噪聲像素點(diǎn)進(jìn)行中值濾波，具體過(guò)程如下：

對(duì)圖像中任意噪聲像素點(diǎn)（i, j），以該點(diǎn)為中心，采用大小為N×N的濾波窗口，將窗口內(nèi)所有非噪聲像素點(diǎn)的集合定義為集合 Ei,j。

對(duì) Ei,j中灰度值進(jìn)行升序排列，得到集合 FI,J，并提取集合 FI,J中像素點(diǎn)的排序中值 Mi,j。濾波后Ei,j中不存在噪聲像素點(diǎn)，只保留了非噪聲像素，降低了噪聲像素對(duì)濾波結(jié)果的影響。

3 濾波算法實(shí)驗(yàn)結(jié)果與分析

本文以matlab 為實(shí)驗(yàn)工具，對(duì)兩副600 MW熱力發(fā)電鍋爐火焰電視采集的煤粉燃燒火焰圖像的灰度圖像為實(shí)驗(yàn)對(duì)象進(jìn)行濾波處理，其中 3.1中實(shí)驗(yàn)圖像為鍋爐啟動(dòng)階段時(shí)采集且經(jīng)灰度變換后圖像，其中爐膛內(nèi)燃燒器為濃縮型 EI-XCL燃燒器，與高負(fù)荷運(yùn)行的鍋爐相比，該圖像可清晰看到點(diǎn)火槍的燃燒狀態(tài)；3.2中實(shí)驗(yàn)圖像為爐膛低負(fù)荷運(yùn)行狀態(tài)（圖像采集時(shí)的爐膛負(fù)荷為396 MW）時(shí)采集且經(jīng)灰度變換后圖像，相比較 3.1中圖像，該圖像提供的圖像細(xì)節(jié)較為簡(jiǎn)單。

3.1 爐膛啟動(dòng)階段煤粉燃燒灰度圖像濾波實(shí)驗(yàn)

本章以 600MW 熱力發(fā)電鍋爐啟動(dòng)階段火焰電視采集的煤粉燃燒火焰圖像的灰度圖像為實(shí)驗(yàn)對(duì)象，分別對(duì)該圖像添加噪聲方差0.3與0.6的椒鹽噪聲，利用自適應(yīng)中值濾波算法與ABDND算法分別對(duì)添加噪聲后圖像進(jìn)行濾波處理。

對(duì)圖1（a）中啟動(dòng)階段煤粉燃燒灰度圖像添加噪聲方差 0.3的椒鹽噪聲，并利用自適應(yīng)中值濾波算法與ABDND算法對(duì)其進(jìn)行濾波處理，實(shí)驗(yàn)結(jié)果如圖2所示。

圖1 啟動(dòng)階段煤粉燃燒灰度圖像&其灰度值分布直方圖Fig.1 Pulverized coal combustion gray image during start-up & Gray value distribution histogram

圖2 添加噪聲方差0.3的椒鹽噪聲后圖像&濾波后圖像濾波結(jié)果&灰度值分布直方圖Fig.2 Add a noise variance of 0.3 after the image & Algorithm filtering results & Each Gray value distribution histogram

實(shí)驗(yàn)結(jié)果顯示，在噪聲方差為 0.3時(shí)，中值濾波算法濾波后圖像噪聲像素明顯減少，但相比原圖像仍存在噪聲像素，ABDND算法濾波后圖像基本不存在噪聲像素，濾波后圖像的灰度值分布直方圖基本與原圖灰度值分布直方圖基本吻合。

對(duì)圖1（a）中啟動(dòng)階段煤粉燃燒灰度圖像添加噪聲方差 0.6的椒鹽噪聲，并利用自適應(yīng)中值濾波算法與ABDND算法對(duì)其進(jìn)行濾波處理，實(shí)驗(yàn)結(jié)果如圖3所示。

實(shí)驗(yàn)結(jié)果顯示，當(dāng)噪聲方差增大到0.6時(shí)，自適應(yīng)中值濾波算法濾波后圖像仍存在大量噪聲像素，濾波后圖像灰度值分布直方圖在最大灰度值像素的數(shù)量分布上與原圖直方圖相差較大，而ABDND算法濾波后圖像濾波效果遠(yuǎn)優(yōu)于自適應(yīng)濾波算法，濾波后圖像灰度值分布直方圖與原圖基本吻合。

3.2 爐膛運(yùn)行階段煤粉燃燒灰度圖像濾波實(shí)驗(yàn)

本章以 600MW 熱力發(fā)電鍋爐運(yùn)行階段階段火焰電視采集的煤粉燃燒火焰圖像的灰度圖像為實(shí)驗(yàn)對(duì)象，分別對(duì)該圖像添加噪聲方差0.3與0.6的椒鹽噪聲，利用自適應(yīng)中值濾波算法與ABDND算法分別對(duì)添加噪聲后圖像進(jìn)行濾波處理。

圖3 添加噪聲方差0.6的椒鹽噪聲后圖像&濾波后圖像濾波結(jié)果&灰度值分布直方圖Fig.3 Add a noise variance of 0.6 after the image & Algorithm filtering results & Each Gray value distribution histogram

圖4 運(yùn)行階段煤粉燃燒灰度圖像&灰度值分布直方圖Fig.4 Pulverized coal combustion gray image during furnace operation & Gray value distribution histogram

對(duì)圖4（a）中啟動(dòng)階段煤粉燃燒灰度圖像添加噪聲方差 0.3的椒鹽噪聲，并利用自適應(yīng)中值濾波算法與ABDND算法對(duì)其進(jìn)行濾波處理，實(shí)驗(yàn)結(jié)果如圖5所示：

圖5 添加噪聲方差0.3的椒鹽噪聲后圖像&濾波后圖像濾波結(jié)果&灰度值分布直方圖Fig.5 Add a noise variance of 0.3 after the image & Algorithm filtering results & Each Gray value distribution histogram

實(shí)驗(yàn)結(jié)果顯示，在噪聲方差為 0.3時(shí)，中值濾波算法濾波后圖像噪聲像素明顯減少，但相比原圖像仍存在噪聲像素，且噪聲像素分布明顯，ABDND算法濾波后圖像基本不存在噪聲像素，濾波后圖像的灰度值分布直方圖與原圖灰度值分布直方圖基本吻合。

對(duì)圖4（a）中啟動(dòng)階段煤粉燃燒灰度圖像添加噪聲方差 0.6的椒鹽噪聲，并利用自適應(yīng)中值濾波算法與ABDND算法對(duì)其進(jìn)行濾波處理，實(shí)驗(yàn)結(jié)果如圖6所示。

實(shí)驗(yàn)結(jié)果顯示，當(dāng)噪聲方差增大到 0.6時(shí)，自適應(yīng)中值濾波算法濾波后圖像仍存在大量噪聲像素，但由于原圖灰度值分布直方圖中灰度值分布較密集，濾波后圖像灰度值分布直方圖與原圖直方圖雖存在差異，但不明顯，而ABDND算法濾波后圖像濾波效果遠(yuǎn)優(yōu)于自適應(yīng)濾波算法，算法濾波后圖像雖存在噪聲像素，但噪聲像素點(diǎn)較少，濾波后圖像灰度值分布直方圖與原圖基本吻合。

本章通過(guò)對(duì)兩幅不同運(yùn)行狀態(tài)的爐膛煤粉燃燒圖像的灰度圖像添加兩種不同方差的椒鹽噪聲，并對(duì)噪聲圖像利用自適應(yīng)濾波算法與ABDND算法進(jìn)行濾波處理，直觀地反映了兩種算法的濾波效果，實(shí)驗(yàn)結(jié)果顯示，相比較自適應(yīng)濾波算法，ABDND算法在火焰灰度圖像噪聲濾波中效果良好

4 標(biāo)準(zhǔn)差數(shù)值對(duì)比

本章通過(guò)計(jì)算原灰度圖像、中值濾波算法濾波后圖像與ABDND算法濾波后圖像的標(biāo)準(zhǔn)差來(lái)客觀評(píng)定兩種算法，濾波后圖像的標(biāo)準(zhǔn)差越接近原圖像，則說(shuō)明算法濾波效果越好，具體計(jì)算方式如下：

圖6 添加噪聲方差0.6的椒鹽噪聲后圖像&濾波后圖像濾波結(jié)果&灰度值分布直方圖Fig.6 Add a noise variance of 0.6 after the image & Algorithm filtering results & Each Gray value distribution histogram

其中I（i,j）為圖像 I（x, y ）中的某像素點(diǎn)，M與N為圖像的長(zhǎng)和寬，W為圖像均值，б為圖像標(biāo)準(zhǔn)差。經(jīng)過(guò)計(jì)算統(tǒng)計(jì)后，兩種濾波算法的濾波后圖像標(biāo)準(zhǔn)差如表1與表2所示：

表1 爐膛啟動(dòng)階段圖像標(biāo)準(zhǔn)差&濾波后圖像標(biāo)準(zhǔn)差Tab.1 Standard deviation of pulverized coal combustion gray image during start-up & Filtered image standard deviation

從表1與表2的標(biāo)準(zhǔn)差計(jì)算數(shù)據(jù)中可以對(duì)比看出，對(duì)于爐膛啟動(dòng)階段煤粉燃燒灰度圖像，當(dāng)噪聲方差大小為 0.2時(shí)，兩種算法濾波后圖像標(biāo)準(zhǔn)值接近原圖標(biāo)準(zhǔn)值，但當(dāng)噪聲方差增大時(shí)，中濾波算法濾波后圖像的標(biāo)準(zhǔn)差異與原圖標(biāo)準(zhǔn)差的差值開(kāi)始增大且增大幅度較大，而 ABDND算法濾波后圖像標(biāo)準(zhǔn)差仍接近原圖標(biāo)準(zhǔn)差，當(dāng)噪聲方差增大到0.5時(shí)，濾波后圖像標(biāo)準(zhǔn)差開(kāi)始增大但幅度較?。粚?duì)于爐膛運(yùn)行階段煤粉燃燒灰度圖像，由于圖像較為簡(jiǎn)單灰度分布較為集中，中值濾波算法濾波后圖像的標(biāo)準(zhǔn)差相比原圖像差值較大，濾波結(jié)果不理想，ABDND算法濾波后圖像在噪聲方差大小為 0.3及小于 0.3時(shí)的圖像標(biāo)準(zhǔn)差接近原圖，但當(dāng)噪聲方差從 0.3開(kāi)始增大后其濾波后圖像標(biāo)準(zhǔn)差與原圖相比差值開(kāi)始增大，但增大幅度較小。

圖7 爐膛啟動(dòng)階段圖像%濾波后圖像標(biāo)準(zhǔn)差折線(xiàn)分布圖Fig.7 Sandard deviation of pulverized coal combustion gray image during start-up &Filtered image standard deviation line profile

表2 爐膛運(yùn)行階段圖像標(biāo)準(zhǔn)差&濾波后圖像標(biāo)準(zhǔn)差Tab.2 Standard deviation of pulverized coal combustion gray image during furnace operation & Filtered image standard deviation

圖8 爐膛運(yùn)行階段圖像%濾波后圖像標(biāo)準(zhǔn)差折線(xiàn)分布圖Fig.8 Standard deviation of pulverized coal combustion gray image during furnace operation & Filtered image standard deviation

5 總結(jié)

本文通過(guò)列舉經(jīng)典濾波算法以及ABDND算法的具體實(shí)現(xiàn)方法，對(duì)比其算法實(shí)現(xiàn)過(guò)程，并利用matlab對(duì)應(yīng)用廣泛的自適應(yīng)濾波算法以及ABDND算法進(jìn)行編譯實(shí)現(xiàn)，通過(guò)對(duì)兩幅不同運(yùn)行階段的爐膛煤粉燃燒灰度圖像進(jìn)行添加噪聲與濾波處理，經(jīng)過(guò)濾波效果圖直觀對(duì)比與標(biāo)準(zhǔn)差數(shù)據(jù)的客觀對(duì)比后得出結(jié)論：ABDND算法可良好完成爐膛火焰圖像在高噪聲密度環(huán)境下（噪聲方差小于等于 0.5）的圖像濾波工作，為爐膛監(jiān)控系統(tǒng)提供清晰的實(shí)時(shí)圖像，且圖像細(xì)節(jié)保留完整，可保證爐膛的穩(wěn)定安全運(yùn)行。

[1] 張欣, 火電廠煤粉鍋爐燃燒狀態(tài)只能監(jiān)測(cè)與評(píng)判研究[D].保定: 華北電力大學(xué), 2009.

[2] 劉芳, 基于數(shù)字圖像處理的爐膛火焰檢測(cè)及燃燒診斷系統(tǒng)研究[D]. 西安: 陜西科技大學(xué), 2008.

[3] 李保根, 600MW火力發(fā)電機(jī)組培訓(xùn)教材(鍋爐篇)電廠鍋爐設(shè)備及系統(tǒng) [M]. 云南: 滇東能源, 2008.

[4] 張曉青, 田書(shū)榕, 肖蘆洋, 等. 基于Simulink 的外調(diào)制微波光纖延遲線(xiàn)噪聲分析[J]. 軟件, 2013, 34(1): 79-81.

[5] 王佳飛, 張強(qiáng), 彭向偉. 一種基于ZSP800 核改進(jìn)型LMS自適應(yīng)濾波算法的實(shí)現(xiàn)[J]. 軟件, 2013, 34(12): 139-141.

[6] 王賽男, 邢冬梅. 基于MATLAB 的小波閾值濾波與應(yīng)用[J]. 軟件, 2015, 36(10): 102-104.

[7] 劉聞, 別紅霞. 基于蟻群算法的噪聲圖像邊緣檢測(cè)[J]. 軟件, 2013, 34(12): 256-259.

[8] 陳國(guó)強(qiáng). 基于GPU 的圖像處理算法研究[J]. 軟件, 2014,35(2): 135-136.

[9] 林王兵, 許燕, 韓飛, 呂超賢, 龔佳俊等. 實(shí)驗(yàn)室環(huán)境下的邊緣檢測(cè)技術(shù)研究綜述[J]. 軟件, 2015, 36(9): 29-32.

[10] 鐘彩. 邊緣檢測(cè)算法在圖像預(yù)處理中的應(yīng)用[J]. 軟件,2013, 34(1): 158-159.

[11] Chang C C, Hsiao j Y, A adaptive median filter for image denoising[J]. IEEE. Washington DC:IEEE Computer society,2008:346-350.

[12] Huang T S, Zhang G T, A Fast Two Dimentional Median Foltering Algorithm[J]. EE. Trans. Acoust. Speech, Sig.Proce, 1979: 13-18.

[13] Nallaperumal, K., Varghese, J., Saudia, S.: Salt&peppr impulse noise removal using adaptive switching median filter.Proc. IEEE Oceans, Asia Pacific, Singapore, 2006.

[14] Ng, P. E., Ma, K. K.: ‘A Switching median filter with boundary discriminative noise detection for extremely corrupted images’, IEEE Trans. Image Process. , 2006, 15, (6),pp. 1506–1516.

[15] A. K. Tripathi U. Ghanekar S. Mukhopadhyay, (Y)Switching median filter: advanced boundary discriminative noise detection algorithm IEY Image Processing, 2010.