[发明专利]一种基于稳定映射点夜间灯光影像间相互辐射校正方法

权利要求书

1.一种基于稳定映射点夜间灯光影像间相互辐射校正方法，其特征在于包括以下步骤：

1)读取待相互辐射校正的夜间灯光影像；

2)根据读入的夜间灯光影像确定参考地区；

3)确定校正参考影像；

4)提取影像中的稳定映射点作为相互校正中模型参数的回归解算点；

5)将回归解算点中的参考影像DN值与各待校正影像映射的DN值之间建立有理函数校正模型，计算得出校正参数，其中，有理函数校正模型为：

式中，DN_(C，i)为稳定映射点的参考影像DN值数据，DN_(n，i)为待校正第n年影像稳定映射点对应的DN值数据，p1、p2、p3、q1、q2为待计算的校正参数；

6)采用步骤5)得出的校正参数对各待校正影像进行校正。

2.如权利要求1所述的一种基于稳定映射点夜间灯光影像间相互辐射校正方法，其特征在于，步骤4)提取影像中的稳定映射点作为相互校正中模型参数的回归解算点，具体过程为：

4.1)统计每幅待校正影像与参考影像对应像元DN值的映射数目；

4.2)选取每一幅待校正影像与参考影像间的稳定映射点作为该待校正影像校正参数的回归解算点，具体选取步骤如下：

4.21)从参考影像DN值映射到待校正影像的映射函数表示为：

DN_C＝F_(n，C)(DN_n)+δ_(n，C)，

式中，DN_C为待校正影像中的DN值数据，DN_n为待校正影像第n年影像对应的DN值数据，F_(n，C)为待校正影像对应参考影像的映射关系函数，δ_(n，C)为该映射关系下的噪声及误差，分别选取参考影像DN值从1～62中映射数目最多的点作为稳定映射点的候选点；

4.22)从校正影像DN值映射到参考影像的映射函数表示为：

DN_n＝F_(C，n)(DN_C)+δ_(C，n)，

式中，F_(C，n)为参考影像对应待校正影像的映射关系函数，δ_(C，n)为该映射关系下的噪声及误差，分别选取校正影像DN值从1～62中映射数目最多的点作为稳定映射点的候选点；

4.23)将上述两个稳定映射点的候选点中映射关系一致的点作为最终的稳定映射点，并将最终的稳定映射点作为回归解算点。

3.如权利要求1或2所述的一种基于稳定映射点夜间灯光影像间相互辐射校正方法，其特征在于，6)采用步骤5)得出的校正参数对各待校正影像进行校正具体为：

6.1)根据有理函数模型对不同传感器不同时间的待校正影像进行校正：

式中，DN′_(n，i)为校正后的第n年第i像元的DN值结果，DN_(n，i)为待校正的第n年影像中第i像元的DN值；

6.2)对不同传感器获取的同一年影像进行连续性校正：

将校正后的不同传感器获取的同一年影像对应像元取平均值得到最终校正后的结果：

式中，DN″_(n，i)为DMSP/OLS非辐射定标第n年影像中第i个像元的最终校正结果，和分别为模型校正后的第n年两个不同传感器第i像元的DN值。