一种拓展的半物理时空融合算法及其初步应用

Landsat 5卫星较低的时间分辨率(16天)使得其很难获得大区域的、时相一致的清晰影像数据集。本文发展了一种基于半物理模型的时空融合算法-即乘性调制融合算法,并借助多时序的MODIS反射率数据来生成多时相的Landsat TM/ETM+反射率合成影像,经镶嵌后得到区域尺度的高时空分辨率地表反射率数据集(Landsat TM/ETM+)。本文利用吉林省2006年—2011年的Landsat 5 TM地表反射率数据以及500 m的MOD09A1反射率产品来生成3个时相的Landsat 5 TM反射率合成数据,从而获得研究区在上述时相下地表反射率数据的镶嵌图。初步分析表明,所生成的Landsat 5 TM反射率数据的光谱分布特征与MOD09A1反射率数据较为一致,且图像在整体上光谱特征的连续性较好。     

Spatial and Temporal Variation Analysis of Snow Cover Using MODIS over Qinghai-Tibetan Plateau during 2003–2014

The Qinghai-Tibetan Plateau (QTP) snow cover information acquisition of the high precision spatial and temporal characteristics is of great significance for the research on its land surface atmosphere coupled system and global climate change effects. The Moderate Resolution Imaging Spectro-radiometer (MODIS) daily snow cover products (MOD10A1 and MYD10A1) have been widely used in long time series of spatial and temporal variation analysis, but they are limited to be used because of high cloud cover ratio. In this paper, a 7-day rolling combination algorithm was presented to eliminate cloud obscuration, and the whole cloud amount falls below 7 %. The ground station in situ measurements verify that the overall precision is more than 90 %. The presented algorithm guaranteed the same spatial resolution and temporal resolution, and has higher precision than products MOD10A1 and MYD10A1. The MODIS 7-day rolling combination snow cover datasets products were obtained between 2003 and 2014 in the QTP, and the snow cover area of spatial and temporal variation was analyzed. The change characteristics of snow cover duration was also studied combining with the Digital Elevation Model data. Results show that the snow cover area of the whole QTP has a slowly decreased trend, but increases in autumn. Thus, the snow cover proportion of annual periodic and unstable in different elevations has the highest correlation with area of the elevation.     

A physically based approach to model LAI from MODIS 250 m data in a tropical region

A time series of leaf area index (LAI) has been developed based on 16-day normalized difference vegetation index (NDVI) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) at 250 m resolution (MOD250_LAI). The MOD250_LAI product uses a physical radiative transfer model which establishes a relationship between LAI, fraction of vegetation cover (FVC) and given patterns of surface reflectance, view-illumination conditions and optical properties of vegetation. In situ measurements of LAI and FVC made at 166 plots using hemispherical photography served for calibration of model parameters and validation of modelling results. Optical properties of vegetation cover, summarized by the light extinction coefficient, were computed at the local (pixel) level based on empirical models between ground-measured tree crown architecture at 85 sampling plots and spectral values in Landsat ETM+ bands. Influence of view-illumination conditions on optical properties of canopy was simulated by a view angle geometry model incorporating the solar zenith angle and the sensor viewing angle. The results revealed high compatibility of the produced MOD250_LAI data set with ground truth information and the 30 m resolution Landsat ETM+ LAI estimated using the similar algorithm. The produced MOD250_LAI was also compared with the global MODIS 1000-m LAI product (MOD15A2 LAI). Results show good consistency of the spatial distribution and temporal dynamics between the two LAI products. However, the results also showed that the annual LAI amplitude by the MOD15A2 product is significantly higher than by the MOD250_LAI. This higher amplitude is caused by a considerable underestimation of the tropical rainforest LAI by the MOD15A2 during the seasonal phases of low leaf production.     

MODIS数据在积雪检测中的应用

积雪作为影响环境的一个因素,是非常重要的。自1999年Terra卫星升空以来,MODIS数据在环境监测的各个方面得到了广泛的应用。由于MODIS数据的高光谱、高空间分辨率、高时间分辨率等特征,越来越多地应用到积雪监测方面。本文就MODIS数据的积雪检测算法进行了探讨,对森林中雪的检测以及云和雪的区分进行了大量的研究。结果显示：MODIS数据对积雪检测非常有效。     

基于多光谱综合的MODIS数据云检测研究

