关于MODIS卫星数据的几何校正方法.docVIP

关于MODIS卫星数据的几何校正方法 遥感信息匝用技 关于MODIS卫星数据的几何校正方法 郭广猛 (中科院地理科学与资源研究所北京,100101) 摘要:MODIS(MODerate--resolutionImagingSpectroradiometer)LIB产品中存在着bowtieeffect,即”双眼皮”现象,几 何校正之后仍无法去除,因此必须在几何校正之前加以去除.本文使用builtGLT和exportGCP两种方法对MODIS数据进行 几何校正,结果表明两种方法的效果都比较准确.但是前一种方法运算耗时太多,因此建议采用后一种方法校正. 关键词:MODIS;几何校正;Bowtieeffect 中图分类号;TP751文献标识:A文章编号:1OOO一3177(2002)一67—026—03 1999年12月美国的TERRA卫星发射成功,星 上搭载的中分辨率成像光谱仪MODIS(MODerate — resolutionImagingSpectroradiometer)是该星上 唯一实行全球免费接收的传感器.MODIs具有36 个波段(0.405～14.385um),拥有250m,500m, 1000m三个不同的空间分辨率,扫描宽度 2330km,只需2～3轨即可覆盖全国,相对于现在常 用的NOAA—AVHRR数据(5个波段,分辨率 1.1km)有了极大改进,这些数据将在地球科学的综合 研究和海洋,陆地,大气的专门研究中发挥重要作用. 1去除处理”双眼皮”现象(bowtieeffect) MODIS数据经辐射校正之后生成的L1B产品 中存在着独有的”双眼皮”现象(bowtieeffect),使得 MODIs的边缘数据无法使用,影响了数据的实际 应用,因几何校正亦无法去除”双眼皮”现象,因此必 须在几何校正之前就加以去除.”双眼皮”现象表现 为相邻两个扫描行之间有部分数据相同,越向边缘 重复数据越多,在线状地物附近表现尤为明显,见图 1所示,河流不仅左右错开,而且上下有部分重叠. 图2是bowtieeffect示意图l】],图中显示有3个扫描 行,每行有10个像元,每个像元在星下点处分辨率 为lkm,所以每个扫描行宽度为10km.由于每扫描 一 次地球会自转一定距离,因此图上显示3个扫描 行之间有一个小的错动.MODIS扫描角度为±55., 对应地面宽度为2330km.由于地球曲率的影响和 扫描角度的增大,越向边缘像元尺寸越大,在最边缘 处一个扫描行的宽度已增大到20km,相邻的扫描行 之间已有10km的重叠,于是就产生了”双眼皮”现 象.对于lkm分辨率的MODIS数据而言,在角度大 约25度时开始有重叠,分辨率250m的数据在17度 时就有了重叠. 图1MODIS250m数据的bowtieeffect ‘11B5km 图2MODISlkm数据左半部分bowtieeffect示意图 (BowtieEffectSketchMapofLeftHalfMODISlkmData) MODIS数据采用HDF格式(HierarchicalData Format),是美国国家高级计算应用中心研制的新 型数据格式.它将多种类型的数据,例如图像数 收稿日期:2.O2O5 作者简介:郭广猛,(1976～),男,1998年毕业于青岛海洋大学,2001年毕业于中国地质大学(武汉),现就读于中科院地理与资 源研究所,主要从理],极大地方便了普通用户.本文即利用IDL中 的函数对HDF文件进行处理以去除”双眼皮”现象. 下面以250m的MODIS数据为例说明数据处 理思路.250m的数据每扫描一次产生5416列,40 行数据,每相邻的两个40行之间存在有”双眼皮”现 象.一幅250m的MODIS数据有5416列,4080行. 首先计算相邻的两个40行之间存在多少行重复数 据.以第1列为例: (1)第1列到第3列,第38行到第40行数据产 生一个3*3矩阵A,以第1列到第3列,第41行到第 43行数据产生矩阵B,计算A和B的相关系数. (2)将矩阵B下移一行,重新计算相关系数. (3)重复步骤2,直到第80行,这样就得到40个 相关系数. (4)取4O个相关系数中最大的一个,它所对应 的行数减去40,再加上3,即为重复的行数.最后算出 对于第1～3列数据每两个扫描行有20行是重复的. (5)去除每个扫描行的20行重复数据,这样即 可去除”双眼皮”现象.但是原来的3列,4080行变 成了3列,2040行,这时需重采样将其变为原来大 小. (6)将矩阵A向图像中心移动,计算每列上的 重复,结果应该是逐渐减小,利用这一规律可以对某 些异常值进行取舍.图3是与图1对应的去除”双眼 皮”现象之后的结果,河流边缘平滑,且周围地物没 有重叠现象. 图3去除bowti