摘要
基于DVBS(Digital Video Broadcast by Satellite)系统转发的MODIS 1B影像,依托ENVI/IDL软件,集成了一套MODIS 1B数据的预处理技术流程,提出了合成晴空影像的新思路。详细论述了条带去除、几何校正、大气校正、标准晴空影像合成4个技术环节,为利用MODIS影像开展景观尺度土地覆被变化研究奠定了技术支撑体系。
研究表明,MODIS 1B影像存在横向和纵向条带,横向条带可分为单行条带、宽条带和多行条带3种类型。尽管包括:频率域去噪法、空间域去噪法和统计量调整法等多种去除条带方法,但由于MODIS影像有36个波段,单独使用某种方法,都不能得到较为理想的条带去除效果。基于均衡化曲线补偿技术,提出了一种快速检测和去除条带的算法,给出了条带检测判定参数,即:如果A帧和B帧的10组传感器校正系数Ci值分别位于0的两侧,且ABS(Ci)>0.008,则存在宽条带;如果10组传感器的校正系数MAX(Ci)>0.05时,则存在单行条带。去除宽条带时,需以当前列为首列,从其后连续的100列影像中提取各传感器相对当前列的校正系数。试验表明,条带检测模型能有效识别单行条带和宽条带,多行条带需目视识别,单行条带和宽条带去除效果明显好于NASA网站公布数据,多行条带虽未彻底去除,但处理后的影像质量得到明显改善。
基于ENVI/IDL软件,总结了一套针对HDF格式的MODIS 1B影像的几何校正和大气校正的技术处理方案。采用ENVI中的Georeference MODIS进行去“双眼皮”(bow-tie)和几何校正;采用FLAASH模块对几何校正的短波辐射亮度波段进行大气校正,用热红外大气校正模块(Thermal Atm Correction)对几何校正的长波发射率波段进行大气校正。预处理后的影像质量明显提高,为MODIS产品的反演计算创造了有利条件。研究表明,影像预处理流程应根据不同产品的反演算法来定。一般来讲,几何校正后的数据可用于反演大气产品,大气校正后的影像可用于开发陆地和海洋产品。在反演陆地和海洋产品时,如果指标或算法具有减小大气影响的特性,也可直接使用几何校正后的影像数据作为这些模型和算法的输入数据。
将一七波段比值法合成的影像与NASA公布的8日合成的MOD09GA数据比较,得出一七波段比值法适用于提取晴空陆地像元的结论;将最大饱和度法合成的影像与NASA公布的8日合成的MOD09GA数据比较,得出最大饱和度合成法适用于水体像元的结论。结合一七波段比值法和最大饱和度法,将陆地表面B1/B7最小的像元定为晴空像元;将水体表面饱和度最大的像元定为晴空像元,提出适合华北地区的晴空影像合成方法。基于地表反射率影像合成晴空影像的过程中,逐像元建立了标准晴空的时间索引,生成了HC_time影像,以HC_time为合成影像索引,拓展合成了表观反射率、辐射亮度、发射率、地表发射率等系列晴空影像。利用本方法,面向不同时间跨度的数据集,可合成候晴空影像、旬晴空影像、月晴空影像。晴空合成算法改善了京津冀地区土地覆被变化时空演化规律研究的数据基础。
Based on the MODIS 1B images transponded by the DVBS (Digital Video Broadcast by Satellite) system, relying on the ENVI / IDL software, integrate MODIS 1B image pre-processing technique flow, propose a new idea for clear sky image synthesis. Pre-processing includes three main components of destriping, geometric correction and atmospheric correction. In destriping, improve equalization curve compensation, obtain good results.
There are two kinds of stripes in MODIS 1B images transponded by DVBS(Digital Video Broadcast by Satellite)system: horizontal and vertical stripes, and horizontal stripes have three types: single stripes, wide stripes and multi-line stripes. The destriping methods mainly include the frequency domain method, the spatial domain method and the statistic adjustment method. MODIS sensor has 36 spectral bands, using only one destriping method can not reach desired effect. Based on compensation and elimination of equalization curve method, this article proposes a new set of algorithm on rapidly detecting and removing stripes. Given a definite stripe detection judgment parameter,namely: the correction factor Ci of 10 groups of sensors,if the value Ci of A-frame and B-frame are separated by 0, moreover ABS (Ci)> 0.008, then there exists wide stripes; if the maximum of correction factor Ci ' is greater than 0.05, then there exists single stripes. In the wide destriping model, the current column correction factor of each sensor can be extracted from the data that the current column as the first column to continuous 100 columns of image. The experimental results show that: the striping detection model can effectively identify the single stripes and the wide stripes, multi-line stripes need visual identification .The removal effect of single and wide stripes is better than NASA website’s data, although multi-line stripes are not completely removed, image quality is significantly improved after the treatment.
Based on ENVI software, this article summarizes a set of pre-processing process for MODIS 1B image, including three main components of destriping, geometric correction and atmospheric correction. In this paper, Use Georeference MODIS module in ENVI to wipe off bow-tie and do geometric correction; Use FLAASH module to do atmospheric correction for VNIR radiation intensity bands that after geometric correction, use thermal infrared atmospheric correction module (Thermal Atm Correction) to do atmospheric correction for thermal infrared band emissivity that after geometric correction. Image quality is markedly improved after the pre-processing , advantage for retrieving the MODIS products. Image pre-processing procedures should be based on retrieval a?lgorithm of different products , in general, the data after geometric correction can be used to retrieve atmospheric products, image after atmospheric correction can be used to develop the land and ocean products. When retrieve the land and ocean products, if the indicator or algorithm has the characteristic that can reduce the atmosphere impact, then can directly use image data after geometric correction.
Compared the images composed by b1/b7 ration method with MOD09GA released by NASA,the conclusion can be leaded to that this method is good at the land part clear image composition; by contrast, the max saturation method was sutible for the water part. Combined these two methods, the special method was proposed,which regarded the land piexl of the minimum value of b1/b7 in a certain days as the clear piexl;considered the water piexl of the maximum value of saturation as the clear piexl of the water. According to the index image, producted during the compound progress of the land surface reflectance, named as HC_time, the compound clear images of Apparent reflectance, radiance, emitivity or land suface emitivity could be producted. As the time span change, the 5 days clear image,10 days clear image and month clear image could be produced. This algrithm of compound clear image improved the data basis for the research of land-use spatio-temporal evolution rule around Beijing ,Tianjin and Hebei province.
引文
[1] Parkinson, C.L.,A. Ward and M.D. King, Earth Science Reference Handbook (A Guide to NASA’s Earth Science Program and Earth Observing Satellite Missions). 1999.
[2]Corporation,L.M.,Direct Access System User's Guide for the EOS-AM Spacecraft.1998.
[3]Nishihama, M., R.Wolfe and D. Solomon,MODIS Level 1A Earth Location: Algorithm Theoretical Basis Document Version 3.0.1997.
[4]刘闯,葛成辉.美国对地观测系统(EOS)中分辨率成像光谱仪(MODIS)遥感数据的特点与应用[J].遥感信息,2000(3):45-48.
[5]王正兴,刘闯,HUETEAlfredo.植被指数研究进展:从AVHRR—NDVI到MODIS—EVI[J].生态学报,2003.23(5):979-987.
[6]刘玉洁,杨忠东.MODIS遥感信息处理原理与算法[M].2001,北京:科学出版社.
[7]郭广猛.非星历表法去除MODIS图像边缘重叠影响的研究[J].遥感技术与应用,2003. 18(3):172-175.
[8]张杰,王介民,郭铌.应用6S模式对EOS-MODIS可见光到中红外波段的大气订正[J].应用气象学报,2004(06):651-657.
[9]刘劲松,卢纪临,王卫等.DVBS系统的地学遥感数据基础平台的构建[J].河北师范大学学报(自然科学版),2008,118(02):253-256.
[10]刘劲松,卢纪临,王卫等.DVBS系统局地文件与ARC/INFO软件平台的技术整合方案[J].河北师范大学学报(自然科学版),2008,119(03):404-408.
[11]王玲,龚健雅.基于HDF文件的组织方式与影像提取[J].2003(4):35-37.
[12]相云.MODIS1B资料处理方法研究与软件实现.中国农业大学硕士学位论文.2005.
[13]Berthold K P Horn, R.J.W.,Destriping LANDSAT MSS images by Histogram Modification[J].COMPUTEB GRAPHICS AND IMAGE PROCESSING,1979(10):69-83.
[14]陈劲松,邵芸,朱博勤.中分辨率遥感图像条带噪声的去除[J].遥感学报,2004.8(3): 227-233.
[15]郭小方,王润生.基于小波分析的成像光谱图像随机点噪声消除[J].遥感学报,1999, 3(3):183-186.
[16]郭子祺,卢刚,王超.海洋SAR图像小波Speckle滤波及边缘信息提取[J].遥感学报,2001,5(06):428-434.
[17]刘光达,赵立荣.基于小波分析的医学CR影像随机噪声消除[J].光学精密工程, 2000,8(5):428-431.
[18]彭玉华.基于离散正交小波变换的图象去噪方法[J].中国图象图形学报:A辑, 1999, 4(8):677-679.
[19]杨忠东,张文健,李健.应用小波收缩方法剔除MODIS热红外波段数据条带噪声[J].遥感学报,2004,8(1):23-30.
[20]张文江,许晓东.基于小波阈值优化和边缘检测的SAR影像斑点噪声滤除[J].遥感学报,2003,7(1):41-46.
[21]黄晓园,周汝良,罗辉.MODIS影像条带噪声去除邻域插值法研究[J].地理空间信息, 2008,6(1):101-103.
[22]蒋耿明,牛铮,阮伟利. MODIS影像条带噪声去除方法研究[J].遥感技术与应用,2003,18(6):393-398.
[23]C.H.Li,P.K.S.T.,A global energy approach to facet model and its minimization using weighted least-squares algorithm[J].Pattern Recognition,2000,33:281 -293.
[24]Preesan Rakwatin,W.T.Y.Y.,Stripe Noise Reduction in MODIS Data by Combining Histogram Matching With Facet Filter[C].IEEE Transactions on geoscience and remote sensing,2007,45(6):1844-1856.
[25]吴军,张万昌.MODIS影像条带噪声去除的自相关插值法[J].遥感技术与应用,2006, 21(3):253-258.
[26]杨金红,顾松山,程明虎.插值法在去除MODIS遥感影像条带噪声中的应用[J].气象科学,2007,27(6):604-609.
[27]di Antonelli,P.B.M.,Destriping MODIS data using IFOV overlapping[J].Proc IGARSS,2004,7:4568–4571.
[28]陈劲松,邵芸,朱博勤.一种改进的矩匹配方法在CMODIS数据条带去除中的应用[J].遥感技术与应用,2003,18(5):313-316.
[29] Corsini G, D.M.W.T.,Striping removal in MOS-B data[C].IEEE Transactions on geoscience and remote sensing,2000,38(3):1439-1446.
[30]Gadallah, F.L., F. Csillag and E.J.M. Smith, Destriping multisensor imagery with moment matching[J].International Journal of Remote Sensing,2000,21(12): 2505-2511.
[31]Kautsky J, N.N.K.A.,Smoothed histogram Modification for Image Processing[J]. Computer vision,graphics,and image processing,1984,26(3):271-291.
[32]刘正军,王长耀,王成.成像光谱仪图像条带噪声去除的改进矩匹配方法[J].遥感学报,2002,6(4):279-284.
[33]牛生丽,唐军武,郭茂华.基于均衡化曲线的条带消除研究[J].遥感信息,2007(3):44-48.
[34]石光明,王晓甜,张犁.基于方向滤波器消除遥感图像孤立条带噪声的方法[J].红外与毫米波学报,2008,27(3):214-218.
[35]孙凌,唐军武,张杰.我国“海洋1”号卫星(HY—1)CCD图像数据定量化条带消除研究[J].海洋学报,2002,24(6):20-33.
[36]王杰生.TM热红外图像的横纹条带噪声及消除[J].遥感技术与应用,1995,10(01): 53-56.
[37] M, W.,Destriping multiple sensor imagery by improved histogram matching[J]. International Journal of Remote Sensing,1990,11(5):859-875.
[38]郭广锰.关于MODIS卫星数据的几何校正方法[J].遥感信息,2002(3):26-28.
[39]徐萌,郁凡,李亚春.6S模式对EOS/MODIS数据进行大气校正的方法[J].南京大学学报(自然科学版),2006(06):582-589.
[40]亓雪勇,田庆久.光学遥感大气校正研究进展[J].国土资源遥感,2005(04):1-6.
[41]宋晓宇,王纪华,刘良云.基于高光谱遥感影像的大气纠正:用AVIRIS数据评价大气纠正模块FLAASH[J].遥感技术与应用,2005,20(04):393-398.
[42]李国砚,张仲元,郑艳芬.MODIS影像的大气校正及在太湖蓝藻监测中的应用[J].湖泊科学,2008(02):160-166.
[43]ENVI4.6.1Help Document.2008.Standard FLAASH Input Parameters.
[44]杨春燕,陈圣波,包书新.长春地区云检测及云相态反演研究[J].遥感技术与应用,2010,25(2):245-250.
[45]Vermote, E.F., S.Y. Kotchenova and J.P. Ray, MODIS Surface Reflectance User’s Guide.2008.
[46]唐海蓉.基于光谱时间变化统计特征的MODIS09时间序列图像异常检测方法.2010.
[47]Levy, R.C., et al., Second-generation operational algorithm:Retrieval of aerosol properties over land from inversion of Moderate Resolution Imageing Spectroradiometer spectral reflectance[J].JOURNAL OF GEOPHYSICAL RESEARCH, 2007,112:D13211.doi :10. 1029/ 2006JD007811.
[48]王玲,田庆久,李姗姗.利用MODIS资料反演杭州市500米分辨率气溶胶光学厚度[J].遥感信息,2010,3:50-54.
[49]Kaufman, Y.J., et al., The MODIS 2.1μm Channel--Correlation with Visible Reflectance for Use in Remote Sensing of Aerosol[C].IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING,1997,35(5):1286-1300.
[50]Kutzner and Kendy, Processing MODIS Data for Fire Detection in Australia. 2001.
[51]C., R., Gonzalez and R.E. Woods,数字图象处理(第二版)(英文版)[M].2007,北京:电子工业出版社.
[52]F, V.E. and V. A, Atmospheric Correction Algorithm Spectral Reflectances. 1999.
[53]Dorothy K. Hall, G.A. and V.V.S. Riggs, Algorithm Theoretical Basis Document (ATBD) for the MODIS Snow and Sea Ice-Mapping Algorithms.2001.
[54]Zhengming, W., MODIS Land-Surface Temperature Algorithm Theoretical Basis Document.1999.
[55]Brown,O.B.and P.J.Minnett,MODIS Infrared Sea Surface Temperature Algorithm Algorithm Theoretical Basis Document.1999.