Arcgis提取河道纵剖面

1. 提取出Stream

from arcpy import env  
from arcpy.sa import *  
env.workspace = PATH  
pj_fill = Fill("pj_dem")  
pj_flowdir = FlowDirection("pj_fill")  
pj_flowacc = FlowAccumulation("pj_flowdir")  
pj_stream = Raster("pj_flowacc") > 3000  
pj_sorder = StreamOrder("pj_stream","pj_flowdir")  
pj_slink = StreamLink("pj_stream","pj_flowdir")  
pj_basin = Basin("pj_flow")  
pj_net = SnapPourPoint("pj_slink", "pj_flowacc", 0, "VALUE")  
pj_watershed = Watershed("pj_flowdir","pj_net")  
pj_river = StreamToFeature("pj_stream", "pj_flowdir", env.workspace, "NO_SIMPLIFY")

2. 插值为polyline z

3D Analyst Tools→Functional Surface→Interpolate Shape

（如果3D Analyst Tools无法打开，需要在Customize→Extensions中加载）

Interpolate Shape工具如下：Input Surface选择DEM，Input Feature Class选择polyline（为上面的river）。Output Feature Class就是输出路径。Method可选Linear或BILINEAR。

由此可以提取出Polyline z，即空间多段线。

如果插值不成功也可以尝试Adding surface information。可能下一步也没法做，但是可以直接用Attribute Table自己导出数据画图。

3. 提取剖面

加载3D Analyst工具条：Customize→Toolbars→3D Analyst

选中之前的Polyline Z，然后选择上面的红框中的工具，弹出剖面。

原文:大专栏  Arcgis提取河道纵剖面

原文地址：https://www.cnblogs.com/petewell/p/11585196.html