数据来源：http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html
非辐射定标的夜间灯光影像包括3 种全年平均影像：无云观测频数影像、平均灯光影像和稳定灯光影像。数据缺少星上的辐射校正
8 期辐射定标的夜间灯光影像。与稳定灯光影像相比辐射定标的夜间灯光影像解决了像元DN值饱和的问题（下载）
本文选择了F16 卫星传感器获取的2006 年的影像。
存在2 个问题需在校正过程中解决：
（1）原始影像数据集中的影像是非连续性的；
（2）数据集中的每一期影像都存在着像元DN值饱和的现象。
校正方法：相互校正、饱和校正和影像间的连续性校正
传统方法：一种不变目标区域法的全球范围的相互校正方法，选择意大利西西里岛作为不变目标区域，F121999 的夜间灯光频数影像作为参考影像，再用其他的夜间灯光频数影像分别对参考影像做回归分析得到相对应的二次校正方程。
预处理：
1. 将所有影像的投影坐标系都转换为兰伯特等面积投影坐标系并将影像中的网格重采样为1 km；
2. 论文中通过分析社会经济统计数据，选择黑龙江鹤岗市市辖区作为不变目标区域，DN值从高到底分布广。3）鹤岗市F12 卫星传感器1992 年的夜间灯光影像与F18 卫星传感器2012 年的夜间灯光影像有很好的线性相关关系。
3. F162006 辐射定标的夜间灯光影像的对应区域分别作为待校正影像和参考影像。
相互校正和饱和校正
幂数回归的相关系数在总体上相对最高，说明幂数方程拟合的待校正影像和参考影像的相关性最好，因此，选择幂数方程建立校正模型
影像间的连续性校正
最后插张2013年全国灯光数据
