[发明专利]一种基于自适应卡方变换的多时相遥感影像变化检测方法

权利要求书

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。