第6章GIS空間數(shù)據(jù)分析

1、-1-2022年5月25日13時(shí)34分地理信息系統(tǒng)原理地理信息系統(tǒng)原理 胡守成胡守成 昆明理工大學(xué)信自學(xué)院計(jì)科系昆明理工大學(xué)信自學(xué)院計(jì)科系 第六章第六章 GISGIS空間分析空間分析-2-2022年5月25日13時(shí)34分第六章第六章 GIS空間分析空間分析 空間數(shù)據(jù)分析是指一切用以提取地理空間信息以及時(shí)空分布、組合、聯(lián)系和發(fā)展的知識(shí)、涉及空間位置要素和空間關(guān)系的分析。 GIS空間分析主要包括：空間幾何分析、網(wǎng)絡(luò)分析、空間統(tǒng)計(jì)分析、數(shù)字地形分析等。-3-2022年5月25日13時(shí)34分 6.1數(shù)字高程模型數(shù)字高程模型 6.1.1 數(shù)字高程模型的概念數(shù)字高程模型的概念 1、數(shù)字地形模型、數(shù)字地形模

2、型 數(shù)字地形模型（數(shù)字地形模型（Digital Terrain Model DTM）是描述地面特征空間分布的有序數(shù)值陣列，其坐標(biāo)是描述地面特征空間分布的有序數(shù)值陣列，其坐標(biāo)可以是可以是X、Y或經(jīng)、緯度，是對(duì)區(qū)域地理空間數(shù)據(jù)描或經(jīng)、緯度，是對(duì)區(qū)域地理空間數(shù)據(jù)描述的基本形式和手段之一，是進(jìn)行地理空間分析的述的基本形式和手段之一，是進(jìn)行地理空間分析的基礎(chǔ)數(shù)據(jù)。基礎(chǔ)數(shù)據(jù)。 2、數(shù)字高程模型、數(shù)字高程模型 當(dāng)將數(shù)字地形模型的地面特征用于描述地面高當(dāng)將數(shù)字地形模型的地面特征用于描述地面高程時(shí)，程時(shí)，DTM杯稱為數(shù)字高程模型（杯稱為數(shù)字高程模型（Ditatal Elevation、ModelDEM）-4-2

3、022年5月25日13時(shí)34分 6.1數(shù)字高程模型數(shù)字高程模型3、DTM的數(shù)學(xué)表達(dá) 從數(shù)學(xué)的觀點(diǎn)來看，DTM是定義在區(qū)域D上的n維向量的有限序列 DTM =Vi i=1、2n 其中Vi是一個(gè)m維向量。 Vi =Vi1，Vi2，Vi3，Vim4、DEM的數(shù)學(xué)表達(dá) 數(shù)字高程模型是DIM的一個(gè)地形分量，它能比較準(zhǔn)確地表示地形的表面形態(tài)，表示區(qū)域D上地形的三維向量的有限序列： DEM =Vi=（Xi,Yi,Zi）Xi,YiD i=1,2,nXi,Yi表示平面位置，Zi表示相應(yīng)的高程。-5-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型 6.1.2 DEM數(shù)據(jù)采集方法 建立DEM首先必須

4、量測(cè)一些點(diǎn)的三維坐標(biāo)，這就是DEM數(shù)據(jù)采集或數(shù)據(jù)獲取，這些具有三維坐標(biāo)的點(diǎn)稱為數(shù)據(jù)點(diǎn)或參考點(diǎn)。 根據(jù)采集方式的不同，可分為選點(diǎn)采集、隨機(jī)采集、沿等高線采集、沿?cái)嗝娌杉?。根?jù)采集方式的不同，數(shù)據(jù)采集方法可分為： （1）直接從地面測(cè)量，如用GPS、全站儀、野外測(cè)量等。 （2）根據(jù)航空貨航天影像，通過攝影測(cè)量途徑獲取。 （3）從現(xiàn)有地圖地圖采集，利用數(shù)字化儀對(duì)已有地圖上的信息進(jìn)行數(shù)字化。-6-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型 6.1.3 DEM空間內(nèi)插方法 利用各種方法采集的DEM數(shù)據(jù)點(diǎn)數(shù)據(jù)必須經(jīng)過數(shù)據(jù)處理才能應(yīng)用，常用的數(shù)據(jù)處理包括：數(shù)據(jù)格式的轉(zhuǎn)換、坐標(biāo)系的變換、數(shù)

