code: def _interpolate(self, kind='linear'):
          
          # drop any existing non-dimension coordinates along the resampling
          # dimension
          dummy = self._obj.copy()
          for k, v in self._obj.coords.items():
              if k != self._dim and self._dim in v.dims:
                  dummy = dummy.drop(k)
          return dummy.interp(assume_sorted=True, method=kind,
                              kwargs={'bounds_error': False},
                              **{self._dim: self._full_index})
docs: Apply scipy.interpolate.interp1d along resampling dimension.