【ArcGIS|空间分析】利用DEM提取山脊、山谷线（水文分析、表面分析）
参考汤国安、杨昕等编著的《地理信息系统空间分析实验教程第二版》
利用水文分析（流量）和表面分析（坡向变率）得到山脊和山谷线
文章目录
- 1、基本概念
- 正地形
- 负地形
- 正负地形的影响
- 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 操作步骤
1、基本概念
正地形
是相对高于邻区或新构造上升地区的地形。山地、高原、丘陵都是正地形。世界上许多石油矿床与高地、丘陵、特别是穹隆或背斜构造形成的正地形有关。天然气、金属矿床常常形成或储存在上述正地形中。
负地形
是相对低于邻区或新构造下沉地区的地形。洼地、盆地都是负地形。负地形是沉积物堆积的有利条件，也是冲刷微弱的场所。煤、石油、铝土、铁、泥炭、盐类和锰结核等沉积矿床多形成在盆地、凹地、平原和洋盆等负地形中。
正负地形的影响
负地形的地貌格局决定了一些地区不断地接纳、汇聚、堆积，在外力作用下，由四周正地形区搬运迁移的陆相碎屑物。“夷平填洼”是不可改变的自然地理过程中外营力作用的运动规律，只要这种正负地形差异格局的存在，就将持续地进行下去，直至夷平或填平。
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。
【ArcGIS|空间分析】利用DEM提取山脊、山谷线（水文分析、表面分析）相关推荐
- matlab 山脊 提取,ArcGIS中利用水文分析提取山脊线山谷线
1 流程图 利用水文分析提取山脊线及山谷线,山脊线相当于分水线,山谷线相当于山谷线.分水线是水流的起源点,这些栅格的水流方向只存在流出方向而不存在流入方向,所以汇流累积量为零.通过对零值的提取就可以得 ...
- 利用水文分析提取山脊线山谷线
1 流程图 利用水文分析提取山脊线及山谷线,山脊线相当于分水线,山谷线相当于山谷线.分水线是水流的起源点,这些栅格的水流方向只存在流出方向而不存在流入方向,所以汇流累积量为零.通过对零值的提取就可以得 ...
- 基于matlab山脊线,教你如何利用水文，分析提取山脊线山谷线
原标题:教你如何利用水文,分析提取山脊线山谷线 利用水文分析提取山脊线及山谷线,山脊线相当于分水线,山谷线相当于山谷线.分水线是水流的起源点,这些栅格的水流方向只存在流出方向而不存在流入方向,所以汇流 ...
- ArcGIS实战：利用LandSat8 提取水体并计算面积
ArcGIS实战:利用LandSat8 提取水体并计算面积 1.数据下载 LandSat 8 遥感影像属于多波段遥感影像,利用不同波段可以做一些分析,比如NDWI(水体提取),NDVI(绿地提取)等. ...
- 生信分析-利用TBtools提取序列
准备文件 gtf/gff3 全基因组fasta 背景知识 CDS=coding sequence,是编码区,是可翻译成蛋白质的exon的集合 cDNA比CDS多了5'-UTR和3'-UTR区域,是所有 ...
- 【ArcGIS|空间分析】数字地形分析
文章目录数字地形分析原理基于ArcGIS的数字地形分析操作DEM的建立1 栅格表面的创建2 TIN的创建3 等高线的创建4 Terrain(地形数据集)的建立基本因子分析1 坡度(栅格表面与TIN表面 ...
- 01 ArcGIS空间分析--水文分析--利用水文分析方法提取山谷山脊线
ArcGIS空间分析的一个特别好的分析教程是汤国安老师的空间分析的内容. 实例与练习 练习1. 利用水文分析方法提取山脊.山谷线 1． 背景:山脊线.山谷线是地形特征线,它们对地形.地貌具有一定的控制 ...
- [空间分析] DEM提取山脊线与山谷线
文章目录 基于图像处理技术的原理 基于地形表面几何形态分析原理的算法 基于地形表面流水物理模拟分析原理的算法 基于地形表面几何形态分析和流水物理模拟分析相结合 平面曲率与坡位组合法 水文分析法 [山脊 ...
- ArcGIS山脊线、山谷线和山顶点的提取(附练习数据下载)
特征地形要素,主要是指对地形在地表的空间分布特征具有控制作用的点.线或面状要素.特征地形要素构成地表地形与起伏变化的基本框架.与地形指标的提取主要采用小范围的邻域分析不同的是,特征地形要素的提取更多地 ...
- Arcgis中山脊线，山谷线的提取，以及流域的分割
目录 实验数据下载:https://pan.baidu.com/s/1jfTv5LXBDaabguJsPTFJEw 提取码:4g1b Task 1:利用Ex1中的数据,练习提取不同位置的地形剖面线: ...
最新文章
- linux编译安装jpeg,Linux下JPEG库安装脚本(转)
- Leangoo如何颠覆传统项目管理软件？
- ycsb 测试验证模式的mongodb
- python字典定义方式_Python字典常用方法及汇总
- Ubuntu SSH Algorithm negotiation failed
- 国产中标麒麟Linux部署dotnet core 环境并运行项目 （三） 部署运行WEB API项目
- SSM+Maven+Dubbo+Zookeeper简单项目实战以及易错注意点
- 更新条目时出错。有关详细信息，请参见内部异常。
- rfp计算机,RFP（中英文）.doc
- docker es持久化_docker如何保持数据持久化？其实可以考虑通过数据卷挂载的方式~...
- HCIA-IoT V2.5物联网初级考试总结(附题库,历年真题,刷题软件)
- PHP网站从Apache转移到Nginx后产生404错误的原因和解决办法
- 小萝卜控机大师录制脚本（手机app自动化）
- 省级、县级行政区shapefile下载
- html表单 新页面打开新页面,js提交表单并打开新页面的方法
- java下载图片到本地
- GEE-Python遥感大数据分析、管理与可视化实践技术应用
- 小象学院—小象学院—Python学习联盟09—空气质量指数计算
- 继上篇，关于Rxjava，让你知道什么叫“大海无量”
- Android Settings 应用二 获取应用消耗的流量