EigenCWD: a spatially varying deconvolution algorithm for single metalens imaging

The miniaturization of optics through the use of two-dimensional metalenses has enabled novel applications in imaging. To date, single-lens imaging remains the most common configuration, partly due to the limited focusing efficiency of metalenses. This results in limitations when it comes to wavefront manipulation and, thus, unavoidable aberrations in the formed image that require computational deconvolution to deblur the image. For certain lens profiles, such as the most common hyperbolic one that results in the highest efficiencies, at large fields of view, spatially-varying aberrations such as coma or astigmatism are prominent. These aberrations cannot be corrected by traditional deconvolution methods, such as Wiener filtering. Here, we develop a spatially varying deconvolution algorithm based on eigenvalue column-wise decomposition (eigenCWD). EigenCWD solves a minimization problem of the error between the measured image and the estimated image of the object to be reconstructed through an approximate forward blurring model. This approximate forward model uses an eigendecomposition of the spatially varying point spread functions for fast computation, allowing for efficient scaling to larger image sizes and blurring kernels common in metalens imaging. We demonstrate eigenCWD's ability to correct spatially varying blur and distortions for various lens profiles, surpassing that of the Wiener filter.