城市植被覆盖度遥感信息提取

城市植被覆盖度遥感信息提取摘要：本文基于SPOT5影像，以连云港市为实验区，研究了城市植被信息遥感提取的方法和技术。通过对研究区SPOT5影像的近红外波段、红光波段和绿光波段典型地物光谱信息的统计分析和对比研究，发现NDVI植被指数法用于增强绿色植被的效果最好。研究成果对于连云港地区绿化结构，优化植被空间结构，使城市植被充分、高效地发挥其生态效益和使用功能具有重要的理论意义和实践价值。关键词：遥感；城市植被；SPOT5；植被指数；信息提取1概述城市植被作为城市生态系统的重要组成部分，对于改善城市环境质量，提高居民生活水平具有重要作用。因此，城市植被的研究是人们对城市发展预测至关重要的任务。随着遥感技术的不断发展，遥感技术已经广泛运用于各个领域。在城市规划方面，遥感技术可以实现土地动态监测，空气质量的监督控制和城市环境的规划建设等。近几年国内许多城市将遥感技术应用于植被覆盖度信息提取中，动态掌握植被覆盖的区域，优化植被空间结构，提高城市可持续发展潜能，实现植被的整体规划。利用遥感影像进行城市植被覆盖度信息提取，比传统方式更加优越。它具有视域范围广、信息量多、重复周期短、图像清晰、资料收集方便等优点。城市植被遥感主要研究城市绿化系统分析及规划，是遥感技术应用的重要领域之一，利用遥感技术不仅能够准确判定和量测绿化覆盖面积，且对于判别植被的类型、结构乃至识别植物种类等都显得十分有效。随着遥感影像分辨率的不断提高，国内大部分城市己经采用如TM、SPOT、ETM+等影像数据来进行城市植被的调查及生态规划，从而为城市生态规划及城市建设提供方案及依据。因此，本文基于SPOT5遥感影像对连云港某地区的城市植被进行提取。首先，对原始遥感影像进行裁剪得到连云港地区的影像图。其次，运用多种植被指数法对研究区图像进行增强处理，并通过比较得出NDVI（归一化植被指数）对植被的增强效果最好。最后，在NDVI指数图像上通过统计阈值，提取研究区的植被信息。2城市植被在SPOT5影像中的特征2.1城市植被及其地理特征2.1.1城市植被城市植被指城市范围内的全部植被，包括自然生长的和人工栽培的各种植被类型。城市植被的定义分为广义城市植被和狭义城市植被。广义城市植被指城市规划区范围内的各种植被。包括六大类型：公共植被，即各种公园、休憩林荫带；居住区植被；交通植被；附属植被；生产防护植被；位于市内或城郊的风景区植被，即风景游览区、休养区、疗养区等。狭义城市植被指面积较小、设施较少或没有设施的绿化地段，区别于面积较大、设施较为完善的“公园”，主要包括公共植被、生产和防护植被两类。本文中所说的城市植被指广义城市植被。城市植被不仅为城市环境增添了色彩，改善了城市气候，也是一个区域发展的重要影响因素。2.1.2研究区植被的地理特征连云港全区位于东经118°24′～119°48′和北纬34°～35°07′之间，东西长129公里，南北宽约132公里，水域面积1759.4平方公里。东临黄海，与朝鲜、韩国、日本隔海相望；西与徐州新沂市、宿迁市沭阳县毗邻；南与淮安市涟水、盐城市响水2县相连；北与山东郯城、临沭、莒南、日照等县市接壤。连云港市位于鲁中南丘陵与淮北平原的结合处，地势由西北向东南倾斜，形如一只飞向海洋的彩蝶。境内以平原为主，兼有丘陵、山地、湖泊、滩涂等。有大小山峰214座，云台山主峰玉女峰海拔624.4米，为江苏省的最高峰。境内河网稠密，连云港市有标准海岸线162公里，21个岛屿，其中东西连岛为江苏第一大岛，面积7.57平方公里，基岩海岸是江苏省独有。2.2城市植被的遥感影像特征遥感影像特征主要有光谱特征和空间特征，我们了解遥感影像特征的目的是为了更好的解译遥感影像中的各类地物，为军事、农业、环境、地学和各类决策提供有效的技术支撑。因此研究遥感影像的特征就是一项非常重要而且困难的工作[1]。2.2.1光谱特征图像上的波谱信息表现为已经量化的辐射值，即图像的亮度/灰度值/像元值，它是一种相对的量度。像元值间接反映了地物的波谱特征，不同的地物有着不同的像元值，当然把同谱异物排除。遥感图像解译中识别不同地物的一个重要标志就是图像的像元值差异。光谱特征可具体为提取目标物的颜色、灰度或波谱间的亮度比等。SPOT5影像是由SPOT5卫星搭载3种传感器――除了前几颗卫星上的高分辨率几何装置（HRVIR）和植被探测器（VEGETATION）外，还有一个高分辨率立体成像（HRS）装置获取的遥感影像，包括以下5个波段[2]。PA：0.51～0.73um，全色波段。地面分辨率较高，为2.5m，可用于农林调查和规划，城市规划和较大比例尺专题制图。B1：0.50～0.59um，为绿波段。波段中心位于叶绿素反射曲线最大值，即0.55um处，处于水蒸气衰减最小值的长波端，对于水体混浊度评价以及水深10～20m以内的干净水体的调查是十分有用的。B2：0.61～0.68um，为红波段。位于叶绿素吸收带，受大气散射的影响较小，为可见光最佳波段，用于识别裸露的地表、植被、土壤、岩性地层、地貌现象等。B3：0.78～0.89um，为近红外波段。能够很好地穿透大气层。在该波段，植被表现的特别明亮，水体表现的特别黑。B4：1.58～1.75um，为短红外波段。用于探测植物含水量及土壤湿度，区别云与雪。由于城市植被信息在B1、B2和B3波段的光谱特征明显，所以本文选择这三个波段的数据进行城市植被信息的提取研究。在不同植物的光谱曲线中，反射率涨落最大的是近红外区。不同的植物种类，如阔叶乔、灌木和草地在该波段具有较高的反射率，而针叶树的反射率也会有10%～20%的涨落，一般春末夏初反射较高，而秋季较低。此外植物的健康程度也会在近红外波段有所显示，对于同一种植物而言，受环境污染及病虫害影响的反射率较低，健康植物的反射率较高。植被指数有助于增强遥感影像的解译能力，并已作为一种遥感手段广泛应用于土地利用覆盖探测、植被覆盖密度评价、作物识别和作物预报等方面，并在专题制图方面增强了分类能力。费鲜芸[3]在研究山东省泰安市城区60km2范围内的城市植被时，运用SPOT52.5m全色波段和10米多光谱波段，对该城区的各种城市植被、耕地、居民区和水域等地物进行光谱特征分析，对各波段的灰度值范围及灰度平均值进行统计，统计结果如图1所示。图1几种典型地物的光谱曲线各种植被、耕地和水域在红色（red）波段（第二波段），近红外（nearinfrared）波段（第三波段）与道路和居民区能较好区分，同时裸地与不同植被类型在各波段有不同程度的光谱重叠，利用单一波段难以区分。水域在第一波段和第四波段反射率都较低，可以与其他各类别较好地区分，个别水域在水泥筑底的浅水区域存在少量的与道路有较少的重叠现象，会影响分类精度。本文通过对连云港某地区遥感图像的亮度值进行统计，得到几类典型地物的光谱信息，如图2所示。图2典型地物的光谱曲线从图中可以看出四类地物的光谱曲线变化相似，纵向比较时发现，植被的亮度值与其他类别有较多的交集，与河流也存在部分交集。故在信息提取时会导致信息的误提。2.2.2空间特征空间特征是通过图像的像元值在空间上的变化反映出来的，包括图像上有实际意义的点、线、面或区域的空间位置、纹理、形状、大小、边缘和线性构造等，这些都属于空间特征。①纹理特征纹理又称质地，是遥感影像上重要的信息，是由于像片比例尺的限制，物体的形状不能以个体的形式明显的在影像上表现出来，而是以群体的色调、形状重复所构成的、个体无法辨认的影像特征。不同物体的表面结构特点和光滑程度是不一致的，在遥感影像上形成不同的纹理质地。纹理具有如下特征：a纹理具有局部的随机性和整体上的统计规律性；b复杂纹理是简单纹理的组合；c纹理是图像中某个区域共同的特征，对于单一的像元没有纹理可言；d不同的分辨率下纹理也是不同的。城市植被的纹理特征提取主要运用灰度共生矩阵、自相关法和小波变换等方法，常用的纹理特征值是熵。申广荣[2]等利用SPOT影像，以福建省泉州市鲤城区为例，在详细分析研究区植被不同类别纹理差异的基础上，利用统计和小波分析方法提取纹理信息，并将其引入基于光谱分类的最大似然法进行植被信息的提取分类，结果分类精度比监督分类提高了4个百分点。②边缘特征边缘是指周围像素灰度有阶跃变化或屋顶变化的那些像素的集合，存在于对象和背景、对象和对象、区域和区域之间。边缘总是存在于两个有不同灰度的区域之间，这是由于两个区域之间的灰度是不连续的。地物的轮廓又叫地物的边缘列表，或一条边缘列表的曲线模型。边缘特征可以分割图像，不同的图像对象灰度不同，边界处一般有明显的边缘特征。边缘特征也是纹理分析等其他图像分析的重要信息源和形状特征基础。城市植被的边缘特征是指图像上植被周围与其他地物灰度变化的分界，基本上体现了城市植被的轮廓。③形状特征形状是地物外貌轮廓在影像上的相似记录，任何物体都具有一定的外貌轮廓，在遥感影像上表现出不同的形状，如：游泳池是长方形，足球场则是两端为弧形的长方形，水渠为长条形，公路为蜿蜒的曲线型等。物体在影像上的形状细节显示能力与比例尺有很大关系，比例尺愈大，其细节显示愈清楚；比例尺愈小，其细节就愈不清楚，即地物形状根据比例尺在影像上的表现不同。但应当注意，遥感影像上所表现的形状与我们平常在地面所见的地物形状有所差异：a遥感影像所显示的主要是地物顶部或平面形状，是俯视图；b遥感影像是中心投影，物体形状在影像边缘会发生变形。城市植被在形态上比较规则，有矩形、环形、带状和零星点状等特点，如街道绿化带呈矩形或带状分布。3基于SPOT5影像的城市植被信息遥感提取3.1遥感图像的预处理本文研究的是城市植被覆盖度信息提取的方法，由于原始图像范围过大，包含的其他信息较多，故需要对原始图像进行裁剪，确定研究区范围。原始图像如图3所示，裁剪后图像如图4所示。图3原始图像图4裁剪图像裁剪后的图像是连云港城区的某一区域，从图中可以看出城市植被、道路、河流、建筑物和其他地类。其中，城市植被的分布较规则，多在建筑物周围、广场和道路两旁等地。3.2传统的植被覆盖度信息提取方法传统的植被覆盖度信息提取方法包括监督分类、非监督分类法、计算机分类新算法和目视解译等方法。目视解译法要求解译者运用遥感影像上不同地物所表现出的颜色、色调、形状、大小及纹理特征来判读影像，费时费力且精度不高。下面就只介绍监督分类、非监督分类和计算机分类新算法。3.2.1监督分类与非监督分类监督分类需要从研究区域选取有代表性的训练场地作为样本。根据已知训练区提供的样本，通过选取特征参数（如像元亮度均值、方差等），建立判别函数，据此对样本像元进行分类，依据样本类别的特征来识别非样本像元的归属类别。通过对监督分类的结果图与原图像的对比，可以看出植被信息误提严重，对河流的分类较好，道路提取有一定的信息缺失。非监督分类指在没有先验类别（训练场地）作为样本的条件下，即事先不知道类别特征，主要根据像元间相似度的大小进行归类合并（将相似度大的归为一类）的方法。常用方法有K-均值聚类法（K-mean）和动态聚类法（ISODATA）。通过非监督分类结果图与原图像的对比，可以看出植被、建筑物等信息提取的精度都不高，比监督分类的结果还要差一些。由于遥感图像普遍存在的“同物异谱”及“异物同谱”现象，使得非监督分类的精度很有限，而监督分类对训练场的选择存在很大主观因素，对于训练者不知或因数量太少未被定义的类别，监督分类也不能识别。所以无论监督分类还是非监督分类都存在很大的局限性。3.2.2植被指数法植被指数提取是遥感监测地面植物生长和分布的一种方法。它通过植被的光谱特征来说明植被的生长状况和覆盖范围。由于不同绿色植物对不同波长光的吸收率不同，光线照射在植物上时，近红外波段的光大部分被植物反射，而可见光波段的光大部分被植物吸收。通过对近红外和红光波段反射率进行线性或非线性组合，可以消除地物光谱产生的影响，得到的特征指数就是植被指数。目前，在科学文献中发布了超过150种植被指数模型，在这些植被指数中只有极少数经过系统的实验检验。植被指数虽然有好多种，但常用的植被指数有：归一