Opto-Electronic Science, 2022, 1 (3): 210006, Published Online: May. 5, 2022  

Single-shot mid-infrared incoherent holography using Lucy-Richardson-Rosen algorithm

Author Affiliations
1 Optical Sciences Center and ARC Training Centre in Surface Engineering for Advanced Materials (SEAM), School of Science, Computing and Engineering Technologies, Optical Sciences Center, Swinburne University of Technology, Hawthorn, Melbourne, Victoria 3122, Australia
2 Institute of Physics, University of Tartu, 50411 Tartu, Estonia
3 Infrared Microspectroscopy (IRM) Beamline, ANSTO – Australian Synchrotron, Clayton, Victoria 3168, Australia
4 Tokyo Tech World Research Hub Initiative (WRHI), School of Materials and Chemical Technology, Tokyo Institute of Technology, 2-12-1, Ookayama, Meguro-ku, Tokyo 152-8550, Japan
In recent years, there has been a significant transformation in the field of incoherent imaging with new possibilities of compressing three-dimensional (3D) information into a two-dimensional intensity distribution without two-beam interference (TBI). Most of the incoherent 3D imagers without TBI are based on scattering by a random phase mask exhibiting sharp autocorrelation and low cross-correlation along the depth. Consequently, during reconstruction, high lateral and axial resolutions are obtained. Imaging based on scattering requires an astronomical photon budget and is therefore precluded in many power-sensitive applications. In this study, a proof-of-concept 3D imaging method without TBI using deterministic fields has been demonstrated. A new reconstruction method called the Lucy-Richardson-Rosen algorithm has been developed for this imaging concept. We believe that the proposed approach will cause a paradigm-shift in the current state-of-the-art incoherent imaging, fluorescence microscopy, mid-infrared fingerprinting, astronomical imaging, and fast object recognition applications.

1 Introduction

In the past, holography, and three-dimensional imaging (3D) techniques were mostly powered by coherent light sources since spatially incoherent illumination demanded cumbersome optical architectures such as rotational shearing interferometers, triangle interferometers and conoscopic holography1. The need for such optical architectures is derived from the requirement that two-beam interference (TBI) was needed to compress 3D information in a two-dimensional (2D) hologram. The birth of active optical devices such as spatial light modulators and computer processing methods, eased the constraints of optical architectures of incoherent holography with TBI. One milestone in this research direction was the development of Fresnel incoherent correlation holography (FINCH), which is one of the widely used incoherent holography architectures today2-4. In 1968, deconvolution based 2D coded aperture imaging (CAI) using randomly arranged pinholes was demonstrated by Dicke and Ables to overcome the low-intensity problem when using a single pinhole for X-ray and gamma-ray imaging5, 6. However, the study was not explored beyond 2D imaging. In 2008, a 3D imaging technique based on CAI was demonstrated with wavelength as the third dimension7. Later, the CAI method was extended to 3D imaging and four-dimensional (4D) imaging with depth as the third dimension8-10, and depth and wavelength as the third and fourth dimensions, respectively11, 12. Therefore, CAI offers multidimensional imaging capabilities without TBI and with a minimum number of optical components, which makes them more attractive than FINCH-based techniques.

This research project began when we attempted the CAI method employing the benchtop Fourier transform infrared microspectroscopy (FTIRm) system at the Australian Synchrotron operating with a conventional GlobarTM IR source. This instrument, which comprises a Bruker Hyperion 3000 IR microscope and V70 FTIR spectrometer, is used in support of many applications at the Infrared Microspectroscopy (IRM) beamline, whereby the identification and analysis of functional groups within a sample enable a better understanding of diverse materials from diseased tissue to composite materials and foods to paintings13-15. We applied the CAI method to the benchtop FTIR instrument at the IRM beamline by manufacturing scattering lenses on calcium fluoride (CaF2) substrates. A matching 15× (NA = 0.4) Cassegrain objective lens (COL) and a 64×64 focal plane array (FPA) imaging detector (pixel size = 40 μm) were used. During this study, we noticed that COL exhibited interesting focal characteristics with sharp autocorrelation functions and low cross-correlation with respect to other planes (see Supplementary information Section 1, Section 2).

