王震明 李領寰 任佳倫 張超

（1浙江省浦江縣林業(yè)生態(tài)工程建設站 322200；2國家林業(yè)和草原局華東調查規(guī)劃設計院 浙江杭州 310019；3浙江農林大學園林設計院有限公司 浙江臨安 311300）

林業(yè)作業(yè)設計是林業(yè)工程項目建設的關鍵環(huán)節(jié)，直接關系到建設工程的質量和投資效果?；A資料（包括林分因子、面積、位置、土壤等屬性）的收集是造林作業(yè)設計中的第一個環(huán)節(jié)，也是造林作業(yè)設計的重要依據。但是，基礎林業(yè)數據的獲取并非易事，如：造林作業(yè)一般采用1∶10000地形圖進行實地區(qū)劃調查，其精度較低，難以精準定位，尤其是小班局部林分長勢不均勻時很難在圖紙上進行準確勾繪標記；若采用實時的高清衛(wèi)星影像進行調查，其昂貴的價格對于造林作業(yè)設計來說是難以承受的，另外受云霧影響，有時也難以獲取到合適的影像數據。由于難以對勘察現(xiàn)場的地形地貌以及林木信息進行精準的定位標示，導致造林作業(yè)設計在基礎資料的收集過程中存在較大誤差，大大降低了作業(yè)設計的質量和施工指導價值。

近年來，無人機遙感技術己成為森林資源調查與監(jiān)測的新途徑，為森林資源精準監(jiān)測以及精細化經營提供了新的技術手段[1]。無人機遙感技術憑借實時性、靈活性、實用性和低成本等優(yōu)勢，成為了快速獲取地理空間數據的有效平臺，其拍攝的高重疊度、多角度的高清影像，很大程度上滿足了獲取精細三維信息的需求。這種以無人機為飛行平臺，同時從垂直和傾斜角度獲取地物信息的技術，即為當前國內外熱門研究的低空攝影測量技術[2-3]。研究旨在將無人機低空攝影測量技術運用于造林作業(yè)設計工作中，通過軟、硬件集成，提出一套完整的技術方案和工作流程。該技術能夠真實還原三維造林場地現(xiàn)狀，顯著提高調查精度和效率，降低作業(yè)設計成本，對推動智慧林業(yè)建設具有重要意義。

1 研究區(qū)概況

浦江縣位于浙江省金華市西北部，地處東經 119°42′～120°07′，北緯 29°21′～29°41′之間，屬中亞熱帶常綠闊葉林北部亞地帶，浙閩山丘甜櫧、木荷林區(qū)。多年平均氣溫16.6℃，極端高溫40.1℃，極端低溫-11.1℃；年均無霜期228天；年平均降水量1412.2毫米；年日照時數 1996.2小時；年蒸發(fā)量1302.3毫米；年平均風速1.6米/秒，常年以東南偏東和東風為主。

研究區(qū)為浦江縣黃宅鎮(zhèn)杭金衢高速公路互通東側林地，以裸巖、宜林地、疏林地和火燒跡地為主，涉及鐘村67號和古塘村1號2個小班，小班面積393畝，工程實施面積260.5畝。

2 數據處理與研究方法

研究采用多旋翼無人機，利用傾斜攝影測量技術，結合差分GPS設備，開展精準造林作業(yè)設計工作。通過軟、硬件集成，提出一套完整的技術方案和工作流程。具體方法包括無人機數據獲取、無人機影像數據處理、制定施工設計方案及施工放線等。

2.1 無人機影像數據獲取

研究所使用的無人機為大疆，型號為PHANTOM4PRO，續(xù)航時間25min，巡航速度28-46km/h，搭載相機鏡頭為FOV84°，2000萬像素，鏡頭畸變＜1.5%，飛行預設航高為120m，航線航向重疊率為80%，旁向重疊率為70%，飛行范圍覆蓋1.15km2。拍攝時間為2017年6月3日上午9：30-11：30，在測量區(qū)域布設像控點，測定像控點三維坐標，使無人機按預設航線飛行；沿南北飛行了8條航線，共拍攝1042張影像，平均飛行速度 9m/s。航拍獲取的無人機影像地面分辨率達到厘米級。飛行時天氣為多云，地面風力2級，陣風2-3級。獲取影像后，篩除成像質量差、航高差大于10m、影像最大旋偏角和航帶彎曲度較大的航片。最后，通過RTK差分GPS采集標識點，將照片組導入Context Capture軟件，輸入控制點坐標對圖片地理信息進行校正，最后獲得TIF格式的高分辨率正射影像圖。

