add nest optional arg in get_healpix_centers

src/gaiaunlimited/utils.py

@@ -6,19 +6,23 @@
 __all__ = ["coord2healpix", "get_healpix_centers"]
 
 
-def get_healpix_centers(order):
+def get_healpix_centers(order,nest=False):
     """Get centers of HEALPix as astropy coordinates.
 
     Args:
         order (int): The order of the pixelisation
 
+    Optional args:
+        nest (boolean): ``False`` for RING healpix scheme, ``True`` for nested.
+            default is ``False``.
+
     Returns:
         astropy.coordinates.SkyCoord: coordinates of the centers.
     """
     nside = hp.order2nside(order)
     npix = hp.order2npix(order)
     ipix = np.arange(npix)
-    ra, dec = hp.pix2ang(nside, ipix, lonlat=True)
+    ra, dec = hp.pix2ang(nside, ipix, lonlat=True, nest=nest)
     return coord.SkyCoord(ra * u.deg, dec * u.deg)