CN113448233A - Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system - Google Patents

Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system Download PDF

Info

Publication number
CN113448233A
CN113448233A CN202110790766.4A CN202110790766A CN113448233A CN 113448233 A CN113448233 A CN 113448233A CN 202110790766 A CN202110790766 A CN 202110790766A CN 113448233 A CN113448233 A CN 113448233A
Authority
CN
China
Prior art keywords
focusing
hologram
reconstruction
sampling
algorithm
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.)
Granted
Application number
CN202110790766.4A
Other languages
Chinese (zh)
Other versions
CN113448233B (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.)
Anhui University
Original Assignee
Anhui University
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 Anhui University filed Critical Anhui University
Priority to CN202110790766.4A priority Critical patent/CN113448233B/en
Publication of CN113448233A publication Critical patent/CN113448233A/en
Application granted granted Critical
Publication of CN113448233B publication Critical patent/CN113448233B/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/0866Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • 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/0443Digital holography, i.e. recording holograms with digital recording means
    • 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/0866Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G03H2001/0883Reconstruction aspect, e.g. numerical focusing

Abstract

The invention discloses an under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and a system, belonging to the technical field of digital holography, comprising the steps of respectively carrying out down-sampling operation on holograms generated by irradiating different positions of an object in an optical structure by a light source to obtain down-sampling holograms at different positions; reconstructing the down-sampling holograms at different positions based on a CS algorithm to obtain amplitude images at different positions; estimating the position of a focal plane according to amplitude images at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm; and reconstructing the down-sampled hologram at the position of the focusing plane based on the TwinT algorithm. The invention provides an EIG-AF-CS self-focusing algorithm for realizing self-focusing under the condition of undersampling, combines the EIG-AF-CS self-focusing algorithm with compression reconstruction without twin images, and effectively overcomes the problems of twin image interference and lack of self-focusing capability in the traditional compression holographic reconstruction while ensuring the reconstruction quality.

