GIS教案：第4章空間數(shù)據(jù)結(jié)構(gòu)

1、第4章 空間數(shù)據(jù)結(jié)構(gòu)空間數(shù)據(jù)結(jié)構(gòu)是指對(duì)空間數(shù)據(jù)邏輯模型描述的數(shù)據(jù)組織關(guān)系和編排方式，對(duì)地理信息系統(tǒng)中數(shù)據(jù)存儲(chǔ)、查詢檢索和應(yīng)用分析等操作處理的效率有著至關(guān)重要的影響。同一空間數(shù)據(jù)邏輯模型往往采用多種空間數(shù)據(jù)結(jié)構(gòu)，例如游程長(zhǎng)度編碼結(jié)構(gòu)、四叉樹結(jié)構(gòu)都是柵格數(shù)據(jù)模型的具體實(shí)現(xiàn)。空間數(shù)據(jù)結(jié)構(gòu)是地理信息系統(tǒng)溝通信息的橋梁，只有充分理解地理信息系統(tǒng)所采用的特定數(shù)據(jù)結(jié)構(gòu)，才能正確有效地使用系統(tǒng)。在地理信息系統(tǒng)中，較常用的有柵格數(shù)據(jù)結(jié)構(gòu)和矢量數(shù)據(jù)結(jié)構(gòu)，除此之外還有混合數(shù)據(jù)結(jié)構(gòu)、鑲嵌數(shù)據(jù)結(jié)構(gòu)和超圖數(shù)據(jù)結(jié)構(gòu)等。空間數(shù)據(jù)結(jié)構(gòu)的選擇取決于數(shù)據(jù)的類型、性質(zhì)和使用的方式，應(yīng)根據(jù)不同的任務(wù)目標(biāo)，選擇最有效和最合適的數(shù)據(jù)結(jié)構(gòu)。

2、 4.1矢量數(shù)據(jù)結(jié)構(gòu)矢量數(shù)據(jù)結(jié)構(gòu)對(duì)矢量數(shù)據(jù)模型進(jìn)行數(shù)據(jù)的組織。它通過(guò)記錄實(shí)體坐標(biāo)及其關(guān)系，盡可能精確地表示點(diǎn)、線、多邊形等地理實(shí)體，坐標(biāo)空間設(shè)為連續(xù)，允許任意位置、長(zhǎng)度和面積的精確定義。矢量數(shù)據(jù)結(jié)構(gòu)直接以幾何空間坐標(biāo)為基礎(chǔ)，記錄取樣點(diǎn)坐標(biāo)。按照這種數(shù)據(jù)組織方式，可以得到精美的地圖。另外，該結(jié)構(gòu)還可以對(duì)復(fù)雜數(shù)據(jù)以最小的數(shù)據(jù)冗余進(jìn)行存貯，它還具有數(shù)據(jù)精度高，存儲(chǔ)空間小等特點(diǎn)，是一種高效的圖形數(shù)據(jù)結(jié)構(gòu)。矢量數(shù)據(jù)結(jié)構(gòu)中，傳統(tǒng)的方法是幾何圖形及其關(guān)系用文件方式組織，而屬性數(shù)據(jù)通常采用關(guān)系型表文件記錄，兩者通過(guò)實(shí)體標(biāo)識(shí)符連接。由于這一特點(diǎn)使得在某些方面有便利和獨(dú)到之處，例如在計(jì)算長(zhǎng)度、面積、形狀和圖形編輯

3、、幾何變換操作中，有很高的效率和精度。矢量數(shù)據(jù)結(jié)構(gòu)按其是否明確表示地理實(shí)體間的空間關(guān)系分為實(shí)體數(shù)據(jù)結(jié)構(gòu)和拓?fù)鋽?shù)據(jù)結(jié)構(gòu)兩大類。4.1.1實(shí)體數(shù)據(jù)結(jié)構(gòu)實(shí)體數(shù)據(jù)結(jié)構(gòu)也稱spaghetti數(shù)據(jù)結(jié)構(gòu)，是指構(gòu)成多邊形邊界的各個(gè)線段，以多邊形為單元進(jìn)行組織。按照這種數(shù)據(jù)結(jié)構(gòu)，邊界坐標(biāo)數(shù)據(jù)和多邊形單元實(shí)體一一對(duì)應(yīng)，各個(gè)多邊形邊界點(diǎn)都單獨(dú)編碼并記錄坐標(biāo)。例如對(duì)圖4-1所示的多邊形A、B、C、D，可以采用兩種結(jié)構(gòu)分別組織。圖4.1 原始多邊形數(shù)據(jù)表4.1 多邊形數(shù)據(jù)文件多邊形ID坐標(biāo)類別碼A（x1，y1），（x2，y2），（x3，y3），（x4，y4），（x5，y5），（x6，y6），（x7，y7），（x8，y8

4、），（x1，y1）A102B（x1，y1），（x8，y8），（x7，y7），（x13，y13），（x12，y12），（x11，y11），（x10，y10），（x9，y9），（x1，y1）B203C（x20，y20），（x25，y25），（x24，y24），（x23，y23），（x22，y22），（x21，y21），（x20，y20）A178D（x5，y5），（x19，y19），（x18，y18），（x17，y17），（x16，y16），（x15，y15），（x14，y14），（x7，y7），（x6，y6），（x5，y5）C523點(diǎn)號(hào)坐標(biāo)1x1，y12x2，y23x3，y34x4，y425x25

5、，y25表4.2 點(diǎn)坐標(biāo)文件表4.3 多邊形文件多邊形ID點(diǎn)號(hào)串類別碼A1,2,3,4,5,6,7,8,1A102B7,8,1,9,10,11,12,13,7B203C20,21,22,23,24,25,20A178D7,13,15,16,17,18,19,5,6,7C523 第一種結(jié)構(gòu)采用表4.1組織，第二種結(jié)構(gòu)同時(shí)采用表4.2組織多邊形頂點(diǎn)坐標(biāo)，在表4.3中記錄多邊形與點(diǎn)的關(guān)系。這種數(shù)據(jù)結(jié)構(gòu)具有編碼容易、數(shù)字化操作簡(jiǎn)單和數(shù)據(jù)編排直觀等優(yōu)點(diǎn)；但這種方法有以下明顯缺點(diǎn)： 相鄰多邊形的公共邊界要數(shù)字化兩遍，造成數(shù)據(jù)冗余存儲(chǔ)，可能導(dǎo)致輸出的公共邊界出現(xiàn)間隙或重疊； 缺少多邊形的鄰域信息和圖形的拓?fù)?/p>

