Combining spatiotemporal fusion and object-based image analysis for improving wetland mapping in complex and heterogeneous urban landscapes

Remote sensing has been proven promising in wetland mapping. However, conventional methods in a complex and heterogeneous urban landscape usually use mono temporal Landsat TM/ETM + images, which have great uncertainty due to the spectral similarity of different land covers, and pixel-based classifications may not meet the accuracy requirement. This paper proposes an approach that combines spatiotemporal fusion and object-based image analysis, using the spatial and temporal adaptive reflectance fusion model to generate a time series of Landsat 8 OLI images on critical dates of sedge swamp and paddy rice, and the time series of MODIS NDVI to calculate phenological parameters for identifying wetlands with an object-based method. The results of a case study indicate that different types of wetlands can be successfully identified, with 92.38%. The overall accuracy and 0.85 Kappa coefficient, and 85% and 90% for the user’s accuracies of sedge swamp and paddy respectively.     

Detection and estimation of mixed paddy rice cropping patterns with MODIS data

In this paper, we developed a more sophisticated method for detection and estimation of mixed paddy rice agriculture from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. Previous research demonstrated that MODIS data can be used to map paddy rice fields and to distinguish rice from other crops at large, continental scales with combined Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) analysis during the flooding and rice transplanting stage. Our approach improves upon this methodology by incorporating mixed rice cropping patterns that include single-season rice crops, early-season rice, and late-season rice cropping systems. A variable EVI/LSWI threshold function, calibrated to more local rice management practices, was used to recognize rice fields at the flooding stage. We developed our approach with MODIS data in Hunan Province, China, an area with significant flooded paddy rice agriculture and mixed rice cropping patterns. We further mapped the aerial coverage and distribution of early, late, and single paddy rice crops for several years from 2000 to 2007 in order to quantify temporal trends in rice crop coverage, growth and management systems. Our results were validated with finer resolution (2.5 m) Satellite Pour l’Observation de la Terre 5 High Resolution Geometric (SPOT 5 HRG) data, land-use data at the scale of 1/10,000 and with county-level rice area statistical data. The results showed that all three paddy rice crop patterns could be discriminated and their spatial distribution quantified. We show the area of single crop rice to have increased annually and almost doubling in extent from 2000 to 2007, with simultaneous, but unique declines in the extent of early and late paddy rice. These results were significantly positive correlated and consistent with agricultural statistical data at the county level (P < 0.01).     

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

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

Improved coastal wetland mapping using very-high 2-meter spatial resolution imagery

Accurate wetland maps are a fundamental requirement for land use management and for wetland restoration planning. Several wetland map products are available today; most of them based on remote sensing images, but their different data sources and mapping methods lead to substantially different estimations of wetland location and extent. We used two very high-resolution (2 m) WorldView-2 satellite images and one (30 m) Landsat 8 Operational Land Imager (OLI) image to assess wetland coverage in two coastal areas of Tampa Bay (Florida): Fort De Soto State Park and Weedon Island Preserve. An initial unsupervised classification derived from WorldView-2 was more accurate at identifying wetlands based on ground truth data collected in the field than the classification derived from Landsat 8 OLI (82% vs. 46% accuracy). The WorldView-2 data was then used to define the parameters of a simple and efficient decision tree with four nodes for a more exacting classification. The criteria for the decision tree were derived by extracting radiance spectra at 1500 separate pixels from the WorldView-2 data within field-validated regions. Results for both study areas showed high accuracy in both wetland (82% at Fort De Soto State Park, and 94% at Weedon Island Preserve) and non-wetland vegetation classes (90% and 83%, respectively). Historical, published land-use maps overestimate wetland surface cover by factors of 2–10 in the study areas. The proposed methods improve speed and efficiency of wetland map production, allow semi-annual monitoring through repeat satellite passes, and improve the accuracy and precision with which wetlands are identified.     

Testing estimation of water surface in Italian rice district from MODIS satellite data

