岳春宇，江萬壽

武漢大學(xué)測繪遙感信息工程國家重點(diǎn)實(shí)驗(yàn)室，湖北武漢430079

1 引 言

合成孔徑雷達(dá)（synthetic aperture radar，SAR）圖像是由接收散射回波信號的相干獲得的，因此會產(chǎn)生隨機(jī)散布的相干斑（speckle），降低了目標(biāo)探測和解譯的能力，增加了圖像處理的難度。傳統(tǒng)的SAR圖像去噪方法如增強(qiáng)Lee濾波［1］、Gamma－MAP濾波［2］等主要利用圖像空域的統(tǒng)計(jì)特征。由于對圖像局部區(qū)域的假設(shè)有時與實(shí)際情況不完全相符，所以在去噪的同時不能較好地保留邊緣和紋理細(xì)節(jié)。因而只通過對圖像域的灰度信息進(jìn)行統(tǒng)計(jì)建模不能充分體現(xiàn)SAR圖像多尺度多分辨等特性，對SAR圖像去噪是遠(yuǎn)遠(yuǎn)不夠的［3］?；谛〔ㄓ蚪y(tǒng)計(jì)特性的方法［4］在多尺度多分辨的統(tǒng)計(jì)模型下分析相干斑噪聲分布特性，能夠保持圖像的細(xì)節(jié)以及輪廓特征。然而小波基只對于具有點(diǎn)狀奇異性的目標(biāo)函數(shù)是最優(yōu)基，小波系數(shù)是稀疏的，對于具有線狀奇異的函數(shù)，小波系數(shù)則不再稀疏［5］。多尺度幾何分析（multiscale geometric analysis，MGA）方法是繼小波變換之后，提出的又一類新的圖像表示方法［6］，能夠在保持小波分析優(yōu)良的時（空）頻局部化特性的同時，通過構(gòu)造具有各向異性的基來檢測小波所不能充分刻畫的幾何結(jié)構(gòu)，很好地保持圖像的奇異性信息，從而有利于圖像去噪質(zhì)量的提高［7］。

輪廓波變換（contourlet transform，CT）［8］變換是一種MGA方法，在SAR圖像去噪中也得到成功應(yīng)用［3］，但由于其缺乏平移不變性，因此在采取收縮法去噪時會引入Gibbs噪聲［6］。文獻(xiàn)［9］在輪廓波變換的基礎(chǔ)上采用非下采樣金字塔分解和非下采樣方向?yàn)V波器組（directional filter banks，DFB）分解，提出了非下采樣輪廓波變換（nonsubsampled contourlet transform，NSCT）。它是一種新的具有完全平移不變性的MGA方法，已在SAR圖像去噪處理中顯示出其優(yōu)勢［7，10－12］。由于非下采樣變換域去噪方法比下采樣變換域去噪方法提供了更加豐富、冗余的圖像細(xì)節(jié)信息，而且具有平移不變性，所以前者的性能一般高于后者［7］。但是目前基于NSCT的方法由于閾值收縮，會產(chǎn)生“裂痕”噪聲，并且大多適用于機(jī)載SAR影像，對星載SAR數(shù)據(jù)處理效果還不是很理想。

筆者根據(jù)SAR圖像數(shù)據(jù)的特征，引入了非對數(shù)加性模型，并在該模型下對SAR圖像NSCT域中的噪聲分布統(tǒng)計(jì)建模，應(yīng)用最大后驗(yàn)（MAP）準(zhǔn)則和在NSCT域中的non－local（NL）約束相結(jié)合的方法解求SAR圖像真實(shí)信號的NSCT系數(shù)，避免了裂痕噪聲，增強(qiáng)了濾波視覺效果。最后使用本文方法與傳統(tǒng)NSCT去噪方法和目前主流去噪方法對真實(shí)星載SAR圖像進(jìn)行去噪試驗(yàn)，證明了本文方法的有效性。

2 SAR圖像的非對數(shù)加性相干斑模型

SAR圖像相干斑模型先驗(yàn)信息對SAR圖像噪聲抑制具有重要意義。在圖像域，完全發(fā)育的相干斑噪聲模型可以表示為