6、關(guān)系； 島只作為一個(gè)單個(gè)圖形，沒有建立與外界多邊形的聯(lián)系。因此，實(shí)體式數(shù)據(jù)結(jié)構(gòu)只適用于簡(jiǎn)單的系統(tǒng)，如計(jì)算機(jī)地圖制圖系統(tǒng)。4.1.2拓?fù)鋽?shù)據(jù)結(jié)構(gòu)拓?fù)潢P(guān)系是一種對(duì)空間結(jié)構(gòu)關(guān)系進(jìn)行明確定義的數(shù)學(xué)方法。具有拓?fù)潢P(guān)系的矢量數(shù)據(jù)結(jié)構(gòu)就是拓?fù)鋽?shù)據(jù)結(jié)構(gòu)。拓?fù)鋽?shù)據(jù)結(jié)構(gòu)是GIS分析和應(yīng)用功能所必需的。拓?fù)鋽?shù)據(jù)結(jié)構(gòu)沒有固定的格式，還沒有形成標(biāo)準(zhǔn)，但基本原理相同的。它們的共同的特點(diǎn)是：點(diǎn)是相互獨(dú)立的，點(diǎn)連成線，線構(gòu)成面。每條線始于起始結(jié)點(diǎn)，止于終止結(jié)點(diǎn)，并與左右多邊形相鄰接。拓?fù)鋽?shù)據(jù)結(jié)構(gòu)最重要的特征是具有拓?fù)渚庉嫻δ?；這種拓?fù)渚庉嫻δ?，不但保證數(shù)字化原始數(shù)據(jù)的自動(dòng)差錯(cuò)編輯，而且可以自動(dòng)形成封閉的多邊形邊界，為由各個(gè)單

7、獨(dú)存儲(chǔ)的弧段組成所需要的各類多邊形及建立空間數(shù)據(jù)庫(kù)奠定基礎(chǔ)。拓?fù)鋽?shù)據(jù)結(jié)構(gòu)包括索引式、雙重獨(dú)立編碼結(jié)構(gòu)、鏈狀雙重獨(dú)立編碼結(jié)構(gòu)等。1. 索引式結(jié)構(gòu)索引式數(shù)據(jù)結(jié)構(gòu)采用樹狀索引以減少數(shù)據(jù)冗余并間接增加鄰域信息，具體方法是對(duì)所有邊界點(diǎn)進(jìn)行數(shù)字化，將坐標(biāo)對(duì)以順序方式存儲(chǔ)，由點(diǎn)索引與邊界線號(hào)相聯(lián)系，以線索引與各多邊形相聯(lián)系，形成樹狀索引結(jié)構(gòu)。圖4.2和圖4.3分別為圖4.1的多邊形文件和線文件樹狀索引圖。組織這個(gè)圖需要3個(gè)表文件，第一個(gè)記錄多邊形和邊界弧段的關(guān)系，第二個(gè)記錄邊界弧段由哪些點(diǎn)組成；第三個(gè)文件記錄每個(gè)頂點(diǎn)的坐標(biāo)，具體的結(jié)構(gòu)見表4.4、表4.5和表4.6。樹狀索引結(jié)構(gòu)消除了相鄰多邊形邊界的數(shù)據(jù)冗余

8、和不一致的問(wèn)題，在簡(jiǎn)化過(guò)于復(fù)雜的邊界線或合并多邊形時(shí)可不必改造索引表，鄰域信息和島狀信息可以通過(guò)對(duì)多邊形文件的線索引處理得到（如多邊形A、B之間通過(guò)公共邊b相鄰接），但是比較繁瑣，因而給鄰域函數(shù)運(yùn)算、消除無(wú)用邊、處理島狀信息以及檢查拓?fù)潢P(guān)系等帶來(lái)一定的困難，而且兩個(gè)編碼表都要以人工方式建立，工作量大且容易出錯(cuò)。fedgdcbebaDCBA圖4.2 多邊形與線之間索引20212223242513141516171819567137191011121312345abcdefg187圖4.3 點(diǎn)與線之間的樹狀索引表4.4 點(diǎn)坐標(biāo)文件點(diǎn)ID坐標(biāo)1x1,y1表4.5 邊文件邊ID組成的點(diǎn)IDa1,2,3

9、,4,5表4.6 多邊形文件多邊形ID組成的邊IDAa,b,c2. 雙重獨(dú)立編碼結(jié)構(gòu)這種數(shù)據(jù)結(jié)構(gòu)最早是由美國(guó)人口統(tǒng)計(jì)系統(tǒng)采用的一種編碼方式，簡(jiǎn)稱DIME（Dual Independent Map Encoding）編碼系統(tǒng)，它是以城市街道為編碼主體，它的特點(diǎn)是采用了拓?fù)渚幋a結(jié)構(gòu)，這種結(jié)構(gòu)最適合于城市信息系統(tǒng)。雙重獨(dú)立編碼結(jié)構(gòu)是對(duì)圖上網(wǎng)狀或面狀要素的任何一條線段，用順序的兩點(diǎn)定義以及相鄰多邊形來(lái)予以定義。例如對(duì)圖4.4所示的多邊形數(shù)據(jù)，利用雙重獨(dú)立編碼可得到以線段為中心的拓?fù)潢P(guān)系表，如表4.7所示。圖4.4 多邊形原始數(shù)據(jù)表4.7 雙重獨(dú)立式編碼線文件結(jié)構(gòu)線號(hào)起點(diǎn)終點(diǎn)左多邊形右多邊形a16QAb

10、21QAc32QAd43BAe54CAf65CAg67QCh78QCi89QCj94BCk910QBl103QBm1113ADn1312ADo1211AD表中第一行表示線段a的方向是從結(jié)點(diǎn)1到結(jié)點(diǎn)6，其左側(cè)面域的多邊形是Q，右側(cè)面域多邊形是A。在雙重獨(dú)立式數(shù)據(jù)結(jié)構(gòu)中，結(jié)點(diǎn)與結(jié)點(diǎn)或者多邊形與多邊形之間為鄰接關(guān)系，結(jié)點(diǎn)與線段或者多邊形與線段之間為關(guān)聯(lián)關(guān)系。利用這種拓?fù)潢P(guān)系可以來(lái)組織數(shù)據(jù)，可以有效地進(jìn)行數(shù)據(jù)存儲(chǔ)正確性檢查（如多邊形是否封閉），同時(shí)便于對(duì)數(shù)據(jù)進(jìn)行更新和檢索。因?yàn)橥ㄟ^(guò)這種數(shù)據(jù)結(jié)構(gòu)的格式繪制圖形，當(dāng)多邊形的起始結(jié)點(diǎn)與終止結(jié)點(diǎn)相一致，并且按照左側(cè)面域或右側(cè)面域自動(dòng)建立一個(gè)指定的區(qū)域單元時(shí)，則

11、空間點(diǎn)的左邊應(yīng)當(dāng)自行閉合。如果不閉合，或者出現(xiàn)多余線段，則表示數(shù)據(jù)存儲(chǔ)或編碼有誤，這樣就可以達(dá)到數(shù)據(jù)自動(dòng)編輯的目的。同樣利用該結(jié)構(gòu)可以自動(dòng)形成多邊形，并可以檢查線文件數(shù)據(jù)的正確性。除線段拓?fù)潢P(guān)系文件外，雙重獨(dú)立編碼結(jié)構(gòu)還需要點(diǎn)文件和面文件，其結(jié)構(gòu)同表4.4和表4.5。DIME編碼結(jié)構(gòu)尤其適用于城市地籍宗地的管理，在宗地管理中，界址點(diǎn)對(duì)應(yīng)于點(diǎn)、界址邊對(duì)應(yīng)于線段，面對(duì)應(yīng)于多邊形，各種要素都有惟一的標(biāo)識(shí)符。3. 鏈狀雙重獨(dú)立編碼結(jié)構(gòu)鏈狀雙重獨(dú)立式數(shù)據(jù)結(jié)構(gòu)是DIME數(shù)據(jù)結(jié)構(gòu)的一種改進(jìn)。在DIME中，一條邊只能用直線兩端點(diǎn)的序號(hào)及相鄰的多邊形來(lái)表示，而在鏈狀數(shù)據(jù)結(jié)構(gòu)中，將若干直線段合為一個(gè)弧段（或鏈段）

