[发明专利]一种基于光谱斜率差异检测地表覆盖变化的方法

权利要求书

1.一种基于光谱斜率差异检测地表覆盖变化的方法，所述方法用于对同一地点的两个不同时刻的Landsat图像进行地表覆盖变化检测，从而测定所述地区发生地表覆盖变化的区域范围，所述两个不同时刻分别为T1、T2时刻，其特征在于，所述方法包括如下步骤，

A、将同一地点的Landsat图像的六个波段的光谱数据依据波谱连接成五个光谱段，分别针对所述T1、T2时刻的Landsat图像的每个点位，计算每个光谱段的斜率，光谱段（n,n+1）的斜率计算公式如下：

k(n,n+1)=ΔrefΔλ=refn+1-refλn+1′-λn′]]>

其中，Δref是波段n+1和波段n的反射率差值，Δλ是波段n+1和波段n的波段长度之差，λ′_n是归一化后无量纲的第n波段波长，λ′_n的计算公式如下：

λn′=λn-min(λ)max(λ)-min(λ)]]>

其中，λ_n是第n波段的波长，min(λ)是所述Landsat图像的所有波段中波长的最小值，max(λ)是所述Landsat图像的所有波段中波长的最大值，λ_n、min(λ)、max(λ)的单位均为微米；

B、根据步骤A中的斜率计算公式，分别计算所述T1、T2时刻的Landsat图像的每个点位的各波段的斜率，形成斜率向量，T1时刻的斜率向量是K_i=(_ki_(1,2),k_i(2,3),...，k_i(n，n+1))^T，T2时刻的斜率向量是k_j=(k_j(1，2),k_k(2,3)，...,k_j(，n+1))^T，

通过斜率向量的绝对距离计算斜率变化强度，计算公式如下：

dK=Σn=15|Δk(n,n+1)|=Σn=15|ki(n,n+1)-kj(n,n+1)|=Σn=15|refi,n+1-refi,nλn+1′-λn′-redj,n+1-refj,nλn+1′-λn′|]]>

=Σn=15|(refi,n+1-refj,n+1)-(refi,n-refj,n)λn+1′-λn′|]]>

其中，ref_i，n+1和ref_i,n是T1时刻第n+1波段和第n波段的反射率，ref_j，n+1和ref_j,n是T2时刻第n+1波段和第n波段的反射率，对应波段的归一化波长分别是λ_n+1和λ_n；d_K是两斜率向量的变化强度；

C、根据步骤B中计算出的所述T1、T2时刻的Landsat图像的每个点位的斜率变化强度计算结果d_K，通过设置阈值确定所述Landsat图像中的地表覆盖变化区域和地表覆盖不变区域。

2.根据权利要求1所述的方法，其特征在于，所述方法进一步包括如下步骤D，用于确定地表覆盖变化区域的地物变化类别，所述步骤D为，根据已知的各种地物特征点及其所对应变化类型的斜率变化特征数据，利用地表覆盖变化区域的像素和上述已知的变化类型斜率变化特征数据的比对，确定出所述像素的变化类型，从而确定出步骤C中所述Landsat图像中的地表覆盖变化区域的每个像素的地物变化类型。