式中，（x，y）為圖像坐標(biāo)；I為觀測的SAR圖像強(qiáng)度；X為真實(shí)后向散射強(qiáng)度；R為產(chǎn)生相干斑噪聲的過程，是一個均值為1，方差與等效視數(shù)成反比，服從Γ分布的二階隨機(jī)過程。R和X是兩個相互獨(dú)立的隨機(jī)過程。

在乘性噪聲模型下，通常去噪的思路是首先對圖像進(jìn)行對數(shù)變換，使圖像轉(zhuǎn)換為加性模型，假設(shè)加性對數(shù)噪聲為零均值高斯白噪聲進(jìn)行去噪。但是這樣處理不能很好地保持原圖像的輻射特性，變換后的相干斑噪聲并不是零均值的，這會使去噪后的結(jié)果圖像的均值偏離原始圖像的均值，圖像在視覺上會變暗［13］。同時對數(shù)模型計(jì)算的時間代價也非常大。

因此，對引入SAR圖像非對數(shù)加性模型［14］，乘性模型式（1）改寫為

N（x，y）＝［R（x，y）－1］X（x，y）可認(rèn)為是和實(shí)際真實(shí)反射強(qiáng)度X（x，y）相關(guān)的噪聲。文獻(xiàn)［10］證明了平穩(wěn)區(qū)域噪聲N（x，y）的NSCT系數(shù)符合零均值高斯分布。

3 非下采樣輪廓波變換

輪廓波變換是一種圖像多尺度幾何表示方法［8］。由拉普拉斯塔式分解（Laplacian pyramid，LP）和方向?yàn)V波器組（DFB）濾波兩部分實(shí)現(xiàn)。先由LP把原始圖像分解為低頻子帶和高頻子帶，然后高頻子帶經(jīng)過DFB分解為多個方向子帶。對低頻子帶重復(fù)上述過程即實(shí)現(xiàn)了對原始圖像的多分辨率、多方向分解。圖1（a）為CT變換的濾波器組結(jié)構(gòu)。CT變換具有良好的各向異性，基的支撐區(qū)間是長寬比隨尺度變換的“長方形”結(jié)構(gòu)，能夠沿圖像中的輪廓邊緣用最少的系數(shù)來逼近奇異曲線。但是由于下采樣機(jī)制，在LP分解時，對圖像進(jìn)行了隔行隔列下采樣，這在較大程度地降低變換冗余性的同時導(dǎo)致了CT變換不具有平移不變性，限制了其在一些圖像處理領(lǐng)域里的應(yīng)用。

因此，文獻(xiàn)［9］又提出了NSCT變換，保留了CT變換的頻率分割結(jié)構(gòu)，去掉原CT變換中的下采樣步驟，采用非下采樣塔式分解（NSPL）和構(gòu)造非下采樣濾波器組（NSDFB）。圖1（b）為NSCT變換的濾波器組結(jié)構(gòu)。NSCT具有更加豐富的基函數(shù)集，提供了更好的頻率選擇性和正則性，有利于更好地捕捉圖像中的細(xì)節(jié)信息，同時保證了平移不變性。

圖1 Fig.1

4 NSCT域最大后驗(yàn)和非局域約束噪聲抑制

圖像中的特征一般具有方向性，經(jīng)過NSCT變換后，具有方向性的空間域的圖像特征只分布在少數(shù)幾個子帶的NSCT系數(shù)上。而噪聲一般不具有方向性，經(jīng)過NSCT變換后會均勻分布在所有方向子帶。因而可以在NSCT域?qū)D像統(tǒng)計(jì)建模，去噪處理。

4.1 SAR圖像NSCT域MAP多層閾值估計(jì)［10］

對SAR圖像的一條子帶的NSCT變換在非對數(shù)加性模型下可以表示為

式中，CI、CX、CN分別為觀測圖像、真實(shí)反射強(qiáng)度、噪聲的NSCT系數(shù)；NSCT（·）為NSCT變換過程。

圖像可看成是由各個很小的同質(zhì)區(qū)域組成［15］，文獻(xiàn)［16］指出在同質(zhì)區(qū)域內(nèi)的地物真實(shí)后向散射強(qiáng)度是一個常數(shù)，這種小的同質(zhì)區(qū)域的極限情況就是單個像素為一同質(zhì)區(qū)域。與圖像空域同質(zhì)區(qū)域相對應(yīng)的是NSCT域高頻子帶中的平穩(wěn)區(qū)域，文獻(xiàn)［10］證明了在NSCT域高頻子帶中的平穩(wěn)區(qū)域CX、CN分別可用零均值的高斯分布和零均值的Laplace分布建模。CX的估計(jì)值根據(jù)最大后驗(yàn)（MAP）準(zhǔn)則由CX和CI計(jì)算得到

