康哥定量遥感实验报告

定量遥感实验报告班级：地信1302：高永康学号：12实验一:辐射信息获取与大气校正一实验目的与意义初步了解目前主流的遥感图象处理软件ERDAS，ENVI的主要功能模块。掌握LandsatETM遥感影像数据，数据获取手段。掌握Erdas遥感影像辐射信息获取。加深对遥感理论知识理解，掌握遥感大气校正方法。实验容与遥感影像数据LE097EDC002003年4月7日LandsatEMT+影像对其进行大气校正处理实验步骤准备阶段：对遥感图像先定标，BasicTools-Preprocessing-CalibrationUtilities-LandsetCalibration并编辑头文件3.1简化黑暗象元法大气校正打开待校正图像文件。在主菜单中，选择BasicTools->Preprocessing->GeneralPurposeUtilities-〉DarkSubtract，在文件选择对话框中选择待校正图像文件，单击0K按钮，打开DarkSubtraction(3)在DarkSubtractionParameters面板中，确定黑暗像素值包括三种方法(SubtractionMethod):波段最小值(BandMinimum)ROI的平均值(RegionOfInterest)自定义值(UserValue)

(4)在OutputResultto中选择File以与相应的输出路径和文件名，单击OK执行操作。(1)在ENVI主菜单中，选择以下方式启动lBasicTools-〉Preprocessing-〉CalibrationUtilities->QUickAtmosphericCorrectionlSpectral->QUickAtmosphericCorrectionlSpectral->Preprocessing-〉CalibrationUtilities->QUickAtmosphericCorrection(2)打开QUickAtmosphericCorrectionParameters面板在SensorType中选择相应的传感器类型，选择文件名和路径输出。1•设置辐射率转换因子相应参数在MTL文件中查询。并选择Multispectralsetting选实验小结因为我们从卫星上获得的影像实际上已经受到了各种因素的干扰，如地表朗伯面反射、大气性质不一与大气多次散射辐照作的用和邻近像元漫反射的作用，这些都使得我们获得的图像有一定的偏差，因此在使用这些图像之前，需要进行校正处理，如大气校正，辐射校正等。黑暗像元法：是一种古老、简单的经典大气校正方法。它的使用前提是在假设待校正的遥感图像上存在黑暗像元、地表朗伯面反射和大气性质均一，并忽略大气多次散射辐照作用和邻近像元漫反射作用下，反射率很小（近似0）的黑暗像元由于大气的影响，使得这些像元的反射率相对增加，可以认为这部分增加的反射率是由于大气影响产生的。这样，将其他像元减去这些黑暗像元的像元值，就能减少大气（主要是大气散射）对整幅影像的影响，达到大气校正的目的。整个过程的关键是寻找黑暗像元以与黑暗像元增加的像元值。快速大气校正工具（简称QUAC）:其自动从图像上收集不同物质的波谱信息，获取经验值多光谱的快速大气校正（图4）。它得到结果的精度近似FLAASH或者其他基于辐射传输目前它支持的多光谱和高光谱波谱围是（0.4〜2.5“m）。Flaash大气校正：FLAASH采用了M0DTRAN4+辐射传输模型，该算法精度高。任何有关影像的标准MO型和气溶胶类型都可以直接使用。通过影像像素光谱上的特征来估计大气的属性，不依赖遥感成像时同步测量的大气参数，可以有效地去除水蒸气/气溶胶散射效应，同时基于像素级的校正，矫正目标像元和邻射的“邻近效应”。并且对由于人为抑止而导致波谱噪声进行光谱平滑处理。作为结果，除了真实地表反射率到整幅图像的能见度、卷云与薄云的分类影像、水气含量数据。5.5.1遥感影像大气校正为什么要定标?定标是将传感器所得的测量值变换为绝对亮度或变换为与地表反射率、表面温度等物理量有关的相对值的处理过程，目的是消除传感器本身的误差。大气校正就是将辐射亮度或者表观反射率转换为地表实际反射率，目的是消除大气散射、吸收、反射引起的误差。辐射定标可以说是大气校正的基础，一切都是为了获得更精确的遥感影像。5.2遥感大气校正主要影响因素？1遥感所利用的各种辐射能与地球大气层发生的相互作用辐射能的散射、吸收，而使能量衰减。大气的对不同波长的光有选择性的衰减作用。太阳-目标-遥感器之间的几何关系不同，使其穿越的大气路径长度不同，从而使图像中不同地区地物的像元灰度值所受大气影响程度不同。同一地物的像元灰度值在不同获取时间所受大气影响程度不同。实验二:地物识别和定量反演一•实验目的理解基于遥感光谱曲线的地物识别原理，掌握遥感定量反演模型和方法。学习ENVIBandmath工具进行地表温度反演的过程。加深定量遥感反演知识理论，增强遥感软件平台操作能力。二•实验容.地物识别与定量反演基础.Bandmath基础.地表温度反演实验数据研究区NDVI数据，ETM数据经数据预处理后并重采样为60米分辨率的NDVI数据TM-NDVI-60m.img研究区第六波段辐射数据，ETM6经过传感器定标，几何校正，大气校正，裁剪后的辐射数据TM6-rad-subset-jz-xiangfan.img

