一种基于夜光遥感数据的城市发展重心轨迹的获取方法技术

本发明专利技术公开了一种基于夜光遥感数据的城市发展重心轨迹的获取方法，包括如下步骤：(1)筛选夜光遥感数据；(2)对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据，得到各年份的DMSP/OLS或NPP/VIIRS影像；(3)对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；(4)对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹。本发明专利技术运用DMSP/OLS及其NPP/VIIRS夜间灯光遥感数据来计算城市发展重心，为理解区域发展规划和经济可持续发展具有重要意义。

A method of obtaining the center of gravity of urban development based on night light remote sensing data

【技术实现步骤摘要】
一种基于夜光遥感数据的城市发展重心轨迹的获取方法
本专利技术涉及遥感
，尤其是一种基于夜光遥感数据的城市发展重心轨迹的获取方法。
技术介绍
随着遥感技术的进步，遥感影像分辨率的不断提高，遥感在各行各业都将有着广阔的应用前景。由美国国防气象卫星搭载的可见光成像线性扫描业务系统(DMSP/OLS)和国家极轨卫星搭载的可见光近红外成像辐射仪(NPP/VIIRS)获取的夜间灯光影像是监测人类社会经济活动的主要数据源。目前所能获得的夜光影像有DMSP/0LS的年产品和NPP/VIIRS所得到的月产品和年产品，并且发布的这两种夜光遥感影像全都为全球合成产品，不需要进行大气校正、镶嵌等处理。城市化是发展中国家经济发展的重要推动力之一，同时城市化也深刻影响了全球和区域气候变化。因此，对城市化进行监测和评估有利于理解全球社会经济发展和气候变化。由于城市在夜间发出灯光，因此夜光影像可以用来准确提取和分析建成区范围和空间聚集现象(城市群)，随着城市化的加速，城市建成区逐渐扩大，研究城市的发展重心，准确获取区域的发展信息对于指导区域发展规划和经济可持续发展具有重要意义，对于理解全球和区域城市化进程可以发挥重要作用，也能为进一步科学地保护土地资源，改善生态环境的重要决策奠定基础。
技术实现思路
本专利技术所要解决的技术问题在于，提供一种基于夜光遥感数据的城市发展重心轨迹的获取方法，运用DMSP/OLS及其NPP/VIIRS夜间灯光遥感数据来计算城市发展重心，为理解区域发展规划和经济可持续发展具有重要意义。r>为解决上述技术问题，本专利技术提供一种基于夜光遥感数据的城市发展重心轨迹的获取方法，包括如下步骤：(1)筛选夜光遥感数据；(2)对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据，得到各年份的DMSP/OLS或NPP/VIIRS影像；(3)对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；(4)对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹。优选的，步骤(2)中，预处理包括夜光遥感影像背景噪声去除及图像之间相互校正；首先建立2012年DMSP/OLS与NPP/VIIRS两种数据源的对数回归模型，再利用所建立的模型对1992年至2011年的DMSP/OLS数据进行校正，这样就得到了1992年至2018年长时间序列一致的夜间灯光影集，并移除DN值0以下的点进行背景噪声去除；对数回归模型的式子如下：Y＝a×log10(b×X)其中Y为DMSP/OLS数据经过转换至NPP/VIIRS尺度上的像元DN值，X为DMSP/OLS数据像元值，a、b为对数回归得到的校正参数。优选的，步骤(4)中，对所得数据进行处理计算中，计算公式如下：其中X、Y分别为t年中n区域的夜间灯光重心的经度坐标和纬度坐标；li为第i个像素的灯光值；Mn为区域n中的像元总数；xi、yi分别为第i个像素点的经度坐标和纬度坐标。本专利技术的有益效果为：本专利技术运用DMSP/OLS及其NPP/VIIRS夜间灯光遥感数据来计算城市发展重心，该夜光遥感数据具有获取上的便捷性、时间上的动态性和地理空间上的客观性的特性；结合了这些特性遥感影像得出的重心轨迹，可为后期研究全国甚至世界范围的城市化扩张及发展轨迹提供了可靠有效的支撑数据，为理解区域发展规划和经济可持续发展具有重要意义。附图说明图1为本专利技术的方法流程示意图。图2(a)为本专利技术南京市地区矢量图层示意图。图2(b)为本专利技术南京市2018年夜间灯光影像示意图。图3为本专利技术南京市1992年到2018年发展重心移动轨迹示意图。具体实施方式如图1所示，一种基于夜光遥感数据的城市发展重心轨迹的获取方法，包括如下步骤：(1)筛选夜光遥感数据；(2)对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据，得到各年份的DMSP/OLS或NPP/VIIRS影像；(3)对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；(4)对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹。以下选取南京城市圈为例将结合附图和实施例的流程详细说明本专利技术技术方案。实施例提供的基于夜光遥感数据的城市发展重心轨迹具体实现步骤如下：步骤一，筛选夜光遥感数据；因为NPP/VIIRS对夜间灯光有更高的灵敏度，能够更准确反映人类经济社会活动。所以选择2012年到2018年的NPP/VIIRS影像及DMSP/OLS的从1992年到2012年的影像。在NOAA(NationalOceanicandAtmosphericAdministration，美国国家海洋和大气管理局)网站下载DMSP/OLS及NPP/VIIRS影像。步骤二，对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据；由于夜光影像中并未去除火光、气体燃烧、火山和极光等，相应的背景噪声也未过滤，并且这两种数据在长时间序列上不具有一致性和可比性，难以直接应用于多源夜间灯光影像集的相关研究。具体实施例中需要对步骤二中的遥感影像进行相互校正，去除背景噪声；首先建立2012年DMSP/OLS与NPP/VIIRS两种数据源的对数回归模型，再利用所建立的模型对1992年至2011年的DMSP/OLS数据进行校正，这样就得到了1992年至2018年长时间序列一致的夜间灯光影集，并移除DN值0以下的点进行背景噪声去除。对数回归模型的式子如下：Y＝a×log10(b×X)其中Y为DMSP/OLS数据经过转换至NPP/VIIRS尺度上的像元DN值，X为DMSP/OLS数据像元值，a、b为对数回归得到的校正参数。在NaturalEarthData网站下载南京市区域边界矢量图如图2(a)；使用ArcGIS10.3软件中的SpatialAnalyst工具-提取分析-按掩膜提取，提取出南京市从1992年到2018年的夜光遥感影像，如图2(b)。步骤三，对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；针对步骤三，实施例利用MATLAB程序进行统计分析出每张图即南京市从1992到2018年每年的夜光遥感图中的灯光总值和各像素点的灯光值及经纬度坐标。步骤四，对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹；针对步骤四，实施例利用计算公式处理数据，计算公式如下：其中X、Y分别为t年中n区域的NTL重心的经度坐标和纬度坐标；li为第i个像素的灯光值；Mn为区本文档来自技高网...

【技术保护点】
1.一种基于夜光遥感数据的城市发展重心轨迹的获取方法，其特征在于，包括如下步骤：/n(1)筛选夜光遥感数据；/n(2)对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据，得到各年份的DMSP/OLS或NPP/VIIRS影像；/n(3)对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；/n(4)对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹。/n

【技术特征摘要】
1.一种基于夜光遥感数据的城市发展重心轨迹的获取方法，其特征在于，包括如下步骤：
(1)筛选夜光遥感数据；
(2)对筛选出的夜光遥感数据进行预处理，并裁剪出目标城市相应历史年份的区域夜光影像数据，得到各年份的DMSP/OLS或NPP/VIIRS影像；
(3)对所得区域夜光遥感数据进行分析，统计得到各历史年份的灯光总值和其各像素点的灯光值及经纬度坐标；
(4)对所得数据进行处理计算，就能得到其夜间灯光重心即区域发展的重心，统计出各历史年份的重心，即可生成城市发展重心轨迹。


2.如权利要求1所述的基于夜光遥感数据的城市发展重心轨迹的获取方法，其特征在于，步骤(2)中，预处理包括夜光遥感影像背景噪声去除及图像之间相互校正，建立DMSP/OLS与NPP/VIIRS两种...

【专利技术属性】
技术研发人员：陈嘉琪，王巍，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32