由于面状水系可能存在多条中间线，因此批量提取时需要使用 ArcGIS 中的 Feature To Line 工具结合 Python 循环和游标来完成。
以下是代码：
import arcpy
import os
# 设置输入输出路径和文件名
input_folder = r"C:\data\river_polygons"
output_folder = r"C:\data\river_midlines"
# 创建输出文件夹
if not os.path.exists(output_folder):
os.makedirs(output_folder)
# 遍历输入文件夹中所有 Shapefile 文件
for _, _, files in os.walk(input_folder):
for file in files:
if file.lower().endswith(".shp"):
# 获取当前文件路径和名称
input_shapefile = os.path.join(input_folder, file)
# 在提取线之前，删除现有的缓存文件 (如果已存在)
arcpy.Delete_management("in_memory")
# 根据面状水系区域生成几何网络
arcpy.Create geometric network（流体模拟）
network_path = "in_memory/network"
arcpy.FeatureClassToFeatureClass_conversion(input_shapefile, "in_memory", "polygons_for_network")
arcpy.CreateGeometricNetwork_management(network_path, "River", "SIMPLE_JUNCTION", "#", "#", "#", "#")
arcpy.AddFeatureClassToGeometricNetwork_management(network_path, "River", "polygons_for_network", "#", "#")
arcpy.BuildGeometricNetwork_management(network_path)
# 取得面状水系中心线
centerlines = []
with arcpy.da.SearchCursor(input_shapefile, ["OID@", "SHAPE@"]) as cursor:
for row in cursor:
# 使用 Feature To Line 工具计算中心线
input_polyline_feature = "in_memory/polyline_feature"
arcpy.FeatureToLine_management(row[1], input_polyline_feature, "#", "NO_ATTRIBUTES")
arcpy.MergeDividedRoads_cartography(input_polyline_feature, "#", input_polyline_feature + "_merged")
# 向列表追加中心线
centerlines.append((row[0], input_polyline_feature + "_merged"))
del cursor
# 将所有中心线融合成一个图层并保存到输出文件夹中
output_shapefile = os.path.join(output_folder, "midlines_" + file)
arcpy.Merge_management([i[1] for i in centerlines], output_shapefile)
with arcpy.da.UpdateCursor(output_shapefile, ["PolyID"]) as cursor:
for row in cursor:
# 更新新的线性要素 ID 与原始面状水系多边形相关
row[0] = [i[0] for i in centerlines if i[1].endswith("_" + str(row[0]))][0]
cursor.updateRow(row)
del cursor
上述代码假设输入为包含多个面状水系区