Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Appl Opt. Author manuscript; available in PMC 2011 February 23.
Published in final edited form as:
Appl Opt. 2011 January 10; 50(2): 170–176.
PMCID: PMC3043200

Analysis of diffracted image patterns from volume holographic imaging systems and applications to image processing

Diffracted image patterns from volume holograms that are used in volume holographic imaging systems (VHIS) are investigated. It is shown, that in VHIS, prior information about the shape and spectral properties of the diffracted patterns is important not only to determine the curvature and field of view of the image, but also for image registration and noise removal. A new methodology to study numerically and analytically the dependence of VHIS diffraction patterns with the hologram construction parameters and the read-out wavelength is described. Modeling and experimental results demonstrate that in most cases VHIS diffracted shapes can be accurately represented by hyperbolas.

1. Introduction

Highly selective spherical volume holograms are key elements in multimode multiplex spectroscopy [1-2] and 3D imaging applications [3-10]. For spectroscopy applications, volume holograms can be used as a spectral diversity filters (SDFs) [1] to map the spectrum of a diffuse source on a detector array using a spatially encoded pattern [2]. For the 3D imaging volume holograms are used to image simultaneously different depths of a sample on a CCD camera [3-10].

Volume holographic imaging systems (VHIS) can use a single wavelength laser or broadband sources such as LEDs. Using broadband source, a new technique denominated Spatial-Spectral VHIS (S2-VHIS) has the advantage of eliminating the need of any mechanical scanning [3-10]. Therefore, video rate images from different depths of a sample can be obtained.

A crescent shape, is a common shape produced by the transmitted or diffracted light from planar or curved volume hologram optical element (VHOE) under certain illumination conditions. It can be classified as a dark crescent when they are observed in the transmission beam or bright crescent when observed in the diffracted beam. For S2-VHIS the crescent shape is important since it provides information about the curvature of the field of view (FOV) and the spectral content of the image (4D imaging) [3-10]. Moreover, for in-vivo biomedical imaging, prior information about the crescent shape is essential for image registration and noise removal.

Previous work has utilized the Born approximation to predict the crescent shape [1,11]. However, this method assumes a weak interaction between the incident and diffracted light [7]. In the case of thick volume holograms, where the incident light is strongly depleted, another approach is needed in the optical design process.

Kogelnik's two-beam approximated coupled-wave analysis (ACWA) [12-13] and rigorous coupled-wave analysis (RCWA) [14-15] are well suited to analyze these strong and highly selective gratings. In fact, these have been used extensively to analytically model the diffraction effects of uniform gratings formed by two spatially planar waves. Among these methods, RCWA provides the best accuracy, at the expense of significantly increasing computation time. In case of thick volume holograms where one diffracted order is dominant, a modified ACWA which includes the effects of arbitrary conical incidence provides accurate and fast results. This capability coupled with short computation time allows to be implemented in ray-tracing packages such as Zemax® [10].

In this paper, we develop a model to quantify the properties of crescent shaped images formed in volume holographic imaging systems. This model, based on ACWA, provides more accurate results than models based on the Born approximation. In most of the cases, the crescent shapes used in VHIS are not circular or elliptical but hyperbolic. These results are confirmed with experiments. An important application of this result is to pre-determine the image field to set post-processing boundaries. We also describe this process in the paper.

The remainder of this paper is divided in five sections. In the second Section the basic operation of VHIS is described. The third section describes the ACWA and the relationship between the crescent curvature and the Bragg detuning condition. In Section fourth experimental validation is presented. The last Section is devoted to discussion and summary.

2. Basic operation of a S2-VHIS

The layout depicted in Figure 1 serves as an example of the S2-VHIS operation. A more detailed description of the S2-VHIS operation can be found in several references [3-10]. The basic configuration of a S2-VHIS consists of a 4F system with a holographic element placed in the Fourier plane. The holographic optical element is composed of thick angle-multiplexed planar and spherical wave gratings with high angular and spectral selectivity. For each plane in the object space, a mapping is produced to the image space by two properties of volume holograms [12-17] the spatial degeneracy and the spectral-angular dispersion. Due to the first property, the wavefront of a point source at any position along the y-axis (cylinders in Fig. 1) satisfies the Bragg phase matching condition of the VHOE and therefore is diffracted to the collector lens. This diffraction is responsible for the y-axis field of view. Current prototypes have achieved a y-field of view of ~1 mm using an objective lens with NA=0.55, focal length of 3.6 mm and relays lenses as indicated in [6-10]. This y-field of view value corresponds to a full width half maximum (FWHM) angular FOV of 8°.

