Updating Landsat-based forest cover maps with MODIS images using multiscale spectral-spatial-temporal superresolution mapping

With the high deforestation rates of global forest covers during the past decades, there is an ever-increasing need to monitor forest covers at both fine spatial and temporal resolutions. Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat series images have been used commonly for satellite-derived forest cover mapping. However, the spatial resolution of MODIS images and the temporal resolution of Landsat images are too coarse to observe forest cover at both fine spatial and temporal resolutions. In this paper, a novel multiscale spectral-spatial-temporal superresolution mapping (MSSTSRM) approach is proposed to update Landsat-based forest maps by integrating current MODIS images with the previous forest maps generated from Landsat image. Both the 240 m MODIS bands and 480 m MODIS bands were used as inputs of the spectral energy function of the MSSTSRM model. The principle of maximal spatial dependence was used as the spatial energy function to make the updated forest map spatially smooth. The temporal energy function was based on a multiscale spatial-temporal dependence model, and considers the land cover changes between the previous and current time. The novel MSSTSRM model was able to update Landsat-based forest maps more accurately, in terms of both visual and quantitative evaluation, than traditional pixel-based classification and the latest sub-pixel based super-resolution mapping methods The results demonstrate the great efficiency and potential of MSSTSRM for updating fine temporal resolution Landsat-based forest maps using MODIS images.     

Net primary productivity of forest stands in New Hampshire estimated from Landsat and MODIS satellite data

Background

A simulation model that relies on satellite observations of vegetation cover from the Landsat 7 sensor and from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate net primary productivity (NPP) of forest stands at the Bartlett Experiment Forest (BEF) in the White Mountains of New Hampshire.

Results

Net primary production (NPP) predicted from the NASA-CASA model using 30-meter resolution Landsat inputs showed variations related to both vegetation cover type and elevational effects on mean air temperatures. Overall, the highest predicted NPP from the NASA-CASA model was for deciduous forest cover at low to mid-elevation locations over the landscape. Comparison of the model-predicted annual NPP to the plot-estimated values showed a significant correlation of R² = 0.5. Stepwise addition of 30-meter resolution elevation data values explained no more than 20% of the residual variation in measured NPP patterns at BEF. Both the Landsat 7 and the 250-meter resolution MODIS derived mean annual NPP predictions for the BEF plot locations were within ± 2.5% of the mean of plot estimates for annual NPP.

Conclusion

Although MODIS imagery cannot capture the spatial details of NPP across the network of closely spaced plot locations as well as Landsat, the MODIS satellite data as inputs to the NASA-CASA model does accurately predict the average annual productivity of a site like the BEF.     

Evaluating MODIS snow products for modelling snowmelt runoff: Case study of the Rio Grande headwaters

Snow-covered area (SCA) is a key variable in the Snowmelt-Runoff Model (SRM) and in other models for simulating discharge from snowmelt. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM + ) or Operational Land Imager (OLI) provide remotely sensed data at an appropriate spatial resolution for mapping SCA in small headwater basins, but the temporal resolution of the data is low and may not always provide sufficient cloud-free dates. The coarser spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) offers better temporal resolution and in cloudy years, MODIS data offer the best alternative for mapping snow cover when finer spatial resolution data are unavailable. However, MODIS’ coarse spatial resolution (500 m) can obscure fine spatial patterning in snow cover and some MODIS products are not sensitive to end-of-season snow cover. In this study, we aimed to test MODIS snow products for use in simulating snowmelt runoff from smaller headwater basins by a) comparing maps of TM and MODIS-based SCA and b) determining how SRM streamflow simulations are changed by the different estimates of seasonal snow depletion. We compared gridded MODIS snow products (Collection 5 MOD10A1 fractional and binary SCA; SCA derived from Collection 6 MOD10A1 Normalised Difference Snow Index (NDSI) Snow Cover), and the MODIS Snow Covered-Area and Grain size retrieval (MODSCAG) canopy-corrected fractional SCA (SCA_MG), with reference SCA maps (SCA_REF) generated from binary classification of TM imagery. SCA_MG showed strong agreement with SCA_REF; excluding true negatives (where both methods agreed no snow was present) the median percent difference between SCA_REF and SCA_MG ranged between −2.4% and 4.7%. We simulated runoff for each of the four study years using SRM populated with and calibrated for snow depletion curves derived from SCA_REF. We then substituted in each of the MODIS-derived depletion curves. With efficiency coefficients ranging between 0.73 and 0.93, SRM simulation results from the SCA_MG runs yielded the best results of all the MODIS products and only slightly underestimated discharge volume (between 7 and 11% of measured annual discharge). SRM simulations that used SCA derived from Collection 6 NDSI Snow Cover also yielded promising results, with efficiency coefficients ranging between 0.73 and 0.91.In conclusion, we recommend that when simulating snowmelt runoff from small basins (<4000 km²) with SRM, we recommend that users select either canopy-corrected MODSCAG or create their own site-specific products from the Collection 6 MOD10A1 NDSI.     

