羅曉梅12索志勇3劉且根24

（1.西安電子科技大學(xué)綜合業(yè)務(wù)網(wǎng)理論及關(guān)鍵技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710071；2.南昌大學(xué)信息工程學(xué)院，江西南昌 330031；3.西安電子科技大學(xué)雷達(dá)信號(hào)處理國家重點(diǎn)實(shí)驗(yàn)室，陜西西安710071；4.中國科學(xué)院勞特伯生物醫(yī)學(xué)成像研究中心，廣東深圳 518055）

采用自適應(yīng)字典學(xué)習(xí)的In SAR降噪方法

羅曉梅1，2，索志勇3，劉且根2，4

（1.西安電子科技大學(xué)綜合業(yè)務(wù)網(wǎng)理論及關(guān)鍵技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710071；2.南昌大學(xué)信息工程學(xué)院，江西南昌 330031；3.西安電子科技大學(xué)雷達(dá)信號(hào)處理國家重點(diǎn)實(shí)驗(yàn)室，陜西西安710071；4.中國科學(xué)院勞特伯生物醫(yī)學(xué)成像研究中心，廣東深圳 518055）

提出了一種基于字典學(xué)習(xí)的干涉合成孔徑雷達(dá)相位降噪算法.首先利用字典學(xué)習(xí)，建立了干涉相位濾波的優(yōu)化模型.鑒于該模型非凸難以求解，采用分裂技術(shù)和增廣拉格朗日框架，獲得松弛后的基于l1范數(shù)正則化的優(yōu)化模型，然后引入交替方向乘子法對松弛后的問題求解，獲得最終的相位濾波結(jié)果.通過InSAR復(fù)相位數(shù)據(jù)訓(xùn)練字典，從稀疏表達(dá)式重建所需的復(fù)相位圖像.對仿真數(shù)據(jù)和實(shí)測數(shù)據(jù)的處理顯示這種新的InSAR相位降噪方法在殘點(diǎn)數(shù)、均方誤差和邊緣完整性保持方面優(yōu)于現(xiàn)有的經(jīng)典濾波方法.

InSAR；相位降噪；字典學(xué)習(xí)；l1范數(shù)正則化；交替方向乘子法

干涉合成孔徑雷達(dá)（Interferometric Synthetic Aperture Radar，InSAR）已經(jīng)廣泛應(yīng)用在許多領(lǐng)域，例如:數(shù)字地形模型（Digital Elevation Model，DEM）的獲取和地表形變的監(jiān)測.相位濾波是干涉處理的關(guān)鍵步驟之一，其目的是為了降低或消除各種相干因素引入的相位噪聲，為后續(xù)的相位解纏繞和DEM生成奠定基礎(chǔ)［1］.根據(jù)干涉相位的特點(diǎn)，人們提出了許多有效的相位濾波算法［2-6］.Lee自適應(yīng)濾波器［2］能夠有效地濾除干涉相位圖中的噪聲，其本質(zhì)是對干涉相位的一種加權(quán)處理，然而，該濾波器需要進(jìn)行局部相位解纏繞來確定濾波方向窗和加權(quán)系數(shù)，而存在噪聲的相位解纏繞結(jié)果會(huì)對濾波產(chǎn)生不可忽略的影響.Goldstein濾波器［3］具有很好的噪聲抑制能力，并且是以“塊”濾波的方式進(jìn)行的，具有很高的執(zhí)行效率，目前應(yīng)用比較廣泛.但由于在處理數(shù)據(jù)過程中參數(shù)的不變性，該濾波器在噪聲抑制和邊緣完整性保持性能方面往往不能兼顧.文獻(xiàn)［4］在應(yīng)用區(qū)域增長法以相位變化最小為準(zhǔn)則自適應(yīng)選擇濾波窗口，取得較好濾波效果，但未能考慮條紋數(shù).文獻(xiàn)［5］通過增加45°方向和135°方向的梯度計(jì)算，以及根據(jù)干涉圖的相干系數(shù)確定條紋邊緣的尺度函數(shù)兩方面的改進(jìn)改善權(quán)系數(shù)的選擇，從而改進(jìn)復(fù)數(shù)空間自適應(yīng)濾波效果.但是相干系數(shù)確定的尺度函數(shù)并不是最優(yōu)的.文獻(xiàn)［6］給出了一種利用相干系數(shù)加權(quán)的最優(yōu)干涉相位濾波方法，濾波效果較好.但由于選擇的樣本必須滿足獨(dú)立同分布的條件，而且由于降噪是逐個(gè)像素進(jìn)行的，因此比經(jīng)典的以塊為基礎(chǔ)的降噪濾波器時(shí)間長，例如:Goldstein濾波器.

1 相位濾波建模