5、據(jù)編輯、數(shù)據(jù)分塊、數(shù)據(jù)內(nèi)插等。 內(nèi)插時(shí)DEM的核心問題，它貫穿在DEM的生產(chǎn)、質(zhì)量控制、精度評(píng)定和分析應(yīng)用等各個(gè)環(huán)節(jié)。DEM內(nèi)插就是根據(jù)參考點(diǎn)的高程求出其它待定點(diǎn)的高程。 DEM內(nèi)插分為整體內(nèi)插、分塊內(nèi)插和逐點(diǎn)內(nèi)插3類。常見的內(nèi)插方法有距離加權(quán)、移動(dòng)擬合、雙線性多項(xiàng)市內(nèi)插、最小二乘配置方法等。-7-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型 6.1.4 地學(xué)分析與應(yīng)用 DEM的應(yīng)用十分廣泛，在測(cè)繪中，它可以用來繪制等高線、坡度圖、坡向圖、立體透視圖，制作正射影像圖、立體景觀圖，進(jìn)行地圖的修側(cè)以及進(jìn)行三維仿真等。 1、地形因子的自動(dòng)提取 1）坡度和坡向分析 地面上某點(diǎn)的坡度是

6、表示地表面在該點(diǎn)的傾斜程度的一個(gè)量。因此是一個(gè)既有大小又有方向的矢量。 坡度和坡向計(jì)算分為兩種，一種為微分表示，另一種為矩陣3X3表示。坡度和坡向的計(jì)算通常在3X3個(gè)DEM矩形格網(wǎng)中進(jìn)行。 坡度 slope =tanP = (z/x)2+ (z/y)21/2 坡向 Dir = (-z/y)/(z/x)-8-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型 2）面積和體積的計(jì)算 （1）表面積的計(jì)算：如果是規(guī)則矩形格網(wǎng)的DEM，則將每個(gè)格網(wǎng)分解為三角形，使用海倫公式計(jì)算三角型的表面積。 S= P(P-D1)(P-D2)(P-D3) P =1/2(D1+D2+D3) Di = Xi2+

7、 Yi2+ Zi2 (1i 3) 式中：Di表示第i對(duì)三角型兩頂點(diǎn)之間的表面距離；S表示三角型的表面積；P表示三角型的周長的一半，整個(gè)DEM的表面積則是每個(gè)三角型表面積之和。 （2）投影面積的計(jì)算：投影面積是指任意多邊形在水平面上的投影面積，可以直接采用海倫公式進(jìn)行計(jì)算，只要將公式中的距離改為平面上兩點(diǎn)的距離即可。 -9-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型（3）體積的計(jì)算：DEM體積可由四棱柱和三棱柱的體積今昔累加得到。四棱柱上表面可用拋物雙面擬合，三菱住上表面可用斜平面擬合，下表面均為水平面或參考平面，計(jì)算公式如下： V3 =(Z1+Z2+Z3)/3S3 V4 =

