[发明专利]一种基于辐射校正的长时间序列夜光遥感图像的获取方法

权利要求书

1.一种基于辐射校正的长时间序列夜光遥感图像的获取方法，其特征在于：包括如下步骤：

S1：对现有年份的NPP夜光图做去噪处理，以及临时灯光的去除；

S2：采用改进的脊回归法对DMSP夜光图做相对辐射校正；

S3：利用DMSP和NPP重合年份的年度图像构造出绝对辐射校正参数；

S4：将待校正的DMSP夜光图分别乘以绝对校正参数得到校正后的DMSP的夜光图像，最终获得长时间序列的夜光遥感图像；

所述步骤S3中绝对辐射校正的具体过程为：将与DMSP重合年份的NPP夜光图做重采样，用NPP夜光图分别比上对应年份的DMSP夜光图，得到两幅商图s1和s2，取s1和s2的几何平均，得到绝对辐射校正参数s，即：

2.根据权利要求1所述的一种基于辐射校正的长时间序列夜光遥感图像的获取方法，其特征在于：所述步骤S1中去噪处理具体为：对于现有年份的NPP夜光遥感图像中DN值小于零的值赋值为0。

3.根据权利要求1所述的一种基于辐射校正的长时间序列夜光遥感图像的获取方法，其特征在于：所述步骤S1中临时灯光的去除具体为：选取稳定年份的NPP夜光遥感图像为参考图像，对临时灯光中DN值异常的点赋值为0。

4.根据权利要求1所述的一种基于辐射校正的长时间序列夜光遥感图像的获取方法，其特征在于：所述步骤S2中相对辐射校正的具体过程为：选取参考图像DMSPF152000，遍历整个参考图像，求出每个DN值的像素点个数sum_DN并记录下像素点位置，然后遍历整个待校正的夜光图，记下对应位置的像素点的DN值并统计每个DN值的个数num_DN，接着采用加权求和的方法得到校正后的DN值，即：

在得到64个修正后的DN值后，采用二次曲线拟合的方法，具体公式为f(x)＝ax²+bx+c，得到参考图像到修正图的一个映射P_修正＝f(P_待修正)，其中x是待校正的DN值，f(x)是校正后的DN值，a、b、c表示拟合出的二次函数的系数。