12、，每個(gè)弧段可以有許多中間點(diǎn)。在鏈狀雙重獨(dú)立數(shù)據(jù)結(jié)構(gòu)中，主要有四個(gè)文件：多邊形文件、弧段文件、弧段點(diǎn)文件、點(diǎn)坐標(biāo)文件。多邊形文件主要由多邊形記錄組成，包括多邊形號(hào)、組成多邊形的弧段號(hào)以及周長(zhǎng)、面積、中心點(diǎn)坐標(biāo)及有關(guān)“洞”的信息等，多邊形文件也可以通過(guò)軟件自動(dòng)檢索各有關(guān)弧段生成，并同時(shí)計(jì)算出多邊形的周長(zhǎng)和面積以及中心點(diǎn)的坐標(biāo)，當(dāng)多邊形中含有“洞”時(shí)則此“洞”的面積為負(fù)，并在總面積中減去，其組成的弧段號(hào)前也冠以負(fù)號(hào)；弧段文件主要有弧記錄組成，存儲(chǔ)弧段的起止結(jié)點(diǎn)號(hào)和弧段左右多邊形號(hào)；弧段坐標(biāo)文件由一系列點(diǎn)的位置坐標(biāo)組成，一般從數(shù)字化過(guò)程獲取，數(shù)字化的順序確定了這條鏈段的方向。結(jié)點(diǎn)文件由結(jié)點(diǎn)記錄組成，存

13、儲(chǔ)每個(gè)結(jié)點(diǎn)的結(jié)點(diǎn)號(hào)、結(jié)點(diǎn)坐標(biāo)及與該結(jié)點(diǎn)連接的弧段。結(jié)點(diǎn)文件一般通過(guò)軟件自動(dòng)生成，因?yàn)樵跀?shù)字化的過(guò)程中，由于數(shù)字化操作的誤差，各弧段在同一結(jié)點(diǎn)處的坐標(biāo)不可能完全一致，需要進(jìn)行匹配處理。當(dāng)其偏差在允許范圍內(nèi)時(shí)，可取同名結(jié)點(diǎn)的坐標(biāo)平均值。如果偏差過(guò)大，則弧段需要重新數(shù)字化。對(duì)圖4.1所示的矢量數(shù)據(jù)，其鏈狀雙重獨(dú)立式數(shù)據(jù)結(jié)構(gòu)需要多邊形文件、弧段文件、弧點(diǎn)文件、點(diǎn)坐標(biāo)文件，見表4.8、表4.9、表4.10和表4.11。表4.8 多邊形文件多邊形ID弧段號(hào)屬性（如周長(zhǎng)、面積等）Aa,b,eBc,d,bCgDf,e,d,-g表4.9 弧段文件弧段ID起始點(diǎn)終結(jié)點(diǎn)左多邊形右多邊形a51QAb71ABc113Q

14、Bd137DBe75DAf135QDg2525DC表4.10 弧段點(diǎn)文件弧段ID點(diǎn)號(hào)弧段ID點(diǎn)號(hào)a5,4,3,2,1e7,6,5b7,8,1f13,14,15,16,17,18,19,5c1,9,10,11,12,13g25,20,21,22,23,24,25d13,7表4.11 點(diǎn)坐標(biāo)文件點(diǎn)號(hào)坐標(biāo)點(diǎn)號(hào)坐標(biāo)1（x1，y1）14（x14，y14）2（x2，y2）15（x15，y15）12（x12，y12）25（x25，y25）13（x13，y13）國(guó)際著名GIS軟件平臺(tái)開發(fā)商美國(guó)ESRI公司的ARCGIS產(chǎn)品中的COVERAGE數(shù)據(jù)模型就是采用鏈狀雙重獨(dú)立編碼數(shù)據(jù)結(jié)構(gòu)的。4.2 柵格數(shù)據(jù)結(jié)構(gòu)以規(guī)

15、則柵格陣列表示空間對(duì)象的數(shù)據(jù)結(jié)構(gòu)稱為柵格數(shù)據(jù)結(jié)構(gòu)。陣列中每個(gè)柵格單元上的數(shù)值表示空間對(duì)象的屬性特征。即柵格陣列中每個(gè)單元的行列號(hào)確定位置，屬性值表示空間對(duì)象的類型、等級(jí)等特征。每個(gè)柵格單元只能存在一個(gè)值。點(diǎn)線面圖4.5 實(shí)體在柵格數(shù)據(jù)結(jié)構(gòu)中的表示柵格數(shù)據(jù)結(jié)構(gòu)表示的地表是不連續(xù)的，是量化和近似離散的數(shù)據(jù)。在柵格數(shù)據(jù)結(jié)構(gòu)中，地理空間被分成相互鄰接、規(guī)則排列的柵格單元，一個(gè)柵格單元對(duì)應(yīng)于小塊地理范圍。在柵格數(shù)據(jù)結(jié)構(gòu)中，點(diǎn)用一個(gè)柵格單元表示；線狀地物則用沿線走向的一組相鄰柵格單元表示，每個(gè)柵格單元最多只有兩個(gè)相鄰單元在線上；面或區(qū)域用記有區(qū)域?qū)傩缘南噜彇鸥駟卧募媳硎荆總€(gè)柵格單元可有多于兩個(gè)的相鄰

16、單元同屬一個(gè)區(qū)域，如圖4.5所示。柵格數(shù)據(jù)結(jié)構(gòu)的顯著特點(diǎn)是：屬性明顯，定位隱含，即數(shù)據(jù)直接記錄屬性的指針或?qū)傩员旧恚谖恢脛t根據(jù)行列號(hào)轉(zhuǎn)換為相應(yīng)的坐標(biāo)給出，也就是說(shuō)定位是根據(jù)數(shù)據(jù)在數(shù)據(jù)集中的位置得到的。同時(shí)具有數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單、數(shù)學(xué)模擬方便的優(yōu)點(diǎn)；但也存在著缺點(diǎn)：數(shù)據(jù)量大、難以建立實(shí)體間的拓?fù)潢P(guān)系、通過(guò)改變分辨率而減少數(shù)據(jù)量時(shí)精度和信息量同時(shí)受損等。4.2.1柵格單元的確定1. 柵格數(shù)據(jù)的參數(shù)一個(gè)完整的柵格數(shù)據(jù)通常由以下幾個(gè)參數(shù)決定：（1）柵格形狀。柵格單元通常為矩形或正方形。特殊的情況下也可以按經(jīng)緯網(wǎng)劃分柵格單元。列行西南角格網(wǎng)坐標(biāo)（XWS，YWS）格網(wǎng)分辨率格網(wǎng)方向圖4.6 柵格數(shù)據(jù)的坐標(biāo)