Description

Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system
Technical Field
The invention relates to the technical field of digital holography, in particular to an under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system.
Background
Digital Holography (DH) is a holographic technique in which an intensity image is recorded by a CCD, and then the original object information is reconstructed by computer numerical calculations. Its disadvantage is the high storage capacity and bandwidth required, which greatly limits the wide-spread use of digital holography. In 2006, Cand, Donoho, Tao et al proposed the Compressed Sensing (CS) theory, which indicated that if the signal is sparse or compressible, only limited measurements of the sparse signal are needed to reconstruct the original signal accurately, saving resources and time greatly, breaking through the limitation that the sampling rate specified by the traditional Nyquist sampling theorem must be at least twice the signal bandwidth.
In 2009, Brady et al combine compressed sensing and digital holography together, propose a compressed holographic method, build a bridge between the compressed sensing theory and the holographic display technology, and by using the compressed sensing theory, an original object can be accurately reconstructed from a 2D hologram by a small amount of measurement sample data, and a Brady group experiment realizes holographic layer chromatography reconstruction of two layers of dandelions with a long distance; subsequently, the group also successfully completed several studies, such as Compressive Holographic Tomography (CHT), Millimeter Wave (MWH), and blur (DH) holograms. In addition, in 2014, Stern et al in israel proposed a Compressive Fresnel Holography (CFH) based on a variable down-sampling strategy, and provided conditions required to be satisfied for realizing the actual Compressive Fresnel Holography, and then the extension of the CFH is applied to imaging fields such as multi-dimensional imaging, tomography, reconstruction of partially-occluded objects, and the like.
The traditional compression holographic reconstruction method mainly faces the interference of a twin image and a zero-order image, and the twin image and the zero-order image are usually eliminated by adopting a multi-exposure method. In 2018, Zhang et al, university of Qinghua, based on the sparsity difference between the reconstructed image and the twin image caused by the physical symmetry of the hologram, proposed a CS method to successfully realize the twin-image-free holographic reconstruction from the single-exposure hologram.
The method proposed by Zhang et al can realize twinning-free image reconstruction, but only aims at a 2D-2D full-sampling imaging model, lacks a compression sampling process, and has no self-focusing capability because the position of a focal plane must be known during reconstruction.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and realize the twinning-free image self-focusing compression holographic reconstruction from a single undersampled hologram.
To achieve the above object, in one aspect, the present invention employs an under-sampled hologram compression holographic multi-scale self-focusing reconstruction method, including:
carrying out down-sampling operation on the hologram generated by the object through the illumination of the light source to obtain a down-sampling hologram;
and reconstructing the down-sampling hologram based on a TwinST algorithm to obtain a reconstruction result.
Further, the down-sampling operation is performed on the hologram generated by the object illuminated by the light source, so as to obtain a down-sampled hologram, and the formula is as follows:
Figure BDA0003160796980000021
wherein D isΩRepresenting down-sampling, G representing forward transformation,. phiΩRepresenting the term 2Re { DΩG ρ } versus ρ, ρ representing the object density, e representing the model error.
Further, the down-sampling operation is performed on the hologram generated by the object illuminated by the light source to obtain a down-sampled hologram, and the method further includes:
and respectively carrying out down-sampling operation on the holograms generated by the irradiation of the light source at different positions of the object in the optical structure to obtain down-sampled holograms at different positions.
Further, reconstructing the down-sampling hologram based on the TwIST algorithm to obtain a reconstruction result, further comprising:
estimating the position of a focal plane by adopting an evaluation method of a characteristic value algorithm for the downsampled holograms at different positions;
and reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
Further, the estimating the focal plane position of the downsampled hologram at different positions by using an evaluation method of a characteristic value algorithm includes:
reconstructing the down-sampling holograms at different positions based on a CS algorithm to obtain amplitude images at different positions;
and estimating the position of the focal plane according to the amplitude images at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm.
Further, still include: and estimating the position of the focal plane by combining a multi-scale search strategy and adopting an evaluation method of a characteristic value self-focusing algorithm on the downsampled holograms at different positions, wherein the method specifically comprises the following steps:
1) setting a search start point of a first layer scale
Figure BDA0003160796980000031
Searching for an endpoint
Figure BDA0003160796980000032
And search spacing
Figure BDA0003160796980000033
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
2) reconstructing the down-sampling hologram at each position in the current layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
3) the evaluation method adopting the characteristic value self-focusing algorithm is adopted according to different conditionsThe amplitude image of the position is estimated to obtain the search result of the current layer scale
Figure BDA0003160796980000034
4) Taking the search result of the current layer scale as the search result of the previous layer scale, and taking the search result of the previous layer scale as the search result of the previous layer scale
Figure BDA0003160796980000035
As a central point, based on a relational expression
Figure BDA0003160796980000036
i 2,3, …, k obtaining the search starting point of the next layer search space
Figure BDA0003160796980000037
Searching for an endpoint
Figure BDA0003160796980000038
And search spacing
Figure BDA0003160796980000039
And the next layer position sequence, and repeatedly executing the steps 2) -4) until the self-focusing search result of the k-th layer scale is obtained
Figure BDA00031607969800000310
5) Reconstructing a location based on the TwinT algorithm
Figure BDA00031607969800000311
And (4) processing the down-sampling hologram to obtain the reconstruction result.
In another aspect, an under-sampled hologram compression holographic multi-scale self-focusing reconstruction system is adopted, which comprises a down-sampling operation module and a reconstruction module, wherein:
the down-sampling operation module is used for performing down-sampling operation on the hologram generated by the object through illumination of the light source to obtain a down-sampling hologram;
and the reconstruction module is used for reconstructing the down-sampling hologram based on a TwinT algorithm to obtain a reconstruction result.
Furthermore, the down-sampling operation module is further configured to perform down-sampling operations on holograms generated by the illumination of the object by the light source at different positions in the optical structure, so as to obtain down-sampled holograms at different positions.
Further, the reconstruction module further comprises a focal plane position estimation unit and a reconstruction unit, wherein:
the focal plane position estimation unit is used for estimating the focal plane position of the down-sampled holograms at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm;
the reconstruction unit is used for reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
Further, the focal plane position estimation unit is configured to estimate the focal plane position by using an evaluation method of a feature value self-focusing algorithm for the downsampled holograms at different positions in combination with a multi-scale search strategy, and specifically includes:
1) setting a search start point of a first layer scale
Figure BDA0003160796980000041
Searching for an endpoint
Figure BDA0003160796980000042
And search spacing
Figure BDA0003160796980000043
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
2) reconstructing the down-sampling hologram at each position in the current layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
3) the evaluation method of the characteristic value self-focusing algorithm is adopted to estimate and obtain the search result of the current layer scale according to the amplitude images at different positions
Figure BDA0003160796980000044
4) Taking the search result of the current layer scale as the search result of the previous layer scale, and taking the search result of the previous layer scale as the search result of the previous layer scale
Figure BDA0003160796980000045
As a central point, based on a relational expression
Figure BDA0003160796980000051
i 2,3, …, k obtaining the search starting point of the next layer search space
Figure BDA0003160796980000052
Searching for an endpoint
Figure BDA0003160796980000053
And search spacing
Figure BDA0003160796980000054
And the next layer position sequence, and repeatedly executing the steps 2) -4) until the self-focusing search result of the k-th layer scale is obtained
Figure BDA0003160796980000055
5) Reconstructing a location based on the TwinT algorithm
Figure BDA0003160796980000056
And (4) processing the down-sampling hologram to obtain the reconstruction result.
Compared with the prior art, the invention has the following technical effects: aiming at the lack of a compression sampling process in the traditional twin image-free reconstruction process, the invention carries out down-sampling on an under-sampled hologram to supplement the compression sampling process and establishes a compression holographic imaging model based on the under-sampled hologram; and aiming at the problem that accurate self-focusing is difficult to realize after the down-sampling operation is added, an EIG-AF-CS self-focusing algorithm is provided and is combined with the twin-free image compression reconstruction. The reconstruction can be realized by using less original data through down-sampling the hologram, the reconstruction quality can be ensured, the self-focusing reconstruction can be realized by combining an EIG-AF-CS self-focusing algorithm, and the problems of twin image interference and lack of self-focusing capability in the traditional compressed holographic reconstruction are effectively solved.
Drawings
The following detailed description of embodiments of the invention refers to the accompanying drawings in which:
FIG. 1 is a flow chart of a method of compressed holographic reconstruction of an undersampled hologram;
FIG. 2 is a schematic down-sampling of an under-sampled hologram, wherein (a) the down-sampling of the under-sampled hologram is added and (b) the down-sampled hologram is up-sampled prior to two-way propagation;
FIG. 3 is a flow chart of a self-focusing reconstruction method for under-sampled hologram feature values;
FIG. 4 is an overall flow chart of an under-sampled hologram feature value self-focusing reconstruction method;
FIG. 5 is a flow chart of a method for multi-scale self-focusing reconstruction of an undersampled hologram.
Detailed Description
To further illustrate the features of the present invention, refer to the following detailed description of the invention and the accompanying drawings. The drawings are for reference and illustration purposes only and are not intended to limit the scope of the present disclosure.
As shown in fig. 1 to fig. 2, the present embodiment discloses a method for reconstructing an under-sampled hologram by compressing a hologram in a multi-scale self-focusing manner, which includes the following steps S1 to S2:
s1, carrying out down-sampling operation on the hologram generated by the object through the illumination of the light source to obtain a down-sampled hologram;
s2, reconstructing the down-sampling hologram based on the TwinT algorithm to obtain a reconstruction result.
It should be noted that, for the problem that the existing twinning-free holography method (TIFH-CS) lacks a compressive sampling process and does not have a self-focusing capability, the present embodiment overcomes the disadvantage that the existing twinning-free holography method lacks a compressive sampling process, as shown in fig. 2, a downsampling operation formula is shown as follows:
Figure BDA0003160796980000061
wherein D isΩRepresenting down-sampling, e.g. uniform random down-sampling, G representing forward transformation, phiΩRepresenting the term 2Re { DΩG ρ } versus ρ, ρ representing the object density, e representing the model error,
Figure BDA0003160796980000062
representing a down-sampled hologram.
As a further preferred solution, the reconstruction needs to be known from
Figure BDA0003160796980000063
To solve ρ, this embodiment adopts a Two-step Iterative shrinkage/Threshold (TwIST-step Iterative shrinkage/Threshold) algorithm to solve the problem:
Figure BDA0003160796980000064
wherein the content of the first and second substances,
Figure BDA0003160796980000065
represents the density of the object finally optimized by the TwinT algorithm, i.e. the result of the optimization of the TwinT algorithm, and tau represents the residual error l2Relative weight between norm and estimated TV norm.
It should be noted that, after down-sampling is added to the hologram in the twin image free holography method, it is difficult to achieve accurate self-focusing, and the evaluation method of the Eigenvalue (EIG) self-focusing algorithm of this embodiment is improved, and a TwIST algorithm is used instead of back propagation to reconstruct the down-sampled holograms at different positions to obtain amplitude images at different positions, so as to achieve the self-focusing of the twin image free holography. As shown in fig. 3 to 4, in order to implement the self-focusing function, the present embodiment provides an under-sampled hologram eigenvalue self-focusing reconstruction method (EIG-AF-CS), which includes the following steps:
1) respectively carrying out down-sampling operation on holograms generated by the irradiation of the light source at different positions of the object in the optical structure to obtain down-sampled holograms at different positions;
2) reconstructing the down-sampling holograms at different positions based on a CS algorithm to obtain amplitude images at different positions;
3) estimating the position of a focal plane according to amplitude images at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm;
4) and reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
As a further preferred technical solution, the estimating method using the EIG algorithm estimates the position of the focal plane according to the amplitude images at different positions, specifically: from a distance z1=zsStarting to perform a hologram reconstruction with a reconstruction distance interval zstepLet M × N and NzRepresenting the size of the hologram and the total number of reconstructed images, respectively, using the set of reconstructed images to perform automatic focal plane detection, comprising the steps of:
at a given reconstruction distance z, the corresponding amplitude image Az(x, y) is normalized by its energy as:
Figure BDA0003160796980000071
the average of the normalized amplitude images is calculated as:
Figure BDA0003160796980000072
covariance matrix QzComprises the following steps:
Figure BDA0003160796980000073
wherein the content of the first and second substances,
Figure BDA0003160796980000074
(·)Trepresenting a matrix transpose, comprising QzColumns of characteristic valuesVector EzComprises the following steps:
Ez=eig(Qz)
where the eigenvalue calculation function in MATLAB is represented. Taking into account the eigenvalues in the vector EzIn ascending order, the metric L for detecting the focal plane is defined as:
Figure BDA0003160796980000081
where κ is a parameter that may be customized.
In the present case, the covariance matrix QzContains most of the information about spatial variations in the reconstructed image. Therefore, all these main eigenvalues should be utilized in the calculation of the metric L. However, the degree of spatial variation in the reconstructed image is almost monotonic with reconstruction distance, except in the focal plane. To compensate for this effect, a certain number of dominant eigenvalues need to be discarded in the computation of the autofocus metric. This number is given by the kappa parameter. In general, a suitable value of κ depends on the size of the reconstructed image. Due to QzThere are M eigenvalues, so the k value can be chosen as a percentage of M. Finally, z ═ z when the selection metric L takes the maximum/minimum valuefThe value serves as the position of the focal plane.
As shown in fig. 5, in order to improve the accuracy of self-focusing, the present embodiment combines a multi-scale search strategy to position a focal plane layer by layer, and implements multi-scale self-focusing (multi-scale EIG-AF-CS) under the condition of undersampling, which specifically includes the following steps:
1) respectively carrying out down-sampling operation on holograms generated by the irradiation of the light source at different positions of the object in the optical structure to obtain down-sampled holograms at different positions;
2) setting a search start point of a first layer scale
Figure BDA0003160796980000082
Searching for an endpoint
Figure BDA0003160796980000083
And search spacing
Figure BDA0003160796980000084
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
3) reconstructing the down-sampling hologram at each position in the first layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
4) evaluation method using the eigenvalue autofocus algorithm
Figure BDA0003160796980000085
Estimating and obtaining a search result of a first-layer scale according to amplitude images at different positions
Figure BDA0003160796980000086
5) Search results at a first hierarchical scale
Figure BDA0003160796980000087
As a central point, based on a relational expression
Figure BDA0003160796980000091
i
2,3, …, k obtaining the search starting point of the next layer search space
Figure BDA0003160796980000092
Searching for an endpoint
Figure BDA0003160796980000093
And search spacing
Figure BDA0003160796980000094
And a second layer position sequence, which reduces the search space layer by layer until the self-focusing search result of the k-th layer scale is obtained
Figure BDA0003160796980000095
6) Reconstructing a location based on the TwinT algorithm
Figure BDA0003160796980000096
And (4) processing the down-sampling hologram to obtain the reconstruction result.
The embodiment provides an EIG-AF-CS self-focusing algorithm for realizing self-focusing under the undersampling condition, and combines the EIG-AF-CS self-focusing algorithm with twin-free image compression reconstruction. Experiments prove that reconstruction can be realized by using less original data for down-sampling the hologram, the reconstruction quality can be ensured, self-focusing reconstruction can be realized by combining an EIG-AF-CS self-focusing algorithm, and the problems of twin image interference and lack of self-focusing capability in the traditional compressed holographic reconstruction are effectively solved.
The embodiment also discloses an under-sampling hologram compression holographic multi-scale self-focusing reconstruction system, which is characterized by comprising a down-sampling operation module and a reconstruction module, wherein:
the down-sampling operation module is used for performing down-sampling operation on the hologram generated by the object through illumination of the light source to obtain a down-sampling hologram;
and the reconstruction module is used for reconstructing the down-sampling hologram based on a TwinT algorithm to obtain a reconstruction result.
As a further preferable technical solution, the down-sampling operation module is further configured to perform down-sampling operations on holograms generated by the illumination of the object by the light source at different positions in the optical structure, respectively, so as to obtain down-sampled holograms at different positions.
The reconstruction module further comprises a focal plane position estimation unit and a reconstruction unit, wherein:
the focal plane position estimation unit is used for estimating the focal plane position of the down-sampled holograms at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm;
the reconstruction unit is used for reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
As a further preferred technical solution, the focal plane position estimating unit is configured to estimate the focal plane position by using an evaluation method of an EIG algorithm for the downsampled holograms at different positions in combination with a multi-scale search strategy, and specifically includes:
1) setting a search start point of a first layer scale
Figure BDA0003160796980000101
Search endpoint
Figure BDA0003160796980000102
And search spacing
Figure BDA0003160796980000103
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
2) reconstructing the down-sampling hologram at each position in the current layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
3) the evaluation method of the characteristic value self-focusing algorithm is adopted to estimate and obtain the search result of the current layer scale according to the amplitude images at different positions
Figure BDA0003160796980000104
4) Taking the search result of the current layer scale as the search result of the previous layer scale, and taking the search result of the previous layer scale as the search result of the previous layer scale
Figure BDA0003160796980000105
As a central point, based on a relational expression
Figure BDA0003160796980000106
i
2,3, …, k obtaining the search starting point of the next layer search space
Figure BDA0003160796980000107
Searching for an endpoint
Figure BDA0003160796980000108
And search spacing
Figure BDA0003160796980000109
And the next layer position sequence is repeatedly executedStep 2) -4) until the self-focusing search result of the k-th layer scale is obtained
Figure BDA00031607969800001010
5) Reconstructing a location based on the TwinT algorithm
Figure BDA00031607969800001011
And (4) processing the down-sampling hologram to obtain the reconstruction result.
The system provided by the embodiment of the present invention is used for executing the above method embodiments, and for details of the process and the details, reference is made to the above embodiments, which are not described herein again.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (10)

