def curvelength(x0, y0, x1, y1, x2, y2, x3, y3, n=20): """Returns the length of the spline. Integrates the estimated length of the cubic bezier spline defined by x0, y0, ... x3, y3, by adding the lengths of lineair lines between points at t. The number of points is defined by n (n=10 would add the lengths of lines between 0.0 and 0.1, between 0.1 and 0.2, and so on). The default n=20 is fine for most cases, usually resulting in a deviation of less than 0.01. """ length = 0 xi = x0 yi = y0 for i in range(n): t = 1.0 * (i+1) / n pt_x, pt_y, pt_c1x, pt_c1y, pt_c2x, pt_c2y = \ curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3) c = sqrt(pow(abs(xi-pt_x),2) + pow(abs(yi-pt_y),2)) length += c xi = pt_x yi = pt_y return length
def curvelength(x0, y0, x1, y1, x2, y2, x3, y3, n=20): """Returns the length of the spline. Integrates the estimated length of the cubic bezier spline defined by x0, y0, ... x3, y3, by adding the lengths of lineair lines between points at t. The number of points is defined by n (n=10 would add the lengths of lines between 0.0 and 0.1, between 0.1 and 0.2, and so on). The default n=20 is fine for most cases, usually resulting in a deviation of less than 0.01. """ length = 0 xi = x0 yi = y0 for i in range(n): t = 1.0 * (i + 1) / n pt_x, pt_y, pt_c1x, pt_c1y, pt_c2x, pt_c2y = \ curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3) c = sqrt(pow(abs(xi - pt_x), 2) + pow(abs(yi - pt_y), 2)) length += c xi = pt_x yi = pt_y return length
def point(path, t, segments=None): """Returns coordinates for point at t on the path. Gets the length of the path, based on the length of each curve and line in the path. Determines in what segment t falls. Gets the point on that segment. When you supply the list of segment lengths yourself, as returned from length(path, segmented=True), point() works about thirty times faster in a for-loop, since it doesn't need to recalculate the length during each iteration. Note that this has been deprecated: the Bezier now caches the segment lengths the moment you use them. >>> path = Bezier(None) >>> point(path, 0.0) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.moveto(0, 0) >>> point(path, 0.0) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.lineto(100, 0) >>> point(path, 0.0) Curve(LINETO, ((0.0, 0.0),)) >>> point(path, 0.1) Curve(LINETO, ((10.0, 0.0),)) """ from ..gfx.bezier import Curve if len(path) == 0: raise DeviceError, "The given path is empty" i, t, closeto = _locate(path, t, segments=segments) x0, y0 = path[i].x, path[i].y p1 = path[i+1] if p1.cmd == CLOSE: x, y = linepoint(t, x0, y0, closeto.x, closeto.y) return Curve(LINETO, ((x, y),)) elif p1.cmd == LINETO: x1, y1 = p1.x, p1.y x, y = linepoint(t, x0, y0, x1, y1) return Curve(LINETO, ((x, y),)) elif p1.cmd == CURVETO: x3, y3, x1, y1, x2, y2 = p1.x, p1.y, p1.ctrl1.x, p1.ctrl1.y, p1.ctrl2.x, p1.ctrl2.y x, y, c1x, c1y, c2x, c2y = curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3) return Curve(CURVETO, ((c1x, c1y), (c2x, c2y), (x, y))) else: raise DeviceError, "Unknown cmd for p1 %s" % p1
def point(path, t, segments=None): """Returns coordinates for point at t on the path. Gets the length of the path, based on the length of each curve and line in the path. Determines in what segment t falls. Gets the point on that segment. When you supply the list of segment lengths yourself, as returned from length(path, segmented=True), point() works about thirty times faster in a for-loop, since it doesn't need to recalculate the length during each iteration. Note that this has been deprecated: the Bezier now caches the segment lengths the moment you use them. >>> path = Bezier(None) >>> point(path, 0.0) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.moveto(0, 0) >>> point(path, 0.0) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.lineto(100, 0) >>> point(path, 0.0) Curve(LINETO, ((0.0, 0.0),)) >>> point(path, 0.1) Curve(LINETO, ((10.0, 0.0),)) """ from ..gfx.bezier import Curve if len(path) == 0: raise DeviceError, "The given path is empty" i, t, closeto = _locate(path, t, segments=segments) x0, y0 = path[i].x, path[i].y p1 = path[i + 1] if p1.cmd == CLOSE: x, y = linepoint(t, x0, y0, closeto.x, closeto.y) return Curve(LINETO, ((x, y), )) elif p1.cmd == LINETO: x1, y1 = p1.x, p1.y x, y = linepoint(t, x0, y0, x1, y1) return Curve(LINETO, ((x, y), )) elif p1.cmd == CURVETO: x3, y3, x1, y1, x2, y2 = p1.x, p1.y, p1.ctrl1.x, p1.ctrl1.y, p1.ctrl2.x, p1.ctrl2.y x, y, c1x, c1y, c2x, c2y = curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3) return Curve(CURVETO, ((c1x, c1y), (c2x, c2y), (x, y))) else: raise DeviceError, "Unknown cmd for p1 %s" % p1
def insert_point(path, t): """Returns a path copy with an extra point at t. >>> path = Bezier(None) >>> path.moveto(0, 0) >>> insert_point(path, 0.1) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.moveto(0, 0) >>> insert_point(path, 0.2) Traceback (most recent call last): ... DeviceError: The given path is empty >>> path.lineto(100, 50) >>> len(path) 2 >>> path = insert_point(path, 0.5) >>> len(path) 3 >>> path[1] Curve(LINETO, ((50.0, 25.0),)) >>> path = Bezier(None) >>> path.moveto(0, 100) >>> path.curveto(0, 50, 100, 50, 100, 100) >>> path = insert_point(path, 0.5) >>> path[1] Curve(LINETO, ((25.0, 62.5), (0.0, 75.0), (50.0, 62.5)) """ i, t, closeto = _locate(path, t) x0 = path[i].x y0 = path[i].y p1 = path[i+1] p1cmd, x3, y3, x1, y1, x2, y2 = p1.cmd, p1.x, p1.y, p1.ctrl1.x, p1.ctrl1.y, p1.ctrl2.x, p1.ctrl2.y if p1cmd == CLOSE: pt_cmd = LINETO pt_x, pt_y = linepoint(t, x0, y0, closeto.x, closeto.y) elif p1cmd == LINETO: pt_cmd = LINETO pt_x, pt_y = linepoint(t, x0, y0, x3, y3) elif p1cmd == CURVETO: pt_cmd = CURVETO pt_x, pt_y, pt_c1x, pt_c1y, pt_c2x, pt_c2y, pt_h1x, pt_h1y, pt_h2x, pt_h2y = \ curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3, True) else: raise DeviceError, "Locate should not return a MOVETO" new_path = Bezier(None) new_path.moveto(path[0].x, path[0].y) for j in range(1, len(path)): if j == i+1: if pt_cmd == CURVETO: new_path.curveto(pt_h1x, pt_h1y, pt_c1x, pt_c1y, pt_x, pt_y) new_path.curveto(pt_c2x, pt_c2y, pt_h2x, pt_h2y, path[j].x, path[j].y) elif pt_cmd == LINETO: new_path.lineto(pt_x, pt_y) if path[j].cmd != CLOSE: new_path.lineto(path[j].x, path[j].y) else: new_path.closepath() else: raise DeviceError, "Didn't expect pt_cmd %s here" % pt_cmd else: if path[j].cmd == MOVETO: new_path.moveto(path[j].x, path[j].y) if path[j].cmd == LINETO: new_path.lineto(path[j].x, path[j].y) if path[j].cmd == CURVETO: new_path.curveto(path[j].ctrl1.x, path[j].ctrl1.y, path[j].ctrl2.x, path[j].ctrl2.y, path[j].x, path[j].y) if path[j].cmd == CLOSE: new_path.closepath() return new_path
def insert_point(path, t): """Returns a path copy with an extra point at t. >>> path = BezierPath(None) >>> path.moveto(0, 0) >>> insert_point(path, 0.1) Traceback (most recent call last): ... NodeBoxError: The given path is empty >>> path.moveto(0, 0) >>> insert_point(path, 0.2) Traceback (most recent call last): ... NodeBoxError: The given path is empty >>> path.lineto(100, 50) >>> len(path) 2 >>> path = insert_point(path, 0.5) >>> len(path) 3 >>> path[1] PathElement(LINETO, ((50.0, 25.0),)) >>> path = BezierPath(None) >>> path.moveto(0, 100) >>> path.curveto(0, 50, 100, 50, 100, 100) >>> path = insert_point(path, 0.5) >>> path[1] PathElement(LINETO, ((25.0, 62.5), (0.0, 75.0), (50.0, 62.5)) """ i, t, closeto = _locate(path, t) x0 = path[i].x y0 = path[i].y p1 = path[i + 1] p1cmd, x3, y3, x1, y1, x2, y2 = p1.cmd, p1.x, p1.y, p1.ctrl1.x, p1.ctrl1.y, p1.ctrl2.x, p1.ctrl2.y if p1cmd == CLOSE: pt_cmd = LINETO pt_x, pt_y = linepoint(t, x0, y0, closeto.x, closeto.y) elif p1cmd == LINETO: pt_cmd = LINETO pt_x, pt_y = linepoint(t, x0, y0, x3, y3) elif p1cmd == CURVETO: pt_cmd = CURVETO pt_x, pt_y, pt_c1x, pt_c1y, pt_c2x, pt_c2y, pt_h1x, pt_h1y, pt_h2x, pt_h2y = \ curvepoint(t, x0, y0, x1, y1, x2, y2, x3, y3, True) else: raise NodeBoxError, "Locate should not return a MOVETO" new_path = BezierPath(None) new_path.moveto(path[0].x, path[0].y) for j in range(1, len(path)): if j == i + 1: if pt_cmd == CURVETO: new_path.curveto(pt_h1x, pt_h1y, pt_c1x, pt_c1y, pt_x, pt_y) new_path.curveto(pt_c2x, pt_c2y, pt_h2x, pt_h2y, path[j].x, path[j].y) elif pt_cmd == LINETO: new_path.lineto(pt_x, pt_y) if path[j].cmd != CLOSE: new_path.lineto(path[j].x, path[j].y) else: new_path.closepath() else: raise NodeBoxError, "Didn't expect pt_cmd %s here" % pt_cmd else: if path[j].cmd == MOVETO: new_path.moveto(path[j].x, path[j].y) if path[j].cmd == LINETO: new_path.lineto(path[j].x, path[j].y) if path[j].cmd == CURVETO: new_path.curveto(path[j].ctrl1.x, path[j].ctrl1.y, path[j].ctrl2.x, path[j].ctrl2.y, path[j].x, path[j].y) if path[j].cmd == CLOSE: new_path.closepath() return new_path