8、(Z1+Z2+Z3+Z4)/4S4 式中：S3與S4分別是三棱柱和四棱柱的底面積。 3）地表粗糙度的計(jì)算 地表粗糙度一般定義為地表單元的曲面面積與其在水平面上的投影面積之比。但根據(jù)此定義，對(duì)表面光滑而傾角不同的斜面所求出的粗糙度卻不準(zhǔn)確。實(shí)際應(yīng)用中，常以格網(wǎng)頂點(diǎn)空間對(duì)角線L1和L2的中點(diǎn)的高差D來表示地表的粗糙度。其計(jì)算公式如下： D = |(Zi+1,j+1+Zij)/2 (Zi,j+1+Zi+1,j)/2| = 1/2|Zi+1,j+1+Zij Zi,j+1-Zi+1,j|-10-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型 2、地形剖面的繪制 剖面是一個(gè)垂直于海拔零平面

9、的假想平面，它與地形表面相交并延伸其地表與海拔零平面之間的部分。 例子: 已知兩點(diǎn)A和B，求經(jīng)過這兩點(diǎn)的剖面圖的原理是： 首先內(nèi)插出A，B兩點(diǎn)的高程值，然后求出AB連線與DEM格網(wǎng)的所有交點(diǎn)，插值出各交點(diǎn)的坐標(biāo)和高程，并將交點(diǎn)以離起始點(diǎn)的距離為序進(jìn)行排列，最后選擇一定的垂直比例尺和水平比例尺，以各點(diǎn)的高程和距始點(diǎn)的距離為縱坐標(biāo)繪制剖面圖。-11-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型3、地形特征點(diǎn)的提取 見p199頁4、立體透視圖 繪制立體透視圖是DEM的一個(gè)極其重要的應(yīng)用,立體圖是表現(xiàn)物體三維模型最直觀形象的圖形。-12-2022年5月25日13時(shí)34分6.1數(shù)字高程

10、模型數(shù)字高程模型-13-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型-14-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型-15-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型-16-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型5、可視性分析 可視性分析屬于對(duì)地形進(jìn)行最優(yōu)化處理的范疇，比如道路選擇、航海導(dǎo)航等。 可視性分析的基本因子有兩個(gè): 一個(gè)是兩點(diǎn)之間的通視性，另一個(gè)是可視域，即給定的觀察點(diǎn)所覆蓋的區(qū)域。 1）判斷兩點(diǎn)之間的可視性算法： （1）確定過觀察點(diǎn)和目標(biāo)點(diǎn)所在的線段與XY平面垂直的平面S （2）求出地形模型中與S相

11、交的所有邊。 （3）判斷相交的邊是否位于觀察點(diǎn)和目標(biāo)點(diǎn)所在的線段上，如果有一條邊在其上，則觀察點(diǎn)和目標(biāo)點(diǎn)不可視。 -17-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型5、可視性分析2）計(jì)算可視域的算法 （1）基于規(guī)則格網(wǎng)DEM的可視域的算法： 沿著視線的方向，從視點(diǎn)開始到目標(biāo)格網(wǎng)點(diǎn)，計(jì)算與視線相交的格網(wǎng)單元（邊或面），判斷相交的格網(wǎng)單元是否可視，從而確定視點(diǎn)與目標(biāo)點(diǎn)之間是否可視。 （2）基于TIN地形模型的可視域的算法： 一般通過計(jì)算地形中單個(gè)的三角型面元可視的部分來實(shí)現(xiàn)。 -18-2022年5月25日13時(shí)34分6.1數(shù)字高程模型數(shù)字高程模型6.1.5 DEM水文分析 基于

12、DEM地表水文分析的主要內(nèi)容是利用水文分析工具提取水流方向、匯流累計(jì)量、水流長度、河流網(wǎng)絡(luò)、河網(wǎng)分級(jí)以及流域分割等。 1、基于DEM進(jìn)行流域分析的基本原理 河網(wǎng)分級(jí)是對(duì)一個(gè)線性的河流網(wǎng)絡(luò)以數(shù)字標(biāo)識(shí)的形式劃分級(jí)別。 2、ArcGIS中的Hydrology功能介紹-19-2022年5月25日13時(shí)34分6.2 空間關(guān)系與空間量算 在GIS中存儲(chǔ)了空間分布位置信息、屬性信息和拓?fù)潢P(guān)系信息?？臻g位置、關(guān)系與度量在GIS中起著至關(guān)重要的作用。 6.2.1 空間關(guān)系的類型 廣義的空間關(guān)系有三種基本類型：拓?fù)潢P(guān)系、方向關(guān)系和度量關(guān)系。 GIS空間關(guān)系主要是指地理實(shí)體之間的位置關(guān)系，即點(diǎn)、線、面3類基本數(shù)據(jù)之

