arcgis提取河网的阈值计算

基于DEM提取河网中的阈值研究

[摘要]本研究以大兴安岭地区的塔河流域为研究对象，基于DEM数据，利用ArcGIS10软件自动提取水系信息。探讨不同阈值对提取结果的影响，通过实地考察数据和所提取的河网矢量数据对照，最终确定一个最佳的阈值，完成塔河流域的河网提取。从提取的实验结果来看，表明用该方法提取研究流域特征是可行并有效的。

[关键词]DEM ArcGIS10 水系信息 提取

[中图分类号] S276.3 [文献码] B [文章编号] 1000-405X（2014）-7-208-2 0引言

数字流域水系提取系统是数字水文研究中的重要组成部分，高质量的水系提取模型能够大大提高整个水文模型在运作过程中的有效性和准确性，从而为提高区域水文模拟的精度提供了重要的保证[1]。本文运用ArcGIS10的水文分析工具，利用大兴安岭地区的DEM数据，对新林县境内的塔河流域进行河网的提取，讨论汇流面积阈值设置的大小与生成河网总长度存在的函数关系。 1研究区域概况与研究方法 1.1研究区域概况

大兴安岭境内的呼玛河系黑龙江一级支流，塔河是呼玛河的最大支流。呼玛河发源于大兴安岭主峰伊勒呼里山北麓，流向大致成“Ω

基于DEM的ArcGIS水文分析

—河网和流域的提取

一、实验背景

水文分析是DEM 数据应用的一个重要方面。而利用DEM生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型研究与地表水流有关的各种自然现象例如洪水水位及泛滥情况，划定受污染源影响的地区，预测当某一地区的地貌改变时对整个地区将造成的影响等。

二、实验目的

通过本实验，使读者理解基于DEM数据进行水文分析的基本原理， 掌握利用ArcGIS提供的水文分析工具进行水文分析的基本方法和步骤，并利用DEM数据提取出河网及流域。

三、实验数据

某地区栅格数据DEM，数据来源于随书光盘（…\\Chp9\\Ex2）。

四、实验要求

根据DEM利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。

五、实验流程图

六、实验内容及步骤

1．无洼地DEM生成

DEM 是比较光滑的地形表面模型，但由于DEM 误差以及一些真实地形或特殊地形的影响，使得DEM 表面存在一些凹陷的区域。

在进行水流方向计算时，由于这些区域的存在，往往得到不合理的甚至错误的水流方向。 因此，在进行水流方向的计算之前，

arcgis中提取线交点方法

提示：在做下列处理时先做好数据备份！！！

一下操作针对shapefile文件操作，不是shapefile的先转换格式。 1、预处理数据

在arcmap中编辑shapefile文件，工具全部选中要素，点击map topology，建立地图拓扑关系（1），然后点击点击打断线工具，在交点处打断线（2），保存结果，我这里保存到了test.shp。

1

2

2、在arccatalog中右键点击刚才保存的文件，选新建网络数据集，一路next和yes，别管内容。

2.jpg

这样就生成了结合点文件test_ND_Junctions.shp，但这还不是我们要的线交点文件，因为里

边有悬挂结点和自相交结点。

4.jpg

3、用空间连接，将test_ND_Junctions.shp和test.shp连接，生成文件

test_Junctions_SpatialJoin.shp。打开此文件属性表，发现多了一个字段join_count:

5.jpg

这个字段的数值就是点和多少条线相连的数。

4、编辑test_Junctions_SpatialJoin.shp，按属性选择join_count小于2的点，删除，剩下的点就是我们要的线的交

arcgis中提取线交点方法

提示：在做下列处理时先做好数据备份！！！

一下操作针对shapefile文件操作，不是shapefile的先转换格式。 1、预处理数据

在arcmap中编辑shapefile文件，工具全部选中要素，点击map topology，建立地图拓扑关系（1），然后点击点击打断线工具，在交点处打断线（2），保存结果，我这里保存到了test.shp。

