精品基于遥感影像提取水体信息用于水文研究可行性分析可研报告最新.docx

精品基于遥感影像提取水体信息用于水文研究可行性分析可研报告最新.docx

《精品基于遥感影像提取水体信息用于水文研究可行性分析可研报告最新.docx》由会员分享，可在线阅读，更多相关《精品基于遥感影像提取水体信息用于水文研究可行性分析可研报告最新.docx（27页珍藏版）》请在冰点文库上搜索。

精品基于遥感影像提取水体信息用于水文研究可行性分析可研报告最新

（此文档为word格式，下载后您可任意编辑修改！

）

基于遥感影像提取水体信息用于水文研究可行性分析研究报告

（课程作业）

一、研究目的

熟练地掌握遥感图像处理技术，结合遥感图像与水文学的联系，解释和解决一些水文现象和水文问题。

通过水体在遥感图像上表现的光谱特征提取遥感图像中水体，并与目视解译下的遥感影像中的水体及用数字高程模型-提取流域的进行对比，分析遥感影像提取水体的精度，从而论证从遥感影像中直接提取水体应用于水文研究的可行性。

对比遥感影像解译、遥感影像数据水体信息提取和根据数字高程模型提取的影像之间的对比。

分析各种水体之间的差异和精确度分析。

二、重点难点

（1）重点：

从遥感影像提取水体：

选择适当的方法对水体的提取，关系到提取水体的准确性，因此选择正确恰当的方法对水体进行提取是本文的重点。

（2）难点：

晋江流域边界的切割。

流域范围是指河流、湖泊的集水区，在地形图上一般表现为山脊或分水岭。

如何在遥感影像中切出研究流域的范围是本文的难点。

（3）DEM数据提取河道信息也是新增的重点。

三、创新点

河流、湖泊等水体的演变和发展规律是水文学研究的重要内容，地形，植被、降水、土壤等因素均对水文要素的产生深刻的影响。

目前对遥感影像的信息提取主要依赖于目视解释。

对遥感影像中单一的水体研究而言，采用地理信息技术对遥感影像水体的直接提取不仅可以节约大量目视解释的劳动力成本，也可以提高水体精度，减少一些不必要的误差。

在对水文、水系演化的研究中采用地理信息技术对水体的提取具有独一无二的优势。

四、原始数据和遥感影像

（1）团队已有的TM2006年11月份遥感影像（图1）

（2）从国际科学数据平台下载遥感影像:

由于2006年没有查询到条带号为120042的遥感影像，因此下载在2004年的相关影像进行替换。

（3）晋江附近各地区30数字高程模型：

图

（1）TM2006年11月份遥感影像

五、主要操作流程

（1）对遥感图像的初步处理。

（2）流域边界的提取：

利用晋江流域数字高程模型根据分水岭提取流域边界。

（3）晋江流域遥感图像水体的提取：

利用水体在遥感图像上表现的光谱空间特征和纹理特征对其进行信息提取。

（4）基于DEM影像提取水系信息。

（5）建立解译标志，解译晋江流域土地利用状况。

（6）将晋江流域遥感图像水体提取以及水系信息与目视解译水体进行对比，分析对流域内水体提取的精度。

具体方法思路可见图2（线路流程图）

图

（2）研究线路流程图

六、主要操作步骤

（一）晋江流域水体信息的提取

1.1研究方法

归一化水体指数法

水体的波谱特征是在TM2具有较强的透射性，TM5处于水的吸收带内，反映含水量敏感，因此采用TM2波段与TM5波段数值之差和这两个波段数值之和的比值作为归一化指数（NDMI），当NDMI>0时提取的地物为水体。

NDMI=（TM2-TM5）（TM2+TM5）

1.2.组合多波段数据

方法一、采用PCI软件将多波段遥感影像的进行融合图（3），并将融合影像转化为.img格式，在转化为.img格式时出现以下提示图（4），Pix文件不能直接转化为.img格式。

因此只能先转化.tif文件再通过ENVI软件转化为.img文件。

图（5）为融合之后的遥感影像。

图（3）

图（4）

