• 路径教程
    • 贝塞尔示例
    • 复合路径

    路径教程

    原文:Path Tutorial

    译者:飞龙

    协议:CC BY-NC-SA 4.0

    位于所有matplotlib.patch对象底层的对象是Path,它支持movetolinetocurveto命令的标准几个,来绘制由线段和样条组成的简单和复合轮廓。 路径由(x,y)顶点的(N,2)数组,以及路径代码的长度为 N 的数组实例化。 例如,为了绘制(0,0)(1,1)的单位矩形,我们可以使用这个代码:

    1. import matplotlib.pyplot as plt
    2. from matplotlib.path import Path
    3. import matplotlib.patches as patches
    4. verts = [
    5. (0., 0.), # left, bottom
    6. (0., 1.), # left, top
    7. (1., 1.), # right, top
    8. (1., 0.), # right, bottom
    9. (0., 0.), # ignored
    10. ]
    11. codes = [Path.MOVETO,
    12. Path.LINETO,
    13. Path.LINETO,
    14. Path.LINETO,
    15. Path.CLOSEPOLY,
    16. ]
    17. path = Path(verts, codes)
    18. fig = plt.figure()
    19. ax = fig.add_subplot(111)
    20. patch = patches.PathPatch(path, facecolor='orange', lw=2)
    21. ax.add_patch(patch)
    22. ax.set_xlim(-2,2)
    23. ax.set_ylim(-2,2)
    24. plt.show()

    路径教程 - 图1

    下面的路径代码会被接受:

    代码 顶点 描述
    STOP 1 (被忽略) 标志整个路径终点的标记(当前不需要或已忽略)
    MOVETO 1 提起笔并移动到指定顶点
    LINETO 1 从当前位置向指定顶点画线
    CURVE3 2 (一个控制点,一个终点) 从当前位置,以给定控制点向给定端点画贝塞尔曲线
    CURVE4 3 (两个控制点,一个终点) 从当前位置,以给定控制点向给定端点画三次贝塞尔曲线
    CLOSEPOLY 1 (点自身被忽略) 向当前折线的起点画线

    贝塞尔示例

    一些路径组件需要以多个顶点来指定:例如CURVE3是具有一个控制点和一个端点的贝塞尔曲线,CURVE4具有用做两个控制点和端点的三个顶点。 下面的示例显示了CURVE4贝塞尔曲线 - 贝塞尔曲线将包含在起始点,两个控制点和终点的凸包中:

    1. import matplotlib.pyplot as plt
    2. from matplotlib.path import Path
    3. import matplotlib.patches as patches
    4. verts = [
    5. (0., 0.), # P0
    6. (0.2, 1.), # P1
    7. (1., 0.8), # P2
    8. (0.8, 0.), # P3
    9. ]
    10. codes = [Path.MOVETO,
    11. Path.CURVE4,
    12. Path.CURVE4,
    13. Path.CURVE4,
    14. ]
    15. path = Path(verts, codes)
    16. fig = plt.figure()
    17. ax = fig.add_subplot(111)
    18. patch = patches.PathPatch(path, facecolor='none', lw=2)
    19. ax.add_patch(patch)
    20. xs, ys = zip(*verts)
    21. ax.plot(xs, ys, 'x--', lw=2, color='black', ms=10)
    22. ax.text(-0.05, -0.05, 'P0')
    23. ax.text(0.15, 1.05, 'P1')
    24. ax.text(1.05, 0.85, 'P2')
    25. ax.text(0.85, -0.05, 'P3')
    26. ax.set_xlim(-0.1, 1.1)
    27. ax.set_ylim(-0.1, 1.1)
    28. plt.show()

    路径教程 - 图2

    复合路径

    所有在 matplotlib,Rectangle,Circle,Polygon 等中的简单补丁原语都是用简单的路径实现的。通过使用复合路径,通常可以更有效地实现绘制函数,如hist()bar(),它们创建了许多原语,例如一堆Rectangle,通常可使用复合路径来实现。bar创建一个矩形列表,而不是一个复合路径,很大程度上出于历史原因:路径代码是比较新的,bar在它之前就存在。虽然我们现在可以改变它,但它会破坏旧的代码,所以如果你需要为了效率,在你自己的代码中这样做,例如,创建动画条形图,在这里我们将介绍如何创建复合路径,替换bar中的功能。

    我们将通过为每个直方图的条形创建一系列矩形,来创建直方图图表:矩形宽度是条形的宽度,矩形高度是该条形中的数据点数量。首先,我们将创建一些随机的正态分布数据并计算直方图。因为 numpy 返回条形边缘而不是中心,所以下面的示例中bins的长度比n的长度大 1:

    1. # histogram our data with numpy
    2. data = np.random.randn(1000)
    3. n, bins = np.histogram(data, 100)

    我们现在将提取矩形的角。 下面的每个leftbottom等数组长度为len(n),其中n是每个直方图条形的计数数组:

    1. # get the corners of the rectangles for the histogram
    2. left = np.array(bins[:-1])
    3. right = np.array(bins[1:])
    4. bottom = np.zeros(len(left))
    5. top = bottom + n

    现在我们必须构造复合路径,它由每个矩形的一系列MOVETOLINETOCLOSEPOLY组成。 对于每个矩形,我们需要 5 个顶点:一个代表MOVETO,三个代表LINETO,一个代表CLOSEPOLY。 如上表所示,closepoly的顶点被忽略,但我们仍然需要它来保持代码与顶点对齐:

    1. nverts = nrects*(1+3+1)
    2. verts = np.zeros((nverts, 2))
    3. codes = np.ones(nverts, int) * path.Path.LINETO
    4. codes[0::5] = path.Path.MOVETO
    5. codes[4::5] = path.Path.CLOSEPOLY
    6. verts[0::5,0] = left
    7. verts[0::5,1] = bottom
    8. verts[1::5,0] = left
    9. verts[1::5,1] = top
    10. verts[2::5,0] = right
    11. verts[2::5,1] = top
    12. verts[3::5,0] = right
    13. verts[3::5,1] = bottom

    剩下的就是创建路径了,将其添加到PathPatch,将其添加到我们的轴域:

    1. barpath = path.Path(verts, codes)
    2. patch = patches.PathPatch(barpath, facecolor='green',
    3. edgecolor='yellow', alpha=0.5)
    4. ax.add_patch(patch)

    结果为:

    路径教程 - 图3