2023-02-04 00:46:23 +00:00
|
|
|
from svgpathtools import disvg, Line, CubicBezier
|
|
|
|
|
from scipy.optimize import fminbound
|
2017-03-30 09:22:23 +00:00
|
|
|
|
|
|
|
|
# create some example paths
|
|
|
|
|
path1 = CubicBezier(1,2+3j,3-5j,4+1j)
|
|
|
|
|
path2 = path1.rotated(60).translated(3)
|
|
|
|
|
|
2023-02-04 00:46:23 +00:00
|
|
|
|
2017-03-30 09:22:23 +00:00
|
|
|
def dist(t):
|
|
|
|
|
return path1.radialrange(path2.point(t))[0][0]
|
2023-02-04 00:46:23 +00:00
|
|
|
|
|
|
|
|
|
|
|
|
|
# find minimizer
|
2017-03-30 09:22:23 +00:00
|
|
|
T2 = fminbound(dist, 0, 1)
|
|
|
|
|
|
|
|
|
|
# Let's do a visual check
|
|
|
|
|
pt2 = path2.point(T2)
|
|
|
|
|
T1 = path1.radialrange(pt2)[0][1]
|
|
|
|
|
pt1 = path1.point(T1)
|
|
|
|
|
disvg([path1, path2, Line(pt1, pt2)], 'grb', nodes=[pt1, pt2])
|
