[發明專利]基于廣義似然比的多航過SAR相干變化檢測方法

有效

• 申請號: 201410384298.0
• 申請日: 2014-08-06
• 公開（公告）號: CN104166128B
• 公開（公告）日: 2016-11-09
• 發明（設計）人: 黃鈺林;王園園;武俊杰;劉曉佳;楊建宇
• 申請（專利權）人: 電子科技大學
• 主分類號: G01S7/41
• 分類號: G01S7/41;G01S13/90
• 代理公司: 成都宏順專利代理事務所(普通合伙) 51227
• 代理人: 周永宏

1.基于廣義似然比的多航過SAR相干變化檢測方法，其特征在于：包括以下步驟：?

S1：選取多航過SAR圖像對：選取不同時間對同一地區多次觀測獲得的K幅SAR圖像并進行配準后記為{f₁,f₂,…,f_K}，K≥3；?

S2：選取多航過SAR圖像像素對，依次選取f₁,f₂,…,f_K對應的復矩陣第m行第n列的值，并記為向量假設向量X_mn～CN(0,Γ)，其中Γ＝E[X_mnX_mn^H]為協方差矩陣，設H₀表示目標區域未發生變化，得到協方差矩陣Γ₀，設H₁表示目標區域發生變化，得到協方差矩陣Γ₁；?

S3：對協方差矩陣進行最大似然估計：分別對協方差矩陣Γ₀和Γ₁進行最大似然估計；?

S4：進行似然比假設檢驗；?

S5：依次選取多航過SAR圖像像素對，重復步驟S4，得到變化檢測結果。?

2.根據權利要求1所述的基于廣義似然比的多航過SAR相干變化檢測方法，其特征在于：所述的步驟S2的具體實現方法為：由于圖像存儲為矩陣形式，因此f₁,f₂,…,f_K均為復圖像，依次選取f₁,f₂,…,f_K對應的復矩陣第m行第n列的值分別記為令復向量m＝1,2,…,M，n＝1,2,…,N，其中[]T表示轉置運算，M、N為圖像對應的復矩陣的大小；?

設向量X_mn服從K維圓對稱復高斯分布，即X_mn～CN(0,Γ)，其中協方差矩陣Γ＝E[X_mnX_mn^H]，因此X_mn的概率密度函數表示為：?

其中，E[X_mnX_mn^H]表示求X_mnX_mn^H的均值，X_mn^H表示求X_mn的復共軛轉置，|Γ|表示Γ的行列式，exp表示指數運算；?

設H₀表示目標區域未發生變化，此時對應的向量X_mn～CN(0,Γ₀)，其中，?

其中，j表示復數單元，ρ_ab≈1，Φ_ab≈0°，a＝1，2，…，K，b＝1，2，…，K，a＜b；?

設H₁表示目標區域發生變化，此時對應的向量X_mn～CN(0,Γ₁)，其中，?

ρ′_ab≈0，Φ′_ab≠0°，由于Γ₀、Γ₁未知，因此為復合假設。?