式中,三为混合像元的地表比辐射率;Pv为植被覆盖率;Rv为植被的温度比率;Rm为建筑表面的温度比率;=v表示植被法地表比辐射率,fm表示建筑表面的地表比辐射率;*表示辐射校正项。
根据覃志豪经验公式,匸v=0.986;・讣=0.972
根据Aster提供的常用地物比辐射率光谱库Landsat8数据特点,对于B10波段,匸v=0.98672;m=0.96767。
对于B11波段,匸v=0.98990;Tm=0.977515°
在地表相对较平整的区域,一般可取de=在地表起伏较大的区域,de可根据植被的构成比例估计。
研究表明热辐射相互作用在植被与裸土分别占一半时达到最大,所以可以根据如下经验公式来估计de:
当Pv<0.5时,de=0.0038Pv
当Pv>0.5时,de=0.0038(1-Pv)
当Pv=0.5时,de=0.0019
RV和Rm的确定采用覃志豪等提出的经验公式计算:
RV=0.9332+0.0585*PV
Rm=0.9886+0.1287*PV
PV为植被覆盖率可用下式计算:
PV=(NDVI-NDVIs)/(NDVIv-NDVIs)
其中NDVI是归一化植被指数,NDVIV、NDVIS分别是植被和裸土的NDVI值,由于没有详细的区域植被和土壤光谱资料,采用覃志豪经验值以NDVIV=0.70和NDVIS=0.05来进行植被覆盖度的近似估计。
1.123大气平均作用温度的计算
对于大气平均温度,覃志豪等也提供了一个经验公式,该公式是根据Modtran所提供的标准大气推导出来的,如下表所示,其中To为近地表温度(K)。
3-4大气平均作用凰度与地面附近气温舵关系腳
人气的剖而类醴
7;芍珀拘关垂
哭国1976均大气
7,^2505Ta
热帶屮均丸气
T.=17.9765+0.51715r0
屮纬度旻李半均人气
^=1ClO110+0.9?
«21r0
屮纬度冬季平均大气
I;-19.2704+0.91118rD
般研究区在重庆的,选择中纬度夏季平均大气
1.124大气透射率的计算
由于单窗算法对输入的参数较敏感,因此需要较精确的大气透射率参数来反演地表温度。
影响大气透射率的因素较多,如气压、气温、气溶胶含量、大气水分含量、03、C02、CO、NH4等对热辐射传导均有不同程度的作用,从而使地表的热辐射在大气中的传导产生衰减。
因此,准确的大气透射率求算比较复杂,需要较详细的大气剖面数据,但在实际应用中很难获取详细的大气剖面数据。
研究表明大气水分含量对大气透射率的变化影响较大,而其它因素由于很难在短时间内发生较大变化而对大气透射率的影响较小。
因此,水分含量就成为大气
透射率估计的主要考虑因素。
覃志豪等根据大气中的水汽含量动态变化对大气透射率的影响
程度进行了研究,并对大气透射率进行了估计,如表3.1。
其中:
w是
大气水汽含量。
以重庆为研究区的研究一般选用是气温较高时的估算公式。
大气透射率估计方程
大气剖面水分含量大气透射率估计方相关系数标准误
r口w/(g•cm-2)程平方(R2)差(SEE)
t6=0.974290-
0.08007W
t6=1.031412-
0.11536W
041.6
0.99611|0.002368
高气温
1.6-3.0
0.9982710.002539
低气温
041.6
1.6-3.0
t6=0.982007-
0.09611W
t6=1.053710-
0.14142W
0.994630.003340
0.998990.002375
接下来只需要计算出大气水分含量便可计算得到大气透射率。
学
者杨景梅根据地面、高空气象要输资料,提出利用地面湿度计算整层
大气可降水量,由此计算得出整层大气水汽含量的方法(杨景梅邱金恒,用地面湿度参量计算我国整层大气可降水量及有效水汽含量的方法的研究,大气科学,2002.01)公式如下:
&)=cQ+q职
式中,3为大气水汽含量(g/cm2),W为整层大气可降水量(cm);
COC1为经验系数,以重庆为研究对象,一般可选取成都的经验值来
计算。
CO二0.1274;C仁0.6878.计算大气可降水量公式如下:
W-坷总
式中,W为整层大气可降水量(cm),e为地面水汽压(hPa),a0al为经验系数,可从以下公式计算得出。
0_03+2.74J7+0.15)(卩>33°)
0.04+£(^<33°)
=\
0.20—£了<3y)
J0.05
d.;
(p-25.0)J+0.25
fo9>2(y)
乩-\
■-0.9(^<204)
0.066
((p-33)2十4.41
式中,©为地理纬度(°,H为海拔高度。
重庆市主城区纬度取29.5,平均高程取0.35km,则有。
3=0.19604e+0.04691
其中:
e是绝对水汽压(hPa)。
采用WMO在1996年公布的《气象仪器和观测方法指南》中公布的
饱和水汽压E公式。
对于纯水面用公式计算,适用于温度-45~60C。
E=6.112exp[17.62*To/(243.12+To)]
RH=e/E*100
推算出:
e=6.112exp[17.62*To/(243.12+To)]*RH
式中,RH为相对湿度,可以通过气象观测数据获得,To为气温「C)。
1.2地表温度反演步骤
▼
1.2.1影像预处理
(1)辐射定标
选择RadiometricCorrection/RadiometricCalibration。
在FileSelection对话框中,选择数据LC81230322013132LGN02_MTL_Thermal,单击SpectralSubset选择ThermalInfrared1(10.9),打开RadiometricCalibration面板。
(2)大气校正
本教程选择Flaash校正法。
FLAASHAtmosphericCorrection,双
击此工具,打开辐射定标的数据,进行相关的参数设置进行大气校正。
注意:
如果在多光谱数据辐射定标时Scalefactor值忘记设置,可在本步骤中打
开辐射定标数时设置singlescalefaceor值为0.1,若已设置,则默认值为1即
1)InputRadianeeImage打开辐射定标结果数据;
2)设置输出反射率的路径,由于定标时候;
3)设置输出FLAASH校正文件的路径,最优状态:
路径所在磁盘空间足够大;
4)中心点经纬度SceneCenterLocation自动获取;
5)选择传感器类型:
Landsat-8OLI;其对应的传感器高度以及影像数据的分辨率自动读取;
6)设置研究区域的地面高程数据;
7)影像生成时的飞行过境时间:
在layermanager中的Lc8数据图层
右键选择ViewMetadata,浏览time字段获取成像时间;
注:
也可以从元文件“LC81230322013132LGN02_MTL.tXt中找到,具体名称:
DATE_ACQUIRED=2013-05-12;SCENE_CENTER_TIME=02:
55:
26.6336980Z;
8)大气模型参数选择:
Sub-ArcticSummer(根据成像时间和纬度信息选择);
9)气溶胶模型AerosolModel:
Urban,气溶胶反演方法Aerosol
Retrieval:
2-band(K-T);
10)其他参数按照默认设置即可。
11)多光谱参数设置中,
K-T反演选择默认模式:
Defaults->Over-LandRetrievalstandard
(600:
2100)
波谱响应函数:
默认指向..
'ProgramFiles\Exelis\ENVI51\classic\filt_func\landsat8_oli.sli
把它重新指向:
..'Program
Files\Exelis\ENVI51\resource\filterfuncs\landsat8_oli.sli
注:
这是因为ENVI5.1版本的一个小bug,即Classic中的L8的波
谱响应函数不正确,另外一个一劳永逸的方法是:
将
“\Exelis\ENVI51\resource\filterfuncs中的”Landsat8_oli.sli和
Landsat8_oli.hdr两个文件拷贝覆盖:
“...\ENVI51\classic\filt_func中
的两个文件。
否则SWIR1波段大气校正后的结果全为0。
见下图
12)高级参数设置:
根据内存大小设置TileSize(Mb):
100(8g物理
内存),其他参数默认即可,详细见下图
经过FLAASH校正的影像基本去除了空气中水汽颗粒等因子的影响,植被的波谱曲线趋于正常。
1.2.1.1地表比辐射率计算
(1)植被覆盖度计算
第一步,利用多光谱遥感影像,计算研究区域归一化植被指数NDVI,ENVI软件中有工具集成,使用即可。
第二步,利用波段运算工具,根据覃志豪经验公式,计算研究区植被覆盖度。
(b1gt0.7)*1+(b