Recent changes in rice crop management within Northern Italy rice district led to a reduction of seeding in flooding condition, which may have an impact on reservoir water management and on the animal and plant communities that depend on the flooded paddies. Therefore, monitoring and quantifying the spatial and temporal variability of water presence in paddy fields is becoming important. In this study we present a method to estimate dynamics of presence of standing water (i.e. fraction of flooded area) in rice fields using MODIS data. First, we produced high resolution water presence maps from Landsat by thresholding the Normalised Difference Flood Index (NDFI) made: we made it by comparing five Landsat 8 images with field-obtained information about rice field status and water presence. Using these data we developed an empirical model to estimate the flooding fraction of each MODIS cell. Finally we validated the MODIS-based flooding maps with both Landsat and ground information. Results showed a good predictability of water surface from Landsat (OA = 92%) and a robust usability of MODIS data to predict water fraction (R² = 0.73, EF = 0.57, RMSE = 0.13 at 1 × 1 km resolution). Analysis showed that the predictive ability of the model decreases with the greening up of rice, so we used NDVI to automatically discriminate estimations for inaccurate cells in order to provide the water maps with a reliability flag. Results demonstrate that it is possible to monitor water dynamics in rice paddies using moderate resolution multispectral satellite data. The achievement is a proof of concept for the analysis of MODIS archives to investigate irrigation dynamics in the last 15 years to retrieve information for ecological and hydrological studies.     

Mapping Wetland Areas Using Landsat-Derived NDVI and LSWI: A Case Study of West Songnen Plain,Northeast China

Increasing interest in wetlands for environmental management requires an understanding of the location, spatial extent, and configuration of the resource. The National Wetlands Inventory is the most commonly used data source for this information. However, its accuracy is limited in some contexts, such as agricultural and forested wetlands. An large number of studies have mapped wetlands worldwide from the perspective of land use and land cover change. However, information on the actual wetland planting areas annually is limited, which greatly impacts ongoing research. In this case study of the West Songnen Plain, we developed a simple algorithm for the quick mapping of wetlands by utilizing their unique physical features, such as annual display of phenological land-cover change of exposed soils, shallow flooding water, and plants from multi-temporal Landsat images. Temporal variations of the Normalized Difference Vegetation Index (NDVI) and Land Surface Water Index (LSWI) derived from Landsat images in 2010 for wetlands at different growth stages were analyzed. Results show that during the ante-tillering phase, the NDVI value (above zero) is lower than the LSWI value of paddies because of flooding of shallow water; during the reproductive and ripening phases, the NDVI value is higher than the LSWI value (above zero); and during the post-harvest wetland planting phase, the NDVI value is still higher than the LSWI value, but the LSWI value is negative. Wetland areas can be detected using one or two images in the optimum time window. The algorithm based on the difference of NDVI and LSWI values derived from Landsat images was used to extract the actual wetland planting area. Validated alongside statistical data, the algorithm showed high accuracy. Therefore, this algorithm highlights the unique features of wetlands and can help in mapping the actual wetland area annually on a regional scale. Results further indicate that the new method has a classification accuracy of 92 %. In comparison, two traditional methods based on Landsat-7/ETM registered accuracy rates of only 83 % and 87 % respectively.     

Rice crop phenology mapping at high spatial and temporal resolution using downscaled MODIS time-series

Satellite data holds considerable potential as a source of information on rice crop growth which can be used to inform agronomy. However, given the typical field sizes in many rice-growing countries such as China, data from coarse spatial resolution satellite systems such as the Moderate Resolution Imaging Spectroradiometer (MODIS) are inadequate for resolving crop growth variability at the field scale. Nevertheless, systems such as MODIS do provide images with sufficient frequency to be able to capture the detail of rice crop growth trajectories throughout a growing season. In order to generate high spatial and temporal resolution data suitable for mapping rice crop phenology, this study fused MODIS data with lower frequency, higher spatial resolution Landsat data. An overall workflow was developed which began with image preprocessing, calculation of multi-temporal normalized difference vegetation index (NDVI) images, and spatiotemporal fusion of data from the two sensors. The Spatial and Temporal Adaptive Reflectance Fusion Model was used to effectively downscale the MODIS data to deliver a time-series of 30 m spatial resolution NDVI data at 8-day intervals throughout the rice-growing season. Zonal statistical analysis was used to extract NDVI time-series for individual fields and signal filtering was applied to the time-series to generate rice phenology curves. The downscaled MODIS NDVI products were able to characterize the development of paddy rice at fine spatial and temporal resolutions, across wide spatial extents over multiple growing seasons. These data permitted the extraction of key crop seasonality parameters that quantified inter-annual growth variability for a whole agricultural region and enabled mapping of the variability in crop performance between and within fields. Hence, this approach can provide rice crop growth data that is suitable for informing agronomic policy and practice across a wide range of scales.     

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.     

