CN113467211B - Holographic encoding method based on gradient reduction of spectral loss function - Google Patents

Holographic encoding method based on gradient reduction of spectral loss function Download PDF

Info

Publication number
CN113467211B
CN113467211B CN202110703230.4A CN202110703230A CN113467211B CN 113467211 B CN113467211 B CN 113467211B CN 202110703230 A CN202110703230 A CN 202110703230A CN 113467211 B CN113467211 B CN 113467211B
Authority
CN
China
Prior art keywords
hologram
frequency spectrum
phase
signal
spectrum
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.)
Active
Application number
CN202110703230.4A
Other languages
Chinese (zh)
Other versions
CN113467211A (en
Inventor
李海峰
张隽怿
刘旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202110703230.4A priority Critical patent/CN113467211B/en
Publication of CN113467211A publication Critical patent/CN113467211A/en
Application granted granted Critical
Publication of CN113467211B publication Critical patent/CN113467211B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/08Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0841Encoding method mapping the synthesized field into a restricted set of values representative of the modulator parameters, e.g. detour phase coding
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/08Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0841Encoding method mapping the synthesized field into a restricted set of values representative of the modulator parameters, e.g. detour phase coding
    • G03H2001/085Kinoform, i.e. phase only encoding wherein the computed field is processed into a distribution of phase differences

Abstract

The invention discloses a holographic coding method based on gradient decline of a spectrum loss function, which divides a frequency spectrum of a pure phase hologram into a signal frequency spectrum region and a noise frequency spectrum region, extracts the signal frequency spectrum region by filtering the frequency spectrum of a complex amplitude hologram, and then uses a gradient decline method to enable the signal frequency spectrum of the pure phase hologram to be approximately equal to the signal frequency spectrum of the complex amplitude hologram, wherein the noise frequency spectrum region of the pure phase hologram can be regarded as an optimized free variable and ignored. Compared with the common bi-phase encoding method and GS algorithm, the signal-to-noise ratio of the method is higher than that of the bi-phase encoding method, and the method is more suitable for holographic display of three-dimensional objects than the GS algorithm.

Description

Holographic encoding method based on gradient reduction of spectral loss function
Technical Field
The invention belongs to the field of three-dimensional display technologies, and particularly relates to a holographic encoding method based on gradient reduction of a spectrum loss function.
Background
Three-dimensional display technology is an important research and development direction of display technology. In recent years, with the popularization and application of virtual reality technology and augmented reality technology, research on head-mounted display devices has been increasing. Among them, holographic near-eye displays are receiving a lot of attention.
Holographic near-eye displays can carry depth information and have potential advantages in future near-eye display applications. To display the hologram, it is preferable to adjust the amplitude and phase of the light at the same time. But the spatial light modulator can only show amplitude or phase. A feasible method in experiments is to calculate the pure phase hologram that can be displayed in the phase-type spatial light modulator by using an encoding method. Currently, there are two commonly used holographic encoding algorithms: the double-phase encoding (double-phase) algorithm and the GS (Gerchberg-Saxton) algorithm.
Bi-phase encoding is the encoding of a complex amplitude value into two pure phase values. For the display of the three-dimensional object, firstly, a chromatography method or a point cloud method is used for firstly calculating the complex amplitude hologram of the three-dimensional object, and then a pure phase hologram is directly calculated by using double-phase encoding. However, bi-phase encoding suffers from coding noise, which causes the displayed image to appear noisy. The GS algorithm is an iterative method with the amplitude of the display image as the optimization target and the phase of the display image as the free variable. But such algorithms are difficult to compute holograms of three-dimensional objects. Moreover, when the phase of the displayed image is taken as a free variable, each pixel value of the actually reconstructed image is inevitably partially overlapped, and the phase values are different, so that interference occurs to cause noise.
Therefore, it is an important research direction in the field of holographic display to improve the display quality by proposing a new encoding method.
Disclosure of Invention
The invention aims to solve the problem of display quality caused by an encoding algorithm in a holographic near-eye and is suitable for displaying a three-dimensional object.
The purpose of the invention is realized by the following technical scheme: a holographic coding method based on gradient descent of a spectrum loss function divides the frequency spectrum of a pure phase hologram into a signal frequency spectrum region and a noise frequency spectrum region, extracts the signal frequency spectrum region by filtering the frequency spectrum of a complex amplitude hologram, and then enables the signal frequency spectrum of the pure phase hologram to be approximately equal to the signal frequency spectrum of the complex amplitude hologram by using a gradient descent method, wherein the noise frequency spectrum region of the pure phase hologram can be regarded as an optimized free variable and ignored; a filter is needed to be used for filtering the noise spectrum in an actual optical path; assuming a pure phase hologram, a complex amplitude hologram and an original size of m x n pixels (typically 1920 x 1080 pixels), the method comprises the steps of:
1) calculating the frequency spectrum F of a complex amplitude hologram by fast Fourier transform for a two-dimensional or three-dimensional original image light field0,F0Is a matrix of m x n, and F0Expanding the vector into an m x n dimensional vector;
2) filtering the frequency spectrum of the complex amplitude hologram, namely selecting a rectangular region with the length and the width of the center position of the frequency spectrum being half of the original size respectively, and extracting the rectangular region as a signal frequency spectrum region; representing the signal spectrum as
Figure BDA0003131049080000021
Wherein
Figure BDA0003131049080000022
Representing a Hadamard product, wherein P represents an m-n-dimensional vector, and if a certain element of P is in a signal spectrum region, the value of the certain element is 1, otherwise, the value of the certain element is 0;
3) for the initial pure phase hologram, performing fast Fourier transform to calculate a frequency spectrum F, selecting a rectangular region with a certain position length and width respectively being half of the original size as a signal frequency spectrum region, and taking the rest region as a noise frequency spectrum region; representing the signal spectrum as
Figure BDA0003131049080000023
4) The loss function is calculated for the signal spectrum of a phase-only hologram, since the spectrum is complex, the loss function C is defined as the form of the real part Re plus the imaginary part Im:
Figure BDA0003131049080000024
5) calculating the gradient of the loss function according to the loss function in the step 4), and updating the pure phase hologram by using the following gradient descent formula:
Figure BDA0003131049080000025
wherein theta iskA vector of m x n dimensions, C, which is expanded for the matrix formed by the phase values of the current phase-only hologramkFor the current loss function, η is the step size of the iteration, θk+1A vector of m x n dimensions is expanded for the matrix formed by the phase values of the updated phase-only hologram, k representing the kth iteration,
Figure BDA0003131049080000026
is a vector of dimensions m x n, the formula is:
Figure BDA0003131049080000027
wherein
Figure BDA0003131049080000031
Is a matrix with m x n rows and columns, F(k)Is the frequency spectrum F after the number k of iterations,
Figure BDA0003131049080000032
the FFT is a fast fourier transform, which is used to convert the optical field into a frequency spectrum,
Figure BDA0003131049080000033
is formed by thetakObtaining a m x n-dimensional complex vector;
6) and repeating the steps 4) and 5) for iteration until the loss function is converged to obtain the final pure phase hologram.
Furthermore, the amplitude value of the light field of the original image is obtained by calculating the gray value of the original image, and the phase value of the light field is set to be a constant value at the same depth, so that the noise caused by interference due to different phases of adjacent pixels of the reconstructed image can be avoided.
Further, the initial phase-only hologram may be generated using random values, or may be generated using a bi-phase hologram as the initial phase-only hologram, which may reduce the number of iterations.
Further, for practical purposes, the spatial light modulator used to display phase-only holograms can be formulated to have errors if crosstalk effects are present
Figure BDA0003131049080000034
Adding correction to crosstalk effect; the crosstalk effect can be expressed as a point spread model, and the phase values on the phase-only hologram can be rolled up into a point spread function in iteration; the formula for calculating the frequency spectrum can be rewritten as:
Figure BDA0003131049080000035
wherein
Figure BDA0003131049080000036
Representing a convolution, and a represents a point spread function of a pixel subject to crosstalk effects.
Further, if the size of the point spread function is in the order of one pixel, then the phase-only hologram needs to be oversampled.
Further, if off-axis holography is to be realized, a phase grating cannot be added in the last step like a common method, because the added phase grating is also influenced by the crosstalk effect, the position of a signal frequency spectrum region of the pure-phase hologram can be deviated from the center of a frequency spectrum in an iteration process instead of being arranged in the center of the frequency spectrum, and therefore off-axis holography is realized.
The invention has the beneficial effects that: compared with the common bi-phase coding method and the GS algorithm, the signal-to-noise ratio of the method is higher than that of the bi-phase coding method, and the method is more suitable for holographic display of three-dimensional objects than the GS algorithm.
Drawings
FIG. 1 is a schematic diagram of a signal spectrum region and a noise spectrum region according to an embodiment of the present invention, in which the abscissa and ordinate fx,fyWhen the spatial distance of the spectral plane in the actual experiment is shown, f is the focal length of the lens for Fourier transform, λ is the wavelength of light, and p is the pixel size of the spatial light modulator, the spectral region is
Figure BDA0003131049080000041
The spectral region of the signal can be selected as
Figure BDA0003131049080000042
Figure BDA0003131049080000043
Fig. 2 is a schematic diagram of simulation results of a reconstructed image of bi-phase encoding according to an embodiment of the present invention, (a) is an original image which is a two-dimensional picture of A, B, C letters, (b) is a simulation result of a reconstructed image of bi-phase encoding, and (c) is a simulation result of a reconstructed image of a method of the present invention;
fig. 3 is a schematic diagram showing simulation results of a reconstructed image by GS encoding according to an embodiment of the present invention, (a-c) shows an original image composed of 3-layer two-dimensional pictures, wherein the 1 st layer (a) is an alphabetic picture, the 2 nd and 3 rd layers (b-c) are random snowflake pictures, respectively, (d) shows simulation results of a reconstructed image by GS encoding, and (e) shows simulation results of a reconstructed image by the method of the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples.
Example 1
As shown in fig. 2, the present embodiment provides a comparison with the simulation result of a bi-phase encoded reconstructed image, where (a) is the original image, which is a two-dimensional picture of A, B, C letters. (b) The peak signal-to-noise ratio PSNR is 20.579, which is a simulation result of a reconstructed image of bi-phase encoding. (c) Is the simulation result of the reconstruction graph of the method of the invention, and the peak signal-to-noise ratio PSNR is 24.236. The initial pure phase hologram of the method of the invention adopts a biphase encoded hologram, and the total iteration times are 5 times. The light intensity and gray level conversion formula of holographic coding and reconstruction simulation uses:
Figure BDA0003131049080000044
Figure BDA0003131049080000045
where I is the normalized light intensity and gray is the normalized gray level. The diffraction distance between the reconstructed image plane and the holographic plane is 80 mm. The observation and simulation results show that the method has higher display quality and less noise, and the signal-to-noise ratio of the method is higher compared with that of the double-phase coding.
Example 2
As shown in fig. 3, this example provides a comparison with the simulation result of the reconstructed image of the GS algorithm, where (a-c) is an original, the original is composed of 3-layer two-dimensional pictures, the 1 st layer (a) is an alphabetic picture, the diffraction distance from the hologram surface is 80mm, the 2 nd and 3 rd layers (b-c) are random snowflake pictures, respectively, (b) the diffraction distance from the hologram surface is 81.5mm, and (c) the diffraction distance from the hologram surface is 83 mm. (d) The simulated result of the reconstructed image encoded by GS showed a diffraction distance of 80mm from the hologram surface. (e) The simulation result of the reconstructed image of the method of the invention is that the diffraction distance from the holographic surface is 80 mm. Because the GS algorithm can only optimize two-dimensional original images, 3 holograms of the GS algorithm are respectively calculated for 3 layers of original images, and then the simulated 3 reconstructed images are directly superposed, which is equivalent to simulating a time division multiplexing multi-layer GS algorithm. The GS algorithm uses a total number of iterations of 30. The initial pure phase hologram of the method of the invention adopts a biphase encoded hologram, and the total iteration times are 5 times.
It can be seen that the defocus patterns of the second and third layers of snowflakes in the GS-encoded reconstructed image are very noisy, because the GS-encoded reconstructed image has the phase of the reconstructed image as a free variable, and when the image is out of focus, adjacent pixels are overlapped and interfered, and the phase difference causes noise. The snowflake defocusing patterns in the simulation result of the method have no obvious noise, and the phase of the original image can be set to be a constant value, so that the obvious noise can not appear in a certain defocusing range. Therefore, the method is more suitable for displaying the three-dimensional object than the GS algorithm.
The above-described embodiments are intended to illustrate rather than to limit the invention, and any modifications and variations of the present invention are within the spirit of the invention and the scope of the appended claims.

