Abstract:Large-aperture, lightweight, and high-resolution imaging are hallmarks of major optical systems. To eliminate aberrations, traditional systems are often bulky and complex, whereas the small volume and light weight of diffractive lenses position them as potential substitutes. However, their inherent diffraction mechanism leads to severe dispersion, which limits their application in wide spectral bands. Addressing the dispersion issue in diffractive lenses, we propose a chromatic aberration correction algorithm based on compressed sensing. Utilizing the diffractive lens's focusing ability at the reference wavelength and its degradation performance at other wavelengths, we employ compressed sensing to reconstruct images from incomplete image information. In this work, we design a harmonic diffractive lens with a diffractive order of M=150, an aperture of 40 mm, a focal length f0=320 mm, a reference wavelength λ0=550 nm, a wavelength range of 500–800 nm, and 7 annular zones. Through algorithmic recovery, we achieve clear imaging in the visible spectrum, with a peak signal-to-noise ratio (PSNR) of 22.85 dB, a correlation coefficient of 0.9596, and a root mean square error (RMSE) of 0.02, verifying the algorithm's effectiveness.
engineering, electrical & electronic,chemistry, analytical,instruments & instrumentation
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the application of **chromatic aberration correction** in diffractive lenses. Specifically, due to their inherent diffraction mechanism, diffractive lenses have a serious dispersion phenomenon in a wide spectral range, which limits their application in high - precision imaging. To overcome this problem, the author proposes a chromatic aberration correction algorithm based on Compressed Sensing (CS).
### Background and Problem of the Paper
Traditional optical systems use multiple lenses to eliminate aberrations, but these systems are usually bulky, complex in structure and costly. In contrast, diffractive lenses have the advantages of being ultra - thin and lightweight, but their inherent diffraction mechanism leads to significant dispersion, making it difficult to achieve high - precision imaging in a wide spectral range. Traditional solutions include adding reverse - power diffractive lenses or designing multi - layer diffractive lens structures, but these methods still have limitations.
### Proposed Method
The author proposes a chromatic aberration correction algorithm based on Compressed Sensing. This algorithm utilizes the focusing ability of the diffractive lens at the reference wavelength and the degraded performance at other wavelengths to reconstruct the image from incomplete image information through Compressed Sensing technology. The specific steps are as follows:
1. **Design of Harmonic Diffractive Lens**: A 150 - order harmonic diffractive lens is designed, with an aperture of 40 mm, a focal length of 320 mm, a reference wavelength of 550 nm, a wavelength range of 500 - 800 nm, and contains 7 annular regions.
2. **Image Capture and Reconstruction**: Images are captured through the prototype system, and the Compressed Sensing algorithm is used to perform image reconstruction in other incomplete channels.
3. **Evaluation of Results**: The effectiveness of the algorithm is verified by indicators such as Peak Signal - to - Noise Ratio (PSNR), correlation coefficient and Root Mean Square Error (RMSE).
### Key Formulas
1. **Point Spread Function (PSF)**:
\[
p(x,y;\lambda) = A \int \int P(u,v;\lambda) e^{-j \frac{2\pi}{\lambda z_i} (ux + vy)} du dv
\]
where \( A \) is the amplitude constant, \( z_i \) is the distance from the lens to the imaging plane, \((u,v)\) are the coordinates on the lens plane, and \( P(u,v;\lambda) \) is the aperture function.
2. **Aperture Function**:
\[
P(u,v;\lambda) = \text{Circ}(u,v) e^{j \Phi(u,v)}
\]
where \(\text{Circ}(\cdot)\) is the circular aperture function and \(\Phi(u,v)\) is the phase delay introduced on the lens surface.
3. **Incoherent Imaging Model**:
\[
A(i(x,y;\lambda)) = i(x,y;\lambda) \otimes |p(x,y;\lambda)|^2
\]
where \(\otimes\) represents the two - dimensional convolution operation and \( A(i(x,y;\lambda)) \) represents the image affected by chromatic aberration.
4. **Compressed Sensing Optimization Problem**:
\[
f(x) = \arg \min \left( \frac{1}{2} \| y - Kx \|^2 + \lambda \Phi_{\text{niTV}}(x) \right)
\]
where \( K \) is the linear operation matrix, \(\| \cdot \| \) represents the norm, \(\lambda \in [0, +\infty)\) is the regularization parameter, and \(\Phi_{\text{niTV}}(x)\) is the non - isotropic total variation regularization method.
5. **Non - Isotropic Total Variation Regularization**:
\[
\Phi_{\text{niTV}}(x) = \sum_i \left| \Delta_h x_i \right| + \left| \Delta_v x_i \right|
\]
where \(\Delta\)