基于ArcGIS和蒼穹軟件的土地利用圖的制圖綜合

馮賀珍 ，翁敏

(1.武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院，湖北 武漢 430079; 2.武漢大學(xué)地理信息系統(tǒng)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430079)

1 引 言

受計(jì)算機(jī)技術(shù)及地圖技術(shù)理論的限制，地圖的全自動(dòng)化制圖綜合還不能完全實(shí)現(xiàn)，目前還需要使用人機(jī)交互的辦法。為了提高制圖綜合的質(zhì)量和速度，國(guó)內(nèi)外工作者提出了兩種解決方法:①采用單一的軟件進(jìn)行數(shù)據(jù)的綜合和入庫(kù)，馬晶等人利用ArcGIS 進(jìn)行1∶ 10 000土地利用現(xiàn)狀圖制圖綜合的研究［1］，但效率不高;②采用自主研制的軟件解決制圖綜合的問(wèn)題，艾廷華等人基于VC + +研制了地圖綜合縮編軟件Do-Map［2］，但自主研發(fā)的軟件并沒(méi)有投入市場(chǎng)，廣泛為測(cè)繪行業(yè)人士使用。目前制圖綜合問(wèn)題的難點(diǎn)主要在于以下幾個(gè)方面:①制圖綜合的復(fù)雜性及綜合后信息的不完整;②全自動(dòng)化綜合還不能完全實(shí)現(xiàn);③自動(dòng)綜合的理論不能完全轉(zhuǎn)化為計(jì)算機(jī)可以識(shí)別和執(zhí)行的工作流程;④自動(dòng)綜合使用的軟件雜亂，數(shù)據(jù)格式不統(tǒng)一。

針對(duì)以上問(wèn)題，在現(xiàn)有理論的基礎(chǔ)上，本文采用ArcGIS 10.0 軟件和蒼穹軟件結(jié)合的方法進(jìn)行土地利用現(xiàn)狀圖的編繪和數(shù)據(jù)入庫(kù)的工作。ArcGIS 10.0 軟件可以對(duì)空間數(shù)據(jù)進(jìn)行存儲(chǔ)、編輯、更新、管理、分析等操作。蒼穹(KQ)空間數(shù)據(jù)庫(kù)具有對(duì)空間數(shù)據(jù)的顯示、查詢、修改、輸出等功能，能夠進(jìn)行破碎圖斑的融合以及成果數(shù)據(jù)的入庫(kù)［3］，兩種軟件的結(jié)合可以高效地制作土地利用圖。ArcGIS 的數(shù)據(jù)格式主要有E00 、Shape、Geodatabase、Coverage，其中E00 是ArcGIS 的外部交換格式。ArcGIS 通過(guò)E00 數(shù)據(jù)格式和其他軟件相互轉(zhuǎn)化，在轉(zhuǎn)化的過(guò)程中一般存在屬性和圖形數(shù)據(jù)的丟失問(wèn)題，但是ArcGIS 的shape 數(shù)據(jù)格式能夠和蒼穹空間數(shù)據(jù)庫(kù)無(wú)縫連接，在數(shù)據(jù)轉(zhuǎn)換的過(guò)程中可以直接把ArcGIS 導(dǎo)出的shape 數(shù)據(jù)導(dǎo)入到蒼穹軟件已建好的數(shù)據(jù)庫(kù)中［4，5］。

2 地類圖斑制圖綜合的理論

2.1 土地分類規(guī)范

根據(jù)《第二次全國(guó)土地調(diào)查技術(shù)規(guī)程》的規(guī)定，土地利用現(xiàn)狀采用二級(jí)分類體系，分為一級(jí)類和二級(jí)類，其中一級(jí)類12 個(gè)，二級(jí)類57 個(gè)，二級(jí)類是在一級(jí)類的基礎(chǔ)上劃分的。12 個(gè)一級(jí)類依次為:耕地、園地、林地、草地、商服用地、工礦倉(cāng)儲(chǔ)用地、住宅用地、公共管理與公共服務(wù)用地、特殊用地、交通運(yùn)輸用地、水域及水利設(shè)施用地和其他土地［6］。在進(jìn)行與土地利用規(guī)劃有關(guān)的工作時(shí)，必須計(jì)算出地類圖斑的面積，土地利用現(xiàn)狀分類是準(zhǔn)確統(tǒng)計(jì)各地類圖斑橢球面面積和凈面積的基礎(chǔ)。第二次全國(guó)土地調(diào)查給出各地類面積之間的關(guān)系為:

