Phenological metrics-based crop classification using HJ-1 CCD images and Landsat 8 imagery

Crop type data are an important piece of information for many applications in agriculture. Extracting crop type using remote sensing is not easy because multiple crops are usually planted into small parcels with limited availability of satellite images due to weather conditions. In this research, we aim at producing crop maps for areas with abundant rainfall and small-sized parcels by making full use of Landsat 8 and HJ-1 charge-coupled device (CCD) data. We masked out non-vegetation areas by using Landsat 8 images and then extracted a crop map from a long-term time-series of HJ-1 CCD satellite images acquired at 30-m spatial resolution and two-day temporal resolution. To increase accuracy, four key phenological metrics of crops were extracted from time-series Normalized Difference Vegetation Index curves plotted from the HJ-1 CCD images. These phenological metrics were used to further identify each of the crop types with less, but easier to access, ancillary field survey data. We used crop area data from the Jingzhou statistical yearbook and 5.8-m spatial resolution ZY-3 satellite images to perform an accuracy assessment. The results show that our classification accuracy was 92% when compared with the highly accurate but limited ZY-3 images and matched up to 80% to the statistical crop areas.     

基于HJ-1A/B卫星CCD数据的江汉平原农田物候探测

利用HJ-1A/B卫星CCD数据,提取2013—2015年三年江汉平原农田的归一化植被指数NDVI,构建时间序列曲线,利用小波变换对HJ-1A/B卫星所得的NDVI数据进行平滑降噪处理,结合地面调研资料,提取江汉平原农作物的物候信息。研究结果表明,HJ-1A/B卫星可用于农田物候监测,对于小区域的农田作物长势监测具有独特的优势。     

Landsat8和MODIS融合构建高时空分辨率数据识别秋粮作物

本文利用Wu等人提出的遥感数据时空融合方法 STDFA(Spatial Temporal Data Fusion Approach)以Landsat 8和MODIS为数据源构建高时间、空间分辨率的遥感影像数据。以此为基础,构建15种30 m分辨率分类数据集,然后利用支持向量机SVM(Support Vector Machine)进行秋粮作物识别,验证不同维度分类数据集进行秋粮作物识别的适用性。实验结果显示,不同分类数据集的秋粮作物分类结果均达到了较高的识别精度。综合各项精度指标分析,Red+Phenology数据组合对秋粮识别效果最好,水稻识别的制图精度和用户精度分别达到91.76%和82.49%,玉米识别的制图精度和用户精度分别达到85.80%和74.97%,水稻和玉米识别的总体精度达到86.90%。     

鄱阳湖HJ-1A/B卫星CCD影像大气校正研究

以准同步的Terra/MODIS反演的气溶胶为辅助,采用FLAASH模型对2009-10-24鄱阳湖HJ-1A/B卫星CCD影像进行大气校正处理。结果表明,大气影响可以被有效去除,在水体遥感反射率较高的红、绿波段,大气校正精度较高,平均相对误差分别为13.4%和9.8%;而在水体遥感反射率较低的近红外、蓝波段,大气校正精度较低,这可能与波段不同的信噪比和陆地邻近像元效应有关。     

Water body mapping method with HJ-1A/B satellite imagery

This paper proposes an integrated water body mapping method with HJ-1A/B satellite imagery, the CCD (charge coupled device) data of the Chinese environmental satellites that were launched on September 6th, 2008. It combines the difference between NDVI and NDWI (NDVI–NDWI) with SLOPE and near-infrared (NIR) band. The NDVI–NDWI index is used to enhance the contrast between water bodies and the surrounding surface features; the topographic SLOPE is used to eliminate the mountain shadow; and the NIR band is used to reduce the effects of artificial construction land. The objectives are evaluating the potential of the HJ-1A/B imagery on water body monitoring, and proposing ideally mapping method. The test study results indicated that the NDVI–NDWI index is superior to the single index of NDVI and NDWI to enhance the contrast between water bodies and the rest of the features. On the basis of the accurately mapped water bodies in the HJ-1A/B CCD images of the study area, we conclude that the HJ-1A/B multi-spectral satellite images is an ideal data source for high spatial and temporal resolution water bodies monitoring. And the integrated water body mapping method is suitable for the applications of HJ-1A/B multi-spectral satellite images in this field.     

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.     

Generation of high spatial and temporal resolution NDVI and its application in crop biomass estimation

Abstract

