冬小麦田空间变异性信息和采样策略研究.docxVIP

冬小麦田空间变异性信息和采样策略研究 采样点空间布局研究 随着国内外详细农业实践的发展，对空间信息的空间变换的研究越来越受到重视。人们利用仪器直接测试和化学分析的方法获取农田信息,利用地理信息系统绘制农田信息空间分布图,分析农田信息的空间分布情况,以便进行农田的精细灌溉、精细施肥和其他的精细农业生产管理。然而,要获得较准确的农田信息空间分布图,必须首先获取准确的采样信息。 精细农业描述田间信息空间变异性,一般采取田间规则网格定位采样,当网格较小时,采样数多,采样较为费钱费力。采样方式对田间信息的平均值和变异程度影响很大,因此,在满足一定精度的前提下,选择合理的采样点分布方式和确定合理的采样点个数非常重要。目前,国内外研究人员只是对采样数据进行直接的统计分析,求出采样个数,而对采样点空间布局的研究还不多。 为准确获取农田信息,对按规则栅格采集的数据点,按50m倍数的样点间隔抽取采样(不失一般性),并进行空间插值,提取对应原63个采样点位置的插值数据,然后进行统计分析,确定较准确又经济的采样方案。 1 试验条件 1.1 土壤理化性质的变化 试验田是一块冬小麦田,面积约13.3hm2,土壤为褐土,土质为沙壤土。土壤肥沃,通透性好,易耕作,抗涝。地势平坦,东部略高,西部略低。南部略高,北部略低。水资源丰富。生产过程已基本实现机械化。主要作物是小麦,近几年实行一年一茬小麦平作。 在经济许可的条件下,本着尽可能多采样的原则,获取土壤耕作层的主要土壤特性信息和产量的有关信息。试验采用50m×50m设置规则栅格,利用GPS接收机定位采样,试验于2002年3月18日进行。 麦田土壤的测试信息包括:土壤有机质、全氮、速效氮、速效磷、速效钾等土壤化学性质,以及土壤含水量、密度等土壤物理性质。所有土壤特性均在现场采样后,在实验室处理测定。 麦田产量数据的测试信息包括:籽粒产量、生物产量、千粒质量等。 1.2 采样点采样方案描述 按照50m×50m的规则网格进行田间信息采样,采样点的布局如图1a所示,共有63个采样点。然后按照50m的倍数抽取采样点,各种不同的采样方案见图1。 图1b为原数据按50m×100m的栅格采样,取原采样点中的35个。图1c为原数据按100m×100m的栅格采样,取原采样点中的20个。图1d、1e、1f分别为原数据按150m×150m、150m×150m左、150m×150m右的栅格采样,取原采样点中的9个。图1g为原数据按100m×200m的栅格采样,取原采样点中的12个。 2 学习方法 2.1 分布特征值选取 在一定的显著水平(α)和抽样允许误差Δ范围内,要求的合理采样数目n计算公式为 n=λα,f(Cv/Δ)2(1) 式中 λα,f——t分布特征值 Cv——变异系数 Δ可取5%、10%或15%。如果计算所得样本数n大于总体样本容量N的10%,则应采用不重复抽样公式 n′=n/(1+n/N) (2) 2.2 克立格插值法 空间插值方法有线性插值法、距离反比法、样条法和克立格插值法。本文采用地理信息系统中常用的克立格插值法进行地统计分析,绘制每一种田间信息不同采样率的空间分布图,提取有关的插值数据。克立格插值法是一种局部估计的加权平均,对各采集点的权重的确定是通过半方差图获得的,然后根据统计学上无偏和最优的要求,利用拉格朗日极小化原理,推导出权重与半方差之间的公式。空间插值采用点状克立格插值法。 2.3 采样方案a:lsd 当两样本均值差的绝对值超过最小显著差数,就可以认为两样本有差异,对应的新采样方案精度不高,不能通过最小显著差数检验(LSD),该采样方案不可取。计算公式为 Lα=tn?k,α2Mni???√(3)Lα=tn-k,α2Μni(3) 式中 ni——各样本的容量数 M——样本均方误差 本研究对不同的田间信息不同的采样方案,按照克立格法插值后,提取与原63个点位置对应的63个插值数据,组成LSD检验样本,而且以原63个采样点插值后的数据样本作为比较的参考样本。 3 计算与分析 3.1 小麦产量和平均ph的变异强度 按照经典统计方法,土壤有机质、全氮、速效氮、速效磷、速效钾、土壤密度的统计特征值列于表1。其中全氮和速效氮略偏低,土壤有机质、速效磷为中等肥力,速效钾比较丰富。这5种养分的变异系数都属于中等变异强度,其中速效磷和速效钾的变异强度比较大。土壤含水量的最大值与最小值相差较大,变异系数为0.151,属中等变异强度,造成这种情况的原因在于农田局部小区不平,灌溉时有径流。土壤密度的变异系数为0.066,变异强度较弱,这表明土壤的质地比较均匀,这与耕作制度有关。小麦产量的统计特征值也列于表1,从籽粒产量看,均值为6 426kg/hm2,属于高产量水平,该地块为高产田。籽粒产量的最大值为9 153.45