基于林地“一張圖”的小班自動編號探究

普紹林

(易門縣十街鄉林業工作站，云南 易門 651100)

ArcGIS軟件已廣泛應用于林業行業，在繪制好林業調查面層矢量數據后，一般按照從上到下，從左到右進行小班編號?，F階段小班編號的方法較多，也有許多小班編號工具，但對于極少數小班周界橫坐標最小值與縱坐標最大值相等時或者縱坐標相等時會產生小班編號錯誤，目前對小班編號不準確的研究相對較少。本文針對小班編號不準確問題結合小班排序編號實例進行探究，以某縣林地一張圖(Shapefile格式，以下簡稱一張圖)為數據，以ArcGIS 10.2浮動版為平臺，用Python語言編寫簡單的語句，求算小班周界折點橫坐標最小值與縱坐標最大值，利用查找相同的工具制作表比較橫坐標最小值、縱坐標最大值是否相等，采用Excel 2010函數countif自動進行小班編號，以VBScript語言編寫字符串函數left表達式格式化小班號進行探究。

1 查找相同的記錄工具的定義及用法

如果選定的輸入字段值與這些記錄值相同，則這些記錄為相同記錄。可對輸入數據集中多個字段的值進行比較。如果指定了多個字段，則記錄將首先按第一個字段中的值進行匹配，然后按第二個字段的值進行匹配，依此類推。

如果在輸出表中僅需要重復的記錄，則選中僅輸出重復記錄參數。如果未選中此參數(默認)，則輸出的記錄數與輸入數據集的記錄數相同。

輸出表含有2個字段：IN_FID 和 FEAT_SEQ。注意輸出表只能保存在地理數據庫中。

IN_FID 字段可用于將輸出表的記錄連接到輸入數據集。

相同的記錄具有相同的 FEAT_SEQ 值，而不相同的記錄具有順序值。FEAT_SEQ 值與輸入記錄 ID 無關。

2 Excel 2010函數countif的定義及用法

countif(要計數的單元格范圍，計算條件)：對區域中滿足某個指定條件的單元格進行計數。計算條件，其形式可以是數字、條件或文本。

用圖層屬性表Excel進行小班編號時，用Excel函數countif可實現按林班分組自動編小班號，編完一個林班再編另一個林班。在編號前按林班號升序、y降序、x升序、順序值升序等排序。小班自動編號時在編號(小班號)所在列第二單元格輸入公式：“=countif (c$2:c2,c2)”。輸入引號中的內容，回車。雙擊或下拉復制公式，統計第幾次重復次數，即為編號(小班號)自動編號的值。式中c$2為林班號所在列第2單元格，相對行引用，c2，c2為林班號所在列第2單元格。

華覺明:在國內,有一種將工藝美術代替傳統工藝的傾向。其實,工藝美術只是傳統工藝的一部分,在傳統工藝的十五個部類中,有八類是工藝美術。在陶瓷行業中,生產得更多的是日用品,土布生產出來在日常生活中使用,這些都不是藝術品。一直以來，由工藝美術行業代表中國參加國際手工藝理事會,這種情況應當有所改變。工藝美術行業制作的是帶有藝術性的產品,但不是手工藝產品的全部。

3 小班編號原理

面層小班按從上到下、從左到右編號，可按面層小班周界(折點)橫坐標最小值及縱坐標最大值排序編號。在Arcmap字段計算器窗口中用Python語言分別輸入x=“!Shape.extent.xmin!”，y=“!Shape.extent.ymax!”語句(引號中內容)，分別求算小班周界折點的橫坐標最小值(x)和縱坐標最大值(y)，再用查找相同的工具制作表查看表中是否存在橫坐標最小值與縱坐標最大值相等的情況。存在3種情況：(1)橫坐標最小值與縱坐標最大值不相等；(2)相鄰小班橫坐標最小值與縱坐標最大值相等；(3)不相鄰小班橫坐標最小值與縱坐標最大值相等。

橫坐標最小值與縱坐標最大值相等的相鄰小班是用同一個折點求算的橫坐標最小值、縱坐標最大值，橫坐標最小值與縱坐標最大值相等的不相鄰小班不是用同一個折點求算的。如果表中存在橫坐標最小值與縱坐標最大值相等的情況，排序時按橫坐標最小值、縱坐標最大值等排序，小班編號時相鄰幾個小班可能不會嚴格按照從左到右的順序編號。相鄰小班橫坐標最小值與縱坐標最大值相等時，可將圖層與表連接(連接字段選擇FID與IN_FID)，在橫坐標最小值與縱坐標最大值相等的相鄰小班屬性表順序值字段中輸入從左到右順序值，如1、2、3等，排序時用橫坐標最小值、縱坐標最大值、順序值字段等參與排序，再編小班號，可以解決極少數小班號編號未嚴格按照從左到右、從上到下編號的問題。不相鄰小班橫坐標最小值與縱坐標最大值相等時和不存在橫坐標最小值與縱坐標最大值相等時，不需要將圖層與表連接，不需要在順序值字段輸入順序值。

