[發明專利]基于相關加權聯合子空間投影的InSAR干涉相位估計方法

無效

• 申請號: 201210003415.5
• 申請日: 2012-01-09
• 公開（公告）號: CN102565798A
• 公開（公告）日: 2012-07-11
• 發明（設計）人: 李海;吳仁彪;賈瓊瓊
• 申請（專利權）人: 中國民航大學
• 主分類號: G01S13/90
• 分類號: G01S13/90;G01S7/36
• 代理公司: 天津才智專利商標代理有限公司 12108
• 代理人: 龐學欣

1.一種基于相關加權聯合子空間投影的InSAR干涉相位估計方法，其特征在于：所述的InSAR干涉相位估計方法包括按順序進行的下列步驟：

1)用傳統的相關法以主圖像作為參考對副圖像進行粗配準，得到粗配準圖像；

2)利用主圖像和粗配準圖像構造相關加權觀測矢量；

3)計算相關加權聯合協方差矩陣和加權聯合相關函數矩陣；

4)計算噪聲子空間和信號子空間；

5)用噪聲子空間投影來構造代價函數，并通過代價函數估計干涉相位。

2.根據權利要求1所述的基于相關加權聯合子空間投影的InSAR干涉相位估計方法，其特征在于：在步驟2)中，所述的構造相關加權觀測矢量的方法是：

a.分別計算主圖像中像素a與粗配準圖像中像素m的相關系數：

r31(m,a)=|Σk=-KKs3(m+k)s1*(a+k)|Σk=-KK|s3(m+k)|2Σk=-KK|s1(i+k)|2m=a-5,a-4,a-3,a-1,a,a+1,a+3,a+4,a+5a=i-1,i,i+3,i+4]]>

b.用所得到的相關系數r₃₁(m，a)分別對粗配準圖像中的像素m進行加權，得到元素s_3W(a，w)，公式如下：

s3W(a,w)=waTVa]]>

其中：

V_a＝[s₂(a-5)，s₂(a-4)，s₂(a-3)，s₂(a-1)，s₂(a)，

????????s₂(a+1)，s₂(a+3)，s₂(a+4)，s₂(a+5)]^T

??????????????????????????(a＝i-1，i，i+3，i+4)

w_a＝[r₃₁(a-5，a)，r₃₁(a-4，a)，r₃₁(a-3，a)，r₃₁(a-1，a)，

???????r₃₁(a，a)，r₃₁(a+1，a)，r₃₁(a+3，a)，r₃₁(a+4，a)，r₃₁(a+5，a)]^T

??????????????????????????(a＝i-1，i，i+3，i+4)

c.用主圖像中的像素a(a＝i-1，i，i+3，i+4)和得到的元素s_3W(a，w)構造相關加權觀測矢量：

js(i,w)=[s1(i-1),s2W(i-1),s1(i),s2W(i),s1(i+3),s2W(i+3),s1(i+4),s2W(i+4)]T]]>

=[s1(i-1),wi-1TVi-1,s1(i),wiTVi,s1(i+3),wi+3TVi+3,s1(i+4),wi+4TVi+4]T]]>

其中：s₁和s₃分別為主圖像和粗配準圖像所對應的復圖像數據，

i表示當前要計算的主圖像和粗配準圖像的像素標號，

2K+1為從相鄰像素中獲得的獨立同分布樣本數，

上標T表示轉置操作，上標*表示共軛操作。