Off-axis illumination direct-to-digital holography

What is claimed is:

1. A method of recording an off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprising,

reflecting a reference beam from a reference mirror at a non-normal angle;

reflecting an object beam from an object at an angle with respect to an optical axis defined by a focusing lens;

focusing the reference beam and the object beam at a focal plane of a digital recorder to form the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis;

digitally recording the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis;

Fourier analyzing the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes by transforming axes of the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined as an angle between the reference beam and the object beam;

applying a digital filter to cut off signals around an original origin; and then

performing an inverse Fourier transform.

2. The method of claim 1, further comprising refracting the object beam with an object objective before reflecting the object beam from the object at an angle with respect to the optical axis defined by the focusing lens; and refracting the object beam with the object objective after reflecting the object beam from the object at an angle with respect to the optical axis defined by the focusing lens.

3. The method of claim 1, further comprising fusing the off-axis illuminated spatially heterodyne hologram with at least one hologram selected from the group consisting of an on-axis illuminated spatially heterodyne hologram and another off-axis illuminated spatially heterodyne hologram to compute a single reconstructed image.

4. The method of claim 1, wherein the off-axis illuminated spatially heterodyne hologram is an off-axis illuminated spatially low-frequency heterodyne hologram.

5. The method of claim 1, wherein a region of the true spectrum of the object observed by the off-axis illuminated spatially heterodyne hologram is represented by

G.sub.k (q)=.mu..sub.k exp(j.gamma..sub.k).multidot.F(q)W.sub.k (q),

where q can represent a sample in the discrete Fourier domain or a coordinate vector in the continuous Fourier plane, .mu..sub.k and .gamma..sub.k are the fringe contrast and phase offset, respectively, for the off-axis illuminated spatially heterodyne hologram k, W.sub.k (q) is a window function that represents the region of the spectrum observed by the off-axis illuminated spatially heterodyne hologram k, F(q) is the true spectrum of the object, and j is the square root of negative one.

6. The method of claim 5, wherein W.sub.k (q) is modeled with circularly symmetric Butterworth function.

7. The method of claim 1, further comprising fusing the off-axis illuminated spatially heterodyne hologram with at least one hologram selected from the group consisting of an on-axis illuminated spatially heterodyne hologram and another off-axis illuminated spatially heterodyne hologram to compute a single reconstructed image,

wherein a region of the true spectrum of the object observed by a hologram k selected from the group consisting of the off-axis illuminated spatially heterodyne hologram and the on-axis illuminated spatially heterodyne hologram is represented by

G.sub.k (q)=.mu..sub.k exp(j.gamma..sub.k).multidot.F(q)W.sub.k (q),

where q represents a sample in the discrete Fourier domain, .mu..sub.k and .gamma..sub.k are the fringe contrast and phase offset, respectively, for the hologram k, W.sub.k (q) is a window function that represents the region of the spectrum observed by the hologram k, F(q) is the true spectrum of the object, and j is the square root of negative one, and

an estimate of the computed single reconstructed image is represented by ##EQU9##

where ##EQU10##

where c(q) represents a regularization parameter.

8. The method of claim 1, wherein region of the true spectrum of the object observed by the off-axis illuminated spatially heterodyne hologram is represented by

G.sub.k (q)=F(q)W.sub.k (q)

where q represents a sample in the discrete Fourier domain, W.sub.k (q) is a window function that represents the region of the spectrum observed by the off-axis illuminated spatially heterodyne hologram k, and F(q) is the true spectrum of the object.

9. The method of claim 8, wherein W.sub.k (q) is a function, at least in-part of a Butterworth filter.

10. The method of claim 1, further comprising fusing the off-axis illuminated spatially heterodyne hologram with at least one hologram selected from the group consisting of an on-axis illuminated spatially heterodyne hologram and another off-axis illuminated spatially heterodyne hologram to compute a single reconstructed image,

wherein a region of the true spectrum of the object observed by a hologram k selected from the group consisting of the off-axis illuminated spatially heterodyne hologram and the on-axis illuminated spatially heterodyne hologram is represented by