小班編號可采用導出圖層屬性表，將圖層屬性表排序后，用Excel函數countif將圖層屬性表Excel按林班自動分組編小班號，編完一個林班再編另一個林班，再將編好小班號的圖層屬性表Excel用FID與圖層(FID字段)連接，用這種方式編小班號字段類型為長整型，小班號如“1”，“10”，“100”等，可用VBScript 語言“字段計算器”中字符串函數left格式化小班號，如“00001”，“00010”，“00100”等。由于橫坐標最小值與縱坐標最大值相等的情況極少，也可以用橫坐標最小值、縱坐標最大值參與排序，在編好小班號后，用查找相同的工具制作表，打開表查看是否存在橫坐標最小值與縱坐標最大值相等的情況，如果存在橫坐標最小值與縱坐標最大值相等的情況(表中IN_FID與圖層中FID值相等)，打開圖層屬性表，很快可以找到橫坐標最小值與縱坐標最大值相等的相鄰小班，可標注好小班號后再手動修改小班號。

4 操作步驟

1)新建字段

打開Arcmap，加載一張圖，添加x(雙精度型，默認)字段，y(雙精度型，默認)字段，順序值(長整型，默認)字段，編號(長整型，默認)字段，原已有小班號(字符串)字段。

2)計算橫坐標最小值、縱坐標最大值

打開圖層屬性表，右鍵單擊x字段，打開字段計算器窗口，單擊選擇字段計算器，勾選Python，勾選整型，英文狀態輸入!Shape.extent.xmin!，確定，計算小班周界橫坐標最小值(圖1)。

同理，右鍵單擊y字段，輸入!Shape.extent.ymax!，計算小班周界縱坐標最大值。

3)用查找相同的工具制作表查看橫坐標最小值、縱坐標最大值是否相等

在數據管理工具常規中打開查找相同的工具，輸入數據集選擇一張圖，輸出數據集保存在地理數據庫中，選擇x、y字段，選擇僅輸出重復記錄(圖2)。

在內容列表中右鍵單擊查找相同的工具生成的表，選擇打開，查看是否存在橫坐標最小值與縱坐標最大值相等的情況(圖3，圖4)。

4)一張圖與表連接，順序值字段賦值

相鄰小班橫坐標最小值與縱坐標最大值相等時，用一張圖(FID字段)與用查找相同的工具生成的表(IN_FID字段)連接，再導出圖層。在導出的一張圖橫坐標最小值與縱坐標最大值相等的相鄰小班屬性表順序值字段中輸入順序值，在橫坐標最小值與縱坐標最大值相等的不相鄰小班屬性表順序值字段中不輸入順序值。

圖1 計算小班周界橫坐標最小值界面

圖2 查找相同的工具制作表界面

圖3 橫坐標最小值與縱坐標最大值不相等界面

圖4 橫坐標最小值與縱坐標最大值相等界面

5)在一張圖屬性表Excel表中用Excel 函數countif自動進行小班編號。

導出一張圖屬性表Excel，先排序，再用Excel 函數countif自動進行小班編號。

6)一張圖與一張圖屬性表Excel連接

用導出的一張圖(FID字段)與一張圖屬性表Excel(FID字段)連接，字段計器，林地一張圖.編號= [林地一張圖_TableToExcel$.編號]。

7)格式化小班號

本操作步驟中編號格式為“1”，“10”，“100”等，而實際林業調查小班號為文本型，如果小班號的固定位數是5位，那么小班號應為“00001”，“00010”，“00100”等。序號前面補的“0”位數是變化的?？捎肰BScript語言 “字段計算器”中字符串函數left實現，即可用字段計算器格式化小班號。其表達式為小班號=left("00000",5-len([編號]))& ([編號])，式中編號(長整型)，小班號(字符串)(圖5)。

8)手動修改小班號錯誤

由于小班號錯誤的情況極少，可手動修改小班號錯誤?？捎脵M坐標最小值、縱坐標最大值參與排序編小班號。編好小班號后，標注小班號，用查找相同的工具制作表查看是否存在橫坐標最小值與縱坐標最大值相等，如果存在相等，橫坐標最小值與縱坐標最大值相等的相鄰小班編號可能會出現錯誤，根據表IN_FID值和圖層FID值在圖層屬性表中選中，根據小班號標注判斷小班編號是否錯誤，再手動修改編號錯誤的小班號。

5 結論

用ArcGIS10.2 Python簡單的語句求算面層小班周界橫坐標最小值、縱坐標最大值，在查找相同的工具制作表中，圖層與表連接，在橫坐標最小值與縱坐標最大值相等的相鄰小班順序值字段中輸入從左到右順序值用于排序，導出圖層屬性表，用Excel函數countif自動進行小班編號，圖層與Excel表連接，用VBScript字符串函數left表達式格式化小班號，對橫坐標最小值與縱坐標最大值相等出現錯誤的小班編號，手動修改小班號錯誤。其操作簡單，容易理解，實用，可解決極少數小班號編號不準確的問題，使小班自動編號準確，符合林業調查規劃小班編號規則，在實際工作中是實用的辦法。

圖5 格式化小班號界面

6 討論

1)在Python窗口中輸入語法時，要注意空格縮進及字母大小寫等語法格式。

2)在小班數量超過7萬個時，用Excel編小班號容易造成Excel卡死。

3)要實現更快更便利操作，需進一步對Python軟件進行研討。