基于多时相中分辨率时间过程特征的江汉平原棉花种植面积变化监测

doi:10.3724/SP.J.1047.2016.01141

摘要/Abstract

摘要：

棉花是中国重要的经济作物,快速、准确地提取棉花的种植面积和分布信息,对于优化棉花种植空间格局、科学指导棉花生产及提高其管理水平具有十分重要的意义。为了探讨多时相中高分辨率影像在棉花种植面积监测中的可行性,本文以江汉平原为研究区,根据棉花物候特征,选取2012年、2014年江汉平原棉花生长关键期的多时相HJ-1A/1B卫星数据,通过分析研究区棉花不同生育期的光谱特征和归一化植被指数（NDVI）时序变化特征,对分类影像进行阈值分割、掩膜处理,最后利用决策树算法提取研究区2012年、2014年棉花种植面积。通过计算混淆矩阵评价分类精度的方法和提取面积精度方法对棉花提取结果进行评价,总体精度达到95.96%,Kappa系数为0.93,以农业局统计数据为参考,2012年、2014年HJ数据提取的棉花种植面积精度分别达到了97.91%、91.27%。因此,在不受云和降水等因素的影响下,基于江汉平原区域关键时相HJ卫星CCD影像数据,可利用该方法进行棉花种植面积监测。

关键词: HJ卫星, 棉花, 种植面积, 江汉平原

Abstract:

Cotton is an important economic crop in our country, and also it is one of the main sources of income to the farmers of Jianghan plain in Hubei Province. Therefore, it is very important to obtain the information about cotton planting area, and accurately and timely get the spatial distribution for optimizing the cotton planting spatial pattern, scientifically guiding the cotton production and improving the management level. To explore the feasibility of middle and high resolution remote sensing images in monitoring the cotton plantation area, we took Jianghan plain as the study area in this paper. According to the growth period duration and the phonological characteristics of cotton, we selected the multi-temporal HJ-1A/1B images during the critical period of cotton growth of Jianghan plain in 2012 and 2014 as the data source. By analyzing the spectral property and normalized difference vegetation index (NDVI) time series variation during the cotton growing period in the study area, we use the threshold configurations and the mask processes to build the decision trees for the estimation of cotton planting area. Finally, the planting areas of cotton were extracted using the decision trees method in 2012 and 2014. With respect to the confusion matrix calculations, the overall accuracy of the study area reached 95.96% and the Kappa coefficient was 0.93. Then, by taking the statistics data provided by the Agricultural Bureau as a reference, the accuracies of area extraction results from HJ data are 97.91% and 91.27% in 2012 and 2014, respectively. Results indicated that this method can accurately reflect the distribution of the cotton area in Jianghan plain. In addition, it is concluded that the cotton planting area can be accurately extracted as long as we have the HJ satellite CCD data for the critical period of cotton growth and the image data is not affected by factors such as clouds and precipitation.

Key words: HJ satellite, cotton, planting area, Jianghan plain

尤慧, 高华东, 苏荣瑞, 刘凯文, 肖玮钰. 基于多时相中分辨率时间过程特征的江汉平原棉花种植面积变化监测[J]. 地球信息科学学报, 2016, 18(8): 1141-1149.DOI:10.3724/SP.J.1047.2016.01141

YOU Hui,GAO Huadong,SU Rongrui,LIU Kaiwen,XIAO Weiyu. Monitoring the Changes of Cotton Plantation Area Based on the Multi-temporal Middle Resolution Features of Temporal Process in Jianghan Plain[J]. Journal of Geo-information Science, 2016, 18(8): 1141-1149.DOI:10.3724/SP.J.1047.2016.01141

图/表 9

表1

图1

图2

图3

表2

图4

表3

表4

图5

参考文献 18

