空間分析之疊加分析第一頁，共三十三頁，編輯于2023年，星期一2023/6/31一、空間疊置概念空間疊置分析(SpatialOverlayAnalysis)是指在統(tǒng)一空間參照系統(tǒng)條件下，每次將同一地區(qū)兩個地理對象的圖層進行疊置，以產生空間區(qū)域的多重屬性特征，或建立地理對象之間的空間對應關系?？臻g疊置分析目的：尋找和確定同時具有幾種地理屬性的地理要素的分布，或是按照確定的地理指標，對疊加后產生的具有不同屬性級的多邊形進行分類或分級第二頁，共三十三頁，編輯于2023年，星期一2023/6/3二、疊置分析的分類將同一地區(qū)，同一比例尺的兩組或更多的多邊形要素的數據層進行疊置，會出現兩種成果形式，一種是根據兩組多邊形邊界的交點來建立具有多重屬性的多邊形，稱為合成疊置；另一種是進行多邊形范圍的屬性特征的統(tǒng)計分析，稱為統(tǒng)計疊置。合成疊置得到一張新的疊置圖，產生了許多新的多邊形，每個多邊形內都具有兩種以上的屬性，通過區(qū)域多重屬性的模擬，尋找和確定同時具有幾種地理屬性的分布區(qū)域。統(tǒng)計疊置得到的是新的屬性數據，得到的新屬性值可能與原數據層的屬性意義不一致，因為得到的結果可能是名稱、分類值，也可能是分析后的結果。第三頁，共三十三頁，編輯于2023年，星期一2023/6/3疊置分析的數學基礎（空間邏輯運算）疊加過程往往是對空間信息和對應的屬性信息作集合的交、并、差、余運算，也可再進一步對屬性作其他的數學運算。為討論方便將空間圖層A，B，C定義為二值圖象

3、空間邏輯差運算；A-B=XX∈A且

X∈B1、空間邏輯并（或）運算；A∪B=XX∈A或X∈B2、空間邏輯交（與）運算；A∩B=XX∈A且X∈B4、空間包含；AB第四頁，共三十三頁，編輯于2023年，星期一2023/6/3例：土壤厚度（大于50厘米）和土壤類型（紅壤和其他類型）兩個二值化圖層，不同的邏輯運算結果如下：AND關系：結果是將土層厚度大于50厘米，且土壤為紅壤的土壤單元顯示出來；OR關系：結果將土層厚度大于50厘米，或者土壤為紅壤的土壤單元顯示出來；XOR：結果將土層厚度小于50厘米，或者土壤不是紅壤的土壤單元顯示出來；NOT:如結果是將土層厚度大于50厘米，但土壤不是紅壤的土壤單元顯示出來；邏輯關系運算例第五頁，共三十三頁，編輯于2023年，星期一2023/6/3柵格數據的算術運算1)算術運算指兩層以上的對應網格值經加、減運算，而得到新的柵格數據系統(tǒng)的方法。

（4）疊加分析的數學運算疊加法第六頁，共三十三頁，編輯于2023年，星期一2023/6/3算術運算——以灰度(亮度)賦值柵格為例（乘除運算）10505010020100100200×2＝（圖像反差增大）陽光砂巖頁巖花崗巖陽坡陰坡波段MSS4MSS5MSS4/MSS5陽坡34510.67陰坡22340.65砂巖對陽光的反射情況乘法案例除法案例0.670.65第七頁，共三十三頁，編輯于2023年，星期一2023/6/3無數據置換①算術運算——賦值柵格的算術條件疊加比較后取最大值比較后取最小值第八頁，共三十三頁，編輯于2023年，星期一2023/6/3①算術運算——賦值柵格的算術條件疊加數據置換第九頁，共三十三頁，編輯于2023年，星期一2023/6/32)函數運算：指兩個以上層面的柵格數據系統(tǒng)以某種函數關系作為復合分析的依據進行逐網格運算，從而得到新的柵格數據系統(tǒng)的過程。這種復合疊置分析方法被廣泛地應用到地學綜合分析、環(huán)境質量評價、遙感數字圖像處理等領域中。利用土壤侵蝕通用方程式計算土壤侵蝕量時，就可利用多層面柵格數據的函數運算復合分析法進行自動處理。一個地區(qū)土壤侵蝕量的大小是降雨(R)、植被覆度(C)、坡度(S)、坡長（L）、土壤抗蝕性(SR)等因素的函數

在ArcView中，使用MapCalculator可以很方便地實現柵格圖層的復合（疊加）運算土壤侵蝕多因子函數運算復合分析示意圖函數運算露點第十頁，共三十三頁，編輯于2023年，星期一2023/6/31、視覺信息疊加：將不同專題的內容疊加顯示在結果圖件上，視覺信息疊加之后，參加疊加的平面之間沒發(fā)生任何邏輯關系，仍保留原來的數據結構。面狀圖、線狀圖和點狀圖之間的疊加；面狀圖區(qū)域邊界之間或一個面狀圖與其他專題區(qū)域邊界之間的疊加；遙感影像與專題地圖的疊加；專題地圖與數字高程模型疊加顯示立體專題圖；遙感影像與DEM疊加生成真三維地物景觀。

疊加分析分為以下五類：視覺信息疊加；點與多邊形疊加；線與多邊形疊加；多邊形疊加；柵格圖層疊加

面狀圖、線狀圖和點狀圖之間的疊加面狀圖與專題區(qū)域邊界之間的疊加遙感影像與專題地圖的疊加專題地圖與數字高程模型疊加遙感影像與DTM疊加第十一頁，共三十三頁，編輯于2023年，星期一2023/6/32、點與多邊形的疊加疊加圖層：將一個含有點的圖層（目標圖層）疊加在另一個含有多邊形的圖層（操作圖層）上，以確定每個點落在哪個區(qū)域內。例如：一個中國政區(qū)圖（多邊形）和一個全國礦產分布圖（點），二者經疊加分析后，并且將政區(qū)圖多邊形有關的屬性信息加到礦產的屬性數據表中，然后通過屬性查詢，可以查詢指定省有多少種礦產，產量有多少；而且可以查詢指定類型的礦產在哪些省里有分布等信息。例如將水井與規(guī)劃區(qū)圖層相疊置，可確定每口井所屬的規(guī)劃區(qū)范圍。*1*3*4*2ADBC點屬性1234多邊形屬性1屬性2ABCD點多邊形點屬性面屬性1面屬性21A2D3C4B第十二頁，共三十三頁，編輯于2023年，星期一2023/6/33、線與多邊形的疊加疊加圖層：將線的圖層（目標圖層）疊加在多邊形的圖層（操作圖層）上，以確定一條線落在哪個多邊形內。疊加原理：與前面不同的是，往往一個線目標跨越多個多邊形，這時需要先進行線與多邊形的求交，并將線目標進行切割，形成一個新的空間目標（新的線目標）的結果集例如：當線狀圖層為河流，疊加的結果是多邊形將穿過它的所有河流打斷成弧段，可以查詢任意多邊形內的河流長度，進而計算它的河流密度等；如果線狀圖層為道路網，疊加的結果可以得到每個多邊形內的道路網密度，內部的交通流量，進入、離開各個多邊形的交通量，相鄰多邊形之間的相互交通量。第十三頁，共三十三頁，編輯于2023年，星期一2023/6/3例如：確定某一行政區(qū)內各種等級道路的里程數線與多邊形的疊合的目的是確定某一線狀圖層上的弧段落在另一多變性圖層上的哪個多邊形內，以便為圖層的每條弧段建立新的屬性。路網密度？——疊加！各種道路長=830KM區(qū)域面積=2797KM2結果？第十四頁，共三十三頁，編輯于2023年，星期一2023/6/34、多邊形與多邊形的疊加多邊形與多邊形的疊置是指將兩個不同圖層的多邊形要素相疊加，根據兩組多邊形的交點來建立多重屬性的多邊形或進行多邊形范圍內的屬性特征的統(tǒng)計分析。多邊形疊置過程分幾何求交過程和屬性確定過程，算法的核心是多邊形求交。①對兩個多邊形進行邊界求交和弧段分割運算，并以新弧段為單位重建拓撲關系；②判斷重建多邊形落在原始多邊形層的哪個多邊形內，從而建立新疊置多邊形與原始多邊形的關系，并抽取屬性。第十五頁，共三十三頁，編輯于2023年，星期一2023/6/3多邊形的不同疊加方式取本底多邊形取上疊多邊形兩個多邊形相交？兩個多邊形是包含關系？多邊形求交，重建拓撲關系屬性抽取，建立關聯連接上疊多邊形結束？本底多邊形結束？ynnn建立包含關系屬性抽取yyy結束多邊形的疊置流程圖第十六頁，共三十三頁，編輯于2023年，星期一2023/6/3新多邊形的屬性多邊形之間的疊置ID屬性101AID屬性1X2Y3Z新多邊形ID層1多邊形屬性層2多邊形屬性10A20X3AX40X5A06A07AY8AZ90Y100Z110Y12Z0AXYZ379101456112層1層2新層12坡度圖層肥力圖層疊加？第十七頁，共三十三頁，編輯于2023年，星期一2023/6/3多邊形疊加產生碎屑多邊形由于矢量有限精度原因，幾何對象不可能完全匹配，出現碎屑多邊形，通常設模糊容限消除第十八頁，共三十三頁，編輯于2023年，星期一2023/6/35、基于柵格數據的疊置分析