A=T－S－L－K

式中:A 為地類圖斑的面積;T 為圖斑的面積;S 為零星地物的面積;L 為線狀地物的面積;K 為扣除地類的面積［7］。

2.2 制圖綜合的原則

空間數(shù)據(jù)綜合的問(wèn)題雖然復(fù)雜，但歸納起來(lái)只有3個(gè):是否需要綜合、如何綜合、綜合的實(shí)施［2］。其中是否需要綜合的問(wèn)題會(huì)因工作性質(zhì)的不同、空間尺度的不同和相關(guān)操作人員對(duì)問(wèn)題理解的不同而發(fā)生相應(yīng)的變化;如何綜合的問(wèn)題會(huì)隨著地圖學(xué)理論及相關(guān)理論的發(fā)展不斷完善;綜合的實(shí)施受計(jì)算機(jī)技術(shù)的影響很大。制圖概括的目的是將地表最重要、最基本的事物或現(xiàn)象反映在地圖上，舍去某些非本質(zhì)的、次要的事物或現(xiàn)象［8］，并可以延長(zhǎng)地圖的時(shí)效性，避免地圖很快地失去作用。總體來(lái)說(shuō)，制圖綜合的實(shí)質(zhì)是對(duì)地圖要素進(jìn)行合理的選取和概括。選取是從制圖要素中選擇出重要的物體表現(xiàn)在地圖上;概括是對(duì)物體的形狀、數(shù)量和質(zhì)量的化簡(jiǎn)。選取和概括既有區(qū)別又緊密相連，概括通常是在選取的基礎(chǔ)上進(jìn)行的［9，10］，具有科學(xué)性和藝術(shù)性。

2.3 土地利用圖制圖綜合的特點(diǎn)

土地利用現(xiàn)狀圖是表達(dá)土地資源的類型、面積、結(jié)構(gòu)、關(guān)系等要素的專題地圖，具有全覆蓋、無(wú)縫隙、無(wú)重疊的特點(diǎn)［7］。其內(nèi)容包括小圖斑的綜合、大圖斑的合并、狹長(zhǎng)圖斑的骨架化和線狀圖斑的處理。圖斑合并時(shí)需要考慮兩種鄰近方式:語(yǔ)義鄰近和空間鄰近。土地利用圖和普通地圖制圖綜合的不同之處表現(xiàn)在:①普通地圖的制圖綜合主要依據(jù)比例尺的變化而變化，土地利用圖的制圖綜合主要受地圖內(nèi)容和表示方法的限制。②普通地圖要素的制圖綜合程度一致，土地利用圖的專題要素主要包括地類界線、地類圖斑、線狀地物及其注記，其綜合的程度也不一致。③普通地圖綜合后不需要過(guò)多考慮面積的平衡，土地利用圖綜合后要保證要素之間無(wú)重疊、無(wú)縫隙，還要使綜合后各地類面積保持平衡。

3 地類圖斑自動(dòng)綜合的技術(shù)方法及試驗(yàn)分析

土地利用現(xiàn)狀圖的編制以及數(shù)據(jù)的入庫(kù)是一項(xiàng)規(guī)范性、系統(tǒng)性很強(qiáng)的工作，能夠?yàn)閲?guó)土資源等相關(guān)部門提供決策、規(guī)劃。為了避免傳統(tǒng)純手工綜合和利用單一軟件綜合造成的工作效率不高、工作量大的問(wèn)題，本次工作在青海省同仁縣已有土地利用現(xiàn)狀數(shù)據(jù)的基礎(chǔ)上，借助ArcGIS 10.0 軟件、蒼穹(KQ)數(shù)據(jù)庫(kù)軟件以及數(shù)據(jù)質(zhì)量檢查軟件，采用人機(jī)交互的方法制作土地利用現(xiàn)狀圖，最終得到數(shù)據(jù)庫(kù)成果和各種圖形、圖表成果。技術(shù)路線如圖1所示。

