專利名稱：基于mcmc的小麥品種特征參數(shù)估算方法
技術(shù)領(lǐng)域：
本發(fā)明屬于精準農(nóng)業(yè)中小麥品種特征定量評價領(lǐng)域，涉及小麥品種特征參數(shù)估算方法，具體涉及融合了馬爾可夫鏈蒙特卡洛方法(Markov Chain Monte Carlo, MCMC)的小麥品種特征參數(shù)估計方法。
背景技術(shù)：
隨著農(nóng)業(yè)技術(shù)的發(fā)展進步，新的小麥品種不斷出現(xiàn)，科學(xué)家通過對不同品種的生理生態(tài)特性研究，提出了不少可以量化作物品種特性的參數(shù)，如株高、千粒重、生理春化天數(shù)、灌漿時間、發(fā)生分蘗的能力等，這些參數(shù)被稱為作物品種特征參數(shù)，它主要反映作物基因型特征。近半個多世紀以來，國際上出現(xiàn)了不少作物生長模型，如美國的CERES模型、菲律賓國際水稻所的ORYZA模型、澳大利亞的APSIM模型、中國的WheatGrow模型等，它們通過解析“氣象_ 土壤_技術(shù)措施”與作物生理生態(tài)過程的機理關(guān)系，對作物的生長發(fā)育過程進行定量的表達，將反映基因型差異的作物品種特征參數(shù)融入到具體作物生長發(fā)育過程當(dāng)中，再結(jié)合環(huán)境因子作為輸入變量，達到對作物生長發(fā)育過程的動態(tài)模擬?？茖W(xué)家在廣泛利用模型進行作物生長預(yù)測的同時，也在考慮如何逆向地利用作物生長模型進行品種特征參數(shù)的反演，以定量評估作物的品種特性，從而指導(dǎo)作物育種。He等利用GLUE方法和 CERES-Maize來反演玉米的品種特征參數(shù)，金之慶等利用“試錯法”和CERES-wheat來反演小麥品種特征參數(shù)。但這些方法都存在反演效率低且無法搜尋到全局最優(yōu)變量等缺點。
發(fā)明內(nèi)容
本發(fā)明的目的在于克服上述缺陷，將MCMC方法和小麥生長模型(WheatGrow)相結(jié)合，通過構(gòu)建生育期和產(chǎn)量的目標函數(shù)，提出一種小麥品種特征參數(shù)估算方法，為定量研究小麥品種特性提供了一種新的途徑。本發(fā)明選用了馬爾可夫鏈蒙特卡洛(Markov Chain Monte Carlo, MCMC)方法，其最早被廣泛應(yīng)用于圖像處理、統(tǒng)計物理學(xué)領(lǐng)域。近年來，隨著計算機技術(shù)的迅猛發(fā)展，MCMC 方法在根系水分吸收模型、農(nóng)業(yè)生態(tài)模型等諸多復(fù)雜模型中得到了廣泛應(yīng)用。該方法是基于貝葉斯理論框架，首先建立平衡分布為η (χ)的馬爾可夫鏈，并對其平衡分布進行采樣， 然后通過不斷更新樣本信息而使馬爾可夫鏈能充分搜索模型參數(shù)空間，最終收斂于高概率密度區(qū)；該方法能將一些復(fù)雜的高維問題轉(zhuǎn)化為一系列簡單的低維問題，因此適用于復(fù)雜模型的貝葉斯計算。同時本發(fā)明選用了國內(nèi)建立并具有較好廣適性的WheatGrow生長模型，該模型涉及到9個小麥品種特征參數(shù)，分別是春化天數(shù)(VD)、光周期敏感因子(PS)、基本早熟因子(IE)、溫度敏感因子(TS)、灌漿期因子(FD)、葉熱間距(LT)、比葉面積(SLA)、收獲指數(shù)(HI)、小麥分蘗能力(TA)。其中，春化天數(shù)(VD)反映了小麥品種必須經(jīng)歷一定天數(shù)的春化，才能進行花芽分化，極冬性的小麥品種為60，極春性的小麥品種為O。光周期敏感因子(PS)，該參數(shù)反應(yīng)了小麥能否開花受到的日照長短的影響，用光周期敏感因子反映不同品種對日照長度的敏感性不同?；驹缡煲蜃?IE)，小麥出苗后，即使?jié)M足了溫光條件，也不會進行花芽分化，這段時間進行營養(yǎng)生長，不同品種表現(xiàn)不同，用基本早熟因子反映該生育階段的差異。溫度敏感因子(TS)，反映了小麥發(fā)育受溫度高低影響的程度。灌漿期因子(FD)，反映了小麥開花以后不同品種灌漿時間的長短不同。葉熱間距(LT)，反映了連續(xù)兩個葉原基出現(xiàn)的時間間隔。比葉面積(SLA)，單位重量的小麥葉片所包含的葉面積。收獲指數(shù)(HI)，反映了小麥干物質(zhì)向穗部分配的比例。小麥分蘗能力(TA)，不同小麥品種產(chǎn)生分蘗的能力不同，受品種基因控制，用分蘗能力反映小麥產(chǎn)生分蘗多少的特征。具體方案為一種基于MCMC的小麥品種特征參數(shù)估算方法，包括如下步驟1)小麥生長模型即WheatGrow模型的數(shù)據(jù)獲取通過自動氣象站記錄小麥生長季田間逐日氣象數(shù)據(jù)，包括日最高氣溫、日最低氣溫、日照時數(shù)、日降水量；利用五點取樣法， 將耕層土壤從上到下分4層進行田間土壤取樣，然后通過室內(nèi)物理化學(xué)分析，測定土壤容重、飽和含水量、田間持水量、萎焉含水量、實際含水量、土壤全氮、氨態(tài)氮和硝態(tài)氮；記錄小麥生長過程中的田間栽培管理數(shù)據(jù)；實測小麥主要生育期(包括播種期、拔節(jié)期、開花期、 成熟期)和收獲時的籽粒產(chǎn)量；2)在小麥品種特征參數(shù)范圍內(nèi)選擇一套初始品種特征參數(shù)組合θ to)，在未知參數(shù)分布的情況下，先假定各個參數(shù)的先驗分布q^rve1-1)為均勻分布，ef分別代表九個品種特征參數(shù)(i = 1,2,3,4,5,6,7,8,9)；3)在θ丨―1基礎(chǔ)上提出一個候選樣本θ:0"ew = θ廣1 + r χ (max(0,) - min(0i))/Dr是分布在0-1之間的一個隨機數(shù)；max ( θ J和min( θ J分別是參數(shù)θ ,的上下限；D等于5，控制著推薦步長的大?。?)將田間觀測的η年(η= 1，2，3...Ν)氣象數(shù)據(jù)、土壤數(shù)據(jù)、田間栽培管理數(shù)據(jù)以及兩套小麥品種特征參數(shù)θη 、θ (η)輸入到WheatGrow模型中，運行模型η遍，得到模擬的η年小麥主要生育期和產(chǎn)量結(jié)果；5)通過田間實測的η年小麥生育期及產(chǎn)量數(shù)據(jù)和模擬的η年小麥生育期及產(chǎn)量結(jié)果，計算生育期和產(chǎn)量的釋然函數(shù)π ( θ )，計算公式如下
權(quán)利要求
1. 一種基于MCMC的小麥品種特征參數(shù)估算方法，其特征在于，包括如下步驟1)小麥生長模型即WheatGrow模型的數(shù)據(jù)獲取通過自動氣象站記錄小麥生長季田間逐日氣象數(shù)據(jù)，包括日最高氣溫、日最低氣溫、日照時數(shù)、日降水量；利用五點取樣法，將耕層土壤從上到下分4層進行田間土壤取樣，然后通過室內(nèi)物理化學(xué)分析，測定土壤容重、飽和含水量、田間持水量、萎焉含水量、實際含水量、土壤全氮、氨態(tài)氮和硝態(tài)氮；記錄小麥生長過程中的田間栽培管理數(shù)據(jù)；實測小麥主要生育期(包括播種期、拔節(jié)期、開花期、成熟期)和收獲時的籽粒產(chǎn)量；2)在小麥品種特征參數(shù)范圍內(nèi)選擇一套初始品種特征參數(shù)組合θ((|)，在未知參數(shù)分布的情況下，先假定各個參數(shù)的先驗分布q(erve”)為均勻分布，θ〖分別代表九個品種特征參數(shù)(i = 1,2,3,4,5,6,7,8,9)；3)在θ”基礎(chǔ)上提出一個候選樣本θ廠 r是分布在0-1之間的一個隨機數(shù)；max( θ J和min( θ J分別是參數(shù)θ ,的上下限；D 等于5，控制著推薦步長的大?。?)將田間觀測η年(η=1，2，3...Ν)氣象數(shù)據(jù)、土壤數(shù)據(jù)、田間栽培管理數(shù)據(jù)以及兩套小麥品種特征參數(shù)θ ne\ θ (η)輸入到WheatGrow模型中，運行模型η遍，η = 1，2，3. . . N， 得到模擬的η年小麥主要生育期和產(chǎn)量結(jié)果；5)通過田間實測的η年小麥生育期及產(chǎn)量數(shù)據(jù)和模擬的η年小麥生育期及產(chǎn)量結(jié)果計算生育期和產(chǎn)量的釋然函數(shù)η ( θ )，計算公式如下 其中，Z(t)為觀察值，X(t)是模擬值，c0V(et)是et的協(xié)方差矩陣，θ為品種參數(shù)；6)按照如下公式計算釋然比a
2.根據(jù)權(quán)利要求1所述的基于MCMC的小麥品種特征參數(shù)估算方法，其特征在于，上述步驟2)所述的品種特征參數(shù)是與小麥生育期相關(guān)的小麥品種特征參數(shù)，包括了春化天數(shù) VD、對光周期反應(yīng)的敏感性PS、基本營養(yǎng)生長時間IE、對溫度高低的敏感性TS、灌漿時間的長短FD;與產(chǎn)量相關(guān)的特征參數(shù)，包括了相鄰兩個葉原基出現(xiàn)的時間間隔LT，單位重量的葉片的表面積SLA,品種發(fā)生分蘗的能力TA，光合干物質(zhì)向穗部分配的比例HI。
3.根據(jù)權(quán)利要求1所述的基于MCMC的小麥品種特征參數(shù)估算方法，其特征在于，上述步驟7)按照M-H準則，將計算出來的釋然比ap、ay和均勻分布在W，1]之間的隨機變量U 進行比較，對于生育期參數(shù)如果ap>U，設(shè)定θ〖 = θ廠(i = 1，2，3，4，5)，否則設(shè)定0〖二0廣1; 對于產(chǎn)量參數(shù)，如果ay彡U，設(shè)定時=0廠(i = 6，7，8，9)，否則設(shè)定0丨=0廣1。
4.根據(jù)權(quán)利要求1所述的基于MCMC的小麥品種特征參數(shù)估算方法，其特征在于，上述步驟12)馬爾可夫鏈收斂的條件是根據(jù)中心極限定理通過M-H準則的馬兒可夫鏈最終收斂于靜態(tài)分布，用一種定量收斂判斷指標GRi對終止計算進行判斷，該GRi指標的計算基于馬爾可夫鏈鏈內(nèi)和鏈間方差，如下式
全文摘要
本發(fā)明公開了一種基于MCMC的小麥品種特征參數(shù)估算方法，其方法主要是通過田間實測小麥品種的生育期及產(chǎn)量與小麥生長模型模擬的生育期及產(chǎn)量構(gòu)建概率密度函數(shù)，利用MCMC方法反演得到小麥品種特征參數(shù)。具體過程首先通過測試運算得到小麥品種特征參數(shù)的先驗概率分布，根據(jù)初始參數(shù)和先驗分布提出候選參數(shù)，計算產(chǎn)量和生育期的概率密度函數(shù)和釋然比，根據(jù)M-H準則判斷是否接受新的參數(shù)，最終得到了品種各特征參數(shù)的后驗概率密度分布。本發(fā)明所述方法估計的結(jié)果既準確又高效，在同類模型小麥品種特征參數(shù)估算當(dāng)中具有普遍適用性。
文檔編號G06F19/00GK102495948SQ201110368758
公開日2012年6月13日 申請日期2011年11月21日 優(yōu)先權(quán)日2011年11月21日
發(fā)明者倪軍, 劉小軍, 呂尊富, 姚霞, 曹衛(wèi)星, 朱艷, 湯亮 申請人:南京農(nóng)業(yè)大學(xué)