While data like HJ-1 CCD images have advantageous spatial characteristics for describing crop properties, the temporal resolution of the data is rather low, which can be easily made worse by cloud contamination. In contrast, although Moderate Resolution Imaging Spectroradiometer (MODIS) can only achieve a spatial resolution of 250 m in its normalised difference vegetation index (NDVI) product, it has a high temporal resolution, covering the Earth up to multiple times per day. To combine the high spatial resolution and high temporal resolution of different data sources, a new method (Spatial and Temporal Adaptive Vegetation index Fusion Model [STAVFM]) for blending NDVI of different spatial and temporal resolutions to produce high spatial–temporal resolution NDVI datasets was developed based on Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). STAVFM defines a time window according to the temporal variation of crops, takes crop phenophase into consideration and improves the temporal weighting algorithm. The result showed that the new method can combine the temporal information of MODIS NDVI and spatial difference information of HJ-1 CCD NDVI to generate an NDVI dataset with both high spatial and high temporal resolution. An application of the generated NDVI dataset in crop biomass estimation was provided. An average absolute error of 17.2% was achieved. The estimated winter wheat biomass correlated well with observed biomass (R ² of 0.876). We conclude that the new dataset will improve the application of crop biomass estimation by describing the crop biomass accumulation in detail. There is potential to apply the approach in many other studies, including crop production estimation, crop growth monitoring and agricultural ecosystem carbon cycle research, which will contribute to the implementation of Digital Earth by describing land surface processes in detail.     

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.     

全球30m分辨率多光谱影像数据自动化处理的实践与思考

全球地表覆盖遥感制图与关键技术研究项目要求对两个基准年度(2000年、2010年)全球30 m分辨率的多光谱遥感数据进行辐射处理和几何精纠正处理,为地表覆盖制图完成数据准备。数据以Landsat TM/ETM+为主,HJ-1A/B CCD数据为补充,共计2万多景影像需要进行辐射处理,有1000多景HJ-1A/B CCD影像需要几何精纠正。如此大规模的数据处理,自动化处理是必然的选择。本文介绍了HJ-1A/B CCD图像几何精纠正自动化实现中关键问题的解决方法和精度评价结果,Landsat TM/ETM+和HJ-1A/B CCD图像自动化辐射校正中关键问题的解决方法和精度评价结果,以及大规模的数据处理活动引发的一些思考。     

遥感影像时—空融合的“点”—“线”—“面”质量评价

遥感影像时—空融合可集成多源数据高空间分辨率和高时间分辨率互补优势,生成时间连续的高空间分辨率影像,在遥感影像的动态监测与时序分析等方面具有重要应用价值。然而,现有多数研究往往基于单一数据产品对时—空融合算法进行评价,而在实际生产应用中,需要验证算法在多种遥感产品数据的融合表现;此外,目前研究大多基于"单点时刻"进行评价,忽略了时—空融合在"时间线"上的有效验证。本文提出遥感影像时—空融合的"点"—"线"—"面"多角度综合质量评价策略,基于Landsat TM和MODIS影像,建立了时—空融合系列数据集,包括地表反射率、植被指数和地表温度,并在此基础上从单时相("点")、时间序列("线")、多种数据产品("面")多个角度对4种典型融合算法进行定性和定量的综合评价。结果表明,基于不同产品类型的数据集更能充分验证算法性能,且结合单点时刻和时间序列的评价更加客观。     

统计数据总量约束下全局优化阈值的冬小麦分布制图

大范围、长时间和高精度农作物空间分布基础农业科学数据的准确获取对资源、环境、生态、气候变化和国家粮食安全等问题研究具有重要现实意义和科学意义。本文针对传统阈值法农作物识别过程中阈值设置存在灵巧性差和自动化程度低等弱点,以中国粮食主产区黄淮海平原内河北省衡水市景县为典型实验区,首次将全局优化算法应用于阈值模型中阈值优化选取,开展了利用全局优化算法改进基于阈值检测的农作物分布制图方法创新研究。以冬小麦为研究对象,国产高分一号(GF-1)为主要遥感数据源,在作物面积统计数据为总量控制参考标准和全局参数优化的复合型混合演化算法SCE-UA (Shuffled Complex Evolution-University of Arizona)支持下,提出利用时序NDVI数据开展阈值模型阈值参数自动优化的冬小麦空间分布制图方法。最终,获得实验区冬小麦阈值模型最优参数,并利用优化后的阈值参数对冬小麦空间分布进行提取。通过地面验证表明,利用本研究所提方法获取的冬小麦识别结果分类精度均达到较高水平。其中冬小麦识别结果总量精度达到了99.99%,证明本研究所提阈值模型参数优化方法冬小麦提取分类结果总量控制效果良好;同时,与传统的阈值法、最大似然和支持向量机等分类方法相比,本研究所提阈值模型参数优化法区域冬小麦作物分类总体精度和Kappa系数分别都有所提高,其中,总体精度分别提高4.55%、2.43%和0.15%,Kappa系数分别提高0.12、0.06和0.01,这体现出SCE-UA全局优化算法对提高阈值模型冬小麦空间分布识别精度具有一定优势。以上研究结果证明了利用本研究所提基于作物面积统计数据总量控制以及SCE-UA全局优化算法支持下阈值模型参数优化作物分布制图方法的有效性和可行性,可获得高精度冬小麦作物空间分布制图结果,这对提高中国冬小麦空间分布制图精度和自动化水平具有一定意义,也可为农作物面积农业统计数据降尺度恢复重建和大范围区域作物空间分布制图研究提供一定技术参考。     

