ArcGIS地形起伏度提取及出图设置

地形起伏度

1.1介绍

地形起伏度是指特定区域内,最高点海拔与最低点海拔的差值,它是描述一个区域地形特征的宏观性指标。利用领域分析中的焦点统计来计算地形起伏度,其中统计类型选择RANGE— 计算邻域内像元的范围(最大值和最小值之差)。得到地形起伏的结果图见图6-4。地形起伏度可表示为

RFi=Hmax-Hmin

式中,RFi为地形起伏度;Hmax分析窗口内的最大高程值;Hmin为分析窗口内的最小闻程值。地形的起伏是导致水土流失的最直接因素地形起伏度能够直观地反映地形起伏特征。

1.2 DEM的提取

 

打开全国县级行政规划地图

 

打开其属性表,在其中按属性选择(双流区代码为510122)双流区

 

在选中的情况下,打开toolbox-分析工具-提取分析-裁剪

 

其中输入要素和输出要素都是刚才已经选中的双流区的图层,其他默认

 

得到的裁剪后的图层如图所示

### 使用ArcGIS计算和显示地形起伏 #### 创建地理数据库并准备环境 为了便于管理和存储处理后的数据,在特定文件夹中新建一个地理数据库用于存放结果。确保该地理数据库的位置方便访问,并且有足够的空间来保存即将生成的数据集。 #### 准备DEM数据 确保已经拥有高质量的数字高程模型(Digital Elevation Model, DEM)数据,这是后续所有分析的基础。对于重庆市的研究区而言,应获取覆盖整个研究区域的DEM数据[^2]。 #### 利用邻域统计工具计算地形起伏 通过ArcToolbox中的`Spatial Analyst Tools -> Neighborhood - **输入栅格层**:选择之前准备好的DEM数据作为输入源。 - **邻域形状**:定义要应用于每个像元周围区域的形状;通常可以选择圆形或矩形窗口。 - **统计数据类型**:选择`Range`选项,这会自动计算指定邻域内的最大值与最小值之间的差异,即代表了局部地区的地形起伏[^3]。 执行上述命令后即可获得反映不同位置之间高变化情况的地——这就是所谓的“地形起伏”。 #### 结果可视化 完成计算之后,可以进一步调整颜色方案和其他符号化属性以便更好地展示所得的信息。例如,采用渐变色彩映射方法能够直观地区分高低起伏明显的地带以及相对平坦之处。 ```python import arcpy from arcpy.sa import * arcpy.CheckOutExtension("Spatial") # 设置工作空间 workspace = r"C:\path\to\your\geodatabase.gdb" arcpy.env.workspace = workspace # 输入DEM数据路径 dem_raster = Raster(r"path_to_dem.tif") # 定义邻域大小 (此处假设为3*3矩阵) neighborhood = NbrRectangle(3, 3) # 应用FocalStatistics函数计算RANGE terrain_roughness = FocalStatistics(dem_raster, neighborhood, "RANGE") # 输结果至当前工作空间下的新栅格文件 output_path = os.path.join(workspace, "TerrainRoughness") terrain_roughness.save(output_path) ```
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值