基于TM影像的水体信息提取实习报告

1、精选优质文档-倾情为你奉上目录实习目的1、熟悉ENVI软件的具体安装和操作；2、理解并能掌握水体信息的提取方法、步骤及其结果评价方法；3、掌握遥感专题图的制图输出；4、通过实习，加深对所学遥感专业知识的理解与掌握，使理论知识融会贯通，提高实践动手能力和综合解决实际问题的能力。实习内容一、操作环境软件：ENVI 4.7影像：某地区的TM影像 二、常用的水提取方法1、非监督分类法也称聚类分析，是一种“盲目”的分类；仅凭遥感影像地物的光谱特征的分布规律，即自然聚类的特性，进行分类；2、监督分类法根据已知训练区提供的水体样本，通过选择特征参数，求出特征参数作为决策规则，建立判别函数以对各待分类影像进行

2、的图像分类；3、单波段阈值法 单波段阈值法主要利用对水陆界线反映较好的TM影像的第5波段（中红外波段）, 根据影像的灰度特征经过数据采样确定其阈值。进行水体的提取。提取模型如式1所示, 其中TM5表示第五波段的灰度值, T为水体提取的灰度阈值: TM 5 TM4+TM5 (2)5、归一化水体指数法水体指数法的基本原理是通过波段之间的比值运算增强某些地物之间的反差。考虑水体的波谱特征是在TM2具有较强的透射性,TM5处于水的吸收带内, 反映含水量敏感, 因此采用TM2波段与TM5波段数值之差和这两个波段数值之和的比值作为规一化水体指数(NDWI) , 当满足NDWIC时所提取的地物为水体。归一化

3、水体指数计算式如式3所示:NDWI = (TM2- TM5) / (TM2+TM5) (3)三、具体操作步骤及效果比较1、非监督分类法 ClassificationUnsupervisedIsoData选择实验影像弹出对话框 ISODATA Parameter（选择最多分10类，循环10次）生成非监督分类图像（通过类别的属性设置，使水体高亮显示 黑色为其他地物，白色为水体，绿色为河岸、人造水塘等浅水体） 原图像 局部 非监督分类图像局部2、监督分类法 1）建立模板： 原图像的Overlay菜单 Region of Interest 弹出对话框 ROI Tool，选择样本 2）评价分类模板：对话

4、框 ROI Tool Option Compute ROI Separability 弹出Select Input File对话框（选择原图影像） 弹出ROI Separability Cacilation对话框（选择要评价的ROI） 生成 ROI Separability Report。 注:大于1.8即为合格 3）进行监督分类 ClassificationSupervisedMaxinum Likelihood选择实验影像生成监督分类图像 监督分类图像 局部图像3、单波段阈值法3.1步骤1）确定阈值 原图像的Overlay菜单 Region of Interest 弹出对话框 ROI To

5、ol 选择样本 确定阈值确定阈值为21；2）进行提取 Basic Tools Band Math （输入提取模型 b5 b4+b5 ） 生成谱间关系法提取的图像。 谱间关系法图像 局部4.2 分析杨存建提出的水体具有(TM2+TM3)(TM4+TM5)的特征的方法易操作、稳定性好。本文采用此方法经试验，得到如图所示的谱间关系法水体提取效果图。用ENVI 对水体提取结果影像和原影像进行叠加判读，谱间关系法提取了几乎所有的水体区域，但是，此方法存在把阴影误提为水体的问题。5、水体指数法5.1 步骤1）建立NDWI模型： Basic Tools Band Math （输入 （b2-b4）/（b2+b

6、4） 生成NDWI的图像。 NDWI图像 局部2）选择水体样本：确定水体的NDWI范围 NDWI图像的Overlay菜单 Region of Interest 弹出对话框 ROI Tool 选择样本 确定范围确定范围为： max=0.49 min=0.223）波段运算输出提取图像： Basic Tools Band Math（输入（b10.22） 生成NDWI的图像 水体指数图像某局部1 某局部 2 某局部 35.2 分析改进的归一化差异水体指数(NDWI)法提取漏提了部分水体（如局部3），主要也是因为水体被污染（如山洪爆发使水体混入了泥沙）后和非水体在影像上表现为差别不大，由于确定的范围问题

7、，部分水体（如局部2）中间部分有“椒盐噪音”现象，可以通过分类后图像处理解决（如腐蚀、膨胀）。四、水体提取的综合 1、水体及背景地物的光谱值分析我们对易混淆的地物进行了采样统计，做出水体与易混淆地物的光谱特征曲线图，分析了以上方法产生误提和漏提的原因，并观察地物的光谱特征规律，研究如何综合各方法的优势，取得好的提取效果。首先，在遥感影像上，就每种地物选取一定的研究样本，测定各波段的光谱亮度值，并从中抽取一些典型而具有代表性的样本。统计后平均值数据列于表1中，并对同一地物同一波段的采样点亮度值做均值统计，做出如图所示的地物光谱特征曲线图。TM 2TM 3TM 4TM 5水体72601729河岸等

8、浅水体74765261阴影43352134表1地物光谱特征曲线图从地物光谱特征曲线图中可以看出由于在第五波段阴影和水体的亮度值较接近，单波段阈值法无法将水体单独提取出来；水体、河岸及浅水体和阴影在第2、3、4、5波段上均有(TM2+TM3)(TM4+TM5)的光谱特性，所以用光谱间关系法提取水体时误提了很多东西出来；河岸及浅水体TM2和TM5采样点的亮度值比较接近，所以用改进归一化差异水体指数法会提取出部分河岸及浅水体。同时，阴影的亮度值也具有TM2TM5的特征，导致阴影被误提为水体。2、水体提取的解决方案分析了上述3种方法的提取效果，并结合地物的光谱特征曲线图，我们得出的解决方法如下 (1)

9、 如何解决阴影的误提： (TM2+TM3)(TM4+TM5)特征对消除易混淆的阴影的效果较好，但对部分阴影，(TM2+TM3)(TM4+TM5)影像仅起到了增加水体和阴影亮度差值的作用。用谱间关系法和改进的归一化差异水体指数(NDWI)法都不易消除阴影的影响。因此，可以用T0(TM2+TM3) -(TM4+TM5) T1 (T1为阈值)的方法来解决。(2) 如何解决河岸及浅水体的误提：光谱特征曲线图第五波段水体和阴影的亮度值差异不大，而河岸及浅水体的亮度值与水体和阴影的亮度值差异很大。因此可采用TM5T2 (T2为阈值)来消除河岸及浅水体的影响。通过上述分析得到水体提取综合方法，对影像进行处理，流程如图所示，且得出水体提取效果图，如提取效果图所示。T0(TM2+TM3) -(TM4+TM5) T1TM 5 T2水体提取效果图流程图取T0为61，T1为133，