Fisheye warping

[Stéfan van der Walt]

Hi all,

After yesterday's post on warping, another teaser:

http://mentat.za.net/refer/lena_warp.png

From:

def fisheye(xy):
    center = np.mean(xy, axis=0)
    c = xy - center
    r = np.sqrt(np.sum(c**2, axis=1))
    theta = np.arctan2(c[:, 1], c[:, 0])

    r = 0.87 * np.exp(r**(1/2.5) / 1.75)

    return np.column_stack((
        r * np.cos(theta), r * np.sin(theta)
        )) + center


Stéfan