為了簡單起見，考慮大小為N×N的InSAR圖像，若將圖像中的像素按列排成列向量，則該圖像在復(fù)域內(nèi)表示成加性噪聲模型為

其中，k表示像素的序列號(hào)，y（k）和u（k）分別表示含噪測量值和理想像素值（復(fù)形式），n（k）表示加性噪聲.引入字典，一幅圖像的每一子塊yi可表示為

其中，yi∈CM2×1，yi表示從原圖像中抽取大小為M×M的含噪圖像子塊按列排列形成的列向量；xi∈CK×1，xi表示對應(yīng)圖像子塊yi的稀疏表示系數(shù)；字典D∈CM2×K，D表示含有K個(gè)原子di∈CM2×1的復(fù)字典；wi∈CM2×1，wi表示均值為零、方差為σ21+jσ22的加性高斯白噪聲，j=（-1）1/2.

眾所周知，由于InSAR相位表現(xiàn)出高的自相似性，因此，可以采用合適的字典基對復(fù)相位數(shù)據(jù)進(jìn)行稀疏表示.通常訓(xùn)練采樣值自相似性越高，字典的稀疏表示能力就越強(qiáng)［11］基于此，以塊為基礎(chǔ)的降噪優(yōu)化模型為

其中，n=（N-M+1）2，表示圖像分塊的塊數(shù).

由于約束條件中存在耦合項(xiàng)，式（4）是非凸的，因而可以采用塊坐標(biāo)輪回法［12］求解.由于字典的引入，采用標(biāo)準(zhǔn)的優(yōu)化軟件CVX［13］直接求解每個(gè)優(yōu)化子問題是非常耗時(shí)的甚至是不能收斂的.因此，設(shè)計(jì)有效的算法求解優(yōu)化問題式（4）是非常必要的.

2 基于ADMM的干涉相位濾波算法

2.1ADMM-DL算法

優(yōu)化問題式（4）是帶有等號(hào)約束的優(yōu)化問題，因此采用ADMM算法［14］求解該問題.采用增廣Lagrangian函數(shù)框架［15］，對式（4）優(yōu)化變量xi和zi（i=1，2，…，n）的追蹤可以等效于求解如下的i=1，2，…，n個(gè)優(yōu)化子問題:

而對于優(yōu)化變量D的追蹤等效于求解如下的優(yōu)化問題:

其中，vi∈CM2×1，vi表示拉格朗日乘子以及β＞0是一個(gè)罰參數(shù).X∈CK×n，X表示由列向量｛xi｝∈CK×1組成的矩陣；Z∈CM2×n，Z表示由列向量｛zi｝∈CM2×1組成的矩陣；Y∈CM2×n，Y表示由列向量｛yi｝∈CM2×1組成的矩陣；V∈CM2×n，V表示由列向量｛vi｝∈CM2×1組成的矩陣.

引入Gaussian-Seidel迭代［12，15］，并采用ADMM技術(shù)，提出的相位濾波算法具體實(shí)現(xiàn)如下:

（1）計(jì)算松弛變量zk+1i，i=1，2，…，n，

（2）計(jì)算稀疏表達(dá)式xk+1i，i=1，2，…，n，

（3）計(jì)算字典Dk+1，

（4）計(jì)算拉格朗日乘子向量vk+1i，i=1，2，…，n，

式（7）～（10）優(yōu)化子問題是凸的，其閉式解或者迭代形式分別為

其中，SShrink表示收縮算子（Shrinkage Operator），定義為

由于篇幅的原因，式（11）和式（12）的證明略.

2.2ADMM-DL算法的計(jì)算復(fù)雜度

由于式（9）～（11）每個(gè)子問題具有閉式解或者簡單的迭代形式，很容易應(yīng)用到真實(shí)的InSAR圖像.通過下一節(jié)對仿真數(shù)據(jù)和實(shí)測數(shù)據(jù)的處理可以獲得驗(yàn)證.ADMM-DL算法的整個(gè)計(jì)算復(fù)雜度分析如下:

由于式（10）～（13）的計(jì)算復(fù)雜度都是O*（n M2K），其中，O*表示低階項(xiàng)可以忽略.總之，ADMM-DL算法的計(jì)算復(fù)雜度是O*（4n M2K）.

3 性能分析

本節(jié)通過對仿真數(shù)據(jù)和實(shí)際數(shù)據(jù)的處理分析筆者所提算法的有效性，并與其他算法進(jìn)行處理性能的對比，對比的指標(biāo)包括濾波后干涉相位圖的殘點(diǎn)數(shù)量［4］，剩余相位的均方差以及相位條紋邊緣的保持性等指標(biāo).

3.1仿真數(shù)據(jù)實(shí)驗(yàn)

