[发明专利]一种遥感器衰变的定量分析方法及系统

权利要求书

1.一种遥感器衰变的定量分析方法，其特征在于，包括：

将不同时间获得的同一地理区域的卫星遥感数据投影到同一个地理栅格，构建同一地理区域的不同时相的卫星数据集；

计算所述卫星数据集中每个卫星遥感数据所对应的表观反射率；

基于IR-MAD变换，根据预设的不变像元判断阈值，获取所述同一地理区域的不变像元；

基于所述不变像元，对各个光谱通道的所述表观反射率进行正交回归，获取卫星遥感器在各时间段内各个光谱通道的衰变值序列；

根据所述衰变值序列，建立遥感器衰变监测模型。

2.根据权利要求1所述的遥感器衰变的定量分析方法，其特征在于，所述基于IR-MAD变换，根据预设的不变像元的判断阈值，获取所述同一地理区域不变像元，包括：

获取所述卫星数据集中两个不同时刻的卫星遥感数据构建图像对；

根据每个卫星遥感数据所对应的表观反射率，对所述图像对进行IR-MAD变换，获取所述图像对的所有MAD分量；

根据每个所述MAD分量对应的的标准差对所述MAD分量进行归一化；

对所述归一化的MAD分量进行平方求和，根据所述预设的不变像元判断阈值，确定所述图像对上的不变像元。

3.根据权利要求2所述的遥感器衰变的定量分析方法，其特征在于，所述根据所述预设的不变像元判断阈值，确定所述图像对上的不变像元，包括：

其中，N为光谱通道数，MAD_i为第i个光谱通道的MAD分量，为MAD_i的标准差，k为不变像元判断阈值；为自由度为N的卡方分布χ²在概率为prob处的上分位点值。

4.根据权利要求2所述的遥感器衰变的定量分析方法，其特征在于，所述获取所述卫星数据集中两个不同时刻的卫星遥感数据构建图像对，包括：

从所述卫星数据集中选择具有接近相同的地面轨道的两个卫星遥感数据构建图像对；

其中，所述具有接近相同的地面轨道的两个卫星遥感数据为卫星遥感器过境时与赤道交点经度在±1°内获取的。

5.根据权利要求2所述的遥感器衰变的定量分析方法，其特征在于，所述获取所述卫星数据集中任意两个时刻的卫星遥感数据构建图像对，还包括：

所述任意两个时刻的卫星遥感数据的时间间隔为：365.25*y±20天之内，其中，y＝0,1,2,3,…years。

6.根据权利要求2所述的遥感器衰变的定量分析方法，其特征在于，在对所述图像对进行IR-MAD变换，获取所述图像对的所有MAD分量之前，还包括：

对所述图形对中的两个卫星遥感图像进行无效像元掩模处理以及配准处理；

所述对所述图像对中的两个卫星遥感图像进行无效像元掩模处理，包括：

基于阈值判别法，根据所述图像对的各个光谱通道信息，剔除所述图像对中的无效像元；

利用海洋陆地掩模版，剔除所述图像对中的水体及海洋区域的像元；

根据与所述图像对相对应的卫星的轨道高度，剔除所述图像对中观测角超预设范围的像元；

所述对所述图形对中的两个卫星遥感图像进行配准处理，包括：

以所述图像对中的任一卫星遥感图像作为基准图像，在不同方向上移动另一幅卫星遥感图像，并获取各对应位置上的Pearson相关系数；

确定所述Pearson相关系数取最大值时，所述图像对中的两个卫星遥感图像的相对位置为配准位置。

7.根据权利要求2所述的遥感器衰变的定量分析方法，其特征在于，所述计算所述卫星数据集中每个卫星遥感数据所对应的表观反射率，包括：

其中，α_i为第i个光谱通道的DN值的定标斜率，β_i为第i个光谱通道的DN值的截距，其中DN值为卫星遥感数据对应的辐射观测值，d²为日地距离修正因子，θ_s为太阳天顶角，ρ_i为不变像元在第i个光谱通道的表观反射率。