徐劍松，楊文總，李果明

(1.廣西壯族自治區(qū)林業(yè)勘測(cè)設(shè)計(jì)院,廣西 南寧 530011；2.廣西裕博林業(yè)勘測(cè)規(guī)劃設(shè)計(jì)有限公司,廣西 南寧 530011；3.廣西蒙山縣林業(yè)局,廣西 梧州 546700)

1 林業(yè)數(shù)據(jù)處理中合并圖斑的應(yīng)用場(chǎng)景

1.1 林權(quán)不動(dòng)產(chǎn)調(diào)查

針對(duì)位置彼此相鄰，且享有相同屬性值的圖斑。從幾何上看，他們享有共同邊，享有相同的屬性。按照林權(quán)不動(dòng)產(chǎn)調(diào)查的有關(guān)技術(shù)規(guī)程，幾個(gè)相鄰的不動(dòng)產(chǎn)單位，并且屬于同一林權(quán)權(quán)利人，應(yīng)劃為同一宗地[1]。這就要求將類圖斑按照此規(guī)則進(jìn)行合并。以往的林權(quán)調(diào)查，通常以小班作為調(diào)查單元，記載各小班的林權(quán)狀況，小班就相當(dāng)于不動(dòng)產(chǎn)調(diào)查中的不動(dòng)產(chǎn)單元。現(xiàn)有的林權(quán)數(shù)據(jù)大都為2009年林權(quán)制度改革的成果，林權(quán)證要逐步銜接到不動(dòng)產(chǎn)證，就要將不動(dòng)產(chǎn)單元?dú)w整到宗地[2]。

1.2 森林資源管理“一張圖”

在繪制森林資源管理“一張圖”小班后，會(huì)產(chǎn)生大量的碎斑(小面積圖斑)，需要對(duì)這些碎斑進(jìn)行合并操作[3]。雖然可以使用ArcGIS軟件中的消除工具進(jìn)行合并，但這個(gè)工具只能進(jìn)行有限的條件控制，往往會(huì)將一些本應(yīng)保留的圖斑也合并了。若要精準(zhǔn)合并圖斑，必須將此工具進(jìn)行功能拓展，而運(yùn)用C#語言進(jìn)行功能拓展是解決這一問難題的有效方案。

2 具體思路

2.1 找出具有鄰接關(guān)系的圖斑

根據(jù)圖斑的空間位置信息判斷圖斑的位置關(guān)系，通過循環(huán)語句(搜索游標(biāo))遍歷矢量數(shù)據(jù)的所有圖斑，再利用ArcGIS Engine提供的判斷圖斑之間位置關(guān)系的接口--IRelationalOperator，該接口支持判斷的空間關(guān)系類型包括Contains(包含)、Crosses(相交)、Equals(相等)、Overlaps(重疊)、Touches(相接)等[4]，采用該接口的Touch方法，輸入比對(duì)要素的幾何信息(IFeature.Shape)作為參數(shù)，可判斷兩個(gè)圖斑是否具有相鄰位置。

2.2 判斷鄰接圖斑是否符合合并條件

識(shí)別出當(dāng)前圖斑的所有相鄰圖斑后，遍歷讀取該圖斑屬性表相應(yīng)字段的值，通過比較運(yùn)算符確定該鄰接圖斑是否符合合并條件。使用流程控制語句中的分支語句，使符合合并的條件圖斑進(jìn)入下一處理流程。當(dāng)符合合并條件的相鄰圖斑為2個(gè)及以上時(shí)，推薦的做法是創(chuàng)建一個(gè)要素列表保存所有符合合并條件的相鄰圖斑，后續(xù)的合并、刪除要素都可以通過操作列表輕松實(shí)現(xiàn)。

2.3 合并前處理

經(jīng)過以上兩步的位置甄別和屬性值比對(duì)，已經(jīng)獲得符合合并條件的圖斑，合并操作需要包含在編輯業(yè)務(wù)流中，因此在執(zhí)行操作前，要先啟動(dòng)編輯工作流程，需要用到IEngineEditor接口，該接口用于開啟編輯對(duì)象的方法為StarEditing，需要傳入IMap和IWorkspace2個(gè)參數(shù)。

2.4 合并圖斑