1. An under-sampled hologram compression holographic multi-scale self-focusing reconstruction method is characterized by comprising the following steps:
carrying out down-sampling operation on the hologram generated by the object through the illumination of the light source to obtain a down-sampling hologram;
and reconstructing the down-sampling hologram based on a TwinST algorithm to obtain a reconstruction result.
2. The method for reconstructing an under-sampled hologram by compressing a holographic multi-scale self-focusing according to claim 1, wherein the hologram generated by the object illuminated by the light source is down-sampled to obtain a down-sampled hologram, and the formula is as follows:
Figure FDA0003160796970000011
wherein D isΩRepresenting down-sampling, G representing forward transformation,. phiΩRepresenting the term 2Re { DΩG rho to rho, where rho represents object density and e represents modelAnd (4) error.
3. The method for compressed holographic multi-scale self-focusing reconstruction according to claim 1, wherein the down-sampling operation is performed on the hologram generated by the object illuminated by the light source to obtain the down-sampled hologram, further comprising:
and respectively carrying out down-sampling operation on the holograms generated by the irradiation of the light source at different positions of the object in the optical structure to obtain down-sampled holograms at different positions.
4. The method for reconstructing compressed holographic multi-scale self-focusing of under-sampled hologram according to claim 3, wherein said reconstructing said under-sampled hologram based on TwinST algorithm to obtain a reconstruction result further comprises:
estimating the position of a focal plane by adopting an evaluation method of a characteristic value self-focusing algorithm for the downsampling holograms at different positions;
and reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
5. The method for compressed holographic multi-scale self-focusing reconstruction according to claim 4, wherein the estimating the focal plane position of the down-sampled hologram at the different positions by using an evaluation method of a feature value algorithm comprises:
reconstructing the down-sampling holograms at different positions based on a CS algorithm to obtain amplitude images at different positions;
and estimating the position of the focal plane according to the amplitude images at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm.
6. The method of compressed holographic multi-scale self-focusing reconstruction of an undersampled hologram according to claim 4, further comprising: and estimating the position of the focal plane by combining a multi-scale search strategy and adopting an evaluation method of a characteristic value self-focusing algorithm on the downsampled holograms at different positions, wherein the method specifically comprises the following steps:
1) setting a search start point of a first layer scale
Figure FDA0003160796970000021
Searching for an endpoint
Figure FDA0003160796970000022
And search spacing
Figure FDA0003160796970000023
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
2) reconstructing the down-sampling hologram at each position in the current layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
3) the evaluation method of the characteristic value self-focusing algorithm is adopted to estimate and obtain the search result of the current layer scale according to the amplitude images at different positions
Figure FDA0003160796970000024
4) Taking the search result of the current layer scale as the search result of the previous layer scale, and taking the search result of the previous layer scale as the search result of the previous layer scale
Figure FDA0003160796970000025
As a central point, based on a relational expression
Figure FDA0003160796970000026
Obtaining a search starting point of a next search space
Figure FDA0003160796970000027
Searching for an endpoint
Figure FDA0003160796970000028
And search spacing
Figure FDA0003160796970000029
And the next layer position sequence, and repeatedly executing the steps 2) -4) until the self-focusing search result of the k-th layer scale is obtained
Figure FDA00031607969700000210
5) Reconstructing a location based on the TwinT algorithm
Figure FDA00031607969700000211
And (4) processing the down-sampling hologram to obtain the reconstruction result.
7. An under-sampled hologram compression holographic multi-scale self-focusing reconstruction system, comprising a down-sampling operation module and a reconstruction module, wherein:
the down-sampling operation module is used for performing down-sampling operation on the hologram generated by the object through illumination of the light source to obtain a down-sampling hologram;
and the reconstruction module is used for reconstructing the down-sampling hologram based on a TwinT algorithm to obtain a reconstruction result.
8. The under-sampled hologram compression holographic multi-scale self-focusing reconstruction system according to claim 7, wherein the down-sampling operation module is further configured to perform down-sampling operations on the holograms of the object illuminated by the light source at different positions in the optical structure, respectively, so as to obtain the down-sampled holograms at different positions.
9. The under-sampled hologram compression holographic multi-scale self-focusing reconstruction system of claim 8, wherein the reconstruction module further comprises a focal plane position estimation unit and a reconstruction unit, wherein:
the focal plane position estimation unit is used for estimating the focal plane position of the down-sampled holograms at different positions by adopting an evaluation method of a characteristic value self-focusing algorithm;
the reconstruction unit is used for reconstructing the down-sampling hologram at the focusing plane position based on the TwinT algorithm to obtain the reconstruction result.
10. The under-sampled hologram compression holographic multi-scale self-focusing reconstruction system according to claim 9, wherein the focal plane position estimation unit is configured to estimate the focal plane position by using an evaluation method of a feature value self-focusing algorithm for the under-sampled holograms at different positions in combination with a multi-scale search strategy, and specifically includes:
1) setting a search start point of a first layer scale
Figure FDA0003160796970000031
Searching for an endpoint
Figure FDA0003160796970000032
And search spacing
Figure FDA0003160796970000033
Obtaining a first layer position sequence, wherein the total interval number of the first layer is N;
2) reconstructing the down-sampling hologram at each position in the current layer position sequence based on a CS algorithm to obtain amplitude images at different positions;
3) the evaluation method of the characteristic value self-focusing algorithm is adopted to estimate and obtain the search result of the current layer scale according to the amplitude images at different positions
Figure FDA0003160796970000034
4) Taking the search result of the current layer scale as the search result of the previous layer scale, and taking the search result of the previous layer scale as the search result of the previous layer scale
Figure FDA0003160796970000041
As a central point, based on a relational expression
Figure FDA0003160796970000042
Obtaining a search starting point of a next search space
Figure FDA0003160796970000043
Searching for an endpoint
Figure FDA0003160796970000044
And search spacing
Figure FDA0003160796970000045
And the next layer position sequence, and repeatedly executing the steps 2) -4) until the self-focusing search result of the k-th layer scale is obtained
Figure FDA0003160796970000046
5) Reconstructing a location based on the TwinT algorithm
Figure FDA0003160796970000047
And (4) processing the down-sampling hologram to obtain the reconstruction result.
CN202110790766.4A 2021-07-13 2021-07-13 Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system Active CN113448233B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110790766.4A CN113448233B (en) 2021-07-13 2021-07-13 Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110790766.4A CN113448233B (en) 2021-07-13 2021-07-13 Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system

