遥感地学分析报告试验报

1、实验一 植被覆盖度反演一、实验目的 植被覆盖度是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区总面积的百分比。通常林冠 称郁闭度，灌草等植被称覆盖度。它是衡量地表植被覆盖的一个最重要的指标，被覆盖度及其变化是区域 生态系统环境变化的重要指示，对水文、生态、全球变化等都具有重要意义。目前已有许多利用遥感技术 测量植被覆盖度的方法，其中应用最广泛的方法是利用植被指数近似估算植被覆盖度，常用的植被指数为 NDVI ，本次实验完成植被覆盖度反演。二、实验数据 实验选取两景覆盖市的 Landsat8 OLI 影像、土地覆盖类型图以及行政边界矢量数据为数据源。其中， 土地覆盖类型图是作为掩膜文件使用

2、，其目的是为了便于植被覆盖度的估算；行政边界矢量数据是裁剪出 市行政区的围。 Landsat8 OLI影像是从地理空间数据云上下载得到的，其成像时间为2013年10月份。与Landsat7的ETM +成像仪相比， OLI 成像仪获取的遥感图像辐射分辨率达到 12比特，图像的几何精度和数据 的信噪比也更高。 OLI 成像仪包括 9个短波谱段(波段 1波段 9)，幅宽 185km，其中全色波段地面分辨率 为 15m，其他谱段地面分辨率为 30m 。三、实验方法 本文反演植被覆盖度所采用的是像元二分模型方法，像元二分模型是一种简单实用的遥感估算模型， 它假设一个像元的地表由有植被覆盖部分与无植被覆盖

3、部分组成， 而遥感传感器观测到的光谱信息 (S)也由 这2个组分因子线性加权合成，各因子的权重是各自的面积在像元中所占的比率，如其中植被覆盖度可以 看作是植被的权重。 因此，像元二分模型的原理如下： VFC = (S - Ssoil)/ ( Sveg - Ssoil)S为遥感信息， 其中 Ssoil 为纯土壤像元的信息 , Sveg 为纯植被像元的信息。改进的像元二分法 遥感信息选择为 NDVIVFC = (NDVI - NDVIsoil)/ ( NDVIveg - NDVIsoil) 两个参数的求解公式NDVIsoil=(VFCmax*NDVImin- VFCmin*NDVImax)/( V

4、FCmax- VFCmin) NDVIveg=(1-VFCmin)*NDVImax- (1-VFCmax)*NDVImin)/( VFCmax- VFCmin) 当区域可以近似取 VFCmax=100% ， VFCmin=0%VFC = (NDVI - NDVImin)/ ( NDVImax - NDVImin) 当区域不可以近似取 VFCmax=100% ， VFCmin=0% ，当有实测数据的情况下，取实测数据中的植被覆 盖度的最大值和最小值；当没有实测数据的情况下，植被覆盖度的最大值和最小值根据经验估算。其中 , NDVIsoil 为裸土或无植被覆盖区域的 NDVI 值 , 即无植被像元

5、的 NDVI 值；而 NDVIveg 则代表 完全被植被所覆盖的像元的 NDVI 值 , 即纯植被像元的 NDVI 值。四、实验处理步骤1、实验处理流程如下图所示文档2、数据预处理 本文使用的 Landsat8 OLI为 L1T级别数据，不需做几何校正处理。而市需要两景 因此首先要进行图像镶嵌和裁剪，然后进行大气校正等预处理过程。Landsat OLI 数据覆盖，（ 1） 辐射定标 辐射定标是将传感器记录的电压或数字值转换成绝对辐射亮度的过程。实验报告中写出辐射定标的作用。该处理过程在 Envi5.1 中实现，具体操作：在 ENVIToolbox 中，选择 Toolbox/Radiometri

6、c Correction/ Radiometric Calibration ，选择 *_MultiSpectral 多光谱组（ 7个波段），打开辐射定标工具，对两景影像分别做 辐射定标。文档（2）影像镶嵌 因本文所使用的影像数据源是两景 Landsat OLI 影像，因此需进行影像镶嵌，镶嵌的目的是将不同的 影像文件无缝地拼接成一幅完整的包含研究区域的影像。该处理过程在 Envi5.1 中实现，具体操作：在 Toolbox 中，选择 /Mosaicking/Seamless Mosaic ，打开无缝镶嵌工具，然后进行相关参数设置。MB匕f$l F: novsql 箱T t r-3 Ekssi