云检测是卫星遥感数据处理中不可缺少的工作。通过分析云在不同波段中的大气辐射特点，结合MODIS数据的光谱特性，提出 一种多光谱综合的云检测方法。该算法从可见光反射率、红外波段亮温值以及亮温差等方面综合考虑，逐步建立一个云检测掩模。通 过对不同时期不同背景的MODIS数据进行验证和对比分析，结果表明，该模型的云检测效果理想，尤其对可见光波段难以识别的薄卷 云也有很好效果，为有效利用MODIS数据以及进行更加精确的反演提供可靠依据。     

Landsat TM/ETM+与HJ-1A/B CCD数据自动相对辐射处理及精度验证

全球地表覆盖遥感制图与关键技术研究项目要求对两个基准年度(2000年,2010年)的全球覆盖30 m分辨率遥感数据进行辐射处理,转换到地表反射率,数据以Landsat TM/ETM+为主,HJ-1A/B CCD数据为补充。海量数据中有些不适宜进行绝对大气校正,为了保证全球覆盖,对这些数据设计开发了一套自动的相对辐射处理及精度验证流程算法,利用相邻数据重叠区域进行相对辐射校正的方式,将数据由Oigital Number(DN)值直接转换为地表反射率,精度验证以MODIS地表反射率产品MOD09GA作为参考,比较对应波段数据的相对一致性,算法采用了图像分块处理技术及OpenMP加速技术提高效率,实际应用结果表明该算法流程可以满足项目对辐射处理精度、速度及自动化程度的要求。     

GF-4增强型地表反射率库支持法的气溶胶光学厚度反演

针对GF-4等国产卫星气溶胶光学厚度反演算法存在的地表反射率估计困难、云像元污染等问题,本文发展了一种增强型地表反射率库支持的气溶胶光学厚度反演方法,改进了云筛选与地表反射率确定方案,在考虑GF-4逐像元成像角度的情况下,使用6SV模型与MOD09-CMA数据对季度尺度上的GF-4 PMS传感器数据进行大气校正,提出了百分比最小值均值法建立地表反射率库,并据此建立了NDVI与红蓝反射率关系模型,根据地表反射率的分布特点,当NDVI小于0.2的时候使用地表反射率库估计地表反射率,而当NDVI大于0.2时,则使用NDVI来估计地表反射率。使用MOD04气溶胶模式时空分布确定气溶胶参数。在京津冀地区开展气溶胶光学厚度反演实验,使用Aeronet站点数据与MOD04产品对反演结果进行了对比验证,与Aeronet相关系数R为0.964,均方根误差RMSE为0.13,满足±(0.05+0.2τ)的点多于78.9%,相关系数与均方根误差优于MODIS暗目标法产品,满足期望误差线的数量优于MODIS暗目标与深蓝算法产品。     

多样地表和大气状况下的MODIS数据云检测

云覆盖作为天气和气候变化的一个重要因子,对地表-大气能量平衡和水循环有着重要的影响,因此,快速、准确地利用卫星遥感技术检测云覆盖具有重要的实用价值和科学意义。利用卫星遥感数据,尤其是常用的Moderate Resolution Imaging Spectroradiometer(MODIS)影像数据,因其具有较高的光谱和时间分辨率,以及2330 km扫描幅宽,为大范围实时、准确地进行云检测提供了可能。目前,基于MODIS数据发展了大量的云检测方法,但因地表类型的多样性和大气状况(如空气污染和沙尘事件等)的复杂性,目前已有的云检测方法,检测精度通常具有较大的不确定性,且针对不同地表和大气状况缺乏普适性,同时也缺乏对检测精度的定量化评估。因此,本文首先比较了常用的3种云检测算法,并基于前人经验提出了两种改进方法(方法4和方法5),首先区分出云和冰雹,摒弃了不稳定的亮温波段,两种算法均适用于复杂地表和大气状况的云检测算法。结果显示,方法5可以较好地应用于基于MODIS数据的云检测,总体精度达92.6±7%,改进了现有基于MODIS数据的云检测算法;方法4平均总体精度82.9±13%,虽然精度相对较低,但云残留少,适合作为对云敏感度高的研究工作的云检测方法。     

Automatic and improved radiometric correction of Landsat imagery using reference values from MODIS surface reflectance images

