modis影像水体信息提取实验报告.docxVIP

modis影像水体信息提取实验报告 　　基于高分一号卫星多时相数据的洪水监测　　摘要：本文利用两幅高分一号多光谱影像数据，通过软件，采用NDVI对黑龙江地区水体进行了提取，并在图像上展示了水体变化区域，计算了水体变化面积。结果表明：9月9日黑龙江水域面积比8月27日增加了。最后又采用了假彩色合成法展示了水体增加区域。结果表明：两种方法对水体变化信息的提取具有一致性。　　1数据介绍　　本作业获得了两幅高分一号TIF数据，分别是8月27日，9月9日。每幅影像有4个波段，查阅资料得知：1波段波长为，属于蓝、青光，2波段波长为，属于黄、绿光，3波段波长为，属于红光，4波段为，属于近红外。　　图10827影像信息图20909影像信息2研究区域由所给数据的经纬度坐标可知，研究区域为抚远县，其地处黑龙江、乌苏里江交汇的三角地带。地理方位是东经133°40′08″至　　135°5′20″,北纬47°25′30″至48°27′40″。　　图3研究区域的百度卫星地图　　2水体提取方法选择　　单波段：水体在近红外波段的反射率很低，所以可以设置阈值进行提取。　　归一化水体指数NDWI?(?Green??NIR)/(?Green??NIR)　　归一化植被指数NDVI?(?NIR??Red)/(?NIR??Red)　　但单波段方法中阈值的设置需要反复调整，而高分一号数据的1、2波段不完全是蓝、绿光，而3、4波段完全是红、近红外。所以选择归一化植被指数提取水体。-1=0，则DN值为2；其他为0；　　结果如下：　　图7两幅水体二值图的差值　　图中有三种DN值，0表示黑色背景；1表示8月27日存在，而在9月9日不存在的水体；2表示在8月27日不存在，而在9月9日存在的水体。　　4结果分析　　、将两幅水体二值图转化为矢量图，分别加载到8月27日的影像上：　　图8红色为8月27日的水域，蓝色为9月9日的水域　　从上图可以明显看到该流域淹没的范围，主要集中在右图红色矩形框内。　　统计水体变化面积：　　图9两幅水体二值图的差值图的DN统计值　　DN值为2的像元数为个，每个像元代表地面大小为*，即9月9日水体增加面积为。然而还有DN值为2的水体，即相对于8月27日，9月9日减少的水体，该面积为。该部分水体主要集中在图10的标示区域，而该区域本身就是湖泊，所以不能算入洪水淹没面积。　　目录　　实习目的.......................................................................................................................2　　实习内容.......................................................................................................................2　　一、操作环境.........................................................................................................2　　二、常用的水提取方法.........................................................................................2　　1、非监督分类法............................................................................................2　　2、监督分类法................................................................................................2　　3、单波段阈值法............................................................................................2　　4、谱间关系法................................................................................................2　　5、归一化水体指数法................................................................................