地質(zhì)圖空間數(shù)據(jù)庫兩種建庫模式比較

劉福魁，左群超，曹世欣，倪振平

（1.山東省地質(zhì)調(diào)查院，山東 濟(jì)南250013；2.中國地質(zhì)調(diào)查局發(fā)展研究中心，北京100037；3.山東省地質(zhì)環(huán)境監(jiān)測總站，山東 濟(jì)南250014）

建國以來，經(jīng)過山東省廣大區(qū)域地質(zhì)調(diào)查工作者的努力，形成了一大批采用野外人工定點(diǎn)、手工記錄等“傳統(tǒng)填圖”方法取得的填圖成果，這些以紙介質(zhì)為載體的成果涵蓋了全部1∶20萬和大部分1∶5萬、1∶25萬標(biāo)準(zhǔn)圖幅。數(shù)字國土工程開展以來，山東按照中國地質(zhì)調(diào)查局的統(tǒng)一部署，逐步對以往紙介質(zhì)成果地質(zhì)圖開展了地質(zhì)圖空間數(shù)據(jù)庫建庫工作。上述3種比例尺的地質(zhì)圖建庫時(shí)，都基于矢量表達(dá)地質(zhì)和地理內(nèi)容。較早完成建庫的地質(zhì)圖幅中包括了全部1∶20萬和少量1∶5萬地質(zhì)圖，后期為了保持?jǐn)?shù)據(jù)格式的一致性，剩余1∶5萬地質(zhì)圖幅在建庫時(shí)，仍采用了和前期1∶5萬和1∶20萬地質(zhì)圖建庫相同的地理關(guān)系數(shù)據(jù)模型；而稍后才開始建庫的1∶25萬地質(zhì)圖幅，則采用了較先進(jìn)的面向?qū)ο髷?shù)據(jù)模型（又稱地理數(shù)據(jù)庫模型）。

建庫工作中應(yīng)用的地理信息系統(tǒng)（GIS），是采集、貯存、管理、分析和描述整個(gè)或部分地球表面（包括大氣層在內(nèi)）與空間和地理分布有關(guān)的數(shù)據(jù)的空間信息系統(tǒng)［1］。它主要涉及與空間幾何特征有關(guān)的空間數(shù)據(jù)、記錄空間數(shù)據(jù)特征和說明信息的屬性數(shù)據(jù)及相關(guān)拓?fù)潢P(guān)系的組織和管理，具有數(shù)據(jù)量大、應(yīng)用面廣、數(shù)據(jù)模型復(fù)雜的特點(diǎn)?？臻g數(shù)據(jù)庫作為GIS的核心，是GIS技術(shù)發(fā)展的基礎(chǔ)，而空間數(shù)據(jù)模型是空間數(shù)據(jù)庫的基礎(chǔ)和核心?？臻g數(shù)據(jù)模型是尋求一種描述地理實(shí)體的有效的數(shù)據(jù)表示方法，根據(jù)應(yīng)用要求建立實(shí)體的數(shù)據(jù)結(jié)構(gòu)和實(shí)體之間的關(guān)系以便于應(yīng)用［2］。

1 空間數(shù)據(jù)模型的發(fā)展歷程

空間數(shù)據(jù)模型的發(fā)展至今已歷三代：計(jì)算機(jī)輔助設(shè)計(jì)模型的代表，CAD數(shù)據(jù)模型；地理關(guān)系數(shù)據(jù)模型的代表，Coverage數(shù)據(jù)模型；面向?qū)ο髷?shù)據(jù)模型的代表，Geodatabase數(shù)據(jù)模型。

1.1 CAD數(shù)據(jù)模型

出現(xiàn)于20世紀(jì)60～70年代，它用點(diǎn)、線、面描述空間實(shí)體，存儲(chǔ)采取二進(jìn)制文件格式。幾何要素與相關(guān)的顏色、形狀等屬性要素存放在一起。其缺點(diǎn)是只能存儲(chǔ)少量的屬性數(shù)據(jù)，不利于空間信息的描述和 GIS分析［3－6］。

1.2 Coverage數(shù)據(jù)模型