Claims (6)

1. A holographic coding method based on gradient descent of a spectrum loss function is characterized in that a frequency spectrum of a pure phase hologram is divided into a signal frequency spectrum region and a noise frequency spectrum region, the signal frequency spectrum region is extracted by filtering the frequency spectrum of a complex amplitude hologram, then the signal frequency spectrum of the pure phase hologram is approximately equal to the signal frequency spectrum of the complex amplitude hologram by a gradient descent method, and the noise frequency spectrum region of the pure phase hologram is regarded as an optimized free variable and ignored; setting the size of pure phase hologram, complex amplitude hologram and original image as m pixel, the method includes the following steps:
1) calculating the frequency spectrum F of the complex amplitude hologram by using fast Fourier transform for the two-dimensional or three-dimensional original image light field0,F0Is a matrix of m x n, and F0Expanding the vector into an m x n dimensional vector;
2) filtering the frequency spectrum of the complex amplitude hologram, namely selecting a rectangular region with the length and the width of the center position of the frequency spectrum being respectively half of the original size, and extracting the rectangular region as a signal frequency spectrum region; representing the signal spectrum as
Figure FDA0003537797690000015
Wherein
Figure FDA0003537797690000016
Representing a Hadamard product, wherein P represents an m-n-dimensional vector, and if a certain element of P is in a signal spectrum region, the value of the certain element is 1, otherwise, the value of the certain element is 0;
3) for the initial pure phase hologram, performing fast Fourier transform to calculate a frequency spectrum F, selecting a rectangular area with the length and width of a certain position being half of the original size as a signal frequency spectrum area, and taking the rest area as a noise frequency spectrum area; representing the signal spectrum as
Figure FDA0003537797690000017
4) The loss function is calculated for the signal spectrum of a phase-only hologram, since the spectrum is complex, the loss function C is defined as the form of the real part Re plus the imaginary part Im:
Figure FDA0003537797690000011
5) calculating the gradient of the loss function according to the loss function in the step 4), and updating the pure phase hologram by using the following gradient descent formula:
Figure FDA0003537797690000012
wherein theta iskA vector of m x n dimensions, C, which is expanded for the matrix formed by the phase values of the current phase-only hologramkFor the current loss function, η is the step size of the iteration, θk+1A vector of m x n dimensions is expanded for the matrix formed by the phase values of the updated phase-only hologram, k representing the kth iteration,
Figure FDA0003537797690000013
is a vector with dimensions of m x n, and the formula is as follows:
Figure FDA0003537797690000014
wherein
Figure FDA0003537797690000021
Is a matrix with m x n rows and columns, F(k)Is the frequency spectrum F after the number k of iterations,
Figure FDA0003537797690000022
the FFT is a fast fourier transform, which is used to convert the optical field into a frequency spectrum,
Figure FDA0003537797690000023
is formed by thetakObtaining a m x n-dimensional complex vector;
6) and repeating the steps 4) and 5) for iteration until the loss function is converged to obtain the final pure phase hologram.
2. The method of claim 1, wherein the amplitude value of the light field of the original image is calculated by using the gray value of the original image, and the phase value of the light field is set to be constant at the same depth, so as to avoid the noise caused by the interference of adjacent pixels of the reconstructed image due to different phases.
3. The method of claim 1, wherein the initial phase-only hologram or bi-phase hologram is generated using random values as the initial phase-only hologram, and wherein the number of iterations is reduced using the bi-phase hologram.
4. Method according to claim 1, characterized in that for practical situations the spatial light modulator used for displaying the phase-only hologram will cause display errors if there are crosstalk effects, in the formula
Figure FDA0003537797690000024
Adding correction to crosstalk effect; the crosstalk effect is expressed as a point diffusion model, and the phase value of the phase-only hologram is rolled up into a point diffusion function in iteration; the formula for the calculation of the spectrum is rewritten as:
Figure FDA0003537797690000025
wherein
Figure FDA0003537797690000026
Representing a convolution, and a represents a point spread function of a pixel subject to crosstalk effects.
5. The method according to claim 4, characterized in that if the size of the point spread function is in the order of one pixel, a supersampling calculation of the phase-only hologram is required.
6. Method according to claim 4, characterized in that the location of the signal spectral region of the phase-only hologram is displaced in an iterative process not in the center of the spectrum but by a distance, so that off-axis holography is achieved.
CN202110703230.4A 2021-06-24 2021-06-24 Holographic encoding method based on gradient reduction of spectral loss function Active CN113467211B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110703230.4A CN113467211B (en) 2021-06-24 2021-06-24 Holographic encoding method based on gradient reduction of spectral loss function

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110703230.4A CN113467211B (en) 2021-06-24 2021-06-24 Holographic encoding method based on gradient reduction of spectral loss function

