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

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

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

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

全球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图像自动化辐射校正中关键问题的解决方法和精度评价结果,以及大规模的数据处理活动引发的一些思考。     

Improved maize cultivated area estimation over a large scale combining MODIS–EVI time series data and crop phenological information

The accurate and timely information of crop area is vital for crop production and food security. In this study, the Enhanced Vegetation Index (EVI) data from MODerate resolution Imaging Spectroradiometer (MODIS) integrated crop phenological information was used to estimate the maize cultivated area over a large scale in Northeast China. The fine spatial resolution China’s Environment Satellite (HJ-1 satellite) images and the support vector machine (SVM) algorithm were employed to discriminate distribution of maize in the reference area. The mean MODIS–EVI time series curve of maize was extracted in the reference area by using multiple periods MODIS–EVI data. By analysing the temporal shift of crop calendars from northern to southern parts in Northeast China, the lag value was derived from phenological data of twenty-one agro-meteorological stations; here integrating with the mean MODIS–EVI time series image of maize, a standard MODIS–EVI time series image of maize was obtained in the whole study area. By calculating mean absolute distances (MAD) map between standard MODIS–EVI image and mean MODIS–EVI time series images, and setting appropriate thresholds in three provinces, the maize cultivated area was extracted in Northeast China. The results showed that the overall classification accuracy of maize cultivated area was approximately 79%. At the county level, the MODIS-derived maize cultivated area and statistical data were well correlated (R² = 0.82, RMSE = 283.98) over whole Northeast China. It demonstrated that MODIS–EVI time series data integrated with crop phenological information can be used to improve the extraction accuracy of crop cultivated area over a large scale.     

Detecting winter wheat phenology with SPOT-VEGETATION data in the North China Plain

Monitoring phenological change in agricultural land improves our understanding of the adaptation of crops to a warmer climate. Winter wheat–maize and winter wheat–cotton double-cropping are practised in most agricultural areas in the North China Plain. A curve-fitting method is presented to derive winter wheat phenology from SPOT-VEGETATION S10 normalized difference vegetation index (NDVI) data products. The method uses a double-Gaussian model to extract two phenological metrics, the start of season (SOS) and the time of maximum NDVI (MAXT). The results are compared with phenological records at local agrometeorological stations. The SOS and MAXT have close agreement with in situ observations of the jointing date and milk-in-kernel date respectively. The phenological metrics detected show spatial variations that are consistent with known phenological characteristics. This study indicates that time-series analysis with satellite data could be an effective tool for monitoring the phenology of crops and its spatial distribution in a large agricultural region.     

Accurate mapping of forest types using dense seasonal Landsat time-series

An accurate map of forest types is important for proper usage and management of forestry resources. Medium resolution satellite images (e.g., Landsat) have been widely used for forest type mapping because they are able to cover large areas more efficiently than the traditional forest inventory. However, the results of a detailed forest type classification based on these images are still not satisfactory. To improve forest mapping accuracy, this study proposed an operational method to get detailed forest types from dense Landsat time-series incorporating with or without topographic information provided by DEM. This method integrated a feature selection and a training-sample-adding procedure into a hierarchical classification framework. The proposed method has been tested in Vinton County of southeastern Ohio. The detailed forest types include pine forest, oak forest, and mixed-mesophytic forest. The proposed method was trained and validated using ground samples from field plots. The three forest types were classified with an overall accuracy of 90.52% using dense Landsat time-series, while topographic information can only slightly improve the accuracy to 92.63%. Moreover, the comparison between results of using Landsat time-series and a single image reveals that time-series data can largely improve the accuracy of forest type mapping, indicating the importance of phenological information contained in multi-seasonal images for discriminating different forest types. Thanks to zero cost of all input remotely sensed datasets and ease of implementation, this approach has the potential to be applied to map forest types at regional or global scales.     

Comparison and assessment of NDVI time series for seasonal wetland classification

