基于无人机多光谱影像的春玉米田管理分区研究

doi:10.13304/j.nykjdb.2021.0618

摘要/Abstract

摘要：

为了提高大尺度农田管理的针对性，探寻低成本高效的分区方法，以黑龙江省典型黑土区30.8 hm²春玉米田为对象进行分区研究。基于吐丝期无人机多光谱影像，使用多尺度分割与模糊聚类相结合的方法进行分区，同时基于播种前的土壤养分（土壤有机质、速效氮磷钾）、土壤体积含水率、电导率、pH进行模糊聚类分区并作为对照，对分区间春玉米产量和土壤养分进行方差分析，并对分区内变异系数进行比较以评价分区效果。结果表明，基于无人机影像得到4个较优管理分区为M1、M2、M3、M4，各分区产量分别为7 597.53、8 236.35、8 686.98、9 119.93 kg·hm^-2，各分区间产量差异显著，其中M1、M2、M3间土壤养分差异显著（P<0.05），分区内作物产量和土壤养分的变异系数降低；基于土壤数据确定4个分区，即S1、S2、S3、S4，春玉米产量分别为7 754.81、8 173.44、8 860.05、9 153.23 kg·hm^-2，分区间土壤养分、土壤水分差异显著（P<0.05），分区内部土壤均一性提高。综合来看，2种方法的划分结果在空间分布上具有一定的相似性，同级分区的空间重合度分别为40.00%、46.51%、57.45%、59.38%，整体重合度为51.32%。当缺乏土壤数据时，无人机多光谱影像可为农田管理分区提供参考依据。

关键词: 管理分区, 多光谱影像, 无人机, 作物长势, 土壤养分

Abstract:

In order to improve the pertinence of practice in a large scale farm， and explore an efficient method for delineating management zones， a 30.8 hm² typical black farmland of spring maize in Heilongjiang province was taking as object to study management zones. Based on the multispectral images at silking stage of maize obtained by unmanned aerial vehicles（UVA） equipped with Sequoia sensor， a method of multiresolution segmentation combined with fuzzy-C means clustering algorithm was used to delineate the farmland. At the same time， segmentation by fuzzy-C means clustering algorithm based on grid sampling and soil properties， including soil organic matter， available NPK， soil volumetric moisture content， soil conductivity， and pH， was conducted as a contrast. The coefficient of variable （CV） of soil properties and yield were chosen to evaluate effect of the 2 methods. The results showed that the partition based on UAV images was delineated to 4 optimal management zones named as M1， M2， M3， M4， and their yields were 7 597.53，8 236.35，8 686.98，9 119.93 kg·hm^-2， respectively. There were significant differences in maize yield between each zones， and in soil nutrient between M1， M2， M3 zones， and the CV of soil nutrient， crop yield among each zones were much lower than that before the management zones. It was also delineated to 4 optimal zones determined based on soil properties named as S1， S2， S3 and S4， and their yields were 7 754.81，8 173.44，8 860.05，9 153.23 kg·hm^-2. There were significant differences in soil nutrients and soil water content between each zones （P<0.05）， and the homogeneity of soil properties within the zones was improved. The 2 methods had certain similarity in spatial distribution. The spatial coincidence degree under the same level zones were 40.00%， 46.51%， 57.45%， 59.38%， respectively， and the overall coincidence degree was 51.32%. Further，the results based on UAV multispectral images provided basis for farmland management zoning in the absence of soil data.

Key words: management zones, multispectral, UAV, crop growth, soil nutrients

中图分类号:

S126

卢闯, 胡海棠, 覃苑, 淮贺举, 李存军. 基于无人机多光谱影像的春玉米田管理分区研究[J]. 中国农业科技导报, 2022, 24(9): 106-115.

Chuang LU, Haitang HU, Yuan QIN, Heju HUAI, Cunjun LI. Delineating Management Zones in Spring Maize Field Based on UAV Multispectral Image[J]. Journal of Agricultural Science and Technology, 2022, 24(9): 106-115.

图/表 10

参考文献 23

1	CRUZ C, SILVA J T, OCHO A, et al.. Delineation of environmental units by multivariate techniques in the duero river watershed, Michoacan, Mexico [J]. Environ. Model. Assess., 2016, 22(3):1-10.
2	刘刚,闫静雯,谢云,等.黑土坡耕地土壤有机质空间变异及其与土壤侵蚀的关系——以黑龙江省鹤山农场为例[J].地理科学, 2016, 36(11): 1751-1758.
	LIU G, YAN J W, XIE Y, et al.. Spatial variation of soil organic matter on black soil sloping cropland and its relationship with soil erosion: a case study of Heshan farm in Heilongjiang province [J]. Sci. Geographica Sin., 2016, 36(11): 1751-1758.
