docs: def _infer_interval_breaks(coord, kind=None):
          """
          Interpolate the bounds from the data in coord
  
          Parameters
          ----------
          %(CFDecoder.get_plotbounds.parameters.no_ignore_shape)s
  
          Returns
          -------
          %(CFDecoder.get_plotbounds.returns)s
  
          Notes
          -----
          this currently only works for rectilinear grids"""
          if coord.ndim == 1:
              return _infer_interval_breaks(coord)
          elif coord.ndim == 2:
              from scipy.interpolate import interp2d
              kind = kind or rcParams['decoder.interp_kind']
              y, x = map(np.arange, coord.shape)
              new_x, new_y = map(_infer_interval_breaks, [x, y])
              coord = np.asarray(coord)
              return interp2d(x, y, coord, kind=kind, copy=False)(new_x, new_y)