圖1 技術(shù)流程圖

3.1 線狀地類圖斑的制圖綜合

線狀圖斑是指河流、公路、鐵路等呈線狀的地類圖斑［11］。《第二次全國(guó)土地調(diào)查技術(shù)規(guī)程》規(guī)定了基本比例尺下圖斑上圖的最小寬度(城鎮(zhèn)、村及工礦用地為圖上4．0 mm2，耕地、園地為6．0 mm2，林地、草地等其他地類為15．0 mm2)［6］，如果不滿足上圖最小尺寸的要求就需要進(jìn)行制圖綜合。對(duì)于線狀圖斑的綜合，我們參考Bader 和Weibel 提出的骨架線剖分方法［12］，提取線狀圖斑的骨架線，并把線狀圖斑的面積分配到相鄰的圖斑中，當(dāng)幾何鄰近地塊不止一個(gè)時(shí)，需要把骨架線剖分。被骨架線分割的面狀地類圖斑即使在同一個(gè)行政區(qū)域內(nèi)，地類編碼相同也不能合并。相反，如果河流、公路、鐵路等線狀地物的寬度能夠滿足要求，就按照面狀地類圖斑對(duì)待，不需要提取骨架線。如下圖所示:圖2中的河流是合并前的線狀地類圖斑，地塊A 和B 是被河流隔開(kāi)的內(nèi)陸灘涂。圖3是綜合后的線狀河流，可以看出合并前的河流圖斑的面積被分配到與之鄰近的圖斑中，內(nèi)陸灘涂A 和B 由于被河流隔開(kāi)不能合并成一塊圖斑。

圖2 綜合前的河流

圖3 綜合后的河流

3.2 面狀地類圖斑的制圖綜合

除了河流、交通要素等線狀圖斑，其他都是面狀地類圖斑。面狀地類圖斑處理的流程如下:

(1)面狀圖斑的打散處理。由于不同類型的地類圖斑可能連在一起，首先在ArcGIS 10.0 的高級(jí)編輯器下對(duì)圖斑進(jìn)行“拆分多部分要素處理”，打散相連接的地類圖斑。

(2)面狀圖斑的制圖綜合。首先運(yùn)用蒼穹軟件根據(jù)距離最近的原則進(jìn)行微小圖斑的合并，把微小圖斑合并到與之空間鄰近的地類圖斑中，這種做法比在ArcGIS 10.0 軟件下進(jìn)行微小圖斑的合并要快速、簡(jiǎn)單，而且能夠避免手動(dòng)合并過(guò)程中容易出現(xiàn)的跨行政區(qū)劃合并、遺漏和重復(fù)合并等錯(cuò)誤。在ArcGIS 10.0軟件中，微小圖斑的合并是手動(dòng)進(jìn)行的，而蒼穹軟件下的合并是按照要素的地類代碼和權(quán)屬單位代碼根據(jù)距離最近的原則自動(dòng)合并的。如下圖所示:圖4為合并前的圖斑，圖斑個(gè)數(shù)為236 個(gè)，圖5為合并后的圖斑，圖斑個(gè)數(shù)為70 個(gè)，圖斑個(gè)數(shù)減少了166 個(gè)。

圖4 合并前的地類圖斑

圖5 合并后的地類圖斑

