code: def interpolate(self, lon, lat, egy=None, interp_log=True):
          
  
          if self.data.ndim == 1:
              theta = np.pi / 2. - np.radians(lat)
              phi = np.radians(lon)
              return hp.pixelfunc.get_interp_val(self.counts, theta,
                                                 phi, nest=self.hpx.nest)
          else:
              return self._interpolate_cube(lon, lat, egy, interp_log)
docs: Interpolate map values.
  
          Parameters
          ----------
          interp_log : bool
              Interpolate the z-coordinate in logspace.