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、符号化