arcpy从img中提取高程值

以下是使用ArcPy从IMG文件中提取高程值的示例代码：

import arcpy

# 设置工作空间和输入图层路径
arcpy.env.workspace = r"C:\data"
input_img = "elevation.img"

# 定义提取高程值的点坐标
x, y = 100, 200

# 提取高程值
elevation = arcpy.GetCellValue_management(input_img, str(x) + " " + str(y), "1").getOutput(0)
print("坐标 ({}, {}) 的高程值为: {}".format(x, y, elevation))

以上代码将从名为`elevation.img`的IMG文件中提取坐标为(100, 200)的点的高程值，并将其打印输出。`GetCellValue_management`函数用于提取高程值，输入参数为图层路径、点坐标和波段数（此处为1）。提取的高程值将作为字符串类型返回，可以通过`float()`函数将其转换为浮点数类型。

相关问题

arcgis提取各个子流域平均坡度

### 计算各子流域的平均坡度

为了在ArcGIS中计算每个子流域的平均坡度，可以采用以下方法：

#### 准备工作
确保拥有数字高程模型（DEM）数据以及定义好的子流域边界矢量图层。这些数据对于后续操作至关重要。

#### 创建坡度栅格
使用`Spatial Analyst Tools`中的`slope`工具处理DEM文件，生成表示地形表面倾斜程度的地图。此过程会创建一个新的栅格数据集，其中每个像元代表该位置处的地表斜率[^2]。

import arcpy
from arcpy.sa import *

arcpy.CheckOutExtension("Spatial")

# 设置环境参数
arcpy.env.workspace = "C:/data"
demRaster = Raster("elevation.tif")  # 输入 DEM 数据路径
slopeOutput = Slope(demRaster, "DEGREE", 1)

# 输出坡度结果至指定位置
slopeOutput.save("c:/output/slope_degrees.img")

#### 统计子流域内坡度均值
借助Zonal Statistics as Table功能，能够针对每一个单独的子流域区域聚合内部所有像素点上的坡度数值并求取平均数。这一步骤的关键在于设置好zone layer即为之前准备好的子流域范围shp文件；value raster则指向刚刚制作完成的slope raster图像。

inZoneData = r"C:\data\watershed.shp"    # 子流域边界 shapefile 路径
zoneField = "FID"                         # zone 层字段名
inValueRaster = slopeOutput               # 坡度影像作为 value raster 参数传入
outTable = r"c:\output\slope_stats.dbf"

# 执行 zonal statistics 并保存到表格
zStat = ZonalStatisticsAsTable(inZoneData, zoneField, inValueRaster,
                               outTable, "NODATA", "MEAN")

上述脚本实现了自动化批量获取各个子流域对应的平均坡度信息，并将其记录在一个dbase格式的结果集中以便进一步分析应用。