Incoherent digital holographic spectral imaging with high accuracy of image pixel registration∗

2021-05-06 08:54FengYingMa马凤英XiWang王茜YuanZhuangBu卜远壮YongZhiTian田勇志YanliDu杜艳丽QiaoXiaGong弓巧侠CeyunZhuang庄策云JinhaiLi李金海andLeiLi李磊
Chinese Physics B 2021年4期
关键词:王茜李磊金海

Feng-Ying Ma(马凤英), Xi Wang(王茜), Yuan-Zhuang Bu(卜远壮),Yong-Zhi Tian(田勇志), Yanli Du(杜艳丽) , Qiao-Xia Gong(弓巧侠),Ceyun Zhuang(庄策云), Jinhai Li(李金海), and Lei Li(李磊),†

1 School of Physics and Microelectronics,Key Laboratory of Materials Physics of the Ministry of Education,Zhengzhou University,Zhengzhou 450001,China

2School of Mechanical Engineering,Zhengzhou University,Zhengzhou 450001,China

Keywords: incoherent digital holography, high-precision registration, spectral imaging, microspectral imaging

1. Introduction

Holography can produce truly three-dimensional images,and be viewed from multiple angles.[1,2]With the rapid development of computer technology and high-resolution image sensor, digital holography exhibit some advantages of high speed, real time, full field of view, and quantitative phase contrast imaging.[3,4]Limited by the requirements for high coherence light source and the unavoidable laser speckle noise, the reconstructed image of digital hologram is seriously ruined by the noise.[5–7]Incoherent digital holography is free from the dependence of traditional digital holography on coherent light sources, and expands its applications to fluorescence microscopy,[8]color holography,[9–14]and adaptive optics.[15,16]By now, many kinds of incoherent digital holography technologies, such as optical scanning holography,[17,18]triangular holography,[19,20]Michelson interferometer type holography,[21,22]and Mach–Zehnder interferometer-type holography[23]have been intensively investigated. Fresnel incoherent correlation holography(FINCH)[24]was first proposed by Joseph Rosen and Gary Brooker in 2007. In the FINCH system, a spatial light modulator (SLM) and a CCD are coaxially arranged to form an in-line incoherent interferometer. It has the advantages of neither time scanning nor space scanning, high resolution, and easy matching with existing mature optical systems.

Spectral imaging technology,which is known as a revolution in the development history of optical instruments,is a perfect combination of spectral technology and imaging technology. The main application field of spectral imaging is remote sensing detection.[25,26]In recent years, its application in microscopy has gradually become a research hotspot. However,there are two common problems in the existing spectral imaging technology:chromatic difference in magnification of spectral images caused by dispersion effect and incapability of extracting 3D information of the target. Chromatic difference in magnification leads to pixel registration errors in spectral images,and thus reducing spectral reconstruction accuracy.[27,28]Spectral imaging can only provide two-dimensional(2D)spatial information of targets in different wavebands,which cannot fully describe the overall characteristics of targets,resulting in low accuracy of 3D recognition.Since the recording and reconstruction of holograms are no longer dependent on coherent light sources, it is possible to combine incoherent digital holography and spectral imaging.

Meanwhile, extending spectral imaging to biomedical engineering has unparalleled great potential for researchers to obtain more information about organs, tissues, and even cells.[29,30]One of the commonly used optical methods is the spectral diagnostic technique, which can obtain the entire spectrum of a single tissue site within a specific wavelength region. This method is often referred to as the point measure method, but it cannot provide spatial information of samples.According to the dispersion theory of light,the focusing region will move accordingly with wavelength increasing. When the object distance is less than the recording distance,it will cause some single-band images to become blurred with the increase of wavelength. This problem may not be ignored in the application of microspectral imaging.