Investigating relationships between Landsat-7 ETM data and spatial segregation of LULC types under shifting agriculture in southern Cameroon

This paper investigates statistical relationships between land use/land cover (LULC), Landsat-7 ETM⁺ imagery and landscape mosaic structure in southern Cameroon where the conversion of tropical rain forest to shifting cultivation leads to dynamic processes, acting on the spatial aggregation of various LULC types. A Global Positioning System (GPS) was used in the field to identify a total of 171 shifting cultivation patches representing eight LULC types in two sub-areas. Because of the lack of a cloud-free image for the date of field sampling, the ETM⁺ imagery was acquired 2 months after field survey, during which it was assumed that no significant changes in LULC occurred (all dry season). Per pixel correlations were developed between spectral reflectance data, vegetation indices and LULC. As an exploratory study, several statistical methods (analysis of variance, means separations (Tukey HSD), principal component analysis (PCA), geo-statistical analysis, image classification and landscape metrics) were applied on point data and sensor images for evaluating the spatial variability within the landscape. Most variables explained 30–72% of LULC variation in the whole dataset. Those variables with high information content of LULC (infrared bands 4, 5, 7 and derived indices and PC1) also showed long ranges (6 km) spatial dependence as compared to those varying only within 1 km range. The results of these statistical analyses suggested the need to group some LULC types and the application of the Maximum Likelihood Classifier (MLC) for supervised classification provided a LULC map with the highest accuracy (81%) after consolidation of perennial LULC types, such as bush fallow, forest fallow and cocoa plantations. Landscape metrics computed from this map showed a high level of patch diversity and connectivity within the landscape and provided input data that can further be used to simulate predictive maps as substitute to cloud-covered sensor imageries. Landsat-7 ETM⁺ imagery proved to be useful in discriminating (with about 80% accuracy) the most dynamic LULC types such cropped plots and young fallow patches (shifting every season) and the extension front of the agricultural landscape.     

Prediction of Land Cover Change Using Markov and Cellular Automata Models: Case of Al-Ain,UAE, 1992-2030

The UAE has witnessed rapid urban development and economic growth in recent years. With its ambitious vision to become one of the advanced nations by 2021, planners and policy-makers need to know the most likely direction of future urban development. In this study, remotely sensed imagery coupled with cellular automata models were used to predict land cover in Al Ain, the second largest city in the Emirate of Abu Dhabi. Markov and cellular automata models were used for 1992 and 2006 to predict land cover in 2012. Land Use and Land Cover maps for the study area were derived from 1992, 2006, and 2012 Landsat satellite images (TM, ETM⁺). The models achieved an overall accuracy of approximately 80 %. A Markov model was applied for 2006 and 2012 to predict land cover in 2030. The results conformed to the general trend of the Al Ain Master Plan 2030. This study demonstrates that remote sensing, with the availability of free Landsat data, is a viable technology that could be used to help in the prediction process especially in developing countries, where data availability is a problem.     