Radiometric correction is a prerequisite for generating high-quality scientific data, making it possible to discriminate between product artefacts and real changes in Earth processes as well as accurately produce land cover maps and detect changes. This work contributes to the automatic generation of surface reflectance products for Landsat satellite series. Surface reflectances are generated by a new approach developed from a previous simplified radiometric (atmospheric + topographic) correction model. The proposed model keeps the core of the old model (incidence angles and cast-shadows through a digital elevation model [DEM], Earth–Sun distance, etc.) and adds new characteristics to enhance and automatize ground reflectance retrieval. The new model includes the following new features: (1) A fitting model based on reference values from pseudoinvariant areas that have been automatically extracted from existing reflectance products (Terra MODIS MOD09GA) that were selected also automatically by applying quality criteria that include a geostatistical pattern model. This guarantees the consistency of the internal and external series, making it unnecessary to provide extra atmospheric data for the acquisition date and time, dark objects or dense vegetation. (2) A spatial model for atmospheric optical depth that uses detailed DEM and MODTRAN simulations. (3) It is designed so that large time-series of images can be processed automatically to produce consistent Landsat surface reflectance time-series. (4) The approach can handle most images, acquired now or in the past, regardless of the processing system, with the exception of those with extremely high cloud coverage. The new methodology has been successfully applied to a series of near 300 images of the same area including MSS, TM and ETM+ imagery as well as to different formats and processing systems (LPGS and NLAPS from the USGS; CEOS from ESA) for different degrees of cloud coverage (up to 60%) and SLC-off. Reflectance products have been validated with some example applications: time series robustness (for a pixel in a pseudoinvariant area, deviations are only 1.04% on average along the series), spectral signatures generation (visually coherent with the MODIS ones, but more similar between dates), and classification (up to 4 percent points better than those obtained with the original manual method or the CDR products). In conclusion, this new approach, that could also be applied to other sensors with similar band configurations, offers a fully automatic and reasonably good procedure for the new era of long time-series of spatially detailed global remote sensing data.     

Comparison of eMODIS and MOD/MYD13A2 NDVI products during 2012–2014 spring green-up periods in Alaska and northwest Canada

Accurate monitoring of vegetation dynamics is required to understand the inter-annual variability and long term trends in terrestrial carbon exchange in tundra and boreal ecoregions. In western North America, two Normalized Vegetation Index (NDVI) products based on spectral reflectance data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are available. The MOD/MYD13A2 NDVI product is available as a 16-day composite product in a sinusoidal projection as global hdf tiles. The eMODIS Alaska NDVI product is available as a 7-day composite geotif product in a regional equal area conic projection covering Alaska and the entire Yukon River Basin. These two NDVI products were compared for the 2012–2014 late May–late June spring green-up periods in Alaska and the Yukon Territory. Relative to the MOD/MYD13A2 NDVI product, it is likely that the eMODIS NDVI product contained more cloud-contaminated NDVI values. For example, the MOD/MYD13A2 product flagged substantially fewer pixels as “good quality” in each 16-day composite period compared to the corresponding MODIS Alaska NDVI product from a 7-day composite period. During the spring green-up period, when field-based NDVI increases, the eMODIS NDVI product averaged 43 percent of pixels that declined by at least 0.05 NDVI between 2 composite periods, consistent with cloud-contamination problems, while the MOD/MYD13A2 NDVI averaged only 6 percent of pixels. Based on a cloudy Landsat-8 scene, the eMODIS compositing process selected 23 percent pixels, while the MOD/MYD13A2 compositing process selected less than 0.003 percent pixels. Based on the results, it appears that the MOD/MYD13A2 NDVI product is superior for scientific applications based on NDVI phenology in the tundra and boreal regions of northwestern North America.     

Estimation of regional evapotranspiration over the North China Plain using geostationary satellite data

Data from the first operational Chinese geostationary satellite Fengyun-2C (FY-2C) satellite are applied in combination with Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products for the assessment of regional evapotranspiration over the North China Plain. The approach is based on the improved triangle method, where the temperature–vegetation index space includes thermal inertia. Two thermal infrared channels from FY-2C are used to estimate surface temperature (Ts) based on a split window algorithm originally proposed for the MSG-SEVIRI sensor. Subsequently the high temporal resolution of FY-2C data is exploited to give the morning rise in Ts. Combined with the 16 days composite MODIS vegetation indices product (MOD13) at a spatial resolution of 5 km, evaporative fraction (EF) is estimated by interpolation in the ΔTs–NDVI triangular-shaped scatter space. Finally, regional actual evapotranspiration (ET) is derived from the evaporative fraction and available energy estimated from MODIS surface albedo products MCD43. Spatial variations of estimated surface variables (Ts, EF and ET) corresponded well to land cover patterns and farmland management practices. Estimated ET and EF also compared well to lysimeter data collected for the period June 2005–September 2007. The improved triangle method was also applied to MODIS products for comparison. Estimates based on FY-2C products proved to provide slightly better results than those based on MODIS products. The consistency of the estimated spatial variation with other spatial data supports the use of FY-2C data for ET estimation using the improved triangle method. Of particular value is the high temporal frequency of image acquisitions from FY-2C which improves the likelihood of obtaining cloud free image acquisitions as compared to polar orbiting sensors like MODIS.     

