地表反射率温度植被指数文档格式.docx
《地表反射率温度植被指数文档格式.docx》由会员分享,可在线阅读,更多相关《地表反射率温度植被指数文档格式.docx(29页珍藏版)》请在冰豆网上搜索。
在有植被覆盖的情况下,NDVI介于0.1~1.0之间,且随植被覆盖度的增大而增加。
几种典型的地面覆盖类型在大尺度NDVI图像上区分鲜明,植被得到有效的突出。
因此,特别适用于全球或区域植被的宏观动态监测。
3、NDVI是植被生长状态及植被覆盖度的最佳指示因子。
NDVI与叶面积指数(LAI)、绿色植物净第一生产率(NPP)、植被覆盖度、光合作用等指标参数密切相关,并据以完成地表植被物理参数的反演。
NDVI与光合有效吸收辐射(FAPAR)近线性关系;
NDVI与LAI呈非线性相关关系;
NDVI的时间变化曲线可反映季节和人为活动的变化;
NDVI在生长季节内的时间积分与净第一生产力(NPP)相关。
NDVI可以部分消除与太阳高度角、卫星观测角、地形、云/阴影和大气条件有关的辐照度条件变化(大气程辐射)等的影响。
4、地表反射率的计算:
根据遥感影像DN值计算到传感器的各波段辐射量度:
Lλ=“gain”*QCAL+“bias”;
各波段表观反射率的计算:
ρ=Pi*L*d*d/(ESUNλ*cos(θs)),其中ρ=表观反射率或是行星反射率,
Lλ=传感器获得的光谱波段辐射量度,d=日地距离,ESUNλ=大气层顶平均入射太阳辐射,θs=太阳天顶角。
逐波段大气辐射校正并计算表观反射率与地表反射率的线性关系:
Y=aX+b。
5、地表温度的计算(反演):
陆面温度(LST):
根据辐射温度理论,利用热红外波段的遥感数据反演计算获得的地表亮度温度,又成表皮温度;
基于landsat7ETM+数据的陆面温度反演算法T=k2/ln(k1/Lλ+1),其中:
T=卫星接收有效温度,K1,k2=定标系数,L=光谱辐射亮度。
数据准备
1、
用于NDVI计算模式的hrf的img多光谱影像
2、
用于地表温度反演的htm的img第6波段的影像
操作方法及过程
一、
DNVI建模
辐射亮度计算。
根据hrf头文件中的数据,可得到辐射亮度计算亮度计算公式为0.94251966100978102*b+(-5.9425196610097810)。
具体步骤:
BasicTools|BandMath,在BandMath对话框中输入公式,分别计算第3和第4波段的辐射亮度,并分别导出数据b3_1和b4_1。
行星反射率计算。
同样通过hrf头文件中的数据,可得到行星反射率计算公式为!
pi*(0.94251966100978102*b+(-5.9425196610097810))*1.0136^2/(1551*cos(64.799999999999997/180*!
pi))。
BasicTools|BandMath,在BandMath对话框中输入公式,公式中的b分别选取第3和第4波段的辐射亮度,导出结果b3_2和b4_2。
3、
利用5S大气模型计算的地表反射率与行星发射率的拟合函数关系求解。
运行5S拟合模型程序,根据hrf影像的信息,设置拟合模型的各个参数如下,SiteElevation(平均高程)设为2300;
SenorAltitude(传感器高度)设为705000;
MonthandDate(年月)设为May30;
SolarZenithAngle(太阳天顶角)设为25.2;
SolarAzmethAngle<
0-360>
Degree(太阳方位角):
122.2;
standerdAtosphericModels设置为MidlatitudeSummer;
AerosolModels设置为Continental;
其他的默认。
然后在HomogenousReflectanceforband1中设置其参数值为0(该波段行星反射率的值),得出对应的地表反射率—将显示的结果并输入到excel中—每隔0.05显示一个结果—将得到的值统计出来,并将统计结果绘制成散点图—再将其拟合成对应的线性关系。
4、
地表反射率计算。
根据5S模型拟合出来的3、4波段的数学关系式,其关系式分别为y=0.9161x+0.0143和y=0.9401x+0.0065,然后计算出3、4波段的地表反射率。
步骤:
BasicTools|BandMath,在BandMath对话框中输入公式,公式中的x分别选取第3和第4波段的行星反射率,计算第3和第4波段的地表反射率,导出数据b3_3和b4_3。
5、
NDVI值的计算。
根据计算出来的3、4波段的地表反射率,以及DNVI的公式的到NDVI值。
BasicTools|BandMath,在BandMath对话框中输入公式,公式中的b3和b4分别选取第3和第4波段的地表反射率。
然后导出结果。
二、
地表温度反演
1、计算辐射亮度。
加载htm影像,根据头文件中的数据,得到1、2波段的辐射亮度的计算公式0.7667001*b1+(-0.7667001)和0.9868001*b2+(3.1627953249638470),步骤同上,得出辐射量度的计算结果。
2、辐射反演。
利用公式T=k2/ln(k1/Lλ+1)算地物的辐射反演,其中T为开尔温度;
查找参数值:
k1=666.09;
k2=1282.71;
Lλ分别利用步骤1中的波段1和波段2的辐射量度。
3、统计反演后的地物的温度值,并比较其差异。
打开反演后的温度影像,右击影像选择ROITool,统计各种地物值的最大值,最小值,均值,标准差,将其统计到Excel中,比较其差异。
结果与分析
一、DNVI建模
【地表反射率】
第3波段
第4波段
【DNVI】
【3、4波段表观反射率和地表反射率的线性关系】
【表观反射率和地表反射率的线性关系数学表达式】
波段
关系式
1波段
y=0.8933*x+0.0473
4波段
y=0.9401*x+0.0065
2波段
y=0.8801*x+0.0242
5波段
y=0.9399*x+0.001
3波段
y=0.9161*x+0.0143
7波段
y=0.9584*x+0.0004
【部分地物的DNVI值】
地物 DNVI值
min
max
mean
stdev
Reservior
0.057713
0.338587
0.145087
0.038598
Snow
-0.12395
0.152669
0.025088
0.031572
BareLand
0.105628
0.374843
0.192701
0.043621
Urban
-0.356923
0.038094
-0.273288
0.045284
Plant
0.333387
0.786695
0.656094
0.081619
Desert
0.071897
0.155663
0.100783
0.014291
River
0.043469
0.429917
0.127503
0.08131
【结果与分析】:
通过对提取地物的DNVI值的可以发现,绿色的DNVI值比较高,原因是绿色植物叶绿素引起的红光吸收和叶肉组织引起的近红外光反射使得植被在近红外波段和红光波段有很大的差异;
水体和裸地在红光波段和近红外波段反射率相当,因此水库和裸地的NDVI值接近0;
雪地NDVI最低值中出现负值,是由于在近红外波段比可见光波段有较低的反射率;
沙漠中植被很少,因此其近红外波段和红光波段的反射情况和裸地类似,因此其NDVI值接近于0;
河流的NDVI值稍大于由于河流中存在一定的含沙量,使得地物在近红外波段的反射率大于近红外波段。
二、温度反演
【温度反演】
低增益温度反演
高增益温度反演
【第1波段部分地物低增益温度反演数据】
开尔文温度
摄氏温度
反演温度
地物
287.47641
289.289886
288.13127
0.388036
14.32641
16.13989
14.98127
273.154785
293.990417
278.177771
3.788266
0.004785
20.84042
5.027771
295.989319
310.676086
303.445647
2.819391
22.83932
37.52609
30.29565
300.165253
310.928528
307.469228
1.530421
27.01525
37.77853
34.31923
294.278015
305.525879
298.698402
2.333251
21.12802
32.37588
25.5484
302.605286
309.915955
306.491575
1.39902
29.45529
36.76596
33.34158
300.438721
313.922485
305.865796
4.320705
27.28872
40.77249
32.7158
m