3	OLDONI H, TERRA V, TIMM L C, et al.. Delineation of management zones in a peach orchard using multivariate and geostatistical analyses [J]. Soil Tillage Res., 2019, 191: 1-10.
4	GILI A, ALVAREZ C, BAGNATO R, et al.. Comparison of three methods for delineating management zones for site-specific crop management [J]. Comput. Electronics Agric., 2017, 139: 213-223.
5	张泽,吕新,吕宁,等.基于GIS、RS的滴灌棉田土壤养分精确管理分区研究[J].农业机械学报,2014,45(7):125-132.
	ZHANG Z, LV X, LV N, et al.. Defining agricultural management zones using remote sensing and GIS techniques for drip-irrigated cotton fields [J]. Trans. Chin. Soc. Agric. Mach.., 2014, 45(7): 125-132.
6	BREUNIG F M, GALVAO L S, DALAGNOL R, et al.. Delineation of management zones in agricultural fields using cover-crop biomass estimates from PlanetScope data [J/OL]. Int. J. Appl. Earth OBS, 2020, 85: 102004 [2022-05-27]. .
7	KUIAWSKI A C, SAFANELLI J L, BOTTEGA E, et al.. Vegetation indexes and delineation of management zones for soybean [J]. Pesqui. Agrop. Trop., 2017, 47(2): 168-177.
8	乔浪,张智勇,陈龙胜,等.基于无人机图像的玉米冠层叶绿素含量检测与分布研究[J].农业机械学报,2019,50:182-186.
	QIAO L, ZHANG Z Y, CHEN L S, et al.. Chlorophyll content detection and distribution research of maize canopy based on UAV image [J]. Trans. Chin. Soc. Agric. Mach., 2019, 50: 182-186.
9	韩文霆,彭星硕,张立元,等.基于多时相无人机遥感植被指数的夏玉米产量估算[J].农业机械学报,2020,51(1):149-155.
	HAN W T, PENG X S, ZHANG L Y, et al.. Summer maize yield estimation based on vegetation ondex derived from multitemporal UAV remote sensing [J]. Trans. Chin. Soc. Agric. Mach., 2020,51(1):149-155.
10	LEROUX C, JONES H, CLENET A, et al.. A new approach for zoning irregularly-spaced, within-field data [J]. Comput. Electron. Agric., 2017,141:196-206.
11	BENZ U C, HOFMANN P, WILLHAUCK, et al.. Multi-resolution, object-oriented fuzzy analysis of remote sensing data for GIS-ready information [J]. ISPRS-J. Photogramm. Remote Sens., 2004,58(3-4): 239-258.
12	马龙飞,胡乃月,李伟,等.利用无人机多光谱数据监测玉米对不同灌溉模式的响应差异[J]. 植物营养与肥料学报,2022,28(4):743-753.
	MA L F, HU N Y, LI W, et al.. Using multispectral drone data to monitor maize's response to various irrigation modes [J]. J. Plant Nutr. Fert., 2022,28(4):743-753.
13	田庆久.植被指数研究进展[J].地球科学进展,1998,13(4):328-333.
	TIAN Q J. Advances in study on vegetation indices [J]. Adv. Earth Sci., 1998, 13(4):328-333.
14	陈春雷,武刚.面向对象的遥感影像最优分割尺度评价[J].遥感技术与应用,2011,26(1):96-102.
	CHEN C L, WU G. Evaluation of optimal segmentation scale with object-oriented method in remote sensing [J]. Remote Sen. Technol. Appl., 2011, 26(1):96-102.
15	王露,刘庆元.面向对象的高分辨率遥感影像多尺度分割参数及分类研究[J].测绘与空间地理信息,2015,38(3):167-169.
	WANG L, LIU Q Y. The methods summary of optimal segmentation scale selection in high resolution remote sensing images multiscale segmentation [J]. Geomatics Spatial Inf. Technol., 2015, 38(3): 167-169.
16	GAVIOLI A, SOUZA E G, BAZZI C L, et al.. Optimization of management zone delineation by using spatial principal components [J]. Comput. Electron. Agric., 2016. 127: 302-310.
17	FRIDGEN J J, NEWELL R, KITCHEN, et al.. Management zone analyst (MZA): software for subfield management zone delineation [J]. Agron. J., 2004. 96(1): 100-108.
18	刘建刚,赵春江,杨贵军,等.无人机遥感解析田间作物表型信息研究进展[J].农业工程学报,2016,32(24):98-106.
	LIU J G, ZHAO C J, YANG G J, et al.. Review of field based phenotyping by unmanned aerial vehicle remote sensing platforrm [J]. Trans. Chin. Soc. Agric. Eng., 2016, 32(24): 98-106.
