docs: def cutout_centroid(self):
          """
          The ``(y, x)`` coordinate, relative to the `data_cutout`, of
          the centroid within the source segment.
          """
  
          m = self.moments
          if m[0, 0] != 0:
              ycentroid = m[1, 0] / m[0, 0]
              xcentroid = m[0, 1] / m[0, 0]
              return (ycentroid, xcentroid) * u.pix
          else:
              return (np.nan, np.nan) * u.pix