segmentize.py

def segmentize(points, km):
    npoints = len(points)
    if npoints == 0:
        return []
    elif npoints == 1:
        return [points[0]]

    a, b = points[:2]
    rest = segmentize(points[2:], km)
    if haversine(a, b) > km:
        ax, ay = a
        bx, by = b
        c = ((bx + ax) / 2., (by + ay) / 2.)
        left = segmentize([a, c], km)
        right = segmentize([c, b], km)
        return left + right + rest
    else:
        return [a, b] + rest

👀 ☝️

⚫ ⚫ ⚪ ⚫