陳晟 譚瑩 徐軍 胡進剛

（1 浙江省森林資源監(jiān)測中心 浙江杭州 310020；2 浙江省杭州市勘測設計研究院 310012）

森林資源規(guī)劃設計調(diào)查，即二類調(diào)查，是以國有林業(yè)局（場）、自然保護區(qū)、森林公園等森林經(jīng)營單位或縣級行政區(qū)域為調(diào)查單位，以滿足森林經(jīng)營方案、總體設計、林業(yè)區(qū)劃與規(guī)劃設計需要而進行的森林資源調(diào)查。

通用二類調(diào)查內(nèi)業(yè)流程中矢量化采集的圖斑面并不具備屬性信息，其屬性點矢量數(shù)據(jù)記錄了包含圖斑號等內(nèi)容的屬性信息，通過將這些信息賦值到圖斑面中從而形成二類圖斑面成果。在數(shù)據(jù)處理中，屬性點如果采取人工一一賦值的辦法會嚴重影響工作進度。因此，利用GIS 技術，分析屬性點信息的規(guī)律，自動對屬性信息進行編號，能夠大大減少工作步驟，技術人員只需在地圖上單擊即可完成一系列的圖斑屬性編號操作，極大的提高了數(shù)據(jù)采集的效率。

圖1 程序界面 Fig.1 Program interface

1 研究區(qū)概況

研究區(qū)玉環(huán)縣位于浙江省東南沿海黃金海岸線中段，東經(jīng) 121°05`~121°32`，北緯28°01`~28°19`?？h域總面積3419595 畝，其中陸域面積550518 畝，是浙江省第二大島，是浙江省六個海島縣之一。根據(jù)2009年玉環(huán)縣森林資源數(shù)據(jù)，全縣林地面積296245 畝，占土地總面積的53.18%，森林覆蓋率45.1%。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源

本文采用的數(shù)據(jù)主要包括：玉環(huán)縣地理國情普查圖斑數(shù)據(jù)，林業(yè)公益林圖斑面數(shù)據(jù)，“十一五”二類調(diào)查期間可參考的行政界線；玉環(huán)縣2.5m 分辨率衛(wèi)星影像（衛(wèi)片上經(jīng)過外業(yè)調(diào)查，已經(jīng)記錄小班編號）。

2.2 開發(fā)環(huán)境與技術路線

林地二類圖斑自動編號程序其開發(fā)環(huán)境使用了微軟公司可視化軟件Microsoft Visual Basic 6.0（VB），ArcGIS Engine 組件庫進行開發(fā)。VB 擁有圖形用戶界面和快速應用程序開發(fā)系統(tǒng)，可以輕易的連接數(shù)據(jù)庫，輕松的創(chuàng)建ActiveX 控件，程序員可以方便快速的建立一個應用程序。ArcGIS Engine（版本9.1）是ArcGIS 的眾多產(chǎn)品之一，它是用于構建定制應用的一個完整的嵌入式GIS 組件庫。開發(fā)者能夠利用組件庫中的接口、類等獨立進行封裝，構建具備特定功能的應用。ArcGIS Engine具有簡潔、靈活、易用、可移植性強等特點。

程序需要運行在Windows 操作系統(tǒng)下，軟件環(huán)境需要安裝Microsoft Visual Basic 6.0運行庫，而且需要安裝ESRI 公司ArcGIS Engine 產(chǎn)品運行庫，才能保證程序正確運行，程序界面如圖1所示。

2.3 算法流程

林地二類圖斑自動編號算法流程如圖2所示，具體步驟如下。

2.3.1 讀取行政代碼數(shù)據(jù)庫

進行二類圖斑自動編號之前，需要對程序進行一些設置，自動讀取鄉(xiāng)鎮(zhèn)、村行政代碼，以及該村最大小班號，便于操作人員選擇。具體如下：

①讀取行政村數(shù)據(jù)庫，并添加到VB 控件中便于選擇。

strSQL = "select * from county.dbf"

Set pRecordset = ExecuteSQL（App.path,strSQL,Msgstring）

Do While Not pRecordset.EOF

AddCundaima Trim（pRecordset（2）.Value）

pRecordset.MoveNext Loop

②初始小班號自動設定。操作人員依據(jù)遙感衛(wèi)片設定工作的村代碼、林班號，然后程序自動讀取該村的小班號作為工作的初始屬性信息。

strtemp = "村代碼 = '" ＆cundaima ＆"' AND " ＆" 林班號 = '" ＆linbanhao ＆"'"

pQueryfilter.WhereClause = strtemp 查詢條件

Set pFeaturecursor = pFeatureClass.Search （pQueryfilter,False）

Set pFeature = pFeaturecursor.NextFeature

If （pFeature.Value（pFeature.Fields.FindField（"

小班號"）） <901 And pFeature.Value （pFeature.Fields.FindField （"小班號"） >FindXiaobanhao MAX ThenFindXiaobanhaoMAX = pFeature.Value （pFeature.Fields.FindField（"小班號"））獲取該村該林班內(nèi)最大的小班編號信息，并記錄到窗口控件內(nèi)。

2.3.2 圖斑自動編號

經(jīng)過初始設定后，即可在衛(wèi)片上進行單擊，后臺自動存儲屬性信息，然后自動計算下一個圖斑的編號。

