1、 裁剪出出济南市区
2、 分别利用ENVI 、ERDAS 反演地表温度(LST )、NDVI ,对LST 进行彩色显示。
3、 分析LST 、NDVI 的关系。
反演公式
具体流程:
图像的DN 值 辐射亮度 辐射亮温 地表温度。
反演时从图像数值(DN )转换成绝对辐射亮度值时的公式、从辐射亮度值转成辐射亮温时的公式、从亮温转换成地表温度时的公式分别是:
min min max 6255)(L L L DN L tm +-⨯=、 )1/ln(/12+=λL K K T 、 ε
ρλl n )/(1T T T s += 其中:6tm L 为TM 传感器所接收到的辐射亮度(mW .cm -2s r-1.um -1),max L 、min L 分别是传感器所接收到的最大和最小的辐射强度,即对应于DN =255和DN =0时的最大和最小辐射强度。
对于Landsat5的TM 6波段,
1K =60.77mW .cm -2s r-1.um -1,2K =1260.56K 。
S T 为地表温度(K )
;T 为辐射温度(K );λ为有效波谱范围内的最大灵敏度值,λ=11.5um ,ρ=/hc δ=1.438×10-2mk ,其中δ=1.38 ×10-23/J k ,为玻尔兹曼常数,h =6.626×10-34Js ,为Plank′s 常数,c =2.998 ×108/m s ,为光速。
一般地,有植被覆盖的地表取ε=0.95,没有植被覆盖的地表取ε=0.92(Weng ,2004[16])。
min L =0.1238
255
)(min max L L - =0.005632156 )1/ln(/12+=λL K K T 1260.56 / LOG ( 1 + 60.766 / $n8_fu )
$n1_12736l / (1 + (0.0000115 * $n1_12736l /0.01438) * LOG (0.95 ) )。