Fig. 1
Basic layout of S2-VHIS showing the Objective lens (Lo) , the collector lens (Lc) and the hologram (HOE).

The FOV does not follow a straight line indicated by the axis of each cylinder, but a crescent shape as shown in Fig 3 and and4.4. The curvature of the diffractive pattern is not produced by the curvature of the hologram as could be misinterpreted from [1-3]. Both planar and curved hologram fringes result in curved diffraction shapes due to the spatial degeneracy property of the hologram. The diffractive shape, which resembles a crescent pattern, depends on the construction parameters of the hologram and the reconstruction wavelength. An accurate knowledge of this shape is required to correctly determine the space-spectral mapping of this type of imaging system.

Fig. 3
Measured diffractive patterns form a VHOE using two wavelengths. Curve fitting to a hyperbolic shape are indicated in dotted lines. Construction parameters of the VHOE indicated in Section 2
Fig. 4
Modeled diffractive traces compared with curve fitting of Fig. 1.

The spectral angular dispersion produces the field of view along the x-axis. This FOV is represented by several cylinders along the x-axis. Each cylinder has different gray level intensity which indicates that it is being imaged by a different wavelength. The FOV depends on the spectral bandwidth of the read-out source utilized. Current prototypes utilize LEDs with FWHM bandwidth ranging from 20 to 50 nm, producing FOV along the dispersive axis of approximately 150 μm which correspond to a FWHM angular FOV < 2.5°.

Multiplexing several gratings into the same volume can be used to map axial points at multiple depths in the object space to the image or camera plane as explained in [3-10].

3. ACWA Modeling Techniques for Volume Holograms

Modeling techniques to evaluate the performance volume holograms include the Born approximation, approximate coupled wave analysis [12,13] and rigorous couple wave analysis [13-14]. The Born approximation in Ref. [3-5,8-9,11], can provide simplified analytical models for volume holographic gratings. However, this method requires weak interaction between the incident beam and diffracted beam. This is not the case for the holograms used in the S2-VHIS.

The modified RCWA [14-15] provides the most accurate analysis for the diffraction phenomena for gratings with arbitrary reconstruction conditions. However, RCWA is a numerical technique that is computationally intensive and difficult to implement in component design algorithms.

An approximate coupled wave algorithm [12-13] can provide reasonable accuracy and fast estimation of grating diffraction efficiency. Moreover, the analytical expression of ACWA provides a more intuitive description of the diffraction phenomena than RCWA or the Born approximation, which is helpful in the investigation on the shape of the diffractive pattern.

3.1 Approximated coupled wave analysis: general model

In ACWA there are two fundamental parameters that describe the diffraction process: the strength and the direction of the diffracted beam. The direction of the diffracted vector inside the grating medium can be obtained from [12-13]

equation M1

where equation M2 are the incident and diffracted propagation vectors inside the medium. The z direction component is determined using,

equation M3

where equation M4, λ is the free-space wavelength, n2 is the refractive index of the medium, and the grating vector is given by,

equation M5

where λc is the grating construction wavelength, equation M6 and equation M7 are the construction reference and signal beam vectors inside the recording medium . The period of the grating can be computed using equation M8. Utilizing polar coordinate notation the diffracted vector is given by,

equation M9

where θ and ϕ are the zenith and azimuth angle in the holographic medium. In air, the zenith angle is obtained using the Snell's law equation M10.

For a transmission grating, the diffraction efficiency (DE) for the TE wave is given by [12-13,16],

equation M11

where ν is a parameter that determines the maximum diffraction efficiency of the grating for a specific reconstruction wavelength and angle, and ξ indicates the variation of the reconstruction parameters from the Bragg condition. For a lossless transmission grating, these parameter functions are given by,

equation M12
equation M13
equation M14

where Δn and tH are the index modulation and the thickness of the hologram respectively, [theta] is the detuning parameter and cr and cs are the obliquity factors described in [12].

3.2 ACWA: analytical model to predict the crescent shape

The ACWA model described in Section 3.1 can be used to numerically model the complete operation of the VHIS. However, in order to get a good insight of the dependence of crescent shape image properties with wavelength and hologram parameters an analytical model is required.

Here an analytical model from the ACWA equations was developed to provide an accurate description of these conical diffraction phenomena for the FOV used in VHIS. In the layout shown in Fig. 2, the trace produced by the beam diffracted by the HOE can be described using,

equation M15

where d is the distance between the HOE and the screen , θ’ , ϕ are the diffracted angles in air described in (4) and ,

