遥感影像辐射校正.doc

遥感影像辐射校正 实验报告 姓名：杨磊 学号班级：11遥感 指导教师：罗小军 2013/6/13 目录 一、 实验名称： 3 二、 实验概况： 3 三、 实验目的： 4 四、 实验步骤： 4 （一）、绝对大气校正： 4 1、按照实习要求将quickbird_multi_2006_xipu.img的第一波段提取出来 4 2、将第一波段图像以“JPG”格式输出： 6 3、在matlab中进行处理： 7 （二）、相对大气校正： 10 1、把quickbird_multi_2006_xipu.img中图像的第一和第四波段图像提取出来： 10 2、回归分析： 11 4、 在matlab中实现相对大气校正： 12 5、 校正结果： 12 （三）、多时相影像匹配法： 13 1、分别提取spot5_10multi_2003_xipu和spot5_10multi_2005_xipu的第三波段（绿色）的影像： 13 2、多时相影像匹配原理： 14 3、在matlab中实现： 15 4、结果： 16 五、实验体会： 17 实验名称： 遥感影像辐射校正实习： 绝对大气校正； 相对大气校正——回归分析法； 多时相影像匹配法。 实验概况： 影像数据：quickbird_multi_2006_xipu.img spot5_10multi_2003_xipu.img spot5_10multi_2005_xipu.img; 实验所用软件：ERDAS9.2 matlabR2010a； 实验地点：X4245 ； 实验时间：6月9日——6月17日。 实验目的： 复习巩固课堂上所学的对遥感图像的辐射校正，掌握这些校正方法的基本原理； 实际体验对遥感图像进行绝对大气校正、相对大气校正的回归分析法和多时相影像匹配法； 考察了对matlab的基本使用，复习掌握一定的编程能力。 实验步骤： （一）、绝对大气校正： 1、按照实习要求将quickbird_multi_2006_xipu.img的第一波段提取出来： 单击，出现如图1.1 1，在当中选择“Utilities”,在弹出的图1.1 2中选择“Layer stack” 图1.1 1 图1.1 2 图1.1 3 在图1.1 3中的input里面找到quickbird_multi_2006_xipu.img文件，输入output File名quickxipu.img。在Layer一栏选择“1”，点击Add，选择输出无符号整数16位，点击OK，则将第一波段图像提取了出来。 2、将第一波段图像以“JPG”格式输出： 单击，弹出图1.1 4 图1.1 4 在Type中找到JFIF（JPEG）格式，在Media中选择 File，在input File 中找到quickxipu.img，单击OK出现图1.1 5，点击OK 图1.1 5 转换成功。 图1.1 6 3、在matlab中进行处理： 打开matlab，将目录置于图1.2 1 图1.2 1 在明亮窗口输入“m=imread(quickxipu.jpg);”，回车后输入“imshow(m)”，得到图1.2 2 图1.2 2 在以前地物光谱实习中所测过的多个地区的叶子反射率找出来，分别选取了五处的叶子在波长为450~520nm间的反射率，并在Excel中求得各个地方的平均值{[0.0787,0.0739,0.0920,0.0963,0.1076}，在图上分别选取这五个地区的五个相应地物的点，记下其index值{74,87,95,94,97}，选取方法如图1.2 3 图1.2 3 在matlab命令窗口输入“edit”编写一M文件： function m=aac(T,B,P) X=double(T); [p,q]=size(X); a=polyfit(B,P,1);%% ·′é??êoˉêy?í?a￡of=1/a*B-b/a; N=zeros(p,q); for i=1:p for j=1:q N(i,j)=1/a(1)*X(i,j)-a(2)/a(1); if N(i,j)0 N(i,j)=0; else N(i,j)=N(i,j); end end end m=double(N); 在命令窗口输入“a=[0.0787,0.0739,0.0920,0.0963,0.1076]; ”，回车，