CN109060122A - Two-step phase recovery method, equipment and system based on single intensity measurement - Google Patents
Two-step phase recovery method, equipment and system based on single intensity measurement Download PDFInfo
- Publication number
- CN109060122A CN109060122A CN201810743862.1A CN201810743862A CN109060122A CN 109060122 A CN109060122 A CN 109060122A CN 201810743862 A CN201810743862 A CN 201810743862A CN 109060122 A CN109060122 A CN 109060122A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- phase
- spatial domain
- complex amplitude
- domain
- 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
Links
- 238000011084 recovery Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005259 measurement Methods 0.000 title claims abstract description 19
- 230000009466 transformation Effects 0.000 claims abstract description 23
- 230000003287 optical effect Effects 0.000 claims abstract description 9
- 150000001875 compounds Chemical class 0.000 claims description 73
- 230000009977 dual effect Effects 0.000 claims description 18
- 238000003384 imaging method Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 abstract 1
- 238000012634 optical imaging Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000005070 sampling Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 239000000571 coke Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001027 hydrothermal synthesis Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000002424 x-ray crystallography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Image Processing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a two-step phase recovery method, equipment and system based on single intensity measurement, belonging to the technical field of optical imaging1And an initialization phaseAccording to the initialized spatial amplitude A1And an initialization phaseSynthesizing complex amplitudes g (x, y); performing cyclic iteration of Fourier transformation and inverse Fourier transformation on the complex amplitude g (x, y) to obtain an amplitude estimation value of a complex optical field in a spatial domain on a 2D coded aperture plane; and processing the amplitude information of the spatial domain and the frequency domain amplitude information obtained by measurement based on a dual-intensity phase recovery algorithm, and recovering the phase of the complex optical field in the spatial domain on the 2D coded aperture plane. The invention obviously improves the quality of the reconstruction phase and the reconstruction power.
Description
Technical field
The present invention relates to optical image technology field, in particular to a kind of two step phase recovery sides based on single ionization meter
Method, equipment and system.
Background technique
Phase recovery be from the ionization meter of record combine known priori knowledge restore lose phase information, at
As field plays very important effect, such as: X-ray crystallography, optics, astronomical imaging, microscope and biomedicine.
Gerchberg and two people of Saxton propose Phase Retrieve Algorithm --- the GS algorithm of alternating projection earliest
(Gerchberg-Saxton algorithm).The algorithm is mainly the intensity data of use space domain and Fourier to restore
The phase of light field.Then, Fienup demonstrates GS algorithm with apparent error descent property, and proposes error reduction
(Error Reduction, ER) algorithm and Mixed design export (Hybrid Input-Output, HIO) algorithm.Currently, ER is calculated
Method and HIO algorithm are considered as most efficient method in phase recovery field.Since algorithm above is just for direct transform system, and
Any linear transformation system but be nobody shows any interest in.Therefore, Yang Guozhen and Gu Benyuan proposes width in any linear transformation system
It is theoretical to spend phase-detection, i.e. Y-G algorism (abbreviation Y-G algorithm).2015, Guo et al. optimized iterative algorithm, and proposes two
The improved GS iterative phase recovery algorithms of kind --- space phase disturbs Gerchberg-Saxton algorithm and combination GS mixing is defeated
Enter output algorithm.For two kinds of modified hydrothermal process, the square value of square error value is dropped rapidly to acceptable value, in spatial domain
With the phase that can successfully restore to lose in Fourier, it means that two kinds of algorithms can jump out local minimum and restrain
To global minimum.
The GS algorithm proposed earliest is the ionization meter for two planes, but is unable to measure two in some cases and puts down
The intensity in face.Therefore, related scholar proposes the GS algorithm based on single ionization meter, and restores in conjunction with some priori knowledges
Phase.Since the phase recovery method of single ionization meter has that reconstruction quality is poor, success rate is lower.
Summary of the invention
The purpose of the present invention is to provide a kind of two step phase recovery methods based on single ionization meter, to improve phase weight
Build quality.
In order to achieve the above object, the present invention uses a kind of two step phase recovery methods based on single ionization meter, it is used for
The strength information of the 2D compound light field of imaging sensor capture is handled, the back focal plane of lens is arranged in imaging sensor, thoroughly
The code aperture 2D M is arranged in the front focal plane of mirror, this method comprises:
The amplitude and phase of 2D compound light field are initialized, initialization spatial amplitude A is obtained1And intialization phase
According to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
It is flat to obtain the code aperture 2D for the loop iteration that Fourier transformation, inverse Fourier transform are carried out to complex amplitude g (x, y)
Amplitude estimation value of the compound light field in spatial domain on face；
Based on dual intensity Phase Retrieve Algorithm, the frequency domain amplitude information that amplitude information and measurement to spatial domain obtain is carried out
It handles, phase of the compound light field in spatial domain in the recovery code aperture 2D plane.
Further, the loop iteration that Fourier transformation, inverse Fourier transform are carried out to complex amplitude g (x, y), obtains
The amplitude information of spatial domain, comprising:
S101, Fourier transformation is carried out to complex amplitude g (x, y), obtains the complex amplitude G (ξ, η) of frequency domain；
S102, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the amplitude of frequency domain complex amplitude G (ξ, η), after being synthesized
Frequency domain complex amplitude G ' (ξ, η)；
S103, inverse Fourier transform is carried out to frequency domain complex amplitude G ' (ξ, η), obtains the complex amplitude g ' (x, y) of spatial domain；
S104, in spatial domain, the complex amplitude g ' (x, y) of spatial domain and the code aperture M are done into point multiplication operation, to space
The complex amplitude g ' (x, y) in domain is updated, and obtains g (x, y)；
S105, step S101~S104 is repeated to updated spatial domain complex amplitude g (x, y), until updated
Spatial domain complex amplitude g (x, y) convergence；
S106, in the updated spatial domain complex amplitude g (x, y) convergence when, according to updated spatial domain complex amplitude g (x,
Y) it restrains, obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain.
Further, described to be based on dual intensity Phase Retrieve Algorithm, the frequency that amplitude information and measurement to spatial domain obtain
Domain amplitude information is handled, phase of the compound light field in spatial domain in the recovery code aperture 2D plane, comprising:
S201, using compound light field in the plane of the code aperture 2D spatial domain amplitude estimation value and the initialization phase
PositionSynthesize complex amplitude g1(x,y)；
S202, to complex amplitude g1(x, y) carries out Fourier transformation, obtains the complex amplitude G of frequency domain1(ξ,η)；
S203, the frequency domain amplitude A that record will be used in frequency domain2Instead of frequency domain complex amplitude G1The amplitude of (ξ, η), is synthesized
Frequency domain complex amplitude G afterwards1′(ξ,η)；
S204, to frequency domain complex amplitude G1' (ξ, η) carries out inverse Fourier transform, obtains the complex amplitude g of spatial domain1′(x,y)；
S205, in spatial domain, replace sky in the amplitude estimation value of spatial domain with compound light field in the plane of the code aperture 2D
Between domain complex amplitude g1The amplitude of ' (x, y) obtains the updated spatial domain complex amplitude g of iteration1(x,y)；
S206, spatial domain complex amplitude g updated to iteration1(x, y) repeats step S202~S205, until iteration
Updated spatial domain complex amplitude g1The phase convergence of (x, y) obtains in the plane of the code aperture 2D compound light field in spatial domain
Phase.
Further, the initialization spatial amplitude A1For complete 1 amplitude, the intialization phaseFor uniformly with
Machine phase and phase interval are [0, pi/2].
On the other hand, a kind of two step phase recovery equipment based on single ionization meter are provided, comprising: obtain module, initial
Change module, synthesis module, spatial domain intensity estimation module and dual intensity phase restoring module；
Obtain strength information of the module for imaging sensor capture compound light field；
Initialization module obtains initialization spatial amplitude A for initializing to the amplitude and phase of 2D compound light field1With
Intialization phase
Synthesis module is used for according to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
Spatial domain intensity estimation module is used to carry out Fourier transformation, the circulation of inverse Fourier transform to complex amplitude g (x, y)
Iteration obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain；
Dual intensity phase restoring module is used to be based on dual intensity Phase Retrieve Algorithm, amplitude information and measurement to spatial domain
Obtained frequency domain amplitude information is handled, phase of the compound light field in spatial domain in the recovery code aperture 2D plane.
Further, the spatial domain intensity estimation module is for executing following steps:
S101, Fourier transformation is carried out to complex amplitude g (x, y), obtains the complex amplitude G (ξ, η) of frequency domain；
S102, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the amplitude of frequency domain complex amplitude G (ξ, η), after being synthesized
Frequency domain complex amplitude G ' (ξ, η)；
S103, inverse Fourier transform is carried out to frequency domain complex amplitude G ' (ξ, η), obtains the complex amplitude g ' (x, y) of spatial domain；
S104, in spatial domain, the complex amplitude g ' (x, y) of spatial domain and the code aperture M are done into point multiplication operation, to space
The complex amplitude g ' (x, y) in domain is updated, and obtains g (x, y)；
S105, step S101~S104 is repeated to updated spatial domain complex amplitude g (x, y), until updated
Spatial domain complex amplitude g (x, y) convergence；
S106, in the updated spatial domain complex amplitude g (x, y) convergence when, according to updated spatial domain complex amplitude g (x,
Y) it restrains, obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain.
Further, the dual intensity phase restoring module is for executing following steps:
S201, using compound light field in the plane of the code aperture 2D spatial domain amplitude estimation value and the initialization phase
PositionSynthesize complex amplitude g1(x,y)；
S202, to complex amplitude g1(x, y) carries out Fourier transformation, obtains the complex amplitude G of frequency domain1(ξ,η)；
S203, the frequency domain amplitude A that record will be used in frequency domain2Instead of frequency domain complex amplitude G1The amplitude of (ξ, η), is synthesized
Frequency domain complex amplitude G afterwards1′(ξ,η)；
S204, to frequency domain complex amplitude G1' (ξ, η) carries out inverse Fourier transform, obtains the complex amplitude g of spatial domain1′(x,y)；
S205, in spatial domain, replace sky in the amplitude estimation value of spatial domain with compound light field in the plane of the code aperture 2D
Between domain complex amplitude g1The amplitude of ' (x, y) obtains the updated spatial domain complex amplitude g of iteration1(x,y)；
S206, spatial domain complex amplitude g updated to iteration1(x, y) repeats step S202~S205, until iteration
Updated spatial domain complex amplitude g1The phase convergence of (x, y) obtains in the plane of the code aperture 2D compound light field in spatial domain
Phase.
Further, the initialization spatial amplitude A1For complete 1 amplitude, the intialization phaseFor uniformly with
Machine phase and phase interval are [0, pi/2].
On the other hand, a kind of two step phase recovery systems based on single ionization meter, including 2D compound light field, CCD figure are provided
As sensor, fourier lense, the code aperture 2D M and above-mentioned two steps phase recovery equipment；
The code aperture 2D M and ccd image sensor are separately positioned on the front focal plane of fourier lense, back focal plane, 2D compound light
Field is placed on before the M of code aperture and 2D compound light field longitudinal section, the code aperture the 2D section M, fourier lense mirror surface and ccd image
Sensor plane is located in same optical path, and the output end of ccd image sensor is connect with computer.
Compared with prior art, there are following technical effects by the present invention: the present invention is directed to single ionization meter, that is, frequency domain intensity
For, compound light field recycles the frequency obtained by measurement in the amplitude information of spatial domain in the recovery code aperture 2D plane first
Domain intensity and restore obtained spatial domain intensity, restores in the plane of the code aperture 2D compound light field in the phase of spatial domain.This programme
The spatial domain amplitude and frequency domain amplitude that compound light field in the plane of the code aperture 2D can effectively be obtained, overcome certain in conventional method
A little applications are difficult to obtain the dual intensity measurement of spatial domain and frequency domain.Meanwhile it can significantly improve the quality for rebuilding phase
And reconstruction success rate.
Detailed description of the invention
With reference to the accompanying drawing, specific embodiments of the present invention will be described in detail:
Fig. 1 is the flow diagram of the two step phase recovery methods based on single ionization meter；
Fig. 2 is the flow diagram of TSPR algorithm；
Fig. 3 is the optical setup schematic diagram of 2D compound light field imaging process；
Fig. 4 is the structural schematic diagram of the two step phase recovery equipment based on single ionization meter；
Fig. 5 is the structural schematic diagram of the two step phase recovery systems based on single ionization meter；
Fig. 6 is single phase recovery test effect figure；
Fig. 7 is the variation diagram of amplitude and phase with the number of iterations；
Fig. 8 is variation schematic diagram of the Phase Build Out success rate with sample rate；
Fig. 9 be sample rate be 0.3 when 3 kinds of code aperture schematic diagrames；
Figure 10 is the contrast schematic diagram of reconstruction performance under different coding aperture；
Figure 11 is the contrast schematic diagram of Phase Build Out performance under algorithms of different.
Specific embodiment
In order to further explain feature of the invention, reference should be made to the following detailed description and accompanying drawings of the present invention.Institute
Attached drawing is only for reference and purposes of discussion, is not used to limit protection scope of the present invention.
As shown in Figure 1 to Figure 2, present embodiment discloses a kind of two step phase recovery methods based on single ionization meter,
The strength information of 2D compound light field for capturing to imaging sensor is handled, and the rear coke of lens is arranged in imaging sensor
The code aperture 2D M is arranged in the front focal plane in face, lens, and this method comprises the following steps:
S1, the amplitude and phase of 2D compound light field are initialized, obtains initialization spatial amplitude A1And intialization phase
S2, according to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
S3, the loop iteration that Fourier transformation, inverse Fourier transform are carried out to complex amplitude g (x, y), obtain the code aperture 2D
Amplitude estimation value of the compound light field in spatial domain in plane；
S4, it is based on dual intensity Phase Retrieve Algorithm, to the amplitude information of the spatial domain obtained in previous step and measured
To frequency domain amplitude information handled, restore the code aperture 2D plane on compound light field spatial domain phase.
It should be noted that as shown in figure 3, code aperture and CCD are individually positioned in the front focal plane and back focal plane of lens.It is first
First, it is u that compound light field u (x, y) obtains compound light field after code aperture M filteringM(x, y), then be imaged by fourier lense, most
Intensity is recorded in CCD plane afterwards.Expression are as follows:
The main purpose of this embodiment scheme is the ionization meter I from recordCCDThe middle phase information for restoring compound light field.
The present embodiment uses a kind of two step phase recovery (Two-Step Phase based on single ionization meter
Retrieval, TSPR) algorithm carry out the code aperture 2D plane on compound light field the phase of spatial domain recovery, first with
Single intensity phase recovery (Single Intensity Phase Retrieval, SIPR) algorithm restores the amplitude letter of spatial domain
Breath utilizes dual intensity phase recovery (Two further according to compound light field frequency domain intensity and the spatial domain intensity that measurement obtains
Intensity Phase Retrieval, TIPR) algorithm restore the code aperture 2D plane on compound light field spatial domain phase.
In practical applications, measurement intensity information needs the optical measuring apparatus of accurate complexity.Sometimes, it is desirable to while obtaining two
A strength information, but condition may be not achieved in optical detection apparatus.The present embodiment overcomes certain occasions well and is difficult to obtain
The dual intensity measurement for obtaining spatial domain and frequency domain, effectively improves the quality and reconstruction success rate for rebuilding phase.
Further, the above-mentioned loop iteration that Fourier transformation, inverse Fourier transform are carried out to complex amplitude g (x, y), obtains
Compound light field specifically comprises the following steps: in the amplitude estimation value of spatial domain in the plane of the code aperture 2D
S101, Fourier transformation is carried out to complex amplitude g (x, y), obtains the complex amplitude G (ξ, η) of frequency domain, wherein
S102, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the amplitude of frequency domain complex amplitude G (ξ, η), after being synthesized
Frequency domain complex amplitude G ' (ξ, η), wherein G ' (ξ, η)=A2{ j ψ (ξ, η) }, ψ (ξ, η) are the phase of frequency domain complex amplitude G ' (ξ, η)
Position；
S103, inverse Fourier transform is carried out to frequency domain complex amplitude G ' (ξ, η), obtains the complex amplitude g ' (x, y) of spatial domain,
In,
S104, in spatial domain, the complex amplitude g ' (x, y) of spatial domain and the code aperture M are done into point multiplication operation, to space
The complex amplitude g ' (x, y) in domain is updated, and obtains g (x, y)=g ' (x, y) M；
S105, step S101~S104 is repeated to updated spatial domain complex amplitude g (x, y), until updated
Spatial domain complex amplitude g (x, y) convergence；
S106, in the updated spatial domain complex amplitude g (x, y) convergence when, according to updated spatial domain complex amplitude g (x,
Y) restrain, obtain in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain | gest(x, y) |=| g (x, y) |.
Further, it is based on dual intensity Phase Retrieve Algorithm, the frequency domain vibration that amplitude information and measurement to spatial domain obtain
Width information is handled, phase of the compound light field in spatial domain in the recovery code aperture 2D plane, comprising:
S201, using compound light field in the plane of the code aperture 2D spatial domain amplitude estimation value and the initialization phase
PositionIt synthesizes complex amplitude g1 (x, y)；
S202, to complex amplitude g1(x, y) carries out Fourier transformation, obtains the complex amplitude G of frequency domain1(ξ, η), wherein
S203, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the complex amplitude G of frequency domain1The amplitude of (ξ, η), after obtaining synthesis
Frequency domain complex amplitude G1' (ξ, η), wherein G1' (ξ, η)=A2·{jψ1(ξ, η) }, ψ1(ξ, η) is frequency domain complex amplitude G1' (ξ, η's)
Phase；
S204, to frequency domain complex amplitude G1' (ξ, η) carries out inverse Fourier transform, obtains the complex amplitude g of spatial domain1' (x, y),
Wherein
S205, in spatial domain, replace spatial domain multiple in the amplitude estimation value of spatial domain with compound light field in the plane of the code aperture 2D
Amplitude g1The amplitude of ' (x, y) obtains the updated spatial domain complex amplitude g of iteration1(x, y),
WhereinFor spatial domain complex amplitude g1The phase of ' (x, y)；
S206, spatial domain complex amplitude g updated to iteration1(x, y) repeats step S202~S205, until iteration
Updated spatial domain complex amplitude g1The phase convergence of (x, y) obtains in the plane of the code aperture 2D compound light field in spatial domain
Phase.
Further, spatial amplitude A is initialized1For complete 1 amplitude, the intialization phaseFor uniformly randomization phase
Position and phase interval are [0, pi/2].
As shown in figure 4, present embodiment discloses a kind of, the two step phase recovery equipment based on single ionization meter include: to obtain
Module, initialization module, synthesis module, spatial domain intensity estimation module and dual intensity phase restoring module；
Obtain strength information of the module for imaging sensor capture compound light field；
Initialization module obtains initialization spatial domain amplitude A for initializing to the amplitude and phase of 2D compound light field1
And intialization phase
Synthesis module is used for according to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
Spatial domain intensity estimation module is used to carry out Fourier transformation, the circulation of inverse Fourier transform to complex amplitude g (x, y)
Iteration obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain；
Dual intensity phase restoring module is used to be based on dual intensity Phase Retrieve Algorithm, amplitude information and measurement to spatial domain
Obtained frequency domain amplitude information is handled, phase of the compound light field in spatial domain in the recovery code aperture 2D plane.
Wherein, each module in two step phase recovery equipment in the present embodiment based on single ionization meter is for realizing Dan Qiang
Each step in two step phase recovery methods of measurement is spent, is repeated no more at this.The two steps phase recovery equipment can be calculating
The equipment such as machine.
As shown in figure 5, present embodiment discloses a kind of two step phase recovery systems based on single ionization meter, including 2D are multiple
Light field, ccd image sensor, fourier lense, the code aperture 2D M and the described in any item two steps phases of claim 5~8
Restorer；
The code aperture 2D M and ccd image sensor are separately positioned on the front focal plane of fourier lense, back focal plane, 2D compound light
Field is placed on before the M of the code aperture 2D and 2D compound light field longitudinal section, the code aperture the 2D section M, fourier lense mirror surface and CCD figure
It is located in same optical path as sensor plane, the output end of ccd image sensor and the two steps phase recovery equipment (for example count
Calculation machine) connection.
It should be noted that the present embodiment carries out the experiment of single phase recovery, it is therefore an objective to be tested by single phase recovery
To verify the validity of TSPR method.Selected in experiment two width grayscale images (" Lena " and " Cameraman " for 256 × 256 pictures
Element).This two images takes its amplitude and phase with after code aperture (256 × 256 pixel) dot product respectively respectively.Then by amplitude
With phase measurements complex amplitude (256 × 256 pixel).The complex amplitude obtains new complex amplitude after being then encoded aperture filtering.Most
Phase recovery experiment is carried out afterwards.Wherein, code aperture is the mask pattern of 0/1 position random distribution, and 0 expression light does not pass through, 1 table
Show that light passes through, the probability that { 0,1 } value occurs all is 1/2.Sample rate in this section experiment is 0.4, phase recovery result such as Fig. 6 institute
Show.Fig. 6 (a) is the code aperture of 0/1 random distribution, and sample rate is 0.4.Fig. 6 (b)-(c) is located at after code aperture
The amplitude and phase of compound light field.Fig. 6 (d) is the amplitude in frequency domain record.Fig. 6 (e)-(f) is the vibration restored using SIPR algorithm
Width and phase are used as objective metric by signal-to-noise ratio (Signal-to-Noise Ratio, SNR), and SNR is respectively
343.04dB and 19.85dB.Fig. 6 (g)-(h) is the amplitude restored using TSPR algorithm and phase, and SNR is respectively 343.04dB
And 45.58dB.Compare Fig. 6 (f)-(h), easily find, the SNR that TSPR algorithm rebuilds phase in this programme is significantly improved.
Fig. 7 be using TSPR algorithm carry out single phase recovery test when amplitude and phase with the number of iterations variation diagram.
From in the figure it can be found that in first step SIPR algorithm, in the plane of the code aperture 2D compound light field the amplitude of spatial domain can be with
It recovers well, but phase recovery result is poor.With the increase of the number of iterations, the signal-to-noise ratio (SNR) of phase can only
Reach 19dB or so.However, the SNR for restoring phase is up to 45dB, and phase obtains the extensive of high quality in second step TIPR algorithm
Multiple, the SIPR method that compares restores phase with failing.The TSPR method that this programme proposes is using estimating in SIPR step
Amplitude, and known frequency domain amplitude is combined, recycle TSPR method successfully to recover the phase of loss, i.e. SIPR and TIPR
Combination.
Fig. 8 is when carrying out single using SIPR algorithm to rebuild experiment, and the reconstruction success rate of amplitude and phase is with sample rate
Change curve.Uniformly random sampling configuration is selected in code aperture, and the sample rate of code aperture is adopted as the variable element in experiment
Sample rate progressively increases to 0.7 from 0.1, increases by 0.05 every time.Using SIPR method carry out rebuild experiment when, calculate rebuild amplitude with
The SNR of the original amplitude and SNR for rebuilding phase and original phase.Rebuild experiment independent operating 500 times, to it is each rebuild SNR into
Row statistics.SNR threshold value is that 25dB, even reconstruction amplitude or the SNR of phase are greater than 25dB, then is considered as and rebuilds successfully.It calculates and rebuilds
The reconstruction success rate of amplitude and phase rebuilds effect picture as depicted in figure 8.It can be seen that SIPR method weight from the experimental result in figure
The success rate for building amplitude is very high, and the success rate for rebuilding phase is very low, the reason is that SIPR algorithm is only recorded by frequency domain
Single ionization meter restores phase.Therefore the present embodiment proposes TSPR method to solve the phase recovery method of single ionization meter and rebuild
Problem of poor quality, success rate is lower.
Fig. 9 be sample rate be 0.3 when 3 kinds of code apertures, i.e., uniformly random (Uniform Random, UR) sampling, diameter
To line (Radial Line, RL) sampling and variable density (Variable Density, VD) sampling configuration.It is different in the present embodiment
Code aperture refer to the different of 0/1 position random distribution, be multiplied with compound light field, i.e., can change as in Fig. 3 compound light field point
Cloth.
Influence for test different coding aperture to context of methods and SIPR algorithm, the present embodiment use three kinds of notches
Diameter, using the sample rate of code aperture as the variable element in experiment, sample rate is arranged in 0.01 to 0.7 range.It is right respectively
Each code aperture is rebuild using TSPR method, calculates the SNR for rebuilding phase and original phase.Every group of identical parameters weight
Experiment independent operating 500 times is built, each reconstruction SNR is counted.SNR threshold value is 25dB, and the SNR of even reconstruction phase is big
In 25dB, then it is considered as and rebuilds successfully.Calculate the reconstruction success rate for rebuilding phase, 10 institute of reconstruction effect picture in different coding aperture
Show.It can be seen from the figure that the TSPR algorithm that the present embodiment proposes is better than traditional SIPR method.TSPR method for uniformly with
It is best to rebuild effect, and has extraordinary stability for the code aperture of machine sampling configuration.
Figure 11 is under the code aperture tested in uniformly random sampling configuration, with the increase of sample rate, TSPR method,
The performance comparison of TIPR algorithm and SIPR algorithm.The sample rate of code aperture is 0.1 to 0.7, and step-length is set as 0.05.To each
Reconstruction SNR counted.SNR threshold value is 25dB, and the SNR of even reconstruction phase is greater than 25dB, then is considered as and rebuilds successfully.It calculates
The reconstruction success rate of phase is rebuild, the reconstruction effect for comparing algorithms of different is as shown in Figure 8.It can be seen that pair from Fig. 8 experimental result
When algorithms of different carries out phase recovery under different sample rate situations, effect most preferably TIPR algorithm is rebuild.The reason is that
Information known to TIPR method is most, i.e., the strength information of two planes.The effect of TSPR method is taken second place, and SIPR method is worst.
TSPR method is a difference in that TSPR method has used for reference the thought of TIPR method with SIPR method.Estimate first with SIPR method
The amplitude in spatial domain is counted out, then restores phase using TIPR method.In short, TSPR method is better than SIPR method.
The present embodiment first restores in the plane of the code aperture 2D compound light field in sky for single ionization meter, that is, frequency domain intensity
Between domain amplitude information, recycle the spatial domain obtained by the obtained frequency domain intensity (i.e. square of amplitude) of measurement and recovery strong
Degree restores compound light field in the plane of the code aperture 2D and is remarkably improved the quality for rebuilding phase in the phase of spatial domain and is reconstructed into
Power.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (9)
1. a kind of two step phase recovery methods based on single ionization meter, which is characterized in that it is used to capture imaging sensor
The strength information of 2D compound light field handled, the back focal plane of lens is arranged in imaging sensor, and 2D is arranged in the front focal plane of lens
Code aperture M, this method comprises:
The amplitude and phase of 2D compound light field are initialized, initialization spatial amplitude A is obtained1And intialization phase
According to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
The loop iteration that Fourier transformation, inverse Fourier transform are carried out to complex amplitude g (x, y), obtains in the plane of the code aperture 2D
Amplitude estimation value of the compound light field in spatial domain；
Based on dual intensity Phase Retrieve Algorithm, at the frequency domain amplitude information that amplitude information and measurement to spatial domain obtain
It manages, phase of the compound light field in spatial domain in the recovery code aperture 2D plane.
2. the two step phase recovery methods based on single ionization meter as described in claim 1, which is characterized in that described pair of multiple vibration
Width g (x, y) carries out the loop iteration of Fourier transformation, inverse Fourier transform, obtains in the plane of the code aperture 2D compound light field in sky
Between domain amplitude estimation value, comprising:
S101, Fourier transformation is carried out to complex amplitude g (x, y), obtains the complex amplitude G (ξ, η) of frequency domain；
S102, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the amplitude of frequency domain complex amplitude G (ξ, η), the frequency domain after being synthesized is multiple
Amplitude G ' (ξ, η)；
S103, inverse Fourier transform is carried out to frequency domain complex amplitude G ' (ξ, η), obtains the complex amplitude g ' (x, y) of spatial domain；
S104, in spatial domain, the complex amplitude g ' (x, y) of spatial domain and the code aperture M are done into point multiplication operation, to spatial domain
Complex amplitude g ' (x, y) is updated, and obtains g (x, y)；
S105, step S101~S104 is repeated to updated spatial domain complex amplitude g (x, y), until updated space
Domain complex amplitude g (x, y) convergence；
When S106, spatial domain complex amplitude g (x, y) in the updated restrain, received according to updated spatial domain complex amplitude g (x, y)
It holds back, obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain.
3. the two step phase recovery methods based on single ionization meter as described in claim 1, which is characterized in that described based on double
Intensity Phase Retrieve Algorithm, the frequency domain amplitude information that amplitude information and measurement to spatial domain obtain are handled, and are restored 2D and are compiled
Phase of the compound light field in spatial domain in code aperture plane, comprising:
S201, amplitude estimation value and the intialization phase using compound light field in the plane of the code aperture 2D in spatial domainIt synthesizes complex amplitude g1 (x, y)；
S202, to complex amplitude g1(x, y) carries out Fourier transformation, obtains the complex amplitude G of frequency domain1(ξ,η)；
S203, the frequency domain amplitude A that record will be used in frequency domain2Instead of frequency domain complex amplitude G1The amplitude of (ξ, η), after being synthesized
Frequency domain complex amplitude G1′(x,y)；
S204, to frequency domain complex amplitude G1' (x, y) carries out inverse Fourier transform, obtains the complex amplitude g of spatial domain1′(x,y)；
S205, in spatial domain, replace spatial domain in the amplitude estimation value of spatial domain with compound light field in the plane of the code aperture 2D
Complex amplitude g1The amplitude of ' (x, y) obtains the updated spatial domain complex amplitude g of iteration1(x,y)；
S206, to updated spatial domain complex amplitude g1(x, y) repeats step S202~S205, until iteration is updated
Spatial domain complex amplitude g1The phase convergence of (x, y) obtains in the plane of the code aperture 2D compound light field in the phase of spatial domain.
4. the two step phase recovery methods based on single ionization meter as described in claim 1, which is characterized in that the initialization
Spatial amplitude A1For complete 1 amplitude, the intialization phaseFor uniformly randomization phase and phase interval is [0, pi/2].
5. a kind of two step phase recovery equipment based on single ionization meter characterized by comprising obtain module, initialization mould
Block, synthesis module, spatial domain intensity estimation module and dual intensity phase restoring module；
Obtain strength information of the module for imaging sensor capture compound light field；
Initialization module obtains initialization spatial amplitude A for initializing to the amplitude and phase of 2D compound light field1With it is initial
Change phase
Synthesis module is used for according to initialization spatial amplitude A1And intialization phaseIt synthesizes complex amplitude g (x, y)；
Spatial domain intensity estimation module is used to carry out Fourier transformation to complex amplitude g (x, y), the circulation of inverse Fourier transform changes
In generation, obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain；
Dual intensity phase restoring module is used to be based on dual intensity Phase Retrieve Algorithm, obtains to the amplitude information and measurement of spatial domain
Frequency domain amplitude information handled, restore the code aperture 2D plane on compound light field spatial domain phase.
6. a kind of two step phase recovery equipment based on single ionization meter as claimed in claim 5, which is characterized in that the sky
Between domain intensity estimation module for executing following steps:
S101, Fourier transformation is carried out to complex amplitude g (x, y), obtains the complex amplitude G (ξ, η) of frequency domain；
S102, the frequency domain amplitude A that will be recorded in frequency domain2Instead of the amplitude of frequency domain complex amplitude G (ξ, η), the frequency domain after being synthesized is multiple
Amplitude G ' (ξ, η)；
S103, inverse Fourier transform is carried out to frequency domain complex amplitude G ' (ξ, η), obtains the complex amplitude g ' (x, y) of spatial domain；
S104, in spatial domain, the complex amplitude g ' (x, y) of spatial domain and the code aperture M are done into point multiplication operation, to spatial domain
Complex amplitude g ' (x, y) is updated, and obtains g (x, y)；
S105, step S101~S104 is repeated to updated spatial domain complex amplitude g (x, y), until updated space
Domain complex amplitude g (x, y) convergence；
When S106, spatial domain complex amplitude g (x, y) in the updated restrain, received according to updated spatial domain complex amplitude g (x, y)
It holds back, obtains in the plane of the code aperture 2D compound light field in the amplitude estimation value of spatial domain.
7. a kind of two step phase recovery equipment based on single ionization meter as claimed in claim 5, which is characterized in that described double
Intensity phase restoring module is for executing following steps:
S201, amplitude estimation value and the intialization phase using compound light field in the plane of the code aperture 2D in spatial domainSynthesize complex amplitude g1(x,y)；
S202, to complex amplitude g1(x, y) carries out Fourier transformation, obtains the complex amplitude G of frequency domain1(ξ,η)；
S203, the frequency domain amplitude A that record will be used in frequency domain2Instead of frequency domain complex amplitude G1The amplitude of (ξ, η), after being synthesized
Frequency domain complex amplitude G1′(ξ,η)；
S204, to frequency domain complex amplitude G1' (ξ, η) carries out inverse Fourier transform, obtains the complex amplitude g of spatial domain1′(x,y)；
S205, in spatial domain, replace spatial domain in the amplitude estimation value of spatial domain with compound light field in the plane of the code aperture 2D
Complex amplitude g1The amplitude of ' (x, y) obtains the updated spatial domain complex amplitude g of iteration1(x,y)；
S206, spatial domain complex amplitude g updated to iteration1(x, y) repeats step S202~S205, until iteration updates
Spatial domain complex amplitude g afterwards1The phase convergence of (x, y) obtains in the plane of the code aperture 2D compound light field in the phase of spatial domain
Position.
8. a kind of two step phase recovery equipment based on single ionization meter as claimed in claim 5, which is characterized in that described first
Beginningization spatial amplitude A1For complete 1 amplitude, the intialization phaseFor uniformly randomization phase and phase interval be [0,
π/2]。
9. a kind of two step phase recovery systems based on single ionization meter, which is characterized in that passed including 2D compound light field, ccd image
Sensor, fourier lense, the code aperture 2D M and the described in any item two steps phase recovery equipment of claim 5~8；
The code aperture 2D M and ccd image sensor are separately positioned on the front focal plane of fourier lense, back focal plane, and 2D compound light field is put
It sets before the M of the code aperture 2D and 2D compound light field longitudinal section, the code aperture the 2D section M, fourier lense mirror surface and ccd image passes
Sensor plane is located in same optical path, and the output end of ccd image sensor is connect with two step phase recovery equipment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810743862.1A CN109060122B (en) | 2018-07-05 | 2018-07-05 | Two-step phase recovery method, equipment and system based on single intensity measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810743862.1A CN109060122B (en) | 2018-07-05 | 2018-07-05 | Two-step phase recovery method, equipment and system based on single intensity measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109060122A true CN109060122A (en) | 2018-12-21 |
| CN109060122B CN109060122B (en) | 2021-02-12 |
Family
ID=64819613
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810743862.1A Active CN109060122B (en) | 2018-07-05 | 2018-07-05 | Two-step phase recovery method, equipment and system based on single intensity measurement |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109060122B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109859127A (en) * | 2019-01-17 | 2019-06-07 | 哈尔滨工业大学 | Object phase recovery technology based on code aperture |
| CN110047048A (en) * | 2019-04-17 | 2019-07-23 | 清华大学深圳研究生院 | It is a kind of to select excellent phase recovery innovatory algorithm based on MSE |
| CN110309482A (en) * | 2019-05-16 | 2019-10-08 | 中国科学院西安光学精密机械研究所 | A kind of rapid convergence, high-precision phase recovery method |
| CN112486003A (en) * | 2020-12-24 | 2021-03-12 | 四川大学 | Phase hologram generation method based on self-adaptive weight feedback GS algorithm |
| CN112946373A (en) * | 2021-02-01 | 2021-06-11 | 北京邮电大学 | Phase-free measuring method and device based on compact range system |
| CN113376448A (en) * | 2021-04-29 | 2021-09-10 | 北京邮电大学 | Method and device for quiet zone phase recovery in compact range test |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060152508A1 (en) * | 2005-01-10 | 2006-07-13 | Fossum Gordon C | System and method for optimized specular highlight generation |
| WO2009003633A1 (en) * | 2007-06-29 | 2009-01-08 | Eni S.P.A. | Process for the conversion of heavy hydrocarbon feedstocks to distillates with the self- production of hydrogen |
| CN103033260A (en) * | 2012-12-14 | 2013-04-10 | 中国科学院国家天文台南京天文光学技术研究所 | Wave surface separation and defocusing based phase retrieval wavefront analyzer and analytical method thereof |
| CN103559698A (en) * | 2013-10-16 | 2014-02-05 | 中国科学院深圳先进技术研究院 | Coaxial phase contrast imaging phase retrieval method and system based on hybrid iteration |
| US20140114615A1 (en) * | 2011-06-13 | 2014-04-24 | Canon Kabushiki Kaisha | Imaging apparatus and program and method for analyzing interference pattern |
| CN104376526A (en) * | 2014-10-24 | 2015-02-25 | 浙江农林大学 | Image encryption method based on vortex beams and phase recovery algorithm |
-
2018
- 2018-07-05 CN CN201810743862.1A patent/CN109060122B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060152508A1 (en) * | 2005-01-10 | 2006-07-13 | Fossum Gordon C | System and method for optimized specular highlight generation |
| WO2009003633A1 (en) * | 2007-06-29 | 2009-01-08 | Eni S.P.A. | Process for the conversion of heavy hydrocarbon feedstocks to distillates with the self- production of hydrogen |
| US20140114615A1 (en) * | 2011-06-13 | 2014-04-24 | Canon Kabushiki Kaisha | Imaging apparatus and program and method for analyzing interference pattern |
| CN103033260A (en) * | 2012-12-14 | 2013-04-10 | 中国科学院国家天文台南京天文光学技术研究所 | Wave surface separation and defocusing based phase retrieval wavefront analyzer and analytical method thereof |
| CN103559698A (en) * | 2013-10-16 | 2014-02-05 | 中国科学院深圳先进技术研究院 | Coaxial phase contrast imaging phase retrieval method and system based on hybrid iteration |
| CN104376526A (en) * | 2014-10-24 | 2015-02-25 | 浙江农林大学 | Image encryption method based on vortex beams and phase recovery algorithm |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109859127A (en) * | 2019-01-17 | 2019-06-07 | 哈尔滨工业大学 | Object phase recovery technology based on code aperture |
| CN110047048A (en) * | 2019-04-17 | 2019-07-23 | 清华大学深圳研究生院 | It is a kind of to select excellent phase recovery innovatory algorithm based on MSE |
| CN110309482A (en) * | 2019-05-16 | 2019-10-08 | 中国科学院西安光学精密机械研究所 | A kind of rapid convergence, high-precision phase recovery method |
| CN112486003A (en) * | 2020-12-24 | 2021-03-12 | 四川大学 | Phase hologram generation method based on self-adaptive weight feedback GS algorithm |
| CN112486003B (en) * | 2020-12-24 | 2021-12-07 | 四川大学 | Phase hologram generation method based on self-adaptive weight feedback GS algorithm |
| CN112946373A (en) * | 2021-02-01 | 2021-06-11 | 北京邮电大学 | Phase-free measuring method and device based on compact range system |
| CN112946373B (en) * | 2021-02-01 | 2024-02-09 | 北京邮电大学 | Compact range system-based non-phase measurement method and device |
| CN113376448A (en) * | 2021-04-29 | 2021-09-10 | 北京邮电大学 | Method and device for quiet zone phase recovery in compact range test |
| CN113376448B (en) * | 2021-04-29 | 2023-02-28 | 北京邮电大学 | Method and device for quiet zone phase recovery in compact range test |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109060122B (en) | 2021-02-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109060122A (en) | Two-step phase recovery method, equipment and system based on single intensity measurement | |
| Tu et al. | Maxim: Multi-axis mlp for image processing | |
| Xu et al. | Ultrasr: Spatial encoding is a missing key for implicit image function-based arbitrary-scale super-resolution | |
| Cui et al. | Focal network for image restoration | |
| Tao et al. | Spectrum-to-kernel translation for accurate blind image super-resolution | |
| Ma et al. | Extensions of compressed imaging: flying sensor, coded mask, and fast decoding | |
| CN116245726A (en) | Compressed sensing polarization super-resolution imaging method based on deep learning framework | |
| Xu et al. | Moiré pattern removal via attentive fractal network | |
| CN116797676A (en) | Method suitable for coded aperture compression spectrum polarization imaging reconstruction | |
| CN107220945A (en) | The restored method of the pole blurred picture of multiple degeneration | |
| Wen et al. | Reconstruction of underwater image by bispectrum | |
| Liu et al. | Efficient unpaired image dehazing with cyclic perceptual-depth supervision | |
| Ahmed et al. | Dynamic MRI using deep manifold self-learning | |
| Cui et al. | PSNet: Towards efficient image restoration with self-attention | |
| Lyu et al. | Learning spatial-angular fusion for compressive light field imaging in a cycle-consistent framework | |
| Zhang et al. | Terahertz image restoration benchmarking dataset | |
| Li et al. | Subspace modeling enabled high-sensitivity X-ray chemical imaging | |
| Song et al. | SODAS-Net: side-information-aided deep adaptive shrinkage network for compressive sensing | |
| Löfdahl et al. | Two dual-wavelength sequences of high-resolution solar photospheric images captured over several hours and restored by use of phase diversity | |
| Zhang et al. | Super-resolution image reconstruction via adaptive sparse representation and joint dictionary training | |
| Samed et al. | An Improved star detection algorithm using a combination of statistical and morphological image processing techniques | |
| Chen et al. | MoCo‐Flow: Neural Motion Consensus Flow for Dynamic Humans in Stationary Monocular Cameras | |
| Zhang et al. | Computational Super-Resolution Imaging With a Sparse Rotational Camera Array | |
| CN109118474B (en) | Image drawing method for multi-view sparse measurement | |
| CN111397733A (en) | Single/multi-frame snapshot type spectral imaging method, system 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 |