[SciPy-User] Re[SciPy-user] binning to polar coordinates

[arsbbr]

Hi,

I'm trying to rebin some uniform gridded data to polar coordinates. The
original data (Z) is noisy but on a fine grid. 

### snip

# define grid in cartesian coordinates
x = arange(-500, 500)
y = arange(-500, 500)
X, Y = meshgrid(x, y)
Z = X**2 + Y**2
Z += uniform(-5000, 5000, size=Z.shape))

# transform to polar coordinates

def cart2pol(x, y):
    theta = arctan2(y, x)
    rho = sqrt(x**2 + y**2)
    return (theta, rho)  

THETA_xy, RHO_xy = cart2pol(X, Y)

# define new polar grid
theta = arange(-pi, pi, 0.1)
rho = arange(0.1 ,500, 2.0)    

THETA, RHO = meshgrid(theta, rho)

# Tried to use griddata, but it gives of course very jumpy results, because
it 
# does not take the average over the pixels in a patch (dTHETA, dRHO) but an
interpolation at the
# exact new point.

interpolate.griddata((THETA_xy.ravel(), RHO_xy.ravel()), 
   Z.ravel(), (THETA, RHO), method='linear')
        
### snip

Is there a method that rebins the data to a new grid, possibly taking the
average 
of the Z-values inside a patch (dTHETA, dRHO)?

Thanks for any hints!



-- 
View this message in context: http://old.nabble.com/Rebinning-to-polar-coordinates-tp34348219p34348219.html
Sent from the Scipy-User mailing list archive at Nabble.com.