樊景超+丘耘+周國民

摘 要：為解決果樹整形試驗中不可逆、難以快速成型的缺點，以小冠疏層形蘋果樹為研究對象，采用馬爾科夫隨機過程模型建立蘋果樹枝條的隨機生長模型。在此基礎(chǔ)上，利用Visual C++與OpenGL圖形引擎實現(xiàn)了蘋果樹枝條生長過程及光照、陰影在降雨降雪等自然環(huán)境的三維可視化仿真。通過調(diào)控生長參數(shù)來模擬不同生長條件下的果樹枝條發(fā)育過程，為蘋果樹樹形培訓(xùn)和推廣提供了一種可行的技術(shù)手段。

關(guān)鍵詞：OpenGL；果樹培訓(xùn)；虛擬技術(shù)

中圖分類號：S661.1 文獻標識碼：A DOI 編碼：10.3969/j.issn.1006-6500.2015.08.018

我國是果樹種植大國[1]，果樹樹形對果樹的產(chǎn)量和品質(zhì)有著重要的影響[2-3]。科學(xué)高效地推廣先進的樹形成為提高我國蘋果果品質(zhì)量的重要手段。目前果樹樹形的培訓(xùn)和推廣工作，主要存在以下問題：一是口頭講解、教學(xué)光盤等形式不夠直觀形象[4]；二是果樹培訓(xùn)方面的農(nóng)技推廣人員相對短缺[4]；三是果樹整形具有不可逆性，錯誤的整形往往容易帶來不可彌補的損失[4]。虛擬植物技術(shù)通過模擬植物的生長過程，可以替代農(nóng)業(yè)生產(chǎn)過程中相對費時費力的試驗，目前常見的虛擬植物建模方法（如L系統(tǒng)、粒子系統(tǒng)、分形方法）和生成的軟件（如AMAP、Xfrog、Plant Studio）已經(jīng)在多個農(nóng)業(yè)領(lǐng)域展開虛擬植物研究工作。然而與小麥、玉米等一年生作物不同，果樹是一個多年生的、復(fù)雜的人工生態(tài)系統(tǒng)，傳統(tǒng)的基于植物生長機理的建模方法不適于果樹的虛擬生長研究。果樹枝條的分枝過程可以看作是一種開放性高且復(fù)雜的系統(tǒng)，同時亦是一種具有可控概率的系統(tǒng)[5]，非常適合使用馬爾科夫模型預(yù)測果樹的分枝模式[6-7]。

本研究在總結(jié)前人研究的基礎(chǔ)上，以小冠疏層形蘋果樹為對象，利用馬爾可夫模型來描述蘋果樹枝條的生長過程，結(jié)合OpenGL圖形引擎實現(xiàn)蘋果樹枝條生長過程以及果樹生長過程中自然環(huán)境的可視化，對我國果樹樹形的推廣進行有益地探索和嘗試。

1 系統(tǒng)功能設(shè)計

蘋果樹枝條虛擬生長系統(tǒng)擁有3個功能模塊：果樹枝條模擬模塊、自然環(huán)境虛擬模塊、生長控制模塊。蘋果枝條模擬包含了果樹枝干、新梢和葉片的模擬；自然環(huán)境模擬包含氣象數(shù)據(jù)顯示、光照和陰影顯示以及降雨、降雪的模擬；生長控制模塊包含果樹枝條生長參數(shù)的設(shè)置和枝條生長模擬功能。這些功能模塊彼此互相獨立，它們之間的功能組合構(gòu)成了蘋果樹枝條虛擬生長系統(tǒng)，系統(tǒng)功能結(jié)構(gòu)如圖1所示。

2 蘋果枝條模擬

2.1 蘋果枝干、新梢模擬

在OpenGL中，一般使用gluCylinder函數(shù)構(gòu)建一個圓柱體來近似模擬一段枝干，再使用貼圖技術(shù)模擬枝干的外皮。這種做法在普通林木的三維模擬中比較常見且符合人眼直觀感覺，然而經(jīng)過人工整形的果樹枝干多是彎曲生長的，因此不能直接使用該函數(shù)。本研究根據(jù)當前枝條的節(jié)作為自然分割點，將節(jié)間的枝干部分加入彎曲函數(shù)，達到整體枝干自然彎曲的效果。蘋果樹新梢相對筆直，使用gluCylinder函數(shù)直接繪制即可。

