用NDSI归一化雪被指数提取被冰雪覆盖地区
以L5134040_04020051111为数据源，进行冰雪覆盖地区的提取。
在Landsat5中NDSI=(B2-B5)/(B2+B5)
此处用第2、5波段的反射率进行计算，NDSI=(Ρtm2-Ρtm5)/( Ρtm2+Ρtm5)
NDSI>=0.4的区域即为冰雪覆盖区。
1. 把DN值图像转化成亮度值图像
L=（Lmax-Lmin)/(Qmax-Qmin)*(Qdn-Qmin)+Lmin
Lmax,Lmin,Qmax,Qmin分别为最大最小亮度值，最大最小DN值，可从头文件中获取。
2. 把亮度值图像转化成反射率图像
ρ=π*Lλ*d^2/(ESUNλ*COSθ)
d日地距离可由DOY 求得，d=1-0.01674*cos(0.9856*(DOY-4)* π/180)
DOY：Day of year
θ太阳天顶角，90-太阳高度角
3. 在生成的反射率图像上计算NDSI
用BAND MATH 计算NDSI=(float(b2)-float(b5)) /(float(b2)+float(b5))
4. 把NDSI上大于等于0.4的地区提取出来
BAND MATH中输入(b1 ge 0.4)*1+(b1 lt 04)*0,把大于等于0.4的区域赋值1，其他区域赋值0，即可把冰雪覆盖地区提取出来。
左图为反射率图像，给7、4、1波段分别赋予RGB颜色，右边为提取结果白色区域为冰雪覆盖区，可将两图连接查看成果。
5. 将DEM数据加载进去，生成3D影像，更直观的进行观测。
本文来源：https://www.2haoxitong.net/k/doc/735b763c376baf1ffc4fad99.html
文档为doc格式