estimate_flux

xrayvision.mem.estimate_flux(vis, shape, pixel, maxiter=1000, tol=0.001)

Estimate the total flux in the image by solving an optimisation problem.

This function estimates the total flux of an event by solving the problem

\[\underset{\chi}{\mathrm{argmin}}(f) = \sum \left (\Re V(f) - \Re v))^2 + ( \Im V(f) - \Im v)^2 ) / \sigma^2 \right )\]

subject to \(f >= 0\),

the algorithm finds a positive image \(f\) that minimizes the \(\chi\) square function. The estimation of the total flux is then obtained by computing the total fux of \(f\). The method implemented \(f\) the minimization is projected Landweber.

Parameters:

• vis – Input visibilities

• shape – Shape of image

• pixel – Size of pixels

• maxiter (int) – Maximum number of iterations

• tol – Tolerance at which to stop

Returns:

Estimated total flux