1

2

2、在arccatalog中右键点击刚才保存的文件，选新建网络数据集，一路next和yes，别管内容。

2.jpg

这样就生成了结合点文件test_ND_Junctions.shp，但这还不是我们要的线交点文件，因为里

边有悬挂结点和自相交结点。

4.jpg

3、用空间连接，将test_ND_Junctions.shp和test.shp连接，生成文件

test_Junctions_SpatialJoin.shp。打开此文件属性表，发现多了一个字段join_count:

5.jpg

这个字段的数值就是点和多少条线相连的数。

4、编辑test_Junctions_SpatialJoin.shp，按属性选择join_count小于2的点，删除，剩下的点就是我们要的线的交

斜坡单元

地质灾害危险性区划中常用的单元类型有网格单元、地域单元、均一条件单元、子流域单元、斜坡单元等。其中: 网格单元形状较规则，便于实现快速剖分，离散后得到的矩阵形式的数据有利于进一步运算，但是不能完全反映地势起伏，与地质环境条件联系不够紧密; 均一条件单元没有考虑不同区域的地质环境条件差异; 子流域单元适用于泥石流灾害危险性区划，对滑坡、崩塌等则不适用。斜坡单元是滑坡、崩塌等地质灾害发育的基本单元，并且在各类控制或影响因素中，河流和沟谷的发育阶段对滑坡、崩塌的形成具有明显的控制作用，因此采用基于幼年期沟谷划分的斜坡单元作为评价单元，可以与地质环境条件紧密联系，综合体现各类控制或影响因素的作用，使评价结果更贴近于实际。因此，在满足DEM 精度要求的前提下，斜坡单元划分较适用于地质灾害危险性区划【1】。

斜坡单元划分原理

斜坡单元划分的实质是基于DEM 的地表水文分析，包括正反地形无洼地DEM 的生成、水流方向的提取、汇流累积量的计算、河网的生成、集水流域的生成等关键步骤，其基本原理是利用正反地形分别提取山谷线和山脊线( 分别对应于汇水线和分水线)，把生成的集水流域与反向集水流域融合，再经后期处理人工修改不合理的单元，最终得到的由汇水线与分

斜坡单元

地质灾害危险性区划中常用的单元类型有网格单元、地域单元、均一条件单元、子流域单元、斜坡单元等。其中: 网格单元形状较规则，便于实现快速剖分，离散后得到的矩阵形式的数据有利于进一步运算，但是不能完全反映地势起伏，与地质环境条件联系不够紧密; 均一条件单元没有考虑不同区域的地质环境条件差异; 子流域单元适用于泥石流灾害危险性区划，对滑坡、崩塌等则不适用。斜坡单元是滑坡、崩塌等地质灾害发育的基本单元，并且在各类控制或影响因素中，河流和沟谷的发育阶段对滑坡、崩塌的形成具有明显的控制作用，因此采用基于幼年期沟谷划分的斜坡单元作为评价单元，可以与地质环境条件紧密联系，综合体现各类控制或影响因素的作用，使评价结果更贴近于实际。因此，在满足DEM 精度要求的前提下，斜坡单元划分较适用于地质灾害危险性区划【1】。

斜坡单元划分原理

斜坡单元划分的实质是基于DEM 的地表水文分析，包括正反地形无洼地DEM 的生成、水流方向的提取、汇流累积量的计算、河网的生成、集水流域的生成等关键步骤，其基本原理是利用正反地形分别提取山谷线和山脊线( 分别对应于汇水线和分水线)，把生成的集水流域与反向集水流域融合，再经后期处理人工修改不合理的单元，最终得到的由汇水线与分

envi和arcgis结合提取图像上点的坐标信息

2010-06-24 15:36:00| 分类： 默认分类 | 标签： |字号大中小 订阅 有时候我们需要知道一些点的坐标信息，凡是有坐标信息的图像在用专业的软件显示的时候，都能看到每个点的坐标信息，但是我们如果

