学术咨询服务正当时学报期刊咨询网是专业的学术咨询服务平台！
发布时间：2021-06-02 16:57所属平台：学报论文发表咨询网浏览： 次
摘 要：基于NPP/VIIRS夜间灯光遥感数据，采用阈值法结合统计数据提取了城市建成区范围，并通过图斑综合与边界平滑等方法对城市建成区进行了筛选和边界优化。将NPP/VIIRS数据分别与POI核密度、不透水面产品进行耦合性分析发现，虽然NPP/VIIRS空间分布与POI核
摘 要：基于NPP/VIIRS夜间灯光遥感数据，采用阈值法结合统计数据提取了城市建成区范围，并通过图斑综合与边界平滑等方法对城市建成区进行了筛选和边界优化。将NPP/VIIRS数据分别与POI核密度、不透水面产品进行耦合性分析发现，虽然NPP/VIIRS空间分布与POI核密度、不透水面具有一定重叠度，且变化趋势一致性较高，但整体相关性不高。选取扩展速率、扩展强度、紧凑度、分形维数、用地增长弹性系数和标准差椭圆等指标对城市空间结构和变化趋势等进行了定量研究。结果表明，重庆市2013-2019年城市扩张趋势明显，扩展速率和扩展强度呈年际波动趋势;城市以外部扩张为主，重心整体向东南方向偏移;紧凑度水平和不规则程度整体较低，且呈逐年降低的趋势;可能存在人口增长与用地比例失调的情况，但变化趋势逐渐变好。通过对城市空间格局变化的分析，以期能为未来城市建设布局提供一定参考。
关键词：NPP/VIIRS数据;城市建成区;空间发展;重庆市
城市建成区是指实际已成片开发建设、市政公用设施和公共设施基本具备的地区，作为衡量城市发展的重要指标，准确掌握其现状和动态变化特征能优化完善城市空间布局，指导未来城市建设[1]。夜间灯光遥感数据可表达地表在夜间的亮度辐射状况，直观反映城市空间的发展状况[2]。国内外许多学者将夜间灯光遥感数据应用于城市建成区提取、不透水面提取和人口空间化等方面，取得了不错效果[3-8]。在夜间灯光遥感数据中，NPP/VIIRS较DMSP/OLS具有更高的空间分辨率和辐射分辨率，拓展了夜间灯光的应用领域[9]。
城市建设论文范例：中国城市山水人文空间格局研究
利用夜间灯光遥感数据提取城市建成区边界的方法主要包括阈值法和面向对象分类法，但由于数据分辨率限制，需要融合其他高分辨率数据，才能在一定程度上提高城市建成区的提取精度[10-16]。例如，POI数据与夜间灯光数据之间具有相关性，可用于优化城市建成区边界;夜间灯光数据与不透水面范围具有较强相关性，可为城市建成区提取提供验证参考[17-18]。
基于2013-2019年NPP/VIIRS数据，本文首先进行坐标转换、数据裁切、重采样、噪声去除等预处理，然后利用以社会经济数据为参考依据的阈值法提取城市建成区边界范围，再通过破碎图斑综合、边界平滑等方法进行优化，最后形成城市建成区范围成果。本文将POI核密度分析和不透水面数据分别与NPP/VIIRS数据进行耦合性分析，并选取扩展速率、扩展强度、紧凑度、分形维数、用地增长弹性系数和标准差椭圆等指标对重庆市城市建设的空间格局、发展变化趋势等进行定量分析，以期为未来国土空间城市建设布局提供参考。
1研究区概况与数据来源
重庆市位于四川盆地东部，市域面积为8.24万km2，下辖38个区县(26个区、8个县、4个自治县)。根据国土空间规划功能，重庆市可分为中心城区、主城新区、渝东北地区和渝东南地区。2019年全市地区生产总值为23605.77亿元，城镇常住人口为2086.99万人(经济、人口数据来源于《2019年重庆市国民经济和社会发展统计公报》)。本文采用的NPP/VIIRS数据为该月度合成数据产品(https://eogdata.mines.edu/download_dnb_composites.html)，并校正了受天气影响的像元辐射亮度值[7]，数据时间为2013年1月-2019年12月。
本文采用的其他数据包括高分辨率遥感影像数据(WorldView-3，分辨率为0.5m，时间为2019年)、社会经济统计数据、行政区划数据、河流数据、POI数据、不透水面数据和规划相关数据，其中行政区划数据(重庆市界、区县界、区县驻地、乡镇街道界等)、河流数据来源于2019年基础性地理国情监测成果;社会经济统计数据(常住城镇人口、城市建成区等)来源于《重庆市统计年鉴》(https://www.yearbookchina.com/navibookYCQTJ.html)，数据时间为2013-2019年;POI数据为通过网络抓取的高德地图，数据时间为2018年8月;不透水面数据分辨率为30m，现势性为2015年，来源于国家科技基础条件平台—国家地球系统科学数据中心(http://www.geodata.cn);规划相关数据主要用于城市建成区范围参考。
2城市建成区提取
2.1数据预处理
NPP/VIIRS数据预处理包括坐标转换、数据裁切、重采样、噪声去除和对数变换。首先将数据坐标系(WGS84坐标系)统一转换为2000 国家大地坐标系，再利用重庆市界对数据范围进行裁切，然后将数据重采样为0.5km×0.5km的网格。噪声去除包括海洋背景噪声去除、离群辐射值校准和水域赋值噪声去除，其中海洋背景噪声去除选取157.27°～178.33°E，20.89°～337.10°N范围内每月的NPP/VIIRS数据，统计其平均值作为背景噪声值，再将小于该平均值的数据赋值为零[5];离群辐射值校准根据辐射像元值在前一个月大于0则在后一个月也应大于0的原则，将不符合该原则的像元作为离群像元去除;水域赋值噪声去除通常采用高分辨率遥感影像对水域范围进行勾画，再将该范围内的辐射值赋值为0，本文将2019年基础性地理国情监测水域范围数据作为去噪范围。
2.2城市建成区提取
本文采用阈值法提取城市建成区，即假定一个初始亮度阈值对城市建成区进行提取，再不断改变亮度阈值，将与统计数据最接近的阈值作为城市建成区提取标准。
将预处理后的NPP/VIIRS数据按照0.01的间隔进行离散分级，提取城市建成区边界。由于城市建成区是具有一定规模已成片开发建设的区域，需对分散的零碎区域进行剔除。考虑到重庆是典型的山地城市，受地形、地貌与河流等影响，城市以组团型发展为主[20]，整体建设空间较小，因此需要设定合适的规模阈值提取城市建成区。重庆市的38个区县中，城口县位于大巴山区，城市建成区规模最小。
由《重庆市城口县建设用地控制性详细规划方案》可知，规划至2020年中心城区城市建设用地面积为4.0km2，因此结合高分辨率遥感影像，对大于1km2的建设区域进行保留。由于部分建制镇建设面积较大，在进行城市建成区提取时，若以1km2作为城市建成区的提取标准，江津区白沙镇、武隆区平桥镇、石柱县黄水镇、垫江县澄溪镇等部分建制镇将纳入城市建成区范围，因此以各区县相关规划确定的乡镇街道界为范围，将建制镇剔除。按照上述方法，将城市建成区提取结果与统计数据中规模最接近的亮度阈值作为该年度城市建成区的亮度值，以该阈值为该年度城市建成区规模阈值;同时将栅格数据转换为矢量成果，以1km的宽度对数据进行平滑处理，消除边缘锯齿化，得到城市建成区最终拟合成果。
3城市建成区耦合性分析
3.1与POI数据的耦合性分析
1)利用POI提取城市建成区范围。POI数据包括建筑物、基础设施、公共服务设施等多类信息，在城市地区具有集聚效应和规模效益特征，已广泛应用于城市建成区提取，且提取的城市建成区在边缘处具有平滑特征[12-13,21]。本文首先对POI数据进行空间匹配、数据去重、坐标转换等预处理，再通过反复试验选择适合于重庆市特殊地形的搜索半径，最终以500m为搜索半径进行核密度计算(栅格单元为50m)，并基于ArcGIS平台绘制核密度等值线图。
在城市建成区范围内，POI均值随着NPP/VIIRS均值的增加而具有逐渐增加的趋势，但二者的相关性不明显，主要原因在于：它们对城市建成区提取依据的重点特征不同，POI核密度分析提取的城市建成区图斑较零散，且在城市边缘处较平滑;而由于在城市中心区域有灯光溢出现象，NPP/VIIRS数据对于城市建成区的提取表现为集中连片特征，但在城市边缘区域由于分辨率的影响呈现锯齿状特征。存在较多POI低值而NPP/VIIRS值较高的现象，这主要是由于POI数据提取的城市建成区在城市内部存在部分漏洞区域，而NPP/VIIRS数据提取成果在城市内部漏洞区域较少所致。
3.2与不透水面数据的耦合性分析
不透水面数据主要是基于多源遥感数据产品(GlobeLand30地表覆盖产品、VIIRS夜间灯光数据和MODISEVI植被指数产品)，利用随机森林分类等模型方法形成的综合拼接产品(由中国科学院空天信息创新研究院生产，分辨率为30m)。该产品能较为准确地反映人工不透水面的地表空间分布情况[22]。2015年重庆市不透水面的总面积为1479.82km2，而城市规划涉及乡镇街道内不透水面面积为1220.31km2，占NPP/VIIRS数据提取的城市建成区的84.80%。
从空间上来看，NPP/VIIRS数据提取的城市建成区范围内的不透水面面积为719.15km2，占比为49.97%，与POI数据和NPP/VIIRS数据重叠度相比，空间重叠度较低。为实现不透水面与NPP/VIIRS数据的比较分析，本文以城市规划涉及乡镇街道为范围，构建1km×1km的正方形网格，统计网格内不透水面均值和NPP/VIIRS均值，并分析二者的相关性。结果表明，不透水面均值与NPP/VIIRS灯光亮度存在一定的相关性，随着不透水面均值的增加，亮度逐渐增加;但相关性不明显，相关系数仅为0.4272，主要是由于不透水面数据分辨率较高，分布也较离散，城市内部灯光亮度值较高的区域没有不透水面数据分布。
4城市空间演变分析
1)变化特征分析。本文通过扩展速率、扩展强度两个指标对重庆市2013-2019年城市建成区总体变化特征进行分析。其中，扩展速率为相邻两个年度城市建成区的面积之差;扩展强度为相邻两个年度城市建成区面积之差与前一年度的比值。重庆市2013-2019年城市建成区总体扩张面积为373.99km2，年均扩张面积为62.33km2，其中2017-2018年城市建成区扩张面积最大，2013-2014年扩展强度最大。扩展速率和扩展强度均具有波动特征，尤其是2018-2019年扩展速率和扩展强度均为最小，原因可能为：①受全球经济形势整体下滑影响，在城市建设方面有所减少;②由于各项政策影响，特别是房地产和工业用地限制，土地交易市场出现一定降温。
2013-2019年渝中区、大渡口区、江北区、沙坪坝区、九龙坡区、南岸区、北碚区、渝北区等中心城区的城市建成区扩张面积最大，年均变化面积为28.07km2，接近全市扩张总面积的一半;主城新区其次，年均城市建成区面积增长22.81km2;渝东北地区城市建成区扩张以万州区为主要增长点;渝东南面积扩张面积最少，其中黔江区和武隆区年均扩张面积超过1km2。
结果表明，重庆市城市发展仍以中心城区为核心区域，主城新区受中心城区辐射带动作用，城市发展速度较快，中心城区和主城新区的城市建成区扩张占比较高，渝东北地区和渝东南地区为四川盆地盆周区域，山地丘陵较多，受地形地貌限制和功能区定位影响，城市建设面积较小。中心城区城市扩张以渝北区为主要增长极，主要受国家级新区(重庆两江新区)产业发展驱动，近年来渝北区城市扩展主要区域为东部龙兴工业园区。在全市38个区县中，年均城市建成区增长面积超过1km2有21个区县，超过2km2有12个区县，超过3km2有5个区县，分别是渝北区、巴南区、江北区、江津区和永川区。
5结 语
本文利用NPP/VIIRS夜间灯光遥感数据，采用阈值法等方法对城市建成区范围进行了提取与优化;并结合POI核密度分析、不透水面数据，分别分析了二者与NPP/VIIRS数据的耦合性，选取扩展速率、扩展强度、紧凑度、分形维数、用地增长弹性系数和标准差椭圆等指标对城市空间结构及其变化趋势进行了定量研究。
1)POI核密度随NPP/VIIRS夜间灯光亮度增加而不断增加，二者具有一定的相关性，空间重叠度达到64.66%，但相关性不高，主要原因在于POI核密度提取的城市建成区图斑较零散，且在城市边缘处较平滑，而NPP/VIIRS数据对于城市建成区的提取表现为集中连片特征，但分辨率较低，在城市边缘呈锯齿状特征。2)不透水面均值随夜间灯光亮度增加而逐渐增加，相关性也不明显(2=0.4272)，主要是由于不透水面数据分辨率较高，分布也较离散，城市内部灯光亮度值较高区域没有不透水面数据分布所致。
3)重庆市2013-2019年城市扩张趋势明显，扩展速率和扩展强度在年内具有波动趋势，但整体趋势逐年减少，建设重心向东南方向偏移。4)受地形地貌影响，城市以组团式建设为主，空间不规则程度和破碎度呈逐年增加的趋势，城市土地资源利用效率仍需提高，用地增长弹性系数较低，可能存在人口与用地比例失调的情况。当然，本文存在一定的不足之处需要继续深化研究。首先，城市建成区采用简单的阈值法提取，没有考虑多种方法提取效果的对比，以便优化提取方法。其次，对于NPP/VIIRS数据与POI核密度、不透水面的耦合性分析，由于数据来源限制，未采用现势性较高的2019年数据，且没有明确提出二者对NPP/VIIRS数据提取城市建成区边界优化的方法。
参考文献
[1]刘智丽,张启斌,岳德鹏,等.基于Sentinel-2A与NPPVIIRS夜间灯光数据的城市建成区提取[J].国土资源遥感,2019,31(4):227-234
[2]王毓乾,王紫锟,邓志杰,等.兴趣点、夜光遥感和微博签到的空间耦合分析:以北京市为例[J].遥感信息,2019,34(6):18-26
[3]何春阳,史培军,李景刚,等.基于DMSP/OLS夜间灯光数据和统计数据的中国大陆20世纪90年代城市化空间过程重建研究[J].科学通报,2006(7):856-861
[4]陈龙胤.基于DMSP/OLS夜间灯光影像的中国东部沿海地区城市扩张监测与分析[D].长春:东北师范大学,2014
[5]王若曦,李建,李熙,等.DMSP夜间灯光数据与Landsat数据结合的建成区提取研究:以江西省为例[J].华中师范大学学报(自然科学版),2018,52(1):130-136
作者：吴凤敏1，郑稚棚1，梁均军1，陈晓龙1，程宇翔1，李卓锟1
转载请注明来源。原文地址：http://www.xuebaoqk.com/xblw/6560.html
《基于NPP/VIIRS数据的重庆市城市空间发展变化分析》