植被指数、变化检与图像增强的算法实现.pdfVIP

植被指数、变化检测与图像增强的算法实现 武汉大学遥感信息工程学院 20113025900XX XX 【摘要】通过为期一周的实习，笔者思考研究了计算RVI、ND、DVI等植被指数， 基于像素光谱的变化检测，均值滤波、中值滤波、高斯平滑等图像增强以及利用 不同算子进行边缘检测的算法原理与程序实现。通过自主选择实验参数、确定实 验方案、设计程序，解决过程中的问题，直至得出正确的结果和结论，激发了创 新思维，提高了创新能力和实践能力。本文就以上所列算法进行了阐述和展示。 【关键词】植被指数 变化检测 图像增强 边缘检测 算法实现 一、原理介绍 （一）植被指数 植物叶面在可见光红光波段 （R）有很强的吸收特性，在近红外波段 （NIR) 有很强的反射特性，这是植被遥感监测的物理基础，通过这两个波段测值的不同 组合可得到不同的植被指数。 RVI=NIR/R RVI 1 比值植被指数 ，绿色健康植被覆盖地区的 远大于 ，无植 被覆盖的地面（裸土、人工建筑、水体、植被枯死或严重虫害）RVI在1附近。 归一化植被指数NDVI=(NIR-R)/(NIR+R) ，-1=NDVI=1，负值表示地面 0 NIR R 覆盖为云、水、雪等，对可见光高反射； 表示有岩石或裸土等， 和 近似 相等；正值，表示有植被覆盖，且随覆盖度增大而增大。 差值环境植被指数DVI=NIR-R，正值表示植被覆盖，随覆盖度增大而增大。 SAVI=((NIR-R)/(NIR+R+L))(1+L) L 调整土壤亮度植被指数 ， 是根据实际情 况确定的土壤调节系数，取值范围0~1。L=0时，表示土壤背景的影响为零，即 植被覆盖度非常高，土壤背景的影响为零，这种情况只有在被树冠浓密的高大树 木覆盖的地方才会出现。 （二）基于像素光谱的变化检测 变化检测是从不同时期的遥感数据中定量分析和确定地表变化的特征与过 程，先获得两幅通用地点不同时间图像的差异图像，再对差异图像进行处理，将 像素点分成变化和无变化两类。笔者采用了差值和比值两种方式。当差值=0或 比值=1时，表示像素点无变化，否则像素点有变化。 （三）图像增强 均值滤波是一种常用的线性滤波算法，用以去除噪声，它是指在图像上对目 标像素给一个模板，该模板包括了其周围的临近像素 （以目标象素为中心的3*3 5*5 7*7 或 或 个像素，构成一个滤波模板，但不包括目标像素本身），再用模板 植被指数、变化检测与图像增强的算法实现 20113025900XX XX 中的全体像素的平均值来代替原来像素值。 中值滤波与均值滤波类似，区别在于中值滤波的模板包括目标像素本身，且 用模板内全体像素的灰度中值来代替原来像素值。 高斯平滑采用高斯滤波为影像消除噪声，高斯滤波中每一 个像素点的值，都由其本身和邻域内的其他像素值经过加权平 均后得到。笔者采用的是如右图所示的3*3算子。 梯度锐化用来突出图像的边缘和轮廓信息，最简单的梯度算子就是分别求出 目标像素在x和y方向上的灰度值变化率，取其中较大的作为该点梯度的大小与 方向，自定义一个阈值，分梯度大于和小于或等于该阈值两种情况，对原像素值 进行变换。 （四）边缘检测 Sobel Robert Prowitt Laplace Krisch 笔者采用了 、 、 、 、 五种算子对该算法 进行了实现，以Sobel算子为例，取下图所示的一组模板，分别对目标像素进行 计算，取其中较大的作为边缘强度的大小与方向，并自定义一个阈值，若大于该 255 0 阈值，则原像素值变为 ，否则变为 ，这样就生成了一幅二值图像，白色的 边缘部分得到充分突出。