基于GIMMS、VGT和MODIS的中国东部植被指数对比分析

GIMMS NDVI、VGT NDVI和MODIS NDVI/EVI是目前在植被变化有关研究中经常使用的植被遥感数据,它们之间的差异也得到了广泛关注。然而,在分析这些数据之间的差异时,较少有研究注意到植被本身固有的季节循环可能夸大了各数据间的相关关系。应用2000~2006年GIMMS NDVI、VGT NDVI、MODIS NDVI/EVI等不同植被遥感数据,对比了基于这些数据集的中国东部植被年际变化的差异,探讨了植被季节循环对不同遥感数据之间相关性的影响。结果表明:由不同遥感数据提取的植被年际变化特征具有明显的一致性,然而,植被本身固有的季节循环特征掩盖了不同数据集的差异。季节循环去除前,各数据集之间具有显著的相关性;季节循环去除后,各数据集的相关性明显降低,但不同数据集在北部区域依然具有较好的一致性,其差异主要出现在南部区域,差异最明显的是GIMMS与MODIS数据,二者在淮河以南的区域几乎不存在显著相关。     

基于植被分区的中国植被类型分类方法

高精度的土地覆盖分类产品对定量遥感研究及遥感应用等具有非常重要的意义。目前免费的且全球覆盖的土地分类产品已有很多,但这些产品多为国外研究机构和人员所研发,由于对中国区域地形复杂、植被结构特征差异与农作物种植结构差异等没有进行充分的研究,使得这些产品用于中国区域的分类时其精度尤其是植被类型的分类精度较低。因此,生产一种针对中国区域的植被类型分类产品是非常必要的。针对中国区域地形、土壤等信息,并在借鉴现有的植被区划的基础上,发展了一种基于植被分区的中国植被类型分类方法,该分类方法以长时间序列为基础,能以较高的时间分辨率捕捉地表随时间变化的信息,从而利用地物在时间维上的差异提高分类精度,并利用该方法完成了2012年中国土地覆盖分类。此外还通过分层随机采样的方法对分类结果进行了精度评估,发现本分类产品的总体精度和Kappa系数有较大提高,其中本文产品总体精度为90.78%,Kappa系数为0.86;并通过与MODIS土地覆盖数据产品进行比较,发现该产品精度比MODIS土地覆盖数据产品在植被类型上提高了61.38%。     

多云雾地区高时空分辨率植被覆盖度构建方法研究

针对多云雾地区高时空分辨率数据缺乏现状,提出了一套区域尺度高时空分辨率植被覆盖度数据构建方法.首先,通过时空适应反射率融合模型(STARFM)有效地将TM 的较高空间分辨率与MODIS的高时间分辨率融合在一起,构建了研究区植被生长峰值阶段的NDVI数据;然后,以植被生长峰值阶段的NDVI为输入,基于地表覆被类型,综合应用等密度和非密度亚像元模型对研究区的植被覆盖度进行估算.结果表明:①即使数据源存在大量的云雾,且存在一定的时相差异,研究区植被覆盖度的估算结果过渡自然,不存在明显的不接边效应;②以植被生长峰值阶段的NDVI数据为输入进行植被覆盖度估算,有效拉开了同一地表覆被类型不同覆盖度像元的NDVI梯度,提高了亚像元估算模型对输入数据的抗扰动性;③基于地表覆被类型,应用亚像元混合模型,能够提高植被覆盖度的估算精度.经野外实测数据验证,总体约８５％的估算精度表明,针对高时空分辨率遥感数据缺乏的多云雾区域,本研究提出的方法能够实现区域尺度植被覆盖度数据的构建.     

2001~2010年松木希错流域植被动态变化遥感研究

遥感在区域植被变化研究中具有十分重要的作用,能为大面积监测植被状况的演化过程提供技术支持。NDVI在高植被覆盖地区存在过饱和现象,对稀疏地区的植被变化尤其敏感。以古里雅冰帽南部的松木希错流域植被相对稀疏区域为研究区,基于MODIS NDVI数据和逐月气象观测数据,以及RS和GIS平台,对该区域2001~2010年主要植被变化趋势进行了初步研究,并对植被变化与气候驱动因子的关系进行了分析和探讨。结果表明:① 2001~2010年间该区域的植被活动有加强趋势;② NDVI表明研究区植被生长季较短(5~9月),NDVI浮动区间为0.11～0.13,低于全国水平(0.3～0.35),也低于全球稀疏灌丛的平均水平(0.2～0.4);③NDVI与年均气温整体上呈正相关,而与年降水量相关性不强。表明近年来持续升温是影响该区域植被活动加强的最主要原因。     

