Healpix Bug

[AaronParsons]

I ran across a bug awhile ago when running the following code:

        img2d = aipy.img.Img(size=300,res=.5)   #creates a 2D image
        crd = img2d.get_eq(ra=numpy.pi,dec = [some negative value])     #gets equatorial coordinates of the image centered on a specific RA/dec
        x,y,z = crd
        x,y,z = x.flatten(), y.flatten(), z.flatten()

        img3d = aipy.map.Map(fromfits = ' ')      #reads in Haslam map path
        fluxes = img3d[x,y,z]                   #gets values of Haslam map at specified coordinates

The last line generates an error, only when I used RA = exactly pi and DEC = some negative value in the 'get_eq' function. (Dec = 0 works fine).

The error message is always something like this: "index 3312128 is out of bounds for size 3145728."

The Healpix map has a size of 3145728 pixels.

[AaronParsons]

I can confirm that this issues is resolved in the migration to healpy.