① 編號計算。自動讀取屬性點信息數(shù)據(jù)，并計算號碼為：最大小班號+1（即為下一個編入的小班號碼）

ReturnXiaobanhao（Str（FindXiaobanhaoMA X（ReturnCundaima（Trim（m_frmAttributeInput.c boCundaima.Text））,ReturnXiaobanhao（m_frmAttributeInput.txtLinB anhao.Text）） + 1））

② 編號保存到矢量數(shù)據(jù)中。

Set ppoint = New Point；ppoint.PutCoords mapx,mapy 記錄鼠標坐標點

Set pFeature.Shape = ppoint 設屬性點圖層

i = pFeature.Fields.FindField（"村代碼"）

pFeature.Value（i） = ReturnCundaima（Trim（m_frmAttributeInput.cbo Cundaima.Text））保存屬性點村代碼

j = pFeature.Fields.FindField（"小班號"）

pFeature.Value（j） = ReturnXiaobanhao（Trim（m_frmAttributeInput.tx tXiaobanhao.Text））保存屬性點小班號

k = pFeature.Fields.FindField（"林班號"）

pFeature.Value（k） = Trim（m_frmAttributeInput.txtLinBanhao.Text）

保存屬性點林班號信息

其中，pFeature 即為新增的屬性點矢量數(shù)據(jù)，記錄了該點所在的圖斑的屬性信息。

3 結果與分析

作者以玉環(huán)縣坎門鎮(zhèn)為試點，由于采用了本文的方法，加快了二類圖斑面的內(nèi)業(yè)處理作業(yè)，節(jié)省了人工，分析結果見表1。

表1 進度分析表

按照以往借助GIS 軟件人工進行圖斑面的勾繪即人工編號的方法，然后對村代碼、林班號、圖斑號幾個屬性字段逐一賦值的步驟，會產(chǎn)生大量的電腦操作，內(nèi)業(yè)工作量不僅大而且費時，操作人員不停的重復操作，容易造成重碼、漏碼等錯誤。然而，由表1中可以看出，應用本文研究成果不僅克服了步驟繁多，而且提升了圖斑面勾繪工藝步驟，效率是人工編號的一倍以上。并且，圖斑屬性信息只需鼠標點擊即可完成后臺存儲，大大減輕了作業(yè)人員的工作量，加快了內(nèi)業(yè)處理進度。。本文的研究成果應用到了玉環(huán)二類森林調(diào)查內(nèi)業(yè)工作中，各個工作組都得到了很好的應用，并有一致好評。

圖2 林地圖斑自動編號流程圖 Fig.2 Subcompartment numbering automatically flow chart

4 結束語

本文探索基于COM 技術的ArcEngine 組件庫二次開發(fā)，成功研發(fā)出ArcGIS 環(huán)境下林地圖自動編碼程序。本文該技術應用于浙江省“十一五”二類調(diào)查內(nèi)業(yè)處理，目前經(jīng)過程序升級，應用于玉環(huán)縣新一輪二類森林資源調(diào)查工作，如圖3。不僅大大節(jié)省了技術人員操作步驟，提高了效率，而且將復雜的人工操作步驟減少到最小，使得二類圖斑面的勾繪變得更加簡便、易用。

圖3 玉環(huán)縣小班編號應用 Fig.3 Yu Huan sumcompartment numbering automatically

綜上，利用ArcEngine 組件式擴展，可靈活 地把那些人工重復機械操作流程，無縫集成在ArcGIS 中，既降低開發(fā)難度，提高了開發(fā)效率，又增強業(yè)務系統(tǒng)的靈活性和開放性?；贑OM 技術的ArcEngine 組件庫二次開發(fā)方法將會得到各行各業(yè)的更廣泛應用，發(fā)揮出更大的功效。

[1] 浙江省林業(yè)廳.浙江省森林資源規(guī)劃設計調(diào)查技術操作細則[S].2004，6.

[2] 馮克忠，姜遵鋒，徐楊，等．ArcObjects 開發(fā)指南（VB 篇）[M]．北京：電子工業(yè)出版社，2007，1．

[3] 趙萬鋒，劉南，劉仁義，等．基于ArcObjects的系統(tǒng)開發(fā)技術剖析[J]．計算機應用研究，2004，21（3）：130-132．

[4] GB／T14911-2008．測繪基本術語[S]．

[5] GB／T20258．1-2007．基礎地理信息要素數(shù)據(jù)字典第1 部分1：500 1：1 000 1：2 000 基礎地理信息要素數(shù)據(jù)字典[S]．

[6] 鄔倫.地理信息系統(tǒng)原理方法和應用.北京：科學出版社，第八章空間分析：151.

[7] 湯國安.ArcGIS地理信息系統(tǒng)空間分析實驗教程.北京：科學出版社，第七章矢量數(shù)據(jù)的空間分布：196.

[8] ESRI ArcGIS Object Model Diagrams 9.1.

[9] Razavi A H．ArcGIS Developer's Guide for VBA [M]．New York：OnWord Press，2002．

[10] ZeilerM．Exploring ArcObjects[M]．California：Esri Press，2001．

[11] Chang K T．Programming ArcObjects with VBA：A Task-Oriented Approach [M]．London：Taylor ＆Francis,2004.