Satellite-based wetland mapping faces challenges due to the high spatial heterogeneity and dynamic characteristics of seasonal wetlands. Although normalized difference vegetation index (NDVI) time series (NTS) shows great potential in land cover mapping and crop classification, the effectiveness of various NTS with different spatial and temporal resolution has not been evaluated for seasonal wetland classification. To address this issue, we conducted comparisons of those NTS, including the moderate-resolution imaging spectroradiometer (MODIS) NTS with 500?m resolution, NTS fused with MODIS and Landsat data (MOD_LC8-NTS), and HJ-1 NDVI compositions (HJ-1-NTS) with finer resolution, for wetland classification of Poyang Lake. Results showed the following: (1) the NTS with finer resolution was more effective in the classification of seasonal wetlands than that of the MODIS-NTS with 500-m resolution and (2) generally, the HJ-1-NTS performed better than that of the fused NTS, with an overall accuracy of 88.12% for HJ-1-NTS and 83.09% for the MOD_LC8-NTS. Future work should focus on the construction of satellite image time series oriented to highly dynamic characteristics of seasonal wetlands. This study will provide useful guidance for seasonal wetland classification, and benefit the improvements of spatiotemporal fusion models.     

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

大范围、长时间和高精度农作物空间分布基础农业科学数据的准确获取对资源、环境、生态、气候变化和国家粮食安全等问题研究具有重要现实意义和科学意义。本文针对传统阈值法农作物识别过程中阈值设置存在灵巧性差和自动化程度低等弱点,以中国粮食主产区黄淮海平原内河北省衡水市景县为典型实验区,首次将全局优化算法应用于阈值模型中阈值优化选取,开展了利用全局优化算法改进基于阈值检测的农作物分布制图方法创新研究。以冬小麦为研究对象,国产高分一号(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全局优化算法支持下阈值模型参数优化作物分布制图方法的有效性和可行性,可获得高精度冬小麦作物空间分布制图结果,这对提高中国冬小麦空间分布制图精度和自动化水平具有一定意义,也可为农作物面积农业统计数据降尺度恢复重建和大范围区域作物空间分布制图研究提供一定技术参考。     

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.     

Using high resolution QuickBird imagery for crop identification and area estimation

Imagery from recently launched high spatial resolution satellite sensors offers new opportunities for crop assessment and monitoring. A 2.8-m multispectral QuickBird image covering an intensively cropped area in south Texas was evaluated for crop identification and area estimation. Three reduced-resolution images with pixel sizes of 11.2 m, 19.6 m, and 30.8 m were also generated from the original image to simulate coarser resolution imagery from other satellite systems. Supervised classification techniques were used to classify the original image and the three aggregated images into five crop classes (grain sorghum, cotton, citrus, sugarcane, and melons) and five non-crop cover types (mixed herbaceous species, mixed brush, water bodies, wet areas, and dry soil/roads). The five non-crop classes in the 10-category classification maps were then merged as one class. The classification maps were filtered to remove the small inclusions of other classes within the dominant class. For accuracy assessment of the classification maps, crop fields were ground verified and field boundaries were digitized from the original image to determine reference field areas for the five crops. Overall accuracy for the unfiltered 2.8-m, 11.2-m, 19.6-m, and 30.8-m classification maps were 71.4, 76.9, 77.1, and 78.0%, respectively, while overall accuracy for the respective filtered classification maps were 83.6, 82.3, 79.8, and 78.5%. Although increase in pixel size improved overall accuracy for the unfiltered classification maps, the filtered 2.8-m classification map provided the best overall accuracy. Percentage area estimates based on the filtered 2.8-m classification map (34.3, 16.4, 2.3, 2.2, 8.0, and 36.8% for grain sorghum, cotton, citrus, sugarcane, melons, and non-crop, respectively) agreed well with estimates from the digitized polygon map (35.0, 17.9, 2.4, 2.1, 8.0, and 34.6% for the respective categories). These results indicate that QuickBird imagery can be a useful data source for identifying crop types and estimating crop areas.     

北京地区HJ-1卫星CCD数据的气溶胶反演及在大气校正中的应用

与现有的大气卫星传感器相比,环境一号卫星(HJ-1)CCD相机具有较高的空间分辨率(30m),使得在城市地区找到光谱纯像元的机率大大增加。本文提出一种基于纯像元提取的城市地区气溶胶光学厚度(AerosolOpticalDepth,AOD)反演算法,利用像元纯净指数来提取CCD影像上的纯像元,并由HJ-1A星和B星的多时相CCD观测数据结合地表双向反射率模型确定纯像元的地表反射特性,在此基础上反演AOD。与AERONET地基测量数据的对比表明,该算法具有较高精度,相关系数为0.83,线性拟合斜率为1.091,截距为0.053。基于该方法的AOD反演结果作为输入,能较大程度提高HJ-1卫星CCD影像大气校正的精度。     

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%。     

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.     

Examining the Phenological Cycle of Barley (Hordeum vulgare) Using Satellite and in situ Spectroradiometer Measurements for the Detection of Buried Archaeological Remains

The detection of buried archaeological remains using satellite remote sensing is still an open question in archaeological research. This research investigates how the phenological stages of crops can be used support the detection of buried archaeological remains. Ground remote sensing data using the GER-1500 spectroradiometer were obtained from two sites. One site was the Neolithic settlements in central Greece and the other was in Alampra village in Cyprus. For the latter, an archaeological environment was simulated and ground spectroradiometric measurements were systematically acquired over the different phases of the phenological cycle of barley crops. The acquired in situ reflectance measurements have been converted to "in-band" reflectance values of the Landsat TM/ETM+ using the satellite relative spectral responses filters (RSR). Based on the proposed methodology, 97 Landsat MSS, TM, and ETM+ satellite images were acquired (covering a period from 1983 to 2011), for the Thessalian (Greek) site. It has been found that phenological-cycle observations can provide valuable information for identifying buried archaeological remains. Such observations may be used in cases where the spatial resolution of satellite imagery is not high and therefore cannot help support the detection of archaeological remains using standard interpretation techniques.     

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.     

融合Landsat ETM+和MODIS数据估算高时空分辨率地表短波反照率

提出一种通过融合高空间低时间分辨率、低空间高时间分辨率地表短波反照率,来估算高时空分辨率地表短波反照率的方法。首先,利用Landsat ETM+数据,通过窄波段到宽波段的转换得到一景或多景空间分辨率较高的ETM+蓝天空短波反照率;然后,在MODIS短波反照率产品基础上,以天空光比例因子为权重,得到空间分辨率较低的MODIS蓝天空短波反照率;最后,利用STARFM(Spatial and Temporal Adaptive Reflectance Fusion Model)模型融合ETM+短波反照率的空间变化信息和MODIS短波反照率的时间变化信息,得到高时空分辨率的地表短波反照率。针对STARFM模型在异质性区域估算精度降低的问题,通过以MODIS反照率影像各像元的端元(各地类)反照率取代MODIS像元反照率来提取时空变化等信息参与STARFM模型的融合过程,达到提高异质性区域估算精度的目的。结果显示,直接利用STARFM模型估算得到的高空间分辨率地表短波反照率处在合理的精度范围内(RMSE0.02),用改进后的STARFM模型估算得到的异质性区域短波反照率和真实ETM+短波反照率间的相关系数增大。     

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.     

基于时间序列叶面积指数稀疏表示的作物种植区域提取

以华北平原黄河以北地区为研究区域,以时间序列叶面积指数LAI(Leaf Area Index)傅里叶变换的谐波特征作为不同作物识别的数据源,利用稀疏表示的分类方法识别2007年—2016年冬小麦、春玉米、夏玉米等主要农作物种植区域。首先利用上包络线Savitzky-Golay滤波分别对2007年—2016年的时间序列MODIS LAI曲线进行重构,进而对重构的年时间序列LAI进行傅里叶变换,以0—5级谐波振幅、1—5级谐波相位作为作物识别的依据,基于各类地物的训练样本,通过在线字典学习算法构建稀疏表示方法的判别字典,对每个待测样本利用正交匹配追踪算法求解稀疏系数,从而计算对应于各类地物的重构误差,根据最小重构误差判定待测样本的作物类型,并对作物识别结果的位置精度进行验证。结果表明,2007年—2016年作物识别的总体精度为77.97%,Kappa系数为0.74,表明本文提出的方法可以用于研究区域主要作物种植区域的提取。     

资源三号卫星成像在轨几何定标的探元指向角法

从资源三号卫星（ZY-3）影像严格几何模型出发,通过对星敏感器坐标系下各CCD探元指向角的分析,提出一种ZY-3成像在轨几何定标的探元指向角法。经对嵩山和洛阳两个试验区     

资源一号02D高光谱影像内陆水体叶绿素a浓度反演

2019-09-12成功发射的资源一号02D卫星(ZY-102D)搭载了新一代可见短波红外高光谱相机AHSI(Advanced Hyperspectral Imager),其丰富的细分波段和较高的空间分辨率在内陆湖库水质监测方面具有较大潜力,但数据可用性有待分析和验证.本研究以中国华东和华北平原的典型富营养湖库(太湖、...