于 洋，劉吉平

吉林省農(nóng)田土壤速效鉀插值精度研究

于 洋，劉吉平

本文以東北典型黑土區(qū)——吉林省榆樹市弓棚鎮(zhèn)13號(hào)村為研究區(qū)域，研究插值方法對(duì)農(nóng)田土壤速效鉀插值精度的影響。研究結(jié)果表明：在采樣點(diǎn)數(shù)較多的情況下，速效鉀的Kriging插值精度明顯高于BP神經(jīng)網(wǎng)絡(luò)的插值精度。隨著采樣點(diǎn)數(shù)的減少，Kriging插值精度顯著下降，而BP神經(jīng)網(wǎng)絡(luò)的插值精度相對(duì)穩(wěn)定。在采樣點(diǎn)數(shù)較少的情況下，BP神經(jīng)網(wǎng)絡(luò)的插值精度接近Kriging插值精度。

Kriging插值；BP神經(jīng)網(wǎng)絡(luò)；土壤速效鉀；插值精度

一、引言

精準(zhǔn)農(nóng)業(yè)是近年來農(nóng)業(yè)科學(xué)研究的熱點(diǎn)。在精準(zhǔn)農(nóng)業(yè)的實(shí)施過程中，需要采集大量的空間信息樣本，花費(fèi)大量的人力、物力和財(cái)力。因此，研究如何用較少的樣本來反映田間信息的空間變異規(guī)律，再用科學(xué)的插值方法進(jìn)行插值和預(yù)估是精準(zhǔn)農(nóng)業(yè)研究中的一個(gè)關(guān)鍵問題。本文以東北典型黑土區(qū)——吉林省榆樹市弓棚鎮(zhèn)13號(hào)村為研究區(qū)域，分析插值方法對(duì)插值精度的影響，為精準(zhǔn)農(nóng)業(yè)土壤養(yǎng)分插值方法的選取、制定優(yōu)化采樣策略提供科學(xué)依據(jù)。

二、研究方法與數(shù)據(jù)處理

（一）研究方法

1.Kriging插值。首先以SPSS 15.0為平臺(tái)分析土壤速效鉀的基本統(tǒng)計(jì)特征，檢驗(yàn)土壤樣本是否符合正態(tài)分布；然后在GS+3.0下，分析土壤堿解氮的空間變異特性，確定Kriging插值所采用的最優(yōu)模型以及空間變異參數(shù)；最后以ArcGIS9.3為平臺(tái)，完成Kriging插值。

2.BP神經(jīng)網(wǎng)絡(luò)模型。本文神經(jīng)網(wǎng)絡(luò)插值以MATLAB為平臺(tái)，利用newff函數(shù)建立三層網(wǎng)絡(luò)，第1、2層傳遞函數(shù)為tansig，第3層傳遞函數(shù)為purelin；網(wǎng)絡(luò)訓(xùn)練函數(shù)采用動(dòng)量及自適應(yīng)lrBP的梯度遞減訓(xùn)練函數(shù)traincgf；網(wǎng)絡(luò)學(xué)習(xí)函數(shù)為梯度下降動(dòng)量學(xué)習(xí)函數(shù)learngdm；性能函數(shù)采用均方誤差性能函數(shù)mse。

（二）插值精度的檢驗(yàn)

采用均方根誤差（Root mean square error，RMSE） 和平均相對(duì)誤差（）作為插值精度的衡量指標(biāo)。其計(jì)算公式如下：

（三）數(shù)據(jù)處理

采樣區(qū)選在榆樹市弓棚鎮(zhèn)13號(hào)村，農(nóng)田作物為玉米。首先，對(duì)試驗(yàn)地塊進(jìn)行DGPS網(wǎng)格定位，共計(jì)153 hm2，并在ArcGIS軟件上對(duì)其進(jìn)行網(wǎng)格劃分，網(wǎng)格大小為80 m×80 m。采用梅花五點(diǎn)采樣法，共采集386個(gè)樣品，采樣深度為20cm。將采集的樣品裝袋后進(jìn)行標(biāo)記，帶回實(shí)驗(yàn)室分析。

為了比較采樣密度對(duì)模型模擬精度的影響，將研究區(qū)域的386個(gè)采樣點(diǎn)利用ArcGIS軟件Geostatistical Analyst模塊的Create Subsets命令隨機(jī)提取100點(diǎn)、150點(diǎn)、200點(diǎn)、250點(diǎn)、300點(diǎn)、350點(diǎn)作為訓(xùn)練樣本，取樣后在SPSS下對(duì)訓(xùn)練樣本做正態(tài)分布檢驗(yàn)，不符合正態(tài)分布的重新取樣。在全部樣點(diǎn)數(shù)中隨機(jī)抽取50個(gè)作為檢驗(yàn)樣本，本文將對(duì)以上六種方案進(jìn)行比較研究。

三、結(jié)果與分析

（一）基于Kriging插值的農(nóng)田土壤速效鉀空間變異模擬

SPSS下的正態(tài)分布檢驗(yàn)結(jié)果表明，速效鉀在所有尺度下都符合正態(tài)分布。在GS+下對(duì)不同采樣尺度的速效鉀含量進(jìn)行半方差分析后，根據(jù)GS+所求出的最大變異距離和選用的模型作為參數(shù)，將訓(xùn)練樣本在ARCGIS下進(jìn)行Kriging插值。

圖1 不同采樣點(diǎn)下速效鉀的Kriging插值效果圖（單位：mg/kg）

