遥感单窗算法ENVI怎么弄的啊,
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/26 06:29:21
遥感单窗算法ENVI怎么弄的啊,
遥感单窗算法ENVI怎么弄的啊,
遥感单窗算法ENVI怎么弄的啊,
以Landsat TM/ETM影像为例:
一、原理
1.该算法只需要三个参数:大气平均作用温度、大气透射率和地表比辐射率
TM/ETM波段的热辐射传导方程:
B6(T6)=t6(q)[e6B6(Ts)+(1-e6)I6~]+I6_
Ts是地表温度;
T6是TM6的亮度温度;
t6是大气透射率;
e6是地表辐射率.
B6(T6)表示TM6遥感器所接收到的热辐射强度;
B6(Ts)是地表在TM6波段区间内的实际热辐射强度,直接决取于地表温度;
I6~和I6_分别是大气在TM6波段区间内的向上和向下热辐射强度.
化简后最终的单窗体算法模型为:
Ts={a(1-C-D)+[b(1-C+D)+C+D]T6-DTa}/C
式中
C6=t6e6(e6为比辐射率,t6为透射率)
D6=(1-t6)[1+t6(1-e6)]
a =-67.355351,b=0.458606
2.地表比辐射率的估计(地表划分为水体、城镇和自然表面三种地物类型,详细可见http://blog.sina.com.cn/s/blog_764b1e9d0100r484.html)
其中大气平均作用温度Ta的近似估计:
温度换算:T=t+273.15
对于中纬度夏季平均大气Ta=16.0110+0.92621T0,冬季采用Ta=19.2704+0.91118T0
大气透射率t6的估计:
3.关于像元亮度温度计算:
T6=K2/ln(1+K1/L6)
其中:K1=607.76,K2=1260.56(覃志豪,用陆地卫星TM6数据演算地表温度的单窗算法);
L6为遥感器接收的辐射强度.
4.遥感器接收的辐射强度计算
L6=(15.303-1.238)*b1/255.0+1.238 (覃志豪,用陆地卫星TM6数据演算地表温度的单窗算法); b1为第六波段像元灰度值(DN值).
LMAX_BAND6 = 15.303
LMIN_BAND6 = 1.238
二、软件操作:
计算L6
ENVI->Basic Tools->Band Math
在Enter anexpression中输入(15.303-1.238)*b1/255.0+1.238,点击Add to List.
选中刚才输入的公式,点击OK.
选中B1变量,在Available Bands List中选Band6,保存文件名为L6,点击OK.
2.计算T6(同计算L6,只是计算公式不同).
3.计算C6(同计算L6,只是计算公式不同).
4.计算D6(同计算L6,只是计算公式不同).
5.计算Ts
同计算C6相同,在Enter an expression中输入:
(((-67.355351*(1-b1-b2)+(0.458606*(1-b1-b2)+b1+b2)*b3-b2*292.15753)/b1-273.15)),点击Add to List.
选中b1变量,在Available Bands List中选择C6,
选中b2变量,在Available Bands List中选择D6,
选中b3变量,在Available Bands List中选择T6,保存文件名为Ts ,点击OK.