[1]	张成,宋述元,关业虎,等.江汉平原棉区棉花早衰原因调查与预防对策[J].湖北农业科学,2007,46(3):371-373. doi: 10.3969/j.issn.0439-8114.2007.03.017
	[ Zhang C, Song S Y, Guan Y H, et al.Investigation to the reason of premature senility of the cotton and its preventing counter measures[J]. Hubei Agricultural Sciences, 2007,46(3):371-373. ] doi: 10.3969/j.issn.0439-8114.2007.03.017
[2]	马众模,娄径,裴志远,等.长江流域棉花种植面积变化遥感监测——以安庆市和池州市为例[J].合肥工业大学学报:自然科学版,2005,28(7):738-742. doi: 10.3969/j.issn.1003-5060.2005.07.008
	[ Ma Z M, Lou J, Pei Z Y, et al.Monitoring of the change of cotton plantation area in the Changjiang valley by remote sensing-taking Anqing city and Chizhou city as examples[J]. Journal of Hefei University of Technology (Natural Science), 2005,28(7):738-742. ] doi: 10.3969/j.issn.1003-5060.2005.07.008
[3]	湖北省统计局,2014年湖北省国民经济和社会发展统计公报[EB/OL].湖北省人民政府门户网站,2015-02-25.
	[ Bureau of Statistics in Hubei Provincial.2014 national economic and social development statistics bulletin of Hubei Province[J]. People's Government Portal of Hubei Province, 2015-02-25. ]
[4]	刘珺,田庆久,黄彦,等.利用多时相HJ卫星CCD遥感影像提取嘉祥县秋收作物[J].遥感信息,2012,27(2):67-70. doi: 10.3969/j.issn.1000-3177.2012.02.013
	[ Liu J, Tian Q J, Huang Y, et al.Crop extraction based on multi-temporal HJ satellite CCD data in Jiaxiang county[J]. Remote Sensing Information, 2012,27(2):67-70. ] doi: 10.3969/j.issn.1000-3177.2012.02.013
[5]	苏荣瑞,熊勤学,耿一风,等.利用多时相HJ-CCD影像监测江汉平原南部地区棉花和中稻种植面积[J].长江流域资源与环境,2013,22(11):32-36.
	[ Su R R, Xiong Q X, Geng Y F, et al.Cotton and mid-season rice area extraction based on multi-temporal HJ satellite CCD data in southern Jianghan plain[J]. Resources and Environment in the Yangtze Basin, 2013,22(11):32-36. ]
[6]	许文波,田亦陈.作物种植面积遥感提取方法的研究进展[J].云南农业大学学报,2005,20(1):94-98. doi: 10.3969/j.issn.1004-390X.2005.01.020
	[ Xu W B, TianY C. Overview of extraction of crop area from remote sensing[J]. Journal of Yunnan Agricultural University, 2005,20(1):94-98. ] doi: 10.3969/j.issn.1004-390X.2005.01.020
[7]	潘耀忠,李乐,张锦水,等.基于典型物候特征的MODIS-EVI时间序列数据农作物种植面积提取方法—小区域冬小麦实验研究[J].遥感学报,2011,15(3):578-594. doi: 10.11834/jrs.20110066
	[ Pan Y Z, Li L, Zhang J S, et al.Crop area estimation based on MODIS-EVI time series according to distinct characteristics of key phenology phases: a case study of winter wheat area estimation in area[J]. Journal of Remote Sensing, 2011,15(3):578-594. ] doi: 10.11834/jrs.20110066
[8]	苗翠翠,江南,彭世揆,等.基于NDVI时序数据的水稻种植面积遥感监测分析——以江苏省为例[J].地球信息科学学报,2011,13(2):273-280. doi: 10.3724/SP.J.1047.2011.00273
	[ Miao C C, Jiang N, Peng S K, et al.Extraction of paddy land area based on NDVI time-series data: taking Jiangsu province as an example[J]. Journal of Geo-Information Science, 2011,13(2):273-280. ] doi: 10.3724/SP.J.1047.2011.00273
[9]	程乾,王人潮.数字高程模型和多时相MODIS数据复合的水稻种植面积遥感估算方法研究[J].农业工程学报,2005,21(5):89-92.
	[ Cheng Q, Wang R C.Estimation of the rice planting area using digital elevation model and multitemporal moderate resolution imaging spectroradiometer[J]. Transactions of the Chinese Society of Agricultural Engineering, 2005,21(5):89-92. ]
[10]	梁益同,万君.基于HJ-1A/B-CCD影像的湖北省冬小麦和油菜分布信息的提取方法[J].中国农业气象,2012,33(4):573-578. doi: 10.3969/j.issn.1000-6362.2012.04.015
	[ Liang Y T, Wan J.Application of HJ-1A/B-CCD images in extracting the distribution of winter wheat and rape in Hubei Province[J]. Chinese Journal of Agrometeorology, 2012,33(4):573-578. ] doi: 10.3969/j.issn.1000-6362.2012.04.015
[11]	Lennington R k, Sorensen C T, Heydorn R P. A mixture model approach for estimating crop areas from Landsat data[J]. Remote Sensing of Environment, 1984,14(1-3):197-206. doi: 10.1016/0034-4257(84)90015-4
[12]	李卫国,李花,王纪华,等.基于Landsat/TM遥感的冬小麦长势分级监测研究[J].麦类作物学报,2010,30(1):92-95. doi: 10.7606/j.issn.1009-1041.2010.01.019
	[ Li W G, Li H, Wang J H, et al.A study on classification and monitoring of winter wheat growth status by Landsat/TM image[J]. Journal of Triticeae Crops, 2010,30(1):92-95. ] doi: 10.7606/j.issn.1009-1041.2010.01.019
[13]	魏新彩,王新生,刘海,等.HJ图像水稻种植面积的识别分析[J].地球信息科学学报,2012,14(3):382-388. doi: 10.3724/SP.J.1047.2012.00382
	[ Wei X C, Wang X S, Liu H, et al.Extraction of paddy rice coverage based on the HJ satellite data[J]. Journal of Geo-Information Science, 2012,14(3):382-388. ] doi: 10.3724/SP.J.1047.2012.00382
[14]	耿月明,周枝乾,杨光军,等.江汉平原油菜生产状况及发展潜力分析[J].种子世界,2008,11(5):56-59. doi: 10.3969/j.issn.1000-8071.2008.11.003
	[ Geng Y M, Zhou Z Q, Yang G J, et al.Status of rapeseed production and growth potentials in Jianghan plains[J]. Seed World, 2008,11(5):56-59. ] doi: 10.3969/j.issn.1000-8071.2008.11.003
[15]	赵虎,裴志远,马尚杰,等.利用HJ-1-A/BCCD2数据反演冬小麦叶面积指数[J].农业工程学报,2012,28(10):172-176.
	[ Zhao H, Pei Z Y, Ma S J, et al.Retrieving leaf area index of winter wheat using HJ-1-A/B CCD2 data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012,28(10):172-176. ]
[16]	赵丽花,李卫国,杜培军.基于多时相HJ卫星的冬小麦面积提取[J].遥感信息,2011(2):41-45.
	[ Zhao L H, Li W G, Du P J.The area extraction of winter wheat based on multi-temporal HJ remote sensing satellite images[J]. Remote Sensing Information, 2011,2:41-45. ]
[17]	熊勤学,黄敬峰.利用NDVI指数时序特征监测秋收作物种植面积[J].农业工程学报,2009, 25(1):144-148.
	[ Xiong Q X, Huang J F.Estimation of autumn harvest crop planting area based on NDVI sequential characteristics[J]. Transactions of the Chinese Society of Agricultural Engineering, 2009,25(1):144-148. ]
[18]	郑盛,赵祥,张颢,等.HJ-1卫星CCD数据的大气校正及其效果分析[J].遥感学报,2011,15(4):709-721. doi: 10.11834/jrs.20110147
	[ Zheng S, Zhao X, Zhang H, et al.Atmospheric correction on CCD data of HJ-1 satellite and analysis of its effect[J]. Journal of Remote Sensing, 2011,15(4):709-721. ] doi: 10.11834/jrs.20110147

	4月			5月			6月			7月			8月			9月			10月			11月
	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬	下旬	上旬	中旬
双季早稻	播种	三叶		移栽	分蘖期			孕穗	抽穗开花	灌浆期		成熟收获
中稻			播种		三叶		移栽	分蘖期		拔节期		孕穗	抽穗开花		灌浆期		成熟收获
晚稻									播种	三叶		移栽	分蘖期		拔节孕穗		抽穗开花		灌浆期		成熟收获
棉花		播种	出苗		移栽	三叶	五叶		现蕾	开花期						裂铃吐絮期					停止生长	拔杆

影像时间	作物类型	平均值（M）	标准差 (D)	M+D （最大值）	M-D （最小值）
2012-05-12	中稻	0.304714	0.032654	0.337368	0.272060
	棉花	0.297284	0.038867	0.336151	0.258417
	晚稻	0.132896	0.052232	0.185128	0.080664
2012-07-02	中稻	0.814182	0.037437	0.851619	0.776745
	棉花	0.424456	0.062608	0.487064	0.361848
	晚稻	0.205091	0.073252	0.278343	0.131839
2012-07-30	中稻	0.790761	0.030821	0.821582	0.75994
	棉花	0.734289	0.045935	0.780224	0.688354
	晚稻	0.45214	0.046766	0.498906	0.405374
2012-08-08	中稻	0.657455	0.03216	0.689615	0.625295
	棉花	0.714886	0.029092	0.743978	0.685794
	晚稻	0.525734	0.018519	0.544253	0.507215
2012-09-28	中稻	0.348423	0.03674	0.385163	0.311683
	棉花	0.374898	0.067973	0.442871	0.306925
	晚稻	0.666095	0.025196	0.691291	0.640899

地物类型	地表真实像元分类数据
地物类型	水体	棉花	中稻	林地	总计	错分率/（%）	漏分率/（%）	整体精度/（%）	Kappa系数
水体	60	0	0	0	60	0	0	95.96	0.93
棉花	0	332	7	0	339	2.06	2.06
中稻	0	3	145	8	156	7.05	6.45
林地	0	4	3	57	64	10.94	12.31
总计	60	339	155	65	619	-	-

	2012年			2014年
	统计面积/(10³ha）	提取面积/(10³ha)	精度/（%）	统计面积/(10³ha）	提取面积/(10³ha)	精度/（%）	减少面积/(10³ha）
荆门市	10.19	9.23	90.58	9.41	8.51	90.44	0.72
钟祥市	13.93	13.56	97.34	11.94	9.77	81.83	3.79
京山县	11.30	9.58	84.78	3.08	2.28	74.03	7.30
荆州区	9.58	10.14	94.15	5.31	5.79	90.96	4.35
沙市区	3.12	3.85	76.60	2.94	2.73	92.86	1.12
公安县	28.27	32.70	84.33	25.66	23.42	91.27	9.28
监利县	22.44	17.15	76.43	14.37	12.80	89.07	4.35
江陵县	8.05	8.12	99.13	6.03	6.28	95.85	1.84
石首市	12.61	15.20	79.46	11.47	11.14	97.12	4.06
洪湖市	6.91	7.49	91.61	5.60	5.92	94.29	1.57
松滋市	13.83	11.98	86.62	11.78	10.4	88.29	1.58
仙桃市	23.16	19.26	83.16	15.17	13.33	87.87	5.93
潜江市	26.23	26.34	99.58	20.28	19.7	97.14	6.64
天门市	37.97	38.23	99.32	30.01	25.87	86.20	12.36
合计	227.59	222.83	97.91	173.05	157.94	91.27	64.89