1、地類計算與統計一、 數據準備。應用到的數據包括社界（dwg文件）和所求年份的現狀圖（shape或geodatabase）dwg文件的注記的插入點（text為左下角點）要落在對應的社界面上，這樣才能保證數據轉換后注記和面層的一一對應。二、 數據處理。數據的處理主要包括數據的轉換、點面屬性連接、數據相交三個部分。2.1數據轉換2.1.1建立數據集（1）點擊打開arccatalog程序，找合適的路徑然后【右鍵】新建personal geodatabase：（2）新建數據庫后，【雙擊】進入數據庫，【右鍵】新建數據集輸入名稱外，其余使用默認設置，不用修改，直接【下一步】，直到【完成】：輸入數據集名稱坐標

2、系統選unknown(未知),或者用【導入】選dwg文件的投影容限使用默認即可，點擊【完成】2.1.2將dwg數據導入數據集！首先必須確定dwg文件的路徑沒有中文名（d:pssjsj.dwg），否則導入會一片空白（1）【雙擊】進入數據集aaa，【右鍵】【導入】【要素類（多個）】（2）在【input features】中添加dwg文件的注記層（annotation）和面層（polygon）添加注記和面層，然后點【確定】導入2.1.3檢查修改面層的拓撲！面層可能存在裂縫和重疊錯誤，這樣會帶來計算面積的錯誤，因此要進行拓撲的檢查（1） 數據集aaa中，【右鍵】【新建】【拓撲】（2） 前面兩步使用默認

3、直接【下一步】，選擇要素的時候勾選面層然后【下一步】使用默認等級，然后【下一步】添加拓撲規則，規則選擇【不能重疊】和【不能有縫隙】，分兩次添加，然后【下一步】點擊【完成】，然后選【是】驗證拓撲（3） 修改拓撲錯誤打開arcmap，添加aaa_topology，即可看到拓撲檢查結果（紅色部分）可以看出，面層存在重疊，不存在裂縫，修要修改重疊部分，采用挖空的方法【編輯器】【開始編輯】用【選擇工具】（）選中重疊部分：【編輯器】【裁切】（clip，可以將與選擇部分有重疊的所有面擦除）直接【確定】，對所有重疊部分重復以上步驟在arcmap工具欄位置【右鍵】【拓撲】調出拓撲工具欄【驗證全部拓撲】重新驗證拓

4、撲看還有沒有拓撲錯誤修改完沒有拓撲問題后，【編輯器】【保存編輯】【停止編輯】拓撲錯誤已經消除，可以進行下一步操作！如果導入dwg文件的線層，然后用線層構面的話，可以省略掉拓撲檢查和修改這一步，操作會相對簡單些。但拓撲檢查是保證數據質量最常用的檢查方式，因此這里將這部分內容寫上，大家可以對比一下兩種方法，以達到一方面簡化操作步驟增加效率，一方面可以增加對拓撲的理解，保證數據質量的效果。！重疊部分的挖空和縫隙的填補操作是比較常用的編輯方式。挖空用【裁切】命令，填補縫隙則通過【修改要素】修改節點或者【自動完成多邊形】完成，在此不寫具體的操作步驟。2.2點面屬性連接!點面屬性連接是根據點落在面內部的特

5、性，通過空間連接運算，將點的屬性添加到面的屬性里面。具體來說就是標注點（annotation）的refname字段存儲了社名，要將社名這個屬性添加到對應的社界面上。2.2.1添加數據，檢查位置，保留字段（1）在arcmap中添加標注層（sj_dwg_annotation）和面層（sj_dwg_polygon） 原cad圖 對應arcmap圖檢查一下點的位置，確定每個面上對應唯一的注記點，排除異常情況。（2）保留字段保留字段是保留有用的字段參與空間連接，去除不需要的字段，使結果簡潔明了。打開相應圖層的【屬性】切換到【字段】這一欄對于注記層，點【全部清除】，然后勾選【refname】，保留refn

6、ame字段（存儲社名的字段）對于面層，點【全部清除】即可，不需要保留屬性字段2.2.2進行空間相交運算點擊工具欄的【arctoolbox】按鈕，打開arctoolbox工具箱按照【arctoolbox】【analysis tools】【overlay】【spatial join】（空間連接）打開空間連接對話框！早期計算地類采用【intersect】（相交）工具，一樣可以得到結果且有助于理解連接的過程，但操作相對復雜，應用【空間連接】工具操作比較簡單，對結果檢查比較容易對話框中，【目標要素】下拉選擇面層，【連接要素】下拉選擇注記層，【輸出要素類】則自己選擇路徑并為輸出的要素類命名（在此命名為sj

