專利名稱：作物多元數輪作周期病蟲害的遙感診斷方法
技術領域：
本發明涉及農業及相關領域，用于在大地理區域中診斷所有與農作物輪作有關的病蟲害。
背景技術：
通常將在同一地塊上，一種作物種植一定的時間后，再種植另一種作物一定的時間，并且一直交替下去的過程稱為作物輪作，這兩種時間分別稱為從一種作物到另一種作物的輪作時間，也稱為輪作周期。輪作是作物主要的耕作栽培模式之一，具有明顯的生態效益和經濟效益，不僅可以大大改善土壤理化性狀和生態環境，而且還可以有效地破壞病蟲、雜草的生存環境，顯著降低農藥的使用量，減少農藥污染以及提高產量。我國不僅在南方地區水旱以及旱旱輪作廣泛應用，而且在北方地區輪作也廣泛應用，如玉米與大豆的輪作就具有明顯的生態和經濟效益，國家每年投資近億元用于補助實施大豆與玉米輪作的農戶。在世界范圍內，輪作同樣也正在成為提高農業生產的生態和經濟效益的關鍵技術，但是要提高輪作的生態效益和經濟效益，就必須研究與農作物輪作有關的病蟲害的發生情況，包括發生的區域和嚴重的程度，以便采取積極的措施進行治理和控制，因此，與農作物輪作有關的病蟲害的診斷就成為必須要解決的關鍵問題。
通常，將農作物不進行輪作或輪作周期太長所產生的病蟲害稱為農作物輪作病蟲害，由于農作物不進行輪作或輪作周期太長都會使輪作病蟲害加劇，所以及早發現或預測農作物輪作病蟲害，并采取相應的措施，對于提高農業生產的生態效益和經濟效益具有重要的意義。然而，如何在大的地理區域范圍內診斷農作物輪作病蟲害一直是農業領域中面臨的一個難題，這個問題的解決對于制定相應的防治策略，具有重要的價值，國內外的學者對此已進行了大量的探索，但沒有找到解決此問題的方法。
導致對農作物輪作病蟲害研究進展緩慢的主要原因之一是并不知道如何準確的表示和描述農作物輪作周期，而農作物輪作自身的特點已說明，傳統的用單一數值描述農作物之間的輪作周期方法是不科學的，首先農作物之間的輪作存在明顯的對稱性，其次這種農作物之間的輪作是相互的，第三在農作物輪作周期之間存在明顯的不對稱性，即它們數值的大小并不一定相等，而所有這三個關鍵特征根本就無法用單一數值來描述，采用單一數值很難對農作物輪作病蟲害進行客觀診斷，從而制約了用科學的方法對農作物輪作病蟲害進行深入研究，所以必須尋找新的方法。
本發明的目的是采用一種新的方法運用對遙感影像的解譯結果，并通過建立在統計學以及隨機過程基礎上的數學模型，來分析估算大地理區域范圍內的農作物的實際多元數輪作周期，然后再將農作物的實際輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式進行計算，用計算結果來診斷與農作物輪作有關的病蟲害的可能發生情況，以便采取相應的防治措施，該方法具有高效、簡單、易于應用等特點。
發明內容
本發明將從有限多個在年份上兩兩連續的衛星遙感影像上解譯出來的在大地理范圍內分布的輪作區域中的所有同一地塊上不同年份的農作物輪作的基本數據，根據地理范圍的大小不同，分別按市、縣、鄉或村為單位來計算農作物的輪作周期，再根據從一種作物到另一種作物輪作順序不同，他們的輪作周期也不盡相同的結論，運用在輪作為平穩隨機過程基礎上建立的公式來估算不同順序農作物的實際輪作周期，最終估算該農作物的實際多元數輪作周期。然后，再將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，來診斷與農作物輪作有關的病蟲害的可能發生情況。因此，將有限多個在年份上兩兩連續的衛星影像的解譯結果用于高精度估算大地理區域范圍內農作物實際多元數輪作周期以及將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式進行計算，用計算結果來診斷與農作物輪作有關的病蟲害的可能發生情況的方法成為本發明的重要特征。
本發明作物多元數輪作周期病蟲害的遙感診斷方法的技術方案是首先給出農作物多元數輪作周期的基本概念，然后對獲取的覆蓋要研究輪作區域的農作物的遙感影像進行解譯，獲得每個鄉或鎮(市或縣)的農作物輪作數據；再運用農作物輪作周期估算公式對該輪作數據進行分析，獲得農作物的實際輪作周期和實際多元數輪作周期；再將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式的計算結果即可診斷與農作物輪作有關的病蟲害的可能發生情況。
本發明的農作物遙感影像解譯主要包括四個步驟，首先對原始影像進行幾何校正，線性拉伸增強等預處理；其次根據野外實地調查確定輪作農作物以及其他作物的解譯標志；然后再通過人工目視解譯，從影像中提取輪作農作物的信息并生成含輪作農作物圖斑的解譯圖件，之后，對該圖件實地驗證檢查，對有問題的圖斑進行修正，使圖件的解譯準確可靠；最后生成包含輪作農作物的矢量化解譯圖件，并通過地理信息系統對解譯圖件進行空間分析處理，獲得輪作農作物的輪作數據。
研究農作物的多元數輪作周期，通常必須對輪作區域進行多年觀察，并且對所獲得的數據進行統計分析，才能獲得令人信服的結果，因此，本發明利用從有限多個在年份上兩兩連續的遙感衛星影像中獲得的數據來估算農作物的多元數輪作周期。有限多個在年份上兩兩連續的遙感衛星影像是指選取的多年遙感衛星影像在年份上每兩幅影像是連續的，但連續的每兩幅影像之間在年份上可以是間斷的。
本發明設計的農作物多元數輪作周期估算公式，除了利用有限多個在年份上兩兩連續的遙感衛星影像監測外，采用了對輪作區域內所有鄉鎮的農作物輪作水平同時觀察的方法，來估算一個理想化的典型鄉鎮的農作物輪作水平(對市、縣可同樣的討論)，相當于將鄉一級的隨機試驗樣本數量擴大到多個，以便更加準確地推算出該區域的農作物的多元數輪作周期，其依據為各個鄉鎮以及各個農戶之間在輪作及耕作栽培管理上具有明顯的相互獨立性，基本符合統計學上對樣本獨立性的要求；農作物輪作的周期性反映在用隨機過程描述的作物輪作過程時，該過程具有平穩性，在統計學上就相當于每年農作物輪作的數量應大致相同，可根據每年農作物輪作的數量估算出完成整個農作物輪作所需要的時間或周期；有限多個在年份上兩兩連續的多個鄉鎮的獨立樣本在數量上對于研究具有平穩性的鄉一級的農作物輪作過程的統計特征或周期來說，應該是足夠大。
本發明設計的診斷與農作物輪作有關的病蟲害的可能發生的方法適用于所有農作物是本發明的重要特征。
現在，首先以棉花與水稻輪作為例來說明農作物多元數輪作周期估算公式，但本公式適用于所有農作物多元數輪作周期的估算是本發明的重要特征。
在有限多個在年份上兩兩連續的遙感衛星影像監測的相鄰兩年中，用“棉/稻”表示上一年種棉花的地塊在下一年種了水稻的面積。對于給定的地塊，如果“棉i/稻i+1”表示在此地塊上，第i年種棉花而第i+1年種水稻的面積，則下一年種水稻占上一年種棉花面積的比例即“水稻的面積i+1/棉花的面積i”被定義為此地塊i年到i+1年的棉稻輪作因子或簡稱為此地塊i年的棉稻輪作因子，并且用CRRFi(Cotton-Rice Rotation Factor)來表示，CRRFi=NRAiCCAi;]]>其中CCAi為對于給定的地塊在第i年種棉花的面積；NRAi為CCAi中，第i+1年種水稻的面積。
假設僅有棉花與水稻在某一給定的地塊上進行輪作，對此地塊連續監測N年，這一輪作過程可以用具有明顯周期性的平穩隨機過程來描述，而且在此地塊上，棉稻輪作因子為CRRFi(i＝1，......，N-1)，則棉花與水稻的輪作周期CRRP(Cotton-Rice Rotation Period)為CRRP=limN→∞1N-1Σi=1N-11CRRFi;]]>由于實際上觀察的年數N是有限的，因此當N足夠大時，棉稻輪作周期CRRP的近似值CRRP‾=1N-1Σi=1N-11CRRFi.]]>
由于棉花與水稻輪作過程可以用具有明顯周期性的平穩隨機過程來描述，因此每年棉花輪作的數量在統計學上應該大致相同，而且棉花與水稻輪作周期是存在并可計算的。
因為僅有棉花與水稻在給定的地塊上進行輪作，所以根據CRRFi的定義，在給定“上一年棉花面積”CCAi的情況下，如果每年有“棉/稻”面積NRAi改種水稻，則 年可將“上一年棉花面積”改種完，故CRRPi＝1/CRRFi成立。
假設所要求的周期CRRP是隨機變量ξ的數學期望E(ξ)，那么近似確定CRRP的方法是對ξ進行N次重復抽樣，產生相互獨立的ξ值的序列ξ1，......，ξN(ξ1＝1/CRRF1，........，ξN＝1/CRRFN)，則其算術平均值ξN‾=1NΣi=1Nξi;]]>根據強大數定理，(limN→∞ξ‾N=CRRP)=1;]]>因此，當N充分大時，ξN‾≈E(ξ)=CRRP;]]>成立的概率等于1，因此可用 作為所求量CRRP的估算值。
根據上述定理，可以利用對輪作區域的M個鄉鎮連續N年衛星監測數據，計算該輪作區域的棉稻輪作周期，CRRP‾=1N-1Σi=1N-11CRRFi=1N-1Σi=1N-1(Σj=1MCCAij2NRAijΣj=1MCCAij);]]>其中M為鄉鎮總數；NRAij為第j個鄉鎮在第i年種棉花的面積中，在第i+1年種水稻的面積；CCAij為第j個鄉鎮在第i年種棉花的面積。
同理，對于水稻與棉花輪作的情況，可推出計算輪作區域的稻棉輪作周期公式RCRP‾=1N-1Σi=2N1RCRFi=1N-1Σi=2N(Σj=1MCCAij2NCAijΣj=1MCCAij);]]>其中M為鄉鎮總數；RCRFi為i年種棉花的面積中有多少在i-1年種水稻；NCAij為第j個鄉鎮在第i-1年種水稻的面積中，在第i年種棉花的面積；CCAij為第j個鄉鎮在第i年種棉花的面積。
推廣到一般情況，稱用上述方法求得的農作物輪作周期為農作物的統計學輪作周期，簡稱為農作物的統計輪作周期。因此，在大地理范圍內，作物x與y的輪作可以用它們之間的統計輪作周期來描述，根據上述分析，在一般情況下，作物x與作物y的輪作周期與作物y與作物x的輪作周期是不同的，這種在農作物統計輪作周期之間的不對稱性與通常農作物輪作周期之間的不對稱性的物理含義是相同的。因此，農作物x和y之間的輪作周期可以用下述m階多元數來描述(x1，...，xm，y1，...，ym，f1(x1，...，xm，y1，...，ym)，...，fp(x1，...，xm，y1，...，ym))其中xi(i＝1，...，m)為用m種不同的方法所求得的農作物x與農作物y的輪作周期；yi(i＝1，...，m)為用m種不同的方法所求得的農作物y與農作物x的輪作周期；fj(x1，...，xm，y1，...，ym)(j＝1，...，p)為用p種不同的方法所求的反映在輪作周期xi和yi(i＝1，...，m)之間的差異函數，因此可以將多元數輪作周期中的元分為三類，第一類與農作物x與y的輪作有關，第二類與農作物y與x的輪作有關，第三類與第一類和第二類元有關。
現在假定僅用上述稻棉輪作周期與棉稻輪作周期公式求得農作物x與y之間的輪作周期分別為x1和y1，描述x1和y1之間的周期差異函數為f1(x1，y1)＝y1-x1，則描述農作物x與y輪作的一階多元數輪作周期為(x1，y1，y1-x1)。通過與傳統的用單一數值描述的農作物x與y之間的輪作周期相比較，農作物多元數輪作周期提供了更加全面的信息，是一種更為科學的描述方法，并且可以得到許多重要的性質，以一階多元數輪作周期(x1，y1，y1-x1)為例如果y1-x1＞0，則說明輪作使農作物x的收益比作物y的收益更大，這是因為如果對于農作物x與y的輪作來說，輪作是使作物x具有更高的生態和經濟效益，那么農作物x與y輪作周期的長短通常與生態和經濟效益有密切的關系，農作物x與y的輪作周期短些以及農作物y與x的輪作周期長些，通常對農作物x更為有利，即對于兩個多元數輪作周期的比較，在x1相同的情況下，更要考慮差異函數值y1-x1的大小，即對y1值進行比較，y1值越大則輪作的效果越好。
綜上所述，以棉花與水稻輪作為例來說明農作物的一階多元數輪作周期估算公式如下，但本公式適用于所有農作物一階多元數輪作周期的估算是本發明的重要特征，對于農作物的m階多元數輪作周期可作類似分析。農作物x與y輪作的一階多元數輪作周期(x1，y1，y1-x1)其中x1=CRRP‾;y1=RCRP‾;y1-x1=CRRP‾-RCRP‾.]]>在以下分析中，如果不加以特別說明，所有多元數輪作周期均是指一階多元數輪作周期。
本發明設計的診斷農作物病蟲害的方法如下假定農作物i與農作物j進行輪作，用于防治農作物輪作病蟲害k，農作物輪作病蟲害k可能發生的嚴重程度與該農作物多元數輪作周期有關并且用指數αijk來表示αijk＝f(xa，xb，xc)；其中αijk＝f(xa，xb，xc)為求αijk值的函數；變量xa為農作物i的實際多元數輪作周期xa＝(x11，y11，y11-x11)；變量xb為農作物i的最小必須多元數輪作周期xb＝(x12，y12，y12-x12)；變量xc為與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關的量。
f(xa，xb，xc)可以有多種形式，例如可以用多元回歸方法求αijk值，此時有αijk＝f(xa，xb，xc)＝a11x11+b11y11+c11(y11-x11)+a12x12+b12y12+c12(y12-x12)；其中回歸系數a11，...，c12與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關。可以通過在農作物i的輪作區域，根據實際多元數輪作周期的不同，設置若干個樣方，再在每個樣方中測量輪作病蟲害k發生的嚴重程度αijk，并根據相對應的x11，...，(y12-x12)的值，用回歸的方法確定系數a11，...，c12的值。
通過上述多元回歸方程對樣方中的實際數據的處理，并且令xa＝(CRTPijk，RCTPijk，RCTPijk-CRTPijk)，xb＝(CRTNijk，RCTNijk，RCTNijk-CRTNijk)，回歸系數a11，...，c12用系數Cijk和Dijk來取代，可以得到求αijk值的更簡單的公式如下αijk＝Cijk×(CRTPijk-CRTNijk)+Dijk×(RCTPijk-RCTNijk)；其中(1)Cijk和Dijk為與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關的系數。可以在設置的若干個樣方中，調查相關農作物輪作病蟲害實際發生的嚴重程度αijk，再通過回歸的方法來確定系數Cijk和Dijk。
(2)αijk值越大則說明農作物輪作病蟲害k的發生將會越嚴重，αijk＝0說明農作物i與農作物j用于防治農作物輪作病蟲害k的多元數輪作周期符合最低多元數輪作周期的要求，而αijk＜0時|αijk|越大則說明農作物輪作病蟲害k越不容易發生，此時令αijk＝0。
(3)最小輪作多元數周期是指能夠最大限度的發揮農作物i和農作物j輪作的經濟和生態效益同時又能防治農作物輪作病蟲害k所必須的最小多元數輪作周期，可以通過科學試驗、測定，或根據長期積累的經驗來確定。
(4)1≤CRTPijk＜∞，1≤CRTNijk＜∞，1≤RCTPijk＜∞，1≤RCTNijk＜∞。
根據估算的農作物i的實際多元數輪作周期，計算相關的描述農作物輪作病蟲害可能發生程度的指數αijk，就可以判斷出在相關區域農作物輪作病蟲害k可能發生的嚴重程度。
具體實施例方式
實施例1江蘇省蘇北某市為全國著名的農業大市，位于里下河農區，總面積2393平方公里，人口155萬，轄45個鄉鎮，現有耕地13萬公頃，水域面積為全市總面積的四分之一，是國家商品糧和優質棉基地，年產糧食110萬噸、棉花4萬噸，與水稻、棉花同期生長的作物有玉米、大豆、山芋、蔬菜等。研究區域主要為棉花與水稻輪作，此外棉花也與其它作物輪作，但不是主流輪作模式。
研究采用的衛星影像為11937的陸地7號衛星影像和軌道號為11937的陸地5號衛星影像，覆蓋了上述研究區域，影像的每個像元或柵格的尺寸為30米×30米。在研究區域，對稻棉和棉稻輪作進行觀察的最佳時相為7月到8月。因此，我們選擇的四期TM影像的日期為2001年7月26日、2002年7月29日、2003年7月24日、2004年7月26日，影像質量符合解譯要求。
根據上述棉稻輪作周期公式CRRP‾=1N-1Σi=1N-11CRRFi=1N-1Σi=1N-1(Σj=1MCCAij2NRAijΣj=1MCCAij),]]>對遙感影像解譯獲得的數據進行分析，結果如表1所示。
表1蘇北某市2001年到2004年棉稻輪作水平遙感調查(面積單位公頃，周期單位年)
根據表1，2001年到2004年植棉面積波動較大，例如，從2001年到2002年，面積下降21.62％；從2002年到2003年，面積上升6.55％；而從2003年到2004年，面積繼續上升8.05％。同期，各鄉鎮的植棉面積波動也很大，面積增加的鄉鎮，2001年到2002年為10個；2002年到2003年為20個；2003年到2004年為25個。鄉鎮植棉面積增加的個數，基本上反映了蘇北某市植棉面積增減的總趨勢。從周期變化方式的特點來看，根據2001年到2002年的數據計算的棉稻輪作周期為2.01年；根據2002年到2003年的數據計算棉稻輪作周期為3.42年；根據2003年到2004年的數據計算的棉稻輪作周期為3.01年；差異較大，但基本上反映了從2001年以來，特別是2001年，蘇北某市各鄉鎮種植業調整對棉稻輪作的影響。從總體上來看，由連續4年對棉稻輪作周期監測的三個周期數據的算術平均值獲得的平均周期為2.81年，與我國明朝科學家徐光啟的2到3年(平均約為2.5年)的說法僅相差12.40％，比較接近。
根據上述稻棉輪作周期公式RCRP‾=1N-1Σi=2N1RCRFi=1N-1Σi=wN(Σj=1MCCAij2NCAijΣj=1MCCAij),]]>對遙感影像解譯獲得的數據進行分析，結果如表2所示。
表2蘇北某市2001年到2004年稻棉輪作水平遙感調查表(面積單位公頃，周期單位年)
根據表2，稻棉輪作平均周期為2.89年，與棉稻輪作平均周期2.81年十分接近，兩者僅相差2.85％，說明從棉花到水稻的輪作面積與從水稻到棉花的輪作面積基本相同，這一現象符合輪作的基本規律，也符合上述關于輪作過程是具有周期性以及平穩性的隨機過程的假設，并且說明本發明的輪作估算公式是合理的，遙感解譯的精度達到了試驗的要求。
綜上所述，棉花與水稻的一階多元數輪作周期計算如下因為(x1，y1，y1-x1)其中x1=CRRP‾;y1=RCRP‾;y1-x1=CRRP‾-RCRP‾.]]>所以(x1，y1，y1-x1)＝(2.81，2.89，2.89-2.81)＝(2.81，2.89，0.08)即蘇北某市棉花與水稻的多元數輪作周期為(2.81，2.89，0.08)。
在表1和表2的基礎上，可以計算蘇北某市各鄉鎮的多元數輪作周期如表3所示。
表3蘇北某市棉稻多元數輪作周期的遙感估算
由于在棉稻輪作周期和稻棉輪作周期的計算中，采用了不同的解譯面積，前者為當年種棉花下年種水稻的面積，而后者為上年種水稻當年種棉花的面積。如果用棉花將水稻輪作完所需要的時間y1明顯大于用水稻將棉花輪作完所需要的時間x1，那么這就說明在棉花改種水稻后的連作時間要比水稻改種棉花后連作的時間更長，因此可以更加充分利用水稻的浸泡作用殺死植棉土壤中殘存的病蟲害，以達到更好的輪作效果，因此對于兩個多元數輪作周期的比較，在x1相同的情況下，更要考慮差異函數值y1-x1的大小，即對y1值進行比較，y1值越大則輪作的效果越好，根據這一結論，可以對表3進行更加深入的分析，并可得出更多有用的結果。
實施例2棉稻輪作棉花病害遙感診斷表4蘇北某市棉稻輪作枯黃萎病害的遙感診斷
根據農作物輪作病蟲害的遙感診斷方法，對蘇北某市鄉鎮的棉稻輪作棉花枯黃萎病害遙感診斷的結果如表4所示，其中αijk為枯黃萎病的發病率(百分數)，(CRTNijk，RCTNijk，RCTNijk-CRTNijk)＝(2.85，2.89，0.04)為最小必須多元數輪作周期，(CRTPijk，RCTPijk，RCTPijk-CRTPijk)為實際多元數輪作周期，系數Cijk＝16.47，Dijk＝-1.47。注意在本表中，假設在實際多元數輪作周期的CRTPijk不大于最小必須多元數輪作周期的CRTNijk時枯黃萎病不會發生，實際上還可能有輕微的發生，但在這里不予考慮，根據試驗區域的情況枯黃萎病發病率在30％以上時就必須采取防治措施，因此在防治時應該對高于30％發病率的重點鄉鎮予以特別關注。
權利要求
1.本發明涉及作物多元數輪作周期病蟲害的遙感診斷方法，涉及農業及相關領域，診斷的對象為在大地理區域中所有與農作物輪作有關的病蟲害，用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期，再將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式的計算結果來診斷與農作物輪作有關的病蟲害的可能發生情況，以便采取相應的防治措施，該方法具有高效、簡單、易于應用等特點，本發明進一步涉及實現這種方法的技術。
2.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述診斷的對象為在大地理區域中所有與農作物輪作有關的病蟲害是指在同一地塊上，一種作物種植一定的時間后，再種植另一種作物一定的時間，并且一直交替下去的過程稱為作物輪作，所花費的時間稱為從一種作物到另一種作物的輪作時間，也稱為輪作周期，用年份對作物輪作周期研究的結論同樣適用于用其它時間單位對作物輪作周期的研究，在此基礎上，根據作物多元數輪作周期病蟲害的遙感診斷方法可以在大地理區域中用于診斷出所有與農作物輪作有關的病蟲害發生的地理位置和嚴重程度。
3.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指選取用于解譯的多年遙感影像在年份上每兩幅影像是連續的，但連續的每兩幅影像之間在年份上可以是間斷的。
4.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指對從遙感影像上解譯出的數據以相鄰的兩年為基本組合，每一組合包括前一年被輪作作物的面積以及在相同的地塊上前一年被輪作作物在下一年改種或輪作為另一種作物的面積。
5.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指農作物輪作可以用平穩隨機過程來描述，具有周期性并且是可以計算的，這里的周期性是在統計學上計算出來的所有參與輪作的一種作物被另一種農作物輪換種植所需要的年份數。
6.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指參與輪作的甲、乙兩種作物，甲與乙輪作與乙與甲輪作是兩個不同的概念，它們的多元數輪作周期可以不同。
7.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指相對于大地理區域，把自然鄉或村(市或縣)(或更小)符合獨立樣本要求的區域作為計算作物多元數輪作周期的基本單位，將大幅增加計算作物周期的獨立樣本數，減少對遙感影像數量的要求，并可提高估算作物多元數輪作周期的精度。
8.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期是指下述以C農作物與R農作物輪作為例所描述的數學公式、推導過程、計算結果以及應用方法適用于對所有農作物的多元數輪作周期及空間變化的估算，假定已經獲得了C農作物與R農作物輪作區域的M個鄉鎮連續N年遙感衛星影像的解譯數據，那么可用下式來計算C農作物與R農作物輪作周期， 其中M為鄉鎮總數；NRAij為第j個鄉鎮在第i年種C農作物的面積中，在第i+1年種R農作物的面積；CCAij為第j個鄉鎮在第i年種C農作物的面積，同理，假定已經獲得了R農作物與C農作物輪作區域的M個鄉鎮連續N年遙感衛星影像的解譯數據，那么可用下式來計算R農作物與C農作物輪作周期， 其中M為鄉鎮總數；NCAij為第j個鄉鎮在第i-1年種R農作物的面積中，在第i年種C農作物的面積；CCAij為第j個鄉鎮在第i年種C農作物的面積，推廣到一般情況，稱用上述方法求得的農作物輪作周期為農作物的統計學輪作周期，簡稱為農作物的統計輪作周期，因此，在大地理范圍內，作物x與y的輪作可以用它們之間的統計輪作周期來描述，根據上述分析，在一般情況下，作物x與作物y的輪作周期與作物y與作物x的輪作周期是不同的，這種在農作物統計輪作周期之間的不對稱性與通常農作物輪作周期之間的不對稱性的物理含義是相同的，因此，農作物x和y之間的輪作周期可以用下述m階多元數來描述(x1，…，xm，y1，…，ym，f1(x1，…，xm，y1，…，ym)，…，fp(x1，…，xm，y1，…，ym))其中xi(i＝1，…，m)為用m種不同的方法所求得的農作物x與農作物y的輪作周期；yi(i＝1，…，m)為用m種不同的方法所求得的農作物y與農作物x的輪作周期；fj(x1，…，xm，y1，…，ym)(j＝1，…，p)為用p種不同的方法所求的反映在輪作周期xi和yi(i＝1，…，m)之間的差異函數，因此可以將多元數輪作周期中的元分為三類，第一類與農作物x與y的輪作有關，第二類與農作物y與x的輪作有關，第三類與第一類和第二類元有關，現在假定僅用上述稻棉輪作周期與棉稻輪作周期公式求得農作物x與y之間的輪作周期分別為x1和y1，描述x1和y1之間的周期差異函數為f1(x1，y1)＝y1-x1，則描述農作物x與y輪作的一階多元數輪作周期為(x1，y1，y1-x1)。通過與傳統的用單一數值描述的農作物x與y之間的輪作周期相比較，農作物多元數輪作周期提供了更加全面的信息，是一種更為科學的描述方法，并且可以得到許多重要的性質，以一階多元數輪作周期(x1，y1，y1-x1)為例如果y1-x1＞0，則說明輪作使農作物x的收益比作物y的收益更大，這是因為如果對于農作物x與y的輪作來說，輪作是使作物x具有更高的生態和經濟效益，那么農作物x與y輪作周期的長短通常與生態和經濟效益有密切的關系，農作物x與y的輪作周期短些以及農作物y與x的輪作周期長些，通常對農作物x更為有利，即對于兩個多元數輪作周期的比較，在x1相同的情況下，更要考慮差異函數值y1-x1的大小，即對y1值進行比較，y1值越大則輪作的效果越好，綜上所述，以棉花與水稻輪作為例來說明農作物的一階多元數輪作周期估算公式如下，但本公式適用于所有農作物一階多元數輪作周期的估算是本發明的重要特征，對于農作物的m階多元數輪作周期可作類似分析。農作物x與y輪作的一階多元數輪作周期(x1，y1，y1-x1)其中
9.根據權利要求1所述的作物多元數輪作周期病蟲害的遙感診斷方法，其特征在于所述再將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式的計算結果來診斷與農作物輪作有關的病蟲害的可能發生情況是指下述數學公式、推導過程、計算結果以及應用方法適用于對所有農作物輪作病蟲害的遙感診斷，假定農作物i與農作物j進行輪作，用于防治農作物輪作病蟲害k，農作物輪作病蟲害k可能發生的嚴重程度與該農作物多元數輪作周期有關并且用指數αijk來表示αijk＝f(xa，xb，xc)；其中αijk＝f(xa，xb，xc)為求αijk值的函數；變量xa為農作物i的實際多元數輪作周期xa＝(x11，y11，y11-x11)；變量xb為農作物i的最小必須多元數輪作周期xb＝(x12，y12，y12-x12)；變量xc為與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關的量，f(xa，xb，xc)可以有多種形式，例如可以用多元回歸方法求αijk值，此時有αijk＝f(xa，xb，xc)＝a11x11+b11y11+c11(y11-x11)+a12x12+b12y12+c12(y12-x12)；其中回歸系數a11，…，c12與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關，可以通過在農作物i的輪作區域，根據實際多元數輪作周期的不同，設置若干個樣方，再在每個樣方中測量輪作病蟲害k發生的嚴重程度αijk，并根據相對應的x11，…，(y12-x12)的值，用回歸的方法確定系數a11，…，c12的值，通過上述多元回歸方程對樣方中的實際數據的處理，并且令xa＝(CRTPijk，RCTPijk，RCTPijk-CRTPijk)，xb＝(CRTNijk，RCTNijk，RCTNijk-CRTNijk)，回歸系數a11，…，c12用系數Cijk和Dijk來取代，可以得到求αijk值的更簡單的公式如下αijk＝Cijk×(CRTPijk-CRTNijk)+Dijk×(RCTPijk-RCTNijk)；其中(1)Cijk和Dijk為與農作物i、農作物j、農作物輪作病蟲害k、耕作栽培方式、輪作區域的地理條件以及其它因素有關的系數，可以在設置的若干個樣方中，調查相關農作物輪作病蟲害實際發生的嚴重程度αijk，再通過回歸的方法來確定系數Cijk和Dijk，(2)αijk值越大則說明農作物輪作病蟲害k的發生將會越嚴重，αijk＝0說明農作物i與農作物j用于防治農作物輪作病蟲害k的多元數輪作周期符合最低多元數輪作周期的要求，而αijk＜0時|αijk|越大則說明農作物輪作病蟲害k越不容易發生，此時令αijk＝0，(3)最小輪作多元數周期是指能夠最大限度的發揮農作物i和農作物j輪作的經濟和生態效益同時又能防治農作物輪作病蟲害k所必須的最小多元數輪作周期，可以通過科學試驗、測定，或根據長期積累的經驗來確定，(4)1≤CRTPijk＜∞，1≤CRTNijk＜∞，1≤RCTPijk＜∞，1≤RCTNijk＜∞，根據估算的農作物i的實際多元數輪作周期，計算相關的描述農作物輪作病蟲害可能發生程度的指數αijk，就可以判斷出在相關區域農作物輪作病蟲害k可能發生的嚴重程度。
全文摘要
本發明涉及作物多元數輪作周期病蟲害的遙感診斷方法，涉及農業及相關領域，診斷的對象為在大地理區域中所有與農作物輪作有關的病蟲害。用遙感影像中的農作物在空間上變化的解譯結果作為估算農作物多元數輪作周期的數據，通過建立的公式來估算農作物的實際多元數輪作周期，再將農作物的實際多元數輪作周期與最小必須的多元數輪作周期進行比較，根據相關公式的計算結果來診斷與農作物輪作有關的病蟲害的可能發生情況，以便采取相應的防治措施，該方法具有高效、簡單、應用廣泛等特點，本發明進一步涉及實現這種方法的技術。
文檔編號G06Q50/00GK1794290SQ20061003766
公開日2006年6月28日 申請日期2006年1月9日 優先權日2006年1月9日
發明者朱澤生, 孫玲, 朱犁 申請人:江蘇省農業科學院