OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 11, Iss. 6 — Mar. 24, 2003
  • pp: 572–578
« Show journal navigation

Reconstruction of a complex object from two in-line holograms

Yuxuan Zhang and Xinyi Zhang  »View Author Affiliations


Optics Express, Vol. 11, Issue 6, pp. 572-578 (2003)
http://dx.doi.org/10.1364/OE.11.000572


View Full Text Article

Acrobat PDF (437 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

A novel numerical algorithm is proposed to reconstruct a complex object from two Gabor in-line holograms. With this algorithm, both the real and imaginary parts of the complex amplitude of the wave front in the object field can be retrieved, and the “twin-image” noise is eliminated at the same time. Therefore, the complex refractive index of the object can be obtained without disturbance. Digital simulations are given to prove the effectiveness of this algorithm. Some practical experimental conditions are investigated by use of error estimation.

© 2003 Optical Society of America

1. Introduction

X-ray holography is a powerful technique for high-resolution imaging of biological, material, mineral, and fossil specimens. It has many advantages over other microscopy imaging techniques, such as optical microscopy and electron microscopy [1

1. D. Sayer, J. Kirz, R. Feder, D. M. Kim, and E. Spiller, “Potential operating region for ultrasoft X-ray microscopy of biological materials,” Science 196, 1339–1340 (1977). [CrossRef]

3

3. I. McNulty, “The future of X-ray holography,” Nucl. Instrum. Methods A 347, 170–176 (1994). [CrossRef]

]. The resolution of x-ray holography is much higher than that of optical microscopy because the wavelength of x rays is much shorter than that of visible light. It could be several tens of nanometers at present. The penetrating length of x rays is much longer than that of an electron beam, which means that a thicker specimen can be observed by x rays. Soft-x-ray holography has been successfully used to reconstruct many kinds of small biological samples with high resolution [4

4. C. Jacobsen, M. R. Howells, J. Kirz, and S. Rothman, “X-ray holographic microscopy using photoresists,” J. Opt. Soc. Am. A 7, 1847–1861 (1990). [CrossRef]

]. Hard x rays can be used to penetrate a bulk sample such as a piece of fossil or ore several millimeters in size. However, sometimes it is difficult to obtain a high-contrast x-ray transmission image because the difference in the absorption coefficients of the different elements inside the sample is quite small. Thus phase-contrast imaging has attracted more attention in recent years [5

5. P. Spanne, C. Raven, I. Snigireva, and A. Snigirev, “In-line holography and phase-contrast microtomography with high energy x-rays,” Phys. Med. Biol. 44(3), 741–749 (1999). [CrossRef] [PubMed]

7

7. K. Hirano and A. Momose, “Development of an X-ray interferometer for high-resolution phase-contrast X-ray imaging,” Jpn. J. Appl. Phys. 38, L1556–L1558 (1999). [CrossRef]

]. Phase contrast can be much more sensitive than absorption contrast; therefore, by use of phase contrast a much higher signal-to-noise ratio and decreased radiation dose can be reached at the same time. Holography is highly suitable for phase-contrast imaging, because the complex amplitude of the original object wave front can be reconstructed with the digital holography technique. Then both the absorption and the phase shift caused by the object can be retrieved from the complex amplitude, and thus both phase-contrast and absorption-contrast images can be obtained.

2. Theory

Let us consider the formation of an in-line hologram. If the object in plane (X,Y) has a transmittance of t(x,y) and phase shift ϕ(x,y), the complex amplitude of the object wave front Ũ(x,y) becomes

Ũ(x,y)=Ũ0(x,y)t(x,y)eiϕ(x,y)
=Ũ0(x,y)[1-a(x,y)]eiϕ(x,y),
(1)

where Ũ0(x,y) is the incident wave front and a(x,y) is the pure absorption of the object. Let us assume that Ũ0(x,y)=1, which indicates a plane-wave incident condition, and assume that a(x,y)≪1, | φ(x,y)|≪1. These assumptions can be easily satisfied because the complex refractive indices of all elements in high-energy x-ray regions are always near to 1. If the size of the object is not too large (less than 10 mm) and the object does not contain too many heavy atoms, the absorption and phase shift will be quite small. Then we have

Ũ(x,y)[1a(x,y)][1+iϕ(x,y)]
=1a(x,y)+iϕ(x,y)ia(x,y)ϕ(x,y)
1ã(x,y).
(2)

We will neglect the high-order term of a(x,y)ϕ(x,y) in Eq. (2) and in the following equations. Here we define a complex absorption factor ã(x,y)=a(x,y)-iφ(x,y). Thus Re[ã(x,y)]=a(x,y) is the pure absorption and -Im[ã(x,y)]=φ(x,y) is the phase shift.

The complex amplitude distribution of the wave front in the hologram plane can be calculated by the Fresnel–Kichhoff diffraction integral,

Ũh(x,y)=1iλzexp(i2πzλ)Ũ(ξ,η)exp{iπλz[(xξ)2+(yη)2]}dξdη,
(3)

where λ is the wavelength of the x ray and z is the distance between the object and the hologram. The exponential term exp[i2πz/λ] is a constant in a given (x,y) plane, and it will only add a constant phase shift to each pixel; thus we will not take this term into account in the following discussion. Now we define a convolution kernel function hz(x,y) as

hz(x,y)=iλzexp[iπλz(x2+y2)];
(4)

then we have

Ũh(x,y)=Ũ(x,y)hz(x,y),
(5)

where * represents the convolution operation. With Eqs. (2) and (5), the intensity distribution in the hologram plane can be written as

Ih(x,y)=Ũh(x,y)Ũh(x,y)
1ã(x,y)hz(x,y)ã(x,y)hz(x,y).
(6)

If we record two holograms in the distances of z and 2z, we will have two intensity distributions I1(x,y) and I2(x,y):

{I1(x,y)=1ã(x,y)hz(x,y)ã(x,y)hz(x,y)I2(x,y)=1ã(x,y)h2z(x,y)ã(x,y)h2z(x,y).
(7)

From these two equations and the relations of ha(x,y) * hb(x,y)=ha+b(x,y), h*z(x,y)=-hz(x,y), h 0(x,y)=δ(x,y) [10

10. L. Onural and P. D. Scott, “Digital decoding of in-line holograms,” Opt. Eng. 26, 1124–1132 (1987).

], we can obtain the following result,

I1hz+I1hzI2=1(ã+ã).
(8)

Therefore we can obtain the real part of the complex absorption ã(x,y) as

Re[ã]=(I2+1)2I1Re[hz].
(9)

Next, we let ã=ar+iai, and from Eq. (6) we have

I1=1arhzarhz+i(aihz-aihz).
(10)

Now we define I′(x,y) as

I'=i(1I1arhzarhz)
=aihzaihz.
(11)

Since I1(x,y) and ar(x,y) are all known to us, then I′(x,y) becomes known as well. Then we can obtain ai(x,y) from I′(x,y). In fact, from Eq. (11) we can obtain the following relations:

I'hz=aiaih2z,
I'h3z=aih2zaih4z,
I'h5z=aih4zaih6z,
(12)

If we add up all the equations, we can derive

ai=I'hz+I'h3z++I'h(2n1)z+aih2nz.
(13)

The final term in this equation is a residual error. It looks like a hologram with a recording distance of 2nz. If this distance is large enough, the effect from errors is small and can be neglected. Therefore, we can obtain the distribution of ai as

Im[ã]=ai=I'k=1Mh(2k1)z.
(14)

3. Simulations

We tested our new algorithm by computer simulation. First, we consider a simple one-dimensional sample to demonstrate our algorithm. The sample is an annular object with an internal diameter of 0.5 µm and an external diameter of 1.5 µm. The parameters in the simulation are as follows: λ=1 nm; z=1 mm; discrete sampling number N, 256; and total view field width, 16 µm. The absorption of the object is 10%, and the phase shift varies from 0 to 0.07, which is shown in Fig. 1 as curves (a) and (d). The reconstructions using our algorithm and an optical algorithm are shown in Fig. 1, from which we can see that the absorption and the phase shift are clearly retrieved by our algorithm.

Next, we study some factors that might affect the reconstruction’s quality. These factors include parameter of M in Eq. (14) and the position error of z in Eq. (3). In addition, because our algorithm is based on the weak absorption and weak phase-shift assumption according to Eq. (2), we study the reconstruction’s error as a function of the magnitude of the absorption and phase shift as well. Thus we define the rms absorption and phase errors as

Fig. 1. One-dimensional simulation results. Absorption (upper) and phase-shift (bottom) profile curve of the object, (a) and (d) real object, (b) and (e) reconstruction by two-intensity algorithm, (c) and (f) reconstruction by optical algorithm.
Ea=x{Re[ã(x)]a(x)}2N,
Ep=x{Im[ã(x)]ϕ(x)}2N.
(15)

The value of M will affect the result of phase shift Ep only. The relations of M-Ep with different N are shown in Fig. 2(a). From these curves in Fig. 2(a) we can find that the minimum Ep appear at different M values for different N; the optimized value of M is approximately 17 or 45 for N equal to 256 or 512. For large N value, optimized M value also becomes larger. Figures 2(b) and 2(c) are the Ea and Ep with different absorption and phase-shift magnitudes, from which we can find that the errors increase greatly when the magnitude of absorption or phase shift increases. Thus this algorithm is valid only for weak absorption and weak phase-shift situations, as we assumed before in Eq. (2).

In practical experiments, the recording distances of the two holograms must have some errors compared with the ideal value. These errors include distance error Δz and lateral displacement error Δx. The effect of these errors can be studied by digital simulations. First, we keep the recording distance z of the first hologram for 1 mm and vary the recording distance of the second hologram to see the effect caused by distance error. The distance of 2z is varied from 1.7 to 2.3 mm. The reconstructions are calculated as the function of 2z, shown in Fig. 3. The rms errors of the reconstruction are shown in Fig. 4(a). If we want both the rms absorption error and phase error to be less than 1×10-4, the distance of the second hologram should be in the range of 1.92–2.08 mm, which means a position accuracy of ±4%. Figure 4(b) shows the rms errors as the function of lateral displacement Δx, from which we can find that the lateral displacement has a magnificent effect on the reconstructions. The spacing interval of the plotted data in Fig. 4(b) is 1 pixel. Therefore, a lateral position accuracy of ±3 pixels or ±0.2 µm is need to ensure that the rms error is less than 1×10-3.

Fig. 2. rms error: (a) Ep as the function of M, (b) and (c) Ea and Ep as the function of phase shift with different absorption.
Fig. 3. Reconstructions with different recording distance of the second hologram (N=256, M=17): (a) absorption, (b) phase shift.
Fig. 4. rms errors with different positioning error (N=256, M=17): (a) with distance error of 2z, (b) with lateral alignment error Δx.

Next, we test our algorithm with a relatively complicated two-dimensional object, shown in Fig. 5(a). It is an elliptical cell-like object with some inner structures. The absorption and phase shift are proportional to the density of the object. The maximum absorption and phase shift are 10% and 0.1 rad, respectively. There are two small squares in the upper-right and lower-left corners. These two squares act as pure phase objects and do not have any absorption in the formation of the two holograms. The phase shifts that they generate are 0.1 and 0.05 rad, respectively. Figures 5(b) and 5(c) are the holograms in the two different distances.

Figure 5(d) is the reconstruction of the real part of the complex amplitude in the object plane, corresponding to the absorption of the object. The object is clearly reconstructed, and the twin-image noise effect can hardly be seen in this image. The contrast of Fig. 5(d) is reversed compared with Fig. 5(a). This is because the gray level in Fig. 5(a) represents the transmissivity, not the absorption. The two squares can hardly be seen in Fig. 5(d), which is in agreement with our expectations, because the squares are pure phase objects and should not appear in the absorption image. Figure 5(e) shows the imaginary part of the object plane with the parameter M=45. Not only the cell-like object but also the two squares are clearly reconstructed in this figure.

4. Conclusion

In this paper we have proposed a new digital-holography reconstruction algorithm based on two holograms with which an object with complex refractive index can be reconstructed without the twin-image noise. This algorithm is noniterative and therefore does not have the problem of convergence that exists in the iterative phase-retrieval algorithm, and the time consumption is quite small as well. Theoretical analysis and simulation results show that it is an effective means to reconstruct the complex object and eliminate the twin-image noise as well. This algorithm is especially applicable for phase-type objects in x-ray holography, because better contrast and lower dosage can be achieved by this method compared with the absorption-contrast imaging technique.

Fig. 5. Digital simulations of a cell-like sample: (a) original object, (b) and (c) holograms in distance of z and 2z, (d) reconstruction of the real part, (e) reconstruction of the imaginary part (M=45, N=512).

Acknowledgment

This research was supported under project 10105002 by the National Natural Science Foundation of China.

References and links

1.

D. Sayer, J. Kirz, R. Feder, D. M. Kim, and E. Spiller, “Potential operating region for ultrasoft X-ray microscopy of biological materials,” Science 196, 1339–1340 (1977). [CrossRef]

2.

M. R. Howells and C. J. Jacobsen, “X-ray holography,” Synchrotron Radiation News , 3(4), 23–28 (1990). [CrossRef]

3.

I. McNulty, “The future of X-ray holography,” Nucl. Instrum. Methods A 347, 170–176 (1994). [CrossRef]

4.

C. Jacobsen, M. R. Howells, J. Kirz, and S. Rothman, “X-ray holographic microscopy using photoresists,” J. Opt. Soc. Am. A 7, 1847–1861 (1990). [CrossRef]

5.

P. Spanne, C. Raven, I. Snigireva, and A. Snigirev, “In-line holography and phase-contrast microtomography with high energy x-rays,” Phys. Med. Biol. 44(3), 741–749 (1999). [CrossRef] [PubMed]

6.

T. E. Gureyev, C. Raven, A. Snigirev, I. Snigireva, and S. W. Wilkins, “Hard X-ray quantitative noninterferometric phase-contrast microscopy,” J. Phys. D 32, 563–567 (1999). [CrossRef]

7.

K. Hirano and A. Momose, “Development of an X-ray interferometer for high-resolution phase-contrast X-ray imaging,” Jpn. J. Appl. Phys. 38, L1556–L1558 (1999). [CrossRef]

8.

J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt. 21, 2758–2769 (1982). [CrossRef] [PubMed]

9.

G. Liu and P. D. Scott, “Phase retrieval and twin-image elimination for in-line Fresnel holograms,” J. Opt. Soc. Am. A 4, 159–165 (1987). [CrossRef]

10.

L. Onural and P. D. Scott, “Digital decoding of in-line holograms,” Opt. Eng. 26, 1124–1132 (1987).

11.

M. H. Maleki and A. J. Devaney, “Noniterative reconstruction of complex-valued objects from two intensity measurements,” Opt. Eng. 33, 3243–3253 (1994). [CrossRef]

12.

T. Xiao, H. Xu, Y. Zhang, J. Chen, and Z. Xu, “Digital image decoding for in-line X-ray holography using two holograms,” J. Mod. Opt. 45, 343–353 (1998). [CrossRef]

13.

L. Xu, J. Miao, and A. Asundi, “Properties of digital holography based on in-line configuration,” Opt. Eng. 39, 3214–3219 (2000). [CrossRef]

OCIS Codes
(090.0090) Holography : Holography
(100.3010) Image processing : Image reconstruction techniques
(100.5070) Image processing : Phase retrieval
(340.7460) X-ray optics : X-ray microscopy

ToC Category:
Research Papers

History
Original Manuscript: January 17, 2003
Revised Manuscript: March 14, 2003
Published: March 24, 2003

Citation
Yuxuan Zhang and Xinyi Zhang, "Reconstruction of a complex object from two in-line holograms," Opt. Express 11, 572-578 (2003)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-11-6-572


Sort:  Journal  |  Reset  

References

  1. D. Sayer, J. Kirz, R. Feder, D. M. Kim, and E. Spiller, �??Potential operating region for ultrasoft X-ray microscopy of biological materials,�?? Science 196, 1339-1340 (1977). [CrossRef]
  2. M. R. Howells and C. J. Jacobsen, �??X-ray holography,�?? Synchrotron Radiation News, 3(4), 23-28 (1990). [CrossRef]
  3. I. McNulty, �??The future of X-ray holography,�?? Nucl. Instrum. Methods A 347, 170-176 (1994). [CrossRef]
  4. C. Jacobsen, M. R. Howells, J. Kirz, and S. Rothman, �??X-ray holographic microscopy using photoresists,�?? J. Opt. Soc. Am. A 7, 1847-1861 (1990). [CrossRef]
  5. P. Spanne, C. Raven, I. Snigireva, and A. Snigirev, �??In-line holography and phase-contrast microtomography with high energy x-rays,�?? Phys. Med. Biol. 44(3), 741-749 (1999). [CrossRef] [PubMed]
  6. T. E. Gureyev, C. Raven, A. Snigirev, I. Snigireva, and S. W. Wilkins, �??Hard X-ray quantitative noninterferometric phase-contrast microscopy,�?? J. Phys. D 32, 563-567 (1999). [CrossRef]
  7. K. Hirano and A. Momose, �??Development of an X-ray interferometer for high-resolution phase-contrast X-ray imaging,�?? Jpn. J. Appl. Phys. 38, L1556-L1558 (1999). [CrossRef]
  8. J. R. Fienup, �??Phase retrieval algorithms: a comparison,�?? Appl. Opt. 21, 2758-2769 (1982). [CrossRef] [PubMed]
  9. G. Liu and P. D. Scott, �??Phase retrieval and twin-image elimination for in-line Fresnel holograms,�?? J. Opt. Soc. Am. A 4, 159-165 (1987). [CrossRef]
  10. L. Onural and P. D. Scott, �??Digital decoding of in-line holograms,�?? Opt. Eng. 26, 1124-1132 (1987).
  11. M. H. Maleki and A. J. Devaney, �??Noniterative reconstruction of complex-valued objects from two intensity measurements,�?? Opt. Eng. 33, 3243-3253 (1994). [CrossRef]
  12. T. Xiao, H. Xu, Y. Zhang, J. Chen, and Z. Xu, �??Digital image decoding for in-line X-ray holography using two holograms,�?? J. Mod. Opt. 45, 343-353 (1998). [CrossRef]
  13. L. Xu, J. Miao, and A. Asundi, �??Properties of digital holography based on in-line configuration,�?? Opt. Eng. 39, 3214�??3219 (2000). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited