[發明專利]一種基于自適應卡方變換的多時相遙感影像變化檢測方法

有效

• 申請號: 201610814220.7
• 申請日: 2016-09-09
• 公開（公告）號: CN106650571B
• 公開（公告）日: 2019-09-10
• 發明（設計）人: 石愛業;儲艷麗
• 申請（專利權）人: 河海大學
• 主分類號: G06K9/00
• 分類號: G06K9/00;G06T7/30;G06T7/10;G06T5/50
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 許方

1.一種基于自適應卡方變換的多時相遙感影像變化檢測方法，其特征在于，包括以下步驟：

(1)輸入兩時相的高分辨率光學遙感影像，分別記為X₁和X₂；

(2)對X₁和X₂進行影像配準；

(3)利用多元變化檢測方法分別對X₁和X₂進行輻射歸一化校正；

(4)計算多時相差異影像D_X＝X₁-X₂；

(5)自適應估計多時相差異影像D_X高斯分布的均值m和方差矩陣Σ；

(6)計算多時相差異影像D_X上每個點的卡方值；

(7)迭代估計最優檢測閾值，利用最優檢測閾值進行卡方變換，得到變化檢測結果。

2.根據權利要求1所述基于自適應卡方變換的多時相遙感影像變化檢測方法，其特征在于，步驟(5)的具體過程如下：

(A)輸入多時相差異影像D_X、至少一個像素點不是變化點的概率p、變化點占整個像素的比值r、使用非變化點的系數s、確定采樣大小的乘積因子f；

(B)計算用來估計均值和方差矩陣的像素點數目：

上式中，N表示D_X的像素點總數，表示取最接近的整數；

(C)計算迭代次數：

上式中，m＝f(1+b)，b是D_X的波段數目；

(D)開始如下迭代過程，直到達到迭代次數為n為止：

a)從D_X中隨機抽取m個點，估計均值m和方差矩陣Σ；

b)對D_X中所有的像素點計算如下的距離：

$d_{ij}=\sqrt{{\left(D_x\right(i,j)-m)}^T{\mathrm{\Σ}}^{-1}\left(D_x\right(i,j)-m)}$

上式中，D_x(i,j)為D_X上坐標(i，j)處的像素；

c)將所有的d_ij按照從小到大的順序排列；

d)從步驟c)的排列結果中取得前N_max個d_ij對應的像素點坐標；

e)計算步驟d)得到的N_max個像素點坐標對應的均值和方差矩陣；

f)計算步驟e)得到的方差矩陣的行列式值；

g)如果迭代次數小于n，回到步驟a)，否則終止迭代，并找出最小的方差距陣行列式值對應的均值和方差矩陣，即為估計的高斯分布的均值和方差矩陣。

3.根據權利要求2所述基于自適應卡方變換的多時相遙感影像變化檢測方法，其特征在于：在步驟(A)中，p＝0.99，r＝0.15，s＝0.95，f＝2。