实验二 遥感图像的辐射定标

1、实验二遥感图像的辐射定标1实验目的与意义：了解辐射定标原理使用ENVI软件自带的定标工具定标学习波段运算进行辐射定标2为什么要进行辐射定标，定标的原理是什么？目的：消除传感器本身的误差，确定传感器入口处的准确辐射值。原理：辐射定标是将传感器记录的电压或数字量化值(DN)转换为绝对辐射亮度值(辐射率)的过程，或者转换为与地表(表观)反射率、表面温度等物理量有关的相对值的处理过程。3.辐射定标过程一般有两种方式：第一种：利用计算公式，在ENVI中利用bandmath计算福亮度和反射率。第二种：利用ENVI自带的定标工具进行定标，获取福亮度或反射率。第一种方法：用波段运算得到Radiance和Ref

2、lectance(1)表观辅亮度radiance的计算radiance=(Imax-lmin)/(qcalmax-qcalmin)*(qcal-qcalmin)+lmin其中：radiance-表观辐亮度qcalDN(也就是影像数据本身)；Imax和Imin是从参数表中查询；qcalmax是DN值的最大值,对于TM是8bit来说,qcalmax=255；Qcalmin是DN值的最小值，一般为0即Qcal十LMINa(2)表观反射率的计算P=n*L*d2/(ESUN*cos(0)其中P为表观反射率；L为上一步计算出来的表观辐亮度；d为日地距离，这个数据通过下面的表格中获取；ESUN为大气层外的太

3、阳辐射，也可以说是传感器接收处的太阳辐射;0为太阳天顶角。(这个可以通过影像的元数据获取)在本次实验的数据中radiance=(193+1.52)/255*bl-1.52Reflectance=3.14*(bl)*1.0128A2/(1957*0.7381)步骤如下：打开文件L5120036_03620100819_MTL.txt,点击BandMath,输入(193+1.52)/255*bl-1.52,之后即可计算出辐歸度，文件保存为radianceloAvailableBandsListFileOptionsEQ15120036_03620100819TLxt由題MapInfoELS1200

4、36_C：H-TMoTMr-oTM|-oTMh-oTMp-oTM由題MapInfoL51ZOO36_O36ZO100819JrtTLtxt口TMMeta(Band6)(11.4350)MetaMetaMetaMetaMetaFilatA03620100819_MTLtxtCBand1)(Band2)(Band3)(Band4)(Band5)(Band7)tt).4850)(0.S90)(0.300)(0.8400)(1.6760)(22230)QGrayScaleRGBColorSelectedBandTHMeta(Band6)(11.4350)15120036_03620103819Dims

5、8141x7011(Byte)BSQ坦出塑躲汪妁uo一leq=eu-espue夭s=lnUOQeq=euAbDwssodpd舟哑i：FileSpectralSubletTM(Band0485Cl)：L512ClCl3&_03&2ClllJCl819TLtxtSelectBandstoSubsetTMMetaCBand2)TMMeta(BanA3THMeta(Band4)TMMeta(BanA5THMetaBandT(0.5&90):L512003S.03620100819JHL.txt(0.6&00)：L12003&_03&20100819_MTL.It(0.8400):L5120036.03

6、620100819JTL.txt(1.6760)：L12003&_03&20100819_MTL.It0.2230):1512003S.03S20100819JTL.txtFwnberofitemsselected：1Clear|ImportASCII.AddRangeSelectAllOKCancelFreviousLandsatCalibrationInputFileS-electInputFile：123L512003603620100819MTL.txtL5120036J3620100819JiirL.t淤L512003603620100819MTL.txtFile：C:VsersVu

7、serDesktopLT5120035201CDimsSize:ByteFile_SensorTypeByteOrderFrojectioiiPixelDatiuriWavelengthUpperLeftCorner:lz1Description：GEOTIFFFileImportedintoIWIMonApr11l&：09：062016,GROUP=L1JETAIATA.FILE,ACQUISITIOHDATE二2010-08-19,SPACECRAFT_D=LandsatSENSORJD二巩SU1LELEVATI0K二58.8459801,8141Type7011x6ESQ65,535by

8、tes.Lndsa.tMetadataLandsatTMHostffn-tel)UTM,Zone50Korth30MetersWGS-840.485to2.223MicrometersOK|CancelPreviouz点击OK,之后选择输出为radiance文件，命名为radiance2o打开radiancel和radiance2,右击选择link,将两个图像连接起来,然后右击鼠标，选择CursorLocation,出现窗口的内可以看到两种方法得到的辐射度的数值对比。我们可以看到两者有差距，但是很小，这是因为ENVI软件内的数据不是最新的，计算出來会有一定的误差。蛍/?TJ弦序FCursorL

9、ocation/Value796965.00E?3794895.OOHMeters34?51014N,L20?3?2S.85E#1Data：9.935532#2D吐父99.569138FileOptionsDisp#2(584L4T75)Scrn：R:102G:102B102Frojeetion：UTNZone50lUrthMap:LL:DispDisp下面进行反射度的计算，同样的，点击PreprocessingCalibrationUtilitesLandsatCalibration选择所要校正的波段，点击OK,之后选择输出为reflence文件，命名reflence2oXIwnCtoil在

10、BandMath输入3.14*(bl)*1.0092*1.0092/(1957*sin(58.8459801)*3.14计算出辐射度，之后保存为reflencelo打开reflencel和reflence2,将其连接起来,通过CursorLocation可以看到两者的对比。QENVlClassic-cFUeBockTookClocficationTransformFitteSpectralMapVectorToporophkRodorWindowHelpOCursorLocation/V.XFileOptionsDisp#2(2536,4286)Scrn：R:57G:57B:57Firojection；UHitZone50IhrthMap：697635.00E,3809565.0011MetersLL:34?丫311广兀