圖1 研究區(qū)空三點云模型

2.2 無人機影像數據處理

首先，將航拍影像、無人機記錄的POS（Position and Orientation System）數據和像控點導入Context Capture軟件，自動完成空中三角測量計算，生成精度報告；其次，根據地面像控點和POS數據，平臺校驗參數等進行區(qū)域網平差，解算多視影像外方位元素和加密點地面坐標，生成后綴為.LAS的高密度點云數據；最后，自動拼接數據，輸出正射影像（DOM）和數字表面模型（DSM）[4]。

通過8個控制點進行空三加密處理，如圖1所示，空三加密控制點檢查點誤差詳見表1。通過精度指標統(tǒng)計可以看出，在有構架航線、GPS輔助平差的情況下，8個控制點完全可以滿足精度要求，能大大減少外業(yè)控制點的需求；水平誤差值在0.004m-0.09m之間，垂直誤差值在-0.037m-0.046m之間，符合“地形圖航空攝影測量內業(yè)規(guī)范”（GB/T7930- 2008）中1∶1000地形圖測圖精度。

表1 空三加密控制點檢查點誤差統(tǒng)計表

2.3 作業(yè)設計方案確定

利用無人機傾斜攝影測量技術構建高精度三維模型，基于研究區(qū)現(xiàn)狀與種植區(qū)域坐標信息進行造林作業(yè)設計，根據不同類型的立體條件將研究區(qū)分為不同區(qū)塊，采取與之相適宜的種植措施，并反復模擬推敲造林設計方案，提高作業(yè)精度。

2.4 施工放線

使用ArcGIS10.2軟件輸出作業(yè)區(qū)塊主要拐點坐標，利用差分GPS設備完成野外施工放線工作。

3 研究結果與分析

3.1 研究區(qū)植被現(xiàn)狀

根據遙感影像判讀和實地調查數據，研究區(qū)受土壤、水源及人為因素的影響，山腳處林相較好，喬木郁郁蔥蔥；中下部以馬尾松疏林地為主，喬木層間有部分白櫟、黃檀等樹種，灌木層以格藥柃、白檀、白花龍、烏飯樹、杜鵑花為主，間有部分梔子、黃櫨、刺毛越桔、北江蕘花、長葉凍綠等，草本層以芒萁為主，間有蕨、光萼茅膏菜、闊鱗鱗毛蕨等；中上部及山頂由于土壤瘠薄，以灌木林為主，主要有山油麻、烏飯樹、算盤子、白檀、白花龍等，間有部分苦竹、長葉凍綠、格藥柃、梔子、北江蕘花、胡枝子、山礬、黃櫨等，草本層以芒萁、蕨為主，間有刺芒野古草、匙葉鼠麴草、知風草、青葙等；碎石裸巖區(qū)內以草本植物山類蘆為主，間有部分卷柏、刺芒野古草、蘇州薺苧等。

3.2 施工作業(yè)設計

將高精度DOM數據導入平板電腦，結合實地踏查，對造林小班地形地貌以及現(xiàn)有植被信息進行精準的定位標示，設計勾繪不同類型的作業(yè)區(qū)塊。根據土壤、地形及植被狀況，確定8種工程實施類型，如圖2所示，其中研究區(qū)類型F（局部）施工放線拐點分布圖見圖3。

（1）類型A。投影面積18.13畝，坡面面積20.48畝，平均坡度27.3°，平均海拔182米，位于山體中部，為馬尾松火燒跡地，存有大量馬尾松枯立木，立地條件較好。為防止水土流失，應盡快恢復森林景觀。設計將馬尾松枯立木清理后，使用米徑3cm泥球苗，按照64木荷+32香樟+16櫸樹+16楓香+16紅楠+16浙江樟的種植模式，即160株/畝（投影面積）的密度進行景觀林營造。

（2）類型B。投影面積44.96畝，坡面面積51.86畝，平均坡度28.6°，平均海拔191米，位于山體中上部，是立地條件較好的疏林地，地被以芒萁和蕨為主，零星或小組團分布馬尾松、白櫟、楓香、山油麻、烏飯樹、白花龍等喬灌木。設計使用地徑2cm容器苗，按照64木荷+32苦櫧+32青岡+32楓香的種植模式，即160株/畝的密度進行景觀林營造。

（3）類型C。投影面積47.14畝，坡面面積54.09畝，平均坡度29.9°，平均海拔205米，位于山體中上部，是立地條件較好的未成林造林地，以馬尾松、楓香、香樟、木荷等喬木幼樹為主，由于缺少撫育管護，現(xiàn)有密度為83株/畝。設計使用地徑2cm容器苗，按照30木荷+15苦櫧+15青岡+15楓香的種植模式，即75株/畝的密度進行補植造林。