Comprehensive study of the biophysical parameters of agricultural crops based on assessing Landsat 8 OLI and Landsat 7 ETM+ vegetation indices

For three agricultural crop types, winter wheat (Triticum aestivum L.), barley (Hordeum vulgare L.), and canola (Brassica napus L.), we estimated biophysical parameters including fresh and dry biomass, leaf area index (LAI), and vegetation water content, for which we found the equivalent water thickness (EWT), fuel moisture content per fresh weight (FMC_FW), and fuel moisture content per dry weight (FMC_DW). We performed these estimations using data from the newly launched Landsat 8 Operational Land Imager (OLI) sensor, as well as its predecessor the Landsat 7 Enhanced Thematic Mapper Plus (ETM+). Progress in the design of the new sensor (i.e., Landsat 8), including narrower near-infrared (NIR) wavebands, higher signal-to-noise ratio (SNR), and greater radiometric resolution highlights the necessity to investigate the biophysical parameters of agricultural crops, especially compared to data from its predecessor. This study aims to evaluate vegetation indices (VIs) derived from the Landsat 8 OLI and the Landsat 7 ETM+. Both the Landsat 8 OLI and Landsat 7 ETM+ VIs agreed well with in-situ data measurements. However, the Landsat 8 OLI-derived VIs were generally more consistent with in situ data than the Landsat 7 ETM+ VIs. We also note that the Landsat 8 OLI is better able to capture the small variability of the VIs because of its higher SNR and wider radiometric range; in addition, the saturation phenomenon occurred earlier for the Landsat 7 ETM+ than for the Landsat 8 OLI. This indicates that the new sensor is better able to estimate the biophysical parameters of crops.     

Land Cover and Vegetation Change in the Yellow River Delta Nature Reserve Analyzed with Landsat Thematic Mapper Data

Abstract

Coastal wetland is a major part of wetlands in the world. Land cover and vegetation mapping in a deltaic lowland environment is complicated by the rapid and significant changes of geomorphic forms. Remote sensing provides an important tool for coastal land cover classification and landscape analysis. The study site in this paper is the Yellow River Delta Nature Reserve (YRDNR) at the Yellow River mouth in Shangdong province, China. Yellow River Delta is one of the fastest growing deltas in the world. YRDNR was listed as a national level nature reserve in 1992. The objectives of this paper are two fold: to study the land cover status of YRDNR, and to examine the land cover change since it was declared as a nature reserve. Land cover and vegetation mapping in YRDNR was developed using multi‐spectral Landsat Thematic Mapper (TM) imagery acquired in 1995. Land cover and landscape characteristics were analyzed with the help of ancillary GIS. Land use investigation data in 1991 were used for comparison with Landsat classification map. Our results show that YRDNR has experienced significant landscape change and environmental improvement after 1992.     

Landsat系列卫星对地观测40年回顾及LDCM前瞻

