首先，使用FLow Direction得到自然流向的raster。

然后，在ArcMap中打开，用spatial Analyst模块中的Raster to Feature，将需要区域的Raster转换为Point对象，

为了对应流向，先向Point层增加一条新的字段，用如下脚本进行计算：

dim a as double, o as double
a = log([GRID_CODE])/log(2.0)
if a - Int(a) = 0 then
  o =(a/8*360+180) mod 360
else
  o = -1
end if

然后使用Unique Values方式显示Point层，并将－1值的数据显示为一个圆点，将其他数据显示为一个向上的箭头，并用Advanced选项中的Rotation进行字段旋转。选择ArithMetric方式。

这样，就得到了流向图。