（4）類型D。投影面積24.31畝，坡面面積29.26畝，平均坡度28.4°，平均海拔223米，位于山體上部，土壤瘠薄，立地條件較差，以芒萁、蕨、山類蘆、刺芒、野古草等草本植物為主，零星分布馬尾松。由于該區(qū)域景觀視覺敏感度高，立地條件差，因此設計選用耐干旱瘠薄、春季色彩鮮艷的映山紅和紅葉石楠進行景觀林營造，要求地徑2cm容器苗，按照210映山紅+90紅葉石楠的種植模式，即300株/畝的密度進行造林。

（5）類型E。投影面積8.22畝，坡面面積9.48畝，平均坡度28.6°，平均海拔231米，位于山體上部，土壤較瘠薄，苦竹、長葉凍綠、格藥柃等灌木以及馬尾松、楓香等喬木幼苗分布較多。設計使用地徑2cm容器苗，按照38木荷+22楓香+15苦櫧的種植模式，即75株/畝的密度進行補植造林。

（6）類型F。投影面積13.89畝，坡面面積14.51畝，平均坡度13.6°，平均海拔130米，位于山腳，是立地條件較好的造林失敗地，馬尾松幼樹分布較多，以格藥柃、白檀、白花龍、烏飯樹為主。由于周邊有公墓，且多次發(fā)生森林火災，不宜栽植易燃樹種，因此設計使用米徑3cm泥球苗，按照80木荷+48苦櫧+32香樟的種植模式，即160株/畝的密度進行防火景觀林營造。

（7）類型G。投影面積29.70畝，坡面面積33.62畝，平均坡度26.9°，平均海拔187米，位于山體中上部，是立地條件較好的馬尾松疏林地，密度約60株/畝。為提高現(xiàn)有林分防火能力，設計使用米徑3cm泥球苗，按照12木荷+6香樟+3楓香+3櫸樹+3紅楠+3浙江楠的種植模式，即30株/畝的密度進行補植造林。

（8）類型H。投影面積74.12畝，坡面面積94.63畝，平均坡度76.1°，原為馬尾松疏林地，發(fā)生過多次森林火災，以碎石裸巖為主，零星分布山類蘆。設計挑選局部有土壤的地塊，以馬棘、波斯菊和金雞菊種子為主，麻櫟、木荷、苦櫧和青岡種子為輔進行播種，混合比例為4∶4∶2，播種量為2g/m2。

營造林作業(yè)設計方案確定后，利用Arc GIS結合差分GPS設備進行造林施工放線工作。如圖3所示，該區(qū)域是類型F中的局部施工區(qū)域，使用ArcGIS將405號至415號共計11個主要拐點疊加到無人機正射影像圖上，生成含有坐標信息的PDF文件，將文件導入平板電腦，便于施工人員根據影像地物實現(xiàn)野外快速粗略定位，之后利用差分GPS設備根據坐標點開展精準施工放線工作。景觀林營造采用近自然混交的種植方式，同一樹種3-5株組團種植，株間距為1.5-4m，疏密相間；種植點遇有目的樹種或立地條件不佳時，應合理避讓，各地塊設計種植的苗木株數不得減少。

4 結論

研究以浦江縣黃宅鎮(zhèn)杭金衢高速公路互通東側林地造林作業(yè)設計為例，使用多旋翼無人機，利用低空攝影測量技術，結合差分GPS設備，開展精準造林作業(yè)設計工作。首先，使用多旋翼無人機獲取項目區(qū)高分辨率航拍影像，結合地面控制點，使用無人機低空攝影測量技術生成高精度三維場景模型，進而輸出數字正射影像（DOM）和數字表面模型（DSM）；其次，將高精度DOM數據導入平板電腦，結合實地踏查，對造林小班地形地貌以及現(xiàn)有植被信息進行精準的定位標示，設計勾繪不同類型的作業(yè)區(qū)塊；最后，使用三維場景模型反復推敲討論造林設計方案，方案確定后，使用Arc GIS軟件輸出作業(yè)區(qū)塊主要拐點坐標，利用差分GPS設備完成野外施工放線工作。研究將無人機技術與遙感技術、自動化控制技術和林業(yè)實地作業(yè)等方法結合起來互為補充[5]，有助于實現(xiàn)造林作業(yè)設計施工高效、精準的目標，為推動智慧林業(yè)的建設提供了參考。