任俊儒

（華北地質(zhì)勘查局五一四地質(zhì)大隊(duì)，河北 承德 067000）

隨著市場(chǎng)環(huán)境的進(jìn)一步完善，經(jīng)濟(jì)的高速發(fā)展，地質(zhì)勘查工作也變得豐富起來(lái)，從礦藏、石油到地?zé)?、礦泉水，地質(zhì)勘查項(xiàng)目涉及的內(nèi)容越來(lái)越廣泛。近年來(lái)提出了“城市地質(zhì)調(diào)查”的概念，將傳統(tǒng)的地質(zhì)勘查與城市的可持續(xù)發(fā)展相結(jié)合，進(jìn)行基礎(chǔ)性、綜合性、前瞻性的地質(zhì)勘查工作，對(duì)城市的環(huán)境承載力、資源保障力進(jìn)行評(píng)價(jià)、同時(shí)也為城市規(guī)劃、建設(shè)以及管理決策提供基礎(chǔ)數(shù)據(jù)支持。正是由于地質(zhì)勘查類項(xiàng)目的多樣化、細(xì)致化，內(nèi)容進(jìn)一步豐富詳盡，致使各類的地質(zhì)勘查專題地圖的編輯、修飾、輸出等工作量顯著增加。

傳統(tǒng)的制圖方式，是以單幅圖紙為作業(yè)單位，逐個(gè)進(jìn)行處理，需要較大的人力與物力來(lái)保障工作的順利開展。地質(zhì)勘查專題地圖的制作主要分為要素編輯、整飾、輸出、打印等步驟。本文結(jié)合實(shí)際工作項(xiàng)目，對(duì)各步驟進(jìn)行分析，針對(duì)制圖的輸出方式進(jìn)行探索，并結(jié)合Arcpy等相關(guān)功能，實(shí)現(xiàn)了專題圖的批量生成。本文將相關(guān)的技術(shù)方法進(jìn)行總結(jié)、整理，為減少人工處理、提高制圖效率、實(shí)現(xiàn)專題圖批量生成等方面提供參考經(jīng)驗(yàn)。

1 制圖流程

（1）數(shù)據(jù)準(zhǔn)備。專題地圖編制的數(shù)據(jù)準(zhǔn)備主要包括柵格數(shù)據(jù)和矢量數(shù)據(jù)的準(zhǔn)備。柵格數(shù)據(jù)一般為遙感影像，作為底圖使用；矢量數(shù)據(jù)包括基本地形、地物、專題要素等。首先將需要的各類數(shù)據(jù)進(jìn)行收集，并且注明各類數(shù)據(jù)的基本情況，例如影像分辨率、坐標(biāo)系、數(shù)據(jù)生產(chǎn)時(shí)間等，以備后續(xù)的處理和編輯工作有所依據(jù)。

（2）數(shù)據(jù)編輯。數(shù)據(jù)的編輯分為多個(gè)步驟，首先是數(shù)據(jù)格式的轉(zhuǎn)換，在收集的數(shù)據(jù)中可能存在DWG、SHP、WP、WL、WT、HDF等不同格式的數(shù)據(jù)，為了制圖統(tǒng)一化，需要將格式進(jìn)行統(tǒng)一。其次是進(jìn)行坐標(biāo)轉(zhuǎn)換，在工作當(dāng)中會(huì)遇到北京54、西安80、CGCS2000等不同的坐標(biāo)系，一方面為了制圖的統(tǒng)一，另一方面遵循測(cè)繪地理信息行業(yè)的要求，需將數(shù)據(jù)統(tǒng)一轉(zhuǎn)為CGCS2000坐標(biāo)系。再次是數(shù)據(jù)的符號(hào)化，按照規(guī)范要求添加行政界線、道路網(wǎng)、居民點(diǎn)、水系、等高線、高程點(diǎn)等，另外是專題要素的符號(hào)化，例如地質(zhì)類型、水源地、鉆孔、礦山采區(qū)分布等。第四是地圖注記，主要是對(duì)各類數(shù)據(jù)符號(hào)化的補(bǔ)充說(shuō)明，例如村莊、道路、水系的名稱，高程值注記、專題內(nèi)容的注記等。

（3）圖框編輯。圖框的編輯，主要包括圖紙的大小、內(nèi)外圖框、比例尺、指北針、圖名、格網(wǎng)點(diǎn)坐標(biāo)、坐標(biāo)系、作業(yè)單位信息等。一般采用50cm×50cm的圖紙，按照分幅要求對(duì)制圖區(qū)域進(jìn)行分割。在制圖過程中圖框的要求是統(tǒng)一的，主要的問題是處于接邊處的符號(hào)、注記等內(nèi)容可能會(huì)丟失，不能在每幅圖紙上正確顯示，所以要進(jìn)行檢查、修改。

（4）地圖輸出。地圖的輸出是指地理文檔在完成數(shù)據(jù)編輯、圖框編輯之后輸出為圖紙的過程。地圖的輸出之前都是紙質(zhì)的，但隨著技術(shù)的發(fā)展，紙質(zhì)圖紙的利用率降低，現(xiàn)在的地圖輸出一般指電子圖紙，例如矢量圖、柵格圖片等。如果需要紙質(zhì)圖紙可以直接利用打印設(shè)備或者繪圖儀將電子地圖紙質(zhì)化，重要的一點(diǎn)是確保電子地圖能夠達(dá)到柵格采樣分辨率要求[1]。

2 流程分析

地質(zhì)勘查專題圖的制作一般都是采用CAD單幅圖紙編輯，在圖形繪制過程中比較方便，但是在數(shù)據(jù)處理、注記添加、圖框編輯等方面都比較繁瑣，逐個(gè)進(jìn)行編輯使工作量增大，工作內(nèi)容重復(fù)，錯(cuò)誤率提升。在地圖編輯過程中我們采用了ArcGIS軟件，提高數(shù)據(jù)處理以及圖框編輯的工作效率。

除圖紙編輯之外最為耗時(shí)耗工的是圖紙的輸出工作，結(jié)合現(xiàn)階段項(xiàng)目的具體要求，以提供柵格電子地圖為主，例如JPG格式的電子地圖或者地圖冊(cè)。利用CAD進(jìn)行電子地圖輸出，工作效率較低，而且是反復(fù)的機(jī)械操作，為了能夠快速、便捷、批量地完成專題圖的輸出工作，借助ArcGIS的腳本語(yǔ)言Python，對(duì)批量制圖技術(shù)中各工序編寫批處理腳本來(lái)實(shí)現(xiàn)[2]。

3 Arcpy簡(jiǎn)介

Arcpy是Python語(yǔ)言的站點(diǎn)包，其在ArcGIS for Desktop環(huán)境下，通過編寫腳本的方式高效完成地理數(shù)據(jù)分析、轉(zhuǎn)換、管理以及地圖自動(dòng)化等功能。

Python是一種跨平臺(tái)的最初設(shè)計(jì)為編寫自動(dòng)化腳本的計(jì)算機(jī)程序語(yǔ)言，其結(jié)合了解釋性、互動(dòng)性、編譯性和面向?qū)ο蟮忍攸c(diǎn)。該語(yǔ)言定位是“優(yōu)雅”、“明確”、“簡(jiǎn)單”，簡(jiǎn)單來(lái)講Python語(yǔ)言是通過調(diào)用不同的站點(diǎn)包來(lái)實(shí)現(xiàn)功能的，因此Python語(yǔ)言又稱為膠水語(yǔ)言。

基于Python語(yǔ)言的特點(diǎn)，ArcGIS選擇其為腳本語(yǔ)言，自9.0版本開始就集成了python語(yǔ)言包。Arcpy主要包含數(shù)據(jù)訪問、地理分析、自動(dòng)化制圖等模塊，數(shù)據(jù)訪問模塊主要是實(shí)現(xiàn)各版本數(shù)據(jù)編輯的調(diào)用，包括矢量以及屬性域等；地理分析模塊主要是實(shí)現(xiàn)空間分析和運(yùn)行、網(wǎng)絡(luò)分析工具、地統(tǒng)計(jì)分析等功能；自動(dòng)化制圖模塊主要是實(shí)現(xiàn)操作mxd、擴(kuò)展數(shù)據(jù)驅(qū)動(dòng)、自動(dòng)導(dǎo)出PDF和打印等功能。

4 批量生成專題圖實(shí)例

4.1 數(shù)據(jù)整理

數(shù)據(jù)整理是將收集到的數(shù)據(jù)經(jīng)過格式轉(zhuǎn)換、統(tǒng)一坐標(biāo)系等操作之后統(tǒng)一放置在一個(gè)數(shù)據(jù)集中。本文是將數(shù)據(jù)整理到文件地理數(shù)據(jù)庫(kù)(.gdb)中，也可以整理到個(gè)人地理數(shù)據(jù)庫(kù)(.mdb)中，將數(shù)據(jù)按照類型進(jìn)行分類，并依據(jù)屬性將數(shù)據(jù)符號(hào)化。在ArcMap中添加添加整張影像底圖，根據(jù)tfw文件對(duì)影像進(jìn)行配置。至此各類圖形要素均可顯示，然后將行政界線、交通網(wǎng)、水系等，以及專題要素按照特定的符號(hào)顯示出來(lái)。

4.2 圖框制作