(3)在ArcGIS 10.0 軟件中進(jìn)行面狀圖斑的進(jìn)一步綜合。經(jīng)過(guò)蒼穹軟件處理后的地類圖斑需要在Arc-GIS 10.0 軟件中做進(jìn)一步的處理，這時(shí)要根據(jù)語(yǔ)義鄰近和空間鄰近兩方面進(jìn)行圖斑的合并。在同一行政區(qū)域內(nèi)部，語(yǔ)義鄰近要優(yōu)先選擇，如果需要合并的地類圖斑附近沒(méi)有與之地類編碼相同的地類圖斑，就應(yīng)該把該地類圖斑合并到與它距離最近的地類圖斑中。而且要保證合并后各種地類圖斑的面積變化不能過(guò)大。如圖6所示，圖斑A 和B 都是水澆地，圖斑C 是有林地，圖斑D 是天然牧草地。圖斑A 和B、C、D 空間鄰近，圖斑A 和圖斑B 不僅空間鄰近，而且語(yǔ)義鄰近，所以優(yōu)先選擇圖斑A 和圖斑B 合并，如圖7所示。如果把圖斑B 合并到圖斑A 中，合并后圖斑的屬性和圖斑A的屬性相同。相反，如果把圖斑A 合并到圖斑B 中，合并后圖斑的屬性和圖斑B 的屬性相同。另外，如果地類圖斑不屬于同一個(gè)行政區(qū)域，則不能合并。

圖6 合并前的圖斑

圖7 合并后的圖斑

3.3 數(shù)據(jù)質(zhì)量檢查及入庫(kù)

(1)把綜合后的地類圖斑先在ArcGIS 10.0 下進(jìn)行拓?fù)錂z查，檢查線線要素之間是否交叉、懸掛，線面要素之間是否懸掛，面面要素之間是否有空洞、重疊等錯(cuò)誤。

(2)圖形錯(cuò)誤檢查完畢后，把輸出的shape 數(shù)據(jù)導(dǎo)入到蒼穹軟件下已經(jīng)建好的數(shù)據(jù)庫(kù)中并對(duì)屬性數(shù)據(jù)進(jìn)行維護(hù)。主要有數(shù)據(jù)編號(hào)的維護(hù)、ID 值的維護(hù)、圖斑面積和橢球面積的維護(hù)以及圖斑面積的平差結(jié)果等。

(3)數(shù)據(jù)庫(kù)建立完成后可以在專業(yè)的土地調(diào)查質(zhì)檢軟件中建立檢查任務(wù)，對(duì)數(shù)據(jù)庫(kù)進(jìn)行質(zhì)檢。主要包括邏輯錯(cuò)誤、拓?fù)溴e(cuò)誤、報(bào)表錯(cuò)誤和圖件錯(cuò)誤等。

(4)把質(zhì)檢軟件檢查出來(lái)的錯(cuò)誤返回到ArcGIS 10.0 和蒼穹軟件下修改。最終得到無(wú)誤的數(shù)據(jù)庫(kù)和各種圖件成果。

3.4 實(shí)例應(yīng)用及結(jié)果分析

同仁縣位于青海省東南部，東鄰甘肅省夏河縣，西連貴德縣，南接澤庫(kù)縣，北與循化、尖扎縣接壤，地處青藏高原與黃土高原的過(guò)渡地帶，隆務(wù)河縱貫全境南北，形成東西部山區(qū)和中部河谷地區(qū)的地形，全縣總面積3 275 km2。研究區(qū)的土地利用數(shù)據(jù)來(lái)源于青海省第二次土地調(diào)查，原始比例尺1∶ 10 000。經(jīng)過(guò)ArcGIS 10.0 和蒼穹軟件處理后的地類圖斑面積比例基本保持平衡，微小圖斑的處理過(guò)程由手動(dòng)轉(zhuǎn)為自動(dòng)，圖斑的個(gè)數(shù)減少了2 173個(gè)，制圖綜合的結(jié)果令人滿意。最終的土地利用圖如圖8所示。

圖8 同仁縣土地利用圖

4 結(jié) 語(yǔ)