Landsat系列卫星数据凭借其长期连续、全球覆盖、适中的时间空间分辨率和科学的数据存档与分发策略等优势,逐渐成为地表特征和地球系统科学研究中最有效的遥感数据之一,并广泛应用于生态环境、农林地矿、能源资源、教育科研和政府管理等领域。而第8代陆地卫星--陆地卫星数据连续任务卫星(LDCM)于2013年2月发射升空,该卫星携带了运行性陆地成像仪(OLI)和热红外传感器(TIRS)两种传感器。与Landsat 7/ETM+相比,OLI/TIRS在波段设置、辐射分辨性能和扫描方式上都得到很大改进,其中OLI共包括9个波段,新增海岸带(coastal)监测和卷云(cirrus)识别波段,TIRS则设置了两个热红外波段。如果LDCM能够成功升空运行,它将继续承担起长期连续对地观测的使命。     

Land cover classification using Landsat 8 Operational Land Imager data in Beijing,China

The successful launch of Landsat 8 provides a new data source for monitoring land cover, which has the potential to significantly improve the characterization of the earth’s surface. To assess data performance, Landsat 8 Operational Land Imager (OLI) data were first compared with Landsat 7 ETM + data using texture features as the indicators. Furthermore, the OLI data were investigated for land cover classification using the maximum likelihood and support vector machine classifiers in Beijing. The results indicated that (1) the OLI data quality was slightly better than the ETM + data quality in the visible bands, especially the near-infrared band of OLI the data, which had a clear improvement; clear improvement was not founded in the shortwave-infrared bands. Moreover, (2) OLI data had a satisfactory performance in terms of land cover classification. In summary, OLI data were a reliable data source for monitoring land cover and provided the continuity in the Landsat earth observation.     

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.     

Monitoring Mangrove forest landcover changes in the coastline of Bangladesh from 1976 to 2015

Abstract

This study used multi-date Landsat images to quantify mangrove cover changes in the whole of Bangladesh from 1976 to 2015. Images were pre-processed with an atmospheric correction using Dark Object Subtraction (DOS) and Relative Radiometric Normalization (RRN) using Pseudo-Invariant Features (PIFs). Land Use/Land Cover (LU/LC) classification map was generated using Maximum Likelihood (MaxLike) algorithm, indicating the areal extent of mangroves increased by 3.1% between 1976 and 2015, where 1.79% of this increase occurred between 2000 and 2015. Though mangrove areas remained almost constant in the Sundarbans, Chakaria Sundarbans has almost disappeared between 1976 and 1989. The overall accuracy of Landsat MSS, TM, ETM+, and L8 OLI classified images were 80%, 80%, 87%, and 97% respectively. The study also found deforestation, shrimp & salt farm, coastal erosion and sedimentation, and mangrove plantation could be responsible for mangrove changes in Bangladesh.     

利用Landsat 8 OLI进行汉江下游水体浊度反演

综合比较Landsat 8 OLI遥感影像与地面同步水质监测结果,发现Landsat 8的近红外波段与其他波段的组合和水体浊度具有较高的相关性,以此为基础运用OLI的第1、3、5波段组合建立了汉江中下游浊度的遥感反演数学模型。根据该模型生成了2013年4~11月共3幅汉江中下游浊度分布图,并进行了空间分析。精度验证表明,模型相对误差在15%左右,R2=0.71,表明运用Landsat 8 OLI可有效地监测该区域水体浊度分布情况。     

面向对象方法的时间序列MODIS数据湿地信息提取——以洞庭湖流域为例

以洞庭湖流域为研究区,对大范围湿地信息遥感提取方法进行了研究。先基于时间序列MODIS EVI及物候特征参数,通过J-M(Jeffries-Matusita distance)距离分析,构建了MODIS(250 m)最佳时序组合分类数据;其次,通过Johnson指数确定了最佳分割尺度,采用面向对象的遥感分类方法(Random tree分类器)提取了洞庭湖流域的湿地信息,并验证该方法的适用性。研究结果表明,基于时序数据与面向对象的Random tree分类的总体精度和Kappa系数分别为78.84%和0.71,较之基于像元的相同算法的总体分类精度和Kappa系数分别提高了5.79%和0.04。同时,基于面向对象方法的湿地整体的用户精度与生产者精度较基于像元方法分别提高了4.56%和6.21%,可有效提高大区域湿地信息提取的精度。     