在果樹的生長周期內(nèi)，果樹枝條的長度和半徑根據(jù)蘋果物候日期發(fā)生響應(yīng)變化才能達到動態(tài)生長的視覺效果。本研究通過生長間隔設(shè)置枝條長度的動態(tài)增量來實現(xiàn)枝條的增長，而半徑的增長是通過在母枝條和子枝條之間設(shè)定一個響應(yīng)的衰減系數(shù)來同步改變的。隨著各個枝條長度和半徑的連續(xù)變化就能實現(xiàn)枝條形狀連續(xù)生長的模擬。

2.2 蘋果葉片模擬

葉片是果樹的重要功能器官，光合作用和蒸騰作用主要由葉片來完成。蘋果葉片在一年內(nèi)的生長過程中，形狀上呈現(xiàn)由小到大，在顏色上呈現(xiàn)由淺到深的變化，在位置上葉片的著生位置圍繞沿新梢呈螺旋狀生長。同時葉片在動態(tài)生長過程中應(yīng)該呈現(xiàn)一定的飄動效果，使果樹枝條的仿真更加逼真。由于果樹葉片數(shù)量相對較多，為達到實時刷新的目的，本研究使用OpenGL繪制葉片，使用幾何四邊形加上葉片的紋理貼，葉片圖像的四個角對應(yīng)四邊形的四個角。葉片圖像的背景色為白底，在OpenGL中使用glAlphaFunc函數(shù)將白色底色進行透明化處理，使得渲染后的圖片僅顯示葉片部分。該函數(shù)有兩個參數(shù)func和ref，func被允許的符合常量有GL_NEVER、GL_LESS、GL_EQUAL、GL_LEQUAL、GL_GREATER、GL_NOTEQUAL、GL_GEQUAL和 GL_ALWAYS，初始值為GL_ALWAYS。Ref則指明與引入的alpha值比較的參照值，該值屬于區(qū)間[0， 1]，初始值為0。在果樹枝條生長過程中，葉片根據(jù)當前物候日期加載響應(yīng)季節(jié)的蘋果葉片達到動態(tài)生長的效果。為使葉片更加逼真，在調(diào)用繪制葉片函數(shù)前使用glRotatef函數(shù)讓葉片圍繞新梢進行一定角度的隨機擺動，系統(tǒng)在實時繪制果樹葉片時，整株果樹的葉片呈現(xiàn)飄動效果。

3 自然環(huán)境模擬

蘋果園是一個復(fù)雜的人工系統(tǒng)，蘋果果樹生長是由氣象條件和栽培模式共同作用的結(jié)果。雨、雪和光照是植物生長不可或缺的自然環(huán)境，在虛擬場景中加入對自然環(huán)境的模擬，不僅可以增強果樹生長可視化的真實感，同時能夠更準確地重構(gòu)出貼近自然的虛擬場景。

3.1 氣象數(shù)據(jù)顯示

蘋果樹與周圍環(huán)境是一個矛盾的統(tǒng)一體，了解蘋果與環(huán)境之間的相互關(guān)系，特別是對生態(tài)環(huán)境條件的要求與反應(yīng)，有助于揭示其生態(tài)規(guī)律，指導(dǎo)引種、合理區(qū)劃，設(shè)計栽培模式、調(diào)控果園生態(tài)系統(tǒng)，保持相對平衡，以提高蘋果園質(zhì)量、產(chǎn)量和生物量，提高生態(tài)、經(jīng)濟和社會綜合效益。本研究使用蘋果樹生長所在地的年度氣象數(shù)據(jù)包括最高溫度、最低溫度、平均溫度、光照時間、空氣濕度、降雨量和氣象日期數(shù)據(jù)。通過使用OpenGL的DrawString函數(shù)將果樹模擬生長當日的氣象數(shù)據(jù)以文字形式顯示在系統(tǒng)的渲染區(qū)域，使用戶可以直觀看到果樹當前生長的物候時期和當日的氣象情況，同時系統(tǒng)根據(jù)當日的氣象情況動態(tài)控制光照和陰影以及降雨降雪的顯示。

3.2 光照和陰影模擬

光照是蘋果生長不可缺少的能量來源，是果樹進行光合作用的基礎(chǔ)。OpenGL所繪制的三維物體都需要提供光照才能渲染出逼真的材質(zhì)效果。本研究將光源設(shè)置在無窮遠處以模擬太陽相對于地球植物的距離。同時光源需用從所繪制果樹的左上方到右上方的一個拋物線運動來模擬一天內(nèi)太陽軌跡的變化。

