一种确定夜晚灯光数据提取城市建成区阈值的方法技术

本专利涉及一种确定夜晚灯光数据提取城市建成区阈值的方法，使用数据包括夜晚灯光数据和Googleearth影像数据，处理过程为：a、选取样本城市并按照城市形态分类；b、对夜晚灯光数据进行预处理首先对夜晚灯光数据进行预处理，包括辐射校正及去除异常值等；c、在不同假设阈值情况下对灯光数据进行二值化处理，得到假设阈值的城市建成区与非建成区的边界；d、在Googleearth中根据城市建成区在遥感影像中的表现特征采集样点，将采样点与不同假设阈值提取的城市建成区边界进行叠加，以采样点对表示的建成区和非建成区为标准，确定夜晚灯光数据提取城市建成区的最佳阈值；e、利用得到的最佳阈值结果提取不同城市形态的城市建成区。

【技术实现步骤摘要】
一种确定夜晚灯光数据提取城市建成区阈值的方法
本专利技术属于城市地理中遥感信息技术应用领域，重点解决确定利用夜晚灯光数据提取城市建成区最佳阈值的问题，基于城市形态分类实现地形复杂区域内城市建成区的提取，并以GoogleEarth影像作为参照数据降低研究成本,本专利技术提出的方法具有成本低廉、简单易行的特点，既可应用于城市建成区扩张动态研究，也有助于在空间上分析城市社会经济活动的相关问题。
技术介绍
DMSP/OLS夜晚灯光数据具有与其他遥感影像不同的光电放大特性和夜间灯光的获取能力，因此可以作为监测人类活动的数据源，广泛应用于城市建成区提取、人口密度估测和经济发展评估等多个领域。在目前利用夜晚灯光数据提取城市建成区的研究中，使用最佳阈值分割城市建成区和非建成区的方法简单易行且较为可靠，常用的最佳阈值分割法主要有四种：经验阈值法、突变检测法、统计数据验证法和较高分辨率影像比较法。目前利用夜晚灯光数据提取城市建成区的四种方法各有优缺点：（1）经验阈值法具有较强的主观性，得到的结果因研究者的不同而有差异；（2）统计数据比较法需要以研究区的统计数据为研究基础，受到统计数据准确性的限制，并且只能从建成区面积而不能从空间形态上与夜晚灯光数据提取的建成区进行比较；（3）突变检测法具有较强的客观性，但是数据处理较为繁琐，另外对于受地形因素影响、城市建成区外围没有零星斑块的城市，突变检测法需要依据不同城市形态进行方法上的调整；（4）较高分辨率影像数据空间比较法由于需要以研究区域的TM影像作为辅助数据，所以研究成本相对较高，且数据处理过程相对复杂。同时，基于TM影像光谱提取的城市建成区中可能混杂着光谱相似的非建成区地物，比如裸露的山地、沙漠等；城市中成片的公园广场用地和水体可能被划分为非城市建成区，但是从生态和社会服务属性来说，城市公园广场属于城市建成区的范畴。而在利用夜晚灯光数据提取建成区时，由于公园本身具有的灯光和附近辐射灯光的影响，公园绿地和水体被划分成城市建成区，这是其他较高分辨率影像提取城市建成区不精确的因素之一。
技术实现思路
本专利技术的目的是根据两种形态对城市进行分类，确定利用夜晚灯光数据提取不同形态城市建成区的最佳阈值，并且将阈值结果应用于城市建成区的提取。以GoogleEarth影像为标准数据使提取结果更加准确，阈值结果可以用于地形复杂或者具有不同城市形态的区域中提取城市建成区。具体包括以下内容：一种确定夜晚灯光数据提取城市建成区阈值的方法技术流程图见图1，主要包括以下步骤：A、夜晚灯光数据的处理方法a、辐射定标：根据统计数据选择与全国经济发展平均值接近并且发展较为稳定的城市作为辐射定标区域，选取该城市夜晚灯光数据累计像元值最大的年份作为标准年份进行辐射定标。b、辐射校正：利用（其中a、b、c是定标区域各年份之间的校正系数，DN是原始数据，是校正以后的结果）分析各年度灯光数据像元值与定标年份灯光像元值之间的关系，确定系数a、b、c以后，对样本城市灯光数据进行辐射校正。c、去异常值：传感器在接收和传输数据过程中存在不确定性，夜晚灯光数据的像元存在异常值，因此假设随着时间的变化，城市建成区的范围越来越大，灯光亮度值也越来越大，依据此规则对夜晚灯光数据去异常值。d、重采样：为增强两种数据的直观比较以及使灯光影像中城市区域的轮廓更贴近现实情况，对夜晚灯光数据进行三次卷积重采样，使空间分辨率调整到50m。e、二值化：将夜晚灯光数据按照不同的假设阈值进行二值化，将城市建成区赋值为1，非建成区赋值为-1，进而得到城市建成区的轮廓。B、GoogleEarth数据处理方法：利用GoogleEarth影像数据，已有城市的行政边界作为掩膜，依据城市建成区在遥感影像中的显示特征通过目视解译判别城市建成区的边界，在城市建成区和非城市建成区各采一个点形成正负值点对，正值点表示城市建成区，负值点表示非城市建成区，并且两点形成的直线垂直于建成区的边界。由于原始夜晚灯光数据的分辨率为1km，为保证两点落在相邻的栅格内，两点的距离在1.5km到2km之间。C、将GoogleEarth中提取的点对与夜晚灯光数据中的二值化结果叠加，记录符合要求（即正值点落在城市建成区1范围内，负值点落在非城市建成区范围内）的点对数量，统计准确率p=r/t，其中r为符合要求点对数量，t为点对总数。D、按照城市形态分类结果，统计两种城市形态使用夜晚灯光数据提取城市建成区的最佳阈值并且进行应用。附图说明图1是一种确定夜晚灯光数据提取城市建成区阈值的方法技术流程图；图2是宁波市夜晚灯光二值化结果与GoogleEarth点对数据叠加示意图；图3是贵阳市夜晚灯光数据提取城市建成区结果图；图4是兰州市夜晚灯光数据提取城市建成区结果图；图5是拉萨市夜晚灯光数据提取城市建成区结果图；图6是三明市夜晚灯光数据提取城市建成区结果图；图7是西宁市夜晚灯光数据提取城市建成区结果图；图8是长春市夜晚灯光数据提取城市建成区结果图；图9是长沙市夜晚灯光数据提取城市建成区结果图；图10是宁波市夜晚灯光数据提取城市建成区结果图；图11是沈阳市夜晚灯光数据提取城市建成区结果图；图12是厦门市夜晚灯光数据提取城市建成区结果图。具体实施方式本专利技术一种确定夜晚灯光数据提取城市建成区阈值的方法可以通过以下实施例进行说明：A、选取两种城市形态共十个城市作为样本提取城市建成区，将其分为两种城市形态：带状城市和非带状城市，带状城市为三明、贵阳、兰州、拉萨、西宁，非带状城市为宁波、厦门、沈阳、长春、长沙。B、通过《中国城市统计年鉴》取得2007-2012年各城市地区生产总值增长率以及人口密度增长率数据，根据统计数据分析结果确定以潍坊市作为定标区域，分析2007-2012年的灯光像元值确定以2010年为定标年份，其他年份与2010年像元值之间的相关关系如表1，表1：定标区2010年的灯光像元值与其他年份之间的回归分析及系数根据表1对其他城市的夜晚灯光数据进行辐射校正，得到辐射校正后各城市夜晚灯光数据。C、获取2007年至2012年各样本城市夜晚灯光数据，根据公式（1）对各城市夜晚灯光数据去异常值，其中DN（n-1,i），DN（n,i），DN（n+1,i）分别为第n-1，n，n+1年的灯光数据像元值，本文中n=2008,2009,2010,2011D、重采样：对去异常值后的各城市夜晚灯光影像进行三次卷积重采样处理，使空间分辨率调整到50m。E、二值化：将各城市夜晚灯光数据按照假设的不同阈值进行二值化，将二值化结果分为城市建成区1和非建成区-1，进而得到城市建成区的轮廓。F、根据城市建成区在GoogleEarth影像中的显示特征，利用目视解译的方法确定城市建成区与非建成区，并且分别在建成区和非建成区上面采样点形成正负点对，为保证实验的合理性和准确性，每个城市选点100个即50对点并且均匀分布。G、在ArcGIS9.3软件支持下，将GoogleEarth中提取的点对与夜晚灯光数据中提取的城市建成区轮廓叠加，统计叠加数据结果，表2为每个城市对应某一阈值时的准确率p以及每个阈值下所有城市准确率p的方差，选取准确率高并且方差较小的像元值作为最佳阈值。表2为每个城市在40-50假设阈值的情况下对应的准确率以及最佳阈值结果根据不本文档来自技高网...

【技术保护点】
一种确定夜晚灯光数据提取城市建成区阈值的方法，其特征在于以下步骤：A、将样本城市按照城市形态分为两类：带状城市和非带状城市；B、将夜晚灯光数据进行预处理，在Google Earth平台中进行城市建成区和非建成区点对采样，将两种数据的处理结果导入Arcgis9.3软件中叠加显示，并且统计。

【技术特征摘要】
1.一种确定夜晚灯光数据提取城市建成区阈值的方法，其特征在于以下步骤：A、将样本城市按照城市形态分为两类：带状城市和非带状城市；B、将夜晚灯光数据进行预处理，在GoogleEarth平台中进行城市建成区和非建成区点对采样，将两种数据的处理结果导入Arcgis9.3软件中叠加显示，并且统计，具体包括：将经过预处理后的夜晚灯光数据根据一系列不同假设的阈值提取出不同阈值情况下的城市建成区边界；在GoogleEarth平台中，根据城市建成区与非建成区不同的表现特征进行点对采样，正值点表示城...

【专利技术属性】
技术研发人员：李新虎，宋金超，张国钦，
申请(专利权)人：中国科学院城市环境研究所，
类型：发明
国别省市：福建;35