17、系及描述參數(shù)柵格形狀BACO圖4.7 柵格單元屬性值的選?。?）柵格單元大小。也就是柵格單元的尺寸，即分辨率。柵格單元的合理尺寸應(yīng)能有效地逼近空間對(duì)象的分布特征，以保證空間數(shù)據(jù)的精度。但是用柵格來(lái)逼近空間實(shí)體，不論采用多細(xì)小的柵格，與原實(shí)體比都會(huì)有誤差。通常以保證最小圖斑不丟失為原則來(lái)確定合理的柵格尺寸。設(shè)研究區(qū)域某要素的最小圖斑面積為S，柵格單元的邊長(zhǎng)L用如下公式計(jì)算： （4.1）就可以保證最小的圖斑能夠得到反映。（3）柵格原點(diǎn)。柵格系統(tǒng)的起始坐標(biāo)應(yīng)將和國(guó)家基本比例尺地形圖公里網(wǎng)的交點(diǎn)相一致，或者和已有的柵格系統(tǒng)數(shù)據(jù)相一致。并同時(shí)使用公里網(wǎng)的縱橫坐標(biāo)軸作為柵格系統(tǒng)的坐標(biāo)軸。這樣在使用柵格數(shù)據(jù)

18、時(shí)，就容易和矢量數(shù)據(jù)或已有的柵格數(shù)據(jù)相配準(zhǔn)。（4）柵格的傾角。通常情況下，柵格的坐標(biāo)系統(tǒng)與國(guó)家坐標(biāo)系統(tǒng)平行。但有時(shí)候，根據(jù)應(yīng)用的需要，可以將柵格系統(tǒng)傾斜某一個(gè)角度，以方便應(yīng)用。圖4.6 給出了柵格數(shù)據(jù)結(jié)構(gòu)中柵格數(shù)據(jù)的描述參數(shù)示意圖。2. 柵格單元值的選取柵格單元取值是惟一的，但由于受到柵格大小的限制，柵格單元中可能會(huì)出現(xiàn)多個(gè)地物，那么在決定柵格單元值時(shí)應(yīng)盡量保持其真實(shí)性，對(duì)于圖4.7所示的柵格單元，要確定該單元的屬性取值，可根據(jù)需要選用如下方法。（1）中心點(diǎn)法。用位于柵格中心處的地物類型決定其取值。由于中心點(diǎn)位于代碼為C的地物范圍內(nèi)，故其取值為C。這種方法常用于有連續(xù)分布特性的地理現(xiàn)象。（2）

19、面積占優(yōu)法。以占矩形區(qū)域面積最大的地物類型作為柵格單元的代碼。從圖上看，B類地物所占面積最大，故相應(yīng)柵格單元代碼為B。（3）重要性法。根據(jù)柵格內(nèi)不同地物的重要性，選取最重要的地物類型作為相應(yīng)的柵格單元代碼。設(shè)圖中A類地物為最重要的地物類型，則柵格代碼為A。這種方法常用于有特殊意義而面積較小的地理要素，特別是點(diǎn)狀和線狀地理要素。如城鎮(zhèn)、交通線、水系等。在柵格代碼中應(yīng)盡量表示這些重要地物。（4）百分比法。根據(jù)矩形區(qū)域內(nèi)各地理要素所占面積的百分比數(shù)確定柵格單元的取值，如可記面積最大的兩類BA，也可根據(jù)B類和A類所占面積百分比數(shù)在代碼中加入數(shù)字。由于采用的取值方法不同，得到的結(jié)果也不盡相同。逼近原始精

20、度的第二種方法是縮小單個(gè)柵格單元的面積，即增加?xùn)鸥駟卧目倲?shù)，行列數(shù)也相應(yīng)地增加。這樣，每個(gè)柵格單元可代表更為精細(xì)的地面矩形單元，混合單元減少?；旌项悇e和混合的面積都大大減小，可以大大提高量算的精度；接近真實(shí)的形態(tài)，表現(xiàn)更細(xì)小的地物類型。然而增加?xùn)鸥駛€(gè)數(shù)、提高數(shù)據(jù)精度的同時(shí)也帶來(lái)了一個(gè)嚴(yán)重的問(wèn)題，那就是數(shù)據(jù)量的大幅度增加，數(shù)據(jù)冗余嚴(yán)重。4.2.2完全柵格數(shù)據(jù)結(jié)構(gòu)完全柵格數(shù)據(jù)結(jié)構(gòu)（也稱編碼）將柵格看作一個(gè)數(shù)據(jù)矩陣，逐行逐個(gè)記錄柵格單元的值。可以每行都從左到右，也可奇數(shù)行從左到右而偶數(shù)行從右到左，或者采用其他特殊的方法。這是最簡(jiǎn)單最直接的一種柵格編碼方法。通常這種編碼為柵格文件或格網(wǎng)文件。它不采用

21、任何壓縮數(shù)據(jù)的處理，因此是最直觀最基本的柵格數(shù)據(jù)組織方式。完全柵格數(shù)據(jù)的組織有三種基本方式：基于象元、基于層和基于面域，如圖4.8所示。（1）基于像元。以像元為獨(dú)立存儲(chǔ)單元，每一個(gè)象元對(duì)應(yīng)一條記錄，每條記錄中的記錄內(nèi)容包括象元坐標(biāo)及其各層屬性值的編碼；節(jié)省了許多存儲(chǔ)坐標(biāo)的空間，因?yàn)楦鲗訉?duì)應(yīng)象元的坐標(biāo)只需存儲(chǔ)一次。數(shù)據(jù)文件像元1x坐標(biāo)y坐標(biāo)層1屬性值層2屬性值層N屬性值像元2像元N數(shù)據(jù)文件層1層2層N像元1像元2像元Nx坐標(biāo)y坐標(biāo)屬性值數(shù)據(jù)文件層1層2層N多邊形1多邊形2多邊形N屬性值像元1坐標(biāo)像元2坐標(biāo)像元N坐標(biāo)（c） 基于面域方式圖4.8 柵格數(shù)據(jù)組織方式（b） 基于層方式（a）基于象元方式

22、（2）基于層。以層為存儲(chǔ)基礎(chǔ)，層中又以象元為序記錄其坐標(biāo)和對(duì)應(yīng)該層的屬性值編碼。（3）基于面域。也以層作為存儲(chǔ)基礎(chǔ)，層中再以面域?yàn)閱卧M(jìn)行記錄，記錄的內(nèi)容包括：面域編號(hào)、面域?qū)?yīng)該層的屬性值編碼、面域中所有柵格單元的坐標(biāo)；同一屬性的多個(gè)相鄰象元只需記錄一次屬性值?；谙笤臄?shù)據(jù)組織方式簡(jiǎn)單明了，便于數(shù)據(jù)擴(kuò)充和修改，但進(jìn)行屬性查詢和面域邊界提取時(shí)速度較慢；基于層的數(shù)據(jù)組織方式便于進(jìn)行屬性查詢，但每個(gè)象元的坐標(biāo)均要重復(fù)存儲(chǔ)，浪費(fèi)了存儲(chǔ)空間；基于面域的數(shù)據(jù)組織方式雖然便于面域邊界提取，但在不同層中象元的坐標(biāo)還是要多次存儲(chǔ)。4.2.4壓縮柵格數(shù)據(jù)結(jié)構(gòu)1. 游程長(zhǎng)度編碼結(jié)構(gòu)游程長(zhǎng)度（run-length