7、）。其余使用默認。點擊【確定】arcmap自動加載空間連接生成的圖層（sj），打開其【屬性表】找到相應圖層【右鍵】【打開屬性表】【refname】存儲社名，【shape_area】存儲面積，【join_count】字段是空間連接過程中自動生成的，表示連接到這個面的注記點的個數，如果為0則表示這個面上沒有點，如果為2或3則表示這個面上有2個或者3個點，這樣就不符合一一對應原則，要檢查修改。由上圖可知【join_count】全部為1，點與面一一對應，可以進行下一步運算。！到了這一步，社界面積就計算出來了，可以通過access統計每個社的面積。2.3圖層相交求地類面積（1）在上節打開的arcmap中

8、添加2007年的現狀地類shape數據（geodatabase數據也可，如果求2008年地類則添加2008年現狀圖，以此類推）（2）由圖上可以看出社界和現狀圖重疊在一起，相交運算是通過現狀圖的界線將村界切開，并且把現狀圖的屬性放到切開的小塊社界中（可以參看intersect工具的幫助）（3）打開相交工具對話框【arctoolbox】【analysis tools】【overlay】【intersect】（相交）在【輸入要素】中下拉選擇sj（社界）和07mzdl（07現狀圖），【輸出要素類】選擇路徑并輸入名稱（在此命名為sjdl）。其余使用默認，點擊【確定】開始相交運算。運算完成后會自動將sdd

9、l加載到arcmap中。上圖中紅色是原先sj層邊界，淡綠填充層是相交后生成的sjdl層，可以看出sjdl實際上就是07地類界線將sj層切開，或者sj層界線將07地類切開，亦即兩者的公共部分保留。（4）打開sjdl的圖層屬性，在字段對話框中保留【refname】（社名）、【dldm】（地類代碼）、【dlmc】(地類名稱)和【shape_area】（面積）四個字段（參見2.2.1）。（5）打開sjdl的屬性表（參見2.2.2）屬性表里面每一條記錄對應sjdl中的一個面，表明了這塊地所屬的社、地類名稱和代碼以及面積這些信息，比如說臘圃村村委這塊地被07現狀邊界切成了四塊，對應前面四條記錄。這個表是最

10、基本的表，后面的數據都將由這個表整理得出。（6）導出屬性表在屬性表界面中，右下角的【選項】菜單，【選項】【導出】，打開導出表對話框導出數據對話框中，【導出：】欄選all records（所有數據），輸出表里選路徑并命名dbf文件（這里命名為dlsj.dbf）,然后【確定】由上圖可以看出，導出的dbf文件包含兩個：dbf后綴和xml后綴文件。dbf后綴文件存儲所有數據；xml是元數據文檔，起說明作用，可以忽略。用excel【打開】dbf后綴文件，【另存為】xls文件（dlsj.xls）。至此，地類面積已經求出并導出為excel文件，后續統計以導出的excel文件為基礎。三、 地類數據統計（1）打

11、開access（【開始】【程序】下圖）(2)用access打開dlsj.xls(3)上圖中，點擊【打開】后，會出現導入對話框這個不修改，直接【下一步】這里要勾選【第一行包含列標題】，然后【完成】。這里直接點擊【完成】則使用默認名為表名，如果想修改表名則點擊【下一步】，然后再【完成】。如上圖所示，在數據庫的表視圖里可以看到導入的表（dlsj）雙擊表dlsj可以查看導入內容是否正確。（4）地類統計，合并相同地類的面積。由上圖可以看到，每個社都存在有同一地類有多塊地的情況，如村委有兩塊1110的地，布尾社有多塊1310的地。一般情況我們制表只需知道紅線范圍內每個社每個地類有多少面積，所以要進行統計將

12、同一個社的同一地類的所有地塊的面積加起來。點擊【查詢】雙擊【在設計視圖中創建查詢】添加要進行查詢的表，選中dlsj然后【添加】【關閉】在查詢對話框中，分別【雙擊】refname、dldm、dlmc和shape_area，添加四個字段到查詢字段中。在字段欄內部【右鍵】【總計】，調出【總計】選項在【總計】選項中，除了shape_area對應欄選【總計】外，其余欄選【分組】。其他選項不用修改！這個查詢實際上是將同一個分組的相同項進行總計，首先是合并同一個社的相同地類，然后合并同一地類的不同地塊面積。選項設立好后，點擊工具欄的【紅色！號】（運行），則得到查詢結果。查詢結果是每個社的每個地類的總面積【點