多时相MODIS影像水田信息提取研究

水稻种植及其分布信息是土地覆被变化、作物估产、甲烷排放、粮食安全和水资源管理分析的重要数据源。基于遥感的水田利用监测中,通常采用时序NDVI植被指数法和影像分类法分别进行AVHRR和TM影像的水田信息获取。针对8天合成MODIS陆地表面反射比数据的特点和水稻生长特征,选取水稻种植前的休耕期、秧苗移植期、秧苗生长期和成熟期等多时相MODIS地表反射率影像数据,通过归一化植被指数、增强植被指数及利用对土壤湿度和植被水分含量较敏感的短波红外波段计算得到的陆表水指数进行水田信息获取。将提取结果与基于ETM＋影像的国土资源调查水田数据,通过网格化计算处理并进行对比分析,结果表明,利用MODIS影像的8天合成地表反射率数据,进行区域甚至全国的水田利用监测是可行的。     

基于MODIS影像的夜间云检测算法研究

在以往云检测算法的基础上,利用MODIS传感器高光谱和高时间分辨率的特点,建立一套针对MODIS夜间影像的,以单、多时相组合方法为基础的夜间云检测算法.通过对我国南北地震构造带(南北带)影像进行云检测试验证实,该算法对MODIS夜间影像上的各种云类具有较好的检测效果.     

利用改进型MODIS火点探测算法实现河北省秸秆焚烧火点识别

为准确地了解河北省秸秆焚烧火点的空间分布,为秸秆焚烧监测的实现、禁烧工作的开展、环境质量改善提供支持．基于MODIS L1B数据、MODIS标准火点产品MOD14、全国秸秆焚烧火点日报数据为基础,采用改进型MODIS火灾探测算法,并通过IDL语言实现,得到秸秆焚烧火点空间分布信息,并进行空间与定量精度分析．研究表明:火点大部分位于河北省南部的一些地区,其中尤以邢台、石家庄、邯郸火点数量最为突出;该算法运算速度快,获取的秸秆焚烧火点数据具有一定检测精度和可靠性,对秸秆焚烧的监测具有一定的实用价值．      

联合云量自动评估和加权支持向量机的Landsat图像云检测

针对ACCA（云量自动评估）算法难以检测Landsat图像中的半透明云问题,提出了一种ACCA和WSVM（加权支持向量机）相结合的云检测算法．首先根据云在不同波段中的大气辐射特点,结合Landsat ETM+图像数据的光谱特性,利用ACCA算法将图像像元初步分成云像元、非云像元和待定像元,再以云的光谱特性构造特征向量,利用WSVM算法进行待定像元的云层检测,最终获得全部图像的云检测结果．仿真实验结果表明,该方法既具有ACCA算法的云检测优势,还对ACCA算法难以识别的半透明云有很好的检测效果．     

FY-3A/MERSI海上沙尘天气气溶胶光学厚度反演

利用FY-3A/MERSI资料,结合MODIS C005算法的海上气溶胶模型,研究了中国渤海、黄海以及东海海域沙尘天气气溶胶光学厚度反演方法。通过MERSI反演的气溶胶光学厚度与MODIS C005气溶胶光学厚度产品(MOD04)对比,发现MERSI气溶胶光学厚度反演结果存在较大偏差。分析认为,这种偏差是由MODIS C005算法的海上气溶胶模型对MERSI气溶胶光学厚度反演不完全适用造成的。鉴于此,本文研究引进了一种沙尘气溶胶模型,并将其与MODIS C005算法的粗粒子气溶胶模型按照一定比例混合,形成了改进的气溶胶模型。利用改进气溶胶模型再次反演海上沙尘天气气溶胶光学厚度,反演结果与MOD04一致性较好,说明改进气溶胶模型能有效地提高MERSI定量反演沙尘气溶胶的能力。     

PARASOL/POLDER3卫星数据的海洋上空云检测