图（5）组合之后的遥感影像

方法二、由于用PCI组合多波段数据比较复杂，耗时较长，因此后来我发现使用erdas速度较快也更为简洁。

120042影像采用erdas进行各波段组合，主要包括以下步骤：

（1）在ERDAS图标面板工具条中，点击——打开输入输出对话框，如图所示。

并做如下的选择：

（2）为了图像处理与分析，需要将上述转换的单波段IMG文件组合为一个多波段图像文件。

在ERDAS图标面板工具条中，点击Interpreter|Utilities|LayerStack——LayerSelectionandStacking的对话框。

在LayerSelectionandStacking对话框中，依此选择并加载（Add）单波段IMG图像：

点击OK执行并完成波段组合。

1.3三副影像的拼接

由于ERDAS对遥感影像拼接出现故障，因此转由ENVI进行三副影像的拼接,主要分为以下步骤：

（1）选择主菜单File→OpenImageFile，打开要进行拼接的图像

（2）单击主菜单Map→Mosaicking→Georeferenced，打开图像镶嵌窗口；

（3）在MapBasedMosaic窗口中单击Import→ImportFiles命令；

（4）在弹出的MosaicInputFiles对话框中选择待拼接的图像文件

（5）点击OK，把图像加载到了图像镶嵌的窗口中；

（6）在图层列表栏中选择需要调整重叠次序的图层，右键点击选择RaiseImagetoTop或者RaiseImagetoPosition进行重叠次序的调整；

（7）在图层列表栏中选择需要调整重叠次序的图层，右键点击选择EditEntry对图像镶嵌参数进行设置；

（8）点击File菜单下的Savetemplate命令对图像进行虚拟镶嵌；

（9）点击File菜单下的Apply命令，打开镶嵌图像保存对话框，设置输出的像元分辨率，重采样方法以及输出文件名等参数，点击OK完成图像的镶嵌和保存（图2-31）。

（10）点击OK进行图像的镶嵌（图2-32）。

以下为三幅影像的拼接结果：

1.4遥感影像的裁剪

由于拼接之后遥感影像过大，运算时间过长，因此才有先裁剪后提取水体办法，

（1）利用DEM数据确定晋江流域的边界，并在ERDAS里面完成数据的裁剪。

由于这张遥感影像的面积不完全包括晋江流域的面积，因此晋江流域没有被完整的裁剪下来，与北部拼接之后，在西北不出现遥感影像的坏道，因此要与北部和西部的那幅遥感影像将进行拼接才能裁剪出晋江全流域。

图（6）晋江流域平面图图（7）裁剪后的不完整的晋江流域遥感影像

图（8）裁剪后的不完整的晋江流域遥感影像图（9）裁剪后的完整的晋江流域遥感影像

从上图可以看出三副裁剪之后的遥感影像，图7为一张行列号位119042和遥感影像裁来的完整的遥感影像。

1.4遥感影像进行水文信息的提取

使用ERDAS软件构建模型，将已经裁剪好的遥感影像进行水文信息的提取，分析提取精度。

图（10）为归一化水文信息指数模型，通过归一化水文指数的提取可以将遥感影像中的水体提取出来。

图10归一化水文信息指数模型

下列图均为提取结果：

图11可以看出整个流域水体以晋江为主体，其他小块水体零星散落的其周围。

图12、13为提取之后的山美水库，可以很清楚的看出山美水库水体的分布，即使小岛的很清楚，具有较高的准确性，能够较好的反应水体信息。

图11流域提取结果以及与原地物的关系

图12山美水库提取水体与真实水体之间的对比

图（13）提取之后的水体信息与原始图像对比

图（14）小流域水体信息的提取

但是从图13可以看出那些零星的水体在提取后面积明显更小，在较细小的水体提取之后成为断断续续的间断的水体图（14），与现实的水体信息具有较大的差异，因此精度有待提高。

1.5初步结论

1.通过归一化水体指数提取的水体信息，能够较准确的反应遥感中大面积的水体信息。

2.但是在在细小的水体时所提取的水体面积比实际值较小，不能反映部分小地区水体信息，精度有待提高。

