基于像元二分模型的植被覆盖度反演

实习报告题目XX市基于像元二分模型的植被覆盖度反演学院测绘科学与技术学院专业班级测绘工程1903班姓名XXX任课教师XX、XXX222一．实习目的222二．实习地点及时间...............................................................三．实习具体内容................................................................. 2 32.1数据读取与定标 3 4 5 6 63.2掩膜文件制作 8 9 9 3.6植被覆盖度验证和分类 4.植被覆盖度时空变化分析 4.1各期植被覆盖度的基本情况 4.2植被覆盖度的动态变化检测 一．实习目的掌握ENVI软件的基本操作掌握遥感卫星影像的处理流程提高学生分析问题解决问题能力，增强实学会地面测量和遥感估算两种方法。地面测量常用于田间尺度，遥感估算常用于区域尺度。目前已经二．实习地点及时间三．实习具体内容2.数据预处理2.1数据读取与定标示，对数据做全景显示和拉伸，在（1）在ENVIT）olbox中选择Tdolbox/Radiom*_注：此步骤可以使用shp文件对图像进行SpatialSubset，以减少数据量。M（2）辐射定标工具会从元数据文件中自动获取相关的参数信息，包括成像日期、定标参数等，ltiSpecal单击FLAASHSetttings按钮，自动计算ScaleFactor：0.1注：默认OLI定标参数得到的辐射亮度值的单位为W/气校正要求的辐射亮度值的单位为：（卩W）/），t（3）选择输出路径和文件名2014_rad.dat，单击OK执行定标处理结果如下:同样的方法对另一景Landsat8多光谱数据进行辐射定标。选择输出路径和文件名2014_rad.dat。2.2图像裁剪AZITrT-iW曲'EEfK<T3rri_.lit(2)在Toolbox中，选择/Regionsoflnterest/SubsetDatafromR在在"MaskpixelsoutsideofROI?”项中，选择裁图5所示:,单击OK裁剪图像，结果图如(;爭■anHjedhil£=j?jca(]o2.3大气校正在Toolbox中，选择/RadiometricCorrection/FLAASHAtmosphericCorrection，打开FLAASH（1）文件输入与输出信息项目2014_rad_tianshui.dat”。在RadianceScaleFactors对话框中选择Usesinglescalefactorforallbands（Singlescalefactor：1，在定标中已经进行了单位转换，这里直大7)多光谱设置(MultispectralSettings)n羽iwm1ffkdtlikdtlrr选选⑵在NDVI^CalculationParameters对过提示：因为目前NDVI计算工具只能识别原始Landsat8OLI数据，我们这里输入的是经过ta由于大气校正后的结果有部分像元为负值，主要集中在阴在［-1,1］之外，为了便于后面的分析，我们这里统一将这部分像元进行处理，即值大于1的变3.2掩膜文件制作我们利用这五种土地覆盖类型分别为NDVI生成5个掩膜文件。下面以林地为例子介绍classimage.dat。在最大、最小值中同时填入(4)选择输出路径输出，完成掩膜文件的生3.3获取阈值件进行统计，在一定置信度范围内获取每个掩膜文件(也就是土地覆盖类型)对应的最大和(1)在Toolbox中，选择/Statis“NDVI_去除异常值.dat”作为统计文件，以及掩膜文件。表土地覆盖类型3.4生成参数文件这一步是根据上面得到的NDVI阈值分别生成NDVNDVIsoil：b1*0.121569+b2*0.147411+b3*(-0.0352NDVIveg：b1*0.678431+b2*0.214603+b3*0.270588+b4*0.165195+b5*0.427602下:L^w二1'曹[■"w-aw』MSIffiMC•二].HwmLJLFllBIt雌HITLitvLH.THI.ISvalniXlvIlMW*3_4Jr»»ni*i-tUM*1]Ijtf:l*jiv4.lvM«h«i缶KfflOiW*厨.^330]ad.i>uXfatr-arliiaiaa.+!rb»Fillwliiaiaa.+!rrCrCanr^li一-EtdirF-M>!EtdirF-M>!hid匚ZbILSari*Wi'iiMliaMT这一步是根据上面得到的NDVI阈值分别生成NDVTNDVIsoil：b1*0.32765+b2*0.163286+b3*（0.184314）+b4*0.053726+b5*0.129412NDVIveg：b1*0.71365+b2*0.275404+b3*0.505882+b4*0.433725+b5*0.615686赳观Tn赳观TnA巫口□■榊Ci■■JE-IITAI3_I*J_ILB>Ek_杠l*ivi建CLKUM,,—i"•_JJataOiuLtiJat□B□BMLnnaMLd.HrLifflH3.5植被覆盖度我们分析下结果，会发现有一些异常值，即值在[0,1]之外，这些异常值是在NDVI置信度之外的那部分像元产生的（也包括NDVI异常像元）。这些像元数量不多，大约占3.左右。还有背景和水体区域的植被覆盖度的值为-NaN,即无效值，因为分母3.6植被覆盖度验证和分类与地面垂直的角度用数码相机拍摄采样点的地面照片，使用截取长、宽各三分之二的中心地带，采用非监督分类法分为将数码相机得到数字图像被、非植被两类，以此来计算出植被覆盖度。为了保证验证时能够正确定位，野外一般选(2)依次添加5个区间，分别为每个区间设置一定的颜色，单击⑶在图层管理器中，在Slices上右键，选择Exportclasssl4.植被覆盖度时空变化分析4.1各期植被覆盖度的基本情况行了统计分析。如表三、表四0101百分比/%低植被覆盖率全植被覆盖率004.2植被覆盖度的动态变化检测为了进一步分析城市植被覆盖等级的时空变表五2014年一一2018年植被覆盖等级变化监测43210envi-完整的遥感图像处理平台envi(theenvironmentforvisualizingimages)是美国exelisvisualinformationsolutions公司的旗舰产品。它是由遥感领域的科学家釆用交互式数据语言idl(interactivedataIanguage)开发的一套功能强大的遥感图像处理软件。它是快速