code: def estimate_flux(self):
          
  
          from .epsf import _interpolate_missing_data
  
          if np.any(self.mask):
              data_interp = _interpolate_missing_data(self.data, method='cubic',
                                                      mask=self.mask)
              data_interp = _interpolate_missing_data(data_interp,
                                                      method='nearest',
                                                      mask=self.mask)
              flux = np.sum(data_interp, dtype=np.float64)
  
          else:
              flux = np.sum(self.data, dtype=np.float64)
  
          return flux
docs: Estimate the star's flux by summing values in the input cutout
          array.
  
          Missing data is filled in by interpolation to better estimate
          the total flux.