（二）.基于DEM晋江流域河网的提取

利用ArcGis水文分析模块生成河网主要有以下几个步骤：

DEM洼地填平、水流方向矩阵生成、水流累计矩阵生成、设定集水面积阈值形成河网矩阵、河网节点链接信息矩阵、河网分级矩阵、矢量河网（图15）。

主要采用Arctoolbox里面的水文分析工具（Hydrology），由于下载时DEM影像一个经纬度一张图像因此下载范围比较大，

图15基于ArcGis水文分析模块的矢量河网提取

2.1各DEM影像之间进行拼接

由于下载的DEM数据数一个经纬度为一张图，因此要把晋江流域周边的影像进行拼接形成一张统一的图形文件。

由于两块在海上因此DEM为空白。

没有可供下载的影响。

2.2提取水流方向

采用以下操作:

ArctoolBox->SpatialAnalysisTools->FlowDirection工具

结果为：

2.3计算洼地

采用Hydrology下面的WATERSH进行运算

结果为：

洼地

洼地贡献区域

2.4填充洼地计算无洼地的DEM影像：

DEM洼地（水流积聚地）有是洼地和数据精度不够高所造成的洼地。

洼地填平的主要作用是避免DEM的精度不够高所产生的（假的）水流积聚地。

洼地填平使用ArctoolBox->SpatialAnalysisTools->Hydrology->Fill工具。

2.5、计算汇流累计量和生成河网

生成的河网

2.6晋江流域河网的裁剪

通过分析ArcGis软件中的水文分析模块StreamtoFeature工具生成的矢量河网上下节点标识的特点，在由几何网络工具中生成的节点图层中将其计算出来。

这不仅完善了地理数据一致性，也为后续的流域水文地貌过程分析提供了保障。

2.7裁剪结果和分析

下图均为晋江流域裁剪之后的结果，图16可以看出由于晋江流域较大总体上看没有清晰的河网信息，只是一些零散的细点，具体来看图17、18河网脉络清晰可见

此外，由于在DEM影像拼接是还剩下一些缝隙（如图19）为不足之处。

由于提取水系是已经将地表的洼地填充一些集水区如湖泊等已经被填充，无法反映真实的水体分布。

总之，基于DEM提取的水是能够较清楚的放映水系的分布和水的流向，但是所提取的河道，大小均一不能反映实际的水系状态，也没有湖泊等相关的水体信息。

图16晋江流域提取的河网信息图17

图18图19

（3）目视解译

3.1目视解译结果

目视解译为本报告前期准备工作，其包括遥感遥感影像的初步处理、建立解译标志，解译后实地考察等工作。

因此在此不做详细介绍。

本次对晋江流域的解译把地物分为个地类，包括：

水田、旱地、水域、林地、园地建设用地、未利用地。

下图为解译的结果

图20图21

图22晋江流域遥感目视解译结果

3.2解译结果分析

从图20可以看出目视解译可以较为清晰的反应水体信息，信息较为客观精确，不过解译工作量巨大。

图21可以看出对于相同颜色的地区是否为水体，难以判断。

在这种情况下主要取决于解译对地物的主观判断，和解译经验的积累。

从22图可以看出在水库的下游出水口应该有一定量的水体，但是在影像上不明显，解译就会出现水系断裂的情况，不能反映真实的水文系统。

总之，目视解译水体对部分存在一定的主观性，无法判别的地方与解译人的主观意识和解译经验有关。

对于细小的河道是没有办法进行解译，没有显示出具体完整的水系特征。

七结论与存在的问题

1.通过归一化水体指数提取的水体信息，能够较准确的反应遥感中大面积的水体信息。

但是在在细小的水体时所提取的水体面积比实际值较小，不能反映部分小地区水体信息，精度有待提高。

2.基于DEM提取的水是能够较清楚的放映水系的分布和水的流向，但是所提取的河道，大小均一不能反映实际的水系状态，也没有湖泊等相关的水体信息。

3.遥感影像解译，比较精确，部分地区存在主观判断，不同人解译存在人为误差，工作量大。

