def compute_join_intensity(im, px, py):
    tt = linspace(0, 1, 50)
    x = array([round(polyval(px, t)) for t in tt])
    y = array([round(polyval(px, t)) for t in tt])
    return _compute_intensity(im, x, y)
def compute_join_intensity(im, px, py):
    tt = linspace(0, 1, 50)
    x = array([round(polyval(px, t)) for t in tt])
    y = array([round(polyval(px, t)) for t in tt])
    return _compute_intensity(im, x, y)