遥感地学分析第六讲-植被变化动态监测课件

遥感地学分析第六讲植被变化动态监测遥感地学分析第六讲植被变化动态监测1第六讲植被变化动态监测专题概述处理流程介绍第六讲植被变化动态监测专题概述2一、专题概述土默特右旗，简称土右旗，位于内蒙古自治区中南部，黄河北岸，呼和浩特市、包头市之间。总面积2368平方千米。总人口35万人（2004年。土地肥沃，水源丰富，气候适宜。经济以农业为主，主要作物有小麦、莜麦、糜黍，葵花籽、甜菜、胡麻以及青山党参、正北黄芪等名贵药材。一、专题概述土默特右旗，简称土右旗，位于内蒙古自治区中南部，3一、专题概述本专题详细介绍了利用TM影像计算植被覆盖度的完整流程，基于ENVI4.6版本操作，专题涉及数据读取、辐射定标、大气校正、植被覆盖的计算、植被覆盖变化监测等内容。一、专题概述本专题详细介绍了利用TM影像计算植被覆盖度的完整4二、处理流程介绍数据读取和定标暗象元法大气校正波段合并研究区裁剪多时相影像配准归一化植被指数（NDVI）的计算植被覆盖度的计算植被动态变化监测成果后期处理与应用二、处理流程介绍数据读取和定标5（一）数据读取和定标读取两个时相（2002.09.23和2006.09.26）的1、2、3、4波段文件并对其进行辐射定标。（一）数据读取和定标读取两个时相（2002.09.23和206（二）暗象元法大气校正采用暗象元法对两个时相的四个波段的辐射定标文件进行大气校正。（二）暗象元法大气校正采用暗象元法对两个时相的四个波段的辐射7（三）波段合并将两个时相大气校正后的四个波段文件进行合并，得到两个多波段文件。（三）波段合并将两个时相大气校正后的四个波段文件进行合并，得8（四）研究区裁剪将土默特右旗的行政边界文件设为感兴趣区域，再根据感兴趣区域对合并后的多波段文件进行裁剪。（四）研究区裁剪将土默特右旗的行政边界文件设为感兴趣区域，再9（五）多时相影像配准将裁剪后的两个时相的多波段文件进行配准，并修改其尺寸，使两个文件的大小一致。（五）多时相影像配准将裁剪后的两个时相的多波段文件进行配准，10（六）NDVI的计算计算裁剪并配准后的两幅影像的NDVI。（六）NDVI的计算计算裁剪并配准后的两幅影像的NDVI。11（七）植被覆盖度的计算计算植被覆盖度FC采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下：