源于ESRI公司1981年推出的ArcInfo軟件，采用分層結(jié)構(gòu)對空間數(shù)據(jù)進(jìn)行組織管理，用建立了索引的二進(jìn)制文件存儲(chǔ)點(diǎn)、線、面描述的空間實(shí)體數(shù)據(jù)。ESRI公司1996年將以二維表為基礎(chǔ)的關(guān)系數(shù)據(jù)庫管理系統(tǒng)（DBMS）引入數(shù)據(jù)模型中，屬性數(shù)據(jù)就被存儲(chǔ)在DBMS表中；空間數(shù)據(jù)和屬性數(shù)據(jù)通過公共的標(biāo)識(shí)編碼關(guān)聯(lián)，實(shí)現(xiàn)了空間數(shù)據(jù)與屬性數(shù)據(jù)的有機(jī)結(jié)合。不僅如此，它還存儲(chǔ)了矢量數(shù)據(jù)間的拓?fù)潢P(guān)系。因此Coverage數(shù)據(jù)模型提高了對地理空間信息的表達(dá)能力和數(shù)據(jù)的分析能力。

但Coverage數(shù)據(jù)模型將空間數(shù)據(jù)和屬性數(shù)據(jù)存儲(chǔ)的方式，已很難適應(yīng)于以海量數(shù)據(jù)和并發(fā)處理為特征的企業(yè)級和社會(huì)級的GIS應(yīng)用；另外拓?fù)浣Y(jié)構(gòu)不夠靈活，局部的變動(dòng)必須對全局的拓?fù)潢P(guān)系重新建立［3－6］。

1.3 Geodatabase數(shù)據(jù)模型

是ESRI公司在ArcInfo 8之后引入的一個(gè)全新的數(shù)據(jù)模型，融入了面向?qū)ο蠹夹g(shù)，是建立在DBMS之上的統(tǒng)一的、智能化的空間數(shù)據(jù)庫。面向?qū)ο蟮臄?shù)據(jù)模型與人們認(rèn)識(shí)事物的觀點(diǎn)和分類方法很接近，因此直觀且簡單實(shí)用。

在Geodatabase數(shù)據(jù)模型中，為描述由現(xiàn)實(shí)地理世界抽象出的數(shù)據(jù)，引入了要素、對象和關(guān)系的概念，其中非空間實(shí)體被抽象為對象，空間實(shí)體被抽象為要素，數(shù)據(jù)對象間的關(guān)聯(lián)被抽象為關(guān)系。相同類型的對象構(gòu)成對象類，對象類沒有空間特征，但每個(gè)對象類都有其屬性、行為和規(guī)則；相同類型的要素構(gòu)成要素類，一組具有相同空間參考的要素類組成要素?cái)?shù)據(jù)集。數(shù)據(jù)對象之間的關(guān)系規(guī)則構(gòu)成的關(guān)聯(lián)集合構(gòu)成關(guān)系類［3－6］。

2 以Coverage數(shù)據(jù)模型為基礎(chǔ)的1∶5萬地質(zhì)圖的建庫模式

2.1 1∶5萬地質(zhì)圖空間數(shù)據(jù)庫中的數(shù)據(jù)組織

以Coverage數(shù)據(jù)模型為基礎(chǔ)進(jìn)行的1∶5萬地質(zhì)圖空間數(shù)據(jù)庫建設(shè)，實(shí)行水平分幅、垂直分層的數(shù)據(jù)組織管理方式。

在1∶5萬地質(zhì)圖空間數(shù)據(jù)庫中：圖元是圖面上表示空間信息特征的基本單位，分為點(diǎn)、弧段（線）、多邊形（面）三種類型。圖層作為空間數(shù)據(jù)集合，由相同類型的圖元和描述這些圖元特征的屬性數(shù)據(jù)共同構(gòu)成。描述圖層中所有圖元特征的屬性數(shù)據(jù)的集合構(gòu)成屬性表。由于每個(gè)圖層對應(yīng)的空間實(shí)體對象特征不同，與其對應(yīng)的屬性表的數(shù)據(jù)結(jié)構(gòu)也不同；數(shù)據(jù)項(xiàng)是屬性表中最小單元，規(guī)定了明確的數(shù)據(jù)類型，用于存放描述圖元某一方面特征信息的屬性數(shù)據(jù)。