In this work, an incoherent digital holographic spectral imaging system based on liquid crystal tunable filter (LCTF)and SLM is built. A serious of double-lens phase masks,none of whose focal lengths changes with wavelength, is designed and made. For each wavelength of LCTF output,SLM calls the phase masks of the corresponding wavelength, and CCD records the spectral holograms. The spectral images obtained by this method have a constant magnification, which can achieve pixel-level image registration,restrain image registration errors,and improve spectral reconstruction accuracy.Moreover, by replacing specific phase mask loaded on the SLM, we avoid the technical conundrum caused by dispersion theory. The incoherent digital spectral holographic system based on LCTF not only improves the capability of microscopic imaging,but also records and extracts the phase information, amplitude information, and color information of the stained cells.

2. Principle of holographic spectral imaging

The SLM is space-division multiplexed by two diffractive lenses with focal lengths of fd1and fd2. Because of the diffraction effect,the effective focal length of the lens producing chromatic aberration is inversely proportional to the wavelength as shown in the formula f =r2/Nλ.Here r is the radius of diffraction lens and N is the Fresnel number.When λ0is the working wavelength of the mask and the imaging wavelength is λ,the effective focal lengths of the two lenses are λ0fd1/λ and λ0fd2/λ, respectively. Therefore, the reflection function of SLM is expressed as

where

is the reconstructed distance,with

MTis the lateral magnification,c.c. is the complex conjugate of the second term on the right.

Fig.1. Schematic diagram of incoherent digital spectral holography. LCTF:liquid crystal tunable filter; L:collimating lens; SLM:spatial light modulator;CCD:charge-coupled device.

In order to satisfy the optimum interference condition of two beams after being split by SLM,the distance zhbetween CCD and SLM needs to be adjusted continuously with the imaging wavelength. In other words, when the two cones of the spherical waves are perfectly overlapped with each other,the optimal contrast between interference fringes can be achieved.The zhcan be obtained from the triangular similarity relationship in Fig.1 as follows:[31]

The incoherent hologram of a general 3D object g(xs,ys,zs)is

That is to say,in order to accurately adjust the value of zh,the mechanical scanning device needs to be introduced into the system, which will bring some errors inevitably, and change the transverse magnification of the system as shown in the following formula:

The change of transverse magnification of spectral image will lead to image registration errors, resulting in the inaccurate extraction of relative spectral intensity.

In order to obtain 3D images with high spectral reconstruction accuracy,a high-precision spectral imaging technology and method based on FINCH is proposed. A series of dual lens phase masks with constant focal lengths of fd1and fd2at different wavelengths is designed and made. Using the programmable characteristics of SLM, for each wavelength of LCTF output, SLM calls the masks of the corresponding wavelength. Three holograms with different phase constants are recorded sequentially by CCD. Using three-step phase shift technology, three holograms are linearly superimposed to eliminate the zero-order and twin image,thus obtaining the complex value hologram HF(x,y),

Using this method to record the spectral holograms, we can realize that the recorded distance,spectral holograms,reconstruction distance, and spectral reconstruction images are independent of wavelength,and the reconstructed images have the same lateral magnification,which can avoid the mechanical error caused by space scanning and suppress the registration error in spectral image fusion.

Suppose that the spectral power distribution of the incoherent light source is I(λ), the spectral transmittance of the LCTF is T(λ), the spectral reflectance function of the object surface is Ri,j(λ), the spectral sensitivity function of CCD is S(λ),the spectral power distribution of each pixel in the reconstructed image is Ii,j(λ), the exposure coefficient of the camera remains unchanged during the recording of spectral holograms,then the reconstructed image will be obtained by using angular spectrum diffraction algorithm,so the pixel size of the reconstructed image is equal to that of CCD,the reflectance of reconstructed image can be expressed as follows:

The spectral response function of optical instrument is obtained from the instruction manual provided by the instrument manufacturer. Formula(8)is used to extract and calculate the spectral power distribution of pixels in the reconstructed image.

3. Experimental results

In this section, we take the dice, trinkets and biological samples as experimental objects,confirm the feasibility of this method, and realize the high precision microscopic spectral imaging of cells.

