搜索
第２４卷第４期测绘工程Ｖｏｌ.２４，Ｎｏ.４２０１５年４月ＥｎｇｉｎｅｅｒｉｎｇｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇＡｐｒ．，２０１５基于ＥＮＶＩ的高分辨率遥感影像城市绿地信息提取研究陈阳，赵俊三，陈应跃（昆明理工大学国土资源工程学院，云南昆明６５００３３）摘要：近年来众多高分辨率遥感卫星的发射为城市绿地信息提取提供了高效的手段，如何充分利用高分辨率遥感影像对现代城市进行规划和生态环境评价具有重要意义。文中介绍一个完整的遥感工程在ＥＮＶＩ下的实现过程，包含高分辨率数据的处理、高分辨率信息提取流程的合理安排、ＦＸ模块面向对象信息提取等内容，可全面快速掌握城市绿地总量与分布现状，为绿地系统规划与日常管理决策提供参考，对生态城市建设具有重要意义。关键词：城市绿地；ＷｏｒｌｄＶｉｅｗ；信息提取；ＥＮＶＩ中图分类号：ＴＰ７５文献标志码：Ａ文章编号：１００６‐７９４９（２０１５）０４‐００３３‐０４ENVIbasedurbangreenspaceinformationextractionwithhighresolutionremotesensingdataＣＨＥＮＹａｎｇ，ＺＨＡＯＪｕｎ-ｓａｎ，ＣＨＥＮＹｉｎｇ-ｙｕｅ（ＳｃｈｏｏｌｏｆＬａｎｄＲｅｓｏｕｒｃｅｓＥｎｇｉｎｅｅｒｉｎｇ，ＫｕｎｍｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，Ｋｕｎｍｉｎｇ６５００３３，Ｃｈｉｎａ）Abstract：Ｔｈｅｌａｕｎｃｈｏｆｔｈｅｈｉｇｈｒｅｓｏｌｕｔｉｏｎｒｅｍｏｔｅｓｅｎｓｉｎｇｓａｔｅｌｌｉｔｅｐｒｏｖｉｄｅｓａｎｅｆｆｉｃｉｅｎｔｗａｙｆｏｒｕｒｂａｎｇｒｅｅｎｓｐａｃｅｉｎｆｏｒｍａｔｉｏｎｅｘｔｒａｃｔｉｏｎｒｅｓｅａｒｃｈｃｕｒｒｅｎｔｌｙ．Ｔｏｇｅｔｔｈｅｍｏｓｔｏｕｔｏｆｔｈｅｈｉｇｈｒｅｓｏｌｕｔｉｏｎｒｅｍｏｔｅｓｅｎｓｉｎｇｄａｔａｓｅｅｍｓｔｏｂｅｖｅｒｙｉｍｐｏｒｔａｎｔｉｎｍｏｄｅｒｎｕｒｂａｎｐｌａｎｎｉｎｇａｎｄｅｃｏｌｏｇｉｃａｌｅｎｖｉｒｏｎｍｅｎｔｅｖａｌｕａｔｉｏｎ．ＢｙｕｓｉｎｇａｃｏｍｐｌｅｔｅｒｅｍｏｔｅｓｅｎｓｉｎｇｅｎｇｉｎｅｅｒｉｎｇｅｘａｍｐｌｅｉｎＥＮＶＩ，ｉｔｇｉｖｅｓｔｈｅｉｎｔｒｏｄｕｃｔｉｏｎａｂｏｕｔｈｉｇｈｒｅｓｏｌｕｔｉｏｎｒｅｍｏｔｅｓｅｎｓｉｎｇｄａｔａｍａｎａｇｅｍｅｎｔ，ｒｅａｓｏｎａｂｌｅａｒｒａｎｇｅｍｅｎｔｓｏｆｅｘｔｒａｃｔｐｒｏｃｅｄｕｒｅｗｉｔｈｈｉｇｈｒｅｓｏｌｕｔｉｏｎｒｅｍｏｔｅｓｅｎｓｉｎｇｄａｔａ，ｏｂｊｅｃｔ-ｏｒｉｅｎｔｅｄｉｎｆｏｒｍａｔｉｏｎｅｘｔｒａｃｔｉｏｎｉｎＦＸｍｏｄｕｌｅ，ｅｔｃ．Ｔｈｕｓ，ｔｈｅｔｏｔａｌｑｕａｎｔｉｔｙｏｆｔｈｅｕｒｂａｎｇｒｅｅｎｓｐａｃｅａｎｄｔｈｅｐｒｅｓｅｎｔｓｉｔｕａｔｉｏｎｏｆｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｃａｎｂｅｍａｄｅｉｎａｎａｌｌ-ｒｏｕｎｄａｎｄｆａｓｔｗａｙ，ａｎｄｉｔａｌｓｏｐｒｏｖｉｄｅｓａｒｅｆｅｒｅｎｃｅａｂｏｕｔｇｒｅｅｎｓｐａｃｅｓｙｓｔｅｍｐｌａｎｎｉｎｇａｎｄｄａｉｌｙｍａｎａｇｅｍｅｎｔ，ｗｈｉｃｈｈａｓｉｍｐｏｒｔａｎｔｉｍｐｌｉｃａｔｉｏｎｓｆｏｒｔｈｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆｔｈｅｅｃｏｌｏｇｉｃａｌｃｉｔｙ．Keywords：ｇｒｅｅｎｓｐａｃｅ；ＷｏｒｌｄＶｉｅｗ；ｉｎｆｏｒｍａｔｉｏｎｅｘｔｒａｃｔｉｏｎ；ＥＮＶＩ城市绿地作为城市中唯一有生命的基础设施，在改善城市生态环境和人居环境方面起着积极的作用［１］，它是城市的氧源，更是电磁辐射、噪音及有害气体的良好吸收体［２-４］，在城市生态平衡中扮演着重要的角色，准确掌握其分布结构及变化规律，预测其发展走向，对于维护区域生态平衡，保护城市生态环境具有不可忽视的作用［５-６］。城市绿地的规划须先立足于对城市绿地现状的了解，传统的绿地调查采取实地测量与统计相结合的方法，它效率低收稿日期：２０１４‐１１‐０９基金项目：国家自然科学基金资助项目（４１１６１０６２）作者简介：陈阳（１９８９－），男，硕士研究生．下而且统计结果受人为影响。随着遥感技术的发展，近年来众多高分辨率遥感卫星的发射为城市绿地信息提取提供了高效的手段，如何充分利用高分辨率遥感影像对进行现代城市规划和生态环境评价具有重要意义［７-８］。高分辨率遥感技术为城市生态环境的监测方法和评价指标的改进提供了有力的支持。获得的城市绿地覆盖信息与其他资料结合可以进行城市绿地空间分布分析评价，为合理进行城市绿地规划提供决策支持［９］。高分辨率遥感技术具有很多特点：多传感器提供高分辨率全色波段及多光谱数据，通过图像融合方法，可以得到色彩信息丰富的・３４・测绘工程第２４卷高分辨率光谱数据；卫星控制技术，可使重访周期行融合处理，使融合后的影像既保留多光谱特征又控制在数天以内；使地物内部组成要素更丰富，能拥有较高的空间分辨率，如图２所示。ＥＮＶＩ中的更好反映空间信息、地物尺寸形状、地物之间的关融合方法有ＨＩＳ变换、Ｂｒｏｖｅｙ变换、乘积运算、ＰＣＡ系。ＷｏｒｌｄＶｉｅｗ-ＩＩ卫星是Ｄｉｇｉｔａｌｇｌｏｂｅ公司于变换和Ｇｒａｍ-Ｓｃｈｍｉｄｔ，融合方法的选用十分重要，２００９年１０月６日发射升空，运行在７７０ｋｍ高的太不同方法有不同的适用范围，同样的融合方法应用阳同步轨道上的高分辨率商业遥感卫星，它能够为在不同影像中所得结果也往往不一样。本文选择世界各地的商业用户提供满足其需要的高性能图适于最新高空间分辨率影像的Ｇｒａｍ-Ｓｃｈｍｉｄｔ方像产品。ＷｏｒｌｄＶｉｅｗ-ＩＩ在空间分辨率（０.５ｍ），多法，改进了ＰＣＡ中信息过分集中的问题，而且不受光谱成像（１个全色通道、８个多光谱通道）、侧摆速波段限制，可以较好地保持影像的纹理和光谱信息。度、扫描宽度等方面具有显著的优势，能够满足更广领域的遥感用户，为其提供更好的遥感信息源服务，在数字城市建设、自然资源的管理和监控、自然灾害的评估及精细农业等领域中具有应用潜力［１０］。本文介绍了使用ＥＮＶＩ遥感图像处理软件，采用ＷｏｒｌｄＶｉｅｗ-２高分辨率遥感影像，对城市绿地信息提取的方法，其技术路线如图１所示。图1技术路线1实验数据本文选择成像时间为夏季，植被长势较好的ＷｏｒｌｄＶｉｅｗ-２含ＲＰＣ信息的ＬＶ２Ａ级数据，该数据包括０.５ｍ分辨率的全色影像和１.８ｍ分辨率的多光谱影像。其中多光谱影像包含８个波段，除含有蓝、绿、红、近红４个业内标准波段外，还含有海岸、黄、红边、近红４个彩色波段，增加的波段信息为用户提供进行精确变化检测和制图的能力，其中非可见光近红外波段，波长较长，受大气影响较小。因红和近红外波段是对绿地敏感的波段，为了提高工作效率，故本文在图像分类时仅选择４个业内标准波段参与图像分割。2数据预处理2.1图像融合将高空间分辨率的全色影像与多光谱数据进2融合后图像与多光谱图像对比ＥＮＶＩ中分别打开ＷｏｒｌｄＶｉｅｗ-２多光谱图像和全色图像，选择Ｇｒａｍ-ＳｃｈｍｉｄｔＰａｎＳｈａｒｐｅ‐ｎｉｎｇ，选择适当的空间裁剪及波谱裁剪然后进入融合参数面板，该面板会自动识别其传感器类型，重采样方法选择精度最高的三次卷积。2.2自定义坐标系为避免坐标来回转换时产生较大误差，应首先统一控制点和输出正射校正图像的坐标系。本文因ＤＯＭ数据为北京５４坐标系，所以统一坐标系为北京５４坐标系。该步骤在ＥＮＶＩＣｌａｓｓｉｃ中进行，从网上获取北京５４坐标系相应投影、基准面、中央经线长度比等参数，将其添加至软件中并保存。2.3正射校正在卫星影像和航空影像中会有一些几何误差，主要由比例尺变化、传感器姿态的变化、传感器的系统误差引起。对于分辨率较高且具有ＲＰＣ信息的图像可以用正射校正的方法完成几何校正以消除这些误差达到更高的精度要求。在本文中，采取ＧＣＰ＋ＲＰＣ＋ＤＥＭ的方法完成正射校正，其中ＧＣＰ控制点从参考图像ＤＯＭ中选择。首先将参考图像ＤＯＭ和ＤＥＭ数据进行绑定，这样可以在ＤＯＭ上获取x，y及该点的高程值。在选控制点时，不能选择建筑物的顶点，因为传感器若有倾斜角度，建筑物顶点将会偏移，所以确保选点时选择地面上的控制点。选好同名点（链接第４期陈阳，等：基于ＥＮＶＩ的高分辨率遥感影像城市绿地信息提取研究・３５・点）后，在参考图像中右键选择像元位置，可将参数导出至控制点选择面板并添加点。添加点完成之后可查询控制点的误差，若太大可进行删点重选。因图像融合时丢失了ＲＰＣ信息，进行正射校正时需要选择相应ＲＰＣ文件，可以使用全色影像的ＰＲＣ信息进行补充。2.4快速大气校正因提取绿地信息所需精度不用太高，故本文使用大气校正扩展模块中的快速大气校正工具（ＱＵｉｃｋＡｔｍｏｓｐｈｅｒｉｃＣｏｒｒｅｃｔｉｏｎ），可以提高ＮＤＶＩ和光谱属性值的精度，还可以消除部分大气影响。完成大气校正后定位到绿地并分别打开大气校正前后的光谱曲线，发现校正后的光谱曲线更为真实，如图３所示。图3快速大气校正前后绿地波谱曲线对比3面向对象绿地信息提取绿地信息提取使用到特征提取（ＦｅａｔｕｒｅＥｘ‐ｔｒａｃｔｉｏｎ，ＦＸ）这个扩展模块，ＥＮＶＩ５.１中该模块含有３个工具：基于样本选择的特征提取、基于规则的特征提取与分割式特征提取。提取植被选用基于规则的特征提取（ＲｕｌｅＢａｓｅｄＦｅａｔｕｒｅＥｘｔｒａｃｔｉｏｎ）。归一化植被指数（ＮｏｒｍａｌｉｚｅｄＤｉｆｆｅｒｅｎｃｅＶｅｇｅ‐ｔａｔｉｏｎＩｎｄｅｘ，ＮＤＶＩ）计算可以将多光谱数据变换成一个单独的图像波段，用于显示植被分布。较高的ＮＤＶＩ值预示着包含较多的绿色植被。ＥＮＶＩ中的ＮＤＶＩ使用标准算法ＮＤＶＩ＝ＮＩＲ－ＲＥＤＮＩＲ＋ＲＥＤ（ＮＩＲ：近红外波段，ＲＥＤ：红波段）．ＮＤＶＩ值的范围为－１～１［１１］。3.1调整分割尺度与合并尺度创建发现对象时，分割尺度越大则分割的块越大，值越小则分割的块越多，勾选预览可实时预览分割效果图，根据提取地物的大小选择适合的分割尺度。分割时，发现分割的对象相对零碎，使用合并尺度对其进行合并，合并的值越大被合并的块越多。通过预览和经验调整分割尺度值为３５，合并尺度值为６０时效果较为理想。3.2建立规则建立规则时需建立两种规则，一种是把需要的地物提取出来，一种是把不需要的地物去掉。新打开一个ＥＮＶＩ工程文件并打开经过处理过的研究区图像，以ＲＧＢ标准假彩色显示波段使绿地植被显示为红色，利于目视解译找绿地植被。使用ＲＯＩ（感兴趣区域）选择感兴趣区域时，可在研究区域各个地方，选择一些长势不太好的植被进行标注，注意在建筑物阴影中的植被（在图像上显示为近黑色但确定为植被）不要标注，若选择会使结果数值很小，失去研究意义，统计可得绿地对应的ＮＤＶＩ的最小值为０.１４９５１０。回到ＦＸ模块预览窗口设置该值发现有部分显示为紫色屋顶的地物被划分进绿地里面，故再次建立一个ＡＮＤ的规则以去除紫色屋顶。新建感兴趣类型，放大窗口可参考之前预览时被错分的地物目视选择出一些覆盖人造涂料的地物作为样本，统计得出在绿地区域Ｂａｎｄ１的最大值与紫色屋顶的最小值基本接近，为１２８１～１３８１，如图４所示。故使用＜１３８１作为剔除紫色屋顶的阈值。再次预览发现紫色屋顶被剔除。由此获得比较理想的面向对象分类规则［１２］如下：ＳｅｇｍｅｎｔＳｃａｌｅＬｅｖｅｌ：３５.０ＭｅｒｇｅＬｅｖｅｌ：６０.０・３６・测绘工程第２４卷Ｒｅｆｉｎｅ：ＮｏＴｈｒｅｓｈｏｌｄｉｎｇＱＵＡＣ（Ｏｒｔｈｏｒｅｃｔｉｆｉｅｄ（Ｂａｎｄ５）：ＷＶ２＿ＧＳ＿ＡｔｔｒｉｂｕｔｅｓＣｏｍｐｕｔｅｄ：ＯＲ＿ＱＵＡＣ．ｄａｔ）（８３２.５０００）ＳｐｅｃｔｒａｌＣｌａｓｓｉｆｉｃａｔｉｏｎ：Ｒｕｌｅ-ＢａｓｅｄＢａｎｄＲａｔｉｏＲｕｌｅＳｅｔ：ＱＵＡＣ（Ｏｒｔｈｏｒｅｃｔｉｆｉｅｄ（Ｂａｎｄ４）：ＷＶ２＿ＧＳ＿ｂａｎｄｒａｔｉｏ＞０.１４９５ＯＲ＿ＱＵＡＣ．ｄａｔ）（６５８.８０００）ａｖｇｂａｎｄ＿１＜１３８１.００００图4绿地样本与非绿地样本的统计结果4结束语像道路提取［Ｊ］．测绘工程，２０１４，２３（３）：４２-４５．本文介绍了一个完整的遥感工程在Ｅ