Evaluating metrics derived from Landsat 8 OLI imagery to map crop cover

Developing techniques are required to generate agricultural land cover maps to monitor agricultural fields. Landsat 8 Operational Land Imager (OLI) offers reflectance data over the visible to shortwave-infrared range. OLI offers several advantages, such as adequate spatial and spectral resolution, and 16 day repeat coverage, furthermore, spectral indices derived from Landsat 8 OLI possess great potential for evaluating the status of vegetation. Additionally, classification algorithms are essential for generating accurate maps. Recently, multi-Grained Cascade Forest, which is also called deep forest, was proposed, and it was shown to give highly competitive performance for classification. However, the ability of this algorithm to generate crop maps with satellite data had not yet been evaluated. In this study, the reflectance at 7 bands and 57 spectral indices calculated from Landsat 8 OLI data were evaluated for its potential for crop type identification.     

Production of global land cover data – GLCNMO

Abstract

Global land cover is one of the fundamental contents of Digital Earth. The Global Mapping project coordinated by the International Steering Committee for Global Mapping has produced a 1-km global land cover dataset – Global Land Cover by National Mapping Organizations. It has 20 land cover classes defined using the Land Cover Classification System. Of them, 14 classes were derived using supervised classification. The remaining six were classified independently: urban, tree open, mangrove, wetland, snow/ice, and water. Primary source data of this land cover mapping were eight periods of 16-day composite 7-band 1-km MODIS data of 2003. Training data for supervised classification were collected using Landsat images, MODIS NDVI seasonal change patterns, Google Earth, Virtual Earth, existing regional maps, and expert's comments. The overall accuracy is 76.5% and the overall accuracy with the weight of the mapped area coverage is 81.2%. The data are available from the Global Mapping project website (http://www.iscgm.org/). The MODIS data used, land cover training data, and a list of existing regional maps are also available from the CEReS website. This mapping attempt demonstrates that training/validation data accumulation from different mapping projects must be promoted to support future global land cover mapping.     

Primitives as building blocks for constructing land cover maps

Land cover maps play an integral role in environmental management. However, countries and institutes encounter many challenges with producing timely, efficient, and temporally harmonized updates to their land cover maps. To address these issues we present a modular Regional Land Cover Monitoring System (RLCMS) architecture that is easily customized to create land cover products using primitive map layers. Primitive map layers are a suite of biophysical and end member maps, with land cover primitives representing the raw information needed to make decisions in a dichotomous key for land cover classification. We present best practices to create and assemble primitives from optical satellite using computing technologies, decision tree logic and Monte Carlo simulations to integrate their uncertainties. The concept is presented in the context of a regional land cover map based on a shared regional typology with 18 land cover classes agreed on by stakeholders from Cambodia, Laos PDR, Myanmar, Thailand, and Vietnam. We created annual map and uncertainty layers for the period 2000–2017. We found an overall accuracy of 94% when taking uncertainties into account. RLCMS produces consistent time series products using free long term historical Landsat and MODIS data. The customizable architecture can include a variety of sensors and machine learning algorithms to create primitives and the best suited smoothing can be applied on a primitive level. The system is transferable to all regions around the globe because of its use of publicly available global data (Landsat and MODIS) and easily adaptable architecture that allows for the incorporation of a customizable assembly logic to map different land cover typologies based on the user's landscape monitoring objectives     

Land cover mapping applications with MODIS: a literature review

Abstract

Land use/land cover monitoring and mapping is crucial to efficient management of the land and its resources. Since the late 1980s increased attention has been paid to the use of coarse resolution optical data. The Moderate Resolution Imaging Spectroradiometer (MODIS) has features, which make it particularly suitable to earth characterization purposes. MODIS has 10 products dedicated mainly to land cover characterization and provides three kinds of data: angular, spectral and temporal. MODIS data also includes information about the data quality through the ‘Quality Assessment’ product. In this paper, we review how MODIS data are used to map land cover including the preferred MODIS products, the preprocessing and classification approaches, the accuracy assessment, and the results obtained.