landsat遥感影像地表温度反演教程大气校正法

1、基于辐射传输方程的Landsat数据地表温度反演教程一、数据准备Landsa8遥感影像数据一景，本教程以重庆市2015年7月26日的二行列号为(128,049)影像(LC81280402016208LGN00)为例。同时需提前查询影像的基本信息(详见下表)标识日期采集时间中心经度中心纬度LC81280402016208LGN00,016/7/263:26:56106.11288I30.30647注：基本信息在影像头文件中均可查询到，米集时间为格林尼治时间。二、地表温度反演的总体流程/LaridsatSOUniRSUI1轴射定标"OU士气校正*IInum计茸/-I0辐射亮一图/地表比嗣

2、射率十箱，,大气利阳冬教/Jj同温度卜黑体朝射隹度十算地表一度计算三、具体步骤1、辐射定标地表温度反演主要包括两部分，一是对热红外数据，二是多光谱数据进行辐射定标。(1)热红外数据辐射定标选择RadiometricCorrection/RadiometricCalibration0在FileSelection对话框中,选择数据LC81230322013132LGN02_MTL_Thermal,单击SpectralSubse选择ThermalInfraredl(10.9),打开RadiometricCalibration面板。|Radi-:ietricCoirex：ti:m*jApplyGilD

3、iuJnriE.ef3nJItCoreeLH川式cCdTre-tiiinT8s.dkAtnoBEHicTi4CcrrcctluCQOACi禽19hfattA.WEITHIS-rs-s-rrtcfcEllLTilnttlmCsMTKtl:>nBorkSuhsra.-:t!«xn加印RTP*l”用鼬FBnlsslvrliT占lohuPrslilialE<3n£sJivlIj1E'ctl.iillza/tLcKi.BnirfltviiyftrfcrmeCharavei：jErwirlcsiLtriT口1Itisl3ng4：fnpSrl>za.lL-Lij

4、eCnjruteFactits14PlAfFifldCiibrailjTAFRHIfts*C-illUraTioriLciieResiEuiiL:;i:miJiCrsstl口n«1l4metric：CallbratlmIEasterVana：；Ei=ntf匚看tWMFCiTtpiexDana"cor/vmimcririvicU«2rrateC:rczlniti-SysteraStadnc、4mRMIUetaFaleDa”neverj1sta-Ejpe:3fcTtL15ties-11=£trips.EEJQuftieWgmOnuitputHataType(

5、FloatScaleFactor1-00儿口口1了?I_AA.iKScitings计算会报错。保持默认1即可。Scalefactor不能改变，否则后续OutputFiLcnaTie:risni>,Tc5uitOKCincfil(2)多光谱数据辐射定标选择要校正的多光谱数据“LC81230322013132LGN02_MTL_MultiSpectral进行辐射定标。因为后续需要对多光谱数据进行大气校正，可直接单击如下图。ApplyFlaashSettings,RciometricQlibruti0nCallLra.lld：iTvptKadlancgtm;工iciTETlesveE工LOut

6、putmt州Typ尸Float注意与热红外数据辐射定标是的差ScslffacroT0.10别，设置后Scalefactor值为0.1。ApplyFLAASHSettingsdrtpntFilename:iDisrla/resu)10KCanru】2、大气校正本教程选择Flaash校正法。FLAASHAtmosphericCorrection,双击此工具,打开辐射定标的数据，进行相关的参数设置进行大气校正。注意：如果在多光谱数据辐射定标时Scalefactor值忘记设置，可在本步骤中打1）InputRadianceImage打开辐射定标结果数据;开辐射定标数时设置singlescalefaceo

7、r值为0.1,若已设置，则默认值为1即可。2）设置输出反射率的路径，由于定标时候;3）设置输出FLAASH校正文件的路径，最优状态：路径所在磁盘空间足够大;4）中心点经纬度SceneCenterLocation自动获取;5）选择传感器类型：Landsat-8OLI；其对应的传感器高度以及影像数据的分辨率自动读取；6）设置研究区域的地面高程数据；7）影像生成时的飞行过境时间：在layermanager中的Lc8数据图层右键选择ViewMetadata,浏览time字段获取成像时间；注：也可以从元文件“LC81230322013132LGN02_MTL.txt”中找到，具体名称：DATE_ACQU

8、IRED=2013-05-12;SCENE_CENTER_TIME=02:55:26.6336980Z;8）大气模型参数选择：Sub-ArcticSummer（根据成像时间和纬度信息选择）；9）气溶胶模型AerosolModel:Urban,气溶胶反演方法AerosolRetrieval:2-band（K-T）;10）其他参数按照默认设置即可。11)多光谱参数设置中，K-T反演选择默认模式：Defaults->Over-LandRetrievalstandard(600:2100)波谱响应函数：默认指向.ProgramFilesExelisENVI51classicfilt_funcla

