Zernike Polynomials vs Fourier Wavefront Reconstruction: High-Spatial-Frequency Aberrations & Irregular Corneas

Zernike polynomials $Z_n^m(\rho, \theta)$ are mathematically orthogonal over the unit disk:

$$\int_0^1 \int_0^{2\pi} Z_n^m(\rho, \theta) \cdot Z_{n'}^{m'}(\rho, \theta) \, \rho \, d\rho \, d\theta = \frac{\pi}{n + 1} \delta_{nn'} \delta_{mm'}$$

This orthogonality allows each optical aberration (defocus, astigmatism, coma, spherical) to be evaluated and corrected independently. However, truncating the expansion at the 4th, 6th, or 8th order acts as an aggressive low-pass spatial filter, smoothing out sharp localized scars and keratoconic cones.

Fourier Reconstruction treats the measured wavefront slopes as a continuous 2D spatial vector field $\vec{S}(x, y)$, computing the phase via inverse Fast Fourier Transform (FFT):

$$W(x, y) = \mathcal{F}^{-1} \left\{ \frac{i \cdot 2\pi (u \cdot S_x(u, v) + v \cdot S_y(u, v))}{4\pi^2 (u^2 + v^2)} \right\}$$

Fourier transforms do not assume a circular boundary and do not average away local aberrations. They resolve localized corneal indentations, trauma scars, and abrupt surgical transition zones with exceptional fidelity.