Zernike polynomials in NumPy: a telescope's wavefront and Strehl ratio

Afterwards you can fit a measured wavefront with Zernike polynomials in NumPy, name its aberrations, and get the Strehl ratio from the point spread function.