在青海省同仁縣土地利用現(xiàn)狀圖的編繪及成果輸出方面，本研究利用ArcGIS 10.0 軟件、蒼穹(KQ)數(shù)據(jù)庫(kù)軟件、土地調(diào)查質(zhì)檢軟件相結(jié)合的方法，趨利避害，順利完成了這項(xiàng)工作。與在單一軟件下處理的結(jié)果相比，大大減少了勞動(dòng)量，提高了工作效率，也提高了工作質(zhì)量。以上是基于ArcGIS 10.0 和蒼穹軟件相結(jié)合的方法對(duì)土地利用現(xiàn)狀圖制作過(guò)程的一些探討，重點(diǎn)是地類圖斑的制圖綜合。但是要實(shí)現(xiàn)地類圖斑全自動(dòng)化綜合，還需要很大的努力。

隨著計(jì)算機(jī)技術(shù)的發(fā)展，制圖綜合的自動(dòng)化水平也在不斷地提高，但是離全自動(dòng)化綜合還有一定的距離，我們要做的工作分為兩個(gè)方面:一方面是進(jìn)行科學(xué)探究，找到提高制圖綜合水平的技術(shù)和模型，服務(wù)于生產(chǎn)，比如把遺傳算法和專家知識(shí)應(yīng)用于制圖綜合是目前研究的一個(gè)熱點(diǎn);另一方面在現(xiàn)有制圖綜合理論和技術(shù)水平的基礎(chǔ)上，靈活運(yùn)用多種軟件，合理利用各種軟件的優(yōu)點(diǎn)，加快制圖綜合的速度，提高制圖綜合的質(zhì)量，減少人力物力。

［1］馬晶，畢強(qiáng)，崔利等．1:10000 土地利用現(xiàn)狀圖制圖綜合研究［J］．測(cè)繪科學(xué)，2011(5):76 ～78．

［2］艾廷華，國(guó)寶辰，黃亞峰．1:5 萬(wàn)地圖數(shù)據(jù)庫(kù)的計(jì)算機(jī)綜合縮編［J］．武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2005(4):297 ～300．

［3］劉建華．蒼穹(KQ)空間數(shù)據(jù)庫(kù)在第二次農(nóng)村調(diào)查中的應(yīng)用［J］．測(cè)繪與空間地理信息，2011(4):255 ～256．

［4］楊寧，潘敬．淺談多軟件組合運(yùn)用加快地籍詳查成圖的方法［J］．科技創(chuàng)業(yè)月刊，2010(10):179 ～181．

［5］閆琰，李燕，董秀蘭等．ArcGIS 向MapGIS 的數(shù)據(jù)格式轉(zhuǎn)換［J］．測(cè)繪與空間地理信息，2011(6):164 ～167．

［6］TD/T 1014 －2007，第二次全國(guó)土地調(diào)查技術(shù)規(guī)程［S］．

［7］張寶堯，杜震洪，劉仁義等．基于第二次全國(guó)土地調(diào)查數(shù)據(jù)的土地利用現(xiàn)狀圖制圖綜合技術(shù)研究［J］．國(guó)土資源遙感，2012(4):107 ～121．

［8］陳夢(mèng)竹，程雄，張新利等．土地利用現(xiàn)狀圖斑綜合的工程應(yīng)用方法研究［J］．大地測(cè)量與地球動(dòng)力學(xué)，2013(8):83 ～87．

［9］劉國(guó)立，薛紅琳，姚新春．淺談借助ArcGIS 軟件進(jìn)行土地利用現(xiàn)狀圖縮編的方法［J］．測(cè)繪通報(bào)，2011(9):62 ～65．

［10］鄒勝武，鄒時(shí)林，李大軍等．土地利用圖縮編中小圖斑綜合方法的改進(jìn)［J］．地理空間信息，2013(2):138 ～140．

［11］楊德生，于歡，王學(xué)杰等．土地利用現(xiàn)狀數(shù)據(jù)縮編方法探討［J］．測(cè)繪科學(xué)，2012(9):157 ～159．

［12］BADER M，WEIBEL R．Detecting and Resolving Size and Proximity Conflicts in the Generalization of Polygonal Maps［C］．The 18th ICC，Stockholm，1997(3):1525 ～1532．