equation M16
Fig. 2
Experiment layout used to measure crescent shape of the diffracted beam. The Z-axis is normal to the hologram and screen. The angles, θi and θ' , represent the incident and diffracted beam angles respectively. The distance from the hologram ...

For highly selective volume holograms such as the ones used in VHIS light from the object to image space propagate in well defined light paths that satisfy the Bragg condition (7). For the diffracted ray on Bragg condition, equation M17, the detuning equation is equal to zero when,

equation M18

For simplicity and without loss of generality it is assumed here that Ky=0. This can be achieved selecting a proper coordinate system. In this case, the K vector is given by,

equation M19

where equation M20. Using (4) and (12) in (11), the following expression is obtained,

equation M21


equation M22

Dividing (13) by cos θ,

equation M23

Multiplying previous equation by equation M24 and using (10) in (15),

equation M25

Squaring both sides of (16) and after some algebraic manipulation, the following equation is obtained,

equation M26

This equation relates the direction of the diffracted vector in terms of the azimuth and zenith angle in the medium with the hologram construction parameters and reconstructing wavelength (γ and q parameters). As described in Section 2 the S2-VHIS has a small angular FOV ( |θ'| < 10°). In this angular range, equation M27 with error < 2%.

Using this approximation and replacing (9) in (17) produces,

equation M28

Equation (18) describes the diffracted rays trajectories projected on the screen as shown in Fig. 2. Therefore the crescent shape image in terms of wavelength and hologram construction parameters (γ and q parameters) can be estimated.

In the Cartesian coordinate system, the graph of a quadratic equation in two variables (x,y) as shown in (18) is a conic section. The conic discriminant classification described in [18] can be utilized to determine the type of conic produced by the diffraction of the volume hologram. For the conic equation derived in (18), the possible curves are:

equation M29

4. Analysis of the crescent shapes and experimental validation

Previous studies of the diffractive patterns have shown somewhat ambiguous results in estimating the shape of the crescent. For example, the shape has been attributed to the curvature of the hologram fringes [1,3], or it has been assumed exclusively circular or elliptical [1,11]. In some cases this ambiguity can be attributed to the utilization of the paraxial approximation, the effect of other optical components in the imaging system, and the utilization of a reduced FOV.

In order to minimize errors, our initial modeling and experiments use a simple setup indicated in Fig. 2. This setup minimizes the effects of undesired aberrations caused by other optical elements other than the hologram. In this setup, the beams are co-axially aligned before passing through a microscope objective (NA=0.55) which then focuses the beams onto a hologram. This setup permits to fill completely the objective aperture producing a cone of light incident on the hologram larger than the one produced in the imaging system shown in Fig. 1. This allows for a greater portion of the angular bandwidth of the volume grating to be illuminated, therefore more of the resulting diffraction pattern can be observed. The expansion of the FOV minimizes the error in the estimation of the diffractive shape.

Fig. 3 shows bright crescents produced by two beams incident on a hologram using this setup. The light sources producing these images are the 488nm Argon laser and a 632.8 nm HeNe.

The best fit for the measured traces follows a curve given by hyperbolic equation given by

equation M30

For the HeNe wavelength the curve parameters were [A,B,C]=[ 1, 1.08,8900], whereas for the 488 nm beam the curve parameters were [A,B,C]=[1, 0.445, 3900] .

The curve parameters were obtained by minimizing the error between fitting and measured curves (both shown in Fig. 3). In order to simplify comparison, the fitting parameters are normalized by A. More details of this computation are shown in Appendix I.

To model the crescent of this VHOE, the ACWA (described in Section 2) was implemented as a custom dynamic link library (dll) in Zemax™. Results of our general numerical model are shown in Fig 4. In that figure the dotted lines represent the curve fitted to the hyperbolas obtained from previous measurements (Fig. 3). Comparison of modeled and measured traces shows very good accuracy (error < 3%).

However, the numerical model alone does not explain the causes for the obtained hyperbolic shapes. Neither, does it provide a good insight of the dependence of this diffractive pattern on the reconstruction wavelength and the hologram construction parameters. To pre-determine the image field and apply fast image processing techniques to VHIS it is important to understand and predict the diffractive image properties.

The analytical model described in Section 3.2 provides a criterion to determine the type of conic produced by the diffraction based on the hologram construction parameters and reconstruction wavelength.

In this model, the degree of tilt of the grating fringes is given by γ as shown in (12). For the gratings utilized in S2-VHIS, it is desirable to kept γ>80°. The use of smaller angles increases the tilt of the fringes and affects negatively the resolution and FOV of the imaging systems.

Replacing γ=80° in (13) it is found that hyperbolic shapes are produced when,

equation M31

Also, the diffracted angle in the medium has to be lower than the total internal reflection angle of the hologram medium q < sin 90 / n2. These limits indicate that for almost all of the operational spectral range of VHIS, the x2-term coefficient in (18) is positive. Therefore, using (19) it can be found that in the majority of the cases using S2-VHIS diffracted shapes are hyperbolas. The asymptotes slopes of the hyperbolas obtained from (18) are given by,

equation M32

For the gratings diffraction patterns shown in Fig 3 and and4,4, γ=90° the asymptotes slopes are 1.49 for 488 nm and 0.96 for 633 nm. Theoretical and experimental values agree with error< 3% .

In cases where the HOE has large tilt in the K vector, such as cases close to the 90-degree geometry VHIS, γ is very small and it is possible to obtain other shapes such as a circle, ellipse or parabola. However, those grating construction geometries are not very useful for low aberration and high resolution VHIS.

For simplicity, the equations derived in this section were obtained assuming a camera aligned normal to the HOE (Fig 2). However, we have derived the equation for other cases such as the one shown in [11] in which the camera is aligned normal to the in-plane diffracted vector (zero components in the y-axis). In general, this case involves a rotation along the y axis to match the normal of the camera surface to the in-plane diffracted vector. Although this rotation affects the way the curvature of the diffracted pattern is projected in the camera, the discriminant classification shown in (19) still holds for the hyperbolic shape produced by holograms used in S2-VHIS . In general, for a given curve defined by f(x,y), the different orientation of the camera will produce a curve given by f(ax+b,y), where a and b are parameters related with the coordinate system rotation [19] . This means that an ellipse may be transformed in a circle or vice versa. However, there is not a transformation from the hyperbola to an ellipse or circle by simply rotating the camera around the y-axis.

The importance of pre-determining the crescent shape information in S2-VHIS can be appreciated in Fig. 5. This figure shows the image of a tissue sample at two different depths separated by ~ 50 μm. The images were obtained using an LED source with center wavelength of 633 nm and a two-grating multiplexed hologram as described in [7]. In part a) of this figure two depth images are shown without segmentation or any other image processing technique. Part b) shows the segmented images using the hyperbolic shape pattern and spatial filtering optimized for each depth image. The segmentation is necessary to determine the boundaries of each depth section, and to select the best filter for each image. A prior-knowledge of the segmentation area and shape significantly reduce the image processing time. This also improves the SNR since spurious noise out of the segmented regions is avoided.

