第10章 ArcGIS 空間統(tǒng)計分析.doc

第10章 空間統(tǒng)計分析統(tǒng)計分析是空間分析的主要手段，貫穿于空間分析的各個主要環(huán)節(jié)?？臻g統(tǒng)計分析方法不僅僅限于常規(guī)統(tǒng)計方法，還包括利用空間位置的空間自相關分析。本章主要介紹常用統(tǒng)計量、數(shù)據(jù)特征分析（即探索性數(shù)據(jù)分析）、分級統(tǒng)計分析、空間插值和空間回歸分析五方面內容。10.1概述10.1.1基本概念空間統(tǒng)計分析可包括“空間數(shù)據(jù)的統(tǒng)計分析”及“數(shù)據(jù)的空間統(tǒng)計分析”，前者著重于空間物體和現(xiàn)象的非空間特性的統(tǒng)計分析，解決的一個中心議題就是如何以數(shù)學統(tǒng)計模型來描述和模擬空間現(xiàn)象和過程，即將地理模型轉換成數(shù)學統(tǒng)計模型，以便于定量描述和計算機處理，著重于常規(guī)的統(tǒng)計分析方法，尤其是多元統(tǒng)計分析方法對空間數(shù)據(jù)的處理，而空間數(shù)據(jù)所描述的事物的空間位置在這些分析中不起制約作用。如趨勢面擬合被廣泛應用于地理數(shù)據(jù)的趨勢分析中，但在這種分析中，僅考慮了樣本值的大小，而并不考慮這些樣本在地理空間的分布特征及其相互間的位置關系。從這個意義上講，空間數(shù)據(jù)的統(tǒng)計分析在很多方面，與一般的數(shù)據(jù)分析并無本質差別，但是對空間數(shù)據(jù)的統(tǒng)計分析結果的解釋則必然要依托于地理空間進行，在很多情況下，分析的結果以地圖方式來描述和表達的。因此，空間數(shù)據(jù)的統(tǒng)計分析盡管在分析過程中沒有考慮數(shù)據(jù)抽樣點的空間位置，但描述的仍然是空間過程，揭示的也是空間規(guī)律和空間機制。數(shù)據(jù)的空間統(tǒng)計分析則是直接從空間物體的空間位置、聯(lián)系等方面出發(fā)，研究既具有隨機性又具有結構性，或具有空間相關性和依賴性的自然現(xiàn)象。凡是與空間數(shù)據(jù)的結構性和隨機性，或空間相關性和依賴性，或空間格局與變異有關的研究，并對這些數(shù)據(jù)進行最優(yōu)無偏內插估計，或模擬這些數(shù)據(jù)的離散性、波動性，都是數(shù)據(jù)的空間統(tǒng)計分析的研究內容。數(shù)據(jù)的空間統(tǒng)計分析不是拋棄了傳統(tǒng)的統(tǒng)計學的理論和方法，它是在傳統(tǒng)的統(tǒng)計學基礎上發(fā)展起來的。數(shù)據(jù)的空間統(tǒng)計學與經典統(tǒng)計學的共同之處在于：它們都是在大量采樣的基礎上，通過對樣本屬性值的頻率分布、均值、方差等關系及其相應規(guī)則的分析，確定其空間分布格局與相關關系。數(shù)據(jù)的空間統(tǒng)計學區(qū)別于經典統(tǒng)計學的最大特點是：數(shù)據(jù)的空間統(tǒng)計學既考慮到樣本值的大小，又重視樣本空間位置及樣本間的距離?？臻g數(shù)據(jù)具有空間依賴性(空間自相關) 和空間非均質性（空間結構），扭曲了經典統(tǒng)計方法的假設條件，使得經典統(tǒng)計模型對空間數(shù)據(jù)的分析會產生虛假的解釋。經典統(tǒng)計學模型是在觀測結果相互獨立的假設基礎上建立的，但實際上地理現(xiàn)象之間大都不具有獨立性。數(shù)據(jù)的空間統(tǒng)計學研究的基礎是空間對象間的相關性和非獨立的觀測，它們與距離有關，并隨著距離的增加而變化。這些問題為經典的統(tǒng)計學所忽視，但卻成為數(shù)據(jù)的空間統(tǒng)計學的核心。10.1.2主要分析內容 空間統(tǒng)計分析與經典統(tǒng)計學的內容往往是交叉的。空間統(tǒng)計分析使用統(tǒng)計方法解釋空間數(shù)據(jù)，分析數(shù)據(jù)在統(tǒng)計上是否是“典型”的，或“期望”的。同時，它又具有自己獨有的空間自相關分析。主要分析內容包含以下幾點：1. 基本統(tǒng)計量統(tǒng)計量是數(shù)據(jù)特征的反映，也是統(tǒng)計分析的基礎。2. 探索性數(shù)據(jù)分析探索性數(shù)據(jù)分析能讓用戶更深入了解數(shù)據(jù)，認識研究對象，從而對與其數(shù)據(jù)相關的問題做出更好的決策。探索性數(shù)據(jù)分析主要包括確定統(tǒng)計數(shù)據(jù)屬性、探測數(shù)據(jù)分布、全局和局部異常值（過大值或過小值）、尋求全局的變化趨勢、研究空間自相關和理解多種數(shù)據(jù)集之間相關性。3. 分級統(tǒng)計分析分級統(tǒng)計是對數(shù)據(jù)的進一步處理分析，以便于更好得揭示數(shù)據(jù)規(guī)律或在制圖中獲得更好的效果。4. 空間插值基于探索性數(shù)據(jù)分析結果，選擇合適的數(shù)據(jù)內插模型，由已知樣點來創(chuàng)建表面，研究空間分布。5. 空間回歸研究兩個或兩個以上的變量之間統(tǒng)計關系，通過空間關系，包括考慮空間的自相關性，把屬性數(shù)據(jù)與空間位置關系結合起來，更好的解釋地理事物的空間關系。6. 空間分類基于地圖表達，采用與變量聚類分析相類似的方法來產生新的綜合性或者簡潔性專題地圖。包括多變量統(tǒng)計分析，如主成分分析、層次分析，以及空間分類統(tǒng)計分析，如系統(tǒng)聚類分析、判別分析等。10.2 基本統(tǒng)計量平均數(shù)中位數(shù)眾數(shù)分位數(shù)偏度峰度總和比率比例種類集中趨勢描述數(shù)據(jù)特征的統(tǒng)計量離散程度基本統(tǒng)計量其他統(tǒng)計量分布形狀極值極差離差平均離差離差平方和方差標準差變差系數(shù)常用的基本統(tǒng)計量主要包括：最大值、最小值、極差、均值、中值、總和、眾數(shù)、種類、離差、方差、標準差、變差系數(shù)、峰度和偏度等。這些統(tǒng)計量反映了數(shù)據(jù)集的范圍、集中情況、離散程度、空間分布等特征，對進一步的數(shù)據(jù)分析起著鋪墊作用（見圖10.1）。圖10.1 基本統(tǒng)計量10.2.1 代表數(shù)據(jù)集中趨勢的統(tǒng)計量代表數(shù)據(jù)集中趨勢的統(tǒng)計量包括平均數(shù)、中位數(shù)、眾數(shù)，它們都可以用來表示數(shù)據(jù)的分布位置和一般水平。表10.1中，列出了各統(tǒng)計量的含義以及在實際應用的作用。其中，表示數(shù)據(jù)集中的第i個變量，i=1，2，n。表10.1 代表集中趨勢的統(tǒng)計量的含義、特點及作用統(tǒng)計量含義特點及作用平均數(shù)算術平均數(shù)簡單算術平均數(shù)n個數(shù)據(jù)的總和與數(shù)據(jù)的總個數(shù)n的比值平均數(shù)是最常用的表示數(shù)據(jù)集中趨勢的指標，平均數(shù)可分為三種：算術平均數(shù)、幾何平均數(shù)、調和平均數(shù)。其中，前兩者在GIS分析中最常用到。算術平均數(shù)代表了數(shù)據(jù)集的平均水平，不受總體范圍的影響，因此可以作為比較分析的指標，還可作為評價事物的客觀標準。如一個地區(qū)某一年的人均收入可以作為衡量該地區(qū)當年經濟狀況的指標；要比較該地區(qū)近幾年的經濟增長水平，也需要用該地區(qū)近幾年的人均收入來比較，因為人均收入能代表該地區(qū)經濟狀況的平均水平，可以用來做比較分析。另外求取離差、平均離差、離差平方和、方差、標準差、變差系數(shù)、偏度系數(shù)和峰度系數(shù)等時，要先求得算術平均數(shù)；算術平均數(shù)也可用于圖像處理中的平滑運算。加權平均數(shù)與算術平均數(shù)的應用是大致相同的，但加權平均數(shù)要考慮各數(shù)據(jù)點的貢獻作用。加權算術平均數(shù)在求算術平均數(shù)時，考慮到數(shù)據(jù)集中的n個值有時會含有不同的比重，對平均數(shù)的影響也就不同。所以用加權平均法來計算其算術平均數(shù)。權重為出現(xiàn)的頻數(shù)，加權平均數(shù)可以按下式計算： 其中，幾何平均數(shù)n個數(shù)據(jù)的連乘積再開n次方所得的方根數(shù)幾何平均數(shù)用于分析和研究平均改變率、平均增長率、平均定比等，還在偏相關系數(shù)里有應用。中位數(shù)若將數(shù)據(jù)值按大小順序排列，位于中間的那個值就是中位數(shù)或稱中值。當數(shù)據(jù)集中有奇數(shù)個數(shù)據(jù)時，數(shù)據(jù)按大小順序排列，那么第位數(shù)就是中位數(shù)；當有偶數(shù)個數(shù)據(jù)時，中位數(shù)為第項與第項的平均數(shù)。中位數(shù)不受極端數(shù)值的影響，如果數(shù)據(jù)集的分布形狀是左右對稱的，則中位數(shù)等于平均數(shù)；當數(shù)據(jù)集的分布形狀呈左偏或右偏，以中位數(shù)表示它們的集中趨勢比算術平均數(shù)更合理。眾數(shù)是數(shù)據(jù)集中最常出現(xiàn)的，因此一定是數(shù)據(jù)集中的某個值，代表了多少意見，不受極端值的影響，在頻數(shù)分布曲線上位居最高點，即曲線的峰值。眾數(shù)常用于投票選舉。若數(shù)據(jù)集的分布并不是明顯集中在某個數(shù)值上，用眾數(shù)來代表集中情形就沒有多大意義，甚至會有誤導作用。眾數(shù)眾數(shù)是數(shù)據(jù)集中出現(xiàn)頻數(shù)（次數(shù)）最多的某個（或某幾個）數(shù)。10.2.2 代表數(shù)據(jù)離散程度的統(tǒng)計量平均數(shù)、中位數(shù)、眾數(shù)在反映總體一般數(shù)量水平的同時，也掩蓋了總體中各單位的數(shù)量差異。所以，只有這些統(tǒng)計量還不能充分說明一個數(shù)列中數(shù)值的分布情況和波動狀態(tài)。有時雖然兩個數(shù)據(jù)集的平均數(shù)相等，但各數(shù)據(jù)分布在平均數(shù)左右的疏密程度卻不相同，也就是它們的離散程度不一樣，為了把一個數(shù)據(jù)集的離散程度表現(xiàn)出來，就需要研究離散度。代表數(shù)據(jù)離散程度的統(tǒng)計量包括最大值、最小值、分位數(shù)、極差、離差、平均離差、離差平方和、方差、標準差、變差系數(shù)等。離散程度越大，數(shù)據(jù)波動性越大，以小樣本數(shù)據(jù)代表數(shù)據(jù)總體的可靠性越低；離散程度越小，則數(shù)據(jù)波動性小，以小樣本數(shù)據(jù)代表數(shù)據(jù)總體的可靠性越高。表10.2列出了表示離散程度各個統(tǒng)計量的含義及在實際應用的作用。表10.2 代表離散程度的統(tǒng)計量及其特點、作用統(tǒng)計量含義特點及作用最大值與最小值把數(shù)據(jù)從小到大排列，最前端的值就是最小值，最后一個就是最大值通過最大、最小值和極差，可以了解數(shù)據(jù)的取值范圍、分散程度，易于計算，容易理解，但他們都易受極端數(shù)值的影響，漠視了其他值的存在，無法精確地反映所有數(shù)據(jù)的分散情形，因此可能會有誤導作用。在地形分析中，極差主要用于求取一定區(qū)域內的高差。對于兩個不同地區(qū)，雖然它們的平均高程相同，但最高點、最低點及高差不同，說明了這兩個地區(qū)的高程分布狀況有差異。極差一個數(shù)據(jù)集的最大值與最小值的差值稱為極差，它表示這個數(shù)據(jù)集的取值范圍分位數(shù)將數(shù)列按大小排列，把數(shù)列劃分為相等個數(shù)的分段，處于分段點上的值就是分位數(shù)。分位數(shù)剔除了數(shù)據(jù)集中極端值的影響，但計算麻煩，且沒有用到數(shù)據(jù)集中的所有數(shù)據(jù)點。分位數(shù)在數(shù)據(jù)分級中應用較多。離差離差表示各數(shù)值與其平均值的離散程度，其值等于某個數(shù)值與該數(shù)據(jù)集的平均值之差d= x-兩個數(shù)據(jù)集的均值相同，但其離差可以有很大的差別，這說明這兩個數(shù)據(jù)集與各自平均值的離散程度不同。平均離差和離差平方和可以克服恒等于零的缺點，還可以把負數(shù)消除，只剩正值，這樣更易于描述離散程度，而且離差平方和得到的結果較大，使離散程度更明顯。離差平方和用于相關分析中求取相關系數(shù)在回歸分析中，對回歸方程進行顯著性檢驗時，需要對原始數(shù)據(jù)進行離差平方和的分解，即把離差平方和分解為剩余平方和與回歸平方和兩部分，這兩部分的比值可以反映回歸方程的顯著性。在趨勢面分析中，對于趨勢面的擬合程度可以用離差平方和來檢驗，其方法也是將原始數(shù)據(jù)的離差平方和分解為剩余平方和與回歸平方和兩部分，回歸平方和的值越大，表明擬合程度越高。平均離差平均離差是把離差取決對值，然后求和，再除以變量個數(shù)離差平方和離差平方和是把離差求平方，然后求和方差方差是均方差的簡稱。它是以離差平方和除以變量個數(shù)而得到的它們是表示一組數(shù)據(jù)對于平均值的離散程度的很重要的指標，為了應用上的方便，常對方差進行開方，即為標準差。方差和標準差都可應用于相關分析、回歸分析、正態(tài)分布檢驗等，還可用于誤差分析、評價數(shù)據(jù)精度、求取變差系數(shù)、偏度系數(shù)和峰度系數(shù)等。標準差還可用于數(shù)據(jù)分級。標準差對方差進行開方，即為標準差變差系數(shù)變差系數(shù)也稱為離差系數(shù)或變異系數(shù)，是標準差與均值的比值，以表示100%式中為變差系數(shù)，其值為百分率；S是標準差；為平均值。變差系數(shù)是用相對數(shù)的形式來刻畫數(shù)據(jù)離散程度的指標，它可以用來衡量數(shù)據(jù)在時間與空間上的相對變化（波動）的程度。變差系數(shù)可用來求算地形高程變異系數(shù)。10.2.3 代表數(shù)據(jù)分布形態(tài)的統(tǒng)計量分布形態(tài)可以從兩個角度考慮，一是數(shù)據(jù)分布對稱程度，另一個是數(shù)據(jù)分布集中程度。前者的測定參數(shù)稱為偏度或偏斜度，后者的測定參數(shù)稱為峰度。偏度和峰度是衡量數(shù)據(jù)分布特征的重要指標。表10.3 代表分布形態(tài)的統(tǒng)計量及其作用統(tǒng)計量含義作用偏度偏度是刻畫數(shù)據(jù)在均值兩側的對稱程度的參數(shù)，用偏度系數(shù)來衡量。標準偏度系數(shù)（g）：g=當g0時，數(shù)據(jù)的分布情況如下圖：f(x)f(x)f(x)f(x) f(x)f(x)f(x)是數(shù)據(jù)分布的密度函數(shù)，是數(shù)據(jù)的平均值偏度可以表示數(shù)據(jù)分布的不對稱性，刻畫出是向正的方向偏還是向負的方向偏（小于或大于）峰度可以表示數(shù)據(jù)頻數(shù)分布曲線峰形的相對高聳程度或尖平程度。這兩個指標主要用于分析數(shù)據(jù)的頻率統(tǒng)計圖以及評價正態(tài)分布性，當g=0且g=0時，數(shù)據(jù)是標準正態(tài)分布。峰度峰度是刻畫數(shù)據(jù)在均值兩側的集中程度的參數(shù)，用峰度系數(shù)來衡量。標準峰度系數(shù) g按下式計算：g=g0gR時，區(qū)域化變量Z(x)的空間相關性不存在，即當某點與已知點的距離大于變程時，該點數(shù)據(jù)不能用于內插或外推。4. Voronoi圖Voronoi地圖是由在樣點周圍形成的一系列多邊形組成的。某一樣點的Voronoi多邊形的生成方法是：多邊形內任何位置距這一樣點的距離都比該多邊形到其它樣點的距離要近。Voronoi 多邊形生成之后，相鄰的點就被定義為具有相同連接邊的樣點。Voronoi 圖中多邊形值可以采用多種分配和計算方法：l 簡化（Simple）：分配到某個多邊形單元的值是該多邊形單元的值；l 平均（Mean）：分配到某個多邊形單元的值是這個單元與其相鄰單元的平均值；l 模式（Mode）：所有的多邊形單元被分為五級區(qū)間，分配到某個多邊形單元的值是這個單元與其相鄰單元的模式（即出現(xiàn)頻率最多的區(qū)間）；l 聚類（Cluster）：所有的多邊形單元被分配到這五級區(qū)間中，如果某個多邊形單元的級區(qū)間與它的相鄰單元的級區(qū)間都有不同，這個單元用灰色表示，以區(qū)別于其它單元；l 熵（Entropy）：所有單元都根據(jù)數(shù)據(jù)值的自然分組分配到這五級中。分配到某個多邊形單元的值是根據(jù)該單元和其相鄰單元計算出來的熵；l 中值（Median）：分配給某多邊形的值是根據(jù)該單元和其相鄰單元的頻率分布計算的中值；l 標準差（StDev）：分配給某多邊形的值是根據(jù)該單元和其相鄰單元計算出的標準差；l 四分位數(shù)間間隔（IQR）：第一和第三四分位數(shù)是根據(jù)某單元和其相鄰單元的頻率分布得出的。分配給某多邊形單元的值是用第三四分位數(shù)減去第一四分位數(shù)得到的差。圖10.7 簡化（Simple）Voronoi地圖圖10.8 熵（Entropy）Voronoi地圖圖10.7為簡化（Simple）Voronoi地圖，圖10.8為熵（Entropy）Voronoi地圖，顯然不同的多邊形賦值方式，獲取的Voronoi地圖提供信息也不同。簡化（Simple）Voronoi圖可以了解到每個采樣點控制的區(qū)域范圍，也可以體現(xiàn)出每個采樣點對區(qū)域內插的重要性。利用簡化（Simple）Voronoi地圖中就可以找出一些對區(qū)域內插作用不大且可能影響內插精度的采樣點值，可以將它剔除。用聚類和熵的方法生成的Vonoroi圖可用來幫助識別可能的離群值。熵值是量度相鄰單元相異性的一個指標。自然界中，距離相近的事物比距離遠的事物具有更大的相似性，因此，局部離群值可以通過高熵值的區(qū)域識別出來。同樣，一般認為某個特定單元的值至少應與它周圍單元中的某一個的值相近。因此聚類方法也能將那些與周圍單元不相同的單元識別出來。10.3.2檢驗數(shù)據(jù)分布在空間統(tǒng)計的分析中，許多統(tǒng)計分析模型，如地統(tǒng)計分析，都是建立在平穩(wěn)假設的基礎上，這種假設在一定程度上要求所有數(shù)據(jù)值具有相同的變異性。另外，一些克里格插值（如普通克里格法、簡單克里格法和泛克里格法等）都假設數(shù)據(jù)服從正態(tài)分布。如果數(shù)據(jù)不服從正態(tài)分布，需要進行一定的數(shù)據(jù)變換，從而使其服從正態(tài)分布。因此，在進行地統(tǒng)計分析前，檢驗數(shù)據(jù)分布特征，了解和認識數(shù)據(jù)具有非常重要的意義。數(shù)據(jù)的檢驗可以通過直方圖和正態(tài)QQPlot分布圖完成。如果數(shù)據(jù)服從正態(tài)分布，數(shù)據(jù)的直方圖應該呈鐘形曲線，在正態(tài)QQPlot圖中，數(shù)據(jù)的分布近似成為一條直線。10.3.3尋找數(shù)據(jù)離群值數(shù)據(jù)離群值分為全局離群值和局部離群值兩大類。全局離群值是指對于數(shù)據(jù)集中所有點來講，具有很高或很低的值的觀測樣點。局部離群值值對于整個數(shù)據(jù)集來講，觀測樣點的值處于正常范圍，但與其相鄰測量點比較，它又偏高或偏低。離群點的出現(xiàn)有可能就是真實異常值，也可能是由于不正確的測量或記錄引起的。如果離群值是真實異常值，這個點可能就是研究和理解這個現(xiàn)象的最重要的點。反之，如果它是由于測量或數(shù)據(jù)輸入的明顯錯誤引起的，在生成表面之前，它們就需要改正或剔除。對于預測表面，離群值可能引起多方面的有害影響，包括影響半變異建模和鄰域分析的取值。離群值的尋找可以通過三種方式實現(xiàn)：1. 利用直方圖查找離群值離群值在直方圖上表現(xiàn)為孤立存在或被一群顯著不同的值包圍。但需注意的是，在直方圖中孤立存在或被一群顯著不同的值包圍的樣點不一定是離群值。2. 用半變異/協(xié)方差函數(shù)云圖識別離群值如果數(shù)據(jù)集中有一個異常高值的離群值，則與這個離群值形成的樣點對，無論距離遠近，在半變異/協(xié)方差函數(shù)云圖中都具有很高的值。3. 用Voronoi圖查找局部離群值用聚類和熵的方法生成的Voronoi圖可用來幫助識別可能的離群值。熵值是量度相鄰單元相異性的指標。通常，距離近的事物比距離遠的事物具有更大的相似性。因此，局部離群值可以通過高熵值的區(qū)域識別出來。同理，聚類方法也可將那些與它們周圍單元不相同的單元識別出來。圖10.9 直方圖查找離群值圖如圖10.9，直方圖最右邊被選中的一個柱狀條即是該數(shù)據(jù)的離群值。相應地，數(shù)據(jù)點層面上對應的樣點也被刷光。10.3.4全局趨勢分析（a）圖10.10 趨勢面分析透視面圖（b）通常一個表面主要由兩部分組成：確定的全局趨勢和隨機的短程變異?？臻g趨勢反映了空間物體在空間區(qū)域上變化的主體特征，它主要揭示了空間物體的總體規(guī)律，而忽略局部的變異。趨勢面分析是根據(jù)空間抽樣數(shù)據(jù)，擬合一個數(shù)學曲面，用該數(shù)學曲面來反映空間分布的變化情況。它可分為趨勢面和偏差兩大部分，其中趨勢面反映了空間數(shù)據(jù)總體的變化趨勢，受全局性、大范圍的因素影響。如果能夠準確識別和量化全局趨勢，在空間分析統(tǒng)計建模中就可以方便的剔除全局趨勢，從而能更準確地模擬短程隨機變異。透視分析是探測全局趨勢常用方法，準確的判定趨勢特征關鍵在于選擇合適的透視角度。同樣的采樣數(shù)據(jù)，透視角度不同，反映的趨勢信息也不相同。圖10.10（a）為顯示某地區(qū)東西方向（X軸）和南北方向（Y軸）的高程趨勢圖。圖10.10（b）逆時針旋轉45度后，顯示東南-西北方向和西南-東北方向的高程趨勢圖。趨勢分析過程中，透視面的選擇應盡可能使采樣數(shù)據(jù)在透視面上的投影點分布比較集中，通過投影點擬合的趨勢方程才具有代表性，才能有效反映采樣數(shù)據(jù)集全局趨勢。顯然，圖10.10（a）反映的趨勢比圖10.10（b）更為準確。10.3.5空間自相關及方向變異大部分的地理現(xiàn)象都具有空間相關特性，即距離越近的兩事物越相似。這一特性也是空間地統(tǒng)計分析的基礎。半變異/協(xié)方差函數(shù)云圖就是這種相似性的定量化表示?？臻g自相關分析包括全程空間自相關分析和局部空間自相關分析，自相關分析的結果可用來解釋和尋找存在的空間聚集性或“焦點”。 空間自相關分析需要的空間數(shù)據(jù)類型是點或面數(shù)據(jù)，分析的對象是具有點/面分布特征的特定屬性。全程空間自相關用來分析在整個研究范圍內指定的屬性是否具有自相關性。局部空間自相關用來分析在特定的局部地點指定的屬性是否具有自相關性。具有正自相關的屬性，其相鄰位置值與當前位置的值具有較高的相似性。下面介紹兩個常用的分析空間自相關的參數(shù)：MoransI和GearyC。1. 空間權重矩陣地理事物在空間上的此起彼伏和相互影響是通過它們之間的相互聯(lián)系得以實現(xiàn)的，空間權重矩陣是傳載這一作用過程的實現(xiàn)方法。因此，構建空間權重矩陣是研究空間自相關的基本前提之一?？臻g數(shù)據(jù)中隱含的拓撲信息提供了空間鄰近的基本度量。通常定義一個二元對稱空間權重矩陣Wnn來表達n個空間對象的空間鄰近關系，可根據(jù)鄰接標準或距離標準來度量，還可以根據(jù)屬性值xj和二元空間權重矩陣來定義一個加權空間鄰近度量方法?？臻g權重矩陣的表達形式為 （10.1）根據(jù)鄰接標準，當空間對象i和空間對象j相鄰時，空間權重矩陣的元素wij為1，其他情況為0，表達式如下： （10.2）根據(jù)距離標準，當空間對象i和空間對象j在給定距離d之內時，空間權重矩陣的元素wij 為1，否則為0，表達式為 （10.3）如果采用屬性值xj和二元空間權重矩陣來定義一個加權空間鄰近度量方法，則對應的空間權重矩陣可以定義如下： （10.4）2. MoransI參數(shù)MoransI是應用最廣的一個參數(shù)。對于全程空間自相關，MoransI定義是：Morans （10.5）對于局部位置i的空間自相關，MoransI定義是： （10.6）其中：n 是觀察值的數(shù)目；xi 是在位置i的觀察值；Zi 是xi 的標準化變換，。，。是對稱的空間權重矩陣，如果i與j相鄰，取值為，否則取值為。按照行和歸一化后的權重矩陣（每行的和為1），非對稱的空間權重矩陣。MoransI值介于-1到1之間，0為不相關。按照假定的空間數(shù)據(jù)分布可以計算Morans I的期望值和期望方差。對于正態(tài)分布假設， （10.7） （10.8）對于隨機分布假設， （10.9） （10.10）其中， ，是第行權重值之和，是第列權重值之和。原假設是沒有空間自相關。根據(jù)下面標準化統(tǒng)計量參照正態(tài)分布表可以進行假設檢驗。 （10.11）MoransI如果是正的而且顯著，表明具有正的空間相關性。即在一定范圍內各位置的值是相似的，如果是負值而且顯著的，則具有負的空間相關性，數(shù)據(jù)之間不相似。接近于0則表明數(shù)據(jù)的空間分布是隨機的，沒有空間相關性。3. GerayC參數(shù)對于全局空間自相關： （10.12）對于局部位置i的空間自相關： （10.13）其中，Wij是空間權重矩陣。C的值總是正的。假設檢驗是如果沒有空間自相關，C的均值為1。顯著性的低值（0和1之間）表明具有正的空間自相關，顯著性的高值（大于1）表明具有負的空間自相關。10.4分級統(tǒng)計分析分級是對數(shù)據(jù)進行加工處理的一種重要方法，通過分級可以把數(shù)據(jù)劃分成不同的級別，體現(xiàn)數(shù)據(jù)自身的特征，為應用研究及專題制圖提供基礎。10.4.1 分級的概念與目的數(shù)據(jù)分級根據(jù)一定的方法或標準把數(shù)據(jù)分成不同的級別，也就是把一個數(shù)據(jù)集劃分成不同的子集，在此過程中，還可設置分級精度和分級數(shù)目等。數(shù)據(jù)分級之后，僅使原來的數(shù)據(jù)重新歸類，數(shù)據(jù)的屬性沒有發(fā)生改變，研究人員可以根據(jù)分級后的數(shù)據(jù)進行下一步的應用分析。數(shù)據(jù)分級的根本目的在于區(qū)分數(shù)據(jù)集中個體的差別，分級統(tǒng)計的過程就是區(qū)別個體性質的過程。分級的應用目的有兩點：一是為了分級后，圖面制圖效果好，有利于用戶讀圖；二是用不同的分級方法來突出顯示制圖區(qū)域內不同的地貌特征。10.4.2 分級的原則 科學性原則：分級指標的確定要遵循一定的科學規(guī)律； 完整性原則：整個數(shù)據(jù)集中的所有數(shù)據(jù)都應被分到不同的級別中，沒有遺漏，而且同一數(shù)據(jù)集中的每一數(shù)據(jù)只能被分到某一級別中，不能同時分到兩個或多個級別中； 適用性原則：對于一個數(shù)據(jù)集，應該根據(jù)研究或應用的需要選擇合適的分級方法，使得分級結果能較好地滿足目的； 美觀性原則：分級方法及分級數(shù)目的確定不但要依據(jù)研究目的，還要注重制圖效果，制作專題地圖不僅要體現(xiàn)數(shù)據(jù)的空間分布特征，還要使得圖面色彩平衡，特征明顯，易于理解。10.4.3 分級統(tǒng)計的方法分級方法多種多樣，在應用時應根據(jù)研究的需要選擇合適的方法來突出需要的數(shù)據(jù)信息。分級方法的種類也很多，本書主要介紹以下三種： 按使用分級方法的多少可分為單一分級法和復合分級法：單一分級是指對于一個數(shù)據(jù)集只用了一種分級方法；復合分級是指由于數(shù)據(jù)自身的特點，需要對一部分數(shù)據(jù)使用某種分級方法，對另一部分數(shù)據(jù)使用另外一種分級方法，才能更好地滿足研究的需要。如一組坡度數(shù)據(jù)，一部分較?。ㄆ旅嫫骄彛硪徊糠趾艽螅ǖ貏荻盖停?，對這兩部分數(shù)據(jù)，就應選用兩種不同的分級方法，才能更好地突出變化特征。自定義分級模式分級等間距分級分位數(shù)分級等面積分級標準差分級自然裂點法其它分級方法有系統(tǒng)的不等間距分級面積正態(tài)分布分級嵌套平均值分級 按級差是否相等可分為等值分級法和不等值分級法：等值分級法又可以分為等面積分級、等間距分級、分位數(shù)分級等：不等值分級法可以分為自然裂點法、標準差分級、平均值嵌套分級等。 按確定級差的方法可分為自定義分級法和模式分級法。本節(jié)就是按這種分類體系展開討論的，其結構圖如下： 圖10.11 分級方法的類型1. 自定義分級自定義分級即對一個數(shù)據(jù)集，根據(jù)自己的應用目的設定各個級別的數(shù)值范圍來實現(xiàn)分級的方法。這種方法適用于研究者對該數(shù)據(jù)集比較了解，能夠找到合適的分級臨界點。在自定義分級中，臨界點的選擇非常重要，臨界點選擇得好，就能夠增強同一級別區(qū)域間的同質性分級和各級之間的差異性，分級結果就能夠很好地滿足各種分析需求。如在對坡度進行分級的過程中，應根據(jù)應用目標的要求, 確定臨界坡度。2. 模式分級模式分級就是指按固定模式進行分級，在固定模式中，級差由特定的算法自動設定。模式分級分為等間距分級、分位數(shù)分級、等面積分級、標準差分級、自然裂點法分級等。 等間距分級等間距分級是一種最簡單的分級方法，它按某個恒定間隔來對數(shù)據(jù)進行分級。假定數(shù)據(jù)集里有最大值和最小值，那么間距D=。圖10.12 等間距分級示意圖圖10.13 分位數(shù)分級示意圖等間距方法原理簡單、易操作，但當數(shù)據(jù)集中在某一小范圍內時，各分級之間數(shù)據(jù)個數(shù)的差別太大會造成圖面配置不均衡，影響了制圖效果?？梢?，當數(shù)據(jù)具有均勻變化的分布特征時，等間距分級法就簡明實用；若數(shù)據(jù)分布差異過大，將會以響制圖與對統(tǒng)計結果的分析。圖10.12為某一地區(qū)溫度數(shù)據(jù)的等間距分級示意圖。 分位數(shù)分級分位數(shù)分級是把數(shù)列劃分為相等個數(shù)的分段，根據(jù)實際需要選擇四分位、五分位、六分位十分位。為此，要先將數(shù)列按大小排列，從一端開始計算其分位數(shù)，把處于分位數(shù)上的那個值作為分級值。分位數(shù)分級可以使每一級別的數(shù)據(jù)個數(shù)接近一致，往往能產生較好的制圖效果。圖10.13為某地區(qū)溫度數(shù)據(jù)的分位數(shù)分級示意圖。 等面積分級對于矢量數(shù)據(jù)而言，對每個多邊形都可以求出其面積，將面積的值按大小順序排列并累加，把累加面積分為相等的幾段，作為分級范圍，這樣，每個級別中包含的樣本數(shù)目雖然不同，但總面積基本一致。等面積方法使得每一級在圖上占據(jù)的面積相等（或大致相等）。這種方法的特點是在圖面上只反映各級占有相同的面積，制圖效果好，但是沒有充分利用圖面表示級間的差異。對于規(guī)則柵格數(shù)據(jù)而言，一定區(qū)域內的面積可由該區(qū)域內的柵格個數(shù)乘以柵格分辨率得到，所以按等面積分級只需考慮柵格個數(shù)即可。這時可以將數(shù)據(jù)按大小順序排列，將數(shù)據(jù)個數(shù)累加，并把累加的個數(shù)分為相等的幾段，這與分位數(shù)分級法得到的分級結果基本上是一致的。 標準差分級標準差可以反映各數(shù)據(jù)間的離散程度，按標準差分級，首先要保證數(shù)據(jù)的分布具有正態(tài)分布的規(guī)律，才可計算平均值和標準差Std.Dev，然后根據(jù)數(shù)據(jù)波動情況劃分等級。以算術平均值作為中間級別的一個分界點，以一倍標準差參與分級時其余分界點為：Std.Dev，2Std.Dev，3Std.DeviStd.Dev。當然也可以采用1/2倍標準差參與分級，即1/2Std.Dev，2/2Std.Dev，3/2Std.Devi/2Std.Dev的形式，也可以用1/3倍標準差或1/4倍標準差。顯然，分級數(shù)目是由數(shù)據(jù)本身所決定的，且對于同一數(shù)據(jù)集，采用一倍標準差時，分級數(shù)目最少，采用1/4倍標準差時，分級數(shù)目最多。圖10.14為某地區(qū)溫度數(shù)據(jù)的標準差分級示意圖。 自然裂點法分級任何統(tǒng)計數(shù)列都存在一些自然轉折點、特征點，用這些點可以把研究的對象分成性質相似的群組，因此，裂點本身就是分級的良好界限。將統(tǒng)計數(shù)據(jù)制成頻率直方圖、坡度曲線圖、積累頻率直方圖，都有助于找出數(shù)據(jù)的自然裂點。如果頻率最低點與峰值構成一個近似正態(tài)分布曲線，可以把任意兩個正態(tài)分布曲線交點作為分級界線。圖10.15 自然裂點法分級示意圖圖10.14 標準差分級示意圖自然裂點法基本上是基于讓各級別中的變異總和達到最小的原則來選擇分級斷點的。由于自然裂點法是根據(jù)數(shù)據(jù)本身特征找出一些自然裂點進行分級的。圖10.15為某地區(qū)溫度數(shù)據(jù)的自然裂點法分級示意圖。 其他分級方法a. 有規(guī)律的不等間距分級這種方法與等間距分級法的區(qū)別在于它的間距是按一定規(guī)律變化的，而不是一個恒定的間隔。該方法采用的間隔或級差有算術級數(shù)和幾何級數(shù)兩種，每種又都可通過以下六種變化方法來確定各級的分級間隔：按某一恒定速率遞增、按某一加速度遞增、按某一減速度遞增、按某一恒定速率遞減、按某一加速度遞減、按某一減速度遞減。b. 按嵌套平均值分級該方法先計算整個數(shù)據(jù)集的平均值，它將數(shù)據(jù)集分為兩部分，每部分中再計算平均值，又各自把所有的那一部分分成兩段，以此類推，就可以把數(shù)據(jù)集區(qū)分為2個等級，即2的幾何級數(shù)。n是計算中的平均值的嵌套序數(shù)，用這種方法只能得到偶數(shù)個級別，而不可能得到奇數(shù)個級別。c. 按面積正態(tài)分布分級按數(shù)據(jù)的大小排列，累加其面積，然后按正態(tài)分布的規(guī)則使中間級別所占的面積較大，往高端和低端的級別中所占的面積都依次減小，并由此來確定每級的分界線。顯然，這種方法不僅使每個級別中樣本的數(shù)目不相等，而且各級別的累加面積呈正態(tài)分布。總之，關于數(shù)據(jù)的統(tǒng)計分級的研究還很多，其目的都在于改善分級間隔的規(guī)則性、同級之中的同質性和不同級別之間的差異性等等。10.5 空間插值空間數(shù)據(jù)插值是進行數(shù)據(jù)外推的基本方法。常用的插值方法有很多，分類并沒有統(tǒng)一的標準，例如從數(shù)據(jù)分布規(guī)律來講，有基于規(guī)則分布數(shù)據(jù)的內插方法、基于不規(guī)則分布的內插方法和適合于等高線數(shù)據(jù)的內插方法等；從內插函數(shù)與參考點的關系方面，又分為曲面通過所有采樣點的純二維插值方法和曲面不通過參考點的曲面擬合插值方法；從內插曲面的數(shù)學性質來講，有多項式內插、樣條內插、最小二乘配置內插等內插函數(shù)；從對地形曲面理解的角度，內插方法有克立金法、 多層曲面疊加法、加權平均法、分形內插等；從內插點的分布范圍，內插方法分為整體內插、局部內插和逐點內插法?？臻g內插的根本是對空間曲面特征的認識和理解，具體到方法上，則是內插點鄰域范圍的確定、權值確定方法（自相關程度）、內插函數(shù)的選擇等三方面的問題。由于每一種內插方法都有其自身的特點和適用范圍，了解方法的特點是本質所在。本書并不打算對各種內插算法從數(shù)學實現(xiàn)方法上進行分析討論，而是從內插范圍分類方法入手，對每一類內插方法的特點進行簡要的分析歸納，同時為保證內容上的完整性和連續(xù)性，在本節(jié)只介紹概念上的東西，而具體的實現(xiàn)方法與相應內容請查閱相關文檔。10.5.1整體內插整體內插，就是在整個區(qū)域用一個數(shù)學函數(shù)來表達地形曲面，如圖10.16所示。整體內插函數(shù)通常是高次多項式，要求地形采樣點的個數(shù)大于或等于多項式的系數(shù)數(shù)目。當?shù)匦尾蓸狱c的個數(shù)與多項式的系數(shù)相等時，這時能得到一個唯一的解，多項式通過所有的地形采樣點，屬純二維插值；而當采樣點個數(shù)多于多項式系數(shù)時，沒有唯一解，這時一般采用最小二乘法求解，即要求多項式曲面與地形采樣點之間差值的平方和為最小，屬曲面擬合插值或趨勢面插值。從數(shù)學角度講，任何復雜的曲面都可用多項式在任意精度上逼近，但由于以下原因，在空間內插中整體內插并不常用：曲面擬合Z=f (x,y )內插格網點高程圖10.16 整體內插l 整體內插函數(shù)保凸性較差，高次多項式雖然能在任意精度上逼近地形曲面，并且能使數(shù)學曲面與實際地形曲面有更多的重合點，但由于多項式是自變量冪函數(shù)的和的形式，采樣點的增減或移動都需要對多項式的系數(shù)作全面調整，從而采樣點之間會出現(xiàn)難以控制的振蕩現(xiàn)象，致使函數(shù)極不穩(wěn)定，從而導致保凸性較差；l 不容易得到穩(wěn)定的數(shù)值解解，高次多項式的系數(shù)求解一般要解算較高階的線性方程組，計算的舍入誤差和數(shù)據(jù)采樣誤差（平面位置誤差），都有可能引起多項式系數(shù)發(fā)生較大變化，使高次多項式不容易得到穩(wěn)定的數(shù)值解。而且這種微小的數(shù)據(jù)擾動影響，在二元高次多項式中更為嚴重，因為在一元函數(shù)插值時，自變量誤差僅在一個方向上（x軸），而在二維平面上，采樣點的偏移方向卻是無限的。l 多項式系數(shù)物理意義不明顯，在低階多項式中，各個系數(shù)的物理意義非常明確，例如線性多項是中，a、b分別為兩個坐標軸方向的斜率，而在高次多項式中，各個系數(shù)的物理意義一般不明確，容易導致無意義的地形起伏現(xiàn)象。l 解算速度慢且對計算機容量要求較高。l 不能提供內插區(qū)域的局部地形特征；a. 二次趨勢擬合b. 三次趨勢擬合圖10.17 江蘇某地區(qū)GDP全局多項式趨勢擬合整體內插雖然有如上的缺點，但其優(yōu)點也是明顯的，例如整個區(qū)域上函數(shù)的唯一性、能得到全局光滑連續(xù)的空間曲面、充分反映宏觀地形特征等。整體內插函數(shù)常常用來揭示整個區(qū)域內的地形宏觀起伏態(tài)勢。在空間內插中，一般是與局部內插方法配合使用，例如在使用局部內插方法前，利用整體內插去掉不符合總體趨勢的宏觀地物特征。圖10.17（a）和圖10.17（b）分別為江蘇某地區(qū)GDP二次擬合及三次擬合的全局多項式擬合圖。顯然，次數(shù)越低，擬合的表面越粗糙，實際表面擬合的效果越差，大致代表了此區(qū)域的宏觀趨勢；次數(shù)越高，擬合面越光滑，擬合的結果更接近實際的表面。但并不是次數(shù)越高越好，次數(shù)過高使得計算量大大增加而精度提高不大，一般選用到三次即可。10.5.2局部分塊內插利用數(shù)學曲面來模擬實際地形表面，是地形表達的一個常用的手段。例如一階線性平面可模擬具有單一坡度的斜坡地形表面，二次曲面方程可表達山頭、洼地區(qū)域，而三次曲面則能描述較為復雜的地形曲面。然而低階多項式雖然可表達各種地形曲面，但一個地區(qū)確常常包含各種復雜的地貌形態(tài)，簡單的曲面并不能很好的表達這些地形曲面。理論上任何復雜的曲面都可用多項式進行逼近，但高階多項式的上述缺點，也不是理想的地形描述工具。解決這類問題的辦法就是采取分而治之的辦法，即將復雜的地形地貌分解成一系列的局部單元，在這些局部單元內部地形曲面具有單一的結構，由于范圍的縮小和曲面形態(tài)的簡化，用簡單曲面就可較好的描述地形曲面。將地形區(qū)域按一定的方法進行分塊，對每一塊根據(jù)地形曲面特征單獨進行曲面擬合和高程內插，稱為空間分塊內插（圖10.18）。區(qū)域分塊簡化了地形的曲面形態(tài)，使得每一塊都可用不同的曲面進行表達，但隨之而來的是如何進行分塊和如何保證各個分塊之間的曲面的連續(xù)性。一般的可按地形結構線或規(guī)則區(qū)域進行分塊，而分塊大小取決于地形的復雜程度、地形采樣點的密度和分布；為保證相鄰分塊之間的平滑連接，相鄰分塊之間要有一定寬度的重疊，另外一種分塊之間的平滑連接是對內插曲面補充一定的連續(xù)性條件。圖10.18 局部分塊內插方法塊內逐點內插格網點高程分塊已知點格網點不同的分塊單元可用不同的內插函數(shù)，常用的內插數(shù)函數(shù)有線性內插、雙線性內插、多項式內插、樣條函數(shù)、多層曲面疊加法等。l 線性內插和雙線性內插形如的多項式稱為線性平面，它將分塊單元內部的地形曲面視為平面。如果在線性多項式中增加了交叉項，線性內插則變成雙線性內插函數(shù)：，之所以稱為雙線性內插，是因為當y為常數(shù)時，表達的是x方向的線性函數(shù)，而當x為常數(shù)時，則為y方向的線性函數(shù)。線性內插函數(shù)中有三個未知數(shù)，需要三個采樣點才能唯一確定，而雙線性內插函數(shù)中有四個未知數(shù)，需要四個已知點。線性內插和雙線性內插函數(shù)由于物理意義明確，計算簡單，是基于TIN和基于正方形格網分布采樣數(shù)據(jù)的DEM內插和分析應用的最常用的方法。l 二元樣條函數(shù)內插所謂樣條曲面，就是將一張具有彈性的薄板壓定在各個采樣點上，而其它的地方自由彎曲。從數(shù)學上講，就是一個分段的低次多項式，多項式的次數(shù)一般不超過三階。通過樣條函數(shù)，可以獲取在各個采樣點上具有最小曲率的擬合曲面。二元樣條函數(shù)首先對采樣區(qū)域進行分塊，對每一塊用一個多項式進行擬合，為保證各個分塊之間的平滑過渡，按照彈性力學條件設立分塊之間的連續(xù)性條件，即公共邊界上的導數(shù)連續(xù)條件。雖然樣條函數(shù)可適合的任意形狀的分塊單元，但一般還是將其應用在規(guī)則格網分布的采樣數(shù)據(jù)中。與整體內插函數(shù)相比較，樣條函數(shù)不但保留了局部地形的細部特征，還能獲取連續(xù)光滑的DEM。同時樣條函數(shù)在擬合時，由于多項式的階數(shù)比較低，對數(shù)據(jù)誤差的響應不敏感，具有較好的保凸性和逼真性，同時也有良好的平滑性。圖10.19為江蘇某地區(qū)GDP的張力樣條函數(shù)插值結果，圖10.19（A）取權重為0，圖10.19（B）取權重為5，圖10.19（C）為兩次插值結果的差值的絕對值。可以看出，權重越高，表面越粗糙，權重不同，局部地區(qū)插值結果變化比較劇烈。樣條函數(shù)將地表分塊視為彈性剛體，采取具有彈性力學條件的光滑連續(xù)條件。然后地形并不是一個狹義的剛體，也不具備滿足彈性力學光滑性條件。因此雖然樣條函數(shù)具有嚴密的理論基礎，但未必是數(shù)字地形內插的理想數(shù)學模型。圖10.19 張力樣條函數(shù)插值結果A： Weight = 0B： Weight = 5C：