19	CHANG D, ZHANG J, ZHU L, et al.. Delineation of management zones using an active canopy sensor for a tobacco field [J]. Comput. Electron. Agric., 2014, 109: 172-178.
20	LI F, MIAO Y X, FENG G H, et al.. Improving estimation of summer maize nitrogen status with red edge-based spectral vegetation indices [J]. Field Crops Res., 2014, 157: 111-123.
21	赵明明.基于无人机与卫星影像的农田精准管理分区研究[D].哈尔滨:东北农业大学,2019.
	ZHAO M M. Study on site-specific management zones in field using unmanned aerial vehicle and satellite images [D]. Harbin: Northeast Agricultural University, 2019.
22	任佰朝,张吉旺,董树亭,等.生育前期淹水对夏玉米冠层结构和光合特性的影响[J].中国农业科学,2017,50(11):2093-2103.
	REN B C, ZHANG J W, DONG S T, et al.. Effect of waterlogging at early period on canopy structure and photosynthetic characteristics of summer maize [J]. Sci. Agric. Sin., 2017, 50(11): 2093-2103.
23	葛翠萍,赵军,王秀峰,等.东北黑土区坡耕地地形因子对土壤水分和容重的影响[J].水土保持通报,2008,28(6):16-21.
	GE C P, ZHAO J, WANG X F, et al.. Influence of topographic factors on soil water and bulk density in a typical slope land in the black soil area of Northeast China [J]. Bull. Soil Water Conser., 2008, 28(6): 16-21.

土壤指标 Soil properties	样点数目 Sample points No.	平均值Average	最小值Min.	最大值Max.	变异系数CV	偏度Skewness	峰度Kurtosis
有机质Organic matter /%	115	6.52	4.25	9.09	17.41	0.10	0.84
碱解氮Available nitrogen / （mg·kg^-1）	115	238.80	158.11	320.05	19.57	-0.20	-0.04
有效磷/Available phosphorus/ （mg·kg^-1）	115	85.60	39.75	145.90	32.07	0.22	0.38
速效钾/ Available potassium/ （mg·kg^-1）	115	211.00	149.17	298.44	19.53	0.83	0.49
土壤含水率Soil water content/%	115	29.81	19.76	38.47	16.37	0.18	-0.13
土壤电导率 Soil EC/（μs·cm^-1）	115	325.00	193.00	446.67	15.37	-0.09	0.33
pH	115	5.41	5.10	5.60	4.38	-0.56	0.77

土壤指标 Soil properties	样点数目 Sample points No.	平均值Average	最小值Min.	最大值Max.	变异系数CV	偏度Skewness	峰度Kurtosis
有机质Organic matter /%	115	6.52	4.25	9.09	17.41	0.10	0.84
碱解氮Available nitrogen / （mg·kg^-1）	115	238.80	158.11	320.05	19.57	-0.20	-0.04
有效磷/Available phosphorus/ （mg·kg^-1）	115	85.60	39.75	145.90	32.07	0.22	0.38
速效钾/ Available potassium/ （mg·kg^-1）	115	211.00	149.17	298.44	19.53	0.83	0.49
土壤含水率Soil water content/%	115	29.81	19.76	38.47	16.37	0.18	-0.13
土壤电导率 Soil EC/（μs·cm^-1）	115	325.00	193.00	446.67	15.37	-0.09	0.33
pH	115	5.41	5.10	5.60	4.38	-0.56	0.77

植被指数 Vegetation index	叶面积指数 LAI	地上生物量 Above-ground biomass	株高 Plant height	SPAD
NDVI	0.202^*	0.187^*	0.150	0.063 4
OSAVI	0.261^**	0.213^*	0.167	0.115 0
NDREI	0.313^**	0.297^**	0.259^**	0.143 0
RVI	0.151	0.129	0.047	0.051 4
DVI	0.348^**	0.240^*	0.257^**	0.146 0
SAVI	0.295^**	0.231^*	0.218^*	0.135 0
RDVI	0.298^**	0.226^*	0.204^*	0.128 0

植被指数 Vegetation index	叶面积指数 LAI	地上生物量 Above-ground biomass	株高 Plant height	SPAD
NDVI	0.202^*	0.187^*	0.150	0.063 4
OSAVI	0.261^**	0.213^*	0.167	0.115 0
NDREI	0.313^**	0.297^**	0.259^**	0.143 0
RVI	0.151	0.129	0.047	0.051 4
DVI	0.348^**	0.240^*	0.257^**	0.146 0
SAVI	0.295^**	0.231^*	0.218^*	0.135 0
RDVI	0.298^**	0.226^*	0.204^*	0.128 0

主成分Principal component	特征值Eigenvalue	贡献率Contribution rate/%	累积贡献率Accumulative contribution/%
PC1	2.45	34.95	34.95
PC2	1.97	28.12	63.07
PC3	1.77	25.22	88.29
PC4	0.45	6.49	94.78
PC5	0.20	2.91	97.69
PC6	0.11	1.56	99.25
PC7	0.05	0.75	100.00