Let us ask this fundamental question – why is scattering based coded apertures preferred for 3D imaging? The reason is that the average speckle size formed by scattering is approximately the diffraction limited spot size16. Consequently, the autocorrelation function is sharp and about twice that of the size of the focal spot that can be obtained with a lens of the same NA1. Secondly, the speckle patterns change with depth resulting in the capability to discriminate different planes along depth by cross-correlation. Therefore, for 3D imaging, a sharp autocorrelation and low cross-correlation along depth (SALCAD) is needed. The fields generated using scatterers can be classified as random SALCAD fields and the one from COL as a deterministic SALCAD field. Deterministic SALCAD fields require a lower photon budget than the random counterparts and the point spread functions (PSFs) can be calculated, unlike random SALCAD fields where PSFs are recorded9. So, this imaging concept can be considered as correlation holography when the SALCAD condition is satisfied and the object is illuminated by a spatially incoherent source. Recently, a modified approach was employed where sparse randomly arranged focal spots were generated instead of scattering to improve the signal to noise ratio17. However, this method is suitable only for imaging a single plane, and to image multiple planes, many such phase masks are needed to be spatially multiplexed. Fresnel zone aperture based 2D imaging systems have been reported in the past18-20.

2 Results and discussion

The FTIRm consists of collinear dual beams generated from two sources: (i) a standard white light source for bright-field visible light observation of the sample plane and (ii) a broadband GlobarTM IR source that allows the acquisition of mid-IR (MIR) spectral images in the range ~ 899 cm–1–3845 cm–1 as shown in Fig. 1(a). The FTIR single-beam spectrum obtained from the GlobarTM source is shown in Fig. 1(b) in which a total of 765 spectral channels equally separated by an average value of 3.85 cm–1 were used. The MIR beam from the FTIR spectrometer is focused on the sample plane and the scattered signal is recorded by the FPA detector. With 15 × COL, the pixel pitch is ~2.7 μm in the sample plane and the field of view is around ~172 × 172 μm2.

Fig. 1. Simplified schematic of the FTIR microspectroscopic system and the measured spatial and spectral characteristics. (a) Abbreviations: COL – Cassegrain objective lens, IR – Infrared, BS – beam splitter, M – mirror, L – lens, FTIR – Fourier transform infrared spectrometer, FPA – focal plane array, MSP – Motorized sliding plate, MCT - mercury-cadmium-telluride, A – aperture, MIR – mid infrared, PSF – point spread function. The MIR from the FTIR spectrometer is sent into the IR/VISIBLE microscope. The microscope is aligned collinearly for both MIR and reference visible light. There are two modes of operation: (reflection and transmission) and two modes of data recording (mapping using single element MCT detector and imaging using multi-element FPA detector). (b) Spectral profile of the GlobarTM source with 14 spectral channels marked in blue are selected for the study. (ce) 4D plot of the PSFs recorded for an axial aberration Δz = 0, 125 μm and 250 μm, respectively. Spectral cube data PSF(x, y, ν) was formed from the 14 spectral channels each averaged over 100 spectral channels –50 to +50 channels. The non-changing pattern along ν indicates spectral aberrations suppressed by the COL. (fh). Autocorrelation cube data of (c–e), respectively. The thin cylindrical structures obtained for Δz = 125 μm and 250 μm indicates the possibility of imaging in deconvolution mode with a similar resolving power as Δz = 0 μm. (i) Depth cube data PSF(x, y, Δz) from 0 to 250 μm. Unlike a regular lens where the intensity distribution is filled, the COL generates hollow intensity distribution with four distinct lobes causing a sharp autocorrelation function. (j) The spectral and spatial cube data were combined to form 8415 spatio-spectral images corresponding to the 14 spectral channels and aberrations 0 to 250 μm in steps of 25 μm and was cross-correlated with the image corresponding to the central spatio-spectral point (2307 cm–1, 125 μm) and cross-correlation value at the origin (x=0, y=0) is plotted. The FWHM along z was ~50 μm indicating that the intensity of a plane 50 μm apart has a loss of about 50% of intensity which is ideal for 3D MIR imaging applications.

下载图片 查看所有图片

