遥感影像辐射校正实习报告.docVIP

遥感实习报告（4） 遥感影像辐射校正 专 业： 班 级： 姓 名： 学 号： 成 绩： 指导教师： 2013年 月日 目 录 一：实验目的…………………………………………………3 二、影像数据…………………………………………………3 三、实验内容…………………………………………………3 四、实验步骤…………………………………………………4 （一）、绝对大气校正……………………………………4 （二）、相对大气校正——回归分析法…………………7 （三）、多时相影像匹配法………………………………9 五、心得体会…………………………………………………13 六：程序设计…………………………………………………14 一：实习目的： 进一步巩固、掌握遥感影像绝对及相对大气校正基本方法。 二：影像数据：交大犀浦校区200年2005年影像 交大犀浦校区200年影像蓝450-520nm)；波段2—绿（520-660nm）；波段3—红（630-690nm）；波段4—近红外（760-900nm）。 三：实习具体内容：交大犀浦校区200年 基本步骤： （1）、从QuickBird多光谱影像文件中提取出第1波段影像； （2）、从影像中判读出一些典型地物； （3）、从以前实测的光谱物据或光谱数据库中，读出步骤2中判读出的那些地物对应QuickBird第1波段的反射率值； （4）、基于步骤3的反射率值，采用基于ELC的大气校正方法，对交大犀浦校区200年交大犀浦校区200年 以交大犀浦校区200年影像交大犀浦校区200年影像QuickBird多光谱影像文件中提取出第1波段影像 启动ERDAS软件在Viewer #1中打开影像：quickbird_multi_2006_xipu..img 图4.1-1 单击Layer Stack Layer，选择为第一波段单击Add-Data Type 中的Output，设置为无符号十六位，单击OK，如图4.1.2所示： 图4.1-2：多光谱影像文件中提取出第1波段影像 2:利用Matlab从影像中判读出一些典型地物，得到它的灰度值和反射率。 :在Matlab中读取第1波段的影像图片并显示，具体程序语句如下： I=imread(2006-1.tif) %将图像信息读取到矩阵I中 m=min(min(I)) %读取矩阵I灰度的最大值 n=max(max(I)) %读取矩阵I灰度的最小值 imshow(I,n,m) %在matlab中显示读出的图像 得到quickbird_multi_2006_xipu..img的第一波段影像的tif格式图，如图4.1.3: 图4.1-3:quickbird_multi_2006_xipu..img的第一波段影像 （2）：在matlab中读取一些地物特征点的灰度值，正如图4.2.1所示，利用DataCursor工具可得到： 这样读出Index，就是这一点的灰度值，在本次实习 中，我们利用的是树叶、土壤与石头在光谱实习中测得的对应QuickBird第1波段的反射率值，如表4.1-4： 表4.1-4:各地物的光谱反射率 植被（树） 石头 土壤 灰度 302 403 310 反射率（%） 5.817 10.537 4.021 3：采用基于ELC的大气校正方法，对交大犀浦校区200年 （BV为地物点亮度值，为反射率） 即： （二）相对大气校正——回归分析法 以交大犀浦校区200年 是TM4的亮度值的均值，是TM1的亮度均值 利用表4.2.1的数据在matlab中进行计算，得到回归系数： b= 0.3963 a=163.7569 即回归方程为： TM1波段图像的改正值为： 依据回归方程对TM1波段影像的每一点进行改正，可得校正前和校正后影像，如图：4.2-2