式中，P（CI）、P（CX）、P（CN）分別為CI、CX、CN的概率密度函數(shù)。式（4）可以表示為

上式對CX求導(dǎo)可得

令其等于0，求得

由文獻(xiàn)［10］可知sign（CX）與sign（CI）相同，由文獻(xiàn)［11］可確定和σCX。設(shè)閾值對NSCT域各條高頻子帶中的平穩(wěn)區(qū)域進(jìn)行軟閾值去噪。為了下面對全部系數(shù)NL優(yōu)化，保持系數(shù)的獨(dú)立性，在實(shí)際計(jì)算中，取極端情況，每個NSCT系數(shù)看做一個平穩(wěn)區(qū)域。

4.2 NSCT域系數(shù)的非局域約束

對于各個子帶分別軟閾值去噪后重建的結(jié)果存在著裂痕噪聲，并由于閾值法對系數(shù)的“過扼殺”［6］，導(dǎo)致圖像在邊緣處出現(xiàn)模糊。

實(shí)際上，某些NSCT系數(shù)只含有噪聲，而某些NSCT系數(shù)則都為有效信號，對某一子帶采用一個單一的模型建模，設(shè)置同一閾值進(jìn)行收縮，必然會造成去噪不徹底或圖像失真。圖像經(jīng)NSCT由各向異性的基分解后，信號的幾何結(jié)構(gòu)得到了充分保留，實(shí)際信號會表現(xiàn)為很強(qiáng)的相關(guān)性，噪聲因?yàn)槭请S機(jī)分布的表現(xiàn)為弱相關(guān)或者不相關(guān)［17］。因而可以根據(jù)NSCT分解后系數(shù)的相關(guān)性對圖像進(jìn)行去噪后的優(yōu)化處理。

在實(shí)際計(jì)算中，每個NSCT系數(shù)不僅與同子帶系數(shù)之間具有相關(guān)性，同時不同子帶和上下層子帶之間與同一空間位置對應(yīng)的NSCT系數(shù)由于表示同一信號不同層次不同方向的分解，也具有相關(guān)性。在同一空間位置的所有NSCT系數(shù)的關(guān)系可以表示為

文獻(xiàn)［18］中提出的NL濾波方法在對某一信號濾波時同時考慮了本區(qū)域以及相關(guān)鄰域的關(guān)系，能夠考慮到與該信號相似信號對其的影響。對于軟閾值收縮后的各子帶的NSCT某一系數(shù)Ci，其與同子帶其他系數(shù)具有一定相關(guān)性，借鑒文獻(xiàn)［18］的方法，有

式中，Mi和Mj分別是包含Ci和Cj的兩個區(qū)域；權(quán)值w（i，j）由Mi與Mj的相似程度決定。由于NSCT的基是各項(xiàng)異性的，所以直接使用NSCT系數(shù)不能很好的表示他們之間關(guān)系。加入NSCT基α，式（8）可以寫為

令v（Mi）、v（Mj）分別是Mi和Mj的系數(shù)向量，可得

式中，‖v（Mi）－v（Mj）‖為一個加權(quán)歐式距離遞減函數(shù)，a＞0為高斯核的標(biāo)準(zhǔn)差；N（i）＝為正則化系數(shù)；h為濾波尺度系數(shù)。

由式（7）和式（9）可得

對每一子帶NSCT系數(shù)均作NL濾波，對于表示實(shí)際信號的系數(shù)，與其具有相關(guān)性的系數(shù)會提供支持，而噪聲則會進(jìn)一步削弱。

5 算法實(shí)現(xiàn)

步驟1：設(shè)置NSCT變換分解層數(shù)以及各子帶方向數(shù)。對輸入SAR圖像進(jìn)行NSCT變換，得到NSCT系數(shù)。

步驟2：基于MAP的BayesShrink閾值估計(jì)，根據(jù)文中3.1所述方法分層計(jì)算每一個子帶的收縮閾值，分層多閾值對各個子帶軟閾值收縮去噪。