The IR light produced by the GlobarTM source is spatially incoherent and so the light emitted from one object point does not interfere with the light emitted from another point, instead their intensities add up. As a result, this FTIRm system is a shift-invariant imager, linear in intensity, i.e., PSF and the intensity distribution obtained for an object O can be expressed as OH=OPSF+σ, where σ is the noise and ‘ ’ is a 2D convolutional operator (see Supplementary information Section 1, Section 2). Therefore, if the PSF is known for all axial shifts Δz, then the 3D image of O can be reconstructed by solving the above equation. For instance, object O=O(z1,λa)+O(z2,λa)+OT, which absorbs at λa at z1 and z2 while OT refers to non-absorbing layers. The reconstruction of the chemical image of O at plane z1 is given as OR(λa)=OH*PSF(z1)=O(z1,λa)PSF(z1)*PSF(z1)+O(z2,λa)PSF(z2)*PSF(z1)+N, where ‘ *’ is a 2D correlational operator and N is the reconstruction noise. The reconstructed image is O(z1,λa) sampled by the autocorrelation function PSF(z1)*PSF(z1), which is a Delta-like function with a minimum width ~1.22λ/NA added to O(z2,λa) sampled by PSF(z1)*PSF(z2), which is a blurred version of O(z2,λa).

The PSFs were recorded using a 15 μm pinhole for Δz = 0 to 250 μm in steps of 25 μm. The lateral resolution of the system is therefore ~ 15 μm. The spectral data cubes PSF(x,y,ν) for Δz = 0, 125 μm and 250 μm along 14 averaged (~100 channels) spectral channels are shown in Fig. 1(c–e), respectively, whose non-changing images along ν indicate the suppression of spectral aberrations. The data structure of the output generated by the OPUS 8.0 software was transformed into cube data for processing in MATLAB. The COL generates sharp annular intensity distributions with four peaks, which are responsible for the sharp autocorrelation functions as shown in Fig. 1(f–h). Only a slight increase in the width of the autocorrelation functions was noticed for larger values of Δz. The depth cube data PSF(x,y,Δz) was generated from recorded PSFs for Δz = 0 to 250 μm and averaged over all the 765 spectral channels as shown in Fig. 1(i). The spatial and spectral cubes were combined, cross-correlated with the central spatio-spectral point (2307 cm–1, 125 μm) as shown in Fig. 1(j). It is seen that the FWHM along the depth axis is sharp at ~50 μm with mild fluctuations along the spectral axis. The above spatio-spectral characteristics of the FTIRm system are highly desirable for 3D chemical imaging.

