基于ENVI软件的遥感信息处理案例.doc

基于ENVI软件的遥感信息处理 实验报告 一、 遥感数据预处理：辐射亮度与表观反射率的计算 1.实验目的 将传感器记录的数字量化值（DN）转换成绝对辐射亮度值（辐射率）与地表（表观）反射率等物理量有关的相对值的处理过程。 2.实验原理 DN值转为辐射亮度 公式：Lrad =(DN-1)*Unit Conversion Coefficient。 其中：Lrad-----辐射亮度； DN-----遥感影像像元亮度值，记录地物的灰度值（无单位，是一个整数值，值大小与传感器的辐射分辨率、地物发射率、大气透过率和散射率等有关）； Unit Conversion Coefficient------定标系数（表1-1所示）； Aster数据定标系数（表1-1） Aster band Coefficent(W/m2*sr*um)/DN High Gain Normal Low Gain 1 Low Gain 2 1 0.676 1.688 2.25 2 0.708 1.415 1.89 3N 0.423 0.862 1.15 3B 0.423 0.862 1.15 4 0.1087 0.2174 0.2900 0.2900 5 0.0348 0.0696 0.0925 0.4090 6 0.0313 0.0625 0.0830 0.3900 7 0.0299 0.0597 0.0795 0.3320 8 0.0209 0.0417 0.0556 0.2450 9 0.10159 0.0318 0.0424 0.2650 辐射亮度转为表观反射率 公式：RTOA=(pi*Lrad*d2)/(ESUNi*COS(z)) d=(1-0.01672*COS(RADLANS(0.9856*(Julian Day-4)))) 其中：RTOA-----表观反射率； Pi-----3.14159； d-----日地距离； z-----太阳天顶角(z=90-太阳高度角)； Julian Day-----成像日期； ESUNi-----第i个波段的增益（表1-2所示）； Aster数据各个波段增益（表1-2） Aster band Smith： Thome et al(A):ESUNi Thome et al(A):ESUNi B1 1845.99 1847 1848 B2 1555.74 1553 1549 B3N 1119.47 1118 1114 B4 231.25 232.5 225.4 B5 79.81 80.32 86.63 B6 74.99 74.92 81.85 B7 68.66 69.20 74.85 B8 59.74 59.82 66.49 B9 56.92 57.32 59.85 3.数据来源 本次实验所使用数据为2001年宜昌地区遥感影像，该数据为Aster1-9波段，30米分辨率，投影为WGS-84，UL纬度为30UL经度为111 4.实现过程 （1）打开并显示影像文件，选择合适的波段组合加载影像 打开并显示aster影像文件，从envi主菜单中，选择File——Open Image File选择yichang2001aster_打开，选择band3，band2，band1进行组合加载显示。 （2）使用Band Math对原始数据各个波段进行计算，实现DN值转为辐射亮度 envi主菜单选择Basic Tools——Band Math——Enter an expression，输入表达式，即将各个波段的DN值转为辐射亮度的公式依次输入并作用于各个相应的波段，选择File——Save File As——ENVI Stand