G.sub.k (q)=F(q)W.sub.k (q)

where q represents a sample in the discrete Fourier domain, W.sub.k (q) is a window function that represents the region of the spectrum observed by the hologram k, and F(q) is the true spectrum of the object, and

an estimate of the computed single reconstructed image is represented by ##EQU11##

where ##EQU12##

where c(q) represents a regularization parameter.

11. The method of claim 1, wherein the step of digitally recording includes detecting the beams with a CCD camera that defines pixels.

12. The method of claim 1, further comprising storing the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis as digital data.

13. The method of claim 1, further comprising replaying the Fourier analyzed off-axis illuminated spatially heterodyne hologram.

14. The method of claim 1, further comprising transmitting the Fourier analyzed off-axis illuminated spatially heterodyne hologram.

15. A computer program, comprising computer or machine readable program elements translatable for implementing the method of claim 1.

16. A machine readable media comprising data generated by the method of claim 1.

17. An apparatus operable to digitally record an off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprising:

a laser;

a beamsplitter optically coupled to the laser;

a reference beam mirror optically coupled to the beamsplitter;

a focusing lens optically coupled to the reference beam mirror;

a digital recorder optically coupled to the focusing lens; and

a computer that performs a Fourier transform, applies a digital filter, and performs an inverse Fourier transform,

wherein a reference beam is incident upon the reference beam mirror at a non-normal angle, an object beam is incident upon an object at an angle with respect to an optical axis defined by the focusing lens, the reference beam and the object beam are focused by the focusing lens at a focal plane of the digital recorder to form the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis which is recorded by the digital recorder, and the computer transforms axes of the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined by an angle between the reference beam and the object beam and cuts off signals around an original origin before performing the inverse Fourier transform.

18. The apparatus of claim 17, further comprising an object objective optically coupled between the beamsplitter and the object.

19. The apparatus of claim 17, wherein the laser is moveable relative to the beamsplitter.

20. The apparatus of claim 17, wherein the beamsplitter, the reference beam mirror and the digital recorder define a Michelson geometry.

21. The apparatus of claim 17, wherein the beamsplitter, the reference beam mirror and the digital recorder define a Mach-Zehner geometry.

22. The apparatus of claim 17, further comprising a digital storage medium coupled to the computer for performing a Fourier transform, applying a digital filter, and performing an inverse Fourier transform.

23. The apparatus of claim 17, wherein the digital recorder includes a CCD camera that defines pixels.

24. The apparatus of claim 23, wherein the angle between the reference beam and the object beam, and a magnification provided by the focusing lens, are selected in order that the digital recorder may resolve features of the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis and two fringes, each having two pixels per fringe, are provided.

25. A machine readable media comprising data generated using the apparatus of claim 17.

Description Submit all comments and votes

BACKGROUND OF THE INVENTION

1. Field of the Invention

The invention relates generally to the field of direct-to-digital holography (interferometry). More particularly, the invention relates to off-axis illumination for improved resolution in direct-to-digital holography.

2. Discussion of the Related Art

Prior art direct-to-digital holography (DDH), sometimes called direct-to-digital interferometry, is known to those skilled in the art. For instance, FIG. 1 illustrates one simplified embodiment of a DDH system. Light from a laser source 105 is expanded by a beam expander/spatial filter 110 and then travels through a lens 115. Subsequently, the expanded filtered light travels to a beamsplitter 120. The beamsplitter 120 may be partially reflective. The portion of light reflected from the beamsplitter 120 constitutes an object beam 125 which travels to the object 130. The portion of the object beam 125 is that is reflected by the object 130 then passes through the beamsplitter 120 and travels to a focusing lens 145. This light then passes through the focusing lens 145 and travels to a charge coupled device (CCD) camera (not shown).