由于水平分幅，為保證多幅圖拼接后每個(gè)圖形信息及相應(yīng)屬性信息的獨(dú)立性，防止圖層名重復(fù)出現(xiàn)，“1∶5萬區(qū)域地質(zhì)圖空間數(shù)據(jù)庫（分?。┙ㄔO(shè)實(shí)施細(xì)則”對圖層的劃分和命名都作了詳細(xì)的規(guī)定，每個(gè)圖層屬性表的名稱與圖層名相同。如圖1所示。

2.2 MapGIS中1∶5萬地質(zhì)圖空間數(shù)據(jù)庫數(shù)據(jù)組織的實(shí)例

1∶5萬地質(zhì)圖空間數(shù)據(jù)庫的建庫應(yīng)用了MapGIS軟件6.7前的版本。圖2為山東省1∶5萬黃墩幅地質(zhì)圖空間數(shù)據(jù)庫，由圖層文件列表欄可知黃墩幅由30個(gè)圖層文件構(gòu)成。主窗口顯示的局部圖面內(nèi)容以3個(gè)面圖層為主，分別是表達(dá) “面狀水體”的地理圖層 “LEHDF02S”、表達(dá)地質(zhì)圖層“沉積地層單位和火山沉積地層單位”的“DEHDF01D”和 “變質(zhì)巖系地層單位”的“DEHDF01B”。 在 圖 2 所 示 范 圍 內(nèi)，“LEHDF02S”有1個(gè)圖元，“DEHDF01D”有3個(gè)圖元，“DEHDF01B”有4個(gè)圖元。

主窗口內(nèi)因被選中而呈紅色顯示的圖元屬于“DEHDF01B”圖層，“區(qū)屬性編輯”對話窗口展示了該圖層屬性表的部分?jǐn)?shù)據(jù)項(xiàng)及其中所填寫的被選中圖元的屬性。黃墩幅中“DEHDF01B”圖層有24個(gè)圖元，在以“DEHDF01B”命名的DBF格式的屬性表中相應(yīng)的有24條屬性數(shù)據(jù)記錄。屬性數(shù)據(jù)在1∶5萬地質(zhì)圖空間數(shù)據(jù)庫中起著重要作用，是查詢、更新、檢索和分析等操作的基礎(chǔ)［7］。

圖1 構(gòu)成1∶5萬地質(zhì)圖空間數(shù)據(jù)庫的圖層

圖2 MapGIS中1∶5萬地質(zhì)圖空間數(shù)據(jù)庫數(shù)據(jù)組織示意圖

表1展示的是“DEHDF01B”和“DEHDF01D”2個(gè)地質(zhì)圖層屬性表的結(jié)構(gòu)。從中可以看出，不同圖層屬性表的數(shù)據(jù)項(xiàng)個(gè)數(shù)、名稱、數(shù)據(jù)類型都不相同；既使名稱和數(shù)據(jù)類型都相同的數(shù)據(jù)項(xiàng)，其數(shù)據(jù)長度也不一定相同。

3 以Geodatabase數(shù)據(jù)模型為基礎(chǔ)的1∶25萬地質(zhì)圖的建庫模式

3.1 1∶25萬地質(zhì)圖空間數(shù)據(jù)庫中的數(shù)據(jù)組織

以Geodatabase數(shù)據(jù)模型為基礎(chǔ)進(jìn)行的1∶25萬地質(zhì)圖空間數(shù)據(jù)庫建設(shè)，把數(shù)據(jù)組織成關(guān)系型的數(shù)據(jù)對象。在“數(shù)字地質(zhì)圖空間數(shù)據(jù)庫標(biāo)準(zhǔn)”中建立了15個(gè)基本要素類，8個(gè)綜合要素類，12個(gè)對象類，5個(gè)獨(dú)立要素類?？臻g對象以要素類來存儲(chǔ)，非空間對象以表來存儲(chǔ)；關(guān)系以關(guān)系類來存儲(chǔ)，關(guān)系可以存在于空間對象（要素類中的要素）之間、非空間對象（表中的行）之間、空間對象和非空間對象之間。