在卫星海洋遥感中,云作为海气耦合系统最重要的调节器之一,其检测结果对海洋上空云微物理特性的反演精度有较大影响。因此,快速而准确识别海洋上空的云像元是卫星遥感数据处理过程中首要解决的关键问题。以PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar)卫星搭载的POLDER3载荷遥感数据为研究对象,提出一种改进的海洋上空云检测方法。首先剔除海洋耀光;接着利用有云与晴空区近红外反射率差异检验识别有云像元,并利用偏振反射率检验进一步识别低反射率的云像元;然后利用近红外与可见光反射率比值检验识别晴空像元;最后建立多角度云检测结果空间融合规则,重新标记有云、晴空和未定像元。以印度洋海区为例进行实验分析,将云检测结果与Buriez方法进行对比,发现检测精度基本相当,而有云像元的识别速度却平均提高约3倍。结果表明:该方法能有效的检测出海洋上空的云像元,满足业务化数据处理的高精度及时效性要求,为后续云微物理特性反演提供可靠的数据源。     

利用BRDF原型和单方向反射率数据估算地表反照率

地表反照率是影响地表能量收支平衡的决定性参数之一,精确反演地表反照率需要考虑地表各向异性反射特征。本文尝试以双向反射分布函数BRDF原型为地表各向异性反射的先验知识,通过单方向反射率反演地表反照率。首先根据地面实测及MODIS多角度反射率数据对反演方法进行分析和精度评价,然后借助MODIS BRDF产品统计出研究区的主导BRDF原型,并联合环境一号卫星(HJ-1)单方向反射率数据反演30 m地表反照率,最终将结果与地表实测数据进行比较。结果表明:BRDF原型对BRDF的变化进行了约束,且能够适用于几十米尺度的遥感数据反照率的反演;不同级别的各向异性反射特征的分布是不均一的,借助于主导BRDF原型能够使大部分样本的地表反照率满足精度要求;利用研究区MODIS BRDF产品统计得到的主导BRDF原型为先验知识,通过HJ-1数据反演得到的地表反照率与地表实测反照率有较高的一致性,而朗伯假定条件下的反照率高于实测结果。本文算法简单高效,可为产生全国范围的中高分辨卫星反照率产品提供有价值的算法参考。     

MODIS NDVI时间序列数据的去云算法比较

受多重因素的影响,MODIS NDVI数据产品中存在着大量的噪声,需要进行去噪重建.针对目前几种常用的NDVI时间序列数据去云方法,如HANTS法、SPLINE插值法以及Savizky-Golay法,以山东省MODIS NDVI时间序列数据(一年的)作为检验数据,从不同角度比较几种算法的去云能力和使用范围.结果表明:S...     

Optimal MODIS data processing for accurate multi-year paddy rice area mapping in China

ABSTRACT

Researchers, policy makers, and farmers currently rely on remote sensing technology to monitor crops. Although data processing methods can be different among different remote sensing methods, little work has been done on studying these differences. In order for potential users to have confidence in remote sensing products, an analysis of mapping accuracies and their associated uncertainties with different data processing methods is required. This study used the MOD09A1 and MYD09A1 products of the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite, from which the Enhanced Vegetation Index (EVI) and the two-band EVI (EVI2) images were obtained. The objective of this study was to analyze the accuracy of different data processing combinations for multi-year rice area mapping. Sixteen combinations of EVI and EVI2 with two cloudy pixel removal methods (QA/BLUE) and four pixel replacement methods (MO/MY/MOY/MYO) were investigated over the Jiangsu Province of southeast China from 2006 to 2016. Different accuracy results were obtained with different data processing combinations for multi-year rice field mapping. Based on a comparison of the relative performance of different MODIS products and processing method combinations, EVI2_BLUE_MYO was proposed to be the optimal processing method, and was applied to forecasting the rice-planted area of 2017. Study results from 2006 to 2017 were validated against reference data and showed accuracies of rice area extraction of greater than 95%. The mean absolute error of transplanting, heading, and maturity dates were 11.55, 8.10, and 7.78 days, respectively. In 2017, two sample regions (A and B) were selected from places where rice fractional cover was greater than 75%. Rice area extraction accuracies of 85.0% (A) and 92.3% (B) were obtained. These results demonstrated the complementarity of MOD09A1 and MYD09A1 datasets in enhancing pixel spatial coverage and improving rice area mapping when atmospheric influences are significant. The optimal data processing combination indentified in this study is promising for accurate multi-year and large-area paddy rice information extraction and forecasting.