The portion of the light from the lens 115 that passes through the beamsplitter 120 constitutes a reference beam 135. The reference beam 135 is reflected from a mirror 140 at a small angle. The reflected reference beam 135 from the mirror then travels toward the beamsplitter 120. The portion of the reference beam 135 that is reflected from the beamsplitter 120 then travels through the focusing lens 145 and toward the CCD camera (not shown). The object beam 125 from the focusing lens 145 and the reference beam 135 from the focusing lens 145 constitute a plurality of object and reference waves 150 and will interfere at the CCD to produce the interference pattern characteristic of a hologram as noted in U.S. Pat. No. 6,078,392.

In FIG. 1, the object beam 125 is parallel to, and coincident with, the optical axis 127. This type of DDH set-up can be referred to as on-axis illumination.

A limitation of this technology has been that the imaging resolution of the DDH system is limited by the optics of the system. The most notable limitation of the optics is the aperture stop, which is required to prevent degradation of the image quality due to aberrations. With regard to a two-dimensional Fourier plane, only object spatial frequencies within a circle of radius q0 can be transmitted. In the case of on-axis illumination, the aperture with radius q0 appears centered on a zero spatial frequency (q=0). What is needed, therefore, is an approach that permits spatial frequencies outside the circle of radius q0 to be transmitted.

SUMMARY OF THE INVENTION

There is a need for the following aspects of the invention. Of course, the invention is not limited to these aspects.

According to an aspect of the invention, a process of recording an off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprises, reflecting a reference beam from a reference mirror at a non-normal angle; reflecting an object beam from an object at an angle with respect to an optical axis defined by a focusing lens; focusing the reference beam and the object beam at a focal plane of a digital recorder to form the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis; digitally recording the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis; Fourier analyzing the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes by transforming axes of the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined as an angle between the reference beam and the object beam; applying a digital filter to cut off signals around an original origin; and then performing an inverse Fourier transform.

According to another aspect of the invention, a machine operable to digitally record an off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprises: a laser; a beamsplitter optically coupled to the laser; a reference beam mirror optically coupled to the beamsplitter; a focusing lens optically coupled to the reference beam mirror; a digital recorder optically coupled to the focusing lens; and a computer that performs a Fourier transform, applies a digital filter, and performs an inverse Fourier transform, wherein a reference beam is incident upon the reference beam mirror at a non-normal angle, an object beam is incident upon an object at an angle with respect to an optical axis defined by the focusing lens, the reference beam and the object beam are focused by the focusing lens at a focal plane of the digital recorder to form the off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis which is recorded by the digital recorder, and the computer transforms axes of the recorded off-axis illuminated spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined by an angle between the reference beam and the object beam and cuts off signals around an original origin before performing the inverse Fourier transform.

These, and other, aspects of the invention will be better appreciated and understood when considered in conjunction with the following description and the accompanying drawings. It should be understood, however, that the following description, while indicating various embodiments of the invention and numerous specific details thereof, is given by way of illustration and not of limitation. Many substitutions, modifications, additions and/or rearrangements may be made within the scope of the invention without departing from the spirit thereof, and the invention includes all such substitutions, modifications, additions and/or rearrangements.

BRIEF DESCRIPTION OF THE DRAWINGS

The drawings accompanying and forming part of this specification are included to depict certain aspects of the invention. A clearer conception of the invention, and of the components and operation of systems provided with the invention, will become more readily apparent by referring to the exemplary, and therefore nonlimiting, embodiments illustrated in the drawings, wherein identical reference numerals designate the same elements. The invention may be better understood by reference to one or more of these drawings in combination with the description presented herein. It should be noted that the features illustrated in the drawings are not necessarily drawn to scale.

FIG. 1 illustrates a schematic view of a conventional direct-to-digital holography apparatus, appropriately labeled "PRIOR ART."

FIG. 2 illustrates a schematic view of an off-axis illumination direct-to-digital holography apparatus (interferometer) in an on-axis position, representing an embodiment of the invention.

FIG. 3 illustrates a schematic view of the off-axis illumination direct-to-digital holography apparatus (interferometer) of FIG. 2 in an off-axis position.

FIG. 4 illustrates a two-dimensional Fourier plane for an object showing one instance of on-axis illumination and one instance of off-axis illumination, representing an embodiment of the invention.

FIG. 5 illustrates a two-dimensional Fourier plane for an object showing one instance of on-axis illumination and four instances of off-axis illumination, representing an embodiment of the invention.

FIG. 6 illustrates a two-dimensional Fourier plane for an object showing all spatial frequencies that contribute to a fused spectrum (merged image), representing an embodiment of the invention.

FIGS. 7A-7F illustrate the discrete Fourier spectra (7A-7E) obtained from five different holograms and a fused spectrum (7F), representing embodiments of the invention.

FIG. 8 illustrates an object amplitude reconstructed from the on-axis illuminated hologram, representing an embodiment of the invention.

FIG. 9 illustrates an object amplitude reconstructed from a fused result, representing an embodiment of the invention.

FIG. 10 illustrates the intersection of the two apertures W.sub.0 (q) and W.sub.k (q), representing an embodiment of the invention.

DESCRIPTION OF PREFERRED EMBODIMENTS

The invention and the various features and advantageous details thereof are explained more fully with reference to the nonlimiting embodiments that are illustrated in the accompanying drawings and detailed in the following description. Descriptions of well known starting materials, processing techniques, components and equipment are omitted so as not to unnecessarily obscure the invention in detail. It should be understood, however, that the detailed description and the specific examples, while indicating preferred embodiments of the invention, are given by way of illustration only and not by way of limitation. Various substitutions, modifications, additions and/or rearrangements within the spirit and/or scope of the underlying inventive concept will become apparent to those skilled in the art from this disclosure.

Within this application several publications are referenced by Arabic numerals within parentheses. Full citations for these, and other, publications may be found at the end of the specification immediately preceding the claims after the section heading References. The disclosures of all these publications in their entireties are hereby expressly incorporated by reference herein for the purpose of indicating the background of the invention and illustrating the state of the art.

The below-referenced U.S. Patents, and allowed U.S. patent application, disclose embodiments that were satisfactory for the purposes for which they are intended. The entire contents of U.S. Pat. No. 6,078,392, issued Jun. 20, 2000 to C. E. Thomas, L. R. Baylor, G. R. Hanson, D. A. Rasmussen, E. Voelkl, J. Castracane, M. Simkulet and L. Clow, entitled "Direct-to-Digital Holography and Holovision" are hereby expressly incorporated by reference herein for all purposes. The entire contents of allowed U.S. patent application Ser. No. 09/477,267, filed Jan. 4, 2000 by C. E. Thomas and G. R. Hanson, entitled "Improvements To Acquisition and Replay Systems" are hereby expressly incorporated by reference herein for all purposes.

This application contains disclosure that also contained in copending U.S. Ser. Nos. 10/234,043, filed Sep. 3, 2002; and 10/234,042 , filed Sep. 3, 2002, the entire contents of all of which are hereby expressly incorporated by reference for all purposes.

In general, the context of the invention can include obtaining, storing and/or replaying digital data. The context of the invention can include processing digital data that represents an image. The context of the invention can also include transforming data from multiple images into a merged image.

The invention can include a method of acquiring improved resolution holographic imagery from a direct-to-digital holography system using off-axis illumination. The invention can also include an apparatus for acquiring improved resolution holographic imagery with a direct-to-digital holography (DDH) system that uses off-axis illumination.

In general, the object to be observed (imaged) is optically coupled to an illumination source via one or more optical components. As discussed with regard to FIG. 1, the illumination beam is typically passed through the center of the target objective (i.e., lens system) along, and thus parallel to, the optical axis. This type of DDH configuration can be referred to as "on-axis illumination" and allows spatial frequencies (q) of the object to be acquired up to a certain limit (q0), which is determined by the objective aperture.

The invention can include an "off-axis illumination" scenario, where the illumination source is displaced laterally so that the beam will pass through the object objective off-center yet still parallel to the optical axis. The illumination will, due to the focusing effect of the objective, be incident upon the object at some angle to the optical axis. Due to this off-axis illumination, higher spatial frequencies (q>q0) of the object can pass through the objective aperture, and thus be observed, than can with on-axis illumination. This is an important advantage of the invention.

