第6章辐射校正扩展阅读.pdfVIP

辐射校正补充案例 （配第6 章）  测试数据 数据目录：附带光盘下的..\chapter06\data\ 文件名 说明 Landsat5_sd_ref2.dat FLAASH 大气校正后的图像 Landsat5_sd_dem.dat DEM 图像  IDL 程序  程序目录：附带光盘下的..\chapter06\programs\ 文件名 说明 Topo_correction.sav C 地形校正程序  案例背景 在丘陵地带和山区，地形坡度、坡向和太阳光照几何条件等对遥感图像的辐 射亮度影响是非常显著的，朝向太阳的坡面会接受到更多的光照，在遥感图像上 色彩自然要亮一些，背向太阳的阴面由于反射的是天空散射光，在图像上表现得 要暗淡一些。地形校正的目的是消除由地形引起的辐射亮度误差，使坡度不同但 反射性质相同的地物在图像中具有相同的亮度值。ENVI 中没有提供地形校正工 具，但是可以根据半经验C 校正方法采用IDL 编程实现。 半经验C 校正的公式为： cosθ+c L L （6.1） H T cos i +c 式中，LH 指水平面上像元的等效观测值，即校正后的值；LT 指倾斜面上像元的 观测值，即原始像元值；θ指太阳天顶角；i 指太阳入射角；c 为半经验系数。太 阳入射角i 的求解如 ［式（6.2 ）］所示： cos i cosθcosα=+sinθsinαcos(β −φ) （6.2） 式中，α表示坡度；β 表示坡向；ϕ表示太阳方位角。c 没有确切的物理意义，它 的计算是首先利用像元亮度值和太阳入射角余弦值之间的统计回归求取斜率 k 和截距b ，然后取b 和k 的比值，即： 1 c b / k （6.3 ） 所以，地形校正需要输入的数据有遥感图像和相应的DEM 数据，输入的参 数有太阳天顶角和太阳方位角。本案例的测试数据为Landsat5_sd_ref2.dat ，相应 的DEM 数据为Landsat5_sd_dem.dat 。所用的IDL 程序为Topo_correction.sav ， 该程序的原作者为Morton Canty （详见参考文献），该程序利用半经验C 校正方 法对遥感数据进行地形校正。 6.1 添加扩展工具 将Topo_correction.sav 复制到ENVI 安装目录的Exelis\ENVI52\extensions 路 径下。 6.2 启动ENVI 单击开始 所有程序 ENVI 5.2 ENVI ，启动ENVI 5.2.1 。此时，工具箱 中的Extensions 节点下将会自动添加Topo_correction 工具（图6.1 ）。 2 图6.1 Topo_correction 工具 6.3 加载遥感数据 在ENVI 主界面菜单栏中，选择File Open 打开Open 对话框，添加遥感数 据 Landsat5_sd_ref2.dat 和 DEM 数 据 Landsat_sd_dem.dat 。 选 择 Landsat5_sd_ref2.dat 的第5、4 、3 波段进行RGB 组合并进行Linear 2%拉伸显示 （图6.2 ）。