code: def get_sky_coords(self):
          
          if self._ipix is None:
              theta, phi = hp.pix2ang(
                  self._nside, list(range(self._npix)), self._nest)
          else:
              theta, phi = hp.pix2ang(self._nside, self._ipix, self._nest)
  
          lat = np.degrees((np.pi / 2) - theta)
          lon = np.degrees(phi)
          return np.vstack([lon, lat]).T
docs: Get the sky coordinates of all the pixels in this pixelization