Switch to coord_to_idx.
@@ -484,7 +484,7 @@ class Map(object):
msk = self.geom.contains(coords)
vals = np.empty(coords.shape, dtype=self.data.dtype)
coords = coords.mask(msk)
- idx = self.geom.coord_to_pix(coords)
+ idx = self.geom.coord_to_idx(coords)
vals[msk] = self.get_by_idx(idx)
vals[~msk] = np.nan
return vals