地質(zhì)圖中具有空間幾何特征的地質(zhì)和地理實(shí)體，以點(diǎn)、線、面的形式由要素類數(shù)據(jù)集和綜合要素類數(shù)據(jù)集表達(dá)；其中要素類具拓?fù)潢P(guān)系，綜合要素類與要素類的關(guān)系不具拓?fù)潢P(guān)系，是覆蓋關(guān)系。對象類表達(dá)基本要素類和綜合要素類的屬性，用ACCESS表來存儲(chǔ)。獨(dú)立要素類表達(dá)地質(zhì)圖附屬的圖例、柱狀圖和責(zé)任表等圖外內(nèi)容。

表1 兩個(gè)圖層屬性表結(jié)構(gòu)對比

由表2可以看出1∶5萬地質(zhì)圖空間數(shù)據(jù)庫中參與拓?fù)涞摹癉X△△△01D”和“DX△△△05M”等8個(gè)獨(dú)立的圖層，因?yàn)閼?yīng)用面向?qū)ο蠹夹g(shù)的理念，按“1∶25萬區(qū)域地質(zhì)圖空間數(shù)據(jù)庫建設(shè)技術(shù)流程及實(shí)施細(xì)則”中的規(guī)定，被合并成“＿GeoPolygon”一個(gè)圖層。

表2 1∶25萬和1∶5萬地質(zhì)圖空間數(shù)據(jù)庫圖層關(guān)系對比

3.2 RGMapGIS中1∶25萬地質(zhì)圖空間數(shù)據(jù)庫數(shù)據(jù)組織的實(shí)例

21世紀(jì)初，中國地質(zhì)調(diào)查局主導(dǎo)開發(fā)的“數(shù)字地質(zhì)調(diào)查系統(tǒng)（RGMap）”首先在25萬區(qū)調(diào)項(xiàng)目中開展試點(diǎn)，該系統(tǒng)貫穿了從野外地質(zhì)調(diào)查到最終建成1∶25萬地質(zhì)圖空間數(shù)據(jù)庫的整個(gè)過程，為保持?jǐn)?shù)據(jù)格式的一致，對前期人工完成的紙介質(zhì)1∶25萬地質(zhì)圖建庫時(shí)，也采用了RGMap的室內(nèi)綜合整理與數(shù)據(jù)處理子系統(tǒng)（RGMapGIS）［8，9］。

應(yīng)用RGMapGIS建立了山東省1∶25萬淄博市幅地質(zhì)圖空間數(shù)據(jù)庫。在基本要素類“＿GeoPolygon”圖層中的地質(zhì)和地理空間面實(shí)體有5個(gè)需人工錄入的共同屬性數(shù)據(jù)項(xiàng)：①地質(zhì)體面實(shí)體標(biāo)識(shí)號(hào)（ID號(hào)）、②地質(zhì)體面實(shí)體類型代碼（地質(zhì)體符號(hào)）、③地質(zhì)體面實(shí)體名稱（地質(zhì)體中文名稱）、④地質(zhì)體面實(shí)體時(shí)代、⑤子類型標(biāo)識(shí)。子類型標(biāo)識(shí)用1位整型數(shù)值表示，可據(jù)此使用RGMapGIS的“自動(dòng)從地質(zhì)體面實(shí)體提取相應(yīng)對象類”功能，從“＿GeoPolygon”圖層中提取相應(yīng)對象類。如存儲(chǔ)“沉積地層單位和火山沉積地層單位”屬性的對象類數(shù)據(jù)集（＿Strata），其子類型標(biāo)識(shí)是“0”；運(yùn)行自動(dòng)提取功能后，上述5個(gè)共同屬性數(shù)據(jù)項(xiàng)中與“＿Strata”有關(guān)圖元的屬性數(shù)據(jù)會(huì)被自動(dòng)繼承下來。但數(shù)據(jù)項(xiàng)名稱會(huì)依據(jù)對象類數(shù)據(jù)集的不同而改變，如“地質(zhì)體面實(shí)體類型代碼（地質(zhì)體符號(hào)）”更名為“要素分類（地質(zhì)代碼）”。

對象類數(shù)據(jù)集（＿Strata）在繼承下“＿GeoPolygon”圖層5個(gè)共同屬性數(shù)據(jù)項(xiàng)基礎(chǔ)上，又?jǐn)U展出“巖石組合名稱、巖石組合主體顏色、巖層主要沉積構(gòu)造、生物化石帶或生物組合、地層厚度、含礦性”等數(shù)據(jù)項(xiàng)。其它對象類數(shù)據(jù)集也做了類似的數(shù)據(jù)項(xiàng)擴(kuò)展。