环境星归一化植被指数时间序列滤波算法比较

由于环境卫星成像条件和卫星平台的限制,针对环境卫星归一化植被指数(HJ NDVI)时间序列中存在较多噪音的现象,该文比较分析了Savitzky-Golay(S-G)滤波法、非对称高斯函数拟合法(AG)、双逻辑曲线拟合法(DL)和时间序列谐波分析法(Hants)4种滤波算法。建立了南京市典型植被类型区域的HJ NDVI时间序列,对4种滤波方法进行实验。对比纯像元样点、样区的滤波结果以及滤波后5类典型植被的分类精度,评价4种滤波方法的滤波效果,并利用MODIS NDVI时序数据验证结果可靠性。结果表明:滑动窗口大小为5的S-G滤波的滤波效果最佳。该研究结论为基于HJ NDVI时间序列的应用研究滤波方法选择提供参考。     

HJ-1A NDVI与MODIS NDVI时间序列提取植被物候特征对比研究

以HJ-1A和MODIS为数据源,通过动态阈值法提取物候特征参数,对HJ-1A NDVI和MODIS NDVI时间序列进行植被物候特征提取进行定性和定量比较,通过比较结果,提出HJ-1A NDVI数据在该应用中存在的问题,促进国产中空间高时间分辨率影像数据在植被物候信息提取研究中的应用,提高其在生态系统研究中的应用价值。结果表明:在SOS、EOS和LOS以及TOMS几个主要的物候时间点上,MODIS NDVI时间序列的标准差较小,所得物候数据更为集中,偏离度较小,所得物候数据较稳定;而HJ-1A NDVI时间序列所得物候数据的标准差较大,数据偏离程度较大,而在POS、BOS和AOS等表征植被生命周期中生长幅度数据上,其标准差较小,离散程度小。     

青藏高原1982~2015年FPAR时空变化分析

植被吸收利用太阳光合有效辐射比率反映了植被固碳释氧能力,根据青藏高原GIMMS NDVI3g（1982~2015年）和MODIS NDVI（2001~2015年）数据,采用非线性半理论半经验模型进行FPAR反演及时空变化分析。结果表明：①2001~2015年GIMMS NDVI3g和MODIS NDVI反演FPAR在空间分布上具有较高的一致性,相关系数为0.82（P<0.01）,年际变化趋势一致至少6年的区域占80%;②青藏高原FPAR受坡度和海拔影响较大,其中15~35坡度FPAR变化最快,700~2 100 m海拔区间FPAR值最大;不同坡向对应的FPAR除南坡方向偏低外其他方向差异不大。③1982~2015年青藏高原四季FPAR时空变化研究中,冬季FPAR年际变化最明显,约78.5%的区域表现为增长趋势;秋季FPAR下降区域最多,但超过71.5%区域变化不显著;④基于MODIS NDVI和GIMMS NDVI两数据反演的所有植被类型的FPAR都在2012年间出现小幅度下降趋势,且不同植被类型FPAR的年际变化趋势各不相同。     

Logistic函数方法拟合多时序NDVI数据的改进研究

单Logistic函数曲线拟合法是NDVI时间序列重建及物候遥感中关键物候期划分的重要方法之一。虽然该方法不需要设定阈值或经验系数、较适应于不同环境区域的物候遥感监测,但是在山区NDVI噪音较大的情形下,其拟合精度仍会受到较大影响。选取秦岭中部山区为研究区,在分析了多年NDVI时间序列数据特征基础上,利用山区NDVI数据序列最大值相对于最小值更为稳定的特性,对传统单Logistic模型求解方法进行改进,采用更为稳定的参数构建模型以提高NDVI时间序列重建的精度。基于秦岭样区MODIS NDVI多年遥感数据,分别在保持植被生长季特征能力和保留高质量原始真值程度两方面对原方法与改进方法的计算结果进行比较。研究表明改进的方法在上述两个方面都具有更好的效果。在植被指数噪音较大的山区,改进的方法对NDVI时间序列重建表现出更好的适用性,可为复杂的山区物候遥感相关研究提供参考。     

基于MODIS影像的植被覆盖度提取研究

利用线性光谱混合模型对河南省三门峡地区MODIS1B影像进行植被覆盖度（fv）信息提取,将结果与反映植被覆盖度的NDVI比较,并提出在实测资料缺乏的情况下利用同期高分辨率ETM＋图像对选取样本区域进行定量验证的方法。结果表明,对于MODIS数据,线性光谱混合（LSMM）分解方法能有效提取大区域范围的植被覆盖度信息,比NDVI-fv经验统计方法更具有理论意义,为快速、准确、高效的植被监测提供了新思路。     