通過MATLAB中的peaks函數(shù)產(chǎn)生一個(gè)模擬的相位曲面xyphase=appeak（N），其中，a=10，N=900.假設(shè)復(fù)相位模型的噪聲方差是0.6+j0.6.使用ADMM-DL濾波器和經(jīng)典的濾波方法進(jìn)行濾波.文中進(jìn)行性能對比的經(jīng)典濾波方法包括測試均值濾波器、Goldstein濾波器和Lee濾波器.除了Lee濾波器的窗口大小是9×9，其他濾波方法窗口的大小是5×5.濾波后的圖像如圖1所示.

圖1　仿真結(jié)果對比

為了定量分析不同濾波器的濾波性能，對殘點(diǎn)數(shù)和濾波后剩余相位的均方差（Mean Square Error，MSE）進(jìn)行統(tǒng)計(jì)，如表1所示.

表1　不同濾波器的殘點(diǎn)數(shù)和MSE比較

采用ADMM-DL濾波器，按照式（11）～（13），求解優(yōu)化問題式（3），獲得字典D和最稀疏的表達(dá)式X= ｛xi｝，xi是式（3）的優(yōu)化變量，再將二者相乘恢復(fù)出理想的InSAR相位圖.從圖1所示的干涉圖以及表1所列的殘點(diǎn)數(shù)和MSE可知，ADMM-DL濾波器的殘點(diǎn)數(shù)最少并且其殘余相位的均方差最小.與理想的相位圖相比，采用ADMM-DL濾波器濾波后圖像的邊緣失真也是最小的.

3.2實(shí)測數(shù)據(jù)實(shí)驗(yàn)

對帶有噪聲的大小是882×882的實(shí)測InSAR相位圖進(jìn)行降噪.ADMM-DL濾波器的參數(shù)設(shè)定與場景一相同.在實(shí)際中，式（1）中的噪聲可以近似成一個(gè)復(fù)的零均值、互相關(guān)矩陣與相關(guān)值γ有關(guān)的白高斯噪聲，即

這意味著在使用ADMM算法式（7）～（10）之前必須用標(biāo)準(zhǔn)差歸一化式（1）中的噪聲變量n（t）.濾波后的圖像如圖2所示.

圖2　實(shí)測結(jié)果對比

表2　不同濾波器的殘點(diǎn)數(shù)比較

采用ADMM-DL濾波器，按照仿真數(shù)據(jù)實(shí)驗(yàn)相同的方式處理數(shù)據(jù).從表2所示的殘點(diǎn)數(shù)和圖2顯示的干涉相位圖可見，ADMM-DL濾波器不但能濾除更多的噪聲而且能更好地保持邊緣完整性.

4 結(jié)束語

筆者提出了一種新的干涉相位濾波法，以字典學(xué)習(xí)為基礎(chǔ)，對復(fù)InSAR數(shù)據(jù)建立了降噪的優(yōu)化模型.基于Lagrangian/Bregman迭代框架，進(jìn)行松弛變換和采用ADMM技術(shù)，使非凸的難解的優(yōu)化問題被分解成幾個(gè)凸的優(yōu)化子問題，并提出了有效的ADMM-DL算法迭代計(jì)算這些子問題.在抑制噪聲和濾除殘點(diǎn)方面，文中方法的性能都優(yōu)于傳統(tǒng)的濾波方法.

［1］SUKSMONO A B，HIROSE A.Interferometric SAR Image Restoration Using Monte Carlo Metropolis Method［J］. IEEE Transactions on Signal Processing，2002，50（2）:290-298.

［2］郭交，李真芳，劉艷陽，等.一種InSAR干涉相位圖的自適應(yīng)濾波算法［J］.西安電子科技大學(xué)學(xué)報(bào)，2011，38（4）:77-81. GUO Jiao，LI Zhenfang，LIU Yanyang，et al.New Adaptive Noise Suppressing Method for Interferometric Phase Images［J］.Journal of Xidian University，2011，38（4）:77-81.

［3］LEE J S，PAPATHANASSIOU K P，AINSWORTH T L，et al.A New Technique for Noise Filtering of SAR Interferogram Phase Images［J］.IEEE Transactions on Geoscience and Remote Sensing，1998，36（5）:1456-1465.

［4］GOLDSTEIN R M，WERNER C L.Radar Interferogram Filtering for Geophysical Applications［J］.Geophysical Reservation Letters，1998，25（21）:4035-4038.

［5］易輝偉，朱建軍，陳建群，等.一種改進(jìn)的InSAR干涉圖復(fù)數(shù)空間自適應(yīng)濾波［J］.中南大學(xué)學(xué)報(bào)，2013，44（2）:632-641. YI Huiwei，ZHU Jianjun，CHEN Jianqun，et al.An Improved Adaptive Algorithm for Filtering InSAR Interferogram in Complex Plane［J］.Journal of Central South University，2013，44（2）:632-641.

