遥感时间序列数据滤波重建算法发展综述

遥感时间序列数据(MODIS,NOAA/AVHRR,SPOT/VEGETATION等)在植被生长监测、物候信息提取、土地利用类型监测等诸多领域得到了广泛应用,是生产研究的重要数据源之一.由于传感器、云层大气等影响,遥感时间序列数据存在着严重的噪声,应用前必须进行序列滤波重建工作.综述现有各类滤波重建方法,对研究中广为采用的3类主要方法(基于最小二乘的非对称高斯函数拟合、SavitZky-Golay滤波、基于离散傅里叶的系列分析方法)集中阐述其理论基础、应用步骤和优缺点.总结当前遥感时间序列滤波重建方法需要进一步改进之处.     

Biome-level relationships between vegetation indices and climate variables using time-series analysis of remotely-sensed data

ABSTRACT

Climatic factors such as rainfall and temperature play a vital role in the growth characteristics of vegetation. While the relationship between climate and vegetation growth can be accurately predicted in instances where vegetation is homogenous, this becomes complex to determine in heterogeneous vegetation environments. The aim of this paper was to study the relationship between remotely-sensed monthly vegetation indices (i.e. Normalized Difference Vegetation Index and Enhanced Vegetation Index) and climatic variables (temperature and precipitation) using time-series analysis at the biome-level. Specifically, the autoregressive distributed lag model (ARDL1 and ARDL2, corresponding respectively to one month and two month lags) and the Koyck-transformed distributed lag model were used to build regression models. All three models estimated NDVI and EVI fairly accurately in all biomes (Relative Root-Mean-Squared-Error (RMSE): 12.0–26.4%). Biomes characterized by relative homogeneity (Grassland, Savanna, Indian Ocean Coastal Belt and Forest Biomes) achieved the most accurate estimates due to the dominance of a few species. Comparisons of lag size (one month compared to two months) generally showed similarities (Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC) and log-likelihood) with quite high comparability in certain biomes – this indicates the utility of the ARDL1 and ARDL2 model, depending on the availability of appropriate data. These findings demonstrate the variation in estimation linked to the biome, and thus the validity of biome-level correlation of climatic data and vegetation indices.     

采用DBM方法的时间序列LAI建模与估算

运用DBM(Data Based Mechanistic)方法,使用MODIS数据,建立了遥感观测反射率数据与叶面积指数(LAI)在时间序列上的统计关系模型(LAI_DBM模型),并结合部分Bigfoot站点实测LAI数据进行了模型检验。结果显示,LAI_DBM模型能够较好表达时间序列反射率与LAI的动态变化关系。LAI_DBM模型使用遥感观测数据实时估算得到的LAI,在数据质量和时间连续性上比MODISLAI有改进。     

A data mining approach for evaluation of optimal time-series of MODIS data for land cover mapping at a regional level

Optical Earth Observation data with moderate spatial resolutions, typically MODIS (Moderate Resolution Imaging Spectroradiometer), are of particular value to environmental applications due to their high temporal and spectral resolutions. Time-series of MODIS data capture dynamic phenomena of vegetation and its environment, and are considered as one of the most effective data sources for land cover mapping at a regional and national level. However, the time-series, multiple bands and their derivations such as NDVI constitute a large volume of data that poses a significant challenge for automated mapping of land cover while optimally utilizing the information it contains. In this study, time-series of 10-day cloud-free MODIS composites and its derivatives – NDVI and vegetation phenology information, are fully assessed to determine the optimal data sets for deriving land cover. Three groups of variable combinations of MODIS spectral information and its derived metrics are thoroughly explored to identify the optimal combinations for land cover identification using a data mining tool.The results, based on the assessment using time-series of MODIS data, show that in general using a longer time period of the time-series data and more spectral bands could lead to more accurate land cover identification than that of a shorter period of the time-series and fewer bands. However, we reveal that, with some optimal variable combinations of few bands and a shorter period of time-series data, the highest possible accuracy of land cover classification can be achieved.     

植被指数时序数据距离测度方法评价

利用MODIS增强型植被指数(EVI)时序数据,基于中国陆地生态系统55种植被类型上的468个测试点和一个测试区进行了实验,综合比较欧氏距离、光谱信息离散度、光谱角余弦、核光谱角余弦、相关系数、光谱角余弦-欧氏距离6种距离测度方法对遥感植被指数时序数据聚类精度的影响,结果表明:相关系数方法的聚类精度最差;光谱角余弦-欧氏距离方法充分利用了植被指数时序数据的曲线幅度和形状特征,在这6种距离测度方法中表现出了最优的聚类效果;只对光谱亮度敏感的欧氏距离方法或只对曲线形状敏感的光谱角余弦方法,无论是在区分地物类型方面,还是在区域应用上,表现效果均较差;核光谱角余弦虽然在点数据测试上表现较差,但在区域应用上却有较好的表现;光谱信息离散度无论是在点数据测试上还是在区域应用上均表现出了较为适中的效果。     