The invention can include an extended DDH system (apparatus) adapted to digitally capture the on-axis- and one, or more, off-axis-illuminated holograms of the same object. The invention can also include analyzing and/or processing (fusing) the digitally captured data. The resulting, fused image will contain a wider range of spatial frequencies than in any of the original holograms, thus providing a significant increase in the nominal imaging resolution of the system compared to the case where no off-axis-illuminated data is available.

As noted above, the imaging resolution of fundamental DDH systems is limited by the optics, most notably the aperture stop, which is required to prevent degradation of the image quality due to aberrations. The aperture stop is required to prevent aliasing of higher frequencies and subsequent degradation of imaging quality. This means the optics of the DDH system are such that only object spatial frequencies within a circle of radius q0 can be transmitted. In the case of on-axis illumination, the aperture with radius q0 appears centered on a zero spatial frequency (q=0). In the case of off-axis illumination, the aperture with radius q0 appears shifted (e.g., to the left) in the frequency domain. This implies that in the direction in which the aperture is shifted, spatial frequencies with q>q0 are transmitted. On the downside, some spatial frequencies with q close to q0 are "lost" in the opposite direction. By acquiring a second image with the illumination shifted in the opposite direction, the aperture appears shifted (e.g., to the right) and thus the spatial frequencies "lost" from the first image are regained with additional frequencies beyond q0. Fusing the information from the two images results in one image with better resolution. Since DDH records the phase information on the complex image wave, the information from both (or more) images can be fused with surprisingly advantageous results. The invention improves the resolution of generic object structures regardless of orientation.

The invention can include an extension of the fundamental DDH system to automatically capture both on-axis and off-axis illuminated holograms. The invention can also include methods to analyze and fuse the results of these holograms to produce a representation of the observed object with more spatial resolution than available in the prior DDH art.

As evident in FIG. 1, the object beam 125 is parallel to the optical axis 127. As noted above, this set-up can be referred to as on-axis illumination. Off-axis illumination, on the other hand, refers to the case where the object beam 125 is incident upon the object 130 at some angle with respect to the optical axis 127 (an example is illustrated by the object team 215, 305 shown in FIG. 3). There are many methods to achieve off-axis illumination; the approach presented hereafter is intended to serve only as a representative, and therefore non-limiting example.

Referring to FIGS. 2 and 3, an embodiment of an off-axis illumination DDH apparatus is illustrated. In FIGS. 2 and 3, there are two primary modifications from FIG. 1. A first modification is that the laser source 105, the beam expander/spatial filter 110, and the lens 115 are grouped into a computer-controlled, moveable enclosure 205. The enclosure 205 can be movable along an axis that is substantially parallel to the optical axis 127. In more detail, the enclosure 205 can be movable along an axis that is substantially coplanar with a normal to the beamsplitter 120.

Still referring to FIGS. 2 and 3, a second modification is the addition of the object objective 210. In FIG. 2, the laser source enclosure 205 is positioned so that the object beam 125 reflects off of the beamsplitter 120 to pass through the center of the object objective 210. The object beam 125 then leaves the object objective 210 and is incident upon the object 130, centered around the optical axis 127. In this configuration, on-axis illumination is achieved and the system of FIG. 2 is effectively the same as that in FIG. 1.

In FIG. 3, however, the laser source enclosure 205 is shifted (up in this particular configuration) so that the object beam 125 passes through the object objective 210 off-center. Of course, the laser source enclosure 205 can alternatively be shifted down. Because of the focusing properties of the object objective 210, the object beam 215 leaving the object objective 210 is incident upon the object 130 at some angle with respect to the optical axis 127, thereby achieving off-axis illumination. Thus, the object beam 215 can be incident upon the object 130 substantially non-parallel to the optic axis 127. The object beam 305 reflected from the object passes back through the object objective 210 off-axis, but due to the optical properties of the object objective 210 and the focusing lens 150 is still focused on the CCD (not shown). In the off-axis illumination case, the properties of diffraction.sup.(1) imply that the hologram formed at the CCD by the interference of the object beam 305 and the reference beam 135 will contain some spatial frequencies of the object that are not observed using on-axis illumination.