Publications (2)

Publication Number Publication Date
CN113467211A CN113467211A (en) 2021-10-01
CN113467211B true CN113467211B (en) 2022-06-10

Family

ID=77872569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110703230.4A Active CN113467211B (en) 2021-06-24 2021-06-24 Holographic encoding method based on gradient reduction of spectral loss function

Country Status (1)

Country Link
CN (1) CN113467211B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114137812B (en) * 2021-12-03 2022-10-11 北京理工大学 Triple-correlation hybrid holographic encryption method based on metasurface
CN115097709B (en) * 2022-07-05 2023-11-17 东南大学 Holographic coding method based on complex optimizer or complex solver

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738727A (en) * 2019-09-29 2020-01-31 东南大学 complex amplitude hologram algorithm based on light field
CN210666315U (en) * 2019-05-16 2020-06-02 安徽大学 Multi-view three-dimensional display device based on holographic projection
CN111897197A (en) * 2020-08-18 2020-11-06 四川大学 Fourier phase hologram generation method based on double-phase encoding

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10355648A1 (en) * 2003-11-28 2005-06-23 Tesa Scribos Gmbh Computing and generating computer-generated hologram involves using generalized Gerchberg-Saxton algorithm, determining first section of reconstruction independently of second section of reconstruction
CN103645562B (en) * 2013-12-20 2015-10-21 中南大学 The beam shaping method of modulated amplitude and phase place while of a kind of
WO2018021984A1 (en) * 2016-07-28 2018-02-01 Cy Vision Inc. System and method for high-quality speckle-free phase-only computer-generated holographic image projection
CN106210707B (en) * 2016-09-14 2018-05-01 京东方科技集团股份有限公司 A kind of holographic display system and holographic display methods
CN106842880B (en) * 2017-03-27 2018-09-28 深圳市美誉镜界光电科技有限公司 Hologram image generation method, processor and image hologram display device, equipment
GB2579234B (en) * 2018-11-27 2023-07-19 Dualitas Ltd Hologram calculation
US11137719B2 (en) * 2018-12-11 2021-10-05 University Of North Carolina At Chapel Hill Methods, systems, and computer readable media for improved digital holography and display incorporating same

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN210666315U (en) * 2019-05-16 2020-06-02 安徽大学 Multi-view three-dimensional display device based on holographic projection
CN110738727A (en) * 2019-09-29 2020-01-31 东南大学 complex amplitude hologram algorithm based on light field
CN111897197A (en) * 2020-08-18 2020-11-06 四川大学 Fourier phase hologram generation method based on double-phase encoding

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三维纯相位全息显示中的散斑噪声抑制;周婷婷,韩超;《中国激光》;20181231;第45卷(第12期);全文 *
基于空间光调制器的三维相位全息图再现像优化;梁浩聪,蒋晓瑜,牛树来,闫兴鹏,赵锴;《光学技术》;20140531;第40卷(第3期);全文 *