特點：柵格數據的疊置算法，雖然數據存貯量比較大，但運算過程比較簡單。變換方法：（1）點變換（2）區(qū)域變換方法（3）鄰域變換方法基于柵格數據的疊置分析是參與分析的兩個圖層的要素均為柵格數據。柵格疊置的作用包括以下幾種：①類型疊置；②數量統(tǒng)計③動態(tài)分析④益本分析⑤幾何提取第十九頁，共三十三頁，編輯于2023年，星期一2023/6/3（1）點變換點變換只依據參與疊置圖層相應點的屬性值進行新的運算，既與各圖層的鄰域點的屬性無關，也不受區(qū)域內一般特征的影響。運算方法包括：算術運算，指數運算，三角函數運算等244633463377225423643366-=0-21-2-102-200-1-1點變換示意結果特征：運算后得到的新屬性值可能與原圖層的屬性意義完全不同。第二十頁，共三十三頁，編輯于2023年，星期一2023/6/3(2)區(qū)域變換方法指在計算新圖層相應的屬性值時，不僅與原圖層對應的柵格的屬性值有關，而且要顧及原圖層所在區(qū)域的集合特征（區(qū)域長度、面積、周長等）。(3)鄰域變換方法指在計算新圖層相應的屬性值時，不僅考慮原圖層對應的柵格及其屬性，而且還應顧及與該柵格相關聯的鄰域或者影響半徑內的柵格屬性值的影響。第二十一頁，共三十三頁，編輯于2023年，星期一2023/6/3疊加分析與視覺信息疊加的區(qū)別疊置分析是把參與疊加的各數據層經過算術的、幾何的、邏輯的運算生成新的數據層；而視覺信息疊加僅是把不同的數據層同時顯示到一個屏幕上或繪制在一幅圖形里，各數據層間不存在邏輯關系，沒有生成新的數據文件，在同一區(qū)域內，地圖投影與比例尺相同便可進行視覺信息復合。疊置分析和空間集合分析的區(qū)別疊置分析和空間集合分析雖然都可以在基于矢量的地理信息系統(tǒng)和基于柵格的地理信息系統(tǒng)中完成，但由于基于柵格的系統(tǒng)的疊加分析要容易和快得多，而矢量數據在疊置操作時容易產生多邊形誤差且邊界搜索也很困難，運算速度又慢等原因，所以一般用柵格數據進行空間集合分析。第二十二頁，共三十三頁，編輯于2023年，星期一2023/6/3101060601010206030303060101060601010306030303060-=0000001000000空間疊加分析實例1：土地利用變化區(qū)域探測80年遙感影像90年遙感影像點變換后影像點變換后影像分析通過80和90年兩期影像的相減運算后得到變換影像，如果：變換影像值＝0；說明該區(qū)未發(fā)生變化變換影像值≠0；說明該區(qū)已發(fā)生變化10耕地20居民點30水域40草地50未利用地60林地Legend注意：此處的遙感影像可以是分類結果，也可以是原始的遙感影像。在一般應用中，多使用原始的遙感影像，可提高變化探測速度。第二十三頁，共三十三頁，編輯于2023年，星期一2023/6/3信息獲取研究方法變化過程分析機制分析實例2：土地利用變化分析更新土地利用圖社經統(tǒng)計數據動態(tài)變化研究變化信息的解譯與制圖變化數據管理原始遙感信息專題圖件數據預處理空間操作與分析圖像預處理影像分類變化信息探測結果比較變化信息分類高分辨率圖像融合土地利用/覆蓋動態(tài)變化研究框架第二十四頁，共三十三頁，編輯于2023年，星期一2023/6/3信息獲取1999、2003和2007年遙感影像圖第二十五頁，共三十三頁，編輯于2023年，星期一2023/6/3面積統(tǒng)計年份類型耕地園地林地城鎮(zhèn)村與工礦用地交通用地水域未利用地19991135.37258.32678.22192.1230.276.4722.2420031112.65266.32580.97194.31105.6286.4946.462007976.25319.772684.80197.20105.4686.3823.08第二十六頁，共三十三頁，編輯于2023年，星期一2023/6/3土地類型轉移矩陣:根據兩個不同時間（t和t+△t）的土地利用圖計算從一種類型到另一種類型的轉換概率，來分析土地利用變化過程。它依靠GIS技術，將兩個不同時間的土地利用圖進行柵格化處理，計算t時刻上A類有多少格網點轉換成在t+△t時刻的B，C，D等類型，轉換點數占該類型總數的百分比可稱為轉移概率。研究方法使用景觀生態(tài)學方法，計算土地類型轉移矩陣。