Crop classification and acreage estimation in North Korea using phenology features

In North Korea, reliable and timely information on crop acreage and spatial distribution is hard to obtain. In this study, we developed a fast and robust method to estimate crop acreage in North Korea using time-series normalized difference vegetation index (NDVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) data. We proposed a method to identify crop type based on NDVI phenology features using data collected in other areas with similar agri-environmental conditions to mitigate the shortage of ground truth data. Eventually the classification map (MODIScrop) was assessed using the Food and Agriculture Organization (FAO) statistical data and high-resolution crop classification maps derived from one Landsat scene (LScrop). The Pareto boundary method was used to assess the accuracy and crop distribution of the MODIScrop maps. Results showed that acreage derived from the MODIScrop maps was generally consistent with that reported in the FAO data (a relative error <4.1% for rice and <6.1% for maize, and <9.0% for soybean except for in 2004, 2008, and 2009) and the maps derived from the LScrop (a relative error about 5% in 2013, and 7% in 2008 and 2014). The classification accuracy reached 74.4%, 69.8%, and 73.1% of the areas covered by the Landsat images in 2008, 2013, and 2014, respectively. This indicates that features derived from NDVI profiles were able to characterize major crops, and the approaches developed in this study are feasible for crop mapping and acreage estimation in regions with limited ground truth data.     

Mapping Rice Paddies in Complex Landscapes with Convolutional Neural Networks and Phenological Metrics

ABSTRACT

Rice mapping with remote sensing imagery provides an alternative means for estimating crop-yield and performing land management due to the large geographical coverage and low cost of remotely sensed data. Rice mapping in Southern China, however, is very difficult as rice paddies are patchy and fragmented, reflecting the undulating and varied topography. In addition, abandoned lands widely exist in Southern China due to rapid urbanization. Abandoned lands are easily confused with paddy fields, thereby degrading the classification accuracy of rice paddies in such complex landscape regions. To address this problem, the present study proposes an innovative method for rice mapping through combining a convolutional neural network (CNN) model and a decision tree (DT) method with phenological metrics. First, a pre-trained LeNet-5 Model using the UC Merced Dataset was developed to classify the cropland class from other land cover types, i.e. built-up, rivers, forests. Then, paddy rice field was separated from abandoned land in the cropland class using a DT model with phenological metrics derived from the time-series data of the normalized difference vegetation index (NDVI). The accuracy of the proposed classification methods was compared with three other classification techniques, namely, back propagation neural network (BPNN), original CNN, pre-trained CNN applied to HJ-1 A/B charge-coupled device (CCD) images of Zhuzhou City, Hunan Province, China. Results suggest that the proposed method achieved an overall accuracy of 93.56%, much higher than those of other methods. This indicates that the proposed method can efficiently accommodate the challenges of rice mapping in regions with complex landscapes.     

融合时间序列环境卫星数据与物候特征的水稻种植区提取

获取高精度的区域水稻种植面积对于农业规划、配置与决策具有重要意义。区域尺度的水稻面积获取依赖于高时空分辨率影像,但受卫星回访周期和气候影响,难以获取足够时间序列的高时空分辨率影像,从而影响水稻种植面积遥感提取的精度。为此,提出适应于中国南方多雨云天气地区,基于国产环境卫星(HJ-1A/1B)与MODIS融合数据的水稻种植面积提取的新方法。以洞庭湖区为实验区,利用STARFM模型融合环境卫星NDVI数据与MODIS13Q1数据,获取时间序列的环境卫星NDVI数据,利用水稻关键期的NDVI数据结合物候特征参数对水稻种植区域进行提取。结果表明,该方法能有效提取区域水稻种植的面积,水稻种植面积提取的总体精度与Kappa系数分别达到91.71%与0.9024,分类结果明显优于仅采用多光谱影像或NDVI数据。该研究为中国南方多雨云天气地区水稻种植面积提取提供了有效的方法。     