Thus, the invention can include an apparatus operable to digitally record a spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprising: a laser; a beamsplitter optically coupled to the laser; a reference beam mirror optically coupled to the beamsplitter; an object optically coupled to the beamsplitter; a focusing lens optically coupled to both the reference beam mirror and the object; a digital recorder optically coupled to the focusing lens; and a computer for performing a Fourier transform, applying a digital filter, and performing an inverse Fourier transform, wherein a reference beam is incident upon the reference beam mirror at a non-normal angle, an object beam is incident upon the object at an angle with respect to an optical axis defined by the focusing lens, the reference beam and an object beam, which constitute a plurality of simultaneous reference and object waves, are focused by the focusing lens at a focal plane of the digital recorder to form a spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis which is recorded by the digital recorder, and the computer transforms axes of the recorded spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined by an angle between the reference beam and the object beam and cuts off signals around an original origin before performing the inverse Fourier transform. The apparatus can include an object objective optically coupled between the beamsplitter and the object. The apparatus can include an aperture stop coupled between the object and the focusing lens. The beamsplitter, the reference beam mirror and the digital recorder can define a Michelson geometry. The beamsplitter, the reference beam mirror and the digital recorder can define a Mach-Zehner geometry. The apparatus can also include a digital storage medium coupled to the computer for performing a Fourier transform, applying a digital filter, and performing an inverse Fourier transform. The digital recorder can include a CCD camera 350 that defines pixels. The apparatus can include a beam expander/spatial filter 230 optically coupled between the laser and the beamsplitter. The angle between the reference beam and the object beam, and a magnification provided by the focusing lens, can be selected in order that the digital recorder may resolve features of the spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis. So that the digital recorder may resolve a feature, two fringes, each having two pixels per fringe, can be provided. The invention can include a spatially heterodyne hologram produced by the above-described apparatus, embodied on a computer-readable medium.

Accordingly, the invention can include a method of recording a spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis, comprising: splitting a laser beam into a reference beam and an object beam; reflecting the reference beam from a reference mirror at a non-normal angle; reflecting the object beam from an object at an angle with respect to an optical axis defined by a focusing lens; focusing the reference beam and the object beam, which constitute a plurality of simultaneous reference and object waves, with the focusing lens at a focal plane of a digital recorder to form a spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis; digitally recording the spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis; Fourier analyzing the recorded spatially heterodyne hologram including spatially heterodyne fringes by transforming axes of the recorded spatially heterodyne hologram including spatially heterodyne fringes in Fourier space to sit on top of a heterodyne carrier frequency defined as an angle between the reference beam and the object beam; applying a digital filter to cut off signals around an original origin; and then performing an inverse Fourier transform. The method can include refracting the object beam with an object objective before reflecting the object beam from an object at an angle with respect to an optical axis defined by a focusing lens and after reflecting the object beam from an object at an angle with respect to an optical axis defined by a focusing lens. The step of transforming axes of the recorded spatially heterodyne hologram can include transforming with an extended Fourier transform. The step of digitally recording can include detecting the beams with a CCD camera that defines pixels. The off-axis illuminated spatially heterodyne hologram can be an off-axis illuminated spatially low-frequency heterodyne hologram; the phrase low-frequency implies that the fundamental fringe spatial frequency is below the Nyquist sampling limit. The method can also include storing the spatially heterodyne hologram including spatially heterodyne fringes for Fourier analysis as digital data. The method can also include replaying the Fourier analyzed spatially heterodyne hologram. The method can also include transmitting the Fourier analyzed spatially heterodyne hologram. The invention can include a spatially heterodyne hologram prepared by the above-described method(s), embodied on a computer-readable medium.

Fusion Methodology

