code: def get_map_values(self, lons, lats, ibin=None):
          
          pix_idxs = self.get_pixel_indices(lons, lats, ibin)
          idxs = copy.copy(pix_idxs)
  
          m = np.empty_like(idxs[0], dtype=bool)
          m.fill(True)
          for i, p in enumerate(pix_idxs):
              m &= (pix_idxs[i] >= 0) & (pix_idxs[i] < self._npix[i])
              idxs[i][~m] = 0
  
          vals = self.counts.T[idxs]
          vals[~m] = np.nan
          return vals
docs: Return the map values corresponding to a set of coordinates.
  
          Parameters
          ----------
          lons  : array-like
             'Longitudes' (RA or GLON)
  
          lats  : array-like
             'Latitidues' (DEC or GLAT)
  
          ibin : int or array-like
             Extract data only for a given energy bin.  None -> extract data for all bins
  
          Returns
          ----------
          vals : numpy.ndarray((n))
             Values of pixels in the flattened map, np.nan used to flag
             coords outside of map