一种居民地提取方法技术

本发明专利技术公开了一种居民地提取方法，它包括采集研究区的Landsat 8 OLI传感器遥感影像和行政边界数据；对遥感影像进行预处理；计算预处理后遥感影像的地表反照率；根据地表反照率提取水体和阴影；利用单波段阈值法在水体和阴影图层中提取出阴影，从而区分出水体和阴影；对水体提取结果进行闭运算滤波处理；计算预处理后遥感影像的归一化植被指数确定阈值来提取植被；在原始研究区边界图层中裁剪掉上述步骤提取出的水体、阴影和植被图，剩下的图层即为研究区居民地和道路分布图；计算居民地和道路分布图的形状指数并确定形状指数(I)的阈值提取出道路和居民地；解决了现有居民地提取存在的工作强度大、精度低和效率低等问题。

【技术实现步骤摘要】
一种居民地提取方法
本专利技术属于居民地提取技术，尤其涉及一种居民地提取方法。
技术介绍
居民地又称“居民点”或聚落,是人与自然相互作用的产物,是人类从事生产和生活的需要而集聚定居的各种形式的居住场所,是自然景观、建筑景观和文化景观的综合体,是人们生活和生产的基本场所。准确提取居民地，是居民地规划、发展和变化监测的最基本的地理数据来源，对居民地变化以及影响具有深刻的意义，有利于实现区域的可持续发展，并且对研究地区经济社会发展、城镇规划、土地资源保护等相关研究存在重要作用。因此，建立快速、准确的居民地提取方法是非常亟需的。传统的居民地提取方法主要是通过ENVI和ARCGIS等软件，采用监督分类，非监督分类，目视解译，人机交互解译，野外实地调查等方法对遥感影像进行居民地提取。这些方法不仅受不同解译人员经验和主观判断影响大，还受“同物异谱、同谱异物”等问题的影响，而且工作强度大、精度低、效率低。这些方法对居民地信息表征能力存在很大偏差，具有一定的局限性，没有系统地排除干扰居民地提取的其他地物，准确性有待提高。
技术实现思路
本专利技术要解决的技术问题是：提供一种居民地提取方法，以解决现有技术的居民地提取存在的工作强度大、精度低和效率低等技术问题。本专利技术的技术方案是：一种居民地提取方法，它包括：步骤1、数据采集：采集研究区的Landsat8OLI传感器遥感影像和行政边界数据；步骤2、遥感影像数据预处理：对遥感影像进行辐射定标、大气校正、图像融合和图像裁剪；步骤3、计算步骤2预处理后遥感影像的地表反照率(Albedo)；步骤4、根据地表反照率提取水体和阴影；步骤5、利用单波段阈值法在水体和阴影图层中提取出阴影，从而区分出水体和阴影；步骤6、对水体提取结果进行闭运算滤波处理；步骤7、计算步骤2预处理后遥感影像的归一化植被指数(NDVI)；步骤8、为归一化植被指数(NDVI)确定阈值来提取植被；步骤9、在原始研究区边界图层中裁剪掉上述步骤提取出的水体、阴影和植被图，剩下的图层即为研究区居民地和道路分布图；步骤10、计算居民地和道路分布图的形状指数；步骤11、根据形状指数(I)的阈值提取出道路；步骤12、提取居民地。利用行政边界数据裁剪遥感影像得到研究区遥感影像图。步骤3所述地表反照率(Albedo)的计算公式为：α＝0.356*B2+0.130*B4+0.373*B5+0.085*B6+0.072*B7-0.018式中：B2为Landsat8OLI蓝光波段；B4为Landsat8OLI红外波段；B5为Landsat8OLI近红外波段；B6为Landsat8OLI中红外波段；B7为Landsat8OLI中红外波段。步骤4所述根据地表反照率提取水体和阴影的方法为：在ENVI软件中统计地表反照率的直方图，根据直方图上的数值确定阈值的取值范围，然后对每类地物选择训练样本，再统计每种地物的地表反照率并按波段进行排序，根据排序确定出提取地物的阈值；通过ENVI软件里面的波段计算器工具将小于该阈值的区域提取出来，这些区域即水体和阴影。步骤6所述对水体提取结果进行闭运算滤波处理的方法为：利用ENVI软件中的闭运算滤波工具对水体提取结果进行处理，填补水体提取结果中的孔洞和不连续区域。步骤7所述计算步骤2预处理后遥感影像的归一化植被指数(NDVI)的计算公式为：NDVI＝(NIR-Red)/(NIR+Red)式中：NIR为近红外波段的反射值，Red为Landsat8可见光红光波段值。步骤8所述为归一化植被指数(NDVI)确定阈值来提取植被的方法为：利用步骤7中计算出的NDVI分布图中植被这类地物的NDVI分布特征来的分布范围，再在分布范围内选取3个值进行植被的提取，在ENVI软件的波段计算器中将NDVI大于阈值的范围提取为植被，得到3个植被提取结果；利用目视解译将3个提取结果与原始遥感影像对比，与原始遥感影像中植被分布区重合度最好的提取结果作为最终的植被提取结果。步骤10所述的计算居民地和道路分布图的形状指数的公式为：S＝A/P式中S为居民地和道路分布图中图斑的面积，P为居民地和道路分布图中图斑的周长。步骤11所述确定形状指数(I)的阈值提取出道路的方法为：计算面积和周长的比值得到形状指数，并加以统计，将形状指数小于0.1的区域确定为道路。步骤12所述提取居民地的方法为：在ENVI软件中将居民地图层中裁剪掉道路得到居民地最终提取结果图。本专利技术有益效果：本专利技术首次结合单波段阈值法、归一化植被指数(NDVI)、地表反照率(Albedo)、闭运算滤波法以及形状指数这些参数和方法来提取居民地；在确定提取各个地物的阈值时，有针对性得利用不同的方法来确定每一个参数或方法的阈值；经过验证，提取的精度达到91.25％，实现了快速、准确的居民地的提取；解决现有技术的居民地提取存在的工作强度大、精度低和效率低等技术问题。具体实施方式本专利技术提供了一种居民地提取方法，包括如下步骤：第一，数据获取。居民地提取所需要的数据为Landsat8OLI传感器遥感影像和研究区行政边界数据。这些数据可从地理空间数据云网站下载(网址：http://www.gscloud.cn/)。第二，遥感影像数据预处理。数据的预处理主要包括对遥感影像进行辐射定标、大气校正、图像融合和图像裁剪(图像裁剪指利用研究区边界裁剪遥感影像得到研究区遥感影像图)，以便于增加影像的分辨率以及为后续分析和研究奠定基础，该步骤在ENVI软件中实现。第三，计算预处理后遥感影像的地表反照率(Albedo)。使用Liang建立的用Landsat数据估算地表反照率的通用公式来计算，公式为：α＝0.356*B2+0.130*B4+0.373*B5+0.085*B6+0.072*B7-0.018其中，B2为Landsat8OLI蓝光波段；B4为Landsat8OLI红外波段；B5为Landsat8OLI近红外波段；B6为Landsat8OLI中红外波段；B7为Landsat8OLI中红外波段。第四，为步骤三计算得到的地表反照率确定阈值来提取水体。在ENVI软件中统计地表反照率的直方图，根据直方图上的数值确定阈值的取值范围，然后对每类地物选择训练样本，再统计每种地物的地表反照率并按波段进行排序，根据排序确定出提取地物的具体阈值。通过ENVI软件里面的波段计算器工具将小于该阈值的区域提取出来，这些区域即水体和阴影。第五，利用单波段阈值法在水体和阴影图层中提取出阴影，从而区分出水体和阴影。单波段阈值法是基于TM影像提出的一种地物提取算法，主要利用了某种地物在某个波段处的强吸收性和在某波段范围内强反射性的特点，找出某一地物能明显区别于其他地物的波段。首先在ENVI软件中对每种地物选择训练样本，再统计每种地物的光谱值，并按波段和按地物进行两种排序，根据两种排序，找出某一地物能明显区别于其他地物的值(比如光谱值远小于或远大于其他地物)。利用这一阈值筛选出这一地物。利用该方法区分出阴影和水体后，得到分别得水体和阴影分布图。第六，对水体提取结果进行闭运算滤波处理。在第四步提取出水体之后，利用ENVI软件中的闭运算滤波工具对水体提取结果进行处理，填补水体提取结果中的孔洞、不连续等，从而提高水体提取的精度。第七，计算预处理后遥感影像的本文档来自技高网...

【技术保护点】
1.一种居民地提取方法，它包括：步骤1、数据采集：采集研究区的Landsat 8OLI传感器遥感影像和行政边界数据；步骤2、遥感影像数据预处理：对遥感影像进行辐射定标、大气校正、图像融合和图像裁剪；步骤3、计算步骤2预处理后遥感影像的地表反照率(Albedo)；步骤4、根据地表反照率提取水体和阴影；步骤5、利用单波段阈值法在水体和阴影图层中提取出阴影，从而区分出水体和阴影；步骤6、对水体提取结果进行闭运算滤波处理；步骤7、计算步骤2预处理后遥感影像的归一化植被指数(NDVI)；步骤8、为归一化植被指数(NDVI)确定阈值来提取植被；步骤9、在原始研究区边界图层中裁剪掉上述步骤提取出的水体、阴影和植被图，剩下的图层即为研究区居民地和道路分布图；步骤10、计算居民地和道路分布图的形状指数；步骤11、根据形状指数(I)的阈值提取出道路；步骤12、提取居民地。

【技术特征摘要】
1.一种居民地提取方法，它包括：步骤1、数据采集：采集研究区的Landsat8OLI传感器遥感影像和行政边界数据；步骤2、遥感影像数据预处理：对遥感影像进行辐射定标、大气校正、图像融合和图像裁剪；步骤3、计算步骤2预处理后遥感影像的地表反照率(Albedo)；步骤4、根据地表反照率提取水体和阴影；步骤5、利用单波段阈值法在水体和阴影图层中提取出阴影，从而区分出水体和阴影；步骤6、对水体提取结果进行闭运算滤波处理；步骤7、计算步骤2预处理后遥感影像的归一化植被指数(NDVI)；步骤8、为归一化植被指数(NDVI)确定阈值来提取植被；步骤9、在原始研究区边界图层中裁剪掉上述步骤提取出的水体、阴影和植被图，剩下的图层即为研究区居民地和道路分布图；步骤10、计算居民地和道路分布图的形状指数；步骤11、根据形状指数(I)的阈值提取出道路；步骤12、提取居民地。2.根据权利要求1所述的一种居民地提取方法，其特征在于：步骤2所述图像裁剪的方法为：利用行政边界数据裁剪遥感影像得到研究区遥感影像图。3.根据权利要求1所述的一种居民地提取方法，其特征在于：步骤3所述地表反照率(Albedo)的计算公式为：α＝0.356*B2+0.130*B4+0.373*B5+0.085*B6+0.072*B7-0.018式中：B2为Landsat8OLI蓝光波段；B4为Landsat8OLI红外波段；B5为Landsat8OLI近红外波段；B6为Landsat8OLI中红外波段；B7为Landsat8OLI中红外波段。4.根据权利要求1所述的一种居民地提取方法，其特征在于：步骤4所述根据地表反照率提取水体和阴影的方法为：在ENVI软件中统计地表反照率的直方图，根据直方图上的数值确定阈值的取值范围，然后对每类地物选择训练样本，再统计每种地物的地表反照率...

【专利技术属性】
技术研发人员：白晓永，杨钰杰，李朝君，肖建勇，钱庆欢，
申请(专利权)人：中国科学院地球化学研究所，
类型：发明
国别省市：贵州,52