4 1∶5萬和1∶25萬地質(zhì)圖空間數(shù)據(jù)庫2種建庫模式比較

因?yàn)榈刭|(zhì)圖空間數(shù)據(jù)庫建設(shè)依據(jù)的數(shù)據(jù)模型不同，所以1∶5萬和1∶25萬地質(zhì)圖空間數(shù)據(jù)庫的數(shù)據(jù)組織管理方式也不同。

4.1 空間數(shù)據(jù)組織管理的不同

1∶5萬地質(zhì)圖空間數(shù)據(jù)庫，在建庫過程中實(shí)行水平分幅、垂直分層的數(shù)據(jù)組織管理方式。水平分幅與當(dāng)今成果地質(zhì)圖分幅出版、存貯的現(xiàn)狀相適應(yīng)。但由圖2、表2可以看出，垂直分層將一幅地質(zhì)圖中參與統(tǒng)一拓?fù)涞牡刭|(zhì)和地理空間面實(shí)體分成了“DX△△△01D”和“LX△△△02S”等8個(gè)獨(dú)立的圖層。

這樣的方式人為地割裂了所要表達(dá)的地質(zhì)和地理空間，不能直觀反映人們的感知，也不能表達(dá)不同類型和實(shí)體之間的空間拓?fù)潢P(guān)系和語義關(guān)系。地質(zhì)圖空間數(shù)據(jù)庫中參與拓?fù)涞娜我庖粭l弧段，一般都有左右2個(gè)相鄰的多邊形。這時(shí)變動(dòng)某個(gè)圖層中的一條弧段就會(huì)帶來問題，如果該弧段左右相鄰的多邊形不在同一個(gè)圖層中，考慮到還可能涉及不參與拓?fù)涞膱D層，則最少需要在2個(gè)圖層做2次修改，才能保持各空間實(shí)體之間拓?fù)潢P(guān)系的一致性。

1∶25萬地質(zhì)圖空間數(shù)據(jù)庫在建庫時(shí)，仍采用水平分幅，但將“DX△△△01D”等8個(gè)獨(dú)立的圖層用一個(gè)基本要素類圖層“＿GeoPolygon”表達(dá)，避免了上述問題的存在。

4.2 屬性數(shù)據(jù)組織管理的不同

圖3所顯示的屬性數(shù)據(jù)的存儲(chǔ)，上部是1∶5萬地質(zhì)圖空間數(shù)據(jù)庫采用的DBF表，下部是1∶25萬地質(zhì)圖空間數(shù)據(jù)庫采用的ACCESS表。二者除了數(shù)據(jù)存儲(chǔ)格式的不同，1∶25萬地質(zhì)圖空間數(shù)據(jù)庫的屬性數(shù)據(jù)之間在“關(guān)聯(lián)、依賴”關(guān)系之外，還多了“組合、繼承”關(guān)系。

圖3 1∶5萬和1∶25萬地質(zhì)圖空間數(shù)據(jù)庫中屬性數(shù)據(jù)存儲(chǔ)對比

圖3 上部是山東1∶5萬黃墩幅地質(zhì)圖空間數(shù)據(jù)庫“DEHDF01D”圖層中各圖元的屬性。從圖3中可以看出“DEHDF01D”圖層中各圖元都有如表1所示、由9個(gè)數(shù)據(jù)項(xiàng)所組成的屬性數(shù)據(jù)結(jié)構(gòu)，每個(gè)圖元都需要錄入這9個(gè)數(shù)據(jù)項(xiàng)中的數(shù)據(jù)。黃墩幅“DEHDF01D”圖層中有很多特征相同的圖元，它們的屬性數(shù)據(jù)也是相同的。僅在圖3有限可視范圍內(nèi)，“DEHDF01D”屬性表中就有3個(gè)“沂河組”圖元、2個(gè)“萊陽群城山后組”圖元的重復(fù)屬性記錄，形成了很大的數(shù)據(jù)冗余。導(dǎo)致修改某個(gè)圖元的屬性數(shù)據(jù)時(shí)，屬性表中所有同類圖元的相關(guān)記錄都要修改，才能保證圖元屬性的一致。