Object-based image analysis supported by data mining to discriminate large areas of soybean

This research aimed to analyze the possibility to estimate and automatically map large areas of soybean cultivation through the use of MODIS (Moderate-Resolution Imaging Spectroradiometer) images. Two major techniques were used: GEOgraphic-Object-Based Image Analysis (GEOBIA) and Data Mining (DM). In order to obtain the images, the segmentation algorithm implemented by Definiens Developer was used. A decision tree (DT) was created from a training set previously prepared. Time-series of images from the MODIS sensor aboard the Terra satellite were acquired in order to represent the wide variation of the vegetation pattern along the soybean crop cycle. The time-series data were used only for the CEI index. Furthermore, to compare the results obtained from GEOBIA, the slicing technique was used at the CEI level. After the training, the DT was applied to the vegetation indices generating the thematic map of the spatial distribution of soybean. In accordance with the error matrix and kappa parameter analysis, tests for statistical significance were created. Results indicate that the classification achieved by Kappa coefficients is 0.76. In short, the obtained results proved that combining vegetation indices and time-series data using GEOBIA return promising results for mapping soybean plantation on a regional scale.     

Earth observation data based rapid flood-extent modelling for tsunami-devastated coastal areas

Earth observation (EO)-based mapping and analysis of natural hazards plays a critical role in various aspects of post-disaster aid management. Spatial very high-resolution Earth observation data provide important information for managing post-tsunami activities on devastated land and monitoring re-cultivation and reconstruction. The automatic and fast use of high-resolution EO data for rapid mapping is, however, complicated by high spectral variability in densely populated urban areas and unpredictable textural and spectral land-surface changes. The present paper presents the results of the SENDAI project, which developed an automatic post-tsunami flood-extent modelling concept using RapidEye multispectral satellite data and ASTER Global Digital Elevation Model Version 2 (GDEM V2) data of the eastern coast of Japan (captured after the Tohoku earthquake). In this paper, the authors developed both a bathtub-modelling approach and a cost-distance approach, and integrated the roughness parameters of different land-use types to increase the accuracy of flood-extent modelling. Overall, the accuracy of the developed models reached 87–92%, depending on the analysed test site. The flood-modelling approach was explained and results were compared with published approaches. We came to the conclusion that the cost-factor-based approach reaches accuracy comparable to published results from hydrological modelling. However the proposed cost-factor approach is based on a much simpler dataset, which is available globally.     

Correcting superconducting gravity time-series using rainfall modelling at the Vienna and Membach stations and application to Earth tide analysis

We demonstrate the possibility to improve the signal-to-noise ratio of superconducting gravity time-series by correcting for local hydrological effects. Short-term atmospheric events associated with heavy rain induce step-like gravity signals that deteriorate the frequency spectrum estimates. Based on 4D modeling constrained by high temporal resolution rain gauge data, rainfall admittances for the Vienna and Membach superconducting gravity stations are calculated. This allows routine correction for Newtonian rain water effects, which reduces the standard deviation of residuals after tidal parameter adjustment by 10%. It also improves the correction of steps of instrumental origin when they coincide with step-like water mass signals.     

Crop discrimination in Northern China with double cropping systems using Fourier analysis of time-series MODIS data

Crop identification is the basis of crop monitoring using remote sensing. Remote sensing the extent and distribution of individual crop types has proven useful to a wide range of users, including policy-makers, farmers, and scientists. Northern China is not merely the political, economic, and cultural centre of China, but also an important base for grain production. Its main grains are wheat, maize, and cotton. By employing the Fourier analysis method, we studied crop planting patterns in the Northern China plain. Then, using time-series EOS-MODIS NDVI data, we extracted the key parameters to discriminate crop types. The results showed that the estimated area and the statistics were correlated well at the county-level. Furthermore, there was little difference between the crop area estimated by the MODIS data and the statistics at province-level. Our study shows that the method we designed is promising for use in regional spatial scale crop mapping in Northern China using the MODIS NDVI time-series.     

Mapping crop phenology using NDVI time-series derived from HJ-1 A/B data