13、間的關(guān)系，它們兩兩相關(guān)，共有6種關(guān)系。-20-2022年5月25日13時(shí)34分6.2 空間關(guān)系與空間量算6.2.1 空間關(guān)系的類型 1、點(diǎn)點(diǎn)關(guān)系 點(diǎn)和點(diǎn)的關(guān)系主要有：兩點(diǎn)是否相連、兩點(diǎn)之間的距離是多少等。 2、點(diǎn)線關(guān)系 點(diǎn)和線的關(guān)系主要表現(xiàn)在點(diǎn)和線的關(guān)聯(lián)關(guān)系上，如點(diǎn)是否位于線上，點(diǎn)和線之間的距離等。 3、點(diǎn)面關(guān)系 點(diǎn)和面的關(guān)系主要表現(xiàn)在空間包含關(guān)系上，如某個(gè)村子是否位于某個(gè)縣內(nèi)或某縣有多少個(gè)村子等。-21-2022年5月25日13時(shí)34分6.2 空間關(guān)系與空間量算6.2.1 空間關(guān)系的類型 4、線線關(guān)系 線和線的關(guān)系包括兩線是否鄰接、相交等。 5、線面關(guān)系 線和面的關(guān)系表現(xiàn)線是否通過面還是和面

14、關(guān)聯(lián)或包含在面之內(nèi)。 3、面面關(guān)系 面和面之間的關(guān)系主要表現(xiàn)鄰接和包含關(guān)系。-22-2022年5月25日13時(shí)34分6.2 空間關(guān)系與空間量算6.2.2 空間量算 空間量算包括幾何量算、質(zhì)心量算、距離量算和形狀量算。 1、幾何量算 點(diǎn)狀目標(biāo)：主要指點(diǎn)的坐標(biāo)。 線狀目標(biāo)：主要指線的長度、曲率和方向。 面狀目標(biāo)：主要是面的周長和面積。 體狀目標(biāo)：主要指體積和表面積。2、質(zhì)心量算 質(zhì)心是描述地理對(duì)象空間分布的重要指標(biāo)，通常定義為一個(gè)多邊形或面的幾何中心。 -23-2022年5月25日13時(shí)34分6.2 空間關(guān)系與空間量算6.2.2 空間量算 3、距離量算 距離用于描述兩個(gè)地理實(shí)體之間的遠(yuǎn)近程度。歐式

15、距離：即幾何距離。耗費(fèi)距離：考慮到阻力的影響后的距離稱為耗費(fèi)距離。距離表面：一個(gè)地點(diǎn)到其他所有地點(diǎn)的距離，構(gòu)成一個(gè)距離表面。4、形狀量算形狀量算可分為線狀地物形狀量算和面狀地物形狀量算。線狀地物形狀量算主要是指線狀目標(biāo)的起伏程度，可用彎曲度來表示。面狀地物形狀量算主要指多邊形的形狀量算，需要考慮兩個(gè)問題：一是空間一致性問題，即有孔多邊形和破碎多邊形的處理；二是多邊形邊界特征描述問題。多邊形邊界的特征描述常用指標(biāo)：（1）多邊形長、短軸之比 ；（2）周長面積之比；（3）面積長度比；（4）用形狀系數(shù)描述。 -24-2022年5月25日13時(shí)34分6.3 空間查詢分析 空間數(shù)據(jù)的查詢屬于數(shù)據(jù)庫的范疇，

