code: def get_pixel_indices(self, lats, lons):
          
          theta = np.radians(90. - lats)
          phi = np.radians(lons)
          return hp.ang2pix(self.nside, theta, phi, self.nest)
docs: "Return the indices in the flat array corresponding to a set of coordinates