[发明专利]一种两时相遥感影像变化检测方法

权利要求书

1.一种两时相遥感影像变化检测方法，其特征在于，包括：

获取两时相的原始光学遥感影像；

对所述原始光学遥感影像进行影像配准；

利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正；

对辐射归一化校正后的遥感影像，计算每个相应波段的波段差；

对每一个波段差，分别计算扩展属性剖面；

将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌入到FCM模型中，确定影像的变化区域和非变化区域。

2.如权利要求1所述的两时相遥感影像变化检测方法，其特征在于，所述对所述原始光学遥感影像进行影像配准，具体包括：

采用ENVI14.8遥感软件对原始光学遥感影像进行几何粗校正；

利用自动匹配与三角剖分法对几何粗校正后的遥感影像进行几何精校正。

3.如权利要求1所述的两时相遥感影像变化检测方法，其特征在于，所述利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正，具体包括：

将影像配准后的遥感影像的像素灰度值分别调整为从而实现辐射归一化校正，

其中，影像配准后的遥感影像各个波段服从高斯分布，μ_f、σ_f分别为参考影像的像素均值和标准差，g_r、μ_r、σ_r分别为影像配准后的遥感影像的像素灰度值、均值和标准差。

4.如权利要求1所述的两时相遥感影像变化检测方法，其特征在于，所述对每一个波段差，分别计算扩展属性剖面，具体包括：

对每一个波段差D_i，设置其灰度级为L+1，其层级定义为X_l(D_i)＝{p∈D_i|D_i(p)≥l}，对影像D_i在每一灰度级上进行二值化，得到一系列二值图像X_l(D_i)，其中，i＝1,2,…,B，B表示影像的波段数目，p表示像素点，l表示影像的灰度值；

对于D_i上的每一个像素点p，根据不同的属性门限T_λ做属性开变换和属性闭变换，得到p点处的属性开剖面Π((γ^T(D_i))(p))和属性闭剖面其中，

Π((γT(Di))(p))={ΠγTλ:ΠγTλ=(γTλ(Di))(p)},∀λ∈[0,1,...,n]]]>

5.如权利要求4所述的两时相遥感影像变化检测方法，其特征在于，所述将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌入到FCM模型中，确定影像的变化区域和非变化区域，具体包括：

构建FCM结合空间邻域信息的模型为：

Jm=Σj=1CΣk=1Nujkm||D(k)-vj||2+α[Σj=1CΣk=1Nujkm||open(k)-vj||2+Σj=1CΣk=1Nujkm||close(k)-vj||2]]]>

式中，C是聚类数目，C为2，N是样本的总数，表示第k样本对于第j类聚类中心的模糊隶属度，m为隶属度的加权指数，m为2，u_jk∈[0,1]且参数α控制惩罚效果，

D(k)＝[D₁(k)D₂(k)…D_B(k)]，

open(k)＝[Π((γ^T(D₁))(k))Π((γ^T(D₂))(k))…Π((γ^T(D_B))(k))]，

采用式ujk=(||D(k)-vj||2+α(||open(k)-vj||2+||close(k)-vj||2))-1/(m-1)Σj=1C(||D(k)-vj||2+α(||open(k)-vj||2+||close(k)-vj||2))-1/(m-1)]]>更新模糊划分矩阵；

采用式vj=Σk=1Nujkm(D(k)+α(open(k)+close(k)))(1+α)Σk=1Nujkm]]>更新聚类中心；

重复更新模糊划分矩阵和聚类中心，直至相邻两次迭代的聚类中心聚类小于阈值ε；

根据模糊划分矩阵u_jk确定最终的变化检测图，确定影像的变化区域和非变化区域。