3.1. Spectral imaging optical path design based on FINCH

The schematic of the spectral imaging system based on FINCH is shown in Fig.2.

Fig.2. Schematic diagram of spectral imaging system based on LCTF and FINCH.Incoherent light source,xenon lamp;L1,L2,L3: converging lens;P:polarizer;BS:beam splitters.

Fig.3. Single[(a)–(c)focal length of 255 mm]and dual[(d)–(f),focal lengths of 245 mm and 255 mm,random pixel mixing]lens phase masks with 0 phase shift at selected wavelengths of [(a), (d)] 450 nm, [(b), (c)] 550 nm, and [(e), (f)] 650 nm; (g)–(i) and (j)–(l), enlargement of parts in green rectangular boxes shown in panels(a)–(c)and(d)–(f),respectively.

Spatially incoherent illumination with controllable wavelength required for the test, which is generated by an incoherent xenon lamp source(CEL-TCX250,250 W)and LCTF(VariSpec Liquid Crystal Tunable Filters, 400 nm–720 nm,FWHM 10 nm). The focal lengths of L1, L2, and L3 are 60 mm, 60 mm, and 250 mm, respectively. For the convenience of subsequent processing, only pixels in CCD (QIMAGING digital camera RETIGA 6000, pixel size σc=4.54µm,pixel array)are used. The SLM is a series of phase masks only, (Holoeye Pluto, 1920×1080 pixels, 8-µm pixel pitch). The polarization direction of the polarizer is the same as that of the SLM.The other specifications of the experiment are as follows: fd1=245 mm, fd2=255 mm,zh=250 mm,d=150 mm(which is the distance from L3 to the SLM).

Single and dual lens phase masks with 0 phase shift at selected wavelengths are shown in Fig.3. Figures 3(a)–3(c)show single lens phase masks with the same focal length of 255 mm at 450 nm, 550 nm, and 650 nm, respectively. As can be seen,when the focal length and aperture of the diffraction lens are kept unchanged,the concentric rings of the phase mask become more and more sparse with the increase of the wavelength. The same trend is shown in Figs. 3(d)–3(f), in which the dual lens masks using a random pixel mixing space division multiplexing have a focal length of 245 mm and 255 mm at 450 nm,550 nm,and 650 nm,respectively.

3.2. Imaging results

The spectral holograms of two white dices with black dots are recorded from 400 nm to 700 nm in intervals of 10 nm. Figures 4(a)–4(o) show the selected spectral images from 560 nm to 700 nm in intervals of 10 nm. Figure 4(p)shows the spectrum of the point on the dice shown in the inset,and the corresponding CIE1931 chromaticity coordinates(x=0.3552,y=0.3448)are shown in Fig.4(q). For comparison,a cell phone camera picture of the two dices is shown in the inset of Fig.4(q). It can be seen that the incoherent holographic spectral imaging system can accurately reconstruct the color information of the object.

Figure 5 shows color imaging results of a Poker ear stud at wavelengths of 672.8 nm, 562.8 nm, and 452.8 nm. Figures 5(a)–5(c) show the reconstructed images at 672.8 nm,562.8 nm,and 452.8 nm. The color composite image is shown in Fig.5(d). For comparison,a cell phone camera picture and the color spectral image obtained from commercial imaging spectrometer Hyperspec VNIR(400 nm–900 nm)are shown in Figs.5(e)and 5(f). The structural similarity index is 0.68524 measured by this method, and 0.66336 by the commercial spectrometer. The results show that the color reproducibility of the reconstructed image obtained by this method is better than that by the commercial imaging spectrometer.

Fig.4. Spectral imaging results of dice: (a)–(o)selected spectral images from 560 nm to 700 nm in intervals of 10 nm,(p)spectrum of point on dice shown in inset,(q)CIE1931 chromaticity coordinates of point,with inset showing cell phone camera picture of two dices.