步驟3：對各子帶軟閾值收縮后的結(jié)果，根據(jù)4.2中所述方法作non－local平滑約束。

步驟4：對去噪后的NSCT系數(shù)做NSCT逆變換，得到最終的去噪圖像。

6 試驗(yàn)及評價

本文采用兩組星載SAR圖像數(shù)據(jù)驗(yàn)證所提出去噪方法的有效性。第一組為TerraSAR－X的SSC圖像（256×256），主要驗(yàn)證本方法對人工地物和水體的去噪效果；第二組為Cosmo－SkyMed的SCS圖像（256×256），主要驗(yàn)證本方法對田地的去噪效果。試驗(yàn)中采用了增強(qiáng)Lee［1］、BM3D［19］、PPB［20］、基于最大后驗(yàn)的NSCT［10］和本文方法作對比試驗(yàn)。

試驗(yàn)中，增強(qiáng)Lee濾波采用3×3窗口，BM3D濾波采用文獻(xiàn)［19］的方法，PPB濾波采用文獻(xiàn)［20］方法迭代25次，基于最大后驗(yàn)的NSCT濾波采用文獻(xiàn)［10］的方法，基于最大后驗(yàn)的NSCT和本文方法中NSCT變換選擇“max flat” NSPL分解和“dmax flat7”NSDFB濾波器組，為了能夠較好的保存圖像的細(xì)節(jié)，NSCT作4層分解，方向數(shù)為4、8、16、16。試驗(yàn)結(jié)果見圖2、圖3。

圖2 TerraSAR－X圖像去噪結(jié)果Fig.2 The denoising result of TerraSAR－X image

圖3 Cosmo－SkyMed圖像去噪結(jié)果Fig.3 The denoising result of Cosmo－SkyMed image

SAR圖像的相干斑抑制過程是盲的恢復(fù)過程，PSNR（peak signal to noise ratio）等不能用來評價去噪結(jié)果［3］。本文采用等效視數(shù)（equivalent number of looks，ENL）、比值圖像均值［10］，作為去噪結(jié)果客觀的評價指標(biāo)，PSNR反映了濾波后圖像與原始圖像之間的差異，用來評價濾波圖像分辨率的變化和損失。試驗(yàn)數(shù)據(jù)的試驗(yàn)結(jié)果見表1、表2。

表1 不同濾波方法TerraSAR－X的試驗(yàn)結(jié)果Tab.1 The denoising result of all the filter methods for TerraSAR－X image

表2 不同濾波方法Cosmo－SkyMed的試驗(yàn)結(jié)果Tab.2 The denoising result of all the filter methods for Cosmo－SkyMed image

等效視數(shù)用來衡量一幅圖像斑點(diǎn)噪聲相對強(qiáng)度，目前被廣泛接受為SAR相干斑抑制的指標(biāo)。等效視數(shù)越大，去噪效果越好。本文方法去噪結(jié)果的等效視數(shù)高于其他方法，這與目視效果一致。在理想情況下，原始圖像與去噪后結(jié)果圖像的比值圖像為相干斑噪聲圖像。而相干斑噪聲圖像的均值為1。所以比值圖像的均值越接近1，去噪結(jié)果越好，對原始圖像輻射特性保持得越好。本文方法對應(yīng)的比值圖像的均值最接近1，說明對原始圖像輻射特性保持最好。PSNR反映了濾波結(jié)果與原始圖像之間的差異，PSNR值越大則與原始圖像差異越小，因?yàn)槿ピ氲耐瑫r犧牲了圖像的細(xì)節(jié)和分辨率，所以要辯證地看待PSNR值。本文方法PSNR值與增強(qiáng)Lee和PPB方法接近，即在分辨率和細(xì)節(jié)保持上與一些經(jīng)典方法相近，而PSNR值低于某些方法的同時，去噪效果也優(yōu)于這些方法。

從主觀上評價，在兩組試驗(yàn)數(shù)據(jù)中，本文方法對水體和平坦地區(qū)的去噪結(jié)果明顯優(yōu)于其他方法，在同質(zhì)區(qū)域產(chǎn)生的偽吉布斯效應(yīng)是最弱的，同時具有較強(qiáng)的細(xì)節(jié)保持能力，并且改善了基于最大后驗(yàn)的NSCT去噪方法在某些局部區(qū)域會產(chǎn)生劃痕的缺點(diǎn)。由基于最大后驗(yàn)的NSCT濾波方法和本文方法的結(jié)果相比較可以看出，雖然本文方法結(jié)果更加平滑，但由于NL約束的引入，計(jì)算過程中NSCT系數(shù)多作了一次平滑處理，所以對結(jié)果的細(xì)節(jié)紋理保持產(chǎn)生了一定的影響。

