遥感图像处理实习报告

1、*实 习 报 告 书 题 目： 基于遥感的植被信息提取与分析 学 院： 测绘工程学院 专 业： 班 级： 姓 名： 学 号： 2012年6月6日实习报告评阅书学生姓名： 学 号： * 班 级： * 实践教学环节名称： 遥感图像处理实习 教学时间： 2012年6月6日-2012年6月10日 综合成绩： 指导教师评语： 实习报告成绩： 实习日记成绩： 实习表现成绩： 指导教师(签字)： 2012年 6月 日目 录1 引言11.2 实习名称11.3 实习时间11.4 实习地点11.5 实习目的和任务12 实习内容22.1 地物反射光谱的测定22.1.1 地物反射光谱测定原理22.1.2 地物反射光谱

2、的测定32.2 遥感数据的获取和预处理42.2.1 遥感数据的获取42.2.2 遥感数据的几何精校正42.2.3 遥感数据的大气校正52.3 基于遥感的植被信息提取92.3.1 算法介绍92.3.2 基于遥感的植被信息提取92.4 基于ArcGIS的植被信息处理102.4.1植被信息的面积统计与分析102.4.2 植被信息的专题图制作113 实习体会124 实习日记121 引言1.1 实习准备1.1.1动员大会2012年6月6日9点30分在*，测绘工程学院*老师召开了测绘091班遥感实习动员大会。会上 老师介绍了实习的主要内容、要求以及进程安排。此外，会上还进行了人员分组。1.1.2分组情况本

3、次遥感实习是分组进行的。 1.1.3数据准备本次遥感实习所用到的数据有遥感数据：LS5_TM_20091003_000000_000000_120036_GEOTIFF_L4、原始数据；非遥感数据有连云港矢量边界数据；参考资料，其中有ENVI软件操作的教学指导视频，遥感制图指导等；另外还有光谱仪相关资料，其中有FieldSpec HandHeld 波谱仪的相关介绍，其软件安装步骤和其实验操作步骤。1.1.4个人准备在大三第一学期中系统地学习了遥感数据处理这门课程，在大三的第二个学期中又选修了遥感图像处理。通过学习对遥感数据处理有了一定的理论知识的累积，对于ENVI和ERDAS两款软件也有一定的

4、了解。在本次遥感实习中，运用的软件是ENVI4.0来进行遥感图像的处理。在实习开始前，复习ENVI软件，熟悉相关操作等，并在实习动员时认真领会，做好相关记录。1.2 实习名称 本次实习的名称为基于遥感的植被信息提取与分析。1.3 实习时间本次实习时间为2012年6月6日2012年6月10日。1.4 实习地点 本次实习地点为测绘楼地理信息系统实验室305室。1.5 实习目的和任务1.5.1实习目的本次实习的目的如下：1.熟练操作遥感图像处理软件ENVI 2.掌握遥感图像处理中的关键技术3.学会使用所学知识，掌握技术应用流程4.掌握遥感图像的处理与信息提取5.得到结果的成图与分析1.5.2实习任务

5、本次实习的任务如下：1.地物光谱反射率的野外或室内测定2.遥感图像的几何精校正3.遥感图像的大气校正4.基于遥感影像的植被信息提取5.基于ArcGIS的植被信息的空间分布范围制图，并统计面积6.撰写实习报告、实习日记并提交实习成果2 实习内容2.1 地物反射光谱的测定2.1.1 地物反射光谱测定原理本次实习使用的光谱测定仪器是FieldSpec HandHeld 便携式近红外光谱仪，该仪器适用于350-1050nm波段。该光谱仪配置构成有：1. 照明：a.太阳：最佳照明，70%以上的工作b.室内光源：少量工作2. 样品：a.被测物b.参考板（白板、灰板等）3. 光谱仪4. 控制计算机和软件该光

6、谱仪的技术参数有：1. 光谱范围350-1050nm2. 最快采集速度：17ms3. 光谱分辨率：3.5nm700nm4. 采样间隔：1.5nm5. 波长精度：+/-1nm6. 波长重复性：+/-0.02nm7. 重复性: 优于0.3%8. 杂散光: 优于0.02%9. 灵敏度调整：全自动10. 重量：1.2公斤 （含电池）   地物光谱反射率的测定原理主要是利用电磁辐射和各地物光谱特征进行测定。所有的地面非成像光谱仪的基本工作原理基本是一致的。主要是由光谱仪通过光导线探头摄取目标光线，经由A/D（模/数）转换成数字信号，进入计算机。整个测量过程由计算机通过操作员控制。便携