Forest cover classification using Landsat ETM+ data and time series MODIS NDVI data

Forest cover plays a key role in climate change by influencing the carbon stocks, the hydrological cycle and the energy balance. Forest cover information can be determined from fine-resolution data, such as Landsat Enhanced Thematic Mapper Plus (ETM+). However, forest cover classification with fine-resolution data usually uses only one temporal data because successive data acquirement is difficult. It may achieve mis-classification result without involving vegetation growth information, because different vegetation types may have the similar spectral features in the fine-resolution data. To overcome these issues, a forest cover classification method using Landsat ETM+ data appending with time series Moderate-resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) data was proposed. The objective was to investigate the potential of temporal features extracted from coarse-resolution time series vegetation index data on improving the forest cover classification accuracy using fine-resolution remote sensing data. This method firstly fused Landsat ETM+ NDVI and MODIS NDVI data to obtain time series fine-resolution NDVI data, and then the temporal features were extracted from the fused NDVI data. Finally, temporal features combined with Landsat ETM+ spectral data was used to improve forest cover classification accuracy using supervised classifier. The study in North China region confirmed that time series NDVI features had significant effects on improving forest cover classification accuracy of fine resolution remote sensing data. The NDVI features extracted from time series fused NDVI data could improve the overall classification accuracy approximately 5% from 88.99% to 93.88% compared to only using single Landsat ETM+ data.     

Evaluation of global land cover maps for cropland area estimation in the conterminous United States

Global land cover data could provide continuously updated cropland acreage and distribution information, which is essential to a wide range of applications over large geographical regions. Cropland area estimates were evaluated in the conterminous USA from four recent global land cover products: MODIS land cover (MODISLC) at 500-m resolution in 2010, GlobCover at 300-m resolution in 2009, FROM-GLC and FROM-GLC-agg at 30-m resolution based on Landsat imagery circa 2010 against the US Department of Agriculture survey data. Ratio estimators derived from the 30-m resolution Cropland Data Layer were applied to MODIS and GlobCover land cover products, which greatly improved the estimation accuracy of MODISLC by enhancing the correlation and decreasing mean deviation (MDev) and RMSE, but were less effective on GlobCover product. We found that, in the USA, the CDL adjusted MODISLC was more suitable for applications that concern about the aggregated county cropland acreage, while FROM-GLC-agg gave the least deviation from the survey at the state level. Correlation between land cover map estimates and survey estimates is significant, but stronger at the state level than at the county level. In regions where most mismatches happen at the county level, MODIS tends to underestimate, whereas MERIS and Landsat images incline to overestimate. Those uncertainties should be taken into consideration in relevant applications. Excluding interannual and seasonal effects, R² of the FROM-GLC regression model increased from 0.1 to 0.4, and the slope is much closer to one. Our analysis shows that images acquired in growing season are most suitable for Landsat-based cropland mapping in the conterminous USA.     

Time-series maps of aboveground biomass in dipterocarps forests of Malaysia from PALSAR and PALSAR-2 polarimetric data

Background

Malaysia typically suffers from frequent cloud cover, hindering spatially consistent reporting of deforestation and forest degradation, which limits the accurate reporting of carbon loss and CO₂ emissions for reducing emission from deforestation and forest degradation (REDD+) intervention. This study proposed an approach for accurate and consistent measurements of biomass carbon and CO₂ emissions using a single L-band synthetic aperture radar (SAR) sensor system. A time-series analysis of aboveground biomass (AGB) using the PALSAR and PALSAR-2 systems addressed a number of critical questions that have not been previously answered. A series of PALSAR and PALSAR-2 mosaics over the years 2007, 2008, 2009, 2010, 2015 and 2016 were used to (i) map the forest cover, (ii) quantify the rate of forest loss, (iii) establish prediction equations for AGB, (iv) quantify the changes of carbon stocks and (v) estimate CO₂ emissions (and removal) in the dipterocarps forests of Peninsular Malaysia.

