Gaussianization

import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns

sns.set_style('white')

R = np.linspace(0,1,1.0)

X = np.asarray([],dtype=np.float)
Y = np.asarray([],dtype=np.float)
for idx,r in enumerate(R):
    x_coords = np.zeros((100,))
    y_coords = np.zeros((100,))
    for i,th in enumerate(np.linspace(0,2*np.pi,100)):
        x_coords[i] = (1+0.05*np.random.randn())*r*np.cos(th)
        y_coords[i] = (1+0.05*np.random.randn())*r*np.sin(th)
    X=np.append(X,x_coords)
    X=np.append(X,x_coords)
    Y=np.append(Y,y_coords)
    Y=np.append(Y,-y_coords)

plt.figure()
plt.plot(X,Y,'.k')