23、）編碼，也稱行程編碼，不僅是一種柵格數(shù)據(jù)無(wú)損壓縮的重要方法，也是一種柵格數(shù)據(jù)結(jié)構(gòu)。它的基本思想是：對(duì)于一幅柵格數(shù)據(jù)（或影像），常常有行（或列）方向上相鄰的若干點(diǎn)具有相同的屬性代碼，因而可采取某種方法壓縮那些重復(fù)的記錄內(nèi)容。其編碼方案是，只在各行（或列）數(shù)據(jù)值發(fā)生變化時(shí)依次記錄該值以及相同值重復(fù)的個(gè)數(shù)，從而實(shí)現(xiàn)數(shù)據(jù)的壓縮，并實(shí)現(xiàn)數(shù)據(jù)的組織。經(jīng)編碼后，原始柵格數(shù)據(jù)陣列轉(zhuǎn)換為（si，li）數(shù)據(jù)對(duì)，其中si為屬性值，li為行程。圖4.9給出了柵格數(shù)據(jù)沿行方向進(jìn)行游程長(zhǎng)度編碼的結(jié)果。0000044400044444004444880044488822448888222488882222888822228

24、888(0,5),(4,3)(0,3),(4,5)(0,2),(4,4),(8,2)(0,2),(4,3),(8,3)(2,2),(4,2),(8,4)(2,3),(8,1),(8,4)(2,4),(8,4)(2,4),(8,4)圖4.9 柵格數(shù)據(jù)的游程長(zhǎng)度編碼及其數(shù)據(jù)結(jié)構(gòu)顯然游程長(zhǎng)度編碼只用了40個(gè)整數(shù)就可以表示，而如果用前述的直接編碼卻需要64個(gè)整數(shù)表示，可見游程長(zhǎng)度編碼壓縮數(shù)據(jù)是十分有效又簡(jiǎn)便的。事實(shí)上，壓縮比的大小是與圖的復(fù)雜程度成反比的，在變化多的部分，游程數(shù)就多，變化少的部分游程數(shù)就少，原始柵格類型越簡(jiǎn)單，壓縮效率就越高。因此這種數(shù)據(jù)結(jié)構(gòu)最適宜于類型面積較大的專題要素、遙感圖像的分

25、類結(jié)構(gòu)，而不適合于類型連續(xù)變化或類型分散的分類圖。游程長(zhǎng)度編碼在柵格加密時(shí)，數(shù)據(jù)量沒有明顯增加，壓縮效率較高，且易于檢索，疊加合并等操作，運(yùn)算簡(jiǎn)單，適用于機(jī)器存貯容量小，數(shù)據(jù)需大量壓縮，而又要避免復(fù)雜的編碼解碼運(yùn)算增加處理和操作時(shí)間的情況。 2. 四叉樹數(shù)據(jù)結(jié)構(gòu)四叉樹（quadtree）數(shù)據(jù)結(jié)構(gòu)也是一種對(duì)柵格數(shù)據(jù)的壓縮編碼方法。其基本思想是將一幅柵格數(shù)據(jù)層或圖像等分為四部分，逐塊檢查其格網(wǎng)屬性值（或灰度）；如果某個(gè)子區(qū)的所有格網(wǎng)值都具有相同的值，則這個(gè)子區(qū)就不再繼續(xù)分割，否則還要把這個(gè)子區(qū)再分割成四個(gè)子區(qū)；這樣依次地分割，直到每個(gè)子塊都只含有相同的屬性值或灰度為止。圖4.10表示了對(duì)柵格數(shù)據(jù)四

26、叉樹的分割過(guò)程及其關(guān)系。這四個(gè)等分區(qū)稱為四個(gè)子象限，按順序?yàn)樽笊希∟W）、右上（NE）、左下（SW），右下（SE），其結(jié)果是一棵倒立的樹。生成四叉樹0102440022002200葉結(jié)點(diǎn)圖4.10 柵格數(shù)據(jù)的四叉樹分割0011444400114444002244440022444422224444000044440000444400004444這種從上而下的分割需要大量的運(yùn)算，因?yàn)榇罅繑?shù)據(jù)需要重復(fù)檢查才能確定劃分。當(dāng)nn的柵格單元數(shù)比較大，且區(qū)域內(nèi)容要素又比較復(fù)雜時(shí)，建立這種四叉樹的速度比較慢。另一種是采用從下而上的方法建立。對(duì)柵格數(shù)據(jù)按如下的順序進(jìn)行檢測(cè)。如果每相鄰四個(gè)網(wǎng)格值相同則進(jìn)行合并，

27、逐次往上遞歸合并，直到符合四叉樹的原則為止。這種方法重復(fù)計(jì)算較少，運(yùn)算速度較快。從圖4.10中可以看出，為了保證四叉樹能不斷的分解下去，要求柵格數(shù)據(jù)的柵格單元數(shù)必須滿足2n2n，n為極限分割次數(shù)，n+1是四叉樹的最大高度或最大層數(shù)。對(duì)于非標(biāo)準(zhǔn)尺寸的圖像需首先通過(guò)增加背景的方法將柵格數(shù)據(jù)擴(kuò)充為2n2n個(gè)單元，對(duì)不足的部分以0補(bǔ)足（在建樹時(shí)，對(duì)于補(bǔ)足部分生成的葉結(jié)點(diǎn)不存儲(chǔ)，這樣存儲(chǔ)量并不會(huì)增加。四叉樹結(jié)構(gòu)按其編碼的方法不同又分為常規(guī)四叉樹和線性四叉樹。常規(guī)四叉樹除了記錄葉結(jié)點(diǎn)之外，還要記錄中間結(jié)點(diǎn)。結(jié)點(diǎn)之間借助指針聯(lián)系，每個(gè)結(jié)點(diǎn)需要用六個(gè)量表達(dá)：四個(gè)葉結(jié)點(diǎn)指針，一個(gè)父結(jié)點(diǎn)指針和一個(gè)結(jié)點(diǎn)的屬性或灰度

28、值。這些指針不僅增加了數(shù)據(jù)貯存量，而且增加了操作的復(fù)雜性。常規(guī)四叉樹主要在數(shù)據(jù)索引和圖幅索引等方面應(yīng)用。線性四叉樹則只存貯最后葉結(jié)點(diǎn)的信息，包括葉結(jié)點(diǎn)的位置、深度和本結(jié)點(diǎn)的屬性或灰度值。所謂深度是指處于四叉樹的第幾層上，由深度可推知子區(qū)的大小。線性四叉樹葉結(jié)點(diǎn)的編號(hào)需要遵循一定的規(guī)則，這種編號(hào)稱為地址碼，它隱含了葉結(jié)點(diǎn)的位置和深度信息。最便于應(yīng)用的地址碼是十進(jìn)制Morton碼（簡(jiǎn)稱MD碼）。十進(jìn)制mortan碼可以使用柵格單元的行列號(hào)計(jì)算（遵循C語(yǔ)言規(guī)范，矩陣的第一行為“0”行、第一列為“0”列），先將十進(jìn)制的行列號(hào)轉(zhuǎn)換成二進(jìn)制數(shù)，進(jìn)行“位”運(yùn)算操作，如圖4.11 所示，即行號(hào)和列號(hào)的二進(jìn)制數(shù)