Publications (2)

Publication Number Publication Date
CN113448233A true CN113448233A (en) 2021-09-28
CN113448233B CN113448233B (en) 2022-08-09

Family

ID=77816014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110790766.4A Active CN113448233B (en) 2021-07-13 2021-07-13 Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system

Country Status (1)

Country Link
CN (1) CN113448233B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113917819A (en) * 2021-10-13 2022-01-11 中国工程物理研究院激光聚变研究中心 Incoherent three-dimensional holographic layered reconstruction method based on Fresnel mask

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080197842A1 (en) * 2007-02-13 2008-08-21 The Board Of Trustees Of The Leland Stanford Junior University K-t sparse: high frame-rate dynamic magnetic resonance imaging exploiting spatio-temporal sparsity
US20140253713A1 (en) * 2011-10-25 2014-09-11 Guangjie Zhai Time-Resolved Single-Photon or Ultra-Weak Light Multi-Dimensional Imaging Spectrum System and Method
EP2816525A1 (en) * 2013-06-18 2014-12-24 Thomson Licensing Method and apparatus for generating a super-resolved image from a single image
CN104407506A (en) * 2014-12-10 2015-03-11 华南师范大学 Compressive sensing theory-based digital holographic imaging device and imaging method
CN111311493A (en) * 2020-02-13 2020-06-19 河北工程大学 Digital holographic image reconstruction method based on deep learning

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080197842A1 (en) * 2007-02-13 2008-08-21 The Board Of Trustees Of The Leland Stanford Junior University K-t sparse: high frame-rate dynamic magnetic resonance imaging exploiting spatio-temporal sparsity
US20140253713A1 (en) * 2011-10-25 2014-09-11 Guangjie Zhai Time-Resolved Single-Photon or Ultra-Weak Light Multi-Dimensional Imaging Spectrum System and Method
EP2816525A1 (en) * 2013-06-18 2014-12-24 Thomson Licensing Method and apparatus for generating a super-resolved image from a single image
CN104407506A (en) * 2014-12-10 2015-03-11 华南师范大学 Compressive sensing theory-based digital holographic imaging device and imaging method
CN111311493A (en) * 2020-02-13 2020-06-19 河北工程大学 Digital holographic image reconstruction method based on deep learning

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113917819A (en) * 2021-10-13 2022-01-11 中国工程物理研究院激光聚变研究中心 Incoherent three-dimensional holographic layered reconstruction method based on Fresnel mask

