日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 >

ArcGIS图解建模工具与DEM提取山脊、山谷线(水文分析、表面分析)

發布時間:2023/12/20 60 豆豆
生活随笔 收集整理的這篇文章主要介紹了 ArcGIS图解建模工具与DEM提取山脊、山谷线(水文分析、表面分析) 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

參考湯國安、楊昕等編著的《地理信息系統空間分析實驗教程第二版》 利用水文分析(流量)和表面分析(坡向變率)得到山脊和山谷線

文章目錄

  • 1、基本概念
  • 1.1 正地形
  • 1.2 負地形
  • 正負地形的影響
  • 2、水文分析方法
  • 2.1 基本思想原理
  • 2.2 技術路線
  • 2.3 操作步驟
  • 2.3.1 正負地形的提取
  • 2.3.2 利用水文分析思想和工具提取研究區域的山脊線
  • 2.3.3 利用水文分析思想和工具提取研究區域的山谷線
  • 2.3.4 利用水文分析提取山脊線和山谷線結果
  • 3、表面分析方法
  • 3.1 基本思想原理
  • 3.2 技術路線
  • 3.3 操作步驟
  • 4、圖解建模方法與思路
  • 4.1 建立概念模型
  • 4.2 構建圖解模型
  • 4.3 模型的使用
  • 4.4 運行結果
  • 4.5 模型文件與參考圖書pdf下載
  • 5、方法評價
  • 5、文章參考與致謝

1、基本概念

1.1 正地形

是相對高于鄰區或新構造上升地區的地形。山地、高原、丘陵都是正地形。世界上許多石油礦床與高地、丘陵、特別是穹隆或背斜構造形成的正地形有關。天然氣、金屬礦床常常形成或儲存在上述正地形中。

1.2 負地形

是相對低于鄰區或新構造下沉地區的地形。洼地、盆地都是負地形。負地形是沉積物堆積的有利條件,也是沖刷微弱的場所。煤、石油、鋁土、鐵、泥炭、鹽類和錳結核等沉積礦床多形成在盆地、凹地、平原和洋盆等負地形中。

正負地形的影響

負地形的地貌格局決定了一些地區不斷地接納、匯聚、堆積,在外力作用下,由四周正地形區搬運遷移的陸相碎屑物。“夷平填洼”是不可改變的自然地理過程中外營力作用的運動規律,只要這種正負地形差異格局的存在,就將持續地進行下去,直至夷平或填平。

2、水文分析方法

2.1 基本思想原理

通過對零值的匯流累積值的柵格的提取就可以得到分水線,也就得到了山脊線。對于山谷線可以利用反地形的特點,可以利用一個較大的數值減去原始的DEM數據得到與原始地形完全相反的地形數據,使得原始的DEM中的山脊變成反地形的山谷,而原始DEM中的山谷在反地形中就變成了山脊。山谷線的提取就可以在反地形中利用提取山脊線的方法進行提取。但是此方法提取出的山脊和山谷位置有些偏差,可以利用正、負地形加以糾正。

2.2 技術路線

2.3 操作步驟

2.3.1 正負地形的提取

1、 在ArcMap中加載樣區的原始DEM數據.
2、 選擇【Spatial Analyst工具】|【鄰域分析】|【焦點統計】,以11×11的窗口計算平均值。
3、 選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,將原始DEM減去平均值DEM,獲得正負地形數據。
4、 選擇【Spatial Analyst工具】|【重分類】|【重分類】,分別進行兩次重分類。一次將大于0的區域賦值為1(正地形),小于0的區域賦值為0,另一次將小于0的區域賦值為1(負地形),大于0的區域賦值為0。

2.3.2 利用水文分析思想和工具提取研究區域的山脊線

1、 在ArcMap中加載樣區的原始DEM數據。
2、 洼地填充:選擇【Spatial Analyst工具】|【水文分析】|【填洼】,輸入表面柵格數據。選擇將所有洼地全部填充,所以不需要計算洼地深度,Z limit為默認值。

3、 基于無洼地的水流方向計算:選擇【Spatial Analyst工具】|【水文分析】|【流向】,輸入填洼后的表面柵格數據。

4、 匯流累積量的計算:選擇【Spatial Analyst工具】|【水文分析】|【流量】,輸入水流方向數據。

5、 匯流累積量為零值的提取:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,輸入計算公式為facc0=(flowaa1==0),facc0為輸出數據,flowaa1為匯流累積量數據。

6、 數據光滑處理:上一步提取數據中很多不是山脊線,選擇【Spatial Analyst工具】|【鄰域分析】|【焦點統計】,以3×3窗口做平均值計算。

7、 輔助判斷山脊線位置:選擇【Spatial Analyst工具】|【表面分析】|【等值線】和【Spatial Analyst工具】|【表面分析】|【山體陰影】。

8、 選擇【Spatial Analyst工具】|【重分類】|【重分類】,以等高線圖和山體陰影圖判斷山脊線位置,調整分級臨界點,最終確定分界閾值為0.5541。將光滑處理后的數據二值化,接近1的部分賦值為1,其余部分賦值為0。

9、 消除存在于負地形區域中的錯誤山脊線:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,將重分類后的數據與正地形數據相乘。

10、 獲得最終山脊線:選擇【Spatial Analyst工具】|【重分類】|【重分類】,將屬性值不為1的賦值為NoData。