7、式计算机控制光谱仪并实时将光谱测量结果显示于计算机屏幕上，有的光谱仪带有一些简单的光谱处理软件，如光谱曲线平滑处理、微分处理等。测得的光谱数据可储存在计算机内。为了测定目标光谱，需要测定三类光谱辐射值：第一类称暗光谱，即没有光线进入光谱仪时由仪器记录的光谱；第二类为参考光谱或称标准板白光，实际上是从较完美漫辐射体-标准板上测得的光谱；第三类为样本光谱或目标光谱，是从感兴趣的目标物上测得的光谱，这也是我们最终需要的光谱。为了避免光饱和或光量不足，依照测量时的光照条件和环境温度需要调整光谱仪的测定时间。最后，感兴趣目标的反射光谱是在相同的光照条件下通过参考光辐射值除目标光辐射值得到。因此目标反射光

8、谱是个相对于参考光谱辐射的比值（光谱反射率）。2.1.2 地物反射光谱的测定手持便携式地物光谱仪软件安装步骤：1把笔记本电脑的语言环境设置为英文（美国），因为，RS3软件必须在英文的语言环境下安装运行，在控制面板/区域与语言里面设置，设置完以后电脑会提示重启。2把RS232toUSB电缆线的USB插入电脑提示找到新硬件，安装驱动UT810。3安装完驱动以后右键单击我的电脑/属性/硬件/设备管理器查看RS232toUSB的com端口，（以com4为例）。一般情况下RS232toUSB电缆线与笔记本的USB口不要改动，如果改动的话需要重新安装驱动和查看RS232toUSB的com端口。4安装RS3

9、软件（放在操作手册里面的闪存里面），先装DISK1,再装DISK2,在安装DISK2的时候注意有一个接口选项，默认是EPP端口，把它改为com口，选择刚才在设备管理器查看RS232toUSB的com端口（以com4为例）安装完以后电脑会提示重启。5安装完RS3软件以后，再安装ViewSpec数据后处理软件，在DISK2文件夹里面。安装完以后找到ViewSpec安装路径C：/Program Files/ViewSpecPro找到快捷图标，发送到桌面快捷方式。6安装完软件以后关闭笔记本电脑。光谱仪和电脑的开关顺序：开机先打开光谱仪电源，再打开笔记本电脑。关机先关闭笔记本再关闭计算机。实习操作步骤：

10、1光谱仪、计算机充电。2连接电池、网线、探头电源、光纤，准备好白板。3打开光谱仪电源，然后打开计算机电源，并启动RS3软件。4在软件上调整光谱平均、暗电流平均和白板采集平均次数。5在软件中选择或填写需要存储数据的路径、名称和其他内容。6开始测量：（1）打开探头电源，探头放在白板上面，点击OPT优化；（2）探头仍然对准白板，点击WR采集参比光谱。此时，软件自动进入反射率测量状态。（3）探头移向被测目标的测量位置，按空格键存储采集到的目标反射光谱。7先关闭计算机再关闭仪器。8分析实测结果：（1）准确绘出地物光谱反射率曲线；（2）根据地物光谱反射率曲线，比较地物光谱曲线特征；（3）分析实习过程中可能

11、引起误差的因素。2.2 遥感数据的获取和预处理2.2.1 遥感数据的获取利用FieldSpec HandHeld 便携式近红外光谱仪采集树叶的光谱反射率，再将得到的数据传输到电脑上，利用手持便携式地物光谱仪软件，处理数据，得出植物光谱反射率的曲线，并且得到在曲线上的反射峰和吸收谷的数据。结果如图1.图1 植被光谱反射率曲线2.2.2 遥感数据的几何精校正 1.打开实验数据：原始数据中的header.dat和L5-TM-120L4.TIF，打开路径如图2.图2 打开原始数据路径 2.打开选点对话框，路径如下map->registration->select

12、 GCPs:image to image。 3.在两个图像中选择header为基准图像，L5-TM-120L4为待校正图像。在两幅图中选择非常明显的黑点或者白点的相对应点。在选点前把zoom窗口放大到12x，这样保证所选的的点在一个像元内。在GCP选点窗口中点击add point添加点，先自己添加4个点，然后利用predict进行精确点。在点击predict后，还要在zoom中进行精确选择，保证选点正确。4.在陆地上选择大概30个点左右后，查看image to image GCP list，将其中RMS值大于1或者偏大的值删除，并重新选择，将GCP选点窗口中的RMS

13、 error值小于0.5，然后再利用predict在海域中，选择10个左右的控制点。这次实验中，我的选点RMS error值结果和image to image GCP list结果如图3和图4. 图3 地面控制点选取 图4 地面控制点情况选点的结果如图5和图6，其中图5是待校准图像，图6是基准图像。 图5 待校正图像控制点选取情况 图6 基准图像控制点选取情况2.2.3 遥感数据的大气校正 1.对经过几何精校正的图像进行辐射定标。操作路径为：Basic tools->preprocessing->calibration utilitis->landset TM。 在TM ca