基于国产遥感数据的秭归县地质灾害区域危险性评价

以位于三峡库区的秭归县为研究区,选择国产2 m全色高分辨率遥感数据进行地质灾害范围信息提取,选用国产HJ-1A/B CCD数据提取NDVI、土地利用分类、土壤含水量等动态地质灾害致灾因子,结合坡度、坡向、降水量、沟壑密度等静态因子,采用单因素叠加的信息量模型进行了区域地质灾害危险性评价,对利用国产遥感数据进行区域地质灾害信息提取和评价进行了有效探索。     

Assessing the effects of land use spatial structure on urban heat islands using HJ-1B remote sensing imagery in Wuhan,China

Urban heat islands (UHIs) have attracted attention around the world because they profoundly affect biological diversity and human life. Assessing the effects of the spatial structure of land use on UHIs is essential to better understanding and improving the ecological consequences of urbanization. This paper presents the radius fractal dimension to quantify the spatial variation of different land use types around the hot centers. By integrating remote sensing images from the newly launched HJ-1B satellite system, vegetation indexes, landscape metrics and fractal dimension, the effects of land use patterns on the urban thermal environment in Wuhan were comprehensively explored. The vegetation indexes and landscape metrics of the HJ-1B and other remote sensing satellites were compared and analyzed to validate the performance of the HJ-1B. The results have showed that land surface temperature (LST) is negatively related to only positive normalized difference vegetation index (NDVI) but to Fv across the entire range of values, which indicates that fractional vegetation (Fv) is an appropriate predictor of LST more than NDVI in forest areas. Furthermore, the mean LST is highly correlated with four class-based metrics and three landscape-based metrics, which suggests that the landscape composition and the spatial configuration both influence UHIs. All of them demonstrate that the HJ-1B satellite has a comparable capacity for UHI studies as other commonly used remote sensing satellites. The results of the fractal analysis show that the density of built-up areas sharply decreases from the hot centers to the edges of these areas, while the densities of water, forest and cropland increase. These relationships reveal that water, like forest and cropland, has a significant effect in mitigating UHIs in Wuhan due to its large spatial extent and homogeneous spatial distribution. These findings not only confirm the applicability and effectiveness of the HJ-1B satellite system for studying UHIs but also reveal the impacts of the spatial structure of land use on UHIs, which is helpful for improving the planning and management of the urban environment.     

Scrutinising MODIS and GIMMS Vegetation Indices for Extracting Growth Rhythm of Natural Vegetation in India

Satellite derived vegetation vigour has been successfully used for various environmental modeling since 1972. However, extraction of reliable annual growth information about natural vegetation (i.e., phenology) has been of recent interest due to their important role in many global models and free availability of time-series satellite data. In this study, usability of Moderate Resolution Imaging Spectro-radiometer (MODIS) and Global Inventory Modelling and Mapping Studies (GIMMS) based products in extracting phenology information about evergreen, semi-evergreen, moist deciduous and dry deciduous vegetation in India was explored. The MODIS NDVI and EVI time-series data (MOD13C1: 5.6 km spatial resolution with 16 day temporal resolution—2001 to 2010) and GIMMS NDVI time-series data(8 km spatial resolution with 15 day temporal resolution—2000 to 2006) were used. These three differently derived vegetation indices were analysed to extract and understand the vegetative growth rhythm over different regions of India. Algorithm was developed to derive onset of greenness and end of senescence automatically. The comparative analysis about differences in the results from these products was carried out. Due to dominant noise in the values of NDVI from GIMMS and MODIS during monsoon period the phenology rhythm were wrongly depicted, especially for evergreen and semi-evergreen vegetation in India. Hence, care is needed before using these data sets for understanding vegetative dynamics, biomass cestimation and carbon studies. MODIS EVI based results were truthful and comparable to ground reality. The study reveals spatio-temporal patterns of phenology, rate of greening, rate of senescence, and differences in results from these three products.     

Extracting Regional Pattern of Wheat Sowing Dates Using Multispectral and High Temporal Observations from Indian Geostationary Satellite