7 結(jié) 論

本文針對SAR圖像對數(shù)加性模型不能夠很好保持原圖像的輻射特性和計(jì)算時間代價大的問題引入了非對數(shù)加性模型?；谖墨I(xiàn)［9］提出的方法，利用NSCT對SAR圖像細(xì)節(jié)信息的刻畫能力［12］對SAR圖像作NSCT分解。然后對SAR圖像同質(zhì)區(qū)域在非對數(shù)加性模型下的噪聲分布應(yīng)用高斯分布建模，并基于MAP準(zhǔn)則計(jì)算每個高頻子帶中真實(shí)信號NSCT系數(shù)的閾值，分層多閾值去噪。針對頻域閾值去噪對真實(shí)信號“過扼殺”的問題，考慮到在頻域真實(shí)信號NSCT之間的相關(guān)性，借鑒NL理論，對每個子帶閾值去噪的結(jié)果進(jìn)行NL濾波，進(jìn)一步增強(qiáng)真實(shí)信號的NSCT系數(shù)，削弱殘存噪聲的影響。

通過試驗(yàn)，與經(jīng)典的增強(qiáng)Lee濾波算法和目前效果較好的PPB去噪方法等相比，在客觀指標(biāo)上本文方法在原圖像輻射特性保持和去噪效果上均優(yōu)于其他方法，在主觀目視效果上，在SAR圖像的同質(zhì)區(qū)域去噪效果明顯，紋理結(jié)構(gòu)和細(xì)節(jié)也得到有效保持。

但是從試驗(yàn)中也可以看到，由于NL方法的引入，在解決了傳統(tǒng)的NSCT閾值去噪方法產(chǎn)生裂痕噪聲和邊緣模糊的同時，即使在同質(zhì)區(qū)域噪聲去除效果良好，算法的計(jì)算時間代價要小于PPB方法，圖像某些細(xì)節(jié)部分不如PPB方法清晰。因此對去噪后圖像細(xì)節(jié)紋理的保持增強(qiáng)還需作進(jìn)一步的研究。

［1］ LEE J S.Digital Image Enhancement and Noise Filtering by Use of Local Statistics［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，1980，PAMI－2（2）：165－168.

［2］ LOPES A，NEZRY E，TOUZI R，et al.Maximum a Posteriori Filtering and First Order Texture Models in SAR Images［C］∥Proceedings of IGARSS’90.Washington，D C：IEEE，1990：2409－2412.

［3］ SHA Yuheng，CONG Lin，SUN Qiang，et al.SAR Image Despeckling Based on Contourlet Domain Hidden Markov Trees Model［J］.Journal Infrared Millim Waves，2009，28（1）：66－71.（沙宇恒，叢琳，孫強(qiáng)，等.基于Contourlet域HMT模型的SAR圖像相干斑抑制［J］.紅外與毫米波學(xué)報，2009，28（1）：66－71.）

［4］ ZHANG Jun，LIU Jian.A Speckle Reduction Algorithm by Soft－thresholding Based on Wavelet Filters for SAR Images［J］.Acta Geodaetica et Cartographica Sinica，1998，27（2）：119－124.（張俊，柳健.SAR圖像斑點(diǎn)噪聲的小波軟門限濾除算法［J］.測繪學(xué)報，1998，27（2）：119－124.）

［5］ JIAO Licheng，TAN Shan.Development and Prospect of Image Multiscale Geometric Analysis［J］.Acta Electronica Sinica，2003，31（12A）：43－50.（焦李成，譚山.圖像多尺度幾何分析：回顧和展望.電子學(xué)報，2003，31（12A）：43－50.）

［6］ DENG Chengzhi.Research on Image Sparse Representation Theory and Its Applications［D］.Wuhan：Huazhong University of Science and Technology，2008.（鄧承志.圖像稀疏表示理論及其應(yīng)用［D］.武漢：華中科技大學(xué)，2008.）