14、libration parameters对话框中进行必要的修改：landset satellite选择5号，采集数据时间为2010年4月29日，太阳高度角为61°，数据含有7个波段，选择radiance，最后选取自己的保存路径，具体情况如图7.图7 辐射定标 2.经过辐射定标之后，需要将定标成果中数据的头文件中的中心波长进行修改。修改路径为：file->edit envi header，在弹出的对话框中选择辐射定标成果，得到图8对话框，在该对话框中选择edit attributes->FWHM，如图9，在FWHM中，将从1波段到7波段的对应的值分别改成0.035,0.04

15、,0.03，0.07,0.1,1.05,0.135，最后点击OK即可。 图8 ENVI头文件修改 图9 修改FWHM值 3.头文件修改之后，进行数据单位的转换。利用basic tools->band math计算得出，计算公式为b1/10，如图10，点击OK，弹出波段加载窗口，如图11，点击map variable to input file，选择辐射定标结果，再选择保存路径，点击OK，就可以得到单位换算之后的图像。 图10 在band math中添加计算公式 图11 实现波段单位转换 4.在改正过单位之后，需要把数据转换一下，操作路径：basic->convert data(BS

16、Q,BIL,BIP),得到格式转换对话框，如图12，选择BIL和输出路径，点击OK即可。图12 格式转换5.选择FLAASH方法进行图像大气校正，操作路径：spectral->FLAASH，即可得到大气校正对话框，如图13，其中参数要进行一下设置：输入格式转换过后的文件，选择输出路径，选择TM5，纬度为34°3643.69，经度为119°3615.40，ground elevation为0.1km，时间为2010年4月29日2:26:20，季节为MLS，区域为urban，可见距离为35km。图13 FLAASH大气校正 6.在图13中大气校正的对话框中，选中multi

17、spectral settings，弹出更多设置，如图14，选择kaufman-tanre aerosol retrieval标签，defaults中选择第一条标签，于是upper channel就对应了第7波段，lower channel对应了第3波段，点击OK，回到大气校正窗口，点击apply即可。图14 大气校正更多设置 完成大气校正后，弹出大气校正结果，如图15.图15 大气校正结果7.打开大气校正后的结果图，在图像窗口右击，选择quick stats，得到统计结果，如图16，图中Max一栏中都小于10000，并且所有的数据都为正，所以此次大气校正较为成功。图16 大气校正quick

18、stats结果 在图像显示窗口菜单中选择路径：tools->profiles->z profile，得到曲线图，图17。图17 大气校正后曲线图2.3 基于遥感的植被信息提取2.3.1 算法介绍在基于遥感的植被信息提取部分，老师给我们提供了三种算法，我是属于第4小组的，于是我们第4小组的提取植被信息算法为：NDVI=(TM4-TM3)/(TM4+TM3)。2.3.2 基于遥感的植被信息提取1.由于数据的单位有差异，因此在进行算法之前，首先要进行数据的单位换算。利用band math编辑公式，就可以将数据换算。用到的公式为float(b1)/10000.0，在添加波段时，由于我们第4

19、小组只需要用到波段3和波段4，所以，只添加了波段3和波段4，这样的话，计算会相对快速些，也减少了运算的数据量。2.在换算过单位之后，就要进行植被信息的提取了。在这里，老师给我们提供了数据，0.3以上的归为植被，0.3以下的归为非植被。操作路径为图像显示窗口菜单overlay->density slice，如图18，在图中，将最小值到0.3的定义为非植被，颜色选择为白色，将0.3到最大值的定义为植被，颜色选择为绿色，点击apply即可。图18 提取植被信息 3.在提取完植被信息之后，得到分类图像，将文件先存为矢量数据，然后再存储为shp格式的文件，这样在接下来的专题地图制作中，可以使用sh

20、p格式文件。其过程如图19和图20。 图19 转换成矢量数据 图20 转换成SHP格式文件2.4 基于ArcGIS的植被信息处理2.4.1植被信息的面积统计与分析 1.打开arcmap软件，点击“+”添加植被信息提取的SHP格式文件，打开data frame properties，在general标签中将map和display的单位都改为meters，如图21.图21 改变单位2.选中植被提取信息文件，右击选择查看其数据信息，在最后一栏area中可以看到其面积，将其全部求和，就可以得到植被的总面积，如图22. 图22 植被总面积2.4.2 植被信息的专题图制作 1.在arctoolbox中选择data management tools->projection and transform->define projection，为植被信息添加一个坐标系统。在弹出的对话框中选择输入文件为植被提取信息，坐标系统为原始数据中的transverse mecator，点击OK即可。 2.添加图框。选择legend wizard，一直点击下一步，最后完成即可。3.添加标题。在insert中选择title，在适当的位置添加标题，标题内容为基于遥感的研究区植被信息空间分布。4.添加比例尺。在insert中选择scale，选择适合的格式，添加到相应位置。5.添加图例。