问题描述：
遥感单窗算法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.4586062.地表比辐射率的估计（地表划分为水体、城镇和自然表面三种地物类型,详细可见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.303LMIN_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.
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.4586062.地表比辐射率的估计（地表划分为水体、城镇和自然表面三种地物类型,详细可见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.303LMIN_BAND6 = 1.238二、软件操作：计算L6ENVI->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.
展开全文阅读