結合局部頻率估計的小波域InSAR干涉相位濾波方法
【專利摘要】本發(fā)明提供了一種結合局部頻率估計的小波域InSAR干涉相位濾波方法。該小波域InSAR干涉相位濾波方法利用局部頻率估計實現了對復干涉相位的小波系數中有用信息子帶和噪聲子帶的區(qū)分，利用通用閾值收縮和鄰域閾值收縮兩種方法分別具有去噪效果好和細節(jié)保持能力強的特點，對有用信息所在子帶的小波系數進行鄰域閾值收縮，而對其它子帶的小波系數則進行通用閾值收縮，從而盡可能的濾除噪聲，同時保持干涉條紋的細節(jié)信息不被破壞，實現高精度的干涉相位濾波，為高精度的干涉測量提供了條件。
【專利說明】結合局部頻率估計的小波域InSAR干涉相位濾波方法

【技術領域】
[0001] 本發(fā)明涉及電子信息技術雷達【技術領域】，尤其涉及一種結合局部頻率估計的小波 域InSAR干涉相位濾波方法。

【背景技術】
[0002] 干涉合成孔徑雷達（Interferometric Synthetic Aperture Radar，InSAR)是利 用合成孔徑雷達（SAR)兩個通道的干涉相位信息提取地表的高程信息或變化信息，將SAR 的測量拓展到三維空間，具有全天時、全天候、高精度的特點，因此在地形測繪、冰川研究、 海洋測繪以及地面沉降監(jiān)測等多個領域都有廣泛的應用。
[0003] 干涉測量的精度和可靠性在很大程度上取決于干涉相位圖的質量。然而，在實際 系統中，受熱噪聲去相干、時間去相干、基線去相干、配準誤差等多種去相干因素的影響，干 涉相位圖不可避免的存在相位噪聲。低質量的干涉相位將會影響后續(xù)的干涉相位解纏及高 程反演的準確性。因此，在相位解纏前必須對干涉相位進行濾波，從而獲取高質量的干涉相 位圖。
[0004] 目前干涉相位的濾波方法可以大致分為空間域濾波和變換域濾波兩類。圓周期均 值濾波或中值濾波是一種最基本的空間域濾波方法（見參考文獻1)，它實現簡單，但濾波 窗口大小不好確定，在條紋密集時容易破壞相位細節(jié)，降低分辨率。變換域濾波方法在實際 中應用更為廣泛，如Goldstein濾波（見參考文獻2)，然而，該方法受分塊大小和濾波參數 的影響較大，在信噪比很低時，濾波效果較差；小波變換由于其良好的時頻分析特性和多分 辨率特性，也可以應用于干涉相位濾波中（見參考文獻3)，該方法通過增大小波系數中的 信號成分來實現濾波，能夠很好地保持干涉條紋的細節(jié)信息，并在一定程度上提高了圖像 的信噪比，但由于對噪聲信息沒有進行有效的抑制，使得去噪效果較差。因此，為滿足InSAR 應用對干涉相位精度的要求，有必要進一步研究能有效去除噪聲并保持相位細節(jié)的變換域 濾波方法。
[0005] 參考文獻：
[0006] [1]R. Lanari. Generation of digital elevation models by using SIR-C/X-SAR multifrequency two-pass interferometry ：The Etna case study. IEEE Transactions on Geoscience and Remote Sensing,1996,34(5) :1097-1114.
[0007] [2]R. M. Goldstein, C. L. Werner. Radar Interferogram filtering for Geophysical Application. Geophysical Research Letters. 1998,25(21) :4035-4038.
[0008] [3] Lopez-Martinez C, Fabregas X. Modeling and reduction of SAR interferometric phase noise in the wavelet domain[J]. IEEE Trans, on Geoscience and Remote Sensing,2002,40(12) :2553-2566


【發(fā)明內容】