Referring to FIG. 4 a two-dimensional Fourier plane is depicted for a hypothetical object, where zero frequency 405 is the center. The zero frequency for the object 405 corresponds to the heterodyne or carrier frequency in the hologram. All of the spatial frequencies of the object are represented by the filled circle 410. On-axis illumination will capture only those frequencies indicated by the sold-line circle 415. One instance of off-axis illumination, however, will capture the frequencies indicated by the dash-line circle 420. Other instances of off-axis illumination can capture other regions of the object spectrum 410.

A DDH apparatus such as that shown in FIGS. 2 and 3 can be used to capture and store multiple holograms--both on-axis and off-axis illuminated--of the same object area. Realizing that each such hologram contains information about difference spatial frequencies of the object, methods to analyze these holograms and compute a single, reconstructed image will be discussed in detail below. The reconstructed image can contain all of the observed spatial frequencies and thus provide a significantly higher resolution image than any single recorded image. The methods presented hereafter are intended to serve as representative examples and are not intended to preclude extensions, modifications, or other approaches.

The invention can include the capture (digital acquisition) of k=(0, . . . , N) holograms. FIG. 5 illustrates the spatial frequencies contained in each of hologram of a set of holograms where N=4. Similar to FIG. 4, FIG. 5 represents the Fourier plane of the object where the zero frequency 505 is the center point. The solid circle 510 represents spatial frequencies observed with on-axis illumination while the dotted-line circles 515, 520, 525, 530 represent spatial frequencies observed by four different off-axis illuminated holograms. By appropriately fusing the information from all five holograms, the resolution (i.e., bandwidth) of the DDH system is effectively increased and all spatial frequencies indicated by the shaded region 610 in FIG. 6 contribute to the final image.

In the Fourier domain, let the true spectrum of the object be given by F(q). The portion of F(q) (region of F(q)), G.sub.k (q), observed by each hologram is then given by

G.sub.k (q)=.mu..sub.k exp(j.gamma..sub.k).multidot.F(q)W.sub.k (q),

where .mu..sub.k and .gamma..sub.k are the fringe contrast and phase offset, respectively, for hologram k and j is the square root of negative one.

The invention can include normalizing and/or compensating for varying fringe contrasts. One such method for fringe contrast normalization can be described as follows. First we note that G.sub.k (q) above represents the Fourier transform of one acquired object wave (or image) g.sub.k (x):

G.sub.k (q)=.Fourier.{g.sub.k (x)}

where .Fourier.{.multidot.} indicates the Fourier transform operation. We now define two additional images

g.sub.k,0 (x)=.Fourier..sup.-1 {G.sub.k (q)W.sub.0 (q)}=.Fourier..sup.-1 {.mu..sub.k exp(j.gamma..sub.k)F(q)W.sub.k (q)W.sub.0 (q)}

and

g.sub.0,k (x)=.Fourier..sup.-1 {G.sub.0 (q)W.sub.k (q)}=.Fourier..sup.-1 {.mu..sub.0 exp(j.gamma..sub.0)F(q)W.sub.0 (q)W.sub.k (q)}

for k=1, . . . , N where k=0 is assumed to be the reference image and where .Fourier..sup.-1 {.multidot.} indicates the inverse Fourier transform operation. The image g.sub.k,0 (x) is constructed by keeping only the frequencies of g.sub.k (x) common to both g.sub.k (x) and the reference image g.sub.0 (x). Similarly, the image g.sub.0,k (x) comprises only the frequencies of g.sub.0 (x) common to both g.sub.0 (x) and g.sub.k (x). These common frequencies can be visualized as the intersection of the two apertures W.sub.0 (q) and W.sub.k (q) as illustrated in FIG. 10. Once g.sub.k,0 (x) and g.sub.0,k (x) have been constructed, a ratio image .chi..sub.k (x) is computed as ##EQU1##

After computing the ratio image, the relative fringe contrast, .mu..sub.k /.mu..sub.0, can be estimated as the sample mean (average) of the magnitude of the ratio image .chi..sub.k (x): ##EQU2##

where P is the total number of pixels in the (digitized) ratio image and .vertline..chi..sub.k (x.sub.i).vertline. represents the magnitude value of the ratio image at a single pixel location.