在現(xiàn)實世界中，陰影無處不在。陰影提供了光線與物體位置的相關(guān)信息，因此陰影對人們理解三維場景的影響是非常巨大的。常用的陰影算法是以光源為視點對象對場景中的景物進行消隱的過程，使用Z-Buffer方法可以快速消隱。但是該方法的計算量隨著場景中三維模型的復(fù)雜度呈線性增長，所需存儲空間也會大大增加。因此本研究引用陳應(yīng)松[8]提出的三維模型投影到平面上的陰影生成方法：先把原物體投影到陰影平面上，然后再把所得的陰影繪制成黑色或其他半透明的暗色，其中把一個物體投影到陰影平面上的過程是利用矩陣變換交換來完成的。構(gòu)造適當?shù)木仃嚥⒊艘援斍暗哪Ｐ鸵晥D矩陣，將所有的頂點投射到陰影平面上，調(diào)用繪制物體的函數(shù)就可以完成其陰影的繪制。用此構(gòu)造的矩陣就可以把復(fù)雜的三維物體投射到二維平面上，從而簡化陰影的繪制工作。

3.3 降雪、降雨模擬

自然雨雪的模擬是計算機圖形學(xué)中具有挑戰(zhàn)性的問題之一[9]。粒子系統(tǒng)是模擬不規(guī)則物體最常用的圖形生產(chǎn)方法之一[10]。本研究使用粒子系統(tǒng)模擬雨雪降落過程。雨雪粒子系統(tǒng)由大量的雨雪粒子組成，對粒子系統(tǒng)主要考慮粒子數(shù)、粒子產(chǎn)生區(qū)域、粒子平均周期；對粒子則要考慮形狀、位置、速度、大小、生存期、顏色、透明度等屬性。圖2和圖3是使用粒子系統(tǒng)模擬降雪和降雨的系統(tǒng)截圖。普通粒子系統(tǒng)是將粒子貼圖之后以隨機角度渲染出來，這就導(dǎo)致粒子圖片容易與用戶的觀察角度成垂直狀態(tài)，從而導(dǎo)致觀察不到雨雪粒子，因此在本研究中粒子系統(tǒng)的圖片旋轉(zhuǎn)角度限定在一定的范圍內(nèi)，這樣方便用戶觀察渲染效果。

4 生長控制

4.1 生長參數(shù)設(shè)置

果樹虛擬生長的參數(shù)主要有主干根部半徑、主干頂部半徑、主干高度，他們用于描述所導(dǎo)入果樹的初始數(shù)據(jù)。衰減系數(shù)用于控制新生枝條與母枝條之間的直徑比例關(guān)系，分枝x軸角度、分枝y軸角度、分枝z軸角度用于控制新生枝條繞母枝條的旋轉(zhuǎn)角度，二級枝數(shù)量用于控制主干上二級枝條的數(shù)量，以符合人工修剪的需要，如圖4所示。

參數(shù)控制面板上有3個按鈕，其中“保存修改”按鈕用于把當前參數(shù)的修改保存到系統(tǒng)的內(nèi)存變量中，能及時觀察所生成果樹是否符合果樹生長的經(jīng)驗值；“同步到數(shù)據(jù)庫”按鈕將當前果樹生長參數(shù)數(shù)據(jù)保存到數(shù)據(jù)庫中永久保存，下次啟動系統(tǒng)時按照當前參數(shù)運行系統(tǒng)；“關(guān)閉”按鈕表示放棄當前修改。

4.2 虛擬生長

果樹枝條發(fā)育過程具有顯著的馬爾可夫鏈性質(zhì)，本研究根據(jù)蘋果樹設(shè)定的初始參數(shù)和生長參數(shù)，按照馬爾科夫[11-12]的分枝概率計算不同年月日的枝干生長變化情況，再由系統(tǒng)根據(jù)生長日期動態(tài)顯示果樹的虛擬生長情況。果樹虛擬生長過程分為3類：第1類是果樹跨年生長的對比截圖，從圖5中可以看到果樹的枝條形態(tài)結(jié)構(gòu)發(fā)生了較大變化；第2類是果樹逐日生長的對比截圖，從圖5中可以看到果樹的枝條和葉片緩慢生長的模擬情況；第3類顯示的為大齡果樹的年內(nèi)變化對比情況，從圖5中可以看到系統(tǒng)對多年生蘋果樹枝條的模擬效果良好，系統(tǒng)響應(yīng)迅速。