9、ndsat8_oli.sli把它重新指向：.ProgramFilesExelisENVI51resourcefilterfuncslandsat8_oli.sli注：这是因为ENVI5.1版本的一个小bug,即Classic中的L8的波谱响应函数不正确，另外一个一劳永逸的方法是：将“ExelisENVI51resourcefilterfuncs”中的Landsat8_oli.sli和Landsat8_oli.hdr两个文件拷贝覆盖：”.ENVI51classicfilt_func”中的两个文件。否则SWIR1波段大气校正后的结果全为00见下图工MuftispectralSettings12)高

10、级参数设置：根据内存大小设置TileSize(Mb):100(8g物理内存),其他参数默认即可，详细见下图经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响，植被的波谱曲线趋于正常。3、地表比辐射率计算(1)植被覆盖度计算计算植被覆盖度Fv采用的是混合像元分解法，将整景影像的地类大致分为水体、植被和建筑，具体的计算公式如下：FV=(NDVI-NDVIS)/(NDVIV-NDVIS)其中，NDVI为归一化差异植被指数，取NDVIV=0.70和NDVIS=0.00,且有,当某个像元的NDVI大于0.70时，FV取值为1;当NDVI小于0.00,FV取值为0。利用ENVI主菜单->

11、BasicTools->BandMath,在公式输入栏中输入：(b1gt0.7)*1+(b1lt0.05)*0+(b1ge0.05andb1le0.7)*(b1-0.05)/(0.7-0.05)b1:选择NDVI图像(2)地表比辐射率计算根据前人的研究，将遥感影像分为水体、城镇和自然表面3种类型。本专题采取以下方法计算研究区地表比辐射率：水体像元的比辐射率赋值为0.995,自然表面和城镇像元的比辐射率估算则分别根据下式进行计算：surface=0.9625+0.0614FV-0.0461FV2building=0.9589+0.086FV-0.0671FV2式中，surface和buil

12、ding分别代表自然表面像元和城镇像元的比辐射率。利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入：(b1le0)*0.995+(b1gt0andb1lt0.7)*(0.9589+0.086*b2-0.0671*b2*b2)+(b1ge0.7)*(0.9625+0.0614*b2-0.0461*b2*b2)b1:NDVI值；b2:植被覆盖度值。得到地表比辐射率数据。4、计算相同温度下黑体的辐射亮度值卫星传感器接收到的热红外辐射亮度值L人由三部分组成：大气向上辐射亮度Lt,地面的真实辐射亮度经过大气层之后到达卫星传感器的能量；大气向下辐射到达地面后反

13、射的能量。卫星传感器接收到的热红外辐射亮度值的表达式可写为(辐射传输方程)：LX=£s)+(1-£)L；-r+LT这里，e为地表辐射率，TS为地表真实温度，B(Ts)为普朗克定律推到得到的黑体在Ts的热辐射亮度，t为大气在热红外波段的透过率。则温度为T的黑体在热红外波段的辐射亮度B(Ts)为：B(Ts)=L-AL1-r-£1)L；/t-£在NASA官网(/)中输入成影时间以及中心经纬度,则会提供上式中所需要的参数。本专题输入的数据是重庆市地区2016年7月26日格林尼治时间时间03:26,平均气温32

14、.601C,气压961.841MP,相对湿度63.689%,Landsat8OLI影像，影像中心的经纬度为：30.30647N,106.11288E宵arzire1&TSfitr卬上ich丁ooishcalculalealmEspnErictransirissLanandLinzradiance重庆属中纬度夏季天气*2S种r：A-jj-Sisur+is-isatnocnhericproflIsforcloseciiniegerIit/1oneLTgeinterpclaledatnoaphiericprofileforsi.venlat/loc选择基于精确信息数Usendd"l

15、atiTudesuaDrierstajidtardspherefarlaspertiro.sphjeric&Tofile上tJsemd-latitudcwinterstandardatmosphericfor叩pciQ(JaeL加dadt-BTIREBand1口mpuctrNlremcinflu匚忖口三Use时iida/tNTEmd£当白巴七tral让巴耳式1!88£亡urvrLandsat8数据,波段10LI用亡Bmd6片口上仁七星上七时psi片七也山肃七OutpLilucilympiuflle,dunolcalcalssteeffsctivsradih出壮片Dp

16、tiaoal：SurfaceConiditicnsJf血或M声二十比厂k疗'he告两力£whs*joc命alerted匚的心十工口而willjtf1msIf曲知e/K白？rcuMz打出i/aJ2t讨端士匕曲营友婚十加良力才(?门?，Altitiade(km);D391FYcssurc(inb)：75BABSTeneraturo(C)：131.02Rei5tiveHunidit第(31)：6849ResultswiLlbrs-nttih-follo&TiTLgaidress：Email：CMrill献CiMrhelfls得到下图参数图大气在热红外波段的透过率弋为0.41,大气向上辐射亮度La5.23W/(m2sr仙m)大气向下辐射亮辐射亮度L功7.52W/(m2sr。仙m)利用ENVI主菜单->BasicTools->BandMath,在公式输入栏中输入：(b2-5.23-0.41*(1-b1)*7,52)/(0.41*b1)b1:60m分辨率的地表比辐射率值；b2：表示热红外波段的辐射定标值。得到了温度为