遥感图像辐射校正实习报告.docx

遥感影像辐射校正实习报告姓名：__________XXX__________学号：_________XXXXXXXX_______班级：________XXXXXXXXXX______ECHODIDTHISFORYOU2013/6/10目录一、实验目的1二、实验内容1（一）绝对大气校正1（二）相对大气校正——回归分析法1（三）多时相影像匹配法1三、实验步骤2（一）绝对大气校正2（二）同时多波段图像的大气校正5（三）多时相影像匹配法8四、实验体会10一、实验目的进一步巩固、掌握遥感影像绝对及相对大气校正基本方法。二、实验内容（一）绝对大气校正以实测或从光谱数据库中查得的光谱数据采用基于ELC的大气校正方法对交大犀浦校区2006年QuickBird多光谱影像中的第1波段影像进行绝对大气校正。基本步骤：1、从QuickBird多光谱影像文件中提取出第1波段影像；2、从影像中判读出一些典型地物；3、从以前实测的光谱物据或光谱数据库中读出步骤2中判读出的那些地物对应QuickBird第1波段的反射率值；4、基于步骤3的反射率值采用基于ELC的大气校正方法对交大犀浦校区2006年QuickBird多光谱影像中的第1波段影像进行绝对大气校正。（二）相对大气校正——回归分析法以交大犀浦校区2006年QuickBird多光谱影像中的第4波段影像为参考采用回归分析法对第1波段影像进行相对大气校正。（三）多时相影像匹配法以交大犀浦校区2003年SPOT5多光谱影像中的第3波段（绿波段）影像为参考采用多时相影像匹配法对交大犀浦校区2005年SPOT5多光谱影像中的第3波段（绿波段）影像进行相对大气校正。三、实验步骤（一）绝对大气校正1.从多光谱影像文件中提取出第1波段影像（1）打开在其中选择Interpreter->Utilities->LayerStack->InputFile:quickbird_multi_2006_xipu.img->OutputFile:2006-1.tif->Layer选择为第一波段->单击Add->DataType中的Output设置为无符号十六位->单击OK如图1所示图1从多光谱影像中提取第1波段的影像的方法图2.获取特征地物的灰度值与反射率（1）在中读取第1波段的影像图片并显示具体语句如下：1).A=imread('2006-1.tif')%将图像信息读取到矩阵A中2).imshow(A[min(min(A))max(max(A))])%在中显示读出的图像如图2所示图2在中显示quickbird_multi_2006_xipu.img的第1波段影像图（2）在图2所示的图形显示对话框中选择按钮后点击地面特征点并从获取其灰度值。在本次实习中我选取了植被、石头、土壤的为特征地物对应光谱实习中获取三者的反射率得到数据如表1所示表1特征地物的反射率与灰度值地物树叶石头土壤反射率（%）6.7312.654.92灰度值3043643273.基于大气校正模型的建立(1)由大气校正模型的表达式式中：为第像元灰度值为特征地物的反射率值为回归系数为常数将原始数据带入上式后可得根据最小二乘原理可得的计算公式为将数据带入后得最终得到1波段遥感图像的大气校正模型为可得校正后各对应地物反射率的计算表达式为经过校正后的图像如图3(a)校正前的图像(b)校正后的图像图3绝对大气辐射校正前后的对比图分析：通过校正前后的图像对比后可以看出经过大气校正后的图像明显要比校正前的图像清晰因此可以认为本次实习基于ELC的大气校正方法建立的大气校正模型有效。（二）同时多波段图像的大气校正1.从多光谱影像文件中提取出第4波段影像(1)打开在其中选择Interpreter->Utilities->LayerStack->InputFile:quickbird_multi_2006_xipu.img->OutputFile:2006-4.tif->Layer选择为第4波段->单击Add->DataType中的Output设置为无符号十六位->单击OK如图4所示图4从多光谱影像中提取第4波段的影像的方法图(2)将多光谱影像第1、4波段影像显示如图5(a)第1波段的影像(b)第4波段的影像图5进行相对大气校正的影像图2.相对大气校正模型的建立（1）根据相对大气校正模型的回归方程式中：为第1波段的灰度均值为第4波段的灰度均值由最小二乘原理得到的计算方法如下（2）在两波段图像中选择从最亮到最暗的一系列目标对每一目标的两波段亮度值进行回归分析。在两波段中选取的目标点的灰度值如表2表2样本点的坐标值及灰度值坐标点位(685772)(681753)(622684)(69478