［6］李錦偉，李真芳，劉艷陽，等.一種相干系數(shù)加權(quán)的最優(yōu)干涉相位濾波［J］.西安電子科技大學(xué)學(xué)報(bào)，2014，40（2）:25-31. LI Jinwei，LI Zhenfang，LIU Yanyang，et al.Coherence-weighted Optimum Interferometric Phase Filtering Method ［J］.Journal of Xidian University，2014，40（2）:25-31.

［7］AHARON M，ELAD M，BRUCKSTEIN A.K-SVD:an Algorithm for Designing Overcomplete Dictionaries for Sparse Representation［J］.IEEE Transactions on Signal Processing，2006，54（11）:4311-4322.

［8］RAZAVIYAYN M，TSENG H W，LUO Z Q.Dictionary Learning for Sparse Representation:Complexity and Algorithm ［C］//Proceedings of the IEEE International Conference on Acoustics，Speech and Signal Processing.Piscataway:IEEE，2014: 5247-5251.

［9］YU G S，SAPIRO G，MALLAT S.Image Modeling and Enhancement via Structured Sparse Model Selection［C］// Proceedings of the IEEE Conference on Image Processing.Piscataway:IEEE Computer Society，2010:1641-1644.

［10］KERVRANN C，BOULANGER J.Local Adaptively to Variable Smoothness for Examplar-based Image Regularization and Representation［J］.International Journal of Computer Vision，2008，79（1）:45-69.

［11］HAO H，BIOUCAS-DIAS J M，KATKOVNIK V.Interferometric Phase Image Estimation via Sparse Coding in the Complex Domain［J］.IEEE Transactions on Geoscience and Remote Sensing，2015，53（5）:2587-2602.

［12］RAZAVIYAYN M，HONG M Y，LUO Z Q.A Unified Convergence Analysis of Block Successive Minimization Methods for Nonsmooth Optimization［J］.SIAM Journal on Optimization，2013，23（2）:1126-1153.

［13］GRANT M，BOYD S.CVX:Matlab Software for Disciplined Convex Programming［EB/OL］.［2014-07-26］.http:// cvxr.com/cvx/.

［14］HE B S，YUAN X M.On the O（1/n）Convergence Rate of the Douglas-Rachford Alternating Direction Method［J］. SIAM Journal on Numerical Analysis，2012，50（2）:700-709.

［15］BERTSEKAS D P.Nonlinear Programming［M］.2nd Edition.Nashua:Athena-Scientific，1999.

［16］DAUBECHIES I，DEFRISE M，De MOL C.An Iterative Thresholding Algorithm for Linear Inverse Problems with a Sparsity Constraint［J］.Communications on Pure and Applied Mathematics，2004，57（11）:1413-1457.

（編輯:王 瑞）

InSAR noise reduction using adaptive dictionary learning

LUO Xiaomei1，2，SUO Zhiyong3，LIU Qiegen2，4
（1.State Key Lab.of Integrated Service Networks，Xidian Univ.，Xi’an 710071，China；2.Department of Electronic Information Engineering，Nanchang University，Nanchang 330031，China；3.National Key Lab.of Radar Signal Processing，Xidian Univ.，Xi’an 710071，China；4.The Paul C.Lauterbur Research Centre for Biomedical Imaging，Shenzhen Key Lab.for MRI，Chinese Academy of Sciences，Shenzhen 518055，China）

We consider the phase noise filtering problem for interferometric synthetic aperture radar （InSAR）based on the dictionary learning technique.Due to the non-convexity of the optimization problem is difficult to solve.By using the splitting technique and employing the augmented Lagrangian framework，we obtain a relaxed nonlinear constraint optimization problem with l1-norm regularization which can be solved efficiently by the alternating direction method of multipliers（ADMM）.Specifically，we firstly train dictionaries from the InSAR complex phase data，and then reconstruct the desired complex phase image from the sparse representation.Simulation results based on simulated and measured data show that this new InSAR phase noise reduction method has a much better performance than several classical phase filtering methods in terms of residual count，mean square error（MSE）and preservation of the fringe completeness. Key Words: interferometric synthetic aperture radar；phase noise reduction；dictionary learning；l1-norm regularization；alternating directional method of multipliers

TN957.51

A

1001-2400（2016）01-0018-06

10.3969/j.issn.1001-2400.2016.01.004

2014-08-02 網(wǎng)絡(luò)出版時(shí)間：2015-04-14

國家自然科學(xué)基金資助項(xiàng)目（61362001，51165033）

羅曉梅（1971-），女，西安電子科技大學(xué)博士研究生，E-mail:xxmluo@gmail.com.

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/61.1076.TN.20150414.2046.004.html