圖形的輸出一般采用50cm×50cm分幅，采用的紙張為620mm的卷紙，在輸出圖件之前要設(shè)置參數(shù)，ArcMap界面下打開“頁(yè)面和打印設(shè)置”，選擇自定義紙張，長(zhǎng)寬均為62cm；布局視圖下的數(shù)據(jù)框的大小設(shè)置為55cm，位置居中；然后在圖框范圍內(nèi)插入指北針、比例尺、圖名、坐標(biāo)系、作業(yè)單位等信息。

4.3 數(shù)據(jù)驅(qū)動(dòng)頁(yè)面

數(shù)據(jù)驅(qū)動(dòng)頁(yè)面是根據(jù)特定的圖層將圖紙輸出為PDF格式的圖集。本例中地質(zhì)勘查專題圖以標(biāo)準(zhǔn)分幅格式進(jìn)行輸出，比例尺為1:10000，選擇結(jié)合圖表作為數(shù)據(jù)驅(qū)動(dòng)頁(yè)面的“索引圖層”，以該數(shù)據(jù)層的Text字段（圖幅號(hào)）作為“名稱字段”和“排序字段”。另一個(gè)參數(shù)是“范圍”，依據(jù)本例的比例尺和圖紙數(shù)據(jù)框設(shè)置，地圖范圍選“最佳大小”，邊距“大小為250m，指定使用地圖單位”。

參數(shù)設(shè)置完成之后，選擇“導(dǎo)出地圖”，設(shè)置分辨率等，由于數(shù)據(jù)驅(qū)動(dòng)頁(yè)面只有在輸出PDF格式時(shí)使用，所以選擇保存格式為PDF，“頁(yè)面”選項(xiàng)中可以確定輸出的范圍以及保存的格式是單個(gè)PDF數(shù)據(jù)集還是多個(gè)PDF文件。

4.4 批量生成JPG

圖件輸出為PDF格式并不能滿足工程項(xiàng)目的需求，最為實(shí)用的是將圖紙輸出為JPG格式，并按圖幅號(hào)命名。在后續(xù)的工作中可以通過坐標(biāo)值按照?qǐng)D幅號(hào)進(jìn)行索引查找，而且可以在沒有專業(yè)軟件的情況下打開圖件查詢相關(guān)信息。

利用Arcpy進(jìn)行腳本程序編寫即可實(shí)現(xiàn)上述相關(guān)內(nèi)容，相關(guān)代碼如下：

圖1 腳本程序編寫代碼

利用Python語(yǔ)言編寫腳本，首先“import”引入arcpy和os，利用“arcpy.mapping.MapDocument”確定工作空間，并用“CURRENT”參數(shù)確定當(dāng)前文檔為工作空間；定義一個(gè)圖層“JHTBWC”作為數(shù)據(jù)驅(qū)動(dòng)頁(yè)面圖層；使用“range”確定圖件的頁(yè)碼“pageNum”的取值范圍，將頁(yè)碼值賦值給數(shù)字驅(qū)動(dòng)頁(yè)面；使用“arcpy.UpdateCursor”依據(jù)結(jié)合圖表圖層迭代更新游標(biāo)“cursor”，通過循環(huán)遍歷每個(gè)圖件；最后使用“arcpy.mapping.ExportToJPEG”批量輸出JPG格式的圖件，通過設(shè)置該函數(shù)的參數(shù)，可以確定輸出路徑、輸出圖件的名稱、圖件的分辨率。在實(shí)際操作中，首先打開編輯好的mxd地理文檔，將Python程序代碼復(fù)制到Python窗口下，運(yùn)行即可得到相應(yīng)結(jié)果。

5 結(jié)語(yǔ)

本文依托地質(zhì)勘查專題圖的制作，通過對(duì)制圖流程的分析，發(fā)現(xiàn)在制圖過程中，尤其是電子圖輸出方面存在自動(dòng)化不夠的情況，為了能夠提高工作效率引入了Arcpy，具體的說(shuō)是利用ArcGIS的腳本語(yǔ)言Python進(jìn)行編程，使用Arcpy站點(diǎn)包中的相關(guān)函數(shù)來(lái)實(shí)現(xiàn)圖件的批量輸出。

實(shí)現(xiàn)數(shù)據(jù)批量處理主要是解決兩個(gè)關(guān)鍵問題，一是遍歷所有要處理的數(shù)據(jù)；二是將遍歷后的結(jié)果用一種數(shù)據(jù)結(jié)構(gòu)（如數(shù)組、鏈表等）存儲(chǔ)起來(lái)[3]。本文僅是結(jié)合Arcpy對(duì)圖紙的批量生成進(jìn)行了分析和探索，為快速制圖提供了實(shí)踐經(jīng)驗(yàn)。Arcpy不僅應(yīng)用于制圖，在GIS數(shù)據(jù)處理中也有廣泛的應(yīng)用。Python語(yǔ)言的應(yīng)用使地理信息工作避免了重復(fù)操作，提高了工作效率，也為GIS行業(yè)的發(fā)展提供了新的方向和切入點(diǎn)。