code: def _xy_idx(self):
          
  
          yidx, xidx = np.indices(self._data.shape)
          return xidx[~self.mask].ravel(), yidx[~self.mask].ravel()
docs: 1D arrays of x and y indices of unmasked pixels in the cutout
          reference frame.