一种基于定性和定量辅助变量的土壤有机质空间分布预测方法
【专利摘要】本发明公开了一种基于定性和定量辅助变量的土壤有机质空间分布预测方法，其特征在于：步骤如下：1）数据来源；2）辅助数据来源与处理；3）研究方法；4）对照方法。由于采用了上述技术方案，与现有技术相比，本发明运用人工神经网络模型，在融合土壤类型、地形因子以及植被指数等定性和定量辅助环境变量的基础上，开展土壤有机质含量的空间分布预测，以期为区域高精度土壤性质空间分布预测提供方法参考。
【专利说明】一种基于定性和定量辅助变量的土壤有机质空间分布预测方法
【技术领域】
[0001]本发明涉及一种基于定性和定量辅助变量的土壤有机质空间分布预测方法。
【背景技术】
[0002]土壤有机质是土壤特性的重要指标，也是植物矿质营养和有机营养的源泉，在区域土壤肥力、农业生产和环境保护等方面均具有重要的作用和意义。由于多种成土因素和生态过程的影响，土壤有机质含量在空间上呈现非均匀分布。通常基于采样获得的样点土壤有机质含量信息，很难满足区域农业生产和环境保护中对面域土壤有机质含量分布特征的实际需求。因此，精确获取土壤有机质含量的空间分布信息，是准确掌握区域土壤肥力状况、科学的土壤养分管理和区域环境保护的实际需要。
[0003]土壤性质的空间分布格局是在各种地表环境因素共同作用下形成的。引入相关地表环境因子作为辅助变量的土壤性质空间分布预测方法，在不同程度上考虑了相关环境因子对土壤性质空间分布的影响，其预测精度较仅基于样点数据进行空间内插的方法有明显的提高。在各类辅助变量中，定量因子如气象、地形特征参数等因子相对易于操作，被广泛用于土壤性质空间分布预测。定性因素如土地利用方式、土壤类型和母质类型等对土壤性质空间变异也具有显著的影响，该类因素甚至会弱化定量因子对土壤性质空间变异的作用。因此，在预测方法中同时引入定性和定量辅助变量是进一步提高土壤性质空间分布预测精度的必然要求和最有效的途径之一。然而，在土壤性质空间分布预测中，定性因素的使用明显偏少。
[0004]在引入环境变量的预测方法中，多元回归模型、协同克里格法以及回归克里格法等方法被广泛用于刻画土壤性质与环境因子间的关系，进而预测土壤性质的空间分布。然而，土壤性质与环境因素间是一种极其复杂的非线性关系，准确表达土壤性质与环境因素间的复杂关系是提高预测精度的关键。人工神经网络模型能较准确地揭示复杂系统中多元环境变量与目标变量间的非线性映射关系，被成功应用于描述土壤与环境因子间的非线性关系，在土壤性质的空间分布预测研究中受到越来越多的关注。
【发明内容】
[0005]本发明要解决的技术问题是克服上述缺陷，提供一种运用人工神经网络模型，在融合土壤类型、地形因子以及植被指数等定性和定量辅助环境变量的基础上，开展土壤有机质含量的空间分布预测的基于定性和定量辅助变量的土壤有机质空间分布预测方法。
[0006]为解决上述问题，本发明所采用的技术方案是:
[0007]—种基于定性和定量辅助变量的土壤有机质空间分布预测方法，其特征在于:步骤如下:
[0008]I)、数据来源:以选取的研究区地形图、土地利用现状图和土壤类型分布图为辅助信息，在兼顾代表性与均匀性原则的基础上综合考虑地形、土壤类型等信息进行土壤采样点的布设，在每个土壤采样点处采用多点混合的方法采集表层土壤样品；每个采样点均以GPS记录其地理坐标和海拔高度，同时详细记录样点所在的环境信息；采集到的样品带回实验室经自然风干后磨碎过Imm筛，采用重铬酸钾容量法测定每个土样有机质含量；
[0009]2)、辅助数据来源与处理:
[0010]研究中定量因子主要选择了地形因子和植被指数；地形因子是土壤性质空间分布预测中最常用的辅助环境变量，包括高程(H)、坡度(S)、平面曲率(Ct)、剖面曲率(Cp)及地形湿度指数(TI);基于研究区1:5万地形图在ArcGIS9.3中生成30m分辨率的数据高程模型，进一步利用ArcGIS9.3的空间分析和栅格计算功能获取研究区30m分辨率的坡度、平面曲率、剖面曲率及地形湿度指数分布图；
[0011]基于遥感影像获取的植被指数能较好地反映区域地表植被的生长状态和植被覆盖度信息，是常用于土壤性质空间分布预测的又一辅助变量；在研究区采用MODIS遥感数据产品中的增强型植被指数作为预测研究区土壤有机质含量空间分布的定量植被因子；
[0012]3)、研究方法:
[0013]研究方法采用的是融合定性和定量辅助变量的径向基函数神经网络方法；
[0014]将各采样点土壤有机质的实测值和环境因素之间的关系表示为:
[0015]Z(xi；k, yJ；k)=m(k) +r (xi；k, yJ；k) (I)
[0016]r (xi；k, yJ； k) =f (t (xi； Yj), v (xi； Yj),...) (2)
[0017]式中:Z(Xi，k，Yj.k)表示属于第k种土壤类型的土样在(Xi，Yj)处的有机质含量值，(Xi, Yj)为采样点坐标，其行列号分别为i和j ；m(k)为第k种土壤类型有机质含量平均值；r(Xi,k) Yj,k)为采样点(Xi，yj)减去该点所属土壤类型有机质含量平均值后的残差。t(xi)yj)为点(Xi，yj)的地形参数值，V(Xil7j)为点(Xi，yj)的植被指数值。即假定特定点位上土壤有机质的变异首先由该点所处的土壤类型决定，而其残差值则由该点所处的地形植被条件等局地环境因素确定；平均值m(k)和残差r(Xi，k，yik)相互独立。
[0018]各点位残差值与局地环境因素间的关系采用径向基函数人工神经网络模型表达；该模型是一种具有单隐层的三层前馈型神经网络模型，预测方法过程如下:
[0019](I)、依据采样点土壤有机质含量值，分别统计研究区中各个土壤类型土壤有机质含量的平均值及每个土壤采样点对应的残差值；
[0020]( 2 )、依据计算出的各个土壤类型有机质含量平均值和研究区土壤类型分布图，制作研究区各土壤类型有机质值含量平均值分布图；
[0021](3)、以各土壤采样点的地形及植被特征参数值作为网络输入，以对应点位归一化后的土壤有机质残差值作为网络输出，构建神经网络模型，在MATLAB中模拟得到研究区土壤有机质的残差分布图；
[0022](4)、最后将研究区各土壤类型有机质值平均值分布图加上神经网络模型预测得到的残差分布图，即得到研究区土壤有机质含量空间分布图；
[0023]4)、对照方法:
[0024]选用普通克里格法、回归克里格法及神经网络模型与普通克里格结合的方法作为对照方法；回归克里格法首先采用多元逐步回归进行土壤有机质与环境因子间的回归预测，再以普通克里法对回归预测结果残差进行插值，最后将回归预测结果和普通克里格法对残差的估计值相加，即得到研究区土壤有机质的空间分布图；神经网络模型与普通克里格结合的方法则是以神经网络模型替代回归克里格法中的回归模型完成对研究区土壤有机质的空间分布预测。
[0025]由于采用了上述技术方案，与现有技术相比，本发明以四川盆地丘陵区的三台县为研究区，运用人工神经网络模型，在融合土壤类型、地形因子以及植被指数等定性和定量辅助环境变量的基础上，开展土壤有机质含量的空间分布预测，以期为区域高精度土壤性质空间分布预测提供方法参考。
【专利附图】
【附图说明】
[0026]图1是研究区土壤有机质含量及其对数转换值的频率分布图；
[0027]图2是研究区土壤有机质各数据项半方差函数图；
[0028]图3是不同预测方法对土壤有机质含量空间分布的预测结果图。
【具体实施方式】
[0029]实施例:
[0030]一种基于定性和定量辅助变量的土壤有机质空间分布预测方法，步骤如下:
[0031]I)、数据来源:以选取的研究区地形图、土地利用现状图和土壤类型分布图为辅助信息，在兼顾代表性与均匀性原则的基础上综合考虑地形、土壤类型等信息进行土壤采样点的布设，在每个土壤采样点处采用多点混合的方法采集表层土壤样品。全区共采集土样2346个采样点，每个采样点均以GPS记录其地理坐标和海拔高度，同时详细记录样点所在的环境信息如土壤类型、耕种作物和农地轮作方式等信息。采集到的样品带回实验室经自然风干后磨碎过1_筛，采用重铬酸钾容量法测定每个土样有机质含量。
[0032]2)、辅助数据来源与处理:
[0033]研究中定量因子主要选择了地形因子和植被指数；地形因子是土壤性质空间分布预测中最常用的辅助环境变量，包括高程(H)、坡度(S)、平面曲率(Ct)、剖面曲率(Cp)及地形湿度指数(TI);基于研究区1:5万地形图在ArcGIS9.3中生成30m分辨率的数据高程模型，进一步利用ArcGIS9.3的空间分析和栅格计算功能获取研究区30m分辨率的坡度、平面曲率、剖面曲率及地形湿度指数分布图。
[0034]基于遥感影像获取的植被指数能较好地反映区域地表植被的生长状态和植被覆盖度信息，是常用于土壤性质空间分布预测的又一辅助变量；在研究区采用MODIS遥感数据产品中的增强型植被指数作为预测研究区土壤有机质含量空间分布的定量植被因子。考虑到土壤有机质的积累特点，本发明中选用了 2001至2006年每年7月和8月16d合成的250m分辨率的EVI数据，在对所下载的EVI数据进行质量检查和处理后取平均值并重采样至30m分辨率，即用2001至2006年生物量最大时的平均植被指数作为定量辅助变量之一。为方便计算，所有定量因子均进行归一化处理。
[0035]研究表明，丘陵区土壤类型较土地利用类型对土壤有机质含量空间分布的影响大，因此研究中选择研究区土壤类型作为本研究的定性辅助变量。研究区土壤类型分布图来源于第二次土壤普查。在对研究区土壤类型图扫描后利用ArcGIS9.3进行矢量化，得到研究区数字土壤类型分布图。
[0036]3)、研究方法:[0037]研究方法采用的是融合定性和定量辅助变量的径向基函数神经网络方法；
[0038]将各采样点土壤有机质的实测值和环境因素之间的关系表示为:
[0039]Z(xuk, yJ；k)=m(k) +r(xi；k, yJ；k) (I)
[0040]r (xi；k, yJ；k) =f (t (xi； Yj), v (xi； Yj),...) (2)
[0041]式中:Z(Xi，k，Yj.k)表示属于第k种土壤类型的土样在(Xi，Yj)处的有机质含量值，(Xi, Yj)为采样点坐标，其行列号分别为i和j ；m(k)为第k种土壤类型有机质含量平均值；r(Xi,k) Yj,k)为采样点(Xi，yj)减去该点所属土壤类型有机质含量平均值后的残差。t(xi)yj)为点(Xi，yj)的地形参数值，V(Xil7j)为点(Xi，yj)的植被指数值。即假定特定点位上土壤有机质的变异首先由该点所处的土壤类型决定，而其残差值则由该点所处的地形植被条件等局地环境因素确定；平均值m(k)和残差r(Xi，k，yik)相互独立。
[0042]各点位残差值与局地环境因素间的关系采用径向基函数人工神经网络模型表达；该模型是一种具有单隐层的三层前馈型神经网络模型，预测方法过程如下:
[0043](I)、依据采样点土壤有机质含量值，分别统计研究区中各个土壤类型土壤有机质含量的平均值及每个土壤采样点对应的残差值；
[0044]( 2 )、依据计算出的各个土壤类型有机质含量平均值和研究区土壤类型分布图，制作研究区各土壤类型有机质值含量平均值分布图；
[0045](3)、以各土壤采样点的地形及植被特征参数值作为网络输入，以对应点位归一化后的土壤有机质残差值作为网络输出，构建神经网络模型，在MATLAB中模拟得到研究区土壤有机质的残差分布图；
[0046](4)、最后将研究区各土壤类型有机质值平均值分布图加上神经网络模型预测得到的残差分布图，即得到研究区土壤有机质含量空间分布图；
[0047]4)、对照方法:
[0048]选用普通克里格法、回归克里格法及神经网络模型与普通克里格结合的方法作为对照方法；回归克里格法首先采用多元逐步回归进行土壤有机质与环境因子间的回归预测，再以普通克里法对回归预测结果残差进行插值，最后将回归预测结果和普通克里格法对残差的估计值相加，即得到研究区土壤有机质的空间分布图；神经网络模型与普通克里格结合的方法则是以神经网络模型替代回归克里格法中的回归模型完成对研究区土壤有机质的空间分布预测。
[0049]针对上述基于定性和定量辅助变量的土壤有机质空间分布预测方法预测精度评价方法采用独立验证样点集对各方法预测结果进行评价，即随机从2346个土壤样点提取20%(469个)作为验证点，其余80%(1877个)作为建模点，以平均绝对误差(MAE)、均方根误差(RMSE)和平均相对误差(MRE)对验证点的预测值和实际观测值进行对比分析，得出精度评价结果。
[0050]上述基于定性和定量辅助变量的土壤有机质空间分布预测方法结果与分析:
[0051]1、描述性统计分析。
[0052]根据建模点统计分析表明，如表I所示，研究区土壤有机质质量分数的平均值为17.97g/kg ;变异系数处于10?100%之间，为中等程度的变异性。不同土壤类型土壤有机质含量差异较大，水稻土有机质质量分数平均值达到22.49g/kg，远高于其它3个土壤类型；其次是潮土和紫色土，平均值分别为15.87和14.09g.kg-1 ;黄壤土最小，平均值为13.66g.kg-1。从变异程度来看，水稻土和紫色土有机质含量变异程度略高于另2个土类，但均属于中等程度变异。
[0053]从数据的分布形态来看，原始数据呈明显的偏态分布，通过对数转换后其偏度值和峰度值以及分布频率图明显更接近正态分布，如表1、图1所示。因此，对数转换后的数值用于研究中各种预测方法的计算。
[0054]表1研究区不同土壤类型有机质含量描述性统计特征
【权利要求】
1.一种基于定性和定量辅助变量的土壤有机质空间分布预测方法，其特征在于:步骤如下: 1)、数据来源:以选取的研究区地形图、土地利用现状图和土壤类型分布图为辅助信息，在兼顾代表性与均匀性原则的基础上综合考虑地形、土壤类型等信息进行土壤采样点的布设，在每个土壤采样点处采用多点混合的方法采集表层土壤样品；每个采样点均以GPS记录其地理坐标和海拔高度，同时详细记录样点所在的环境信息；采集到的样品带回实验室经自然风干后磨碎过Imm筛，采用重铬酸钾容量法测定每个土样有机质含量； 2)、辅助数据来源与处理: 研究中定量因子主要选择了地形因子和植被指数；地形因子是土壤性质空间分布预测中最常用的辅助环境变量，包括高程(H)、坡度(S)、平面曲率(Ct)、剖面曲率(Cp)及地形湿度指数(TI);基于研究区1:5万地形图在ArcGIS9.3中生成30m分辨率的数据高程模型，进一步利用ArcGIS9.3的空间分析和栅格计算功能获取研究区30m分辨率的坡度、平面曲率、剖面曲率及地形湿度指数分布图； 基于遥感影像获取的植被指数能较好地反映区域地表植被的生长状态和植被覆盖度信息，是常用于土壤性质空间分布预测的又一辅助变量；在研究区采用MODIS遥感数据产品中的增强型植被指数作为预测研究区土壤有机质含量空间分布的定量植被因子； 3)、研究方法: 研究方法采用的是融合定性和定量辅助变量的径向基函数神经网络方法； 将各采样点土壤有机质的实测值和环境因素之间的关系表示为:
Z(xi；k, yJ;k`)=m(k) +r (xi；k, yJ；k) (I)
r (xi；k, yJ； k) =f (t(xi, Yj), v (Xi, Yj),...) (2) 式中:Z(Xi，k，yj，k)表示属于第k种土壤类型的土样在(Xi，yj)处的有机质含量值，(Xi, Yj)为采样点坐标，其行列号分别为i和j ；m(k)为第k种土壤类型有机质含量平均值；r(Xi,k) Yj,k)为采样点(Xi，yj)减去该点所属土壤类型有机质含量平均值后的残差。t(xi)yj)为点(Xi，yj)的地形参数值，V(Xil7j)为点(Xi，yj)的植被指数值。即假定特定点位上土壤有机质的变异首先由该点所处的土壤类型决定，而其残差值则由该点所处的地形植被条件等局地环境因素确定；平均值m(k)和残差r(Xi，k，yik)相互独立； 各点位残差值与局地环境因素间的关系采用径向基函数人工神经网络模型表达；该模型是一种具有单隐层的三层前馈型神经网络模型，预测方法过程如下: (1)、依据采样点土壤有机质含量值，分别统计研究区中各个土壤类型土壤有机质含量的平均值及每个土壤采样点对应的残差值； (2)、依据计算出的各个土壤类型有机质含量平均值和研究区土壤类型分布图，制作研究区各土壤类型有机质值含量平均值分布图； (3)、以各土壤采样点的地形及植被特征参数值作为网络输入，以对应点位归一化后的土壤有机质残差值作为网络输出，构建神经网络模型，在MATLAB中模拟得到研究区土壤有机质的残差分布图； (4)、最后将研究区各土壤类型有机质值平均值分布图加上神经网络模型预测得到的残差分布图，即得到研究区土壤有机质含量空间分布图； 4)、对照方法:选用普通克里格法、回归克里格法及神经网络模型与普通克里格结合的方法作为对照方法；回归克里格法首先采用多元逐步回归进行土壤有机质与环境因子间的回归预测，再以普通克里法对回归预测结果残差进行插值，最后将回归预测结果和普通克里格法对残差的估计值相加，即得到研究区土壤有机质的空间分布图；神经网络模型与普通克里格结合的方法则是以神经网络模型替代回归克里格法中的回归模型完成对研究区土壤有机质的空间分布预测； 5)、预测精度评价方法 采用独立验证样点集对各方法预测结果进行评价，即随机从2346个土壤样点提取20% (469个)作为验证点，其余80% (1877个)作为建模点(图1b)，以平均绝对误差(MAE)、均方根误差(RMSE)和平均相对误差(MRE)对验证点的预测值和实际观测值进行对比分析，得出精度评价结果。`
【文档编号】G06F19/00GK103529189SQ201310268497
【公开日】2014年1月22日 申请日期:2013年6月28日 优先权日:2013年6月28日
【发明者】李启权, 王昌全, 李冰, 张新, 高雪松, 张毅, 袁大刚 申请人:四川农业大学