7、MbT nf - !-= EV tlT文档（3）影像裁剪 因本文所使用的影像数据包含了市行政区划以外的部分地区，因此需进行影像裁剪，以将研究区裁剪 出来，并且减小了数据量，加快了数据处理速度，本文使用行政边界矢量裁剪图像。过程在 Envi5.1 中的 具体操作如下： 在 Toolbox 中，选择 /Regions of Interest/Subset Data from ROIs ，打开裁剪工具：文档（4）Flaash 大气校正ENVI中的 FLAASH模型是基于 MODTRAN4辐+射传输模型，通过参数查找表来进行大气校 正的商业化软件。实验报告中需写出为什么用进行大气校正。在 Toolbo

8、x 中打开 FLAASH 工具 /Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction报告中需要对大气校正前后同一地物的光谱曲线进行对比。文档3、植被覆盖度估算（ 1）计算 NDVI本文选取 NDVI 值为参数，采用像元二分模型对植被覆盖度进行反演，根据植被覆盖度的计算公式可 知，要求取植被覆盖度，首先需要计算 NDVI。在 Envi5.1 中的具体操作如下：在 Toolbox 中，选择 Spectral/Vegetation/NDVI ， NDVI Calculation Inpu

9、t File 面板中，选择大气校正后的 图像，求算 NDVI，如下：由于大气校正后的结果有部分像元为负值， 主要集中在阴影地区， 这部分区域计算得到的 NDVI 在-1,1 之外，为了便于后面的分析，这里统一将这部分像元进行处理，即 NDVI 值大于 1 的变为 1，小于 -1 的变 成-1。使用 Bandmath 工具，（写出运算公式） ，得到去除异常值文件。文档(b1 lt (-1)*(-1) + (b1 gt (1)*1 + (b1 le(1) and b1 ge(-1)*b1（ 2） 生成掩膜文件 该过程主要是为了计算 NDVI的最大值、最小值，根据土地利用分类图（共 5 类，林地、农

10、业用地、 城市用地、水体与其他）制作各种土地利用类型的掩膜文件，在 Envi5.1 中的具体操作如下：在 Toolbox 中选择 /Raster Management/Masking/Apply Mask ，打开制作掩膜工具：文档得到林地、农业用地、城市用地、水体与其他的掩膜文件。应用掩膜文件文档（3）获取阈值 计算 NDVImax 和NDVImin 值，使用获取的掩膜文件分别对 NDVI 图像文件进行统计，在一定置信度围 获取每个掩膜文件（也就是土地覆盖类型）对应的最大和最小 NDVI 值。在Toolbox中，选择 /Statistics/Compute Statistics ，进行统计，然

11、后在统计结果中，取一定的置信度获取，选择 NDVImin=0.3804 ，NDVImax=0.8667 。 水体没有植被（水藻不属于植被） ，认为这部分最大和最小的 NDVI 值。林地覆盖区域的统计结果（如下图） 同样的方法得到其他地物覆盖类型的 NDVI 阈值，其中 区域的植被覆盖度为 0，如下表：文档土地覆盖类型NDVImin(NDVIsoil )NDVImax(NDVIveg)林地-0.0039220.937255农业用地-0.0039220.701961城市用地-0.0039220.301961水体00其他-0.0039220.3490204)生成参数文件植被覆盖度的计算公式： VFC

12、 =(NDVI - NDVIsoil ) (NDVIveg -NDVIsoil)，该过程是根据上面得到的NDVI阈值分别生成 NDVIsoil和NDVIveg参数文件，也即 NDVImin 与NDVImax 。该过程主要使用 Envi5.1 的 bandmath 工具，并且：NDVIsoil：b1*0.128627+b2*0.090196+b3*0.011765+b4*0+b5*0.003922其中， b1 ：林地掩膜文件， b2 ：农业用地掩膜文件， b3 ：城市用地掩膜文件， b4 ：水体掩膜文 件， b5：其他用地掩膜文件文档NDVIveg：b1*0.992157+b2*0.992157

13、+b3*0.568627+b4*0+b5*0.639216其中， b1:林地掩膜文件， b2：农业用地掩膜文件， b3 ：城市用地掩膜文件， b4：水体掩膜文件， b5：其他用地掩膜文件 最终，生成的参数文件。4）植被覆盖度估算利用上一步得到的 NDVIsoil 和NDVIveg 参数文件带入公式： VFC =（NDVI - NDVIsoil ） （NDVIveg-NDVIsoil） ，该过程也是利用 Envi5.1 中的 Bandmath工具来实现（写出运算公式）。分析下结果，会发现有一些异常值，即值在0,1 之外，这些异常值是在 NDVI置信度之外的那部分像元 产生的（也包括 NDVI异常

14、像元）。这些像元数量不多，大约占 3.7%左右。还有背景和水体区域的植被覆盖文档度的值为 -NaN，即无效值， 因为分母为 0造成的。 第一种异常值可以将小于 0 的值变成 0，大于 1的值变成 1， 用 bandmath 工具即可， Bandmath 表达式为： 0.0b1Compute ROI Separability, 计算样本的可分离性。其中 Jeffries-Matusita, Transformed Divergence 参数表示，这两个参数的值在 02.0 之间，大于 1.9 说 明样本之间可分离性好，属于合格样本；小于 1.8，需要重新选择样本；小于 1，考虑将两类样本 合成一

15、类样本。文档08 年的09 年的文档2、土地利用分类有以下两个分类流程方法，选择其中一种方法即可。1）采用监督分类方法包括平行六面体、最小距离、马氏距离、最大似然，基于神经网络模型，支持向量机、模糊分类等， Toolbox 选择 Classification Supervised 多种分类模型，实验挑选 种分类模型进行土地利用分类（报告中标出采用何种分类模型），并对分类结果的处理，使用 Classification-Post Classification 工具中方法。最大似然法08 年：文档2) Toolbox 中，打开 /Classification/Classification Workf

16、low ，采用图像分类流程化工具。3、分类精度评定使用 Toolbox 中 Classification-Post Classification-Confusion Matrix ，得到分类精度评价混淆 矩阵和 Kappa 系数，对分类结果进行评价分析。文档4、土地利用变化监测使用 Thematic Change Workflow 分类后比较法工具， 分析两期分类数据的变化情况， 并在流 程中Smoothing 和 Aggregation 中设置合适的值去除噪声和合并小斑块。 结果以图像输出， 并统计文档变化的结果，分析各种地类变化的情况。文档四、实验心得 应用遥感技术进行土地利用变化监测研究

17、是土地调查研究的重要组成部分，本子实验要注意 店颜色的选取和样点选取数量要多，全图选取这样劫过才不会出现大误差。文档实验三：基于规则的道路信息提取一、实验目的道路作为基础地理信息，为保持数据的现势性，需要对该类数据实时或准实时更新，在现实 中，基于高分辨率的数据道路数据的提取是更新道路信息数据的有效方法。本次实验使用基于规 则的面向对象信息提取方法，完成自动提取道路信息的过程。二、实验容实验容包括高分辨率数据的处理， 高分辨率信息提取流程， ENVI FX 模块面向对象信息提取， 矢量后处理的方法。实验数据采用 Quickbird 影像数据( QB-CJ)，应用 ENVI FX 扩展模块中的

18、Feature Extraction 工具，提取影像中道路信息。三、实验步骤1、启动 Rule Based FX 工具在Toolbox 中，选择 /Feature Extraction/Rule Based Feature Extraction Workflow。关注 Custom Bands 面板，有两个自定义波段，包括归一化植被指数或者波段比值、HSI 颜色空间，这些辅助波段可以提高图像分割的精度，如植被信息的提取等自定义的属性。选择Spectral，Band 下面选择Normalized Difference 。在第一步自定义波段中选择的波段是红色和近红外波段， 所以在此计算 的是NDV

19、I ，参数设置为 band1：band2， band2：band4。2、设置影像分割阈值( 20,30,40)、合并阈值 (70,85)文档3、制定提取规则根据规则进行特征提取，在规则分类界面。每一个分类有若干个规则（Rule）组成，每一个规则有若干个属性表达式来描述。道路提取的规则如下（根据数据可以动态调整阈值的大小，能 较好的反映道路特征）剔除植被： NDVI 的 Spectral mean设小于 0.01 剔除高亮度的建筑物和空地，第一波段的 Spectral mean阈值小于 310 剔除面积较小的建筑物和空地，参数 Spatial-Area，阈值大于 2500 剔除延伸性小于 3的地

20、物，参数 Spatial-Elongation ，阈值大于 34、结果输出与处理对提取的结果叠合与处理，得到道路信息特征，如果效果不好，调整制定的提取规则文档四、实验心得 道路信息作为一种重要的基础地理信息，可以作为提取其他地物目标的线索和参考系，具有很强的现 实意义。从遥感影像自动提取人工地物,特别是线状地物 （主要是道路 ）,不仅是摄影测量与遥感领域的难题也是计算机视觉与图像理解研究的重点之一。本次试验中，主要是 使用基于规则的面向对象信息提取方法，完成自动提取道路信息的过 程。实验四 太湖水体叶绿素浓度反演一、实验目的根据环境小卫星 CCD数据，反演太湖的叶绿素 a 浓度，要求掌握环境小

21、卫星的数据读取、几何校正、 大气校正、反演模型的建立、遥感反演过程、反演结果验证等操作技能及原理，熟悉遥感水质反演过程。二、实验要求根据环境小卫星 CCD数据特点及太湖水质反演模型要求，先对环境小卫星CCD数据进行数据预处理，大气校正、太湖区裁剪， 利用波段比值法对实测的叶绿素 a 浓度数据建立反演模型， 将模型应用于太湖水 面区域影像，反演出整个太湖区的叶绿素 a 浓度 ,并验证模型精度。三、实验步骤（一）数据预处理1、安装环境小卫星数据读取和定标补丁，复制到HomeProgram Files ENVI51Extensions 目录下。本次实验中已完成小卫星数据转化和定标，得到 HJ1B-C

22、CD1-20091006-Cal-sub 数据，用于流域叶绿素反演。 2、几何校正（图像配准），以 TM 作为基准影像对环境小卫星图像进行图像配准。（ 1）打开基准影像 TM_baseimage.img 。（2）选择 Tool中Registration -Automatic Registration ： Image to Image，选择基准影像第 4波段作为匹配波 段（ TM_baseimage.img ）。（3）选择被配准影像 HJ1B-CCD1-20091006-Cal-sub.img 第 4波段作为被配准波段。文档（ 4）选择控制点，并调整控制点，直到总的RMS Error 小于 1个

23、像素时，完成控制点的选择。点击GroundControl Points Selection 上的 File-Save GCPs to ASCII ，保存控制点。调整控制点后的 RMS 错误5）在 Ground Control Points Selection 面板上，选择 Options-Warp File （As Image Map ），选择校正文件。 在校正参数面板中，投影参数默认。（6）在 X 和 Y 的像元大小输入 30米，按回车，图像输出大小自动更改重采样方法选择 Nearest Neighbor （保持图像辐射值不变），输出为 registratio文档3、大气校正nm），后面四列

24、分别表示 4 个波环境小卫星提供了波谱响应函数，以文本形式提供，第一列表示波长（ 段对应波长的波谱响应值。为了做大气校正，需要制作波谱曲线来描述波谱响应函数，用于大气校正。（1）制作波谱曲线 打开 Display-New Plot Window 面板，在波谱绘制窗口中，选择 ImportASCII ，导入 “ 681_HJ1ACCD1.txt” 文本文件。 在绘制窗口生成了 4 条曲线， 选择 Edit-Data Parameters ，编辑每条线的名称为 b1，b2，b3，b4， 便于区分。将数据 Export-Spectral Library ，在 Output Spectral Libr

25、ary 面板中，有输出曲线相关参数设置， 将波谱曲线保存为波谱库文件 “HJ1B-CCD1 光谱响应 .sli 。”(2) FLAASH 大气校正1) 主菜单 Spectral-FLAASH打开 FLAASH 大气校正模块；首先采用 Convert Data( BSQ、BIL 、BIP ）工具，转换为 BIL 格式的数据。文档2）大气模型选择 Mid-Latitude Summer ，气溶胶模型选择 Rural ，气溶胶反演方法选择 None，能见度给 40km 。3）大气模型选择 Mid-Latitude Summer ，气溶胶模型选择 Rural ，气溶胶反演方法选择 None，能见度给

26、40km 。4）大气校正完成后，检查大气校正的结果，分别加载校正前后的图像，将两幅影像进行地理，移动到植 被区域（植被的波谱曲线比较特殊），在影像上右键，选择Z Profile （ Spectrum）打开光谱曲线窗口，显文档 示两幅图像同一位置的光谱曲线图。4、太湖区裁剪使用 shp 格式太湖围数据，对影像数据进行裁剪(二)叶绿素反演实验中选择波段比值法( BNIR/B RED )建立模型。1、采样点实测数据处理 需要结合实地调查数据，将水面调查点与实测叶绿素浓度在空间上相对应，其中实地调查数据中包括水面 调查点的经纬度、叶绿素含量，实验数据使用 “叶绿素实测数据 ”。2、获取采样遥感影像上的数据1) 采用波段计算器，输入 float(b4)/b3 ，计算得到比值图像。文档2) 打开 Envi classic ，Basic tool-Region of Interest-ROI Tool 工具，选择 ROI_Type-Input Points from ASCII ， 选择文本格式的 “反演点 .txt ”。注意参数选择， x：选择经度； y：选择纬度； These point comprise：Individual Points。设置投影信息地理经纬度。将实地调查的点位置信息加载到图像中。3) 在 ROI Tool 中，选择 File-Output ROIs