Summarystats: Wavefront summaries
In matwey/Rfringe: Rfringe

Description Usage Arguments Value Note Author(s) See Also Examples

Estimate the RMS or P-V wavefront error over the pupil given in the argument.

1
2
3

pupilrms(pupil)
pupilpv(pupil)
strehlratio(rms)

`pupil`	`pupil` is the matrix created by the call to `pupil`
`rms`	The rms wavefront error

An estimate of the RMS or P-V error of the wavefront, or Strehl ratio.

The function pupilrms simply returns the standard deviation of the defined values in pupil, which is a crude but usually good enough approximation to the properly defined integral over the aperture.

For a wavefront defined entirely in terms of a vector zcoef of Zernike coefficients rms <- sqrt(crossprod(zcoef)) is faster and more accurate.

pupilpv does the obvious. There is no analytical solution in general for P-V.

strehlratio computes Mahajan's approximation to the Strehl ratio.

M.L. Peck mpeck1@ix.netcom.com

pupil.

# A random vector of Zernike coefficients

zcoef <- rnorm(length(zlist.qf$n), mean=0, sd=0.01)
wf <- pupil(zcoef=zcoef)
image(wf, col=topo.colors(256), asp=1)
contour(wf, add=TRUE)
pupilrms(wf)
sqrt(crossprod(zcoef)) #should be the same to about 4 digits
pupilpv(wf)
strehlratio(sqrt(crossprod(zcoef))) #probably around 0.8