The next important aspect of imaging is the reconstruction technique. It is well-established from previous studies that the PSF is not the optimal reconstruction function by cross-correlation21, and a non-linear reconstruction (NLR) method was introduced by Rosen in which the magnitudes of the spatial frequency spectrum of the PSF and OH are tuned in powers of α and β respectively in R=|F1exp[iarg(OH˜)]}| until minimum entropy was obtained (see Supplementary information Section 3)22. Another approach is the Lucy-Richardson algorithm (LRA), which consists of a forward convolution between the initial guess solution and the PSF, which is compared with the recorded OH and the ratio is backward cross-correlated with the PSF and the result is multiplied to the initial guess. This process is iterated until the maximum likelihood object function is reconstructed23, 24. The (n+1)th reconstructed image is given as Rn+1=Rn{OHRnPSFPSF'}, where PSF' is the 180 degrees rotated version and the complex conjugate of PSF. In this study, it was noticed that neither NLR nor LRA were optimal; however, when the linearity was broken of the backward cross-correlation of LRA by the NLR approach, a rapid convergence and significantly better estimation was obtained. The schematic of the Lucy-Richardson-Rosen algorithm (LRRA) is shown in Fig. 2(a) (see Supplementary information Section 3).

Fig. 2. Reconstruction results using Lucy-Richardson algorithm, Non-linear reconstruction and Lucy-Richardson-Rosen algorithm. (a) Lucy-Richardson-Rosen algorithm, O – object, OH – object hologram, PSF – point spread function, OTF – optical transfer function, × - element by element product, - 2D convolution, - Fourier transform, * - complex conjugate operation following a Fourier transform, 1- inverse Fourier transform, ~ Fourier transform of a variable, Rn and R(n+1) are the nth and (n+1)th solutions and OH was used as the initial guess solution R1, ML(O) – maximum likelihood solution of the object. (b) Direct imaging of the cross object using a visible light channel of the microscope. Average image of all the 765 spectral channels of the cross object at (c) Δz = 0. (d) Δz = 100 μm, (e) The PSF recorded for Δz = 100 μm. Reconstruction results of the cross object using (f) LRA (150 iterations), (g) NLR (α = 0, β = 0.6) and (h) LRRA (20 iterations, α = 0.2, β = 1). (i) Direct MIR imaging result of four pinholes Δz = 0. (j) Intensity pattern recorded for four pinholes for Δz = 150 μm and (k). its reconstruction result. (l) Intensity pattern recorded for four pinholes for Δz = 200 μm and (m). its reconstruction result. (n) Bar plots of RMSE and entropy for LRA, NLR and LRRA. The performance of LRRA is found to be significantly higher than LRA and NLR with RMSE ratio of 13 and 5 times, respectively. The LRRA has a rapid convergence with 7 times faster than LRA for the cross object. It must be noted that the above values are for only the cross object. The performances vary with the complexity of the object with the minimum difference for simple objects and significant difference for complicated objects. The scale bar is 50 μm.

下载图片 查看所有图片

The spatio-spectral aberrations of the system have been studied25. A cross shaped object (150 μm × 150 μm) and four random pinholes (50 μm in diameter) ablated on chromium-gold layer coated calcium fluoride substrate were used as test objects. The visible light image and IR image when Δz = 0 are shown in Fig. 2(b) and 2(c) respectively. The IR image when Δz = 100 μm is shown in Fig. 2(d). The PSF recorded for Δz = 100 μm is shown in Fig. 2(e). The reconstruction results with LRA (150 iterations) (Supplementary Video 1), NLR (α = 0, β = 0.6) and LRRA (20 iterations, α = 0.2, β = 1) (Supplementary Video 2) are shown in Fig. 2(f–h), respectively. The recorded IR images for Δz = 0, 150 μm and 200 μm are shown in Fig. 2(i), 2(j)and 2(l), respectively, and the reconstructions using LRRA for Δz = 150 μm and 200 μm are shown in Fig. 2(k) and 2(m), respectively. The RMSE and entropy for LRA, NLR and LRRA are compared for the cross object and LRRA was found to perform significantly better than both LRA and NLR as shown in Fig. 2(n).

The experiment was then repeated on a bundle of silk fibers. The absorption spectrum of silk is shown in Fig. 3(a) and the chemical image (Δz = 0) is shown in Fig. 3(b). The reconstruction results of Fig. 3(b) using LRRA with PSFs recorded at Δz = 25 μm, 50 μm and 75 μm are shown in Fig. 3(c–e), respectively. The image directly recorded at Δz = 50 μm is shown in Fig. 3(f). Comparing Fig. 3(d) and 3(f), a good overlap can be observed. A single strand of silk fiber was imaged as shown in Fig. 3(g) whose varying blur indicates varying depth. It was reconstructed using a PSF recorded at Δz = 50 μm as shown in Fig. 3(h) and the corresponding direct image is shown in Fig. 3(i). A polymer test pattern, prepared on a CaF2 substrate, using a chrome-on-glass United States Air Force (USAF) 1954 target as a lithography mask was examined next, whose absorption spectrum is shown in Fig. 3(j). The direct imaging and reconstruction of object ‘3’ from the USAF target using LRRA are shown in Fig. 3(k) and 3(l), respectively. In some cases, it was necessary to adjust the magnification of the PSFs before reconstruction to achieve optimal reconstruction and for some cases, synthetic PSFs generated from recorded PSFs were used for reconstruction (see supplementary information Section 4).

Fig. 3. Reconstruction results of thick and strong IR absorbing materials. (a) Normalised absorbance of the silk fiber sample. The characteristic absorption peaks are noted at 3297 cm–1 and 1600 cm–1. The spectral matrices for Δz = 0 from 1708 cm–1 to 1558 cm–1 were averaged to obtain Fig. 3(b). (b) Direct imaging of a bundle of silk threads in 3D space, showing focused and out-of-focus objects. Reconstruction result using LRRA (α = 0.5, β = 1, iterations = 18) using PSFs recorded at (c) Δz = 25 μm, (d) Δz = 50 μm and (e), Δz =75 μm. (f) Direct imaging of a bundle of silk fibers for Δz ≈ 50 μm. The magnified regions of the direct imaging and reconstruction results indicate that the LRRA could obtain a good reconstruction even for thick objects and strong IR absorbing objects. (g) Direct imaging of a single silk fiber in 3D space. (h) Reconstruction results using PSF of Δz = 50 μm. (i) Direct imaging of the single silk fiber for Δz ≈ 50 μm. (j) Normalised absorbance of the polymer material of USAF target with characteristic peaks at 1500 cm–1 and 1238 cm–1. (k) Direct imaging of the USAF object ‘3’ located at Δz ≈ 50 μm. (l) Reconstruction result of the USAF object ‘3’ using LRRA (α = 0.6, β = 1, iterations = 14). The scale bar is 50 μm. The lateral resolution of the imaging system is ~15 μm.

下载图片 查看所有图片

3 Conclusion

In conclusion, we demonstrated 3D imaging without TBI using deterministic fields. Unlike scattering based SALCAD fields, where the PSFs have to be recorded at all axial planes, with deterministic SALCAD fields it is possible to calculate the PSFs from the complex amplitude of the optical modulator. Consequently, the reconstruction in this case mimics that of conventional incoherent holography. Furthermore, we have invented a new reconstruction method by introducing non-linearity into the LRA using the NLR approach, which performs significantly better than both LRA and NLR. The method and reconstruction have been demonstrated in the FTIRm chemical imaging system by introducing necessary data structure conversions (see Supplementary information Section 5). We believe that the approach can be easily adapted to power sensitive areas such as astronomical, fluorescence and biomedical imaging. Above all, we also believe that this research work will benefit the users of synchrotron based FTIRm technique to be able to investigate multiple samples during the limited beamtime in the future which was the motivation that led to this research work. The proposed method opens a new direction where deterministic fields and structured light can be engineered for rapid 3D imaging with a low photon budget, high lateral, and axial resolutions.

The evolution of incoherent holography over the years from the complicated architectures such as rotational shearing interferometer26, conoscopic holography27, and FINCH2, 3, to this version of holography with deterministic SALCAD fields using Lucy-Richardson-Rosen algorithm is interesting. We believe that the proposed method will compete with the existing coherent holography methods such as MIR digital holography and holographic interferometry using quantum cascade laser28, single shot Raman holography29, and label-free second harmonic phase imagers30.

This research was undertaken on the IRM beamline at the Australian Synchrotron (Victoria, Australia), part of ANSTO (Proposal ID. 15775, Reference No. AS1/IRM/15775 and Proposal ID. M17333, Reference No. AS2/IRM/17333). This work was performed in part at the Swinburne's Nanofabrication Facility (Nanolab). Funded by European Union’s Horizon 2020 research and innovation programme under grant agreement No. 857627 (CIPHR).

The authors declare no competing financial interests.

Supplementary information for this paper is available at https://doi.org/10.29026/oes.2022.210006


[1] Rosen J, Vijayakumar A, Kumar M, Rai MR, Kelner R et alRecent advances in self-interference incoherent digital holographyAdv Opt Photonics20191116610.1364/AOP.11.000001

[2] Rosen J, Brooker GDigital spatially incoherent Fresnel holographyOpt Lett20073291291410.1364/OL.32.000912

[3] Rosen J, Brooker GNon-scanning motionless fluorescence three-dimensional holographic microscopyNat Photonics2008219019510.1038/nphoton.2007.300

[4] Siegel N, Lupashin V, Storrie B, Brooker GHigh-magnification super-resolution FINCH microscopy using birefringent crystal lens interferometersNat Photonics20161080280810.1038/nphoton.2016.207

[5] Ables JGFourier transform photography: a new method for X-ray astronomyPubl Astron Soc Aust1968117217310.1017/S1323358000011292

[6] Dicke RHScatter-hole cameras for X-rays and gamma raysAstrophys J1968153L101L10610.1086/180230

[7] Wagadarikar A, John R, Willett R, Brady DSingle disperser design for coded aperture snapshot spectral imagingAppl Opt200847B44B5110.1364/AO.47.000B44

[8] Lee K, Park YExploiting the speckle-correlation scattering matrix for a compact reference-free holographic image sensorNat Commun201671335910.1038/ncomms13359

[9] Vijayakumar A, Rosen JInterferenceless coded aperture correlation holography–a new technique for recording incoherent digital holograms without two-wave interferenceOpt Express201725138831389610.1364/OE.25.013883

[10] Antipa N, Kuo G, Heckel R, Mildenhall B, Bostan E et alDiffuserCam: lensless single-exposure 3D imagingOptica201851910.1364/OPTICA.5.000001

[11] Park JH, Park J, Lee K, Park YDisordered optics: exploiting multiple light scattering and wavefront shaping for nonconventional optical elementsAdv Mater202032190345710.1002/adma.201903457

[12] Anand V, Ng SH, Maksimovic J, Linklater D, Katkus T et alSingle shot multispectral multidimensional imaging using chaotic wavesSci Rep2020101390210.1038/s41598-020-70849-7

[13] Nasse MJ, Walsh MJ, Mattson EC, Reininger R, Kajdacsy-Balla AHigh-resolution Fourier-transform infrared chemical imaging with multiple synchrotron beamsNat Methods2021841341610.1038/nmeth.1585

[14] Ong L, Pax AP, Ong A, Vongsvivut J, Tobin MJ et alThe effect of pH on the fat and protein within cream cheese and their influence on textural and rheological propertiesFood Chem202033212732710.1016/j.foodchem.2020.127327

[15] Osmond G, Boon JJ, Puskar L, Drennan JMetal stearate distributions in modern artists' oil paints: surface and cross-sectional investigation of reference paint films using conventional and synchrotron infrared microspectroscopyAppl Spectrosc2012661136114410.1366/12-06659

[16] Dainty JC. Stellar speckle interferometry. In Dainty JC. Laser Speckle and Related Phenomena 255–280 (Springer, 1975)

[17] Rai MR, Rosen JResolution-enhanced imaging using interferenceless coded aperture correlation holography with sparse point responseSci Rep202010503310.1038/s41598-020-61754-0

[18] Lohman AW. OpticalInformationProcessing (Physikalisches Institut der Universitat, 8520 Erlangen, EDR, 1978).

[19] Wu JC, Zhang H, Zhang WH, Jin GF, Cao LC et alSingle-shot lensless imaging with Fresnel zone aperture and incoherent illuminationLight Sci Appl202095310.1038/s41377-020-0289-9

[20] Wu JC, Cao LC, Barbastathis GDNN-FZA camera: a deep learning approach toward broadband FZA lensless imagingOpt Lett20214613013310.1364/OL.411228

[21] Horner J L, Gianino PDPhase-only matched filteringAppl Opt19842381281610.1364/AO.23.000812

[22] Rai MR, Vijayakumar A, Rosen JNon-linear adaptive three-dimensional imaging with interferenceless coded aperture correlation holography (I-COACH)Opt Express201826181431815410.1364/OE.26.018143

[23] Richardson WHBayesian-based iterative method of image restorationJ Opt Soc Am197262555910.1364/JOSA.62.000055

[24] Lucy LBAn iterative technique for the rectification of observed distributionsAstron J19747974575410.1086/111605

[25] Anand V, Ng SH, Katkus T, Maksimovic J, Klein AR et alExploiting spatio-spectral aberrations for rapid synchrotron infrared imagingJ Synchrotron Radiat2021281616161910.1107/S1600577521007104

[26] Murty MVRK, Hagerott ECRotational-shearing interferometryAppl Opt1966561561910.1364/AO.5.000615

[27] Sirat GYConoscopic holography. I. Basic principles and physical basisJ Opt Soc Am A19929708310.1364/JOSAA.9.000070

[28] Ravaro M, Locatelli M, Pugliese E, Di Leo I, de Cumis MS et alMid-infrared digital holography and holographic interferometry with a tunable quantum cascade laserOpt Lett2014394843484610.1364/OL.39.004843

[29] Shi KB, Li HF, Xu Q, Psaltis D, Liu ZWCoherent anti-stokes Raman holography for chemically selective single-shot nonscanning 3D imagingPhys Rev Lett201010409390210.1103/PhysRevLett.104.093902

[30] Shaffer E, Moratal C, Magistretti P, Marquet P, Depeursinge CLabel-free second-harmonic phase imaging of biological specimen by digital holographic microscopyOpt Lett2010354102410410.1364/OL.35.004102

Vijayakumar Anand, Molong Han, Jovan Maksimovic, Soon Hock Ng, Tomas Katkus, Annaleise Klein, Keith Bambery, Mark J. Tobin, Jitraporn Vongsvivut, Saulius Juodkazis. Single-shot mid-infrared incoherent holography using Lucy-Richardson-Rosen algorithm[J]. Opto-Electronic Science, 2022, 1(3): 210006.

引用该论文: TXT   |   EndNote



关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。