29、兩兩交叉，得到以二進(jìn)制數(shù)表示的MD碼，再將其轉(zhuǎn)換為十進(jìn)制數(shù)。圖4.11 MD碼的“位”運(yùn)算生成In In-1 In-2 . I3 I2 I1Jn Jn-1 Jn-2 . J3 J2 J1In Jn In-1 Jn-1 In-2 Jn-2 . I3 J3 I2 J2 I1 J1MD=行號(hào)I=列號(hào) J=例如下圖第二行和第三列對(duì)應(yīng)的柵格單元，其二進(jìn)制的行列號(hào)分別為：I=0010，J=0011 ；得到的MD碼為：MD=(00001101)2=(13)10 ；用類似的方法，也可以由MD碼反求柵格單元的行列號(hào)。對(duì)于88柵格單元，MD碼的排列順序如圖4.12 所示。 01234567800145161720

30、21641236718192223662891213242528297231011141526273031744323336374849525396534353839505154559864041444556576061104742434647585962631068128129132133144145148149192YX圖4.12 柵格單元的MD碼順序行方向列方向按以上MD順序，對(duì)圖4.10所示的柵格數(shù)據(jù)按線性四叉樹進(jìn)行編碼，可得到線性四叉樹數(shù)據(jù)文件，其結(jié)構(gòu)如圖4.13所示。0140222224000000MD碼屬性值004180122164322332340350362372380390

31、400440484圖4.13 按MD碼建立的線性四叉樹結(jié)構(gòu)3. 二維行程編碼結(jié)構(gòu)在生成線性四叉樹之后，仍然存在前后葉結(jié)點(diǎn)的值相同的情況，因而可以進(jìn)一步壓縮數(shù)據(jù)，將前后值相同的葉結(jié)點(diǎn)合并，形成一個(gè)新的線性表列。如圖4. 14（a）所示的線性四叉樹的線性表，是按Morton碼的大小順序排列的，可以看出，在這個(gè)表中還有屬性值相同而又相鄰排列的情況，將值相同的葉結(jié)點(diǎn)合并后的編碼表見圖4. 14（b）。這種記錄方式類似游程編碼，但是所合并的不是柵格單元，而是合并了代表范圍大小不一的葉結(jié)點(diǎn)，所以稱它為二維行程編碼。比較兩個(gè)表可以看出，二維行程編碼又進(jìn)一步壓縮了數(shù)據(jù)。MD碼屬性值00418012216432

32、2332340350362372380390400440484二維行程碼屬性值004180122164322340362380484（a）（b）圖4.14 二維行程編碼結(jié)構(gòu)及其生成二維行程編碼采用了線性四叉樹的地址碼（Morton碼），并按照碼的順序完成編碼，但卻是沒有結(jié)構(gòu)規(guī)律的四叉樹。二維行程編碼比規(guī)則的四叉樹更節(jié)省存貯空間，而且有利于以后的插入、刪除和修改等操作。它與線性四叉樹之間的相互轉(zhuǎn)換也非常容易和快速，因此可將它們視為相同的結(jié)構(gòu)概念。4.2.4鏈碼結(jié)構(gòu)32140567圖4.15 Freeman方向碼 起始點(diǎn)起始點(diǎn)65656774556701222圖4.16 線、面的鏈?zhǔn)骄幋a鏈碼數(shù)據(jù)結(jié)

33、構(gòu)首先采用弗里曼（Freeman）碼對(duì)柵格中的線或多邊形邊界進(jìn)行編碼，然后再組織為鏈碼結(jié)構(gòu)的文件。鏈?zhǔn)骄幋a將線狀地物或區(qū)域邊界表示為：由某一起始點(diǎn)和在某些基本方向上的單位矢量鏈組成。單位矢量的長(zhǎng)度為一個(gè)柵格單元，每個(gè)后續(xù)點(diǎn)可能位于其前繼點(diǎn)的8個(gè)基本方向之一（如圖4.15）。圖4.16所示的線實(shí)體和面實(shí)體可編碼為表4.12所示的方式。具體編碼過(guò)程是：起始點(diǎn)的尋找一般遵從從上到下、從左到右的原則。當(dāng)發(fā)現(xiàn)沒有記錄過(guò)的點(diǎn)，而且數(shù)值也不為零時(shí)，就是一條線或邊界線的起點(diǎn)。記下該地物的特征碼及起點(diǎn)的行列數(shù)；然后按順時(shí)針方向?qū)ほE，找到相鄰的等值點(diǎn)，并按八個(gè)方向編碼。如遇不能閉合的線段，結(jié)束后可以返回到起始點(diǎn)再

34、開始尋找下一個(gè)線段。已經(jīng)記錄過(guò)的柵格單元，可將屬性代碼置零，以免重復(fù)編碼。表4.12 鏈碼結(jié)構(gòu)文件特征碼起點(diǎn)行起點(diǎn)列鏈碼2146,5,6,5,6,7,77284,5,5,6,7,0,1,2,2,2鏈碼結(jié)構(gòu)可以有效地壓縮柵格數(shù)據(jù)，特別是對(duì)計(jì)算面積、長(zhǎng)度、轉(zhuǎn)折方向和凹凸度等運(yùn)算十分方便。缺點(diǎn)是對(duì)邊界做合并和插入等修改，編輯比較困難。這種結(jié)構(gòu)有些類似矢量結(jié)構(gòu)，但不具有區(qū)域的性質(zhì)，因此對(duì)區(qū)域空間分析運(yùn)算比較困難。圖4.17 金字塔數(shù)據(jù)結(jié)構(gòu)X0X1X2Xn4.2.5影像金字塔結(jié)構(gòu)影像金字塔是指在統(tǒng)一的空間參照下，根據(jù)用戶需要以不同分辨率進(jìn)行存儲(chǔ)與顯示，形成分辨率由粗到細(xì)、數(shù)據(jù)量由小到大的金字塔結(jié)構(gòu)。影像

35、金字塔結(jié)構(gòu)用于圖像編碼和漸進(jìn)式圖像傳輸，是一種典型的分層數(shù)據(jù)結(jié)構(gòu)形式，適合于柵格數(shù)據(jù)和影像數(shù)據(jù)的多分辨率組織，也是一種柵格數(shù)據(jù)或影像數(shù)據(jù)的有損壓縮方式。在金字塔結(jié)構(gòu)里，圖像被分層表示。在金字塔結(jié)構(gòu)的最頂層，存儲(chǔ)最低分辨率的數(shù)據(jù)；隨著金字塔層數(shù)的增加，數(shù)據(jù)的分辨率依次降低；在金字塔結(jié)構(gòu)的底層，則存儲(chǔ)能滿足用戶需要的最高分辨率的數(shù)據(jù)。每一層相當(dāng)于降低分辨率的圖像估計(jì)。影像金字塔有多種結(jié)構(gòu)，其中最簡(jiǎn)單的是M-金字塔結(jié)構(gòu),如圖4.17所示。4.3 矢柵一體化數(shù)據(jù)結(jié)構(gòu)4.3.1 柵格結(jié)構(gòu)與矢量數(shù)據(jù)結(jié)構(gòu)的比較柵格數(shù)據(jù)結(jié)構(gòu)類型具有“屬性明顯、位置隱含”的特點(diǎn)，它易于實(shí)現(xiàn)，且操作簡(jiǎn)單，有利于基于柵格的空間信息