三•实验步骤（一）地表比辐射率计算植被覆盖度计算ENVI:主菜单-BasicTools-BandMathle公式：(blgt0.7)*l+(bllt0.)*0+(blge0andblle0.7)*((bl-0.0)/(0.7-0.0))注意：bl选择NDVI图像操作如下：namet计尉』tnamet计尉』tW1N7(C;)kPiograrnFilei(ka6]zhibeiRjgazhiberfijgifiiiduidu.hdrlJ收审夹■臺㈢已最近访月的進忘仁逗[叫：2hibeilugaidu结果如下：联1联1BandMath((blgtd-7)siH-CblltO.^O<-(blgetl..宴F-ihsOverlayEnhanceTools-Window2•地表比辐射率计算ENVI:主菜单-BasicTools-BandMath公式：(blle0)*0.995+(blgt0andbllt0.7)*(0.9589+0.086*b2-0.0671*b2辽)+(b1ge0.7)*(0.9625+0.0614*b2-0.0461*b2J)注意：b1选择NDVI图像；b2选择植被覆盖度操作如下：⑨BandMath[—^■1⑨BandMath[—^■1FreviousBand ExTruwTons；(bl“0)*0，旳5+(bl寸QblIt0,T)*(结果如下：计算一样温度下黑体辐射亮度值根据辐射传输方程和普朗克定律，温度为T的黑体在热红外波段辐射亮度B(Ts)表示为：B(TS}=[LJL-LX大气向上和向下辐射亮度；地表比辐射率地表真实温度；普朗克定律推到得到的黑体在Ts的热辐射亮度；大气在热红外波段透过率；热红外波段亮度值。•参数可得到t为0®川为3・39叭肿》：• 为5/2操作如下：：：》BandMathFreviouzE；iiniMathEspressicmz：[bl£.闕-Cl.时〔1-吃）怖.12]旳.帥垃OKC：=lTl<:el结果如下：（三）反演地表温度获取温度为Ts的黑体热红外波段辐射亮度后，根据普朗克公式反函数，求得地表真实温度T5二K2/ln(K17*(7i)+L)对于ETM匕K1=666.09 -sr-firn)K2=1282.71K注意’自然对数IDL表达式创og计算得到地表真实温度，并要转换为摄氏温度操作如下：

专题图输出•显示温度反演结果；•利用ArcGIS生成专题图，并根据温度围划分4个区间：-39C以上，红色-35-39C，黄色-30-35C，绿色-低于30C，蓝色寰樊市地衣温度反演专题图17.HBIt：料国人t地忖】盹炳应除

7^!^31^^200212

4倆时心auiMi*1Ji17.HBIt：实验小结热辐射是物体由于具有温度而辐射电磁波的现象。一切温度高于绝对零度的物体都能产生热辐射，温度越高，辐射出的总能量越大，短波成分也就越大。物体在向外辐射的同时，还会吸收从其它物体辐射来的能量，物体辐射或吸收的能量与它的温度，表面积，黑度等因素有关。思考：(1)地表温度反演利用的是什么方法，有哪些优缺点？1.单窗算法。其适用于只有一个通道的热红外遥感数据。其将大气平均作用温度和大气下行平均作用温度合二为一，简化了热辐射传输模型。缺点是与大气传输模型一样，在参数化过程中使用的数据仍然是标准化的，并没有使用实施大气轮廓数据。分裂窗算法。分裂窗算法是目前应用最广泛的温度反演方法，分裂窗方法是利用大气窗口10〜13fim里的两个相邻通道（11“m附近和12im附近）的大气吸收作用不同，通过两个通道测量值的组合来消除大气影响。3.温度，比辐射率分离算法。温度比辐射率分离算法（TES）针对高光谱分辨率的热红外测量数据，利用地物热红外光谱的一些共性特点作为先验知识或者约束条件，从一次观测中同时反演温度和光谱发射率。它集合了三种已有方法（NEM,RAT,MMD）的优点和一些新的特点。且第一次提出发射率的无偏和精确估计改进了陆面温度的估算。（2）哪些方法可以有效的提高温度反演精度？单窗算法，分裂窗算法，TES算法等。实验三:遥感反演模型•实验目的与意义初步了解ERDAS主要功能模块。掌握ErdasModeler遥感建模方法，遥感定量反演算法设计。加深对ETM影像植被反演、温度反演方法学习。实验容与遥感影像数据LE097EDC002003年4月7日LandsatEMT+影像

Flaash大气校正影像数据三实验步骤1•数据准备1.1Erdas数据转换：ManageDataTmportData将Envi数据转换为IMG数据格式N】mportFomn冃N】mportFomn冃t;InputFile:OutputFile:EN训（Di用utR曲d] *|嫌|訝作业F定量遥感/新建文件夹/I71123039J〒虚b30.inngOK.Clo?eDataView...Help点击OK即可。Rorf Rorf FlocasTMeStalePiOCiffiS1eW^pace2ImageCommandDONE-LickDisniR5boRflmmrE1呼■3ImpuCcftinHntlDONE-OickDiyinb^toRerwc 4InpoilinaIMAGINEImawD抽DONE•□冰DiitiiisstoRerrovaion^5InpailrigIMAGINEImage-DeJa□LIME-匚1匚kDismEEtaflera<Eions一6M«Wer•lumingmedstI^BinrrgtprrdlDONE-Click少帥馭也Removi?Cjncd□o粧 |HelpCjncd1.2Erdas多波段合成Raster-Spectral-LayerStackLayerStackLayerStackLnyeiStackLnyeiStackStazkmultiple^.allysirtghband]W:m-aa-sasban.ds/^usra-in^taaurgteoirtpitteuitti-bapdimag-fsh.ExarnofwbsHthisu:ommonrlyj-s-d:n-zhdeKON-05muJth-pectulm-agcryd-lr«™dass-epar-at-sTIFForN1TFfiles[onebandpsrfile]-artozomblriemufltjp]-denvateimaos■m±wre,itextirrc,mdspendsnt-cmportn±saetc.)into-asuxgkmu.^ti-handimaj-toimprcvsTatsrf'zaton-acc-rr-BCj1.2.影像裁剪Raster-Subset&Chip-CreateSubsetImage从AOI裁剪RasterVector T町個ini T暫口Ibex HelpPariiCPassifica-tionCF▼ DetCLASS.N>呈MTFCopy=3selectsdporticin,ors-ubK-t,ofminputdatafileinto*anioutputdat-afils.Etabmere1/triincBtadbasedanbomndin^-boK,laysraand^atatjrpe•mdrtaitretchingisdappled.I±»l^ss^ahrfyoui'atshistransformth■■imaceckrtatype-结果如下3.各类植被指数计算3.1未校正下Raster—ClassificationNDVI或Indices计算NDVI和比值植被指数RVIDrawing3・2大气校正(1)转成BIL格式File-Ba^i匚T口口Is匚iBssiFicatrnnTran^larmFResizeData(Spathl/Spectral)SubsetDataviaROIsRotate/FlipDotaiLayerStackingStretchDsrfcaConvertFileParametersenv渚嘩支件StatisticsSpatialSaati^ticBChangeDefectionMeasurementToolStretchDsrfcaConvertFileParametersenv渚嘩支件StatisticsSpatialSaati^ticBChangeDefectionMeasurementToolBandMathSpectralMath£eqmentationImageRegionOflEure^tMosaickingMaskingPrenp-rocessing(2)FLAASH大气校正T「石n^FpriTi1]FilterSpectralMapVtImage-Sha『ipuningBandIRiotiosPrincipalComponents-IndependentCompDn&ntsMNFRotationColorTrsnsfaimiEDe匚o-rreistionStretchPhotographicStretchSaturationStretchSyntheticColorImage「NDVITasseledCapt-0fL&&ih□PLUSH(tajur田md]：L?;卜□PLUSH（L呼■「ipKTiJ1:匚7;卜□PLMSH田KTiJ]:匚:;卜□PLM!ZH田md]:L?;卜□PLUk^H(Bu.d】:LTI卜□PLUl^H(L*y・rCBKrid】:LTFileOplian-s建模操作（同时计算出