遥感应用模型遥感反演土地变化检测

遥感应用模型实习报告

学院：

班级：

学号：

姓名：

日期：

指导老师：

Part1大冶研究区土地使用情况分析

1.监督分类最大似然法

2009年影像分类结果：

分了4类：建筑，水体，裸地，植被

2010年影像分类结果：

2.栅格转矢量

3.变化检测-叠加分析

（1）先提取2009年和2010年分类图的建筑用地：

2009年 2010年

（2）变化检测（叠加分析）

Chang=2010-2009

4.空间查询

变化图层Chang 与建设用地红线进行查询

5.提取结果效果图

：批而未用 :

用而未尽 ：正规使用

：越界开发

： 未批先用

将变化矢量图与2010年影像叠加显示：

Part2：遥感反演与建模

1.数据预处理

1.安装环境卫星数据处理补丁

将ENVI_HJ1A1B_Tools.sav补丁放在…\ITT\IDL\IDL80\products\envi48\save_add目录下。

2.数据读取和定标

主菜单->File->Open External File->HJ-1A/1B Tools，打开环境卫星数据处理补丁后，选择CCD，Input path选择环境卫星数据文件夹，点Search，设置输出路径，勾选“Calibration”“Layer Stacking”，点Apply，如下图：

3.工程区裁剪

由于整景影像范围太大了，进行几何校正之前，裁剪出我们需要的太湖及其周边区域。

1.打开上一步处理好的数据：HJ1B-CCD1- Calbrated_LayerStacking.img

2.完成太湖及其周边区域的裁剪

主菜单->File->Save File As->ENVI Standard，弹出New File Builder面板，单击Import File，弹出 Create New File Input File面板，选中Select Input File 面板中的数据，单击Spatial Subset，弹出Select Spatial Subset面板，单击image 弹出Subset by Image对话框，在其中裁剪出太湖及其周边区域，在几个对话框中单击ok，在New File Builder面板中，单击choose，设置保存路径，输出文件名

HJ1B-CCD1-Cal-sub.img

4.几何校正

(1)打开基准影像TM_baseimage.img，选择Map-> Registration ->Automatic

Registration：Image to Image，选择基准影像TM_baseimage.img的波段4作为匹配波段

(2)选择被配准影像HJ1B-CCD1-20091006-Cal-sub.img ，选择band4作为匹配波段，在提示是否手动选择同名点时，选择否后，弹出Automatic Registration Parameters 面板

(3)在Automatic Registration Parameters面板中，设置下图参数单击ok执行

基于像元值自动寻找同名点。

对于RMS高的点直接删除，水域中的点由于没有固定参考物都删掉，湖岸线的点RMS高的点可以直接删掉，或者在两个影像的ZOOM窗口上将十字光标定位到正确位置，再点击Update进行微调，在缺少控制点的地方手动添加一些控制点。调整控制点，直到总的RMS小于1时，完成控制点的选择，点击Ground Control Points Selection上的File->Save GCPs to ASCII，保存控制点。

(4)在Ground Control Points Selection面板上，选择Options->Warp File

（As ImageMap）选择校正文件。在校正参数面板中，投影参数默认。在X和Y的像元大小输入30，按回车，图像输出大小自动更改。选择输出路径和文件名HJ1B-CCD1-20091006-Cal-sub-jz.img，点击ok进行Image Registration。

(5)打开基准影像和校正影像，在显示校正后影像的窗口中，右键选择Geographic Link命令，选择显示图像的两个窗口，打开十字光标查看校正结果。

5.大气校正

环境小卫星提供了波谱响应函数，以文本形式提供第一列表示波长（nm），后面4列分别表示4个波段对应波长的波谱响应函值。为了做大气校正，需要制作波谱曲线来描述波谱响应函数，用于大气校正（波谱响应函数与CCD型号相对应）（1）制作波谱响应函数

主菜单Window->Start New Plot Window，打开ENVI Plot Window面板，选择File->Input Data->ASCII，导入681_HJ1BCCD1文本文件，波长单位选择Nanometers，单击ok。

在绘制窗口出现了4条曲线，选择Edit->Data Parameters，编辑每条线的名称为b1，b2，b3，b4，便于区别

选择File->Save Plot As->Spectral Library，在Output Plots to Spectral Library窗口中，单击Select All Items，单击ok。

在Output Spectral Library面板中，选择输出路径和文件名“HJ1BCCD1光谱响应”，单击ok，将波谱曲线保存为波谱库文件。

（2）FLASSH大气校正

选择主菜单Basic Tools->Convert Data（BSQ 、BIL、 BIP ）选择校正后的影像HJ1B-CCD1-20091006-Cal-sub-jz.img，在Convert File Parameters中选择BIL,选择Convert In Place：yes，单击ok。