Fig. 5
Sample images: A) Before segmentation and without any imaging processing. B) After segmentation and spatial filtering. In both cases, images from holograms with curved fringes (z=-50μm) are on the left side.

For 4D imaging, each hyperbolic line can be related to a specific wavelength of the illumination source. Therefore, an accurate knowledge of the diffractive shape at each position of the FOV allows partial recovery of spectral information.

It can be seen in this figure, that both planar and curved HOE produced hyperbolic traces. The hyperbolic shape is related to the construction parameters of the grating as described in the preceding sections of this paper. An example showing the derivation of the hyperbolic traces fitting parameters from the analytical model is shown in Appendix I.

5. Conclusion

Light focused on a planar or curved fringe hologram diffract light into crescent shaped images. Prior knowledge of the crescent shape allows selecting the image regions for post processing to improve the image contrast. Previous studies used the Born approximation to predict the shape and width of the crescent image. However, this model is not accurate for the gratings strength parameters and selective characteristics of gratings used in volume holographic imaging systems.

In this paper, we have studied the diffractive patterns of a hologram using approximate couple wave analysis. It has been shown that ACWA accurately predict the diffracted beam shapes for planar and spherical gratings to within 3%. The model also explains the image shape dependence on wavelength and holographic grating construction parameters. Results of our analysis demonstrate that for the gratings used in S2-VHIS, the diffracted shape is a hyperbola and not a circles or ellipse as previously assumed.

The techniques shown in this paper to pre-determine the diffraction pattern, greatly decrease the time it takes to segment and process the acquired images as described in previous section.


The authors gratefully acknowledge the support from the National Institutes of Health (NIH-RO1CA134424).

Appendix I

I.1 Example showing the derivation of the fitting parameters from the analytical model

For an unslanted planar grating (γ=90°), the K vector can be obtained using (3),

equation M33

where θi is the construction angle in air and λc the construction wavelength For this example, using θi =34 ° and λc =488 nm, equation M34. The period of the grating is given by, equation M35 Using a reconstruction wavelength of λ=632 nm, a medium with refractive index n2=1.5 the q parameter is computed using (14),

