利用ARCGIS進(jìn)行地類計算與統(tǒng)計

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

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

3、直接【下一步】，選擇要素的時候勾選面層然后【下一步】使用默認(rèn)等級，然后【下一步】添加拓?fù)湟?guī)則，規(guī)則選擇【不能重疊】和【不能有縫隙】，分兩次添加，然后【下一步】點擊【完成】，然后選【是】驗證拓?fù)洌?） 修改拓?fù)溴e誤打開arcmap，添加aaa_topology，即可看到拓?fù)錂z查結(jié)果（紅色部分）可以看出，面層存在重疊，不存在裂縫，修要修改重疊部分，采用挖空的方法【編輯器】【開始編輯】用【選擇工具】（）選中重疊部分：【編輯器】【裁切】（clip，可以將與選擇部分有重疊的所有面擦除）直接【確定】，對所有重疊部分重復(fù)以上步驟在arcmap工具欄位置【右鍵】【拓?fù)洹空{(diào)出拓?fù)涔ぞ邫凇掘炞C全部拓?fù)洹恐匦买炞C拓

4、撲看還有沒有拓?fù)溴e誤修改完沒有拓?fù)鋯栴}后，【編輯器】【保存編輯】【停止編輯】拓?fù)溴e誤已經(jīng)消除，可以進(jìn)行下一步操作！如果導(dǎo)入dwg文件的線層，然后用線層構(gòu)面的話，可以省略掉拓?fù)錂z查和修改這一步，操作會相對簡單些。但拓?fù)錂z查是保證數(shù)據(jù)質(zhì)量最常用的檢查方式，因此這里將這部分內(nèi)容寫上，大家可以對比一下兩種方法，以達(dá)到一方面簡化操作步驟增加效率，一方面可以增加對拓?fù)涞睦斫?，保證數(shù)據(jù)質(zhì)量的效果。！重疊部分的挖空和縫隙的填補(bǔ)操作是比較常用的編輯方式。挖空用【裁切】命令，填補(bǔ)縫隙則通過【修改要素】修改節(jié)點或者【自動完成多邊形】完成，在此不寫具體的操作步驟。2.2點面屬性連接!點面屬性連接是根據(jù)點落在面內(nèi)部的特

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

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

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

8、添加2007年的現(xiàn)狀地類shape數(shù)據(jù)（geodatabase數(shù)據(jù)也可，如果求2008年地類則添加2008年現(xiàn)狀圖，以此類推）（2）由圖上可以看出社界和現(xiàn)狀圖重疊在一起，相交運(yùn)算是通過現(xiàn)狀圖的界線將村界切開，并且把現(xiàn)狀圖的屬性放到切開的小塊社界中（可以參看intersect工具的幫助）（3）打開相交工具對話框【arctoolbox】【analysis tools】【overlay】【intersect】（相交）在【輸入要素】中下拉選擇sj（社界）和07mzdl（07現(xiàn)狀圖），【輸出要素類】選擇路徑并輸入名稱（在此命名為sjdl）。其余使用默認(rèn)，點擊【確定】開始相交運(yùn)算。運(yùn)算完成后會自動將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）屬性表里面每一條記錄對應(yīng)sjdl中的一個面，表明了這塊地所屬的社、地類名稱和代碼以及面積這些信息，比如說臘圃村村委這塊地被07現(xiàn)狀邊界切成了四塊，對應(yīng)前面四條記錄。這個表是最

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

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

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

13、擊】左上角(紅色圈)全選數(shù)據(jù)，然后【右鍵】【復(fù)制】，將復(fù)制的數(shù)據(jù)粘貼到excel表中，可以用于制作各種表格。！理解查詢操作可以有助于各種面積的統(tǒng)計，比如說要統(tǒng)計紅線范圍內(nèi)各個地類的面積，則在查詢字段里不選refname字段即可，得到結(jié)果如下：！整個地類計算過程中，要求前后面積一致，這樣方能保證計算的過程不會出錯。一般要確保下面等式：紅線面積 = 拓?fù)涮幚砗蟮纳缃缈偯娣e = 與現(xiàn)狀圖相交后的地類總面積 = 查詢后的地類總面積計算過程會有系統(tǒng)差，但一般前后面積相差不到1個平方米，如果范圍比較大可以將誤差適當(dāng)放大，但是超過10個平方米則應(yīng)當(dāng)認(rèn)為計算過程有錯誤，可以分環(huán)節(jié)進(jìn)行檢查。在arcgis中查看

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

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

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

17、一格進(jìn)行粘貼：可以看到，缺失的0已經(jīng)補(bǔ)上?！揪庉嬈鳌俊颈４婢庉嫛俊就Ｖ咕庉嫛?，關(guān)閉sj屬性表。打開sj層屬性對話框，對sj層進(jìn)行標(biāo)注：鼠標(biāo)在sj層上【右鍵】【屬性】【標(biāo)注】：先勾選【標(biāo)注這個圖層中的要素】，然后點擊【表達(dá)式】在【表達(dá)式】中輸入：mu1 &畝，其中&是連接符號，將數(shù)字與字符“畝”連接起來。點擊【確定】，在圖層屬性對話框也點擊【確定】！不同的表達(dá)式可以取得不同的標(biāo)注效果，比如說現(xiàn)狀圖中的（如），可以通過輸入：tbbh & vbnewline & & vbnewline & dldm 來完成。標(biāo)注結(jié)果如上圖，此時要將圖層【縮放】一下，以調(diào)整標(biāo)注文字的大小，一般以標(biāo)注文字填滿最小面積的

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