[0009] (一）要解決的技術問題
[0010] 為解決上述問題，本發(fā)明提供了一種結合局部頻率估計的小波域InSAR干涉相位 濾波的方法，以克服現有的小波域干涉相位濾波方法不能兼顧去噪和細節(jié)保持的缺點，從 而提高干涉相位的精度。
[0011](二）技術方案
[0012] 根據本發(fā)明的一個方面，提供了一種結合局部頻率估計的小波域InSAR干涉相位 濾波方法。該小波域InSAR干涉相位濾波方法包括：步驟A :將InSAR的干涉相位Φ變換 到復數域&_Φ，分別取復數域干涉相位&_φ的實部和虛部；步驟B :對干涉相位Φ進行局部 頻率估計，得到干涉相位所在的頻率范圍；步驟C :對復數域干涉相位的實部和虛部，分別 進行尺度為s的小波分解，得到不同子帶的小波系數的頻率范圍，其中，m， n為小波系數 的位置，i為分解的尺度，其范圍為1?s ;步驟D :對復數域干涉相位的實部和虛部，分別根 據干涉相位所在的頻率范圍和不同子帶小波系數的頻率范圍，分別確定有用信息和噪聲所 在的子帶；步驟E :對復數域干涉相位的實部和虛部，分別對噪聲所在子帶的小波系數分別 進行通用閾值收縮處理；步驟F :對復數域干涉相位的實部和虛部，分別對有用信息所在子 帶的小波系數分別進行鄰域閾值收縮處理；步驟G :對復數域干涉相位的實部和虛部，分別 將通用閾值收縮處理處理后的噪聲所在子帶的小波系數和鄰域閾值收縮處理后的有用信 息所在子帶的小波系數共同進行小波重構，得到濾波后的復數域干涉相位的實部和虛部； 以及步驟H :由濾波后的復數域干涉相位的實部和虛部得到InSAR濾波后的干涉相位。
[0013] (三）有益效果
[0014] 本發(fā)明利用局部頻率估計實現了對復干涉相位的小波系數中有用信息子帶和噪 聲子帶的區(qū)分，利用通用閾值收縮和鄰域閾值收縮兩種方法分別具有去噪效果好和細節(jié)保 持能力強的特點，對有用信息所在子帶的小波系數進行鄰域閾值收縮，而對其它子帶的小 波系數則進行通用閾值收縮，從而盡可能的濾除噪聲，同時保持干涉條紋的細節(jié)信息不被 破壞，實現高精度的干涉相位濾波，為高精度的干涉測量提供了條件。

【專利附圖】

【附圖說明】
[0015] 圖1為根據本發(fā)明實施例結合局部頻率估計的小波域InSAR干涉相位濾波方法的 流程圖；
[0016] 圖2為意大利Etna火山實測的星載InSAR干涉相位圖；
[0017] 圖3A?圖3D為分別利用圓周期均值濾波,Goldstein濾波，Lopez-Martinez C提 出的小波濾波以及本實施例方法濾波后的干涉相位圖

【具體實施方式】
[0018] 為使本發(fā)明的目的、技術方案和優(yōu)點更加清楚明白，以下結合具體實施例，并參照 附圖，對本發(fā)明進一步詳細說明。需要說明的是，在附圖或說明書描述中，相似或相同的部 分都使用相同的圖號。附圖中未繪示或描述的實現方式，為所屬【技術領域】中普通技術人員 所知的形式。另外，雖然本文可提供包含特定值的參數的示范，但應了解，參數無需確切等 于相應的值，而是可在可接受的誤差容限或設計約束內近似于相應的值。
[0019] 本發(fā)明將將干涉相位變換到小波域后，利用局部頻率估計判斷出小波系數中包含 有用信息的子帶，對有用信息所在子帶的小波系數利用鄰域閾值進行收縮處理，而對其他 子帶的小波系數利用通用閾值進行收縮處理，從而盡可能的濾除噪聲，同時不損失干涉條 紋的細節(jié)信息，從而為InSAR干涉相位濾波提供了一種有效的方法。
[0020] 根據本發(fā)明的一個方面，提供了一種結合局部頻率估計的小波域InSAR干涉相位 濾波方法。圖1為根據本發(fā)明實施例結合局部頻率估計的小波域InSAR干涉相位濾波方法 的流程圖。如圖1所示，本實施例結合局部頻率估計的小波域InSAR干涉相位濾波方法包 括：
[0021] 步驟A :將InSAR預處理后得到的干涉相位Φ變換到復數域e# = cos<i)+jsin<i) 中，分別取復數域干涉相位6']<*>的實部1^{6' ]<*>}=〇08<1)，虛部11]1{6']<* >}=8；[11?。?br> [0022] 此處的預處理包括：雙通道InSAR數據的成像處理，配準處理、復數據的共軛相乘 等等。這些均為本領域內的公知常識，此處不再詳細說明。
[0023] 步驟B :對復數域干涉相位&_Φ進行局部頻率估計，得到干涉相位所在的頻率范 圍；
[0024] 該步驟B具體包括：
[0025] 子步驟Bl :對干涉相位Φ的每一像素，對其進行局部頻率估計，得到其方位向干 涉相位頻率/3和距離向干涉相位頻率，具體如下：
[0026] 子分步驟Bla :以該像素為中心取（2M+1) X (2N+1)的估計窗口，該估計窗口的相 位模型可表示為：
[0027]

【權利要求】
1. 一種結合局部頻率估計的小波域InSAR干涉相位濾波方法，其特征在于，包括： 步驟A:將InSAR的干涉相位Φ變換到復數域&_φ，分別取復數域干涉相位^1*·的實部 和虛部； 步驟B:對干涉相位Φ進行局部頻率估計，得到干涉相位所在的頻率范圍； 步驟C:對復數域干涉相位的實部和虛部，分別進行尺度為s的小波分解，得到不同子 帶的小波系數4,,的頻率范圍，其中，m，n為小波系數的位置，i為分解的尺度，其范圍為1? s； 步驟D:對復數域干涉相位的實部和虛部，分別根據干涉相位所在的頻率范圍和不同 子帶小波系數的頻率范圍，分別確定有用信息和噪聲所在的子帶； 步驟E:對復數域干涉相位的實部和虛部，分別對噪聲所在子帶的小波系數分別進行 通用閾值收縮處理； 步驟F:對復數域干涉相位的實部和虛部，分別對有用信息所在子帶的小波系數分別 進行鄰域閾值收縮處理； 步驟G:對復數域干涉相位的實部和虛部，分別將通用閾值收縮處理處理后的噪聲所 在子帶的小波系數和鄰域閾值收縮處理后的有用信息所在子帶的小波系數共同進行小波 重構，得到濾波后的復數域干涉相位的實部和虛部；以及 步驟H:由濾波后的復數域干涉相位的實部和虛部得到InSAR濾波后的干涉相位。
2. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟B包 括： 子步驟Bl:對干涉相位Φ的每一像素，對其進行局部頻率估計，得到其方位向干涉相 位頻率尤和距離向干涉相位頻率X; 子步驟B2:對所有像素估計得到的方位向干涉相位頻率分別取最大值和最小值為 .元和，距尚向干涉相位頻率/?分別取最大值和最小值為和，得到整幅干涉相 位圖所在的頻率范圍為i.L,m..Lml·ι./;_.I*
3. 根據權利要求2所述的小波域InSAR干涉相位濾波方法，其特征在于，所述子步驟 Bl中，對干涉相位Φ的每一像素得到其方位向干涉相位頻率￡和距離向干涉相位頻率X 包括： 子分步驟Bla:以該像素為中心?。?M+1)X(2N+1)的估計窗口，該估計窗口的相位模 型為：
其中，爐(w,?)為估計窗口中心點的相位，k、1為窗口內像素相對窗口中心的位移，fa，f； 分別為該窗口沿方位向和距離向的干涉相位頻率； 子分步驟Blb:通過最大化估計窗口的相位模型的代價函數來估計該像素沿方位向的 干涉相位頻率Z和距離向的干涉相位頻率X:

4. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟C 中，不同子帶小波系數的頻率范圍為：
其中，< 為子帶對應干涉相位圖的低頻部分，rff, 和分別為對應干涉相位圖 在各個尺度的垂直方向、水平方向和對角線方向的高頻部分。
5. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟C 中：用于小波分解的函數為Symlets函數、Daubechies函數或Coiflets函數，所述尺度s 介于2至8之間。
6. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟D根 據干涉相位所在的頻率范圍和不同子帶小波系數的頻率范圍，分別確定有用信息和噪聲所 在的子帶包括： 分別判斷不同子帶小波系數的頻率范圍af? , 與 iUjxlUj是否有交集，如果有交集，則判斷該子帶為有用信息所在的子帶；如果 交集為空，判斷該子帶是噪聲所在子帶。
7. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟E 中，按照下式對噪聲所在子帶的小波系數進行通用閾值收縮處理：
其中，為閾值收縮后的小波系數，sgn(·)為符號函數，下標+表示保持正值不變， 將負值置零；T為通用閾值。
8. 根據權利要求1所述的小波域InSAR干涉相位濾波方法，其特征在于，所述步驟F中 對有用信息所在子帶的小波系數分別進行鄰域閾值收縮處理包括： 以當前要處理的小波系數#為中心，選取大小合適的窗口W，令Σ 則閾值收縮后的小波系數P為：c,,=4,,，其中，下標+表示保持正值不變， 將負值置零，T為通用閾值。
9. 根據權利要求7或8所述的小波域InSAR干涉相位濾波方法，其特征在于，所述通用 閾值T按照下式計算： T =^2tog￥ 其中，N為干涉相位圖包含的像素數，σ為噪聲標準差，按下式進行估計： SJedian(yj) σ =-! L-〇 0.6745
10.根據權利要求1至8中任一項所述的小波域InSAR干涉相位濾波方法，其特征在 于，所述步驟H中，按照下式得到InSAR濾波后的干涉相位"
其中，&和?分別為濾波后的復數域干涉相位的實部和虛部。
【文檔編號】G01S7/02GK104459633SQ201410717610
【公開日】2015年3月25日 申請日期:2014年12月1日 優(yōu)先權日:2014年12月1日
【發(fā)明者】李芳芳, 丁赤飚, 雷斌, 林雪, 胡東輝, 仇曉蘭 申請人:中國科學院電子學研究所