土壤特性空间变异性研究.pdfVIP

2005年1O月 农 业 机 械 学 报 第36卷第1O期 土壤特性空间变异性研究 陈翠英 江永真 袁朝春 【摘要】 在R环境中以半方差函数和克里格插值为基础，对镇江市瑞京园种植区土壤的pH值、有效钾、硝态 氮和铵态氮含量进行了空间变异性研究，并绘制了空间分布图。对于使用R语言进行土壤空间变异分析具有一定 的理论和实践意义。 关键词：土壤特性 克里格插值 空间变异 R语言 中图分类号：S159；S126 文献标识码：A StudyonSoilPropertySpatialVariabilityUsingR Language ChenCuiying JiangYongzhen YuanChaochun (JiangsuUniversity) Abstract Basedonthetheoryofsemivariogram andKriging．，inR environment，spatialvariabilityof pH，availablepotassium，nitrateandammonium inthesoilofZhenjiangRuijingParkwere analyzed，and spatialdistributionmapsweredrawn inR．Theresultswereapplicableto soil spatialanalysisusingR． Keywords Soilproperties，Kriging，Spatialvariability，R language 引言 1 R计算环境 田间实际情况表明，即使在土壤质地相同的区 R可 以被认为是 S语言的一个实现，和 S具有 域 内，同一时刻土壤特性在不同空间位置上也具有 很大的兼容性。在R中计算的一般步骤为：①准备 明显差异，这种属性称为土壤特性的空间变异性[1]。 原始数据，以向量的形式存储在变量中。②依据函数 土壤空间变异性对作物生长有直接的影响。在精确 进行计算，结果可存储为能被Excel识别的文本文 农业中，对农作物种植区土壤特性的空间变异性进 件，需要的话可以存储中间结果。③使用绘图函数绘 行描述，提供土壤不同特性的空间分布图是非常必 制出计算结果的图形表示。目前，R的一个很重要的 要的，能够为作物种植规划、施肥决策等提供科学依 发展方向是与GIS结合在一起，已经有R与GIS软 据 。 件GRASS之间的接 口程序[2]。与现有的一些GIS 本文以地统计学原理为基础，在R计算环境中 软件如ArcGIS相比较，使用R语言来进行地学统 对镇江市京口区高科技农业蔬菜种植园区(瑞京园) 计计算具有更大的灵活性，用户可以根据需要编写 中一个地块的土壤特性空间变异性进行了研究，绘 相应的计算程序。 制了不同土壤特性的空间分布图，对该地区制定合 本文在R中编写了相应的半方差计算和克里 理的施肥措施和进行合理的土壤养分管理具有直接 格计算函数。对于需要考虑不同方向半方差的情况， 的现实意义，也为精确农业的研究和实施提供一定 函数默认 4个方 向为 0。、45。、90。和 135。，对应东西、 的参考 。 南北和两个斜向，在半方差拟合时采用的是加权最 收稿 日期：2004—08—02 陈翠英 江苏大学机械工程学院 教授 博士生导师，212013 镇江市 江永真 江苏大学机械工程学院 硕士生 袁朝春 江苏大学机械工程学院 博士生