ArcGIS Engine合并操作，相當(dāng)于用的ArcMap圖形操作界面的編輯器下的合并(Merge)功能，但又有所不同。合并圖斑用到的是ITopologicalOperator2接口，應(yīng)在合并前要預(yù)先創(chuàng)建一個(gè)空要素來保存結(jié)果圖斑，創(chuàng)建要素使用IFeatureClass接口的CreateFeature方法。調(diào)用ITopologicalOperator2接口的Union方法，將兩個(gè)同維度的幾何對(duì)象合并為單個(gè)幾何對(duì)象[5]，其參數(shù)為被合并圖斑的幾何信息，再使用賦值語句將結(jié)果賦值給預(yù)先創(chuàng)建的空要素。

2.5 刪除被合并的圖斑

由于圖斑合并使用的Union方法的功能是產(chǎn)生新的圖斑，再利用參與合并的圖斑邊界生成的新圖斑。該過程本質(zhì)上是新增了一個(gè)圖斑，對(duì)原來的圖斑未做處理。要達(dá)到ArcMap圖形操作界面的合并功能效果，還需要在合并圖斑后，把參與合并的圖斑刪除，否則要素類會(huì)存在重疊。刪除要素使用IFeature接口的Delete方法。

2.6 保存編輯

在編輯操作中，所有的編輯操作在沒有保存之前都被保存在內(nèi)存中，只有執(zhí)行保存操作，才能把所有的編輯結(jié)果保存到文件中[6]。停止編輯、保存編輯內(nèi)容的方法均在IWorkspaceEdit接口下，停止編輯方法為StopEditOperation；保存編輯的方法為StopEditing，該方法需要傳入一個(gè)參數(shù)，“true”表示保存，“false”表示不保存。

3 程序編寫

以Microsoft的可視化開發(fā)環(huán)境Visual Studio作為開發(fā)工具，采用面向?qū)ο蟮能浖_發(fā)語言C#創(chuàng)建一個(gè)項(xiàng)目，并在擴(kuò)展項(xiàng)引入ArcObjects。

關(guān)鍵代碼：

pWorkspaceEdit.StartEditing(true);

pWorkspaceEdit.StartEditOperation();

for (inti = 0; i

{

IFeatureCursorpFeatureCursor = pFeatureClass.Search(null, false);

IFeaturepFeature = pFeatureCursor.NextFeature();

IGeometryCollectiongeometryBag = newGeometryBagClass();

ITopologicalOperator2 topoOperator = null;

IGeometrygeometryMax = null;

IGeometrygeometryOther = null;

if (LstOther.Count> 0)

{

geometryMax = maxFeature.Shape;

IGeometrygeometryNext = null;

IFeatureunionedFeature = pFeatureClass.CreateFeature();

IGeometryunionedGeometry=null;

for (int j = 0; j

{

geometryNext = LstOther[j].Shape;

if (unionedGeometry != null)

{topoOperator = (ITopologicalOperator2)unionedGeometry;

}

unionedGeometry = topoOperator.Union(geometryNext);

pWorkspaceEdit.StopEditOperation();

pWorkspaceEdit.StopEditing(true);

4 結(jié)論與建議

借助C#語言編程及ArcGIS Engine提供的組件，可以讓圖斑之間的合并操作變得簡(jiǎn)單、智能，并可廣泛應(yīng)用于林業(yè)數(shù)據(jù)處理工作中。對(duì)需要合并的600個(gè)圖斑進(jìn)行測(cè)試，使用本文所述方法，耗時(shí)3 min，正確率100%；而采用傳統(tǒng)的手動(dòng)合并方法操作，則需要2 h，且圖斑數(shù)量增多后，容易出現(xiàn)合并錯(cuò)誤。

合并圖斑有許多需要注意的細(xì)節(jié)，建議合并前預(yù)先保存被合并圖斑的屬性值，完成合并后再對(duì)結(jié)果圖斑賦值，使結(jié)果圖斑屬性與合并前某一圖斑屬性一致。若不設(shè)置屬性值，結(jié)果圖斑將只保留圖形，屬性值為空。另外，在合并后要加上保存的操作，才能成功合并圖斑。設(shè)置完屬性值后，也要使用store方法保存。