A Shape-matching Cropping Index (CI) Mapping Method to Determine Agricultural Cropland Intensities in China using MODISTime-series Data

摘要

As an important indicator of cropping intensity, Cropping Index (ci) is defined as the number of peaks in the Vegetation Index (vi) time-series curve in a year. The existing CI mapping algorithms (e.g., the cross-fitting and the second order difference algorithm) are vulnerable to noise contained in the vi time series and need a priori knowledge and some extra constraints which could not be directly derived from the vi time-series data. In this paper, a shape-matching method is developed which can map ci directly from the preprocessed vi time-series data without the de-noising processes. This shape-matching method utilizes a decision-making process to find out the true peaks in the vi time-series curve based on a rank order mathematical morphology algorithm. The processing procedure involves five steps: (a) determination of the temporal moving window size, (b) detection of local maximum/minimum points, (c) exclusion of false maximum/minimum points, (d) determination of the threshold for the minimumgrowth amplitude, and (e) mapping of ci. This shape-matching method only needs two input parameters, the temporal moving window size and the threshold for the minimum growth amplitude, which can be both directly derived from the vi time-series data withsome selected test pixels. Moreover, the response of the shape-matching method is relatively insensitive to the exact values of its design parameters, making it more flexible and effective in adapting to other regions. This new method is applied to mapthe Cis in Jiangsu Province, China, in 2010, based on the Enhanced Vegetation Index (evi) time-series data derived from the Moderate Resolution Imaging Spectroradiometer (modis) product. The overall ci mapping accuracy for the shape-matching method is 80percent, which is much higher than the ci mapping accuracy of 60 percent for the second order difference algorithm. This shape-matching method can be further applied to other regions with a grid-search for its optimal parameters using some test pixels.