With the availability of high frequent satellite data, crop phenology could be accurately mapped using time-series remote sensing data. Vegetation index time-series data derived from AVHRR, MODIS, and SPOT-VEGETATION images usually have coarse spatial resolution. Mapping crop phenology parameters using higher spatial resolution images (e.g., Landsat TM-like) is unprecedented. Recently launched HJ-1 A/B CCD sensors boarded on China Environment Satellite provided a feasible and ideal data source for the construction of high spatio-temporal resolution vegetation index time-series. This paper presented a comprehensive method to construct NDVI time-series dataset derived from HJ-1 A/B CCD and demonstrated its application in cropland areas. The procedures of time-series data construction included image preprocessing, signal filtering, and interpolation for daily NDVI images then the NDVI time-series could present a smooth and complete phenological cycle. To demonstrate its application, TIMESAT program was employed to extract phenology parameters of crop lands located in Guanzhong Plain, China. The small-scale test showed that the crop season start/end derived from HJ-1 A/B NDVI time-series was comparable with local agro-metrological observation. The methodology for reconstructing time-series remote sensing data had been proved feasible, though forgoing researches will improve this a lot in mapping crop phenology. Last but not least, further studies should be focused on field-data collection, smoothing method and phenology definitions using time-series remote sensing data.     

MODIS数据存储格式研究

新型的遥感数据MODIS在近年来得到了比较广泛的应用,该数据有很多优点,在地理科学与资源管理中的应用潜力很大。深入研究该数据采用的HDF文件的数据结构,是开发自主知识产权的遥感图像处理软件的前提工作。本文就:什么是HDF,为什么创造HDF,HDF的高等级应用程序接口,HDF命令行用法和可视化工具,主要的HDF平台,HDF的基本原理,HDF文件格式,用多文件接口对HDF文件的基本操作,确定一个文件是否为HDF文件等问题展开了讨论。     

基于时间序列法的高速公路沉降预测

公路的不均匀沉降会对车辆运行带来极大的安全隐患,高速公路沉降观测和预测是保障高速公路安全运营的必要手段.本文将时间序列法用于高速公路沉降预测,分析时间序列法的建模步骤、参数估计和适应性分析原理.通过程序计算得到预测结果和精度,与传统的双曲线法比较发现,时间序列法具有更高的预测精度和更好的工程实用性.     

Quantifying annual land-cover change and vegetation greenness variation in a coastal ecosystem using dense time-series Landsat data

Land-cover change may affect water and carbon cycles when transitioning from one land-cover category to another (land-cover conversion, LCC) or when the characteristics of the land-cover type are altered without changing its overall category (land-cover modification, LCM). Given the increasing availability of time-series remotely sensed data for earth monitoring, there has been increased recognition of the importance of accounting for both LCC and LCM to study annual land-cover changes. In this study, we integrated 1,513 time-series Landsat images and a change-updating method to identify annual LCC and LCM during 1986–2015 in the coastal area of Zhejiang Province, China. The purpose was to quantify their contributions to land-cover changes and impacts on the amount of vegetation. The results show that LCC and LCM can be successfully distinguished with an overall accuracy of 90.0%. LCM accounted for 22% and 40.5% of the detected land-cover changes in reclaimed and inland areas, respectively, during 1986–2015. In the reclaimed area, LCC occurred mostly in muddy tidal flats, construction land, aquaculture ponds, and freshwater herbaceous land, whereas LCM occurred mostly in freshwater herbaceous land, Spartina alterniflora, and muddy tidal flats. In the inland area, both LCC and LCM were concentrated in forest and dryland. Overall, LCC had a mean magnitude of normalized difference vegetation index (NDVI) change similar to that of LCM. However, LCC had a positive effect and LCM had a negative effect on NDVI change in the reclaimed area. Both LCC and LCM in the inland area had negative impacts on vegetation greenness, but LCC resulted in larger NDVI change magnitude. Impacts of LCC and LCM on vegetation greenness were quantified for each land-cover type. This study provided a methodological framework to take both LCC and LCM into account when analyzing land-cover changes and quantified their effects on coastal ecosystem vegetation.     

MODIS和HJ-1 CCD数据时空融合重构NDVI时间序列