12┄nai+1a11a12┄a1na1+2a21a22┄a2na2+┆┆┆┆

Nan1a32┄annan+a+ja+1a+2┄a+na注：表中ai+代表土地類型I在轉化前的總量；a+j代表土地類型I在轉化后的總量；a是土地總面積類型轉移距陣示例第二十七頁，共三十三頁，編輯于2023年，星期一2023/6/3朱溪河小流域1999-2007土地利用類型轉移矩陣hm2土地利用轉移矩陣耕地園地林地城鎮(zhèn)村及工礦用地交通用地水域未利用地1999年合計1999年轉移量耕地624.5159.42295.474.0245.5129.376.891135.37510.86園地37.73161.5343.2811.846.691.340.12258.396.77林地245.0573.12285.3131.4129.969.210.162678.22392.91城鎮(zhèn)村及工礦用地36.0417.2754.8470.2813.280.330.05192.12121.84交通用地6.563.230.4517.68.010.42030.221.19水域22.650.025.401.242.82076.4733.65未利用地6.1744.22.231.982.157.5122.2414.732007年合計976.25319.572684.68197.2105.4686.3823.084392.92--2007年轉移量351.74158.04399.37126.9297.4543.5615.57--1191.95第二十八頁，共三十三頁，編輯于2023年，星期一2023/6/3分析目的和評價準則：1）估計住宅用地被洪水淹設而造成的損失；2）洪水水位的相對高程為500米；3）損失的大小和居民的財產、地基的穩(wěn)定性有關；應用實例3：洪水淹沒損失分析獲取數據：1）數字化的地塊多邊形地圖。每個地塊均有土地使用、可遭損失的財產狀況（簡稱估計財產）、不同地基類型等屬性。2）地塊多邊形屬性表中有地均財產這一項，地均財產=估計財產/地塊面積。3）對每一類地基，可估計其穩(wěn)定性，并估計房屋倒坍的可能性，稱損失系數（見表）4）數字化的等高線地形圖。第二十九頁，共三十三頁，編輯于2023年，星期一2023/6/3地塊多邊形圖地塊屬性表地形等高線及其