多熟种植是中国重要的种植制度,在农业生产中具有重要的地位。复种指数反映了耕地实际利用强度,及时获取其空间分布信息是国家农业决策的基础。论文在重新理解和定义复种指数的基础上,提出了一种基于NDVI曲线迭代修正的复种指数遥感提取方法,并以SPOT/VGT多时相NDVI数据为基础,提取了1999~2004年中国北方17省市区农用地的复种指数。结果表明:提取的复种指数空间分布符合我国耕作制度区划,基于目视解译的样点检验总体精度为95.24%,Kappa系数为0.9057,与根据统计数据得出的复种指数结果的线性回归斜率为0.9288（R2=0.9159,P<0.001）,表明该方法能够准确快捷地提取复种指数,及时地为农业和土地管理部门提供耕作制度的空间信息。
Multiple cropping system characterized by cropping index is crucial to Chinese food security.Multiple cropping index(MCI)refers to the times of sequential crop planting in the same arable land in one year,reflecting utility degree of arable land to be used at a certain period.It is desired to extract the MCI and its spatial distribution information by remotely sensed data for sustainable agricultural development.Traditionally,the MCI is calculated by statistical data at local administration unit,which is time-lagged,labor-consuming and poor in creditability as well as lack of spatial distribution.Meanwhile,remote sensing technology has been widely applied to agriculture and crop monitoring.The advancement in remote sensing technology provides potenial to obtain the actual MCI information efficiently and reliably.Many researches indicate that crop growth dynamic can be monitored by the time series of Normal Difference Vegetation Index(NDVI) data.It is found that the peak of the time series of the NDVI curve indicates that the ground biomass of crops reaches the maximum,and fluctuates with the crops growing processes such as sowing,seeding,heading,ripeness,and harvesting within one year.Thus,the MCI is defined as the number of peaks of the time series of the NDVI curve.However,since NDVI data are affected by cloud and poor atmospheric conditions,the curve of time series of NDVI turns out to be noisy with a lot of small peaks and valleys in one cycle,which makes it more difficult to extract the actual MCI.This study develops a new method for extracting cropping index based on NDVI time-series by which cloud and other contaminations can be corrected effectively.The details about the method include:(1) smoothing the NDVI time-series by Savitzky-Golay filter to get the long-term growing trend;(2) correcting the growing trend curve by iteration process to get rid of the false peaks;and(3) obtaining the MCI by counting the number of true peaks.Using 1km 10-day Maximum Value Composite SPOT/VGT NDVI time-series,the MCI of 17 provinces of northern China from 1999 to 2004 was extracted by this new method.The results revealed that: the cropping index of Huang-Huai-Hai region was higher than that of the other regions in northern China.The spatial distribution of the extracted MCI was consistent with the actual Chinese cropping system.The total precision of sample validation based on visual identification was 95.24%,the coefficient of Kappa 0.9057,and the slope of linear regression of the MCI between remotely sensed data and statistical data was 0.9288(R2=0.9159,P<0.001),suggesting that this method could provide an effective way to extracting spatial information of the MCI for agriculture and land management.
[1] 封志明,刘宝勤,杨艳昭. 中国耕地资源数量变化的趋势分析与数据重建: 1949～2003 [J]. 自然资源学报, 2005, 20 (1): 35～43. [FENG Zhi-ming, L IU Bao-qin, YANG Yan-zhao. A study of the changing trend of Chinese cultivated land amount and data reconstructing: 1949-2003. Journal of N atural Resources, 2002, 20 (1): 35-43.]
[2] 刘巽浩. 90年代我国耕作制度发展展望[J]. 耕作与栽培, 1992, (2): 1～9. [L IU Xun-hao. The view of Chinese cropping system in 1990 s. Cultivate and Plant, 1992, (2): 1-9.]
[3] 刘巽浩. 中国耕作制度[M]. 北京:农业出版社, 1993112～15. [L IU Xun-hao. Cropp ing System in China. Beijing: Ag2 ricultural Press, 1993. 12-15.]
[4] 刘巽浩. 耕作学[M]. 北京: 中国农业出版社, 199611 ～9. [L IU Xun-hao. Cultivation. Beijing: China Agricultural Press, 199611-9.]
[5] 曾祥熙. 作物栽培学[M]. 西安:陕西人民教育出版社, 1989112. [ZENG Xiang-xi. Crop Cultivation. Xi’an: Shaanxi Peop le’s Education Press, 1989112.]
[6] 马永祥. 复种指数刍议[J]. 作物杂志, 2003, (4): 541[MA Yong2xiang. Multi-cropp ing index discussion. CropM agzine, 2003, (4): 54.]
[7] Reed B C, Brown J F, VanderZee D, et al. Measuring phenological variability from satellite imagery[J]. Journal of Vegetation Science, 1994, 5 (5): 703-714.
[8] 范锦龙,吴炳方. 复种指数遥感监测方法[J]. 遥感学报, 2004, 8 (6): 628 ～636. [FAN Jin-long, WU Bing-fang. A methodology for retrieving cropp ing index from NDV I p rofile. Journal of Rem ote Sencing, 2004, 8 (6): 628-636.]
[9] 闫慧敏,曹明奎,刘纪远,等. 基于多时相遥感信息的中国农业种植制度空间格局研究[J]. 农业工程学报, 2005, 21 (4): 85～90. [YAN Hui-min, CAO Mingkui, L IU Ji-yuan, et al. Characterizing spatial patterns ofmultip le cropp ing system in China from multi2temporal remote sensing images. Transactions of the CSAE, 2005, 21 (4): 85-90.]
[10] 彭代亮,黄敬峰,金辉民. 基于MOD IS2NDV I的浙江省耕地复种指数监测[J]. 中国农业科学, 2006, 39 (7): 1352～1357. [PENG Dai-liang, HUANG Jing-feng, J IN Hui2min. The monitoring for sequential cropp ing index of arable land in Zhejiang Province usingMOD IS-NDV I. Scientia Agriculture S inica, 2006, 39 (7): 1352-1357.]
[11] 刘巽浩. 论我国耕地种植指数(复种)的潜力[J]. 作物杂志, 1997, (3): 1～3. [L IU Xun2hao. Study on potential of multicropp ing index of arable land in China. CropM agzine, 1997, (3): 1-3.]
[12] 卞新民,冯金侠. 多元多熟种植制度复种指数计算方法探讨[J]. 南京农业大学学报, 1999, 22 (1): 11～15. [B IAN Xin-min, FENG J in-xia. Method of calculating cropp ing system index of diversified multip le cropp ing systems. Journal of N anjing Agriculture University, 1999, 22 (1): 11-15.]
[13] Rikie S, Tomoyuki N, Tetsuzo Y. West2east contrast of phenology and climate in northern Asia revealed using a remotely sensed vegetation index[J]. International Journal of B iom eteorol. , 2003, 47: 126-138.
[14] Lee R, Yu F, PR ICE K P, et al. Evaluating vegetation phenological patterns in InnerMongolia usingNDV I time-seriesa-nalysis[J]. International Journal of Rem ote Sensing, 2002, 23: 2505-2512.
[15] Bondiglio A, Cuomo V, LanfrediM, et al. InterfacingNOAA /ANHRR NDV I and soil truth map s formonitoring vegetation phenology at a local scale in a heterogeneous landscape of Southern Italy[J]. International Journal of Rem ote Sensing, 2002, 23: 4181-4195.
[16] Stowe L L, McClain E P, Carey R, et al. Global distribution of cloud cover derived from NOAA /AVHRR operational satellite data[J]. Advances in Space Research, 1991, (3): 51-54.
[17] Chen J, Jonsson P, TamuraM, et al. A simp le method for reconstructing a high2quality NDV I time-series data set based on the Savitzky-Golay filter[J]. Rem ote Sensing of Environm ent, 2001, 91: 332-344.
[18] Savitzky A, GolayM J E. Smoothing and differentiation of data by simp lified least squares p rocedures [J]. Analytical Chem istry, 1964, 36: 1627-1639.
[19] Stropp ianaD, Pinnock S, Pereira JM C, et al. Radiometric analysis of SPOT2VEGETATION images for burnt area detection in Northern Australia[J]. Rem ote Sensing of Environm ent, 2002, 82: 21-37.
[20] Vegetation Program. SPOT Vegetation User Guide[EB /OL]. http://www.spotimage.fr/data/images/vege/vegetat/book_1/e_frame.htm. 1998.
[21] 闫慧敏,刘纪远,曹明奎. 近20 年中国耕地复种指数的时空变化[J]. 地理学报, 2005, 60 (4): 559 ～556. [YAN Hui-min, L IU Ji-yuan, CAO Ming-kui. Remotely sensed multip le cropp ing index variations in China during 1981-2000. Acta Geographica S inica, 2005, 60 (4): 559-556.]
[22] 沈学年,刘巽浩. 多熟种植[M]. 北京:农业出版社, 198319～15. [SHEN Xue-nian, L IU Xun-hao. More Ripeness Cultivation. Beijing: Agriculture Press, 198319-15.]
[23] 中国农业信息网. 全国农业统计提要[EB /OL]. http://www.agri.gov.cn/sjzl/nongyety.htm. [China Agriculture Information Net. The Statistic Abstract of State Agriculture. http://www.agri.gov.cn/sjzl/nongyety.htm.]
[24] 中华人民共和国农业部. 中国农业年鉴(2005) [Z]. 北京:中国农业出版社, 2005. 159. [Ministry ofAgriculture of the Peop le’s Republic of China. China Agriculture Yearbook (2005) . Beijing: China Agricultural Press, 20051159.]