［7］ SUN Qiang，GAO Yong，JIAO Licheng.Image Denoising Based on Spatially Adaptive Bayesian Shrinkage in NSCT Domain［J］.Journal of Computer Applications，2010，30（8）：2080－2084.（孫強(qiáng)，高勇，焦李成.基于空間自適應(yīng)Bayesian縮減的NSCT域圖像去噪方法［J］.計(jì)算機(jī)應(yīng)用，2010，30（8）：2080－2084.）

［8］ DO M N，VETTERLI M.The Contourlet Transform：An Efficient Directional Multiresolution Image Representation［J］.IEEE Transaction on Image Processing，2005，14（12）：2091－2106.

［9］ CUNHA A L，ZHOU Jianping，DO M N.The Nonsubsampled Contourlet Transform：Theory，Design，and Applications［J］.IEEE Transaction on Image Processing，2006，15（10）：3089－3101.

［10］ FENG Hongxiao，HOU Biao，JIAO Licheng，et al.SAR Image Despeckling Based on Local Gaussian Model and MAP in NSCT Domain［J］.Acta Electronica Sinica，2010，38（4）：811－816.（鳳宏曉，侯彪，焦李成，等.基于非下采樣Contourlet域局部高斯模型和MAP的SAR圖像相干斑抑制［J］.電子學(xué)報，2010，38（4）：811－816.）

［11］ YANG Xiaohui，JIAO Licheng，NIU Hongjuan，et al.Image Denoising Method for Nonsubsampled Contourlet Based on Multi－threshold［J］.Computer Engineering，2010，36（4）：200－204.（楊曉慧，焦李成，牛宏娟，等.基于多閾值的非下采樣輪廓波圖像去噪方法［J］.計(jì)算機(jī)工程，2010，36（4）：200－204.）

［12］ CHANG Xia，JIAO Licheng，LIU Fang，et al.SAR Image Despeckling Based on the Estimation of Speckle Variance in Nonsubsampled Contourlet Domain［J］.Acta Electronica Sinica，2010，38（6）：1328－1333.（常霞，焦李成，劉芳，等.基于斑點(diǎn)方差估計(jì)的非下采樣Contourlet域SAR圖像去噪［J］.電子學(xué)報，2010，38（6）：1328－1333.）

［13］ XIE H，PIERCE L E，ULABY F T.Statistical Properties of Logarithmically Transformed Speckle［J］.IEEE Transactions on Geoscience and Remote Sensing，2002，40（3）：721－727.

［14］ ARGENTI F，ALPARONE L.Speckle Removal from SAR Images in the Undecimated Wavelet Domain［J］.IEEE Transactions on Geoscience and Remote Sensing，2002，40（11）：2363－2374.

［15］ LEE J S，HOPPEL K，MANGO S.Unsupervised Estimation of Speckle Noise in Radar Images［J］.International Journal of Imaging Systems and Technology，1992，4（4）：298－305.

［16］ LOPES A，TOUZI R，NEZRY E.Adaptive Speckle Filters and Scene Heterogeneity［J］.IEEE Transactions Geoscience and Remote Sensing，1990，28（6）：992－1000.

［17］ ZHANG Xuemeng，JI Fang.The Study of Image Denoising Based on the Wavelet Coefficients Relevance［J］.Journal of Weifang Educational College，2010，23（2）：90－91.（張學(xué)夢，籍芳.基于小波系數(shù)相關(guān)性的圖像去噪研究［J］.濰坊教育學(xué)院學(xué)報，2010，23（2）：90－91.）

［18］ BUADES A，COLL B，MOREL J M.A Non－local Algorithm for Image Denoising［C］∥Proceedings of IEEE Computer Society Conference on Computer Vision and Pattern Recognition：2.Palma de Mallorca：IEEE，2005：60－65.

［19］ DABOV K，F(xiàn)OI A，KATKOVNIK V，et al.Image Denoising by Sparse 3DTransform－domain Collaborative Filtering［J］.IEEE Transactions on Image Processing，2007，16（8）：2080－2095.

［20］ DELEDALLE C A，DENIS L，TUPIN F.Iterative Weighted Maximum Likelihood Denoising with Probabilistic Patchbased Weights［J］.IEEE Transactions on Image Processing，2009，18（12）：2661－2672.