（二）基于BP神經(jīng)網(wǎng)絡(luò)模型的農(nóng)田土壤速效鉀空間變異模擬

為了方便BP網(wǎng)絡(luò)模型的構(gòu)建和增強(qiáng)訓(xùn)練效果及驗(yàn)證網(wǎng)絡(luò)的泛化能力，需要對(duì)網(wǎng)絡(luò)的輸入和輸出進(jìn)行等價(jià)的歸一化預(yù)處理。本文采用歸一化處理函數(shù)mapminmax對(duì)網(wǎng)絡(luò)輸入和輸出進(jìn)行歸一化處理，并在得到模擬值之后使用mapminmax函數(shù)將模擬值還原到原來的量綱和量級(jí)。

利用仿真函數(shù)sim對(duì)通過訓(xùn)練的神經(jīng)網(wǎng)絡(luò)模型進(jìn)行模擬，得到檢驗(yàn)樣本土壤速效鉀的神經(jīng)網(wǎng)絡(luò)模擬值，見圖2。

圖2 不同采樣點(diǎn)數(shù)速效鉀的BP神經(jīng)網(wǎng)絡(luò)模擬效果（單位：mg/kg）

（三）模型比較

在全部樣本中隨機(jī)抽取50個(gè)樣本對(duì)Kriging模型和BP神經(jīng)網(wǎng)絡(luò)模型的插值結(jié)果進(jìn)行精度檢驗(yàn)，通過對(duì)比采樣密度的大小對(duì)模型模擬效果的影響，來分析兩種模型受采樣密度影響時(shí)模擬精度的變化規(guī)律，結(jié)果見表1。

表1 速效鉀在不同采樣密度下的精度對(duì)比Contrast of precision of soil available K the different sampling densities

檢驗(yàn)結(jié)果表明：在采樣點(diǎn)數(shù)較多的情況下，速效鉀的Kriging插值精度明顯高于BP神經(jīng)網(wǎng)絡(luò)的插值精度。隨著采樣點(diǎn)數(shù)的減少，Kriging插值精度顯著下降，而BP神經(jīng)網(wǎng)絡(luò)的插值精度相對(duì)穩(wěn)定。在采樣點(diǎn)數(shù)較少的情況下，BP神經(jīng)網(wǎng)絡(luò)的插值精度接近Kriging插值精度。

四、結(jié)論

在隨機(jī)采樣的情況下，衡量?jī)煞N模型在不同采樣點(diǎn)數(shù)下插值精度的變化，結(jié)果表明：Kriging模型插值精度隨采樣點(diǎn)數(shù)的減少顯著下降，而BP神經(jīng)網(wǎng)絡(luò)模型的插值精度隨采樣點(diǎn)數(shù)的變化不大；在采樣點(diǎn)數(shù)較多的情況下，Kriging模型插值精度顯著高于BP神經(jīng)網(wǎng)絡(luò)模型的插值精度，隨著采樣點(diǎn)數(shù)的減少，Kriging模型插值精度急劇下降，在采樣點(diǎn)較少的情況下，BP神經(jīng)網(wǎng)絡(luò)模型的插值精度已經(jīng)接近甚至優(yōu)于Kriging模型插值精度。

[1]王建軍,王長(zhǎng)松,殷朝珍.不同采樣密度對(duì)土壤養(yǎng)分空間插值結(jié)果的影響[J].江西農(nóng)業(yè)學(xué)報(bào)2010,22(7)∶87-90.

[2]MARTIN Y,ROSSITER D G,NJOMGANG R.Multi-scale characterization of soil variability within an agricultural landscape mosaic system in southern Cameroon[J].Geoderma,2005,12(5)∶117-143.

[3]YUAN X Y，CHAI X R，GAO R T，et al.Temporaland spatialvariability of soil organic matter in a county scale agricultural ecosystem[J].New Zealand JournalofAgriculturalResearch，2007，50(5)：1157-1168.

[4]LIU X M，ZHAO KL，XU J M，et al.Spatial variability of soil organic matter and nutrients in paddy fields at various scales in southeast China[J].Environmental Geology，2008，53(5)∶1139-1147.

[5]王淑英，路蘋，王建立.不同研究尺度下土壤有機(jī)質(zhì)和全氮的空間變異特征∶以北京市平谷區(qū)為例 [J].生態(tài)學(xué)報(bào)2008，28(10)：4957-4964.

[6]陳光，賀立源，詹向雯.耕地養(yǎng)分空間插值技術(shù)與合理采樣密度的比較研究[J].土壤通報(bào)，2008，39(5)∶1007-1011.

[7]潘瑜春，劉巧芹，閻波杰.采樣尺度對(duì)土壤養(yǎng)分空間變異分析的影響[J].土壤通報(bào)，2008，41(2)：257-262.

責(zé)任編輯：姚 旺

S158

B

1671-6531（2012）04-0042-02

吉林省科技廳科技支撐計(jì)劃項(xiàng)目 (20080207)；國(guó)家農(nóng)業(yè)科技成果轉(zhuǎn)化資金項(xiàng)目 (2009GB2B100095)；吉林師范大學(xué)研究生創(chuàng)新科研計(jì)劃項(xiàng)目 （201111）

于洋/吉林師范大學(xué)旅游與地理科學(xué)學(xué)院助教，碩士（吉林四平 136000）；劉吉平/吉林師范大學(xué)旅游與地理科學(xué)學(xué)院副教授（吉林四平 136000）。