目录

一、对夜间灯光数据进行预处理

<1>图像裁剪(在ENVI处理)

<2>量纲转换(在ENVI中处理)

 <3>去噪处理(这里我们在ENVI中进行中值滤波处理)

 <4>重投影(在Arcmap中进行投影栅格)

 二、提取城市建成区的方法(在Arcmap中操作)

 <1>对上述的数据进行按掩膜提取

<2>转为整型

 <3>栅格转面

 <4>添加“面积”字段

 <5>融合

 <6>表转EXCEL

<7>通过筛选提取城市建成区

---

一、对夜间灯光数据进行预处理

<1>图像裁剪(在ENVI处理)

由于获取到的是全球影像数据，为了节约时间、增加效率，我们首先进行了研究区域的裁剪。

<2>量纲转换(在ENVI中处理)

由于NPP_VIIRS数据单位量级非常小，为了避免数值显示不完全的问题，方便后续的操作，我们利用Band Math对数据进行量纲转换，这里我们将DN值放大10的6次方(根据大家需要)。

 <3>去噪处理(这里我们在ENVI中进行中值滤波处理)

利用中值滤波进行平滑处理。中值滤波会选取数字图像或数字序列中像素点及其周围临近像素点的像素值，将这些像素值排序，然后将位于中间位置的像素值作为当前像素点的像素值，让周围的像素值接近真实值，从而消除孤立的噪声点。通过相关研究, 在进行过滤时选用5×5的中值滤波栅格模块处理结果最优且对城市灯光影响较小。

 <4>重投影(在Arcmap中进行投影栅格)

为使数据更精确，避免网格变形带来的影响以及便于计算影像中灯光亮值像元的面积，通常将夜光遥感影像数据转换成兰伯特方位角等积投影（Asia_Lambert_Conformal_Conic）,该投影比较适合中国地形。同时我们将将栅格大小设置成10×10。

 

 二、提取城市建成区的方法(在Arcmap中操作)

 <1>对上述的数据进行按掩膜提取

得到上海市的夜间灯光数据，必须做这步。

<2>转为整型

通过截断将栅格的每个像元值转换为整型，防止后期报错

 <3>栅格转面

将栅格数据集转换为面要素。字段一般默认为Value字段。

 <4>添加“面积”字段

添加“面积”字段(默认即可)后通过计算几何自动计算出即可，注意单位选择。

 <5>融合

融合字段选择亮度值的字段，统计字段选择面积，统计类型选择SUM。

 <6>表转EXCEL

在EXCEL表中通过对GRIDCODE倒序（从大到小）排列，累加面积大于或等于建成区面积的亮度值即为阈值。表中的阈值即为29824596。

<7>通过筛选提取城市建成区

通过SQL表达式选择大于或等于GRIDCODE的像元即可。

 最终即可提取出城市建成区。

效果图如下：