2-MODIS影像的几何处理算法研究.doc

MODIS影像的几何处理算法研究钱海明作者简介：钱海明(1989-)，男，硕士，研究方向：摄影测量与遥感、遥感图像处理与分析。Email：1070572202，王春林，孙金彦（安徽省（水利部淮河委员会）水利科学研究院,合肥 230088）摘要：针对MODIS影像边缘区域存在明显的重叠状几何畸变数据，及投影转换耗时长、算法效率低等问题，研究了MODIS L1B的数据存储结构，利用自带星下点1公里间隔维度数据，解算扫描方向上的影像重叠度。通过重采样方式消除蝴蝶效应，改正异常值，提高MODIS影像数据的分析及应用水平。实验表明，在计算速度和影像质量上，本文方法具有良好的适应性。关键词：MODIS；蝴蝶结效应；畸变校正；MODIS（中分辨率成像光谱仪）是Terra和Aqua卫星上搭载的主要传感器，共有36个观测通道，波谱范围从0.4m到14m，在环境监测、气象预报等方面有着重要的研究与应用价值。研制MODIS卫星系统旨在建立全球大气、海洋和陆地的动力学模型，并能在地球发生变化之前做出准确的预测。由于地球曲率影响和大扫描角等问题，导致MODIS影像的扫描带之间出现重叠错位现象，简称“蝴蝶结效应”。针对MODIS影像数据的蝴蝶效应问题，黄家洁等阐述了MODIS的成像规律、几何畸变和纠正处理的基本方法，但是没有提出有效的纠正措施；郭广猛等利用相关系数法校正“蝴蝶结效应”；陶卫中等利用MOD03地理数据结合分形插值方法去除蝴蝶结现象；隋明等4利用矢量数据校正MODIS数据。如何进一步使用L1B产品，尚未得到统一且处理得非常彻底的方法。因此，研究消除MODIS数据中的蝴蝶效应的工作显得非常迫切。本文研究了MODIS L1B的数据存储结构，利用自带星下点1公里间隔维度数据，解算扫描方向上的影像重叠度。通过重采样消除蝴蝶效应，改正异常值，提高MODIS影像数据的分析及应用水平。 1蝴蝶结效应MODIS 是一种被动式摆动扫描探测器，摆动扫描角为 55。由于地球曲率以及探测方式的影响，像素的大小随扫描角的增大而增大，扫描带宽实际达到了 2 330 km。每采集一次数据，MODIS 探测沿轨道中心前进10 km，扫描距离20 km。两个相邻扫描带在边缘处存在50%的数据为重复数据，此现象被称为“蝴蝶结效应”。如图1所示。图1 蝴蝶结效应2数据来源以2015年8月的EOS/MODI遥感影像为数据源。所采用的遥感影像是每日实时接收并预处理的中等分辨率成像光谱仪（moderate-resolution imaging spectroradiometer, MOD13Q1）的陆地II级标准数据产品（Terra/MODIS及aqua/MODIS数据），两个波段，空间分辨率为250m。TERRA及AQUA每日过境时间分别为上午10:30和下午1:30。3 “蝴蝶结效应”纠正方法3.1 图像纠正针对空间分辨率为250m的MODIS数据，首先将空间分辨率为1KM的纬度数据重采样到250米，使得地理数据与影像数据保持一致；然后，记录每一列纬度数据的重叠度，在250米影像数据中，每4列影像数据对应一个重叠纬度数据；最后，根据重叠度对影像进行重采样，消除“蝴蝶结效应”。以扫描纬度数据并记录其重叠程度，然后用记录对图像数据去除蝴蝶现象，算法示意图如图2所示。图2 第一列纬度数据蝴蝶结效应去除前后对比3.2 坏点修复由于仪器异常、控制程序出错等原因导致数据存在异常点，即在纬度数据中记录为-999的数据便是数据异常（正常值范围为-90.0+90.0），称为数据坏点。在纬度数据中可能会出现整个扫描行均为数据坏点，这对数据使用带来很大影响。针对以上问题，首先扫描全部纬度数据值，判断是否存在数据坏点；如果存在则进一步判断该列下一个像元纬度数据值是否也同样为数据坏点，依次判断直到检测到不是数据坏点的像元纬度数据值。记录下最初的数据坏点所在列的上一个正常值，同时也记录下检测到的最后一个数据坏点所在列的下一个正常值，而后进行线性插值处理填充数据坏点。采用类似方法处理图像数据，但是图像数据的存贮按列并不是线性分布的，因此对于数据坏点长度较大的列就无法直接线性插值来进行填充。考虑到实际情况，算法最多对宽度达到一个扫描行的数据坏点进行了处理，而对超出这个范围的数据坏点则不进行处理。从图3中的A图与B图的对比中可以看出当异常数据（A中白色像元）的宽度不超过一条扫描行时，算法能很好地修复。由于示例中数据坏点的宽度并不太大，因此进行线性填充时与周围像元仍能较好地融为一体。同时，由于去除坏点的检测是按列进行的，即使出现图3（A）中部那些连续的横向数据坏点也能较好地修复。图3 数据坏点去除前后对比 (图A为坏点去除前，图B为坏点去除后)4实验与分析为了验证蝴蝶结现象的去除效果，在第一波段(Band 1)沿直线采集像素的亮度值，做成剖面图.，如图4所示。以行号作为横坐标，在蝴蝶结现象去除后的曲线不存在与去除前的曲线重复的曲线段。在每一个扫描行（250m分辨率图像为每40个像素内），曲线保留部分（即未重叠部分）的特征得到了很好的保留，并相应地对需要保留的值进行了拉伸。同时，算法去掉了重复部分后，其基本走势仍与原曲线保持一致。图4 蝴蝶结效应去除后数据剖面图(Band 1)图5 去除蝴蝶结效应图（A图为去除前，B图为去除后）从图5可以看出，算法已成功去除MODIS数据中存在的蝴蝶结效应现象。图5（A）的第1波段左上角存在非常明显的蝴蝶结现象，云层重叠非常严重。通过本文算法处理后，图5（B）中云层已经恢复了连续的形状，重叠现象已经基本被去除了。利用直方图对去除蝴蝶结现象的前后进行亮度对比分析，实验结果如6所示，去除前与去除后的直方图曲线在在总体趋势上吻合的比较好，主要特点部位的吻合的相似度非常高。表明了算法具有较好的整体修饰效果。图6 蝴蝶结效应去除前后第一列数据直方图对比为了提高算法的计算效率，通过不同配置的计算机验证算法的可行性，实验数据为250米的MODIS L1B 数据，实验结果如表2所示，对于最大的250m分辨率影像数据，本文算法对计算机的配置问题具有很好的适应性。表2 250m分辨率MODIS 1B影像全图单波段(Band 1)所耗时间对比CPU内存耗时总耗时奔腾IV3.04G1G12秒30秒赛扬1.5G256M17秒42秒相比于利用图像相似性去除蝴蝶结效应，本文算法具有很好的快速响应能力，能较好的满足MODIS地面站对实时性的要求。5 结论本文利用纬度数据计算影像每列的重叠程度，减少了算法的复杂度，提高了计算效率。通过对纬度数据以及影像数据进行快速异常值检查，快速修复影像中存在的异常点。实验表明，本文算法在去除蝴蝶结现象问题上具有很强的适应性。参考文献：1黄家洁,万幼川,刘良明.MODIS的特性及其应用J.地理空间信息. 2003(04)2郭广猛.非星历表法去除MODIS图像边缘重叠影响的研究J. 遥感技术与应用,2003,(03):172-175.3陶卫中. MODIS数据预处理的并行算法设计D.华中科技大学,2007.4隋明,李刚,李国春,蔡冰,全美兰,王阳. 基于数据矢量特性的MODIS数据几何校正J. 气象与环境学报,2013,(04):94-99.5简兴,张振国,邱银国. 基于ENVI的MODIS数据预处理实现J. 牡丹江师范学院学报(自然科学版),2014,(02):15-17.6刘瑜.基于FLAASH模型的MODIS1B数