遥感数据反演高时空分辨率NDVI对监测植被动态变化过程具有重要意义,然而受天气影响,单颗卫星难以提供时间连续的高空间分辨率NDVI数据。以华北平原中东部为实验区,联合HJ-1 CCD数据和MODIS数据,对STARFM算法进行了改进,(1)考虑了不同地物对光谱响应的差异,为减少分类错误利用统计学上()对分类数据进行筛选,按照不同地物类型分别利用线性拟合方法修改光谱距离权重;(2)定义了预测半径,对HJ-1 CCD数据因外界影响而缺失的影像进行了预测。结果表明,与真实影像相比,预测结果呈现了较好的空间一致性,相关系数均达到了极显著相关,改进算法的预测精度要高于原算法。利用该方法将HJ-1 CCD NDVI的空间变化信息与MODIS NDVI时间变化信息有机结合重构了高时空分辨率NDVI序列,有效补充了HJ-1CCD NDVI的缺失数据集。     

Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data

Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha⁻¹ (RMSE% = 46.3%) and R² value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha⁻¹ and 135.9 Mg ha⁻¹ for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha⁻¹ to 112.27 Mg ha⁻¹ and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha⁻¹ to 550 Mg ha⁻¹ of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha⁻¹ to 600 Mg ha⁻¹ of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.     

Big Earth data: disruptive changes in Earth observation data management and analysis?

ABSTRACT

Turning Earth observation (EO) data consistently and systematically into valuable global information layers is an ongoing challenge for the EO community. Recently, the term ‘big Earth data’ emerged to describe massive EO datasets that confronts analysts and their traditional workflows with a range of challenges. We argue that the altered circumstances must be actively intercepted by an evolution of EO to revolutionise their application in various domains. The disruptive element is that analysts and end-users increasingly rely on Web-based workflows. In this contribution we study selected systems and portals, put them in the context of challenges and opportunities and highlight selected shortcomings and possible future developments that we consider relevant for the imminent uptake of big Earth data.     

地球观测数据共享的发展和趋势

地球观测数据共享是地球科学和相关学科科研活动中非常重要的基础性工作,是对地观测信息生命周期中的重要环节。受到由资源提供者、资源消费者和资源加工者组成的社会生态系统发展变化的影响,共享模式经历了无共享、项目共享、部门共享、社会共享等渐进的4个发展阶段,并呈现出区域差异和阶段差异。地球观测数据共享的概念体系包含数据开放、数据共享、数据互联等不同层次的问题,并受到信息技术等使能技术的驱动。其中开放性代表数据在网络中可被访问的状态,共享性是对于数据重复使用的授权和模式,互联性则是强调可共享数据资源在科学含义上的相互理解。而地球观测数据共享的技术体系则包含数据开放技术、数据共享技术和数据出版与引用技术。目前地球观测领域的数据共享正在经历巨大的文化、政策、技术和应用变革,下一代的地球观测数据设施集中体现了数据的共享和协作,并将呈现国际化、多学科化、标准化、设施化、大数据化和公众社会化等新的技术特征,将对相关科学活动产生重大影响。     

基于小波消噪的时序分析改进法在GPS变形监测中的应用

考虑到小波技术在粗差探测、消噪方面有特殊的作用,本文结合实际算例,把小波分析与传统数据处理技术结合作为研究对象,提出一种基于小波消噪的时序分析改进法。并对这种处理方法在GPS监测中的应用作了充分地论证。结果表明,经过小波处理后,时序分析法在预报方面有明显的进步,小波处理后的时序分析法比小波处理前的精度高,从而说明了这种改进方法的可行性。     

一种结合整型最佳谱间预测与SPIHT的MODIS影像无损压缩算法

针对MODIS(Moderate Resolution Imaging Spectroradiometer)影像数据海量并具有重要研究价值的特点,研究MODIS影像的无损压缩算法。采用最佳线性预测方法,通过波段相关性排序确定波段最优预测顺序,并自适应计算"预测波段"与"当前波段"的二阶最佳预测器系数,减少谱间冗余;以多级树集合分裂树(Set Partitioning In Hierarchical Trees,SPIHT)编码算法降低谱内相关。为确保无损压缩,对线性预测系数进行最佳逼近取整操作,并采用基于提升格式的D5/3整数小波变换。实验结果表明本文提出的算法在压缩比上性能较3DSPIHT等算法突出。     

边缘系数建筑物雷达点云数据简化方法

针对基于LiDAR点云数据进行建筑物自动重建中存在的数据冗余问题,该文设计了一种定量描述激光点位于地物边缘区几率大小的指标——边缘系数,并据此提出了基于边缘系数的建筑物LiDAR点云数据简化方法。该方法利用激光点与其邻域点的位置、数量及分布计算该点的边缘系数,通过试验分析确定边缘系数的阈值并对点云数据进行分割,最后保留建筑物边缘区域的点,实现点云数据的简化。实验表明,该方法在对点云数据进行高效压缩的同时有效保留了位于地物边缘处的点云,有助于提高海量点云数据处理能力和建筑物重建效率。