36、模型的分析，如在給定區(qū)域內(nèi)計(jì)算多邊形面積、線密度，柵格結(jié)構(gòu)可以很快算得結(jié)果，而采用矢量數(shù)據(jù)結(jié)構(gòu)則麻煩的多；但柵格數(shù)據(jù)表達(dá)精度不高，數(shù)據(jù)存儲(chǔ)量大，工作效率較低。如要提高一倍的表達(dá)精度（柵格單元減小一半），數(shù)據(jù)量就需增加三倍，同時(shí)也增加了數(shù)據(jù)的冗余。因此，對(duì)于基于柵格數(shù)據(jù)結(jié)構(gòu)的應(yīng)用來(lái)說(shuō)，需要根據(jù)應(yīng)用項(xiàng)目的自身特點(diǎn)及其精度要求來(lái)恰當(dāng)?shù)仄胶鈻鸥駭?shù)據(jù)的表達(dá)精度和工作效率兩者之間的關(guān)系。另外，因?yàn)闁鸥駭?shù)據(jù)結(jié)構(gòu)的簡(jiǎn)單性（不經(jīng)過(guò)壓縮編碼），其數(shù)據(jù)格式容易為大多數(shù)程序設(shè)計(jì)人員和用戶所理解，基于柵格數(shù)據(jù)結(jié)構(gòu)基礎(chǔ)之上的信息共享也較矢量數(shù)據(jù)容易。最后，遙感影像本身就是以象元為單位的柵格結(jié)構(gòu)，所以，可以直接把遙感影像應(yīng)

37、用于柵格結(jié)構(gòu)的地理信息系統(tǒng)中，也就是說(shuō)柵格數(shù)據(jù)結(jié)構(gòu)比較容易和遙感相結(jié)合。矢量數(shù)據(jù)結(jié)構(gòu)類型具有“位置明顯、屬性隱含”的特點(diǎn)，它操作起來(lái)比較復(fù)雜，許多分析操作（如疊置分析等）用矢量數(shù)據(jù)結(jié)構(gòu)難于實(shí)現(xiàn)；但它的數(shù)據(jù)表達(dá)精度較高，數(shù)據(jù)存儲(chǔ)量小，輸出圖形美觀且工作效率較高。兩者的比較如表4.13所示。表4.13 柵格、矢量數(shù)據(jù)結(jié)構(gòu)的對(duì)比優(yōu)點(diǎn)缺點(diǎn)矢量數(shù)據(jù)結(jié)構(gòu)1. 數(shù)據(jù)結(jié)構(gòu)嚴(yán)密，冗余度小，數(shù)據(jù)量??；2. 空間拓?fù)潢P(guān)系清晰，易于網(wǎng)絡(luò)分析；3. 面向?qū)ο竽繕?biāo)的，不僅能表達(dá)屬性編碼，而 且能方便地記錄每個(gè)目標(biāo)的具體的屬性描述 信息；4. 能夠?qū)崿F(xiàn)圖形數(shù)據(jù)的恢復(fù)、更新和綜合；圖形顯示質(zhì)量好、精度高。1. 數(shù)據(jù)結(jié)構(gòu)處理算

38、法復(fù)雜2. 疊置分析與柵格圖組合比較難；3. 數(shù)學(xué)模擬比較困難；4. 空間分析技術(shù)上比較復(fù)雜，需要更復(fù)雜的軟、硬件條件；5. 顯示與繪圖成本比較高。柵格數(shù)據(jù)結(jié)構(gòu)1. 數(shù)據(jù)結(jié)構(gòu)簡(jiǎn)單，易于算法實(shí)現(xiàn)；2. 空間數(shù)據(jù)的疊置和組合容易，有利于與遙感數(shù)據(jù)的匹配應(yīng)用和分析；3. 各類空間分析，地理現(xiàn)象模擬均較為容易；4. 輸出方法快速建議，成本低廉。5. 圖形數(shù)據(jù)量大，用大像元減小數(shù)據(jù)量時(shí)，精度和信息量受損失；6. 難以建立空間網(wǎng)絡(luò)連接關(guān)系；7. 投影變化實(shí)現(xiàn)困難；8. 圖形數(shù)據(jù)質(zhì)量低，地圖輸出不精美。目前，大多數(shù)地理信息系統(tǒng)平臺(tái)都支持這兩種數(shù)據(jù)結(jié)構(gòu)，而在應(yīng)用過(guò)程中，應(yīng)該根據(jù)具體的目的，選用不同的數(shù)據(jù)結(jié)構(gòu)。

39、例如，在集成遙感數(shù)據(jù)以及進(jìn)行空間模擬運(yùn)算（如污染擴(kuò)散）等應(yīng)用中，一般采用柵格數(shù)據(jù)為主要數(shù)據(jù)結(jié)構(gòu)；而在網(wǎng)絡(luò)分析、規(guī)劃選址等應(yīng)用中，通常采用矢量數(shù)據(jù)結(jié)構(gòu)。4.3.2 矢柵一體化數(shù)據(jù)結(jié)構(gòu)數(shù)據(jù)結(jié)構(gòu)矢量和柵格數(shù)據(jù)結(jié)構(gòu)各有優(yōu)缺點(diǎn)，如何充分利用兩者的優(yōu)點(diǎn)，在同一個(gè)系統(tǒng)中將兩者結(jié)合起來(lái)，是GIS中的一個(gè)重要理論與技術(shù)問(wèn)題。為將矢量與柵格數(shù)據(jù)更加有效地結(jié)合與處理，龔建雅（1993）研究提出了矢柵一體化結(jié)構(gòu)。這種數(shù)據(jù)結(jié)構(gòu)中，同時(shí)具有矢量實(shí)體的概念，又具有柵格覆蓋的思想。其理論基礎(chǔ)是：多級(jí)格網(wǎng)方法、三個(gè)基本約定和線性四叉樹編碼。圖4.18 一體化數(shù)據(jù)結(jié)構(gòu)細(xì)分格網(wǎng)多級(jí)格網(wǎng)的方法是將柵格劃分成多級(jí)格網(wǎng)：粗格網(wǎng)、基本格網(wǎng)

40、和細(xì)分格網(wǎng)（圖4.18）。粗格網(wǎng)用于建立空間索引，基本格網(wǎng)的大小與通常柵格劃分的原則基本一致，即基本柵格的大小。由于基本柵格的分辨率較低，難以滿足精度要求，所以在基本格網(wǎng)的基礎(chǔ)上又細(xì)分為256256或1616個(gè)格網(wǎng)，以增加?xùn)鸥竦目臻g分辨率，從而提高點(diǎn)、線的表達(dá)精度。粗格網(wǎng)、基本格網(wǎng)和細(xì)分格網(wǎng)都采用線性四叉樹編碼的方法，用三個(gè)Morton碼（M0、M1、M2）表示，其中M0表示點(diǎn)或線所通過(guò)的粗格網(wǎng)的Morton碼，是研究區(qū)的整體編碼；M1表示點(diǎn)或線所通過(guò)的基本格網(wǎng)的Morton碼，也是研究區(qū)的整體編碼；M2表示點(diǎn)或線所通過(guò)的細(xì)分格網(wǎng)的Morton碼，是基本柵格內(nèi)的局部編碼。以上編碼是基于柵格的，

