基于landsa-8卫星遥感数据的第三次全国国土调查统一时点更新研究.docxVIP

基于landsa-8卫星遥感数据的第三次全国国土调查统一时点更新研究 摘要：以榆中县为样本，研究利用Landsat-8卫星遥感数据，提取和分析不透水面变化信息，用于辅助第三次全国国土调查统一时点更新工作。研究结果表明，使用归一化差值不透水面算法提取的不透水面变化信息准确度高、适用范围大，可以为统一时点更新工作提供指引，减轻人工目视判读寻找变化区的工作量并避免遗漏。 0 引言近日，国务院第三次全国国土调查领导小组办公室印发《关于开展第三次全国国土调查统一时点更新调查的通知》，对第三次全国国土调查(以下简称“三调”)统一时点更新工作进行部署，要求各地严格按照三调《总体方案》和《实施方案》等文件要求，以2019年12月31日为标准时点，全面掌握三调初始时点和统一时点间的地类、面积和权属变化情况，在已完成的三调初始成果基础上更新“三调”数据库，保证调查成果的统一性。可以看出，统一时点更新工作的重点就是发现调查时点与标准时点之间的地类变化图斑，而单纯靠人工目视解译比对，要查清发现几千甚至上万平方千米县域范围内的土地利用现状变化情况，工作量巨大，且难以避免疏漏。目前，基于多时相遥感数据进行变化检测，可以较为可靠地提取土地利用类型的变化信息，已经广泛应用到城市发展变迁、生态环境评价等领域利用遥感数据提取不透水面信息一般使用建模法1 研究区及遥感数据1.1 研究区本文以甘肃省兰州市榆中县作为研究区。榆中县位于甘肃省中部，全县地势南高北低，中部凹，呈马鞍形。南部为石质高寒山区，马衔山最高峰海拔3 670.3 m。北部为黄土丘陵区，最高峰吕家岘海拔2 495 m。南北两山之间是川塬丘陵沟壑区，海拔1 500—2 000 m，土地面积约为3 294 km1.2 遥感数据美国从1972年起至今陆续发射了8颗“陆地卫星”(Landsat)，其最新型号为Landsat-8，于2013年2月11日成功发射。Landsat系列卫星以近极地圆形太阳同步轨道运行，能在同一地方时、同一运行方向通过同一观测点，这就保证了同一研究区每幅影像的观测条件基本一致，非常有利于对不同时期图像进行对比分析。Landsat-8卫星搭载有OLI陆地成像仪和TIRS热红外传感器，同以前的Landsat系列卫星相比新增了一些波段:新增了海蓝波段，主要应用在对海岸带观测;新增了短波红外波段，又称为卷云波段，用于云去除及其他应用。同时，为了避免大气吸收部分特征，Landsat-8对原有波段也进行了部分调整，如缩减了近红外波段(第5波段)的波长范围，排除了0.825μm处水汽吸收特征。除此之外，与Landsat-7相比，Landsat-8最重要的变化是将热红外波段拆分为两个，空间分辨率降低到100 m，但是却对光谱信息进行了细化。Landsat-8数据波段参数见表1。本文研究中选取的遥感数据是2017年12月24日和2019年12月17日覆盖研究区的两幅Landsat-8影像，主要使用其OLI陆地成像仪数据中可见光波段的绿光、近红外波段和中红外波段以及TIRS热红外传感器中的第10波段数据。2 研究方法及结果分析2.1 归一化差值不透水面指数徐涵秋提出的归一化差值不透水面指数(NDISI)模型，可以实现对不透水面信息的快速、准确且自动化提取。模型构建思路是分析多光谱遥感数据中对不透水面辐射最强和最弱的波段组合，并对其进行运算，以进一步增大二者辐射强度差距，使不透水面信息得到显著增强，然后以比值的方式计算得到不透水面信息。该指数模型可以用如下的复合波段运算方程表示:上式中，TIR,NIR和MIR式中，MNDWI即为改进型归一化水体指数:式中，Green为可见光的绿光波段，对应Landsat-8卫星的陆地成像仪第3波段。NDISI为归一化差值指数模型，运算结果介于-1—+1之间，若不经人工判读，直接选取0为阈值，则小于0的为被增强的不透水信息，受抑制的其他地物信息大于或等于02.2 结果分析本文选取的两期Landsat-8卫星数据，获取日期均在12月20日前后，且拍摄时间基本一致，均在下午三点四十分左右。因此，两期影像拥有较为一致的光谱特性，太阳高度角也基本一致。用两期影像的反演结果求差，不仅可以剔除未变化区域和山体阴影等噪声影响，而且可以增大变化区域的影像亮度值。截取榆中县城建成区中心的两期不透水面提取结果如图1、图2所示，图中亮度值较高区域为反演得到的不透水面。可以看出，在2017年12月至2019年12月之间，榆中县城建设规模不断扩大，县城北侧及东侧建筑区都发生了较大变化，将两期反演结果求差，得到不透水面变化区域，如图3所示。为便于使用，选取一定阈值将得到的不透水面变化结果进行二值化处理，最终得到新增