13、擊】左上角(紅色圈)全選數據，然后【右鍵】【復制】，將復制的數據粘貼到excel表中，可以用于制作各種表格。！理解查詢操作可以有助于各種面積的統計，比如說要統計紅線范圍內各個地類的面積，則在查詢字段里不選refname字段即可，得到結果如下：！整個地類計算過程中，要求前后面積一致，這樣方能保證計算的過程不會出錯。一般要確保下面等式：紅線面積 = 拓撲處理后的社界總面積 = 與現狀圖相交后的地類總面積 = 查詢后的地類總面積計算過程會有系統差，但一般前后面積相差不到1個平方米，如果范圍比較大可以將誤差適當放大，但是超過10個平方米則應當認為計算過程有錯誤，可以分環節進行檢查。在arcgis中查看

14、總面積：打開面層的屬性表，在屬性表的面積字段（shape_area）【右鍵】【統計】【總計】一項就是總的面積。在計算過程中多查看面積并進行對比可保證計算過程的準確性。arcgis注記轉化到cad進行空間運算后的結果有時要導出為cad文件，這里以導出社界畝數為例說明如何將arcgis的面層及標注導出到cad中去。打開arcmap，添加社界層（sj）【圖層】【右鍵】【屬性】，打開圖層屬性對話框【常規】【地圖】【分米】（此選項可以控制標注的大小，一般選分米即可）【確定】打開sj層屬性，在字段中只保留shape_area字段（參見地類計算與統計2.2.1）打開sj層屬性表，【選項】【添加字段】（添加一

15、個字段用來存儲畝數）如上圖設置，然后【確認】！常用的字段類型包括text（文本型），用來存儲文字；short integer 和long integer（短、長整型），用來存儲整數，只是長整型可存儲位數更多；float和double（單、雙精度型），用來存儲小數，double型可存儲更長位數。添加字段后如上圖所示在mu字段上【右鍵】【計算字段】，打開計算字段對話框設置如上圖所示，在等號下面輸入round(shape_area*0.0015,3)，亦即mu= round(shape_area*0.0015,3)，其中shape_area*0.0015是根據shape_area字段的平方數計算出畝

16、數，round函數用來對畝數保留三位小數。點擊【確定】，進行計算。計算結果如圖所示，arcgis的double類型默認將小數點前的0去掉，而cad注記則要把0顯示出來（如0.184畝）。在這里采取先導出excel表將0顯示，然后粘貼到新建的text字段中的辦法來解決。導出上圖的屬性表到excel中（參見地類計算與統計2.3之（6）【選項】【導出】，在excel表中打開如下圖：設置【mu】列的單元格格式，保留三位小數，然后選中數據部分【復制】在sj屬性表中新建一個text字段，【選項】【添加字段】，設置如下：點擊【確認】后，結果如下圖：將數據粘貼到mu1中：【編輯器】【開始編輯】，然后在mu1第

17、一格進行粘貼：可以看到，缺失的0已經補上。【編輯器】【保存編輯】【停止編輯】，關閉sj屬性表。打開sj層屬性對話框，對sj層進行標注：鼠標在sj層上【右鍵】【屬性】【標注】：先勾選【標注這個圖層中的要素】，然后點擊【表達式】在【表達式】中輸入：mu1 &畝，其中&是連接符號，將數字與字符“畝”連接起來。點擊【確定】，在圖層屬性對話框也點擊【確定】！不同的表達式可以取得不同的標注效果，比如說現狀圖中的（如），可以通過輸入：tbbh & vbnewline & & vbnewline & dldm 來完成。標注結果如上圖，此時要將圖層【縮放】一下，以調整標注文字的大小，一般以標注文字填滿最小面積的

18、多邊形為標準將圖層標注轉化為注記：選擇sj層，【右鍵】【將標注轉化為注記】記下注記存放的位置和名稱，然后點擊【轉換】 arcmap自動把轉化的注記層（sjanno）加載 數據庫中生成了sjanno要素類到這一步，注記要素在arcgis中已經生成，剩下的就是將其和面層導出為cad格式打開fme workbench（【程序】【fme】【fme workbench】）！fme是加拿大safe軟件公司開發的一款商用數據轉換軟件，在數據格式轉換方面較為常用，功能非常強大。根據向導進行數據轉換直接【ok】這一步是選擇源數據的格式選中esri geodatabase（mdb），然后【ok】格式選好后，【next】進行下一步這一步是對mdb進行定位，點擊“”瀏覽數據庫并將其添加，結果如圖，點擊【next】到下一步這一步是選擇目標格式選擇acad格式，然后【ok】，然后【next】這一步是對目標格式進行設置，比如版本、屬性等。點擊【settings】在release中選擇2000，其他不變，然后【ok】，然后【next】！這個選項規定了輸出的cad版本，默認是2007，由于cad