基于植被物候特征与监督分类的青南高原信息提取

针对大尺度区域的植被信息提取,由于范围广阔、地形复杂、气候迥异,分类精度的提高是个亟待解决的问题;通过对青南高原采用分区处理,利用植被指数的特性,将基于时间序列的NDVI数据所反映的植被物候知识,辅助信息DEM和GIS数据加入监督分类系统,进行植被信息提取,并进行了分类精度评价。研究结果表明,利用该方法对青南高原的3个地区分类后,其分类精度都达到了83.3%以上,达到了较好的分类结果。在监督分类的训练区选取过程中,将植被物候特征作为知识,结合目视解译和DEM辅助知识帮助选取训练区的方法,同时参考GIS土地利用数据,使得训练区的选取更准确可靠,可进一步提高分类精度。     

Land-use land-cover classification analysis of Giba catchment using hyper temporal MODIS NDVI satellite images

Landsat-based land-use land-cover (LULC) mapping studies were previously conducted in Giba catchment, comprising an area of 4019 km². No attempt has been done to map LULC of this catchment through the analysis of Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) time-series data. This article is aimed to see whether time-series MODIS NDVI data set is applicable for LULC mapping of Giba catchment or not. MODIS NDVI data sets of the year 2010 were used for classification analysis. The original data were subjected to MODIS Reproduction Tool and stacking. The re-projected and stacked images were filtered using Harmonic Analysis of Time-Series filtering algorism to remove the effects of cloud and other noises. The MODIS NDVI data sets (16-day maximum value composite) were classified using the ISODATA clustering algorithm available under ERDAS IMAGINE software. A series of unsupervised classification runs were carried out with a pre-defined number of classes (5–24). From this classification, the optimal numbers of classes were determined to be eight after checking for average divergence analysis. The classification result became eight LULC classes namely: bare land, grass land, irrigated land, cultivated land, area closure, shrub land, bush land, and forest land with an overall accuracy of 87.7%. It was therefore concluded that MODIS NDVI time-series image is applicable for mapping large watersheds.     

MODIS time-series imagery for forest disturbance detection and quantification of patch size effects

The Moderate Resolution Imaging Spectroradiometer (MODIS) 250 m single day surface reflectance (MOD09GQK) and 16-day composite gridded vegetation index data (MOD13Q1) were used to detect forest harvest disturbance between 2000 and 2004 in northern Maine. A MODIS multi-date Normalized Difference Vegetation Index (NDVI) forest change detection map was developed from each MODIS data set. A Landsat TM/ETM+ change detection map was developed as a reference to assess the effect of disturbed forest patch size on classification accuracy (agreement) and disturbed area estimates of MODIS. The MODIS single day and 16-day composite data showed no significant difference in overall classification accuracies. However, the 16-day NDVI change detection map had marginally higher overall classification accuracy (at 85%), but had significantly lower detection accuracy related to disturbed patch size than the single day NDVI change detection map. The 16-day composite NDVI data achieved 69% detection accuracy and the single day NDVI achieved 76% when the disturbed patch size was greater than 20 ha. The detection accuracy increased to approximately 90% for both data sets when the patch size exceeded 50 ha. The R² (range 0.6 to 0.9) and slope (range 0.5 to 0.9) of regression lines between Landsat and MODIS data (based on forest disturbance percent of township) increased with the mean disturbed patch size of each township. The 95% confidence intervals of forest disturbance percent estimate for each township were narrow with less than 1% of each township at the mean MODIS forest disturbance level.     

A comparison of MODIS 250-m EVI and NDVI data for crop mapping: a case study for southwest Kansas

Multi-temporal vegetation index (VI) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are becoming widely used for large-area crop classification. Most crop-mapping studies have applied enhanced vegetation index (EVI) data from MODIS instead of the more traditional normalized difference vegetation index (NDVI) data because of atmospheric and background corrections incorporated into EVI's calculation and the index's sensitivity over high biomass areas. However, the actual differences in the classification results using EVI versus NDVI have not been thoroughly explored. This study evaluated time-series MODIS 250-m EVI and NDVI for crop-related land use/land cover (LULC) classification in the US Central Great Plains. EVI- and NDVI-derived maps classifying general crop types, summer crop types and irrigated/non-irrigated crops were produced for southwest Kansas. Qualitative and quantitative assessments were conducted to determine the thematic accuracy of the maps and summarize their classification differences. For the three crop maps, MODIS EVI and NDVI data produced equivalent classification results. High thematic accuracies were achieved with both indices (generally ranging from 85% to 90%) and classified cropping patterns were consistent with those reported for the study area (> 0.95 correlation between the classified and USDA-reported crop areas). Differences in thematic accuracy (< 3% difference), spatially depicted patterns (> 90% pixel-level thematic agreement) and classified crop areas between the series of EVI- and NDVI-derived maps were negligible. Most thematic disagreements were restricted to single pixels or small clumps of pixels in transitional areas between cover types. Analysis of MODIS composite period usage in the classification models also revealed that both VIs performed equally well when periods from a specific growing season phase (green, peak or senescence) were heavily utilized to generate a specific crop map.     