Results

This study found that the annual rate of deforestation within inland forests in Peninsular Malaysia was 0.38% year^?1 and subsequently caused a carbon loss of approximately 9 million Mg C year^?1, which is equal to emissions of 33 million Mg CO₂ year^?1, within the ten-year observation period. Spatially explicit maps of AGB over the dipterocarps forests in the entire Peninsular Malaysia were produced. The RMSE associated with the AGB estimation was approximately 117 Mg ha^?1, which is equal to an error of 29.3% and thus an accuracy of approximately 70.7%.

Conclusion

The PALSAR and PALSAR-2 systems offer a great opportunity for providing consistent data acquisition, cloud-free images and wall-to-wall coverage for monitoring since at least the past decade. We recommend the proposed method and findings of this study be considered for MRV in REDD+?implementation in Malaysia.

     

Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS vegetation continuous fields with lidar-based estimates of error

Abstract

We developed a global, 30-m resolution dataset of percent tree cover by rescaling the 250-m MOderate-resolution Imaging Spectroradiometer (MODIS) Vegetation Continuous Fields (VCF) Tree Cover layer using circa- 2000 and 2005 Landsat images, incorporating the MODIS Cropland Layer to improve accuracy in agricultural areas. Resulting Landsat-based estimates maintained consistency with the MODIS VCF in both epochs (RMSE =8.6% in 2000 and 11.9% in 2005), but showed improved accuracy in agricultural areas and increased discrimination of small forest patches. Against lidar measurements, the Landsat-based estimates exhibited accuracy slightly less than that of the MODIS VCF (RMSE=16.8% for MODIS-based vs. 17.4% for Landsat-based estimates), but RMSE of Landsat estimates was 3.3 percentage points lower than that of the MODIS data in an agricultural region. The Landsat data retained the saturation artifact of the MODIS VCF at greater than or equal to 80% tree cover but showed greater potential for removal of errors through calibration to lidar, with post-calibration RMSE of 9.4% compared to 13.5% in MODIS estimates. Provided for free download at the Global Land Cover Facility (GLCF) website (www.landcover.org), the 30-m resolution GLCF tree cover dataset is the highest-resolution multi-temporal depiction of Earth's tree cover available to the Earth science community.     

Land cover dynamics monitoring with Landsat data in Kunming,China: a cost-effective sampling and modelling scheme using Google Earth imagery and random forests

Changes in forest composition impact ecological services, and are considered important factors driving global climate change. A hybrid sampling method along with a modelling approach to map current and past land cover in Kunming, China is reported. MODIS land cover (2001–2011) data-sets were used to detect pixels with no apparent change. Around 3000 ‘no change points’ were systematically selected and sampled using Google Earth’s high-resolution imagery. Thirty-five per cent of these points were verified and used for training and validation. We used Random forests to classify multi-temporal Landsat imagery. Results show that forest cover has had a net decrease of 14385?ha (1.3% of forest area), which was primary converted to shrublands (11%), urban and barren land (2.7%) and agriculture (2.5%). Our validation indicates an overall accuracy (Kappa) of 82%. Our methodology can be used to consistently map the dynamics of land cover change in similar areas with minimum costs.     

A time-integrated MODIS burn severity assessment using the multi-temporal differenced normalized burn ratio (dNBRMT)