16、一般定義為從數(shù)據(jù)庫中找出所有滿足屬性約束條件和空間約束條件的地理對(duì)象。6.3.1 空間數(shù)據(jù)庫查詢語言 SQL是關(guān)系數(shù)據(jù)庫管理系統(tǒng)的標(biāo)準(zhǔn)查詢語言，適合于查詢屬性數(shù)據(jù)，而空間數(shù)據(jù)比較復(fù)雜，它由點(diǎn)、線、多邊形等混合而成，因此需要對(duì)標(biāo)準(zhǔn)SQL進(jìn)行擴(kuò)展。 SQL的擴(kuò)展 SQL的擴(kuò)展是在標(biāo)準(zhǔn)數(shù)據(jù)查詢語言上加入空間關(guān)系查詢，為此需要增加空間數(shù)據(jù)類型（如點(diǎn)、線、面等）和空間操作算子（如求長度、面積、疊加等）。在給定查詢條件時(shí)也需含有空間概念，如距離、鄰接、疊加等。 -25-2022年5月25日13時(shí)34分6.3 空間查詢分析 6.3.2 空間數(shù)據(jù)的查詢 查詢和定位空間對(duì)象，并對(duì)空間對(duì)象進(jìn)行量算是地理信息系統(tǒng)的

17、基本功能之一，是地理信息系統(tǒng)進(jìn)行高層次分析的基礎(chǔ)。 空間數(shù)據(jù)的查詢檢索就是按一定條件對(duì)空間目標(biāo)的位置和屬性信息進(jìn)行查詢，從而形成一個(gè)新的數(shù)據(jù)集。 圖形與屬性互查是最常用的查詢方法之一，主要有兩種類型： 第一類是按屬性信息的要求查詢定位空間位置； 第二類是根據(jù)對(duì)象的空間位置查詢有關(guān)屬性信息。-26-2022年5月25日13時(shí)34分6.3 空間查詢分析 6.3.2 空間數(shù)據(jù)的查詢 具體查詢方式可分為以下幾種情況： 1、幾何參數(shù)查詢 幾何參數(shù)包括點(diǎn)的位置坐標(biāo)，兩點(diǎn)之間的距離，一段線目標(biāo)的長度，一個(gè)面目標(biāo)的周長或面積等。幾何參數(shù)查詢就是通過查詢屬性庫或者進(jìn)行空間計(jì)算來實(shí)現(xiàn)。 2、空間定位查詢 空間定位

18、查詢用于實(shí)現(xiàn)圖形數(shù)據(jù)和屬性數(shù)據(jù)的雙向查詢。比如給定一個(gè)點(diǎn)或者一個(gè)幾何圖形，然后檢索該圖形范圍內(nèi)的空間對(duì)象及其屬性。-27-2022年5月25日13時(shí)34分6.3 空間查詢分析 6.3.2 空間數(shù)據(jù)的查詢 具體查詢方式可分為以下幾種情況： 3、基于屬性特征的查詢 屬性查詢是給定一個(gè)屬性表和某些屬性值，找出對(duì)應(yīng)的屬性記錄或圖形。 4、基于空間關(guān)系的查詢 空間實(shí)體存在著多種空間關(guān)系，包括：拓?fù)?、距離、順序、方位等。通過空間關(guān)系查詢和定位空間實(shí)體是地理信息系統(tǒng)不同于一般數(shù)據(jù)庫查詢的主要功能之一。 5、地址匹配查詢-28-2022年5月25日13時(shí)34分6.4 緩沖區(qū)分析 6.4.1 緩沖區(qū)分析的概念

