CA2310851A1 - Method and device for hologram synthesis - Google Patents
Method and device for hologram synthesis Download PDFInfo
- Publication number
- CA2310851A1 CA2310851A1 CA002310851A CA2310851A CA2310851A1 CA 2310851 A1 CA2310851 A1 CA 2310851A1 CA 002310851 A CA002310851 A CA 002310851A CA 2310851 A CA2310851 A CA 2310851A CA 2310851 A1 CA2310851 A1 CA 2310851A1
- Authority
- CA
- Canada
- Prior art keywords
- complex
- transform
- image
- dimensional image
- hologram
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000015572 biosynthetic process Effects 0.000 title description 4
- 238000003786 synthesis reaction Methods 0.000 title description 4
- 230000003287 optical effect Effects 0.000 claims abstract description 29
- 238000004519 manufacturing process Methods 0.000 claims abstract description 14
- 230000001131 transforming effect Effects 0.000 claims abstract description 11
- 239000002131 composite material Substances 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 description 28
- 230000002194 synthesizing effect Effects 0.000 description 5
- 230000001427 coherent effect Effects 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H2210/00—Object characteristics
- G03H2210/20—2D object
Abstract
The invention concerns a method which consists in transforming (E1, E2) a two-dimensional digital image defined by a real function into a complex two-dimensional image (IM2) defined by a complex function, oversampling (E3) the complex image (IM2), simulating (E4) the production of a diffracted image (IM4) resulting from the diffraction of an optical wave by the oversampled image (IM3), and adding (E5) a complex field representing a reference optical wave to the resulting diffracted image (IM4) to produce a hologram. The resulting hologram can be used to produce three-dimensional images or in telecommunications.
Description
METHOD AND DEVICE FOR HOLOGRAM SYNTHESIS
The present invention relates generally to synthesizing holograms and more particularly to synthesizing holograms digitally from two-dirnensional images stored in a memory.
Digital synthesis of holograms from two-dimensional images is used in methods of reproducing three-dimensional images, for example.
Respective holograms are computed for two-dimensional digital images representing a three-dimensional object from different viewpoints. These holograms are then combined to produce a hologram of the object which l0 reproduces a three-dimensional image of the object when it is reproduced physically by a spatial light modulator and illuminated by a coherent wave.
There are other applications of digital synthesis of holograms, in particular in telecommunications, radar, X-rays and sonar.
Digital techniques for synthesizing holograms are known in the art.
For example, the article by S. Michelin et al. entitled "Fourier-transform computer generated hologram: a variation on the off-axis principle"
published in SPIE Conferences 1994, Practical Holography VIII, pages 249-254, describes a method of simulating the production of an analog hologram. The method consists of applying a Fourier transform to a two 2 o dimensional image, adding a complex field representing a reference optical wave to the Fourier transform obtained in this way, and then extracting the amplitude information contained in the sum of the complex field and the Fourier transform. Applying the Fourier transform to the two-dimensional image digitally simulates the production of a "diffracted"
2 5 image which results from the diffraction of a fictitious optical wave by the two-dimensional image. The two-dimensional image is also oversampled before the Fourier transform is applied to it. However, the oversampled two-dimensional image obtained in this way is defined by a real intensity distribution which is not always well suited to computing a complex 30 transform such as a Fourier transform.
The present invention aims to provide a method of synthesizing holograms that is more efficient than those of the prior art.
To this end, a method of producing a hologram from a two dimensional image defined by a real function is characterized in that it 35 comprises the following steps:
transforming the two-dimensional image defined by said real function into a complex two-dimensional image defined by a complex function, oversampling the complex image, simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and adding a complex field representing a reference optical wave to the resulting diffracted image in order to produce said hologram.
The method can further comprise the step of encoding values taken by the amplitude of the sum of said complex field and the resulting diffracted image, so that the hologram can be reproduced on a liquid crystal screen or transmitted over a transmission line, for example.
In the present context, a "real or complex function" means a function of two variables, in the form of digital data, and taking real or complex values, respectively. The real function is typically an intensity distribution while the complex function is a distribution of complex numbers each defined by a real amplitude and a real phase.
The step of transforming the given two-dimensional image into a complex image derives from the original two-dimensional image an image which is defined by complex numbers which optimally represent the real optical field and facilitate the computations employed in the simulation step.
The oversampling step increases the number of pixels of the hologram because the computations employed in subsequent steps apply 2 5 to a greater number of image points. This step can consist of inserting the complex image into a larger image in which the intensity of pixels outside the original complex image is made equal to 0. In this case, implementing the step of oversampling the complex image after the steps of transforming the two-dimensional image into a complex image avoids having to calculate the complex function for points of the oversampled image outside the original complex image.
The transform step typically includes the following steps:
determining amplitude values each depending on the square root of a corresponding value taken by said real function, and associating a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
The present invention relates generally to synthesizing holograms and more particularly to synthesizing holograms digitally from two-dirnensional images stored in a memory.
Digital synthesis of holograms from two-dimensional images is used in methods of reproducing three-dimensional images, for example.
Respective holograms are computed for two-dimensional digital images representing a three-dimensional object from different viewpoints. These holograms are then combined to produce a hologram of the object which l0 reproduces a three-dimensional image of the object when it is reproduced physically by a spatial light modulator and illuminated by a coherent wave.
There are other applications of digital synthesis of holograms, in particular in telecommunications, radar, X-rays and sonar.
Digital techniques for synthesizing holograms are known in the art.
For example, the article by S. Michelin et al. entitled "Fourier-transform computer generated hologram: a variation on the off-axis principle"
published in SPIE Conferences 1994, Practical Holography VIII, pages 249-254, describes a method of simulating the production of an analog hologram. The method consists of applying a Fourier transform to a two 2 o dimensional image, adding a complex field representing a reference optical wave to the Fourier transform obtained in this way, and then extracting the amplitude information contained in the sum of the complex field and the Fourier transform. Applying the Fourier transform to the two-dimensional image digitally simulates the production of a "diffracted"
2 5 image which results from the diffraction of a fictitious optical wave by the two-dimensional image. The two-dimensional image is also oversampled before the Fourier transform is applied to it. However, the oversampled two-dimensional image obtained in this way is defined by a real intensity distribution which is not always well suited to computing a complex 30 transform such as a Fourier transform.
The present invention aims to provide a method of synthesizing holograms that is more efficient than those of the prior art.
To this end, a method of producing a hologram from a two dimensional image defined by a real function is characterized in that it 35 comprises the following steps:
transforming the two-dimensional image defined by said real function into a complex two-dimensional image defined by a complex function, oversampling the complex image, simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and adding a complex field representing a reference optical wave to the resulting diffracted image in order to produce said hologram.
The method can further comprise the step of encoding values taken by the amplitude of the sum of said complex field and the resulting diffracted image, so that the hologram can be reproduced on a liquid crystal screen or transmitted over a transmission line, for example.
In the present context, a "real or complex function" means a function of two variables, in the form of digital data, and taking real or complex values, respectively. The real function is typically an intensity distribution while the complex function is a distribution of complex numbers each defined by a real amplitude and a real phase.
The step of transforming the given two-dimensional image into a complex image derives from the original two-dimensional image an image which is defined by complex numbers which optimally represent the real optical field and facilitate the computations employed in the simulation step.
The oversampling step increases the number of pixels of the hologram because the computations employed in subsequent steps apply 2 5 to a greater number of image points. This step can consist of inserting the complex image into a larger image in which the intensity of pixels outside the original complex image is made equal to 0. In this case, implementing the step of oversampling the complex image after the steps of transforming the two-dimensional image into a complex image avoids having to calculate the complex function for points of the oversampled image outside the original complex image.
The transform step typically includes the following steps:
determining amplitude values each depending on the square root of a corresponding value taken by said real function, and associating a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
By averaging the amplitude values of the hologram, associating a phase with each amplitude value avoids peaks of excessively high amplitude in the resulting hologram of the given two-dimensional image.
The simulation step can include computing one of the following complex transforms: Fourier transform, Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform and a transform which is a composite of at least two of the above transforms.
The simulation step advantageously consists of computing a convolutional product, associated with the oversampled complex image, of two components, by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
Until now, the skilled person has regarded the Fourier transform, which is widely used in optics, as the best possible transform for calculating a convolutional product of this kind. However, experiments conducted by the present inventors have shown that using one of the complex transforms mentioned above other than the Fourier transform produces, for a two-dimensional image, a resultant hologram of much better quality, i.e. which, when it is reproduced physically and illuminated by a coherent source, produces an image associated with the two-dimensional image that is finer than those generally produced by prior art systems.
According to another aspect of the invention, a method of producing a hologram from a two-dimensional image defined by a real function is characterized in that it comprises the following steps:
oversampling the two-dimensional image, transforming the oversampled two-dimensional image into a 3o complex two-dimensional image defined by a complex function, simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and adding a complex field representing a reference optical wave to the resulting diffracted image to produce said hologram.
The invention also provides a system for producing a hologram from a two-dimensional image defined by a real function, characterized in that it comprises:
The simulation step can include computing one of the following complex transforms: Fourier transform, Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform and a transform which is a composite of at least two of the above transforms.
The simulation step advantageously consists of computing a convolutional product, associated with the oversampled complex image, of two components, by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
Until now, the skilled person has regarded the Fourier transform, which is widely used in optics, as the best possible transform for calculating a convolutional product of this kind. However, experiments conducted by the present inventors have shown that using one of the complex transforms mentioned above other than the Fourier transform produces, for a two-dimensional image, a resultant hologram of much better quality, i.e. which, when it is reproduced physically and illuminated by a coherent source, produces an image associated with the two-dimensional image that is finer than those generally produced by prior art systems.
According to another aspect of the invention, a method of producing a hologram from a two-dimensional image defined by a real function is characterized in that it comprises the following steps:
oversampling the two-dimensional image, transforming the oversampled two-dimensional image into a 3o complex two-dimensional image defined by a complex function, simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and adding a complex field representing a reference optical wave to the resulting diffracted image to produce said hologram.
The invention also provides a system for producing a hologram from a two-dimensional image defined by a real function, characterized in that it comprises:
transform means for transforming the two-dimensional image defined by said real function into a complex two-dimensional image defined by a complex function, means for oversampling the complex image, simulator means for simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and means for adding a complex field representing a reference optical wave to the resulting diffracted image to produce said hologram.
1 o The system can further comprise means for encoding values taken by the amplitude of the sum of said complex field and the diffracted image.
The transform means can comprise means for determining amplitude values each depending on the square root of the corresponding value taken by said real function and means for associating a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
The simulator means can comprise means for computing one of the following complex transforms: Fourier transform, Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform and a transform which is a composite of at least two of the above transforms.
The simulator means advantageously comprise means for 2 5 computing a convolutional product, associated with the oversampled complex image, of two components, by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
According to another aspect of the invention, a system for 3o producing a hologram from a two-dimensional image defined by a real function is characterized in that it comprises:
means for oversampling the two-dimensional image, transform means for transforming the oversampled two dimensional image into a complex two-dimensional image defined by a 35 complex function, simulator means for simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and means for adding a complex field representing a reference optical wave to the resulting diffracted image to produce said hologram.
Other advantages of the present invention will become apparent 5 on reading the following detailed description with reference to the accompanying drawings, in which:
- Figure 1 is a flowchart of an algorithm according to the invention, - Figure 2 is a block diagram of a computer executing the algorithm shown in Figure 1, - Figure 3 illustrates the production of a hologram from a two-dimensional image, - Figure 4 shows the oversampling of a two-dimensional image by the algorithm shown in Figure 1, and - Figure 5 is a diagram showing geometrical planes used in the algorithm shown in Figure 1.
Figure 1 shows an algorithm of the invention for synthesizing digital holograms, which is executed by a microprocessor MP associated with a memory MM, both shown in Figure 2.
In a preliminary step E0, a two-dimensional digital image IM, shown diagrammatically in Figure 3, is stored in the memory MM
associated with the microprocessor MP in the form of digital data. The two-dimensional image IM is typically defined by a real function of two variables, particularly by a distribution of intensities f(x,y), where (x,y) represent co-ordinates in a two-dimensional system of axes (O,x,y) associated with the image IM.
In a step E1, the two-dimensional image IM is transformed into a transformed two-dimensional image IM1 which is defined by an amplitude distribution by computing for each point of the image IM a value proportional to the square root of the corresponding intensity value.
In the next step E2, a "pseudorandom" diffuser is generated digitally. This diffuser consists of an "image" having the same number of pixels as the two-dimensional image IM and in which each pixel has an intensity value equal. to 1 and a random phase. Each phase of the diffuser is then associated with a corresponding pixel of the transformed two-dimensional image IM1, to transform the image IM1 into a "complex"
image IM2 in which a complex number defined by an amplitude value and a phase value is determined for each pixel. The pseudorandom diffuser prevents the resulting hologram HO, shown diagrammatically in Figure 3, associated with the image IM, having excessive amplitude level disparities by averaging the amplitude values of the hologram.
In a step E3 the complex image IM2 obtained in step E2 is oversampled, i.e. the image is included in a larger image, as shown in Figure 4. An oversampled image IM3 is formed in this way consisting of the complex image IM2 in a central part PC and of pixels whose amplitude is chosen arbitrarily, for example equal to 0, in a complementary periphery part PF. This oversampling of the complex image IM2 increases the number of pixels of the resultant hologram HO and therefore improves its resolution.
In a step E4, the production of a diffracted image IM4 resulting from the diffraction of a fictitious coherent optical wave DIF by the oversampled complex image IM3 is simulated digitally. To this end, first and second parallel and separate geometrical planes P1 and P2 are defined within a three-dimensional system of axes (O',X,Y,Z) as shown in Figure 5. The first plane P1 includes the oversampled complex image IM3 and the second plane P2 constitutes the plane for computing the hologram HO. The production of the diffracted image IM4 can be 2 o simulated in a manner that is known in the art by applying a Fourier transform to the image IM3. In the method according to the invention, the diffracted image IM4 is preferably determined otherwise, namely by computing, in the plane P2, a convolutional product associated with the oversampled complex image IM3. This convolutional product conforms to scalar diffraction theory. For example, using a Rayleigh-Sommerfeld scalar diffraction formulation, the two components of the convolutional product can respectively correspond to a complex field representing the oversampled complex image IM3 and a complex field representing a spherical optical wave with the same wavelength as the optical wave DIF.
The skilled person however knows other types of convolutional product for computing a diffracted image. The convolutional product computed in step E4 uses parameters including the distance D between the geometrical planes P1 and P2 and the wavelength of the coherent optical wave DIF.
In accordance with the invention, the convolutional product is computed by applying a complex transform, also referred to as a fast complex transform, to the two components of the convolutional product, computing the product of the resulting fast complex transforms, and then applying the fast complex transform which is the inverse of said fast complex transform to the aforementioned product of the fast complex transforms.
To be more precise, if CONY denotes the convolutional product, C1 and C2 its two components and T the fast complex transform, then the convolutional product is written:
CONV = C1~C2=T'T(C1~C2) CONV = T~(T(C1 )T(C2)).
1 o In the present context, the expression "fast complex transform"
means a mathematical transform compatible with scalar optical diffraction theory, i.e. whose resulting transformed functions satisfy the conventional scalar diffraction equations. The fast complex transform must also have the property whereby the fast complex transform of a convolutional product of two components is equal to the product of the respective fast complex transforms of each of said two components. The Fourier transform, the orthogonal polynomial transform, the Paley transform, the Hadamar transform, the Walsh transform, the Hankel transform, the Karhunen-Loeve transform, the multiresolution discrete wavelet transform 2 o and the adaptive wavelet transform are all fast complex transforms which meet the above conditions. Other appropriate fast complex transforms are composites of at least two of the aforementioned transforms, such as a composite of the Walsh transform and the Hadamar transform. The application of a composite of two transforms T1 and T2 to any image I is 2 5 defined in standard mathematical terms by the equation:
(T1 ~T2)(I)=T1 (T2(I)).
Each of the aforementioned fast complex transforms can be used in a specific case. In particular, the fast complex transform can be chosen according to the distance D between the planes P1 and P2. A Fourier 3o transform is appropriate for a large distance D. A Walsh transform is more suitable for a smaller distance D. Also, it has been found that using one of the above-mentioned fast complex transforms other than the Fourier transform gives better results in terms of the quality of the hologram HO than those obtained using the Fourier transform.
35 It should be noted that, because the two-dimensional image IM is transformed into a complex image IM2, computing the convolutional product associated with the image IM in step E4 is more practical than in the prior art since the fast complex transform is applied directly to an image IM3 defined by a complex function and not to an image defined by a real function.
At the exit from step E4, the diffracted image IM4 is defined by a complex field made up of a set of complex numbers each of which is associated with a point of the image IM4. Each of these complex numbers also depends on the image IM3 taken as a whole.
In a next step E5 a complex field simulating a reference optical wave REF with the same wavelength as the optical wave DIF and directed l0 towards the hologram computation plane P2 is added, in the plane P2, to the complex field representing the diffracted image IM4. The amplitude information contained in the resulting complex field is then extracted in order to produce an interference field. The addition of the aforementioned two complex fields is performed by adding, at each point of the diffracted image IM4, the complex number associated with that point and the value at the same point of the complex field representing the reference wave REF. The interference field constitutes the hologram HO of the two-dimensional image IM.
A variant of the Figure 1 algorithm dispenses with the steps E1 and E2 of producing the complex image IM2 and/or the oversampling step E3. In another variant, the oversampling step E3 precedes the step E1.
The hologram HO of the two-dimensional image IM obtained in step E5 is a diffractive field, or grating, which is computed for a particular wavelength, namely the wavelength of the optical waves DIF and REF.
This hologram, which is present in virtual form in step E5, i.e. represented by digital data, is such that, if it is reproduced physically by a holographic screen, illuminating said holographic screen with a laser source emitting at the aforementioned wavelength reproduces the original two-dimensional image IM at a given order of diffraction.
The hologram HO obtained in step E5 is defined digitally by a two-dimensional amplitude function A(u,v), where (u,v) designate co-ordinates in the hologram computation plane P2 which correspond to image spatial frequencies when the fast complex transform chosen in step E4 is a Fourier transform, for example. The two-dimensional amplitude function A(u,v) is deduced from the two-dimensional intensity function f(x,y) describing the two-dimensional image IM, as explained above. In practice, the function A(u,v) is computed only for a series of discrete points (u,v) _ (uk, vq), where k and g are integers. The values that the function A(u,v) takes can nevertheless be spread continuously between a minimum amplitude value and a maximum amplitude value.
In a step E6 of the Figure 1 algorithm, the values taken by the function A(u,v) are quantized and encoded, i.e. each value of this function is associated with a discrete value which is encoded digitally, for example on eight bits. To each pair of discrete points (uk, vq) there then corresponds a discrete amplitude value representing one of 256 gray levels. The amplitudes A(u,v) can also be quantized more simply by 1 o allocating to each amplitude value of A(u,v) the discrete value "0" if said amplitude value is below a predetermined threshold or the discrete value "1" if said amplitude value is above the predetermined threshold.
The encoding step E6 enables the hologram HO to be adapted to suit digital display systems, such as a digitally controlled liquid crystal screen, or to facilitate transmission of the hologram via telecommunication systems.
1 o The system can further comprise means for encoding values taken by the amplitude of the sum of said complex field and the diffracted image.
The transform means can comprise means for determining amplitude values each depending on the square root of the corresponding value taken by said real function and means for associating a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
The simulator means can comprise means for computing one of the following complex transforms: Fourier transform, Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform and a transform which is a composite of at least two of the above transforms.
The simulator means advantageously comprise means for 2 5 computing a convolutional product, associated with the oversampled complex image, of two components, by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
According to another aspect of the invention, a system for 3o producing a hologram from a two-dimensional image defined by a real function is characterized in that it comprises:
means for oversampling the two-dimensional image, transform means for transforming the oversampled two dimensional image into a complex two-dimensional image defined by a 35 complex function, simulator means for simulating the production of a diffracted image resulting from the diffraction of an optical wave by the oversampled complex image, and means for adding a complex field representing a reference optical wave to the resulting diffracted image to produce said hologram.
Other advantages of the present invention will become apparent 5 on reading the following detailed description with reference to the accompanying drawings, in which:
- Figure 1 is a flowchart of an algorithm according to the invention, - Figure 2 is a block diagram of a computer executing the algorithm shown in Figure 1, - Figure 3 illustrates the production of a hologram from a two-dimensional image, - Figure 4 shows the oversampling of a two-dimensional image by the algorithm shown in Figure 1, and - Figure 5 is a diagram showing geometrical planes used in the algorithm shown in Figure 1.
Figure 1 shows an algorithm of the invention for synthesizing digital holograms, which is executed by a microprocessor MP associated with a memory MM, both shown in Figure 2.
In a preliminary step E0, a two-dimensional digital image IM, shown diagrammatically in Figure 3, is stored in the memory MM
associated with the microprocessor MP in the form of digital data. The two-dimensional image IM is typically defined by a real function of two variables, particularly by a distribution of intensities f(x,y), where (x,y) represent co-ordinates in a two-dimensional system of axes (O,x,y) associated with the image IM.
In a step E1, the two-dimensional image IM is transformed into a transformed two-dimensional image IM1 which is defined by an amplitude distribution by computing for each point of the image IM a value proportional to the square root of the corresponding intensity value.
In the next step E2, a "pseudorandom" diffuser is generated digitally. This diffuser consists of an "image" having the same number of pixels as the two-dimensional image IM and in which each pixel has an intensity value equal. to 1 and a random phase. Each phase of the diffuser is then associated with a corresponding pixel of the transformed two-dimensional image IM1, to transform the image IM1 into a "complex"
image IM2 in which a complex number defined by an amplitude value and a phase value is determined for each pixel. The pseudorandom diffuser prevents the resulting hologram HO, shown diagrammatically in Figure 3, associated with the image IM, having excessive amplitude level disparities by averaging the amplitude values of the hologram.
In a step E3 the complex image IM2 obtained in step E2 is oversampled, i.e. the image is included in a larger image, as shown in Figure 4. An oversampled image IM3 is formed in this way consisting of the complex image IM2 in a central part PC and of pixels whose amplitude is chosen arbitrarily, for example equal to 0, in a complementary periphery part PF. This oversampling of the complex image IM2 increases the number of pixels of the resultant hologram HO and therefore improves its resolution.
In a step E4, the production of a diffracted image IM4 resulting from the diffraction of a fictitious coherent optical wave DIF by the oversampled complex image IM3 is simulated digitally. To this end, first and second parallel and separate geometrical planes P1 and P2 are defined within a three-dimensional system of axes (O',X,Y,Z) as shown in Figure 5. The first plane P1 includes the oversampled complex image IM3 and the second plane P2 constitutes the plane for computing the hologram HO. The production of the diffracted image IM4 can be 2 o simulated in a manner that is known in the art by applying a Fourier transform to the image IM3. In the method according to the invention, the diffracted image IM4 is preferably determined otherwise, namely by computing, in the plane P2, a convolutional product associated with the oversampled complex image IM3. This convolutional product conforms to scalar diffraction theory. For example, using a Rayleigh-Sommerfeld scalar diffraction formulation, the two components of the convolutional product can respectively correspond to a complex field representing the oversampled complex image IM3 and a complex field representing a spherical optical wave with the same wavelength as the optical wave DIF.
The skilled person however knows other types of convolutional product for computing a diffracted image. The convolutional product computed in step E4 uses parameters including the distance D between the geometrical planes P1 and P2 and the wavelength of the coherent optical wave DIF.
In accordance with the invention, the convolutional product is computed by applying a complex transform, also referred to as a fast complex transform, to the two components of the convolutional product, computing the product of the resulting fast complex transforms, and then applying the fast complex transform which is the inverse of said fast complex transform to the aforementioned product of the fast complex transforms.
To be more precise, if CONY denotes the convolutional product, C1 and C2 its two components and T the fast complex transform, then the convolutional product is written:
CONV = C1~C2=T'T(C1~C2) CONV = T~(T(C1 )T(C2)).
1 o In the present context, the expression "fast complex transform"
means a mathematical transform compatible with scalar optical diffraction theory, i.e. whose resulting transformed functions satisfy the conventional scalar diffraction equations. The fast complex transform must also have the property whereby the fast complex transform of a convolutional product of two components is equal to the product of the respective fast complex transforms of each of said two components. The Fourier transform, the orthogonal polynomial transform, the Paley transform, the Hadamar transform, the Walsh transform, the Hankel transform, the Karhunen-Loeve transform, the multiresolution discrete wavelet transform 2 o and the adaptive wavelet transform are all fast complex transforms which meet the above conditions. Other appropriate fast complex transforms are composites of at least two of the aforementioned transforms, such as a composite of the Walsh transform and the Hadamar transform. The application of a composite of two transforms T1 and T2 to any image I is 2 5 defined in standard mathematical terms by the equation:
(T1 ~T2)(I)=T1 (T2(I)).
Each of the aforementioned fast complex transforms can be used in a specific case. In particular, the fast complex transform can be chosen according to the distance D between the planes P1 and P2. A Fourier 3o transform is appropriate for a large distance D. A Walsh transform is more suitable for a smaller distance D. Also, it has been found that using one of the above-mentioned fast complex transforms other than the Fourier transform gives better results in terms of the quality of the hologram HO than those obtained using the Fourier transform.
35 It should be noted that, because the two-dimensional image IM is transformed into a complex image IM2, computing the convolutional product associated with the image IM in step E4 is more practical than in the prior art since the fast complex transform is applied directly to an image IM3 defined by a complex function and not to an image defined by a real function.
At the exit from step E4, the diffracted image IM4 is defined by a complex field made up of a set of complex numbers each of which is associated with a point of the image IM4. Each of these complex numbers also depends on the image IM3 taken as a whole.
In a next step E5 a complex field simulating a reference optical wave REF with the same wavelength as the optical wave DIF and directed l0 towards the hologram computation plane P2 is added, in the plane P2, to the complex field representing the diffracted image IM4. The amplitude information contained in the resulting complex field is then extracted in order to produce an interference field. The addition of the aforementioned two complex fields is performed by adding, at each point of the diffracted image IM4, the complex number associated with that point and the value at the same point of the complex field representing the reference wave REF. The interference field constitutes the hologram HO of the two-dimensional image IM.
A variant of the Figure 1 algorithm dispenses with the steps E1 and E2 of producing the complex image IM2 and/or the oversampling step E3. In another variant, the oversampling step E3 precedes the step E1.
The hologram HO of the two-dimensional image IM obtained in step E5 is a diffractive field, or grating, which is computed for a particular wavelength, namely the wavelength of the optical waves DIF and REF.
This hologram, which is present in virtual form in step E5, i.e. represented by digital data, is such that, if it is reproduced physically by a holographic screen, illuminating said holographic screen with a laser source emitting at the aforementioned wavelength reproduces the original two-dimensional image IM at a given order of diffraction.
The hologram HO obtained in step E5 is defined digitally by a two-dimensional amplitude function A(u,v), where (u,v) designate co-ordinates in the hologram computation plane P2 which correspond to image spatial frequencies when the fast complex transform chosen in step E4 is a Fourier transform, for example. The two-dimensional amplitude function A(u,v) is deduced from the two-dimensional intensity function f(x,y) describing the two-dimensional image IM, as explained above. In practice, the function A(u,v) is computed only for a series of discrete points (u,v) _ (uk, vq), where k and g are integers. The values that the function A(u,v) takes can nevertheless be spread continuously between a minimum amplitude value and a maximum amplitude value.
In a step E6 of the Figure 1 algorithm, the values taken by the function A(u,v) are quantized and encoded, i.e. each value of this function is associated with a discrete value which is encoded digitally, for example on eight bits. To each pair of discrete points (uk, vq) there then corresponds a discrete amplitude value representing one of 256 gray levels. The amplitudes A(u,v) can also be quantized more simply by 1 o allocating to each amplitude value of A(u,v) the discrete value "0" if said amplitude value is below a predetermined threshold or the discrete value "1" if said amplitude value is above the predetermined threshold.
The encoding step E6 enables the hologram HO to be adapted to suit digital display systems, such as a digitally controlled liquid crystal screen, or to facilitate transmission of the hologram via telecommunication systems.
Claims (14)
1. A method of producing a hologram (HO) from a two-dimensional image (IM) defined by a real function (f(x,y)), characterized in that it comprises the following steps:
transforming (E1, E2) the two-dimensional image defined by said real function into a complex two-dimensional image (IM2) defined by a complex function, oversampling (E3) the complex image (IM2), simulating (E4) the production of a diffracted image (IM4) resulting from the diffraction of an optical wave (DIF) by the oversampled complex image (IM3), and adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) in order to produce said hologram (HO).
transforming (E1, E2) the two-dimensional image defined by said real function into a complex two-dimensional image (IM2) defined by a complex function, oversampling (E3) the complex image (IM2), simulating (E4) the production of a diffracted image (IM4) resulting from the diffraction of an optical wave (DIF) by the oversampled complex image (IM3), and adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) in order to produce said hologram (HO).
2. A method of producing a hologram (HO) from a two-dimensional image (IM) defined by a real function (f(x,y)), characterized in that it comprises the following steps:
oversampling (E3) the two-dimensional image (IM), transforming (E1, E2) the oversampled two-dimensional image into a complex two-dimensional image defined by a complex function, simulating (E4) the production of a diffracted image resulting from the diffraction of an optical wave (DIF) by the oversampled complex image, and adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image to produce said hologram (HO).
oversampling (E3) the two-dimensional image (IM), transforming (E1, E2) the oversampled two-dimensional image into a complex two-dimensional image defined by a complex function, simulating (E4) the production of a diffracted image resulting from the diffraction of an optical wave (DIF) by the oversampled complex image, and adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image to produce said hologram (HO).
3. A method according to claim 1 or claim 2, further comprising the step (E6) of encoding values taken by the amplitude of the sum of said complex field and the resulting diffracted image (IM4).
4. A method according to any one of claims 1 to 3, wherein said transform step comprises the following steps:
determining (E1) amplitude values each depending on the square root of a corresponding value taken by said real function, and associating (E2) a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
determining (E1) amplitude values each depending on the square root of a corresponding value taken by said real function, and associating (E2) a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
5. A method according to any one of claims 1 to 4, wherein said simulation step (E4) includes calculating a complex transform, said complex transform being a Fourier transform.
6. A method according to any one of claims 1 to 4, wherein said simulation step (E4) includes computing one of the following complex transforms: Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform, a transform which is a composite of at least two of the above transforms and a transform which is a composite of at least one of the above transforms and a Fourier transform.
7. A method according to claim 5 or claim 6, wherein said simulation step (E4) consists of computing a convolutional product, associated with the oversampled complex image, of two components, by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
8. A system for producing a hologram (HO) from a two-dimensional image (IM) defined by a real function (f(x,y)), characterized in that it comprises:
transform means (MP, MM) for transforming (E1, E2) the two-dimensional image defined by said real function into a complex two-dimensional image (IM2) defined by a complex function, means (MP) for oversampling (E3) the complex image (IM2), simulator means (MP) for simulating the production of a diffracted image resulting from the diffraction of an optical wave (DIF) by the oversampled complex image (IM3), and means (MP) for adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) to produce said hologram (HO).
transform means (MP, MM) for transforming (E1, E2) the two-dimensional image defined by said real function into a complex two-dimensional image (IM2) defined by a complex function, means (MP) for oversampling (E3) the complex image (IM2), simulator means (MP) for simulating the production of a diffracted image resulting from the diffraction of an optical wave (DIF) by the oversampled complex image (IM3), and means (MP) for adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) to produce said hologram (HO).
9. A system for producing a hologram (HO) from a two-dimensional image (IM) defined by a real function (f(x,y)), characterized in that it comprises:
means (MP, MM) for oversampling (E3) the two-dimensional image (IM), transform means (MP) for transforming (E1, E2) the oversampled two-dimensional image into a complex two-dimensional image defined by a complex function, simulator means (MP) for simulating the production of a diffracted image (IM4) resulting from the diffraction of an optical wave (DIF) by the oversampled complex image, and means (MP) for adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) to produce said hologram (HO).
means (MP, MM) for oversampling (E3) the two-dimensional image (IM), transform means (MP) for transforming (E1, E2) the oversampled two-dimensional image into a complex two-dimensional image defined by a complex function, simulator means (MP) for simulating the production of a diffracted image (IM4) resulting from the diffraction of an optical wave (DIF) by the oversampled complex image, and means (MP) for adding (E5) a complex field representing a reference optical wave (REF) to the resulting diffracted image (IM4) to produce said hologram (HO).
10. A system according to claim 8 or claim 9, further comprising means (MP) for encoding (E6) values taken by the amplitude of the sum of said complex field and the diffracted image (IM4).
11. A system according to any one of claims 8 to 10, wherein said transform means comprise:
means (MP) for determining (E1) amplitude values each depending on the square root of a corresponding value taken by said real function, and means (MP) for associating (E2) a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
means (MP) for determining (E1) amplitude values each depending on the square root of a corresponding value taken by said real function, and means (MP) for associating (E2) a phase with each of said amplitude values so that an amplitude value and a phase value are defined for each point of the complex image.
12. A system according to any one of claims 8 to 11, wherein said simulator means comprise means (MP) for computing a complex transform, said complex transform being a Fourier transform.
13. A system according to any one of claims 8 to 11, wherein said simulator means comprise means (MP) for computing one of the following complex transforms: Walsh transform, Hankel transform, orthogonal polynomial transform, Hadamar transform, Karhunen-Loeve transform, multiresolution discrete wavelet transform, adaptive wavelet transform, a transform which is a composite of at least two of the above transforms and a transform which is a composite of at least one of the above transforms and a Fourier transform.
14. A system according to claim 12 or claim 13, wherein said simulator means comprise means (MP) for computing a convolutional product, associated with the oversampled complex image, of two components (IM3), by applying the transform which is the inverse of said complex transform to the product of the respective complex transforms of said two components.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP97402798.9 | 1997-11-20 | ||
| EP97402798A EP0921450B1 (en) | 1997-11-20 | 1997-11-20 | Method and device for the synthesis of holograms |
| PCT/EP1998/007539 WO1999027422A1 (en) | 1997-11-20 | 1998-11-17 | Method and device for hologram synthesis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2310851A1 true CA2310851A1 (en) | 1999-06-03 |
Family
ID=8229908
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002310851A Abandoned CA2310851A1 (en) | 1997-11-20 | 1998-11-17 | Method and device for hologram synthesis |
Country Status (12)
| Country | Link |
|---|---|
| US (1) | US6344909B1 (en) |
| EP (1) | EP0921450B1 (en) |
| JP (1) | JP2001524697A (en) |
| AT (1) | ATE201775T1 (en) |
| CA (1) | CA2310851A1 (en) |
| DE (1) | DE69705050T2 (en) |
| DK (1) | DK0921450T3 (en) |
| ES (1) | ES2158472T3 (en) |
| GR (1) | GR3036373T3 (en) |
| NO (1) | NO20002591L (en) |
| PT (1) | PT921450E (en) |
| WO (1) | WO1999027422A1 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DK0921449T3 (en) * | 1997-11-20 | 2001-08-06 | Europ Economic Community | Computer controlled holographic method and device |
| US8941904B2 (en) * | 2005-07-04 | 2015-01-27 | Dai Nippon Printing Co., Ltd. | Hologram sheet and hologram observation sheet using same, and blinding device |
| JP4225346B2 (en) * | 2006-12-14 | 2009-02-18 | ソニー株式会社 | Playback device and playback method |
| JP4977886B2 (en) * | 2008-06-19 | 2012-07-18 | 大日本印刷株式会社 | Computer-generated hologram reproduction simulation method, computer-generated hologram production method using the reproduction simulation method, and computer-generated hologram produced using the computer-generated hologram production method |
| EP3175302B1 (en) | 2014-08-01 | 2021-12-29 | The Regents of the University of California | Device and method for iterative phase recovery based on pixel super-resolved on-chip holography |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5347375A (en) * | 1991-11-26 | 1994-09-13 | Kabushiki Kaisha Toshiba | Computer-assisted holographic image formation technique which determines interference pattern data used to form the holographic |
| US5668648A (en) * | 1991-11-26 | 1997-09-16 | Kabushiki Kaisha Toshiba | Computer-assisted holographic display apparatus |
-
1997
- 1997-11-20 DK DK97402798T patent/DK0921450T3/en active
- 1997-11-20 ES ES97402798T patent/ES2158472T3/en not_active Expired - Lifetime
- 1997-11-20 EP EP97402798A patent/EP0921450B1/en not_active Expired - Lifetime
- 1997-11-20 AT AT97402798T patent/ATE201775T1/en not_active IP Right Cessation
- 1997-11-20 DE DE69705050T patent/DE69705050T2/en not_active Expired - Lifetime
- 1997-11-20 PT PT97402798T patent/PT921450E/en unknown
-
1998
- 1998-11-17 CA CA002310851A patent/CA2310851A1/en not_active Abandoned
- 1998-11-17 WO PCT/EP1998/007539 patent/WO1999027422A1/en active Application Filing
- 1998-11-17 JP JP2000522500A patent/JP2001524697A/en active Pending
- 1998-11-17 US US09/530,967 patent/US6344909B1/en not_active Expired - Fee Related
-
2000
- 2000-05-19 NO NO20002591A patent/NO20002591L/en unknown
-
2001
- 2001-08-10 GR GR20010401229T patent/GR3036373T3/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001524697A (en) | 2001-12-04 |
| DE69705050T2 (en) | 2001-09-13 |
| NO20002591D0 (en) | 2000-05-19 |
| WO1999027422A1 (en) | 1999-06-03 |
| EP0921450A1 (en) | 1999-06-09 |
| NO20002591L (en) | 2000-07-19 |
| ES2158472T3 (en) | 2001-09-01 |
| DK0921450T3 (en) | 2001-09-17 |
| PT921450E (en) | 2001-11-30 |
| DE69705050D1 (en) | 2001-07-05 |
| ATE201775T1 (en) | 2001-06-15 |
| GR3036373T3 (en) | 2001-11-30 |
| US6344909B1 (en) | 2002-02-05 |
| EP0921450B1 (en) | 2001-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11460808B2 (en) | Method for generating a head up display for an aircraft using video holograms in real time with the help of sub-holograms | |
| US6621605B1 (en) | Computer-assisted method and device for restoring three-dimensional images | |
| Bryngdahl et al. | I Digital holography–computer-generated holograms | |
| US7636184B2 (en) | Method and device for computing computer-generated video holograms | |
| US7768684B2 (en) | 3D display | |
| Ferri | Visualization of 3D information with digital holography using laser printers | |
| US8218210B2 (en) | Method for generating computer-generated video holograms in real time by means of propagation | |
| US7277209B1 (en) | Computer-assisted holographic method and device | |
| Yoshikawa et al. | Integral holography | |
| US6344909B1 (en) | Method and device for hologram synthesis | |
| Lucente et al. | New approaches to holographic video | |
| Tsang et al. | Modern methods for fast generation of digital holograms | |
| Kondoh et al. | Hidden surface removal in full‐parallax CGHs by silhouette approximation | |
| Shiomi et al. | Fast hologram calculation method using wavelet transform: WASABI-2 | |
| EP1425634A1 (en) | Computation of computer generated holograms | |
| US7161721B2 (en) | Computer generated holograms | |
| WO2002039195A1 (en) | Improved method of producing a computer generated hologram | |
| US20230025687A1 (en) | Method for generating a head up display for an aircraft using video holograms in real time with the help of sub-holograms | |
| Bokor et al. | Computational method for testing computer-generated holograms | |
| Pappu et al. | Generalized pipeline for preview and rendering of synthetic holograms | |
| Rong et al. | Reconstruction efficiency enhancement of amplitude-type holograms by using Single-Scale Retinex algorithm | |
| CN115061348A (en) | Three-dimensional dynamic holographic display method and device based on object points | |
| Roux et al. | Intensity Distribution Transformation For General Beam-Shaping | |
| GB2271194A (en) | Holographic stereogram |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| FZDE | Discontinued |