FC=[(NDVI-NDVISoil)／(NDVIVeg-NDVISoil)]（1）其中，NDVI为归一化植被指数，NDVISoil为完全是裸土或无植被覆盖区域的NDVI值，NDVIVeg则代表完全被植被所覆盖的像元的NDVI值，即纯植被像元的NDVI值。取经验值NDVIVeg=0.70和NDVISoil=0.00，且有，当某个像元的NDVI大于0.70时，FC取值为1；当NDVI小于0.00，FC取值为0。（七）植被覆盖度的计算计算植被覆盖度FC采用的是混合像元分12（七）植被覆盖度的计算利用basictools/BandMath，在公式输入栏中输入：(b1gt0.7)*1+(b1lt0.)*0+(b1ge0andb1le0.7)*((b1-0.0)/(0.7-0.0))b1:选择NDVI图像。（七）植被覆盖度的计算利用basictools/Band13（八）植被覆盖动态监测（1）2006年9月26日植被覆盖区提取：利用basictools/BandMath，在公式输入栏中输入：(b1le0.3)*0+(b1gt0.3)*1b1:选择2006年植被覆盖度图像，0.3为经验值。（八）植被覆盖动态监测（1）2006年9月26日植被覆盖区提14（八）植被覆盖动态监测（2）2002年9月23日植被覆盖区提取：利用basictools/BandMath，在公式输入栏中输入：(b1le0.3)*0+(b1gt0.3)*1b1:选择2002年植被覆盖度图像，0.3为经验值。（八）植被覆盖动态监测（2）2002年9月23日植被覆盖区提15植被变化检测利用basictools/BandMath，在公式输入栏中输入：b1-b2b1:选择2006年9月的植被覆盖区图像。b2:选择2002年9月的植被覆盖区图像;植被变化检测利用basictools/BandMath16（九）成果后期处理与应用植被变化区域图的背景值处理背景值与影像DN值相同，易混淆，需要将二者区分开来。（九）成果后期处理与应用植被变化区域图的背景值处理背景值与影17植被变化区域图的背景值处理此时，需要对背景值进行掩膜处理，改变其DN值。方法：（1）将土默特右旗的行政边界设置为需要进行掩膜处理的影像（植被覆盖图）的感兴趣区域。（2）使用ROI生成掩膜（3）应用掩膜植被变化区域图的背景值处理此时，需要对背景值进行掩膜处理，改18使用ROI生成掩膜在ENVI的菜单下选择BasicTools/Masking/BuildMask。使用ROI生成掩膜在ENVI的菜单下选择BasicTool19使用ROI生成掩膜选择需要进行掩膜的影像使用ROI生成掩膜选择需要进行掩膜的影像20使用ROI生成掩膜在掩膜定义窗口中，在Options下选择ImportROIs，然后选择前面生成的面状ROI。使用ROI生成掩膜在掩膜定义窗口中，在Options下选择I21使用ROI生成掩膜设置ROI的类型，在Options菜单下通过设置SelectedAreas‘off’或是’On’来设定ROI的类型。点击Apply，输入掩膜文件保存的位置，生成掩膜文件。使用ROI生成掩膜设置ROI的类型，在Options菜单下通22应用掩膜在BasicTools下选择Masking

àApplyMask，首先选择需要进行掩膜的文件，然后在选择文件的对话框中选择SelectMaskBand，在弹出的对话框中选择前面生成的掩膜文件。应用掩膜在BasicTools下选择MaskingàA23遥感地学分析第六讲_植被变化动态监测课件24遥感地学分析第六讲植被变化动态监测遥感地学分析第六讲植被变化动态监测25第六讲植被变化动态监测专题概述处理流程介绍第六讲植被变化动态监测专题概述26一、专题概述土默特右旗，简称土右旗，位于内蒙古自治区中南部，黄河北岸，呼和浩特市、包头市之间。总面积2368平方千米。总人口35万人（2004年。土地肥沃，水源丰富，气候适宜。经济以农业为主，主要作物有小麦、莜麦、糜黍，葵花籽、甜菜、胡麻以及青山党参、正北黄芪等名贵药材。一、专题概述土默特右旗，简称土右旗，位于内蒙古自治区中南部，27一、专题概述本专题详细介绍了利用TM影像计算植被覆盖度的完整流程，基于ENVI4.6版本操作，专题涉及数据读取、辐射定标、大气校正、植被覆盖的计算、植被覆盖变化监测等内容。一、专题概述本专题详细介绍了利用TM影像计算植被覆盖度的完整28二、处理流程介绍数据读取和定标暗象元法大气校正波段合并研究区裁剪多时相影像配准归一化植被指数（NDVI）的计算植被覆盖度的计算植被动态变化监测成果后期处理与应用二、处理流程介绍数据读取和定标29（一）数据读取和定标读取两个时相（2002.09.23和2006.09.26）的1、2、3、4波段文件并对其进行辐射定标。（一）数据读取和定标读取两个时相（2002.09.23和2030（二）暗象元法大气校正采用暗象元法对两个时相的四个波段的辐射定标文件进行大气校正。（二）暗象元法大气校正采用暗象元法对两个时相的四个波段的辐射31（三）波段合并将两个时相大气校正后的四个波段文件进行合并，得到两个多波段文件。（三）波段合并将两个时相大气校正后的四个波段文件进行合并，得32（四）研究区裁剪将土默特右旗的行政边界文件设为感兴趣区域，再根据感兴趣区域对合并后的多波段文件进行裁剪。（四）研究区裁剪将土默特右旗的行政边界文件设为感兴趣区域，再33（五）多时相影像配准将裁剪后的两个时相的多波段文件进行配准，并修改其尺寸，使两个文件的大小一致。（五）多时相影像配准将裁剪后的两个时相的多波段文件进行配准，34（六）NDVI的计算计算裁剪并配准后的两幅影像的NDVI。（六）NDVI的计算计算裁剪并配准后的两幅影像的NDVI。35（七）植被覆盖度的计算计算植被覆盖度FC采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下：