Also Published As

Publication number Publication date
CN113448233B (en) 2022-08-09

Similar Documents

Publication Publication Date Title
Ongie et al. Deep learning techniques for inverse problems in imaging
US10671939B2 (en) System, method and computer-accessible medium for learning an optimized variational network for medical image reconstruction
Rajagopalan et al. Depth estimation and image restoration using defocused stereo pairs
CN112132959B (en) Digital rock core image processing method and device, computer equipment and storage medium
Zhang et al. Blind image super-resolution via contrastive representation learning
Wang et al. High-fidelity view synthesis for light field imaging with extended pseudo 4DCNN
FR2763721A1 (en) ELECTRONIC IMAGE PROCESSING DEVICE FOR DETECTION OF DIMENSIONAL VARIATIONS
CN113448233B (en) Under-sampling hologram compression holographic multi-scale self-focusing reconstruction method and system
CN116228912B (en) Image compressed sensing reconstruction method based on U-Net multi-scale neural network
CN115830406A (en) Rapid light field depth estimation method based on multiple parallax scales
KR20180114049A (en) Machine vision systems that form one-dimensional digital representations of low information content scenes
CN115578262A (en) Polarization image super-resolution reconstruction method based on AFAN model
CN111598962A (en) Single-pixel imaging method and device based on matrix sketch analysis
CN114396877B (en) Intelligent three-dimensional displacement field and strain field measurement method for mechanical properties of materials
Huang et al. Light-field reconstruction and depth estimation from focal stack images using convolutional neural networks
Dharejo et al. SwinWave-SR: Multi-scale lightweight underwater image super-resolution
CN116385520A (en) Wear surface topography luminosity three-dimensional reconstruction method and system integrating full light source images
CN115761137A (en) High-precision curved surface reconstruction method and device based on mutual fusion of normal vector and point cloud data
CN115546271A (en) Visual analysis method, device, equipment and medium based on depth joint characterization
Zou et al. DiffCR: A fast conditional diffusion framework for cloud removal from optical satellite images
CN111967292B (en) Lightweight SAR image ship detection method
CN114648608A (en) Tunnel three-dimensional model reconstruction method based on MVSNET
CN113240593A (en) Single image reflecting layer removing method based on position sensing
Yoo et al. Bayesian approach for automatic joint parameter estimation in 3D image reconstruction from multi-focus microscope
CN112017113B (en) Image processing method and device, model training method and device, equipment and medium

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