處方圖的生成過程

1、2精準(zhǔn)施肥處方圖的生成過程1、前言目前多采用空間插值法生成精準(zhǔn)施肥處方圖。插值從有限的樣本點數(shù)據(jù)來預(yù)測柵格數(shù)據(jù)中單元的值。它能夠用來預(yù)測任何地理點數(shù)據(jù)的未知值如高程、降雨量、化學(xué)物濃度、噪聲級等等。精準(zhǔn)施肥處方圖的作業(yè)是在ArcGISDesktop中進行，利用它內(nèi)嵌的空間分析模塊進行插值分析?？臻g分析模塊提供的插值函數(shù)有距離權(quán)重倒數(shù)函數(shù)（IDW）、樣條函數(shù)、克里格插值。距離權(quán)重倒數(shù)插值和樣條法插值是確定性的插值方法，它們是直接基于周圍已知點的值進行計算或是用指定的數(shù)學(xué)公式來決定輸出表面的平滑度的插值方法，因而不能派生出未知點的預(yù)測值，而精準(zhǔn)施肥處方圖的恰是要求派生出未知點的預(yù)測值，而利用克里格

2、插值法可以解決這個問題?？死锔癫逯挡逯捣▋?nèi)插方法公式如下，表達為數(shù)據(jù)的權(quán)重總和。其中，Z（si）是已測得的第i個位置的值；-是在第i個位置上測得值的未知的權(quán)重；s0是預(yù)測的位置；N是已知點（已測得值的點）的數(shù)目。利用克里格方法進行預(yù)測，必須完成以下兩個任務(wù)：（1）揭示相關(guān)性規(guī)則。（2）進行預(yù)測。要完成這兩項任務(wù)，克里格插值方法通過以下兩個步驟完成：（1）生成變異函數(shù)和協(xié)方差函數(shù)，用于估算單元值間的統(tǒng)計相關(guān)（也叫空間自相關(guān)），而變異函數(shù)和協(xié)方差函數(shù)也取決于自相關(guān)模型（擬合模型）。（2）預(yù)測未知點的值。因為前面已經(jīng)說過的兩個明確的任務(wù),因此要用克里格方法對數(shù)據(jù)進行兩次運算:第一次是估算這些數(shù)據(jù)的空

3、間自相關(guān)而第二次是做出預(yù)測?？臻g分析模塊提供了兩種克里格插值方法:普通克里格插值法和全局克里格插值法。普通克里格方法（OrdinaryKriging）是最普通和應(yīng)用最廣的克里格方法。它假設(shè)常數(shù)的均值是未知的。這是一個合理的假設(shè)除非你有一些科學(xué)的理由來否定這些假設(shè)。全局克里格方法假設(shè)數(shù)據(jù)中有主導(dǎo)趨勢（例如盛行風(fēng)），它可以用一個確定性的函數(shù)或多項式來模擬。從原始已知點中減去這一多項式，從隨機誤差中模擬自相關(guān)。在進行預(yù)測運算前，需要先完成從隨機誤差中擬和自相關(guān)的工作，然后將多項式加回到預(yù)測模型以獲得有意義的結(jié)果。全局克里格方法將僅用于知道數(shù)據(jù)的趨勢并能合理而科學(xué)地描述它的情況。2、處方圖的生成過程、

4、利用克里格插值生成一個表面Densij.Cwert單擊SpatialAnalyst模塊的下拉箭頭，選擇InterpolatetoRaster并單擊Kriging。單擊Inputpoints的下拉箭頭并單擊要加入的點數(shù)據(jù)集。單擊Zvaluefield的下拉箭頭并單擊要加入的字段。單擊你所需要的克里格方法。有兩種可供選擇的方法可變搜索半徑的克里格插值和固定搜索半徑的克里格插值單擊Semivariogrammodel下拉箭頭并單擊要使用的模型。單擊Searchradiustype的下拉箭頭并單擊Variable?？蛇x項，改變?nèi)笔〉妮斎朦c的數(shù)目?？蛇x項，指定一個最大距離。可選項，改變輸出數(shù)據(jù)的缺省單元尺寸。可選項，現(xiàn)在是否需要生成預(yù)測的標(biāo)準(zhǔn)誤差。鍵入輸出結(jié)果的名稱或使用缺省值在工作目錄下生成一個臨時文件。單擊0K。94、對生成的柵格數(shù)據(jù)重分類1、單擊SpatialAnalyst模塊的下拉箭頭，選擇重分類，在彈出窗口選擇分類2、修改重分類后的施肥量為合適值，確定后可生成施肥量柵格數(shù)據(jù)最終生成如下圖所示的柵格數(shù)據(jù)：3、單擊SpatialAnalyst模塊的下拉箭頭，選擇轉(zhuǎn)換，柵格到要素把