Monitoring of Agricultural crops using remote sensing data is an emerging tool in recent years. Spatial determination of sowing date is an important input of any crop model. Geostationary satellite has the capability to provide data at high temporal interval to monitor vegetation throughout the entire growth period. A study was conducted to estimate the sowing date of wheat crop in major wheat growing states viz. Punjab, Haryana, Uttar Pradesh (UP), Madhya Pradesh (MP), Rajasthan and Bihar. Data acquired by Charged Couple Detector (CCD) onboard Indian geostationary satellite INSAT 3A have continental (Asia) coverage at 1 km?×?1 km spatial resolution in optical spectral bands with high temporal frequency. Daily operational Normalized Difference Vegetation Index (NDVI) product from INSAT 3A CCD available through Meteorological and Oceanographic Satellite Data Archival Centre (MOSDAC) was used to estimate sowing date of wheat crop in selected six states. Daily NDVI data acquired from September 1, 2010 to December 31, 2010 were used in this study. A composite of 7 days was prepared for further analysis of temporal profile of NDVI. Spatial wheat crop map derived from AWiFS (56 m) were re-sampled at INSAT 3A CCD parent resolution and applied over each 7 day composite. The characteristic temporal profiles of 7 day NDVI composite was used to determine sowing date. NDVI profile showed decreasing trend during maturity of kharif crop, minimum value after harvest and increasing trend after emergence of wheat crop. A mathematical model was made to capture the persistent positive slope of NDVI profile after an inflection point. The change in behavior of NDVI profile was detected on the basis of change in NDVI threshold of 0.3 and sowing date was estimated for wheat crop in six states. Seven days has been deducted after it reached to threshold value with persistent positive slope to get sowing date. The clear distinction between early sowing and late sowing regions was observed in study area. Variation of sowing date was observed ranging from November 1 to December 20. The estimated sowing date was validated with the reported sowing date for the known wheat crop regions. The RMSD of 3.2 (n?=?45) has been observed for wheat sowing date. This methodology can also be applied over different crops with the availability of crop maps.     

高海拔山区草地物候提取方法对比分析

准确量测高海拔山区的植物物候对理解全球变化下的敏感生态系统的响应具有重要意义。利用物候相机和遥感技术开展物候信息的提取和对比,既能准确评估物候相机在山区植物物候提取的性能,又可为山区遥感物候数据反演提供重要参考。利用中国新疆维吾尔自治区天山山区人工观测、物候相机和遥感数据,测试了5种曲线拟合方式与4种物候参数提取方法的20种组合的物候参数提取结果,对比了3种数据在物候信息提取结果的异同。结果表明:（1）植物物候相机能在天山山区草地物候观测中提供高时间分辨率的绿度变化信息,是山区开展物候观测并验证遥感物候数据的有效手段。（2）山区雨雪天气等对相对绿度指数产生较强噪声影响,需要选择合适的滤波器进行去噪。（3）曲线拟合方式和物候提取方法均对物候参数数值产生影响。而提取方法可产生更明显的差异性,其中,阈值法和导数法提取的物候数值相近,开始期与人工观测的返青期一致性较好,停止期与枯黄期一致性较好;而Klosterman方法和Gu方法提取物候数值相近,提取的开始期与人工观测的返青末期一致性较好,停止期与人工观测的枯黄末期一致性较好。（4）20种不同滤波+提取方法的组合形式在山区遥感数据物候信息提取的有效性仅为48%,中分辨率成像光谱仪数据的最有效提取方法为Beck+Derivatives组合,可见光红外成像辐射套件数据的最优提取方法为Beck+Threshold组合和Elmore+Derivatives组合。     

HJ-1A/B卫星CCD影像的武汉市东湖水色三要素遥感研究

以武汉市东湖为研究区域,利用同步的MODIS-Terra气溶胶光学厚度数据为输入参数,采用FLAASH模型对2010年3月11日HJ-1A/B卫星CCD影像进行大气校正处理,并利用多年实测数据建立叶绿素a浓度、悬浮泥沙浓度、黄色物质吸收系数三要素神经网络反演模型,对水色三要素进行反演。通过对反演结果与实测数据的对比分析可知,悬浮泥沙浓度、黄色物质吸收系数和叶绿素a浓度的平均相对误差分别为28.052%、17.628%和35.621%,表明HJ-1A/B卫星CCD传感器基本能满足II类水体水色要素的遥感监测需求。