FC=[(NDVI-NDVISoil)／(NDVIVeg-NDVISoil)]（1）其中，NDVI为归一化植被指数，NDVISoil为完全是裸土或无植被覆盖区域的NDVI值，NDVIVeg则代表完全被植被所覆盖的像元的NDVI值，即纯植被像元的NDVI值。取经验值NDVIVeg=0.70和NDVISoil=0.00，且有，当某个像元的NDVI大于0.70时，FC取值为1；当NDVI小于0.00，FC取值为0。（七）植被覆盖度的计算计算植被覆盖度FC采用的是混合像元分36（七）植被覆盖度的计算利用basictools/BandMath，在公式输入栏中输入：(b1gt0.7)*1+(b1lt0.)*0+(b1ge0andb1le0.7)*((b1-0.0)/(0.7-0.0))b1:选择NDVI图像。（七）植被覆盖度的计算利用basictools/Band37（八）植被覆盖动态监测（1）2006年9月26日植被覆盖区提取：利用basictools/BandMath，在公式输入栏中输入：(b1le0.3)*0+(b1gt0.3)*1b1:选择2006年植被覆盖度图像，0.3为经验值。（八）植被覆盖动态监测（1）2006年9月26日植被覆盖区提38（八）植被覆盖动态监测（2）2002年9月23日植被覆盖区提取：利用basictools/BandMath，在公式输入栏中输入：(b1le0.3)*0+(b1gt0.3)*1b1:选择2002年植被覆盖度图像，0.3为经验值。（八）植被覆盖动态监测（2）2002年9月23日植被覆盖区提39植被变化检测利用basictools/BandMath，在公式输入栏中输入：b1-b2b1:选择2006年9月的植被覆盖区图像。b2:选择2002年9月的植被覆盖区图像;植被变化检测利用basictools/BandMath40（九）成果后期处理与应用植被变化区域图的背景值处理背景值与影像DN值相同，易混淆，需要将二者区分开来。（九）成果后期处理与应用植被变化区域图的背景值处理背景值与影41植被变化区域图的背景值处理此时，需要对背景值进行掩膜处理，改变其DN值。方法：（1）将土默特右旗的行政边界设置为需要进行掩膜处理的影像（植被覆盖图）的感兴趣区域。（2）使用ROI生成掩膜（3）应用掩膜植被变化区域图的背景值处理此时，需要对背景值进行掩膜处理，改42使用ROI生成掩膜在ENVI的菜单下选择BasicTools/Masking/BuildMask。使用ROI生成掩膜在ENVI的菜单下选择BasicTool43使用ROI生成掩膜选择需要进行掩膜的影像使用ROI生成掩膜选择需要进行掩膜的影像44使用ROI生成掩膜在掩膜定义窗口中，在Options下选择ImportROIs，然后选择前面生成的面状ROI。使用ROI生成掩膜在掩膜定义窗口中，在Options下选择I45使用ROI生成掩膜设置ROI的类型，在Options菜单下通