41、因而據(jù)此設(shè)計(jì)的數(shù)據(jù)結(jié)構(gòu)必定具有柵格的性質(zhì)。為了使之具有矢量的特點(diǎn)，對(duì)點(diǎn)狀地物、線狀地物、和面狀地物的作三個(gè)約定：（1）點(diǎn)狀地物僅有空間位置而無(wú)形狀和面積，在計(jì)算機(jī)中僅有一個(gè)坐標(biāo)數(shù)據(jù)；（2）線狀地物有形狀但無(wú)面積，在計(jì)算機(jī)中需要組織一組元子（即柵格單元）填滿的路徑表達(dá)；（3）面狀地物有形狀和面積，在計(jì)算機(jī)中有一組元子表達(dá)的填滿路徑的邊界線和內(nèi)部（空洞外均填滿）的區(qū)域組成。據(jù)此，點(diǎn)狀地物、線狀地物、面狀地物的數(shù)據(jù)組織方式如下：點(diǎn)狀地物：用（M1，M2）代替矢量數(shù)據(jù)結(jié)構(gòu)中的（x，y）；點(diǎn)標(biāo)識(shí)號(hào)M1M2屬性值線狀地物：用Morton碼代替（x，y）記錄原始采樣的中間點(diǎn)的位置；必要時(shí)，還可記錄線目標(biāo)所穿

42、過(guò)的基本網(wǎng)格的交線位置：弧ID起點(diǎn)ID終點(diǎn)ID左域ID右域ID中間點(diǎn)坐標(biāo)（M1，M2）序列面狀地物：除用Morton碼代替（x，y）記錄面狀地物邊界原始采樣的中間點(diǎn)位置，以及它們所穿過(guò)的所有基本網(wǎng)格的交線位置之外，還要用鏈指針記錄多邊形的內(nèi)部柵格。必要時(shí)，還可以記錄邊界所穿過(guò)的所有基本網(wǎng)格的交線位置：面域ID邊界ID序列面域內(nèi)點(diǎn)指針面域內(nèi)點(diǎn)指針面域內(nèi)點(diǎn)坐標(biāo)（M1，M2）序列因此，點(diǎn)狀地物、線狀地物和面狀地物不僅具有如同矢量數(shù)據(jù)結(jié)構(gòu)的位置“坐標(biāo)”，而且還可以有類型編碼、屬性值和拓?fù)潢P(guān)系，因而具有完全的矢量特性。與此同時(shí)，由于用柵格元子表達(dá)了點(diǎn)、填充了線性目標(biāo)、多邊形邊界及其內(nèi)部（空洞除外），實(shí)際

43、是進(jìn)行了柵格化，因而可以進(jìn)行各種柵格操作。 4.4 鑲嵌數(shù)據(jù)結(jié)構(gòu)針對(duì)以正方形和矩形單元進(jìn)行地理空間劃分的規(guī)則鑲嵌數(shù)據(jù)模型，完全可以采用柵格數(shù)據(jù)結(jié)構(gòu)進(jìn)行數(shù)據(jù)的組織。本節(jié)重點(diǎn)討論以Voronoi多邊形和TIN三角網(wǎng)進(jìn)行空間劃分的不規(guī)則鑲嵌數(shù)據(jù)模型的數(shù)據(jù)結(jié)構(gòu)。4.4.1 Voronoi數(shù)據(jù)結(jié)構(gòu)以Voronoi 面塊單元來(lái)組織Voronoi多邊形數(shù)據(jù)。如圖4.19 所示，節(jié)點(diǎn)p1p8 將平面凸多邊形ABCDE區(qū)域劃分為8個(gè)Voronoi 單元，定義相鄰Voronoi 單元之間的垂直平分線的交點(diǎn)為頂點(diǎn)，且各頂點(diǎn)分別賦予頂點(diǎn)ID號(hào)。對(duì)于每個(gè)樣點(diǎn)及其Voronoi 單元，其數(shù)據(jù)結(jié)構(gòu)包括一組文件，在這些文件中

44、分別記錄：（1） 樣點(diǎn)序號(hào)（與Voronoi 單元的ID號(hào)相同）、樣點(diǎn)的坐標(biāo)及其屬性；（2） 各樣點(diǎn)所對(duì)應(yīng)的相鄰樣點(diǎn)的序號(hào)；（3） 生成的Voronoi 單元的頂點(diǎn)坐標(biāo)；（4） Voronoi單元的頂點(diǎn)組成。則其數(shù)據(jù)結(jié)構(gòu)如表4.14 表4.15 所示。在記錄相鄰Vorinoi單元序號(hào)時(shí), 統(tǒng)一采用逆時(shí)針方向。當(dāng)某樣點(diǎn)的Voronoi 單元的邊界就是地理區(qū)域的邊界時(shí)，將Vorinoi邊與地理區(qū)域邊界的交點(diǎn)也作為頂點(diǎn)，但使用特殊標(biāo)識(shí)符予以區(qū)分。這種數(shù)據(jù)結(jié)構(gòu)對(duì)Voronoi多邊形的生成過(guò)程以及基于Voronoi的空間分析將起到重要的作用。表4.15 Voronoi單元鄰接關(guān)系表Voronoi單元ID

45、相鄰Voronoi單元號(hào)12，3，4，5，6，7，821，8，331，2，441，3，551，4，661，5，771，6，881，7，2表4.16 Voronoi頂點(diǎn)信息表Voronoi頂點(diǎn)IDVoronoi頂點(diǎn)坐標(biāo)Voronoi頂點(diǎn)標(biāo)識(shí)v1x1,y1v2x2,y2v3x3,y3v4x4,y4x5,y5AxA,yA邊界點(diǎn)BxB,yB邊界點(diǎn)X8,y8圖4.19 Voronoi多邊形樣點(diǎn)ID樣點(diǎn)坐標(biāo)樣點(diǎn)屬性值1x1,y1A12x2,y2A23x3,y3A34x4,y4A45x5,y5A56x6,y6A67x7,y7A78X8,y8A8表4.14 樣點(diǎn)數(shù)據(jù)AEDB54136278C4.4.2 TIN數(shù)據(jù)結(jié)構(gòu)在數(shù)據(jù)結(jié)構(gòu)上，TIN可以采用類似于多邊形的矢量拓?fù)浣Y(jié)構(gòu)，但不必要描述一般多邊形中的“島嶼”或“洞”的拓?fù)潢P(guān)系。可以采用多種方式來(lái)組織TIN數(shù)據(jù)模型，對(duì)圖4.20所示的TIN三角網(wǎng)，可以三角形作為基本的空間對(duì)象進(jìn)行數(shù)據(jù)組織，見圖4.21。數(shù)據(jù)組織需要兩個(gè)