5 系統(tǒng)實現(xiàn)

本研究使用微軟公司的Visual Studio C++結(jié)合OpenGL進行系統(tǒng)實現(xiàn)。圖6是果樹生長模擬系統(tǒng)的主界面，系統(tǒng)分為菜單區(qū)和渲染區(qū)。菜單區(qū)展示系統(tǒng)所具有的功能，通過菜單來選擇并執(zhí)行系統(tǒng)具體功能。其主要功能是對果樹進行初始參數(shù)設(shè)置和導(dǎo)入果樹的氣象數(shù)據(jù)，再生成果樹枝條的模擬生長數(shù)據(jù)然后交由渲染區(qū)進行繪制。渲染區(qū)主要顯示4部分內(nèi)容：太陽、氣象數(shù)據(jù)、虛擬果樹以及果樹陰影以及雨雪。太陽和陰影與雨雪屬于互斥關(guān)系，不能同時出現(xiàn)，可根據(jù)當前虛擬生長日期的氣象數(shù)據(jù)進行交互繪制。

6 結(jié) 論

本研究是對果樹整形培訓(xùn)推廣的一個初步的嘗試[13-14]。通過參數(shù)化方法構(gòu)建果樹初始三維模型并利用馬爾科夫生長模型對果樹枝條的生長進行了動態(tài)模擬，從果樹枝條生長效果來看，效果逼真、成型快速，有效克服了傳統(tǒng)田間果樹培訓(xùn)中果樹生長不可逆的問題。

參考文獻：

[1] 周國民.我國數(shù)字果園的研究與發(fā)展[J].農(nóng)業(yè)網(wǎng)絡(luò)信息，2012（1）：10-12.

[2] 樊慶忠.樹形改造對紅富士蘋果生長和結(jié)果的影響[D].泰安：山東農(nóng)業(yè)大學(xué)，2005.

[3] 楊飛，盧桂賓，杜俊杰，等.樹形結(jié)構(gòu)對果樹生產(chǎn)能力的影響分析[J].山西林業(yè)科技，2012，41（1）：46-49.

[4] 王劍，周國民.基于虛擬現(xiàn)實技術(shù)的果樹修剪培訓(xùn)系統(tǒng)的研究[J].浙江農(nóng)業(yè)科學(xué)，2009（1）：213-216.

[5] 呂萌萌，郭新宇，陸聲鏈，等.基于Unity3D果樹交互虛擬修剪技術(shù)及其實現(xiàn)[J].農(nóng)機化研究，2015，4（4）：7-11.

[6] 夏寧，李保國，鄧西民.桃樹修剪分枝模式的模擬[J].植物學(xué)報，2004，40（7）：793-802.

[7] 熊瑛，張光年，郭新宇，等.基于馬兒可夫模型的蘋果樹枝條生長仿真[J].農(nóng)機化研究，2009（7）：70-73.

[8] 陳應(yīng)松，肖世德，胡漢春.用OpenGL繪制三維模型投影到平面的陰影方法[J].微計算機應(yīng)用，2003，24（6）：355-358.

[9] 趙春江，郭新宇，陸聲鏈，等.農(nóng)林植物生長系統(tǒng)虛擬設(shè)計與仿真[M].北京：科學(xué)出版社，2010.

[10] 羅維佳，都金康，謝順平.基于粒子系統(tǒng)的三維場地實時模擬[J].中國圖像圖形學(xué)報，2004，9（4）：495-500.

[11] Costers E， Guedon Y. Modelling the sylleptic branching on one-year-old trunks of apple cultivars[J].Journal of the American Society for Horticultural Science， 1997，122（1）：53-62.

[12] Costers E， Sinoquet J H， Kelner J， et al. Exploring within tree architeural development of two apple tree cultivars over 6 years[J].Annals of Botany，2003，91（1）：91-104.

[13] 趙靈芝，陳春園，趙勇，等.果樹立體生態(tài)種植模式的應(yīng)用[J]華北農(nóng)學(xué)報，2006，21（S）：32-36.

[14] 李雄，孫伯筠，李福榮，等.蘋果梨樹適宜負載量及樹體參數(shù)的研究[J].華北農(nóng)學(xué)報，1997，12（4）：111-115.