Hierarchical Land Cover Classification based on MODIS NDVI Time-series

In this paper,we mainly used MODIS NDVI time-series dataset at 16-days temporal resolution and 250-meters spatial resolution to analyze land cover mapping of northeastern China.We used two different filter methods to fit NDVI time-series dataset,and compared their average classes’ separability based on Jeffries-Matusita distance index.In addition,we made use of hierarchical classification method to complete classification,combined with short-wave infrared spectral reflectance data and DEM.We conformed to the principle that separate area hierarchically into several parts first and then classify each part further,and use a single characteristic band first and then multiple feature bands.In the process of classification,we adopted threshold value method,support vector machine,artificial net neural and C5.0 decision tree classification to distinguish each land-cover type hierarchically.Finally,we evaluated the accuracy of the final classification of study area using known land-cover classification data and high-resolution remote sensing imagery,overall accuracy is 84.61%,Kappa coefficient is 0.8262.     

Estimation of wheat production by the integration of MODIS and ground data

This paper presents a methodology capable of combining Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) imagery and ancillary data to estimate durum wheat production in Tuscany (Central Italy). First, the phenological stages of winter wheat are simulated by the use of an agro-meteorological model (Syrius 4.1). Next, MODIS NDVI images at 250 m spatial resolution are utilized to identify fields likely grown with winter wheat. The multitemporal NDVI profiles of these fields are then integrated with Syrius 4.1 outputs through a previously developed procedure in order to simulate wheat above-ground biomass and harvest index. This allows the computation of wheat yield, which, combined with relevant cropped area, provides provincial wheat production estimates. The methodology is tested using ground and MODIS data taken over four Tuscany provinces where winter wheat is widely cultivated. The accuracy of all estimated variables (wheat cropped area, yield and production) is finally evaluated against provincial statistical data. The results of this experiment indicate that the accuracy of wheat cropped area estimation and yield simulation is variable, but interannual production variations are reproduced well for all provinces.     

基于MODIS NDVI时间序列数据的藏北草地类型识别

选取西藏自治区申扎县北部部分区域作为藏北草原的典型代表区进行研究。对区域内2009年全年23期MODIS NDVI数据收集处理,得到23个波段的时间序列数据。通过Savitzky-Golay滤波算法优化时间序列数据并对不同类型的时间序列曲线做对比分析,提取分类的特征波段。最后采用非监督分类法把研究区分为沼泽草甸、高寒草甸、中盖度草原、低盖度草原、裸地和水体6种类型,并结合地形数据分析了其中4种草地类型的空间分布特征。一系列的处理结果及分析内容证明了MODIS NDVI时间序列数据在藏北草地类型识别方面的可行性和适用性。     

Boosted decision tree classifications of land cover over Turkey integrating MODIS,climate and topographic data

This study investigates the impact of using different combinations of Moderate Resolution Imaging Spectroradiometer (MODIS) and ancillary datasets on overall and per-class classification accuracies for nine land cover types modified from the classification system of the International Geosphere Biosphere Programme (IGBP). Twelve land cover maps were generated for Turkey using boosted decision trees (BDTs) based on the stepwise addition of 14 explanatory variables derived from a time series of 16-day MODIS composites between 2000 and 2006 (Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and four spectral bands) and ancillary climate and topographic data (minimum and maximum air temperature, precipitation, potential evapotranspiration, aspect, elevation, distance to sea and slope) at 500-m resolution. Evaluation of the 12 BDTs indicated that the BDT built as a function of all the MODIS and climate variables, aspect and elevation produced the highest degree of overall classification accuracy (79.8%) and kappa statistic (0.76) followed by the BDTs that additionally included distance to sea (DtS), and both DtS and slope. Based on an independent validation dataset derived from a pre-existing national forest map and Landsat images of Turkey, the highest overall accuracy (64.7%) and kappa coefficient (0.58) among the 12 land cover maps was achieved by using MODIS-derived NDVI time series only, followed by NDVI and EVI time series combined; NDVI, EVI and four MODIS spectral bands; and the combination of all MODIS and climate data, aspect, elevation and distance to sea, respectively. The largest improvements in producer's accuracies were observed for grasslands (+50%), barrenlands (+46%) and mixed forests (+39%) and in user's accuracies for grasslands (+53%), shrublands (+30%) and mixed forests (+28%), in relation to the lowest producer's accuracy. The results of this study indicate that BDTs can increase the accuracy of land cover classifications at the national scale.     