将我们感兴趣的点的坐标提取出来呢？

第一个你想到的可能是编程，的确，这样可以提取每个点的信息，但是对于哪些不懂得编程的人员来说，这太难了。因此我们必须寻找其

他的办法来解决这个问题。

第二个方法就是充分利用当前的软件，应该说目前的软件已经能够满足大部分的需求，当然，如果进行批处理的话，用软件进行操作就有点太繁琐了，工作量超大。如果数据量比较少的话，就很方便了。毕

竟，操作的东西还是比较容易学的。

这里我们介绍一种采用envi和arcgis结合的方式来给从遥感图像上提取的点增加坐标信息。当然这个图像要有投影坐标信息才可以 注：当将文件利用envi输出为ascii时，输出的信息是有坐标的,但

是必须在输出格式时选择x y z format 1.从envi中确定感兴趣的点或者区域

打开图像－roi tool－option－band throshold to roi

在选择阈值的时候，首先要确

２０１０．４

理论酮究

遥摩信息

基于阈值分割的影像建筑物角点提取算法优化研究

孙一权①～，冯仲科①－⑦

（①北京林业大学测绘与３Ｓ技术研究中心，北京１０００８３；

②北京林业大学省部共建森林培育与保护教育部重点实验室，北京１０００８３）

摘要：对Ｍｏｒａｖｅｃ和ＳＵＳＡＮ这两种经典的特征提取算法进行了深入研究，以Ｍａｔｌａｂ软件为实验平台，验证了它们的可行性。通过对实验结果的比较，分析了Ｍｏｒａｖｅｃ算子中运用“抑制局部非最大”思想的优势、局限性和适用性；在ＳＵＳＡＮ算子中通过考虑提取过程中的误差影响因素，为算法中几何阈值和灰度阈值的确定和模板选择提供了依据，对其影响范围进行了定量分析。指出在两种算法实现过程中，Ｍｏｒａｖｅｃ算子由于对如何确定判断窗口的大小不能给出可靠的依据，加之缺少对影像的先验知识，无法了解研究区域角点分布状况。因此，仅仅用这种方法很难获得完全准确的角点特征；对于ＳＵＳＡＮ算子，根据分析结果对几何阈值加设一个下限ｄ，在很大程

度上排除了孤立噪声点的干扰；并运用“重心距离法”，即通过判定ＳＵＳＡＮ重心与核心点连线上的像素点的边缘初始值的相近条件，消除了虚假的角点，从而实现了对算法的优化。

关键词ｉ特征提取；兴趣算子；角点检测；抑制局部非最大

地理国情监测云平台

ArcGIS教程：栅格像元的提取

一、按属性提取

摘要

基于逻辑查询提取栅格像元。 插图

用法

· 如果 Where 子句的求值结果是 True，则将为该像元位置返回初始输入值。如果其求值结果是 False，则将为像元位置指定 NoData。

· 对于输出栅格，将舍弃输入栅格中的所有额外项(值和计数除外)。

· 如果在 Where 子句中指定了除输入栅格的值以外的某一项，则将为此像元位置返回初始输入值。 · 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。

默认输出格式是 Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9个字符。

· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层，提取操作将只考虑图层加载(符号化)的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。 · 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。 语法

ExtractByAttributes (in_raster, w

地理国情监测云平台

ArcGIS教程：栅格像元的提取

一、按属性提取

摘要

基于逻辑查询提取栅格像元。 插图

用法

· 如果 Where 子句的求值结果是 True，则将为该像元位置返回初始输入值。如果其求值结果是 False，则将为像元位置指定 NoData。

· 对于输出栅格，将舍弃输入栅格中的所有额外项(值和计数除外)。

· 如果在 Where 子句中指定了除输入栅格的值以外的某一项，则将为此像元位置返回初始输入值。 · 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。

默认输出格式是 Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9个字符。

· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层，提取操作将只考虑图层加载(符号化)的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。 · 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。 语法

ExtractByAttributes (in_raster, w