Burn severity is an important parameter in post-fire management. It incorporates both the direct fire impact (vegetation depletion) and ecosystem responses (vegetation regeneration). From a remote sensing perspective, burn severity is traditionally estimated using Landsat's differenced normalized burn ratio (dNBR). In this case study of the large 2007 Peloponnese (Greece) wildfires, Landsat dNBR estimates correlated reasonably well with Geo composite burn index (GeoCBI) field data of severity (R² = 0.56). The usage of Landsat imagery is, however, restricted by cloud cover and image-to-image normalization constraints. Therefore a multi-temporal burn severity approach based on coarse spatial, high temporal resolution moderate resolution imaging spectroradiometer (MODIS) imagery is presented in this study. The multi-temporal dNBR (dNBR_MT) is defined as the 1-year integrated difference between burned pixels and their unique control pixels. These control pixels were selected based on time series similarity and spatial context and reflect how burned pixels would have behaved in the case no fire had occurred. Linear regression between downsampled Landsat dNBR and dNBR_MT estimates resulted in a moderate-high coefficient of determination R² = 0.54. dNBR_MT estimates are indicative for the change in vegetation productivity due to the fire. This change is considerably higher for forests than for more sparsely vegetated areas like shrub lands. Although Landsat dNBR is superior for spatial detail, MODIS-derived dNBR_MT estimates present a valuable alternative for burn severity mapping at continental to global scale without image availability constraints. This is beneficial to compare trends in burn severity across regions and time. Moreover, thanks to MODIS's repeated temporal sampling, the dNBR_MT accounts for both first- and second-order fire effects.     

European forest cover mapping with high resolution satellite data: The Carpathians case study

The aim of the study was to elaborate a methodology for forest mapping based on high resolution satellite data, relevant for reporting on forest cover and spatial pattern changes in Europe. The Carpathians were selected as a case study area and mapped using 24 Landsat scenes, processed independently with a supervised approach combining image segmentation, knowledge-based rules to extract a training set and the maximum likelihood decision rule. Validation was done with available very high resolution imagery. Overall accuracies per scene ranged from 93 to 96%. The labelling disagreement in overlapping areas of adjacent scenes was 6.8% on average.     

Estimating spatial variations of total evaporation using multispectral sensors within the uMngeni catchment,South Africa

Total evaporation is of importance in assessing and managing long-term water use, especially in water-limited environments. Therefore, there is need to account for water utilisation by different land uses for well-informed water resources management and future planning. This study investigated the feasibility of using multispectral Landsat 8 and moderate resolution imaging spectroradiometer (MODIS) remote sensing data to estimate total evaporation within the uMngeni catchment in South Africa, using surface energy balance system. The results indicated that Landsat 8 at 30 m resolution has a better spatial representation of total evaporation, when compared to the 1000 m MODIS. Specifically, Landsat 8 yielded significantly different mean total evaporation estimates for all land cover types (one-way ANOVA; F_4.964?=?87.011, p < 0.05), whereas MODIS failed to differentiate (one-way ANOVA; F_2.853?=?0.125, p = 0.998) mean total evaporation estimates for the different land cover types across the catchment. The findings of this study underscore the utility of the Landsat 8 spatial resolution and land cover characteristics in deriving accurate and reliable spatial variations of total evaporation at a catchment scale.     

Landsat 8 vs. Landsat 5: A comparison based on urban and peri-urban land cover mapping

An image dataset from the Landsat OLI spaceborne sensor is compared with the Landsat TM in order to evaluate the excellence of the new imagery in urban landcover classification. Widely known pixel-based and object-based image analysis methods have been implemented in this work like Maximum Likelihood, Support Vector Machine, k-Nearest Neighbor, Feature Analyst and Sub-pixel. Classification results from Landsat OLI provide more accurate results comparing to the Landsat TM. Object-based classifications produced a more uniform result, but suffer from the absorption of small rare classes into large homogenous areas, as a consequence of the segmentation, merging and the spatial parameters in the spatial resolution (30 m) of Landsat images. Based exclusively on the overall accuracy reports, the SVM pixel-based classification from Landsat 8 proved to be the most accurate for the purpose of mapping urban land cover, using medium spatial resolution imagery.     