圖3下部是山東1∶25萬淄博市幅地質(zhì)圖空間數(shù)據(jù)庫中“＿Strata”對象類的屬性。由前述可知，在基本要素類“＿GeoPolygon”圖層中各圖元有5個(gè)共同的屬性數(shù)據(jù)項(xiàng)，每個(gè)圖元都需要錄入這5個(gè)數(shù)據(jù)項(xiàng)中的數(shù)據(jù)，但這些屬性數(shù)據(jù)都可以被地質(zhì)圖空間數(shù)據(jù)庫中各個(gè)對象類自動(dòng)繼承下來，不需要重新錄入。另外可以看出，“＿Strata”對象類屬性表中沒有重復(fù)的屬性記錄。也就是說，在“＿GeoPolygon”圖層中所有特征相同的圖元，在與其對應(yīng)的對象類屬性表中只用一條記錄進(jìn)行表達(dá)。

所以1∶25萬地質(zhì)圖空間數(shù)據(jù)庫采用的這種屬性數(shù)據(jù)組織管理方式，直接降低了數(shù)據(jù)冗余、減少了重復(fù)錄入工作量，也給修改屬性數(shù)據(jù)帶來很大便利。

5 結(jié)論

與以地理關(guān)系數(shù)據(jù)模型為基礎(chǔ)建成的1∶5萬地質(zhì)圖空間數(shù)據(jù)庫相比，以面向?qū)ο髷?shù)據(jù)模型為基礎(chǔ)建成的1∶25萬地質(zhì)圖空間數(shù)據(jù)庫，不僅可以將空間數(shù)據(jù)和屬性數(shù)據(jù)有效地融合在一起，并且能直觀反映各空間實(shí)體之間的拓?fù)潢P(guān)系，有效降低數(shù)據(jù)冗余、減輕建庫工作量。這種建庫模式適應(yīng)了以處理海量數(shù)據(jù)和并發(fā)操作為特點(diǎn)的GIS應(yīng)用發(fā)展趨勢，對于處理復(fù)雜的地質(zhì)和地理空間對象具有很大的優(yōu)越性。

如果面向?qū)ο髷?shù)據(jù)模型在根本意義上實(shí)現(xiàn)面向?qū)ο?，必然?huì)為GIS提供更加靈活高效的存儲(chǔ)管理機(jī)制，也將給予GIS更好的發(fā)展前景。

［1］ 龔健雅.地理信息系統(tǒng)基礎(chǔ)［M］.北京：科學(xué)出版社，2001.

［2］ 魏少春，梁煥炎.淺談地理信息系統(tǒng)中的空間數(shù)據(jù)模型［J］.地礦測繪，2006，22（3）：21－23.

［3］ 葉亞琴，左澤均，陳波.面向?qū)嶓w的空間數(shù)據(jù)庫模型［J］.地球科學(xué)：中國地質(zhì)大學(xué)學(xué)報(bào)，2006，31（5）：595－599.

［4］ 李國標(biāo)，莊雅平，王鈺華.面向?qū)ο蟮腉IS數(shù)據(jù)模型—地理數(shù)據(jù)庫［J］.測繪通報(bào)，2001，6（6）：37－39.

［5］ 郭新慧，張和生.第三代地理數(shù)據(jù)模型—Geodatabase及其實(shí)現(xiàn)［J］.太原科技大學(xué)學(xué)報(bào)，2008，29（1）：5－8.

［6］ 陳靜，張樹文.第三代地理數(shù)據(jù)模型及其實(shí)現(xiàn)［J］.三晉測繪，2003，10（2）：7－9.

［7］ 劉?？R莉，蔡青.礦產(chǎn)資源潛力評價(jià)成果圖空間數(shù)據(jù)庫中屬性數(shù)據(jù)質(zhì)量檢查探討［J］.山東國土資源，2013，29（1）：31－34.

［8］ 李超嶺，楊東來，李豐丹，等.中國數(shù)字地質(zhì)調(diào)查系統(tǒng)的基本構(gòu)架及其核心技術(shù)的實(shí)現(xiàn)［J］.地質(zhì)通報(bào)，2008，27（7）：923－944.

［9］ 李超嶺，楊東來，于慶文，等.數(shù)字地質(zhì)調(diào)查系統(tǒng)操作指南（中冊）［M］.北京：地質(zhì)出版社，2011.