基于MODIS遥感数据的宏观土地覆盖特征分类方法与精度分析研究

针对宏观土地覆盖遥感分类的现状,充分利用MODIS相对于AVHRR数据具有的多光谱和分辨率优势,提出了利用MODIS数据进行分类特征选择与提取并结合多时相特征进行宏观土地覆盖分类的分类方法,并在中国山东省进行了分类试验,得出以下结论:①不同比例下的训练样本与验证样本影响着总体分类精度;②从MODIS数据中得到的植被指数EVI、白天地表温度Tday、水体指数NDWI、纹理特征局部平稳Homogeneity等可以作为分类特征配合参与到多波段地表反射率Ref1-7遥感影像中,能明显提高分类精度,而土壤亮度指数NDSI则没有贡献;③提取的分类特征对总体分类精度贡献大小为:EVI贡献最大,提高近6个百分点,其次是Homogeneity、NDWI,均提高近4个百分点,而最少的Tday也贡献了近3个百分点;④各分类特征对不同地物类别具有不同的分离度,在提高某些类别的分离性时,有可能降低了其它类别的分离性。试验结果表明:在没有其它非遥感信息的前提下,仅利用MODIS遥感自身信息对宏观土地覆盖分类就可达到较高精度。     

基于MODIS遥感数据的宏观土地覆盖特征分类方法与精度分析研究

针对宏观土地覆盖遥感分类的现状,充分利用MODIS相对于AVHRR数据具有的多光谱和分辨率优势,提出了利用MODIS数据进行分类特征选择与提取并结合多时相特征进行宏观土地覆盖分类的分类方法,并在中国山东省进行了分类试验,得出以下结论:①不同比例下的训练样本与验证样本影响着总体分类精度;②从MODIS数据中得到的植被指数EVI、白天地表温度Tday、水体指数NDWI、纹理特征局部平稳Homogeneity等可以作为分类特征配合参与到多波段地表反射率Ref1-7遥感影像中,能明显提高分类精度,而土壤亮度指数NDSI则没有贡献;③提取的分类特征对总体分类精度贡献大小为:EVI贡献最大,提高近6个百分点,其次是Homogeneity、NDWI,均提高近4个百分点,而最少的Tday也贡献了近3个百分点;④各分类特征对不同地物类别具有不同的分离度,在提高某些类别的分离性时,有可能降低了其它类别的分离性。试验结果表明:在没有其它非遥感信息的前提下,仅利用MODIS遥感自身信息对宏观土地覆盖分类就可达到较高精度。     

Pattern validation for MODIS image mining of burned area objects

An image mining method was applied to Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data to estimate the area burned by forest fires occurring in Galicia (Spain) between 4 August and 15 August 2006. Five different inputs were considered: post-fire near-infrared reflectance (NIR) band, post-fire Normalized Difference Vegetation Index (NDVI) image, pre-fire and post-fire NDVI difference image and 4-μm and 11-μm thermal bands. The proposed image mining method consists of three steps: a pre-classification step, applying kernel smoothing, based on the fast Fourier transform (FFT), a modelling step applying Gaussian distributions on individual grid cells with deviating values, and a thresholding step classifying the model into burned and unburned classes. Polygons collected in the field with Global Positioning System (GPS) measurements from a helicopter permitted validation of the burned area estimation. A Z-test based on the κ statistic compared the accuracy of this estimation with the accuracies achieved by traditional methods based both on spectral changes in reflectance after the fire and active fire detection. The results showed a significant improvement (7.5%) in the accuracy of the burned area estimation after kernel smoothing. Burned area estimation based on the smoothed difference image between pre-fire and post-fire NDVI image had the highest accuracy (κ = 0.72). We conclude that the image mining algorithm successfully extracted burned area objects and that extraction was best when smoothing was applied prior to classification. Image mining methods based on using the κ statistic thus provide a valuable validation procedure when selecting the optimal MODIS input image for estimating burned area objects.