Phenology from Landsat when data is scarce: Using MODIS and Dynamic Time-Warping to combine multi-year Landsat imagery to derive annual phenology curves

Green-leaf phenology describes the development of vegetation throughout a growing season and greatly affects the interaction between climate and the biosphere. Remote sensing is a valuable tool to characterize phenology over large areas but doing at fine- to medium resolution (e.g., with Landsat data) is difficult because of low numbers of cloud-free images in a single year. One way to overcome data availability limitations is to merge multi-year imagery into one time series, but this requires accounting for phenological differences among years. Here we present a new approach that employed a time series of a MODIS vegetation index data to quantify interannual differences in phenology, and Dynamic Time Warping (DTW) to re-align multi-year Landsat images to a common phenology that eliminates year-to-year phenological differences. This allowed us to estimate annual phenology curves from Landsat between 2002 and 2012 from which we extracted key phenological dates in a Monte-Carlo simulation design, including green-up (GU), start-of-season (SoS), maturity (Mat), senescence (Sen), end-of-season (EoS) and dormancy (Dorm). We tested our approach in eight locations across the United States that represented forests of different types and without signs of recent forest disturbance. We compared Landsat-based phenological transition dates to those derived from MODIS and ground-based camera data from the PhenoCam-network. The Landsat and MODIS comparison showed strong agreement. Dates of green-up, start-of-season and maturity were highly correlated (r 0.86-0.95), as were senescence and end-of-season dates (r > 0.85) and dormancy (r > 0.75). Agreement between the Landsat and PhenoCam was generally lower, but correlation coefficients still exceeded 0.8 for all dates. In addition, because of the high data density in the new Landsat time series, the confidence intervals of the estimated keydates were substantially lower than in case of MODIS and PhenoCam. Our study thus suggests that by exploiting multi-year Landsat imagery and calibrating it with MODIS data it is possible to describe green-leaf phenology at much finer spatial resolution than previously possible, highlighting the potential for fine scale phenology maps using the rich Landsat data archive over large areas.     

The improvement of an object-oriented classification using multi-temporal MODIS EVI satellite data

Abstract

This paper investigates the contribution of multi-temporal enhanced vegetation index (EVI) data to the improvement of object-based classification accuracy using multi-spectral moderate resolution imaging spectral-radiometer (MODIS) imagery. In object-oriented classification, similar pixels are firstly grouped together and then classified; the produced result does not suffer the speckled appearance and closer to human vision. EVI data are from the MODIS sensor aboard Terra spacecraft. 69 EVI data (scenes) were collected during the period of three years (2001–2003) in a mountainous vegetated area. These data sets were used to study the phenology of the land cover types. Different land cover types show distinct fluctuations over time in EVI values and this information might be used to improve object-oriented land cover classification. Two experiments were carried out: one was only with single date MODIS multispectral data, and the other one including also the 69 EVI images. Eight classes were distinguished: temperate forest, tropical dry forest, grassland, irrigated agriculture, rain-fed agriculture, orchards, lava flows and human settlement. The two classifications were evaluated with independent verification data, and the results showed that with multi-temporal EVI data, the classification accuracy was improved 5.2%. Evaluated by McNemar's test, this improved was significant, with significance level p=0.01.     

基于超分辨率重建的多时相MODIS与Landsat反射率融合方法

提出一种基于超分辨率重建的MODIS与Landsat反射率图像融合方法,以STARFM算法与超分辨率重建为基础,使用观测的MODIS和Landsat地表反射率图像预测给定时刻的Landsat合成反射率图像。该方法利用基于稀疏表示的超分辨率重建方法对MODIS图像进行分辨率增强,实验结果表明这一操作能够增加原MODIS图像的空间细节,有助于提高STARFM算法的预测精度;另一方面,考虑输入两个基时刻图像相差较大时原STARFM算法预测的反射率会存在"时间平滑"的问题,限制每次只使用一个基时刻MODIS和Landsat图像对进行STARFM预测,使用逐图像块选择策略,从由两个基时刻图像分别进行预测得到的两组预测图像中选择最优的预测,同样得到了优于STARFM算法的预测结果。     