目视解译水体对部分水体存在一定的主观性，无法判别的地方与解译人的主观意识和解译经验有关。

对于细小的河道是没有办法进行解译，没有显示出具体完整的水系特征。

8、后续工作

从上面工作可以看出无论是目视解译还是归一化指数方法对水体的提取以及和网信息的提取都存在一些问题。

要真正将提取的水体信息运用于水文分析还有一定的差距。

后续还有一些工作需要开展。

初步为一下打算：

1.对于河流下游、湖泊以及较为明显的水体可采用归一化水体指数提取方法提取的水体信息结果。

调整阀值使水体信息能够最大限度的放映实际水体分布。

2.对于一些无法提取和较细小的河流可以采用河网提取的结果，由于网提取的结果是线性不存在面积，因此可以根据地形，降水和地理位置，差异的赋予不同地区的河流一定的缓冲区。

具体方法还有待进一步不学习和探讨。

3.基于于晋江流域水体信息提取的这些研究以及后续研究愿能够真正提取出对水文学真实有用的水体信息，并以论文形式发表。

九、参考文献

[1]王亚飞,刘文军.面向对象的遥感影像水体信息的提取[J].河南科学,2012,30（5）:

于欢，张树清，李晓峰，等．基于TM影像的典型内陆淡水湿地水体提取研究[J]．遥感技术与应用，2008（3）：

310—3l5．

[2]陈树辉，刘忠华，刘翔，等．基于TM遥感影像的南京市水体信息提取[J]．安徽王培培．基于ETM影像的水体信息自动提取与分类研究[J]．首都师范大学学报：

自然科学版，2009（6）：

75-79

[3]胡进刚，张晓东，沈欣，等．一种面向对象的高分辨率影像道路提取方法[J]．遥感技术与应用，2006（3）：

184—188．

[4]汤国安,杨昕.ArcGIS地理信息系统空间分析实验教程[M].北京:

科学出版设,2006:

[5]韩耀军,王乐鹏,秦天保等.数据库系统原理与应用[M].北京:

机械工业出版社,2007:

[6]张行南，叶立华，井立阳．基于DEM的流域流水网对比分析[J]．水利水电科技进展，2004，24（3）：

1-4

[7][2]文康，金管生．地表径流过程的数学模型[M]．北京：

水利电力出版社，1990：

288—289

[8][31任立良．流域水文物理过程的数字模型研究[D]．南京：

河海大学，1999：

90—112

[9]易红伟,汤国安,刘咏梅等.河网径流节点及其基于DEM的自动提取[J].水土保持学报.

毕业设计（论文）原创性声明和使用授权说明

原创性声明

本人郑重承诺：

所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。

尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。

对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。

作者签名：

　　　　　日　期：

指导教师签名：

　　　　　日　　期：

使用授权说明

本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：

按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。

作者签名：

　　　　　日　期：

学位论文原创性声明

本人郑重声明：

所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。

除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。

对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。

本人完全意识到本声明的法律后果由本人承担。

作者签名：

日期：

年月日

学位论文版权使用授权书

本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。

本人授权　　　　大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。

涉密论文按学校规定处理。

作者签名：

日期：

年月日

导师签名：

日期：

年月日

注意事项

1.设计（论文）的内容包括：

1）封面（按教务处制定的标准封面格式制作）

2）原创性声明

3）中文摘要（300字左右）、关键词

4）外文摘要、关键词

5）目次页（附件不统一编入）

6）论文主体部分：

引言（或绪论）、正文、结论

7）参考文献

8）致谢

9）附录（对论文支持必要时）

2.论文字数要求：

理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。

3.附件包括：

任务书、开题报告、外文译文、译文原文（复印件）。

4.文字、图表要求：

1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写

2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。

图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画

3）毕业论文须用A4单面打印，论文50页以上的双面打印

4）图表应绘制于无格子的页面上

5）软件工程类课题应有程序清单，并提供电子文档

5.装订顺序

1）设计（论文）

2）附件：

按照任务书、开题报告、外文译文、译文原文（复印件）次序装订