2.3.3 利用水文分析思想和工具提取研究區域的山谷線

1、 在ArcMap中加載樣區的原始DEM數據。
2、 獲取反地形:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,計算公式為fandem=Abs(dem-2000),fandem為輸出數據,dem為原始DEM數據。

3、 水流方向計算(無需填洼):選擇【Spatial Analyst工具】|【水文分析】|【流向】,輸入填洼后的表面柵格數據。

4、 匯流累積量的計算:選擇【Spatial Analyst工具】|【水文分析】|【流量】,輸入水流方向數據。

5、 匯流累積量為零值的提取:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,輸入計算公式為facc0=(flowaa2==0),facc0為輸出數據,flowaa2為匯流累積量數據。

6、 數據光滑處理:上一步提取數據中很多不是山谷線(反地形中山脊線為原地形中的山谷線),選擇【Spatial Analyst工具】|【鄰域分析】|【焦點統計】,以3×3窗口做平均值計算。

7、 輔助判斷山谷線位置:選擇【Spatial Analyst工具】|【表面分析】|【等值線】和【Spatial Analyst工具】|【表面分析】|【山體陰影】。

8、 選擇【Spatial Analyst工具】|【重分類】|【重分類】,以等高線圖和山體陰影圖判斷山脊線位置,調整分級臨界點,最終確定分界閾值為0.65667。將光滑處理后的數據二值化,接近1的部分賦值為1,其余部分賦值為0。

9、 消除存在于負地形區域中的錯誤山谷線:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,將重分類后的數據與負地形數據(不是反地形數據)相乘。

10、獲得最終山谷線:選擇【Spatial Analyst工具】|【重分類】|【重分類】,將屬性值不為1的賦值為NoData。

2.3.4 利用水文分析提取山脊線和山谷線結果

3、表面分析方法

3.1 基本思想原理

首先利用DEM數據提取地面的平面曲率及地表的正負地形,取正地形上平面曲率的大值即為山脊,負地形上平面曲率的大值為山谷。在實際應用中,由于平面曲率的提取比較繁瑣,而坡向變率(SOA)在一定程度上可以很好地表征平面曲率。

3.2 技術路線

3.3 操作步驟

1、 求DEM坡向變率數據:選擇【Spatial Analyst工具】|【表面分析】|【坡向】,將得到的坡向數據再求一次坡度得到坡向變率數據,記為SOA1。
2、 求DEM最大值數據:選擇【Spatial Analyst工具】|【鄰域分析】|【焦點統計】,以11×11窗口做最大值計算得到DEM最大高程值。
3、 求反地形數據:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,輸入公式為H-DEM(H為DEM最大高程值數據,DEM為原始DEM數據)得到反地形DEM數據。
4、 求反地形坡向變率數據:選擇【Spatial Analyst工具】|【表面分析】|【坡向】,基于反地形DEM數據求反地形DEM坡向變率數據,記為SOA2。
5、 求沒有誤差的DEM坡向變率:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,輸入公式為SOA=(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/ 2 ,得到沒有誤差的DEM坡向變率數據。
6、 求DEM平均值:選擇【Spatial Analyst工具】|【鄰域分析】|【焦點統計】,以11×11窗口做平均值計算得到DEM最平均值。
7、 獲得正負地形數據:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,計算公式為C=DEM-B(C為輸出數據,DEM為原始DEM數據,B為DEM平均值數據)得到正負地形數據。
8、 提取山脊線:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,計算公式為shanji=[C]>0&SOA>70。
9、 提取山谷線:選擇【Spatial Analyst工具】|【地圖代數】|【柵格計算器】,計算公式為shanji=[C]<0&SOA>70。

4、圖解建模方法與思路

4.1 建立概念模型

山脊線是處于正地形上且水流匯集量為0的點的集合;
山谷線是處于負地形上且在DEM反地形上水流匯聚量為0。

4.2 構建圖解模型

4.3 模型的使用

4.4 運行結果

部分山谷線 部分山脊線

4.5 模型文件與參考圖書pdf下載

鏈接:https://pan.baidu.com/s/1-IKBNU7zZeFhJoibaZq7Vw 提取碼:n4lg

5、方法評價

【基本思想是】

按照流水從高至低的自然規律,順序計算每一柵格點上的匯水量,然后按匯水量單調增加的順序,由高到低找出區域中的每一條匯水線。
根據得到的匯水線,通過計算找出各自匯水區域的邊界線,就得到了分水線

【評價】算法采用了DEM的整體追蹤分析的思路與方法,分析結果具有系統性好,還便于進行相應的徑流成因分析
【缺陷】該方法也存在以下兩個明顯的缺陷

  • 由于該算法所計算的匯水量與高程有關,計算的結果必然是高程值大的地形特征線上的點的匯水量小,高程值小的地形特征線上的點的匯水量大。
  • 由于該算法降格匯水區域的公共邊界視為分水線,因此它所確定的分水線均為閉合曲線,這與實際的地形特征線(山脊線)不符
  • 5、文章參考與致謝

    文章參考

    https://blog.csdn.net/qq_40628258/article/details/83473964#_5
    https://blog.csdn.net/summer_dew/article/details/82967666

    感謝參考文章與當時空間分析授課的許老師,也感謝自己的耕耘!^- ^

    總結

    以上是生活随笔為你收集整理的ArcGIS图解建模工具与DEM提取山脊、山谷线(水文分析、表面分析)的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。