GF-2支持下的干旱区稀疏植被区植被覆盖度估算

现有像元二分模型MODIS植被覆盖度模型因其形式简单、适用性较强的特点被广泛应用于区域植被覆盖度（FVC）的估算。然而，研究表明在沙漠和低植被覆盖的西部干旱区，从250 m的影像上很难精准地获取NDVI_veg（全植被覆盖植被指数）和NDVI_soil（全裸土区植被指数）参数。利用常用的直方图累计法获取模型所需参数NDVI_veg和NDVI_soil，估算结果存在普遍高估现象。为此，本文首先引入同期获取的GF-2号卫星数据，从GF-2号影像上提取植被覆盖像元；然后，利用Pixel Aggregate方法重采样至250 m分辨率，获取250 m空间分辨率下纯植被和纯裸土像元；最后，将纯植被和纯裸土像元各自空间位置相对应的MODIS NDVI数据最大值作为模型所需NDVI_veg和NDVI_soil参数，实现研究区内植被覆盖度的估算。试验通过与线性回归法、多项式回归法和直方图累计像元二分模型法估算结果进行精度对比，结果表明：利用GF-2影像辅助的像元二分模型，精准地获取了低植被覆盖区NDVI_veg和NDVI_soil模型参数，提高了干旱区植被覆盖度的估算精度，并有效地抑制了受稀疏植被影响NDVI在干旱区普遍偏高问题导致的FVC高估的现象。     

融合高时空分辨率数据估算植被净初级生产力

植被净初级生产力NPP(Net Primary Production)遥感估算与分析,有赖于高时空分辨率的遥感数据,但目前中高分辨率的遥感数据受卫星回访周期及天气的影响,在中国南方地区难以获取连续时间序列的数据,从而影响了高精度的区域植被净初级生产力的遥感估算。为此,提出一种基于多源遥感数据时空融合技术与CASA模型估算高时空分辨率NPP的方法。首先,利用多源遥感数据,即Landsat8 OLI数据与MODIS13Q1数据,采用遥感数据时空融合方法,获得了时间序列的Landsat8 OLI融合数据;然后,基于Landsat8 OLI时空融合数据,并采用CASA模型,以长株潭城市群核心区为例,进行区域植被NPP的遥感估算。研究结果表明,基于时间序列Landsat融合数据估算的30m分辨率的NPP具有良好的空间细节信息,且估算值与实测值的相关系数达0.825,与实测NPP数据保持了较好的一致性。     

一种高时空分辨率NDVI数据集构建方法-STAVFM

ETM NDVI可以用来在30m的尺度上开展植被的监测,然而在Landsat卫星16天的重访周期和云污染等因素的影响下,常常会在相当长的一段时间内无法获取有效的ETM NDVI数据,给这一尺度下的植被动态监测带来了一定困难。相比之下,MODIS虽然在空间上只有250m分辨率的NDVI产品,却可以每天进行相同区域的监测。针对ETM空间分辨率高和MODIS时间分辨率高的特点,本研究选择实验区,基于对STARFM方法的改进,构建不同时空分辨率NDVI的时空融合模型-STAVFM,使用该模型对ETM NDVI与MODIS NDVI融合,构建了高时空分辨率NDVI数据集。研究结果表明,通过MODIS NDVI时间变化信息与ETM NDVI空间差异信息的有机结合,实现缺失高空间分辨率NDVI的有效预测（3景预测NDVI与实际NDVI的相关系数分别达到了0.82、0.90和0.91）,从而构建高时空分辨率NDVI数据集。所构建的高时空分辨率NDVI数据集在时间上保留了高时间分辨率数据的时间变化趋势,空间上又反映了高空间分辨率数据的空间细节差异。