基于ENVI软件的遥感信息处理

《基于ENVI软件的遥感信息处理》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、基于ENVI软件的遥感信息处理实验报告一、遥感数据预处理：辐射亮度与表观反射率的计算1.实验目的将传感器记录的数字量化值（DN）转换成绝对辐射亮度值（辐射率）与地表（表观）反射率等物理量有关的相对值的处理过程。2.实验原理lDN值转为辐射亮度公式：Lrad=(DN-1)*UnitConversionCoefficient。其中：Lrad-----辐射亮度；DN-----遥感影像像元亮度值，记录地物的灰度值（无单位，是一个整数值，值大小与传感器的辐射分辨率、地物发射率、大气透过率和散射率等有关）；UnitConversionCoefficient------定标系数（表1-1

2、所示）；Aster数据定标系数（表1-1）AsterbandCoefficent(W/m2*sr*um)/DNHighGainNormalLowGain1LowGain210.6761.6882.2520.7081.4151.893N0.4230.8621.153B0.4230.8621.1540.10870.21740.29000.290050.03480.06960.09250.409060.03130.06250.08300.390070.02990.05970.07950.332080.02090.04170.05560.245090.101590.03180.04

3、240.2650l辐射亮度转为表观反射率公式：RTOA=(pi*Lrad*d2)/(ESUNi*COS(z))d=(1-0.01672*COS(RADLANS(0.9856*(JulianDay-4))))其中：RTOA-----表观反射率；Pi-----3.14159；d-----日地距离；z-----太阳天顶角(z=90-太阳高度角)；JulianDay-----成像日期；ESUNi-----第i个波段的增益（表1-2所示）；Aster数据各个波段增益（表1-2）AsterbandSmith：ESUNiThomeetal(A):ESUNiThomeetal(A):ESU

4、NiB11845.9918471848B21555.7415531549B3N1119.4711181114B4231.25232.5225.4B579.8180.3286.63B674.9974.9281.85B768.6669.2074.85B859.7459.8266.49B956.9257.3259.853.数据来源本次实验所使用数据为2001年宜昌地区遥感影像，该数据为Aster1-9波段，30米分辨率，投影为WGS-84，UL纬度为30.82418333，UL经度为111.21726944。4.实现过程（1）打开并显示影像文件，选择合适的波段组合加载影像打开并显

5、示aster影像文件，从envi主菜单中，选择File——OpenImageFile选择yichang2001aster_打开，选择band3，band2，band1进行组合加载显示。（2）使用BandMath对原始数据各个波段进行计算，实现DN值转为辐射亮度envi主菜单选择BasicTools——BandMath——Enteranexpression，输入表达式，即将各个波段的DN值转为辐射亮度的公式依次输入并作用于各个相应的波段，选择File——SaveFileAs——ENVIStandard（注意调整波段顺序），输出数据即为将DN值转换为辐射亮度后的影像。（3）使用

6、BandMath对原始数据各个波段进行计算，实现DN值转为表观反射率envi主菜单选择BasicTools——BandMath——Enteranexpression，输入表达式，即将各个波段的DN值转为表观反射率的公式依次输入并作用于各个相应的波段，选择File——SaveFileAs——ENVIStandard（注意调整波段顺序），输出数据即为将DN值转换为表观反射率后的影像。二、地图投影的转换1.实验目的把UTM投影影像转换成高斯－克吕格投影和阿尔勃斯投影的影像；任选五个坐标点对，做高斯—克吕格投影的正反算，并与ENVI软件比较计算精度。2.实验原理Albers投影，又

7、名“正轴等积割圆锥投影”，圆锥投影的一种。它是假想一个圆锥面与地球椭球旋转轴重合的套在椭球上，按等积条件把地球椭球上的经纬线投影到圆锥面上，然后沿一条母线将圆锥面切开展成平面。Albers投影的参数：参考椭球体KRASOVSKY，X轴偏量0.0000000000，Y轴偏量0.0000000000，第一标准纬线：25°0′0″N，第二标准纬线：47°0′0″N，中央经线：105°0′0″E，00.000/*起始经线，坐标原点：0°0′0″，半球参数：Krasovskya=6378245.0000b=6356863.0