equation M36

For a setup configuration similar to the one shown in Fig. 2, the distance between the hologram and the screen is d~91 mm. Using (18) the coefficients of the conic curve are obtained as follows:

equation M37

These results are divided by A. The normalized coefficients that are then

equation M38

After the operation, the parameter C is proportional to d2. For the utilized setup C= 9072 mm2.

The normalized A and B parameters depends only on the grating properties and not in the distance between the hologram and screen. The slopes of the hyperbola asymptotes can be obtained from these normalized parameters (equation M39). The slope facilitates comparison between the diffractive shapes produced at different wavelengths.

Using the same procedure for another reconstruction wavelength ( λ=488 nm) the values of [A,B,C] are [1, 0.454, 3742 mm2].

Comparison between the analytical model (eq. 18) and the numerical implementation of the full ACWA model (described in Section 3.1) produce errors < 0.5% for gratings with high angular selectivity (< 0.04°) which are commonly used in VHIS.


OCIS codes: (090.4220) Multiplexed hologram; (090.7330) Volume hologram; (110.0110) Imaging systems; (090.2890) Holographic optical elements.


1. Momtahan Omid, Hsieh Chao Ray, Karbaschi Arash, Adibi Ali, Sullivan Michael E., Brady David J. Spherical beam volume holograms for spectroscopic applications: modeling and implementation. App. Opt. 2004;43:6557–6567. [PubMed]
2. Hsieh Chaoray, Momtahan Omid, Karbaschi Arash, Adibi Ali. Compact Fourier-transform volume holographic spectrometer for diffuse source spectroscopy. Optics Letters. 2005 April 15;30(8):836–838. [PubMed]
3. Sinha Arnab, Sun Wenyang, Shih Tina, Barbastathis George. Volume holographic Imaging in Transmission Geometry. Applied Optics. 2004;43(7):1533–1551. [PubMed]
4. Li Z, PSaltis D, Liu W, Johnson WR, Bearman G. Volume Holographic Spectral Imaging. Proc. SPIE. 2005;5694:33.
5. Liu W, Psaltis D, Barbastathis G. Real-time spectral imaging in three spatial dimensions. Opt. Lett. 2002;27:854–856. [PubMed]
6. Luo Y, Gelsinger PJ, Barbastathis G, Barton JK, Kostuk RK. Optimization of multiplexed holographic gratings in PQ-PMMA for spectral-spatial filters. Opt. Lett. 2008;33:566–568. [PubMed]
7. Gelsinger-Austin PJ, Luo Y, Watson JM, Kostuk RK, Barbastathis G, Barton JK, Castro JM. Optical design for a spatial-spectral volume holographic imaging system. Opt. Eng. 2010;49 [PMC free article] [PubMed]
8. Sinha Arnab, Barbastathis George. Broadband volume holographic imaging. Applied Optics. 2004;43(27):5214–5221. [PubMed]
9. Wissmann P, Baek Oh S, Barbastathis G. Simulation and optimization of volume holographic imaging systems in Zemax. Optics Express. 2008;16(10):7516–7524. [PubMed]
10. Lou Yuan, Castro JM, Barton Jennifer, Kostuk Raymond K., Barbastathis G. Simulation and Experiment of non-uniform multiplexed gratings in volume holographic imaging systems. Optics Express. 2010;18(18):19273–19285. [PMC free article] [PubMed]
11. Oh Se Baek, Watson Jonathan M., Barbastathis1 George. Theoretical analysis of curved Bragg diffraction images from plane reference volume holograms. Applied Optics. 2009;48(31) 1,* 1,†. 2. [PubMed]
12. Kogelnik H. Coupled wave theory for thick hologram gratings. Bell Syst. Tech. J. 1969;48:2909–2947.
13. Solymar L, Cooke DJ. Volume Holography and Volume Gratings. Academic Press; 1981.
14. Moharam MG, Gaylord TK. Three-dimensional vector coupled-wave analysis of planar-grating diffraction. J. Opt. Soc Am. OSA. 1983;73(9):1105–1112.
15. Moharam MG, Gaylord TK. Rigorous coupled-wave analysis of planar-grating diffraction. J. Opt. Soc. Am. 1981;71:811–818.
16. Goodman J. Introduction to Fourier Optics. McGraw-Hill; 1996.
17. K Maeda W. Thesis. The University of Arizona, ECE Department; 2005. Edge-illumination gratings in PQ-doped PMMA for OCDMA applications.