1、准备遥感数据：下载的是2017年月度数据，共12组。

 

         下载地址：科罗拉多矿业大学地球观测组(EOG)  网址：VIIRS Nighttime Light (mines.edu)https://eogdata.mines.edu/products/vnl/

2、准备矢量数据：使用的是南京主城区

        下载地址：1）根据所需下载 .json格式DataV.GeoAtlas地理小工具系列 (aliyun.com)http://datav.aliyun.com/portal/school/atlas/area_selector

        2）转化为shpmapshaperhttps://mapshaper.org/

        3)使用工具箱：合并， 将各个区域合并成一个shp

3、按掩膜提取：裁剪夜光数据，(最好保存在数据库中.gbd ，可以避免数据存储不了问题)

4、重投影：投影栅格，将地理坐标系转化为投影坐标系（此处选择的是Albers等面积投影），可以为重采样做铺垫，方便改空间分辨率。

5、重采样：将空间分辨率改为500*500 ， 重复4-5步骤处理完12个月的数据

6、栅格计算器： 遥感数据叠加获取平均值

7、栅格转点：插值需要转化为点

 

8、插值分析：作用是补充未采集到的缺失点数据；主意Z值的选择，此处选择的是IDW反距离权重法

 

10、创建渔网：像元宽高0.005表示分辨率500*500，几何类型选择面PLOYGON

11、按位置选择：裁剪渔网.面

        11.1 发现裁剪无效

 

         11.2 换种方法裁剪，按位置选择

12、裁剪：裁剪渔网点point

 

13、多值提取至点：经过插值的像素点值提取到了栅格点

 

 14、连接：面关联点，提取值

 

 15、符号化