19、緩沖區(qū)是指地理空間目標(biāo)的一種影響范圍或服務(wù)范圍。它是根據(jù)數(shù)據(jù)庫中事物的點(diǎn)、線、面實(shí)體，自動(dòng)建立的其周圍一定寬度范圍內(nèi)的緩沖區(qū)多邊形。 緩沖區(qū)分析是指根據(jù)分析對(duì)象的點(diǎn)、線、面實(shí)體，自動(dòng)建立它們周圍一定距離的帶狀區(qū)，用以識(shí)別這些實(shí)體或主體鄰近對(duì)象的輻射范圍或影響程度，以便為某項(xiàng)分析或決策提供依據(jù)。-29-2022年5月25日13時(shí)34分6.4 緩沖區(qū)分析 6.4.2 緩沖區(qū)的生成 點(diǎn)緩沖區(qū)一點(diǎn)狀地物為圓心，以緩沖區(qū)距離為半徑繪制圓；線狀地物和面狀地物緩沖區(qū)的建立，是以線狀地物或面狀地物的邊界為參考線做平行線，同時(shí)在斷點(diǎn)和拐點(diǎn)等處繪制圓弧，與平行線相切構(gòu)成。 1、角平分線法 在軸線首尾點(diǎn)處，做軸線的

20、垂線并按緩沖區(qū)半徑R截出左右邊線的起止點(diǎn)，在軸線的其他轉(zhuǎn)折點(diǎn)上，用與該線關(guān)聯(lián)的前后兩鄰邊距軸線的距離為R的兩平行線的交點(diǎn)來生成緩沖區(qū)的對(duì)應(yīng)頂點(diǎn)。 -30-2022年5月25日13時(shí)34分6.4 緩沖區(qū)分析 6.4.2 緩沖區(qū)的生成 2、凸角圓弧法 在軸線首尾點(diǎn)處，做軸線的垂線并按雙線和緩沖區(qū)半徑截出左右邊線的起止點(diǎn)，在軸線其他轉(zhuǎn)折點(diǎn)處，首先判斷該點(diǎn)的凸凹性，在凸側(cè)用圓弧彌合，在凹側(cè)則用前后兩鄰邊平行線的交點(diǎn)生成對(duì)應(yīng)頂點(diǎn)。6.4.3 緩沖區(qū)分析舉例（P213）-31-2022年5月25日13時(shí)34分6.5 疊加分析 疊加主要用于多重專題圖層的綜合分析，它是指在統(tǒng)一空間參照系條件下，每次將同一地區(qū)

21、兩個(gè)地理對(duì)象的圖層進(jìn)行疊加，以產(chǎn)生空間區(qū)域的多重屬性特征，或建立地理對(duì)象之間的空間對(duì)應(yīng)關(guān)系。分為空間合成疊加和空間統(tǒng)計(jì)疊加兩種。6.5.1 點(diǎn)與多邊形疊加 點(diǎn)與多邊形疊加就是計(jì)算多邊形對(duì)點(diǎn)的包含關(guān)系，在完成點(diǎn)與多邊形的幾何關(guān)系的計(jì)算后，還要進(jìn)行屬性信息處理。最簡單的方式是將多邊形屬性信息疊加到其包含的點(diǎn)上；也可以將點(diǎn)的屬性疊加到多邊形上。-32-2022年5月25日13時(shí)34分6.5 疊加分析 6.5.2 線與多邊形疊加 線與多邊形的疊加是比較線上坐標(biāo)與多邊形坐標(biāo)的關(guān)系，判斷線是否落在多邊形內(nèi)。通常是計(jì)算線與多邊形的交點(diǎn)，只要相交，就產(chǎn)生一個(gè)節(jié)點(diǎn)，將原線分成多條弧段，并將原線和多邊形的屬性信息一起賦給新弧段。其結(jié)果就是產(chǎn)生一個(gè)新的數(shù)據(jù)層面。6.5.3 多邊形疊加 多邊形疊加是將兩個(gè)或多個(gè)多邊形圖層進(jìn)行疊加，產(chǎn)生一個(gè)新多邊形圖層的操作，其結(jié)果是將原來多邊形要素分割成新要素，新要素綜合了原來兩層或多