Fig.5. Color imaging results of the Poker ear stud: (a)–(c) reconstructed images at wavelengths of 672.8 nm, 562.8 nm, and 452.8 nm; (d)color composite image from panels(a)–(c);(e)cell phone camera picture;(f)imaging result from commercial imaging spectrometer.

Fig.6. Composited color images of ear suds at different reconstruction distances: ((a),(b))reconstructed images in the best focal plane of red mask and Doraemon,respectively;(c),(d)and(e),(f)magnified parts in green rectangular boxes shown in panels(a)and(b),respectively.

Using the same method, we study the color threedimensional spectral imaging characteristics of the system.The tested objects are two ear studs with 2 cm apart. The spectral holograms are recorded from 400 nm to 700 nm, in an interval of 20 nm. The composited color images are shown in Fig.6. Figures 6(a)and 6(b)show the composited color reconstructed images in the best focal plane of the red mask and Doraemon, respectively. Figures 6(c)–6(f) are the magnified parts in rectangular boxes shown in Figs. 6(a) and 6(b). The results show that the FINCH system 3D spectral image with high accuracy of image pixel registration can be achieved by keeping the transverse magnification constant in the spectral imaging process.

3.3. Incoherent digital holographic microspectral imaging based on LCTF

Amicroscopic imaging system based on FINCH is built,a 20×,0.4-NA microscope objective with a working distance of 5.9 mm is placed in front of the polarizer. Figures 7(a)–7(c)show the holograms of the stem transected cells of stained woody dicotyledons with phase factors of 0◦,120◦,and 240◦at 632.8 nm. The complex amplitude and phase diagram are shown in Figs. 7(d) and 7(e), respectively. Figure 7(f) shows the spectrum of the green point on the cells indicated in the inset.

Figures 8(a) and 8(c) show the microscopic imaging results from FINCH microscope and Nikon microscope,respectively. Figure 8(b) shows the enlargement in the green box in Fig.8(a). The green dots in Figs. 8(b) and 8(c) represent the same cell. Although the performance of the spectral microscopy system based on incoherent digital holography needs to be further optimized, compared with commercial optical microscope,it has significant advantages in obtaining 3D spatial information and spectral information of biological samples.

Fig.7. (a)–(c)Holograms at 632.8 nm with 0◦,120◦,and 240◦phase factors;(d),(e)amplitude and phase of complex hologram;(f)spectrum of green point on cells shown in inset.

Fig.8. (a)Spectral fusion image of stem transversely cut cells stained with woody dicotyledonous plant;(b)enlargement of green box in panel(a);(c)image from NIKON microscope.

4. Conclusions and perspectives

In this paper, a high-precision spectral imaging technology and method based on FINCH is proposed and analyzed theoretically in detail. By designing the phase masks loaded on SLM,the system has a constant transverse magnification at different imaging wavelengths, which avoids the mechanical error caused by space scanning and suppresses the registration error in spectral image fusion. The results show that this method can not only obtains the 3D spatial information and spectral information of the object simultaneously,but also has high accuracy of spectral reconstruction and excellent color reproducibility. This method has the potential applications in imaging and detection in label-free biological samples.

Acknowledgment

The authors would like to thank the following collaborators: Ms. Wang X,Mr. Bu Y Z,Prof. Li L,Prof. Tian Y Z,Prof. Du Y L,Prof. Gong Q X,Mr. Zhuang C Y,Mr. Li J H,for their contributions to the present work.

猜你喜欢
王茜李磊金海
MAPS PRESERVING THE NORM OF THE POSITIVE SUM IN Lp SPACES*
青岛金海种苗有限公司
故纸情怀
Effect of nitrogen gas flow and growth temperature on extension of GaN layer on Si*
Effect of non-Markovianity on synchronization
王茜作品赏析
程彦鹏、王茜、刘文作品
百花绽放
金海
High-resolution boosted reconstruction of γ-ray spectra∗