15TechnicalcodeofpracticeforproductivityevaluationofcultivatedI區(qū)農(nóng)牧業(yè)生態(tài)與資源保護(hù)中心、內(nèi)蒙古大學(xué)技術(shù)推廣中心、科右前旗農(nóng)牧業(yè)科學(xué)技術(shù)發(fā)展1黑土區(qū)耕地產(chǎn)能評估技術(shù)規(guī)程本文件規(guī)定了黑土區(qū)基礎(chǔ)條件準(zhǔn)備、耕地估產(chǎn)模型構(gòu)本文件適用于黑土區(qū)縣域尺度耕地產(chǎn)能評估，包括呼耕地產(chǎn)能cultivatedlandproductivity指在目前已有的生產(chǎn)條件和自然環(huán)境下單位面積耕地實(shí)際生產(chǎn)糧食的能力。耕地的產(chǎn)能水平是以4.1.1衛(wèi)星遙感影像4.1.1.1多光譜衛(wèi)星影像），宜時期，從中選擇三個不同時期的影像作為估產(chǎn)模型構(gòu)建。4.1.1.2影像質(zhì)量2耕地范圍內(nèi)影像無云遮擋，傳感器條帶噪聲為0,影像噪點(diǎn)<1%,波段無缺失，邊緣無畸變，空間產(chǎn)樣本地塊面積>1.33hm2,地塊規(guī)整。行混淆矩陣精度驗(yàn)證，精度驗(yàn)證指標(biāo)包括總體精度(OverallAccuracy)和Kappa系數(shù)(Kappa基于地塊的野外調(diào)查樣本點(diǎn)單產(chǎn)數(shù)據(jù)和不同時期衛(wèi)星影像數(shù)據(jù)進(jìn)行作物單產(chǎn)模型構(gòu)建，模型構(gòu)建單產(chǎn)數(shù)據(jù)作為被解釋變量，以代表作物生長期內(nèi)生長發(fā)育狀況的植被指數(shù)指標(biāo)(N波段組合而成)為解釋變量，建立多元回歸模型：NDVI?、NDVI?、NDVI?——分別為作物生長旺盛期內(nèi)的三個關(guān)鍵物候期的NDVI;3采用留一交叉驗(yàn)證法進(jìn)行模型精度評價(jià)。模型的預(yù)測精度采用穩(wěn)定性和預(yù)測能力2個評價(jià)指標(biāo)解釋，穩(wěn)定性用決定系數(shù)R2檢驗(yàn)(優(yōu)于0.5),預(yù)測能力用實(shí)測值與預(yù)測值的均方根誤差RMSE檢驗(yàn)。6耕地產(chǎn)能估算模型構(gòu)建6.1耕地地塊產(chǎn)量標(biāo)準(zhǔn)化將作物分類圖采用ArcGIS軟件掩膜功能做成掩膜數(shù)據(jù)以避免其他地物干擾，以地塊為單位對作物進(jìn)行產(chǎn)量估算，通過空間運(yùn)算得到地塊的產(chǎn)量。并以耕地地塊為基礎(chǔ)進(jìn)行分區(qū)統(tǒng)計(jì)，得到地塊的產(chǎn)量。對每一年的地塊產(chǎn)量進(jìn)行標(biāo)準(zhǔn)化：Ri——為該地塊某一年的產(chǎn)量標(biāo)準(zhǔn)化值；Y(i)——為該地塊產(chǎn)量；Y(i)max——為同一年中產(chǎn)量最大值；Y(i)min——為同一年中產(chǎn)量最小值。6.2耕地產(chǎn)能計(jì)算采用耕地多年產(chǎn)量均值表征耕地產(chǎn)能水平，均值高表示產(chǎn)能水平高，呈正相關(guān)；耕地產(chǎn)能的計(jì)算公FPC——為耕地產(chǎn)能；Ri——為該地塊某一年的產(chǎn)量標(biāo)準(zhǔn)化值；N為參與遙感估產(chǎn)的總年數(shù)。6.3耕地產(chǎn)能穩(wěn)定性計(jì)算采用耕地多年產(chǎn)量標(biāo)準(zhǔn)差反映耕地產(chǎn)能的穩(wěn)定性，標(biāo)準(zhǔn)差高表示年際變異度大，產(chǎn)能水平不穩(wěn)定，呈負(fù)相關(guān)。耕地產(chǎn)能的穩(wěn)定性評估公式：S——為耕地產(chǎn)能的穩(wěn)定性；Ri——為該地塊某一年的產(chǎn)量標(biāo)準(zhǔn)化值；R——為該地塊多年產(chǎn)量標(biāo)準(zhǔn)化后的均值；N——為參與遙感估產(chǎn)的總年數(shù)。7耕地產(chǎn)能評估7.1耕地產(chǎn)能分級4產(chǎn)能等級分級區(qū)間較低中等高穩(wěn)定、較穩(wěn)定、穩(wěn)定