Chinese Optics Letters, 2016, 14 (6): 060902, Published Online: Aug. 3, 2018  

Fast reconstruction of multiple off-axis holograms based on a combination of complex encoding and digital spatial multiplexing Download: 845次

Author Affiliations
1 College of Physics and Electronics, Shandong Normal University, Jinan 250014, China
2 College of Physics and Electronic Engineering, Qilu Normal University, Jinan 250200, China
Abstract
An algorithm is proposed for the fast reconstruction of off-axis digital holograms based on a combination of complex encoding (CE) and spatial multiplexing (SM). In this algorithm, every two off-axis holograms recorded in sequence are first assembled into a CE hologram using the CE method, and then four of the CE holograms are again encoded into one complex spatial multiplexing (CSM) hologram based on the SM algorithm. It is demonstrated that the eight holograms encoded into such one CSM hologram can be quickly reconstructed by performing a two-dimensional (2D) Fourier transform (FT) on the CSM hologram. Using this method, the eight 2D FTs required for the reconstruction of the eight holograms in the conventional spatial filtering methods can be simplified to a process with only one 2D FT, which can largely improve the computation efficiency with the resolution of the reconstructed images nearly unchanged.

Digital holography (DH) has drawn much attention in recent years due to its unique capability for digitally retrieving the complex amplitude distribution transmitted or scattered by a sample. Its applications have been extended to many areas, such as contour measurement[15" target="_self" style="display: inline;">5], digital microscopy[69" target="_self" style="display: inline;">9], image recognition[1012" target="_self" style="display: inline;">–12], x ray imaging[1315" target="_self" style="display: inline;">–15], electron microscopy[1618" target="_self" style="display: inline;">–18], and various noninvasive measurements[1921" target="_self" style="display: inline;">–21]. According to the relative orientations between the object beam and the reference beam in a holographic recording system, DHs are often divided into two categories: on-axis and off-axis. The former is capable of effectively utilizing the spatial frequency bandwidth of image sensors and thus can capture finer sample spatial details, but it suffers from the disturbance of the zero-order autocorrelation and the conjugate image, which need to be eliminated by using phase-shifting techniques or other time-consuming elimination algorithms[2226" target="_self" style="display: inline;">26]. In off-axis DHs, the autocorrelation and the conjugate items are separated from the required image item in the spatial frequency domain and thus can be eliminated simply by digital spatial filtering[2730" target="_self" style="display: inline;">30], although the off-axis geometry will result in an extra demand for the bandwidth of imaging sensors.

In the reconstruction of the recorded wavefront from an off-axis hologram based on spatial filtering, at least a two-dimensional (2D) Fourier transform (FT) is required, which is usually a time-consuming operation, especially for a hologram with a large number of pixels. Although such a reconstruction process could be speeded up to video rate using parallel computing techniques, some special graphic processing units and high programming skill are required. Recently, Girshovitz et al. proposed a novel and effective method for speeding up the reconstruction of off-axis digital holograms based on digital angular multiplexing (AM)[31] and a complex encoding (CE)[32] algorithm, which realizes the quick reconstruction of up to four holograms at the same time. This approach greatly increases the reconstruction efficiency and makes it possible to realize a real-time holographic reconstruction with a personal computer. In this paper, we present a complex spatial multiplexing (CSM) method for the fast reconstruction of off-axis holograms based on a combination of the CE algorithm and the spatial multiplexing (SM) algorithm[33]. Using this method, up to eight holograms can be encoded into one complex hologram and synchronously reconstructed, and so the reconstruction efficiency can be further improved.

Suppose that the interference intensity of an off-axis digital hologram recorded by an image sensor is Ht(x,y)=R2+|Ot(x,y)|2+Ot(x,y)Rexp(i2πξ0x)+Ot*(x,y)Rexp(i2πξ0x),where Ot(x,y) is the complex amplitude of the sample recorded at time t, and R is the amplitude of the reference beam. For the sake of simplicity, we assume that the reference beam is a plane wave tilted in the x direction with a carrier frequency of ξ0. If the bandwidth of the recorded object wave is equal to 2ξm, the size of ξ0 is usually set to one and a half times the bandwidth for an ideal off-axis hologram, that is, ξ0=3ξm. For extracting the sample’s complex amplitude from such an off-axis hologram, it is necessary to first transform Eq. (1) into its spatial frequency domain by a 2D FT. The spatial frequency of the hologram Ht can be expressed as H˜t=F{Ht}=F{R2+|Ot|2}+RF{Ot}δ(ξ+ξ0)+RF{Ot*}δ(ξξ0)=Ct(ξ,η)+RO˜t(ξ,η)δ(ξ+3ξm)+RO˜t*(ξ,η)δ(ξ3ξm),where F{ξ,η} denotes the FT operation, is the convolution operator, Ct(ξ,η) is the autocorrelation item determined by the first two items in Eq. (1), O˜t(ξ,η) and O˜t*(ξ,η) are the FTs of Ot(x,y) and Ot*(x,y). The delta functions of δ(ξ+ξ0) and δ(ξξ0) mean the frequency shifts of O˜t(ξ,η) and O˜t*(ξ,η), respectively, introduced by the carrier frequency ξ0 of the reference beam. Figure 1(a) shows a sketch of the spatial spectrum of the hologram Ht, in which ξ0 is set to three times the maximum frequency of the recorded sample ξm. It can be seen that, in this situation, the spatial spectrum of the sample could be directly extracted from the spatial spectrum of the hologram by spatial filtering[2730" target="_self" style="display: inline;">30], and then the recorded complex amplitude of the sample wave can be reconstructed by a further inverse FT. That is to say, two 2D FTs are needed for the retrieval of the complex amplitude of the sample wave from each hologram based on the conventional spatial filtering methods. However, 2D FTs are often time-consuming operations, especially for a hologram with a large number of pixels.

Fig. 1. Sketch of the spatial spectrum of (a) a conventional ideal off-axis hologram and (b) a CSM hologram encoded by eight off-axis holograms.

下载图片 查看所有图片

To improve the reconstruction efficiency of the off-axis holograms, we propose the CSM algorithm based on a combination of the CE algorithm and the SM algorithm. In this CSM algorithm, we first assemble eight off-axis holograms Hi(x,y)(i=1,2,,8), recorded in sequence into four CE functions as follows: {H1,2=H1+jH2H3,4=H3+jH4H5,6=H5+jH6H7,8=H7+jH8,where the sign j represents the unit’s imaginary number; then, these four CE functions are further synthesized into a CSM function. This function is constructed by adding four complex holograms, each of which is multiplied by a plane wave tilted in the y direction (orthogonal to the tilt direction of the reference beam in a holographic recording). The CSM function can be expressed as H=H1,2exp(j6πξmy)+H3,4exp(j2πξmy)+H5,6exp(j2πξmy)+H7,8exp(j6πξmy).

Thus, all the information of the eight holograms is encoded in one CSM function. Next, we will prove that the object waves recorded in the eight holograms can be retrieved from the spatial spectrum of the CSM function, which means only one 2D FT is required for the reconstruction of the eight holograms.

By substituting Eqs. (1) and (3) into Eq. (4), the CSM function can be further written as H=A1exp[j6πξmy]+S1exp[j2π(ξ0x+3ξmy)]+S2exp[j2π(ξ0x3ξmy)]+A2exp[j2πξmy]+S3exp[j2π(ξ0x+ξmy)]+S4exp[j2π(ξ0xξmy)]+A3exp[j2πξmy]+S5exp[j2π(ξ0xξmy)]+S6exp[j2π(ξ0x+ξmy)]+A4exp[j6πξmy]+S7exp[j2π(ξ0x3ξmy)]+S8exp[j2π(ξ0x+3ξmy)],where Ai(i=1,2,3,4) correspond to the autocorrelation items of the original holograms. Si(i=1,2,,8) are determined by the following formula: {S1=O1+jO2,S2=O1*+jO2*S3=O3+jO4,S4=O3*+jO4*S5=O5+jO6,S6=O5*+jO6*S7=O7+jO8,S8=O7*+jO8*.

Obviously, for retrieving the object waves of the encoded holograms, all the components Si(i=1,2,,8) expressed in Eq. (6) should be firstly separated from the CSM function, which could be realized by spatial filtering in the spatial frequency domain.

The spatial frequency of the CSM function is just equal to the FT of Eq. (5), which can be expressed as H˜=F{H}=[A˜1+RS˜1δ(ξ+ξ0)+RS˜2δ(ξξ0)]δ(η+3ξm)+[A˜2+RS˜3δ(ξ+ξ0)+RS˜4δ(ξξ0)]δ(η+ξm)+[A˜3+RS˜5δ(ξ+ξ0)+RS˜6δ(ξξ0)]δ(ηξm)+[A˜4+RS˜7δ(ξ+ξ0)+RS˜8δ(ξξ0)]δ(η3ξm),where S˜i=F{Si} (I=1,2,8). Figure 1(b) gives a sketch of the spatial spectrum of the CSM function. It can be seen that the spatial spectra S˜i of all the components Si are spatially separated in the spatial frequency domain, and so they can be extracted by spatial filtering. The spatial spectra O˜i of the object waves recorded by the eight holograms can be retrieved according to the following formula: {O˜1(ξ,η)=12[S˜1(ξ,η)+S˜2(ξ,η)*]O˜2(ξ,η)=12j[S˜1(ξ,η)S˜2(ξ,η)*]O˜3(ξ,η)=12[S˜3(ξ,η)+S˜4(ξ,η)*]O˜4(ξ,η)=12j[S˜3(ξ,η)S˜4(ξ,η)*]O˜5(ξ,η)=12[S˜5(ξ,η)+S˜6(ξ,η)*]O˜6(ξ,η)=12j[S˜5(ξ,η)S˜6(ξ,η)*]O˜7(ξ,η)=12[S˜7(ξ,η)+S˜8(ξ,η)*]O˜8(ξ,η)=12j[S˜7(ξ,η)S˜8(ξ,η)*],where S˜i(ξ,η)*=[S˜i(ξ,η)]*. Finally, the complex amplitude Oi of the object waves can be reconstructed by an inverse FT of the corresponding spatial spectrum O˜i.

Figure 2 gives the flowchart of the CSM algorithm described above, which mainly contains the following steps: (A) encode eight holograms Hi (i=1,2,,8) into one CSM function H using the CSM encoding algorithm shown in Eqs. (3) and (4); (B) get the spatial spectrum of the CSM function by calculating a 2D fast FT; (C) extract the components S˜i from the spatial spectrum of the CSM function as shown in Fig. 1(b); (D) retrieve the spatial spectra O˜i of the object waves from the data S˜i based on Eq. (8); (E) if required, take further steps such as inversely transforming O˜i to get the complex amplitudes of the object wave Oi, unwrapping the phase distributions, and enlarging them for display.

Fig. 2. Flowchart of the CSM algorithm for fast reconstruction of off-axis holograms.

下载图片 查看所有图片

Experiments were carried out to demonstrate the method described above. The holographic recording setup used in the experiments is shown in Fig. 3, which is essentially a Mach–Zehnder interferometer. In our experiments, the light source is an He–Ne laser with a wavelength of 632.8 nm. The recorded object is a dynamic refractive index field produced by a flame, which is imaged on the recording plane by an imaging lens and interferes with a tilted plane wave. The holograms were recorded by a CCD image sensor with a pixel size of 6.45μm×6.45μm and pixel number of 1024×1024.

Fig. 3. Off-axis holographic recording setup used in experiments for CSM algorithm. S, source; BE, beam expander; BS, beam splitter; M, mirror; L, lens; CCD, image sensor.

下载图片 查看所有图片

A series of off-axis holograms was recorded in sequence, one of which is shown in Fig. 4(a). For the reconstruction of the holograms using our CSM algorithm described above, we first divided them into groups: each group contains eight holograms. Then we encode the eight holograms in each group into a CSM function according to Eqs. (3) and (4) and calculate the spatial spectrum of the CSM function. Figure 4(b) shows an example of the spatial spectrum of a CSM function of the recorded holograms in experiments. Using the eight parts marked by the red circles in Fig. 4(b), the frequency spectra of the complex amplitude distributions of the object waves recorded in the eight holograms can be synchronously retrieved according to Eq. (8). Next, we can obtain the complex amplitude distributions by carrying the inverse FT operation on each of the eight retrieved frequency spectra. Figures 5(a)5(h) give the phase distributions of the object waves reconstructed from a group of holograms using our CSM algorithm.

Fig. 4. (a) Example of the off-axis holograms recorded in the experiment. (b) Example of the spatial spectrum of a CSM function encoded by eight recorded holograms.

下载图片 查看所有图片

Fig. 5. One group of the reconstructed wrapped phases using the CSM algorithm.

下载图片 查看所有图片

To test the phase reconstruction quality of the CSM algorithm, we compared the phase profile retrieved by using the CSM algorithm with that retrieved using the conventional spatial filtering algorithm (below we abbreviate it to CV algorithm) from the same hologram. In the CSM algorithm, eight different holograms are synchronously reconstructed by one 2D FT, while in the CV algorithm, the holograms are separately reconstructed one by one. One of the compared results is shown in Fig. 6, which shows a phase profile of the retrieved object wave as shown in Fig. 5(a) (the part along the red line). In Fig. 6, the blue line expresses the phase value retrieved based on the CSM algorithm and the red dashed line corresponds to the result using the CV method. It can be seen that the phase reconstructed by the CSM algorithm is consistent with that reconstructed by the conventional algorithm. Meanwhile, the computation efficiency of the CSM algorithm is greatly improved.

Fig. 6. Example of the compared results of the retrieved phases using the CSM algorithm (the blue line) with that using the CV method (the red dashed line) at the same position of the retrieved phase distributions from the same hologram as the part marked by the red line in Fig. 5(a).

下载图片 查看所有图片

To test the influence of the CSM algorithm on the resolution of the reconstructed images when eight holograms are synchronously reconstructed, we changed the phase object to a transmittance USAF resolution target (with group numbers of 2 to 7) and repeated the experiment. In this experiment, the lens L shown in Fig. 3 was also removed, and the distance between the target and the recording plane was set to about 175 mm. Figure 7(a) gives an example of the images reconstructed based on the CSM algorithm. As a comparison, Fig. 7(b) shows the image reconstructed from the same hologram but based on the CV algorithm. It can be seen that the image reconstructed by the CSM algorithm has the same resolution as that reconstructed using the CV algorithm.

Fig. 7. Examples of the images reconstructed by using (a) the CSM algorithm and (b) the CV algorithm when the object is a transmittance USAF resolution target with group numbers of 2 to 7.

下载图片 查看所有图片

Last, we compared the computation efficiency of the CSM algorithm with that of the CV algorithm, the AM algorithm[31], the complex angular multiplexing (CAM) algorithm[32], and the SM algorithm[33] under the same hardware and software environments. The compared processing times were firstly taken from encoding the holograms to getting the complex amplitudes of the object waves recorded in the holograms and then divided by the number of holograms. In order to compare the influence of the hologram size on the computation efficiency, holograms of different sizes were simulated in computer. Table 1 shows the average ratio of the CSM’s processing time (tCSM) to the CV’s processing time (tCV), the AM’s processing time (tAM), the CAM’s processing time (tCAM), and the SM’s processing time (tSM). The sizes of the holograms are set to 512×512, 1024×1024, and 2048×2048. All the calculations involved in Table 1 was performed in a single-core personal computer (Pentium4 2.8 GHz CPU, 2 GB RAM).

Table 1. Ratio of the CSM’s Processing Time to Other Algorithms’ for Reconstructing a Single Off-Axis Hologram

SizetCSM:tCV (%)tCSM:tAM (%)tCSM:tCAM (%)tCSM:tSM (%)
512×51210.739.144.977.0
1024×102411.942.548.577.4
2048×204812.140.845.775.8

查看所有表

It can be seen from Table 1 that the computation efficiency of our CSM algorithm is superior to that of the other four algorithms. For example, the CSM algorithm on average takes only about twelve percent of the processing time cost by the CV algorithm for reconstructing one hologram and less than fifty percent of the CAM’s processing time.

In conclusion, we demonstrate that the off-axis holograms can be quickly reconstructed by the CSM algorithm, in which only one 2D FT is required for getting the spatial spectra of eight holograms by encoding them into a CSM function based on a combination of the CE algorithm and SM algorithm; thus, a great deal of unnecessary calculations can be cut out and the computational efficiency is greatly improved. We think that this improved algorithm may be useful in the implementation of real-time holographic visualizations in a simple personal computer system. In particular, it can be a good choice in the design of a preview system for dynamic wavefront imaging in digital holographic microscopy. In addition, the CSM algorithm can also provide an approach for compressing off-axis holograms in optical information storage and processing.

References

[1] SolísS. M.Hernández-MontesM. S.SantoyoF. M., Biomed. Opt. Express3, 3203 (2012).BOEICL2156-7085

[2] AbdelsalamD. G.MinJ.KimD.YaoB., Chin Opt. Lett.13, 100701 (2015).

[3] NehmetallahG.BanerjeeP. P.FerreeD.KephartR.PraharajS., Chin Opt. Lett.9, 120004 (2011).

[4] AbdelsalamD. G.MinJ.KimD.YaoB., Chin Opt. Lett.13, 100701 (2015).

[5] WangZ.ChenY.JiangZ., Chin Opt. Lett.14, 010008 (2016).

[6] TrujilloC. A.Garcia-SucerquiaJ., Opt. Lett.39, 2569 (2014).OPLEDP0146-9592

[7] Sánchez-OrtigaE.DoblasA.SaavedraG.Martínez-CorralM.Garcia-SucerquiaJ., Appl. Opt.53, 2058 (2014).APOPAI0003-6935

[8] ClarkD. C.KimM. K., Chin Opt. Lett.9, 120001 (2011).

[9] MaJ.YuanC.SituG.PedriniG.OstenW., Chin Opt. Lett.11, 090901 (2013).

[10] AnandA.JavidiB., Chin. Opt. Lett.12, 060012 (2014).CJOEE31671-7694

[11] SuC.XiaX.LiH.LiuX.KuangC.XiaJ.WangB., Chin. Opt. Lett.12, 060007 (2014).CJOEE31671-7694

[12] SeifiM.DenisL.FournierC., J. Opt. Soc. Am. A30, 2216 (2013).JOAOD60740-3232

[13] SchlotterW. F.LüningJ.RickR.ChenK.ScherzA.EisebittS.GüntherC. M.EberhardtW.HellwigO.StöhrJ., Opt. Lett.32, 3110 (2007).OPLEDP0146-9592

[14] GuehrsE.GuntherC. M.KonneckeR.PfauB.EisebittS., Opt. Express17, 6710 (2009).OPEXFF1094-4087

[15] PellicciaD.RigonL.ArfelliF.MenkR. H.BukreevaI.CedolaA., Opt. Express21, 19401 (2013).OPEXFF1094-4087

[16] GanZ.PereaD. E.YooJ.PicrauxS. T.SmithD. J.McCartneyM. R., Appl. Phys. Lett.103, 153108 (2013).APPLAB0003-6951

[17] PantzerA.VakahyA.EliyahouZ.LeviG.HorvitzD.KohnA., Ultramicroscopy138, 36 (2014).ULTRD60304-3991

[18] RöderF.LubkA.WolfD.NiermannT., Ultramicroscopy144, 32 (2014).ULTRD60304-3991

[19] WangY.WangD.ZhaoJ.YangY.XiaoX.CuiH., Chin Opt. Lett.9, 030901 (2011).

[20] LuZ.SunJ.ZhouY.ZhangN.SunZ.LiuL., Chin Opt. Lett.13, 020901 (2015).

[21] YaoL.WuX.YangJ.WuY.GaoX.ChenL.GréhanG.CenK., Chin Opt. Lett.13, 072801 (2015).

[22] XuX. F.CaiL.WangY.MengX.SunW. J.ZhangH.ChengX.DongG.ShenX., Opt. Lett.33, 776 (2008).OPLEDP0146-9592

[23] DengJ.WangH.ZhangD.ZhongL.FanJ.LuX., Opt. Lett.38, 1506 (2013).OPLEDP0146-9592

[24] GuoC. S.ShaB.XieY. Y.ZhangX. T., Opt. Lett.39, 813 (2014).OPLEDP0146-9592

[25] StoykovaE.KangH.ParkJ., Chin. Opt. Lett.12, 060013 (2014).CJOEE31671-7694

[26] LiY.XiaoW.PanF.RongL., Chin Opt. Lett.12, 020901 (2014).

[27] CucheE.MarquetP.DepeursingeC., Appl. Opt.39, 4070 (2000).APOPAI0003-6935

[28] VerrierN.AtlanM., Appl. Opt.50, H136 (2011).APOPAI0003-6935

[29] PhamH.DingH.SobhN.DoM.PatelS.PopescuG., Biomed. Opt. Express2, 1781 (2011).BOEICL2156-7085

[30] FrenklachI.GirshovitzP.ShakedN. T., Opt. Lett.39, 1525 (2014).OPLEDP0146-9592

[31] GirshovitzP.ShakedN. T., Opt. Lett.39, 2262 (2014).OPLEDP0146-9592

[32] GirshovitzP.ShakedN. T., Opt. Express23, 8773 (2015).OPEXFF1094-4087

[33] ShaB.LiuX.GeX. L.GuoC. S., Opt. Express22, 23066 (2014).OPEXFF1094-4087

Bei Sha, Yujie Lu, Yiyan Xie, Qingyang Yue, Chengshan Guo. Fast reconstruction of multiple off-axis holograms based on a combination of complex encoding and digital spatial multiplexing[J]. Chinese Optics Letters, 2016, 14(6): 060902.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!