The invention can include normalizing and/or compensating for the phase offsets. One method for phase offset normalization is to compute the same ratio image described above for fringe contrast normalization. Then the relative phase offset, exp(j(.gamma..sub.k -.gamma..sub.0)), can computed as the (weighted) sample mean of the phase of the ratio image .chi..sub.k (x): ##EQU3##

where .angle..chi..sub.k (x.sub.i) is the phase (angle) of the (digitized) ratio image at a single pixel location and where .alpha..sub.k (x.sub.i) is a weighting factor. In the most straightforward embodiment, the weighting factor is fixed equal to one. In other embodiments, the weighting factor may be related to the magnitudes of the individual images used to compute the ratio image. This may be used because phase data is often inaccurate when the magnitude is very small. One such embodiment for computing the weighting factor is ##EQU4##

Similarly motivated weighting factors may also be employed. Another method for phase offset estimation is described as follows. Often phase offset is due to errors in locating the carrier frequency. Under such a condition, the phase offset of the ratio image takes the form

.angle..chi..sub.k (x)=exp(j(e.sub.1 x.sub.1 +e.sub.2 x.sub.2 +.gamma..sub.k))

where e.sub.1 and e.sub.2 are related to the error in finding the carrier frequency. In this situation, both e.sub.1 and e.sub.2 must be found in addition to .gamma..sub.k. The above equation can also be written as

.angle..chi..sub.k (x)=cos(e.sub.1 x.sub.1 +e.sub.2 x.sub.2 +.gamma..sub.k)+j sin(e.sub.1 x.sub.1 +e.sub.2 x.sub.2 +.gamma..sub.k).

The real and imaginary parts can be considered separate elements of a two-dimensional vector ##EQU5##

The three parameters in question on the right of the above equation--e.sub.1, e.sub.2, and .gamma..sub.k --can then be estimated from the observations (on the left of the above equation) through any standard non-linear optimization technique.

Thus, after fringe contrast and phase offset normalization, the equation for the observed spectrum in hologram k can be simplified to

G.sub.k (q)=F(q)W.sub.k (q)

In the above equations, W.sub.k (q) is essentially a window function that represents the region of the spectrum observed by hologram k. As a simple example, the window functions from FIG. 5 would be equal to one inside and zero outside the respective circles 510, 515, 520, 525, 530. In a more sophisticated embodiment, the window function can be modeled using a circularly symmetric Butterworth function: ##EQU6##

where q.sub.1 and q.sub.2 represent the horizontal and spatial frequency variables of the vector q, c.sub.1 and c.sub.2 represent the center point of the function in the Fourier plane, r represent the radius of the function, and m is the order of the filter.

The goal of the fusion method is to form an estimate of F(q) from the observations G.sub.k (q). One approach to estimate F(q) is to employ a linear estimator based upon the minimum mean square error criterion; this is the well-known linear, minimum mean-square error (LMMSE) estimator. Using a linear estimator, the estimate for F(q), where q represents a sample in the discrete Fourier domain, is given by ##EQU7##

where the coefficients, as determined by the LMMSE criterion, are given by ##EQU8##

where c(q) is a positive number that represents a regularization parameter that can be selected dependent upon the specifics of the system. This regularization parameter is used to compensate for system noise. In a perfect, zero noise case c(q) could be set to zero, in which case the estimate for F(q) reduces to a simple average. In practice, such a technique produces undesirable artifacts in the fused image. The most straightforward approach is to set c(q) to some constant number for all q based on experimental observations. In more sophisticated approaches, c(q) can be set (either to a constant or to vary with q) based on some computational analysis of the observed holograms and images. Many such approaches for computing the regularization parameter may be devised by those skilled in the art.

Simple averaging or LMMSE estimation are perhaps the most straightforward approaches for computing the fused image. A few alternative optimization criteria for computing the fused image, well-known to those skilled in the art of image processing, might include maximum likelihood (ML) estimation, maximum a posteriori (MAP) estimation, and/or total least squares (TLS) estimation. Note that these examples are not exhaustive.