Also Published As

Publication number Publication date
CN113467211A (en) 2021-10-01

Similar Documents

Publication Publication Date Title
Blinder et al. The state-of-the-art in computer generated holography for 3D display
US11137719B2 (en) Methods, systems, and computer readable media for improved digital holography and display incorporating same
CN113467211B (en) Holographic encoding method based on gradient reduction of spectral loss function
TWI409716B (en) A computer device for generating a video image and for expanding Real - time generation of image - like image of 3D rendering drawing pipeline
JP5266223B2 (en) Method for generating computer video holograms in real time using propagation
Liu et al. 4K-DMDNet: diffraction model-driven network for 4K computer-generated holography
US20230205133A1 (en) Real-time Photorealistic 3D Holography With Deep Neural Networks
CN113986163A (en) Method and device for encoding a complex-valued signal for reconstructing a three-dimensional object
KR20140096532A (en) Apparatus and method for generating digital hologram
CN111240177B (en) Holographic speckle noise suppression method based on layered pixel scanning algorithm
CN114387395A (en) Phase-double resolution ratio network-based quick hologram generation method
KR100910642B1 (en) Method for reproducing hologram 3D image by integral imaging scheme and Apparatus thereof
Min et al. Wigner inverse transform based computer generated hologram for large object at far field from its perspective light field
CN104090476A (en) Three-dimensional scene information acquisition method for holographic display
Shiomi et al. Fast hologram calculation method using wavelet transform: WASABI-2
CN112765624B (en) Authenticatable phase-only hologram generation method based on phase optimization and sparse constraint
Ryu et al. Continuous depth control of phase-only hologram with depth embedding block
Park et al. Parallel Synthesis Algorithm for Layer-based Computer-generated Holograms Using Sparse-field Localization
CN115797231A (en) Real-time hologram generation method based on neural network of Fourier inspiration
JP7486381B2 (en) Hologram data generating device and program thereof
CN114764220B (en) Method for improving speckle autocorrelation reconstruction effect based on off-axis digital holography
Dong et al. Vision transformer-based, high-fidelity, computer-generated holography
Zhou et al. Point spread function-inspired deformable convolutional network for holographic displays
Shigematsu et al. Computer-generated holography with ordinary display
Jiao et al. Fast phase-only hologram calculation with a shared phase mask optimized from an image dictionary

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant