CN109472842A - A kind of phase recovery imaging reconstruction method of no lens imaging - Google Patents
A kind of phase recovery imaging reconstruction method of no lens imaging Download PDFInfo
- Publication number
- CN109472842A CN109472842A CN201811541689.3A CN201811541689A CN109472842A CN 109472842 A CN109472842 A CN 109472842A CN 201811541689 A CN201811541689 A CN 201811541689A CN 109472842 A CN109472842 A CN 109472842A
- Authority
- CN
- China
- Prior art keywords
- algorithm
- imaging
- phase
- power spectrum
- ptychography
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The present invention proposes a kind of phase recovery imaging reconstruction method of no lens imaging, fast power spectrogram is captured using coherent diffraction imaging experimental provision, imaging device only needs a collimator and extender laser, one fourier lense and a ccd detector, the power spectrum chart detected is carried out to the reconstruct of three steps, the first step calculates an initial phase by Mixed design output algorithm, then soft aperture scanning is carried out, rather than it is scanned in objective plane and spatial spectrum plane, second step combination Ptychography algorithm enhances Fourier constraint, with phase Iterative restoration algorithm, target image is reconstructed in realization, third step, in order to further decrease normalization mean square deviation root, the algorithm of error reduction is carried out herein, reconstruct target image.Imaging device of the present invention is simple, and algorithm can quickly reconstruct accurate target image from power spectrum chart, realizes the function of image reproducing.
Description
Technical field
The present invention relates to coherent diffraction imaging field more particularly to a kind of phase recovery image reproducing sides of no lens imaging
Method.
Background technique
In a practical situation, image measured directly is usually the intensity distribution of wave field, and the phase distribution of wave field is past
Toward being difficult directly to measure, even it is not possible that, therefore, restoring phase information from intensity measurement data is always Recent study
Hot spot, Phase Retrieve Algorithm be it is this by intensity measurement data export radiation source Amplitude & Phase method.
Currently, there are two ways to common phase recovery, one is (wave specifically can be used using phase difference sensor
Front sensor) direct detection optical system phase difference, another kind be by utilize Gerchberg-Saxton (GS) algorithm, with
Optimal theoretical image is target, is iterated using Fourier transformation and inverse Fourier transform and corrects error, after final convergence
To the phase difference of optical system.
But the detection of optical system phase difference is carried out using Wavefront sensor, it needs to pacify in the specific position of optical system
Wavefront sensor is filled, while also needing to be arranged on imaging surface a point light source as guidance, is brought to the design of optical system
Biggish limitation；And iterative Fourier transform (Iterative Fourier Transform, abbreviation IFT) is considered current
Best Phase Retrieve Algorithm, regrettably IFT algorithm does not ensure that iterative process always converges to correct solution, sometimes even
Some local minimum, which can be stuck in, cannot nearby escape from, to can not restore target image.
Summary of the invention
In view of the problems of the existing technology, the present invention proposes a kind of phase recovery image reproducing side of no lens imaging
Method solves the problems, such as that optical system builds the diffraction limit encountered in the prior art, also solves IFT algorithm in an iterative process
It is stuck in the problem of some local minimum cannot restore target image, realizes that experimental facilities is simple to operation, iterative process
The target that operation time is few, restores clearly effect.
To achieve this purpose, the present invention adopts the following technical scheme:
The embodiment of the invention provides a kind of phase recovery imaging reconstruction methods of no lens imaging, comprising:
1. obtaining the power spectrum chart that target is obtained by fourier lense, specifically apply are as follows: use collimator and extender laser illumination
When target, target placement position is not necessarily to particular focal length f, can be with any position；Obtaining power spectrum chart only need to will be after laser irradiation
Target places ccd detector by fourier lense at focal length f, is not necessarily to imaging len.
2. introducing the solution of Ptychography restructing algorithm in iterative algorithm according to acquired power spectrum chart falls into office
Portion's minimum problem, particular by a kind of soft scanning of the scan aperture of non-physical, by carrying out redundancy to Fourier spectrum
Scanning solves to restrain stagnation problem in coherent diffraction imaging method.
3. combining IFT algorithm and Ptychography algorithm to carry out the iterative phase Restoration model of three steps, according to institute
State the reproduction of model realization target image.
Step 1 obtains initial phase with the Mixed design output algorithm in iterative Fourier transform algorithm；
Step 2 uses Ptychography algorithm and the result of soft scanning step one；
Step 3 reduces normalized mean squared error with the method that the error in iterative Fourier transform algorithm reduces.
The embodiment of the invention provides a kind of phase recovery imaging reconstruction method of no lens imaging, this method utilizes relevant
Diffraction imaging experimental provision captures fast power spectrogram, by the power spectrum chart detected carry out three steps reconstruct, first
Step calculates an initial phase by Mixed design output algorithm, then carries out soft aperture scanning, rather than in objective plane and
Spatial spectrum plane is scanned, in conjunction with the enhancing Fourier constraint of Ptychography algorithm, with phase Iterative restoration algorithm,
Target image is reconstructed in realization, in order to further decrease normalization mean square deviation root, carries out the algorithm of error reduction, weight herein
Structure goes out target image.Imaging device of the present invention is simple, and algorithm can quickly reconstruct accurate target image from power spectrum chart, with
Original image is minimum compared to mean square error, realizes the function of image reproducing.
Detailed description of the invention
Fig. 1 is the schematic diagram of CDI imaging system
Fig. 2 is the algorithm flow chart for reconstructing algorithm steps one and step 3
Fig. 3 is the algorithm flow chart for reconstructing algorithm steps two
Fig. 4 is that schematic diagram is scanned in aperture
Fig. 5 is the power spectrum chart of target ' F '
Fig. 6 is the method for the present invention image reproducing figure
Specific embodiment
Below with reference to the drawings and specific embodiments, the present invention is described in detail.
The present invention is to obtain power spectrum chart based on the simple optical device of CDI, and experimental principle figure is as shown in Figure 1, a branch of swash
Light is by collimation and expands, and irradiates target, by fourier lense, receives the spectrogram of target with camera at focal length f.In reality
The middle laser for using and having collimator and extender is tested, and feux rouges (λ=632.8nm) is selected to irradiate target object, irradiates object
Light source is directional light, and object can be placed in any one position between fourier lense and light source, after laser irradiation target
By fourier lense, target is obtained in the spectrogram of Fourier, is then acquired with the CCD of no imaging lens, in CCD and Fu
The distance of leaf lens is the focal length of fourier lense.
The spectrogram of target is obtained for experiment, reconstructing method of the invention includes three steps, about algorithm flow chart
As shown in Figures 2 and 3, the meaning of parameters in two width figures is defined as follows: F indicates Fast Fourier Transform (FFT)；Arg indicates amplitude；Pk
Indicate target ΘkNew phase；Indicate the intensity of testee；Rpl indicates target ΘkIntensity quiltReplacement；F-1Table
Show inverse Fourier transform；Cons is indicated from O'kObtain OkPositive real number space constraint；Prc indicates the entire mistake of some step
Journey；Scn indicates scanning frequency spectrum；Upd expression continues to update, and step 1 and step 3 have similar iterative process, respectively k1 times
1. 3. iteration and k3 iteration, output result indicate that step 2 is the process of a Ptychography algorithm, the calculation with Y with Y
Method can be used to reconstruct the spectrum I of high qualityh。
Test the power spectrum chart I of measurementmIt is the unique target data of reconstruct, step 1 and step 3 have similar iteration mistake
Journey, as shown in Figure 2.Step 2 flow chart is as shown in Figure 3.A kind of no lens imaging phase recovery imaging reconstruction method, in step
K1 iteration is run to Mixed design-output (Hybrid input-output, HIO) algorithm in one, wherein defining the pact in object domain
Beam γ must be positive real number, the target initial phase P being calculated in step 1k1Step 2 is calculated closer to true result
There is significant contribution.Step 2 is to combine Ptychography restructing algorithm, wherein ImBe size be defined aperture scanning after
Power spectrum chart, bring into the algorithm flow chart of Fig. 3, obtain the reconstruct power spectrum chart of high quality by k2 iteration, based on miss
Difference reduces algorithm, and the output of step 2 calculates the known conditions for replacing the intensity of measurement new as step 3, by k3
Secondary iteration exports reconstruction result figure.
The present invention has passed through experimental demonstration, and the power spectrum chart of 480*480 is obtained by CDI imaging system, and original object is
' F ' in USAF 1951target, the power spectrum chart for obtaining utilize feedback parameter β after step 1 is run k1=20 times
=0.7 quickly calculates accurate initial phase, then, as shown in figure 4, selecting pixel size for 30 scan aperture, carries out
256 scanning (every row sweeps 16 times, totally 16 row), the overlapping area between adjacent apertures is the 2/3 of aperture area, in order to further
Normalized mean squared error is reduced, this method has also carried out k3=10 error and reduced algorithm, during iterative reconstruction, calculates mesh
Mark the mean square deviation E in domainkAs the standard of assessment reconstructed results quality, as shown in formula (1).
The power spectrum chart measure to target ' F ' is as shown in figure 5, target reproduction figure is as shown in Figure 6 by the method for the invention.
Claims (4)
1. a kind of phase recovery imaging reconstruction method of no lens imaging characterized by comprising
It obtains and the power spectrum chart of no lens imaging is obtained by fourier lense；
The solution of Ptychography restructing algorithm is introduced in phase recovery iterative algorithm according to acquired power spectrum chart to fall into
Local minimum problem；
The iterative phase Restoration model that three steps are carried out in conjunction with IFT algorithm and Ptychography algorithm, according to the model
Realize that target image reproduces.
2. obtaining the function of no lens imaging by fourier lense the method according to claim 1, wherein obtaining
Rate spectrogram, comprising:
When with collimator and extender laser illumination target, target placement position is not necessarily to particular focal length f, can be with any position；
Obtain power spectrum chart only need to place ccd detector by the target after laser irradiation by fourier lense at focal length f,
Without imaging len.
3. the method according to claim 1, wherein being calculated according to acquired power spectrum chart in phase recovery iteration
Ptychography restructing algorithm is introduced in method to solve the problems, such as to fall into local minimum, comprising:
Relevant spread out is solved by carrying out redundancy scanning to Fourier spectrum by a kind of soft scanning of the scan aperture of non-physical
It penetrates and restrains stagnation problem in imaging method.
4. the method according to claim 1, wherein carrying out three in conjunction with IFT algorithm and Ptychography algorithm
The iterative phase Restoration model of a step is reproduced according to the model realization target image, comprising:
Step 1 obtains initial phase with the Mixed design output algorithm in iterative Fourier transform algorithm；
Step 2 uses Ptychography algorithm and the result of soft scanning step one；
Step 3 reduces normalized mean squared error with the method that the error in iterative Fourier transform algorithm reduces.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541689.3A CN109472842A (en) | 2018-12-17 | 2018-12-17 | A kind of phase recovery imaging reconstruction method of no lens imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541689.3A CN109472842A (en) | 2018-12-17 | 2018-12-17 | A kind of phase recovery imaging reconstruction method of no lens imaging |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109472842A true CN109472842A (en) | 2019-03-15 |
Family
ID=65676076
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811541689.3A Pending CN109472842A (en) | 2018-12-17 | 2018-12-17 | A kind of phase recovery imaging reconstruction method of no lens imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109472842A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110047048A (en) * | 2019-04-17 | 2019-07-23 | 清华大学深圳研究生院 | It is a kind of to select excellent phase recovery innovatory algorithm based on MSE |
| CN110231310A (en) * | 2019-06-28 | 2019-09-13 | 长沙理工大学 | A kind of imaging method for the target being hidden in opaque scattering medium |
| CN113030024A (en) * | 2021-03-24 | 2021-06-25 | 浙江大学 | Coherent diffraction imaging device and method based on low dynamic range spectrogram |
| CN113203485A (en) * | 2021-04-27 | 2021-08-03 | 浙江大学 | Device and method for realizing axial phase difference wavefront reconstruction through single exposure |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6097856A (en) * | 1998-07-10 | 2000-08-01 | Welch Allyn, Inc. | Apparatus and method for reducing imaging errors in imaging systems having an extended depth of field |
| JP2012238588A (en) * | 2011-04-25 | 2012-12-06 | Hokkaido Univ | Fourier iterative phase retrieval method |
| CN103942805A (en) * | 2014-05-08 | 2014-07-23 | 长沙理工大学 | Rapid image sparse decomposition method based on partial polyatomic matching pursuit |
| US20150036038A1 (en) * | 2013-07-31 | 2015-02-05 | California Institute Of Technology | Aperture scanning fourier ptychographic imaging |
| CN107942523A (en) * | 2017-12-05 | 2018-04-20 | 安徽大学 | Phase recovery system based on light intensity transmission measurement calculation |
| CN108537862A (en) * | 2018-04-11 | 2018-09-14 | 北京理工大学 | A kind of Fourier's Diffraction scans microscope imaging method of adaptive noise reduction |
| CN108550108A (en) * | 2017-09-28 | 2018-09-18 | 武汉大学 | A kind of Fourier's lamination image method for reconstructing minimized based on phase iteration |
-
2018
- 2018-12-17 CN CN201811541689.3A patent/CN109472842A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6097856A (en) * | 1998-07-10 | 2000-08-01 | Welch Allyn, Inc. | Apparatus and method for reducing imaging errors in imaging systems having an extended depth of field |
| JP2012238588A (en) * | 2011-04-25 | 2012-12-06 | Hokkaido Univ | Fourier iterative phase retrieval method |
| US20150036038A1 (en) * | 2013-07-31 | 2015-02-05 | California Institute Of Technology | Aperture scanning fourier ptychographic imaging |
| CN103942805A (en) * | 2014-05-08 | 2014-07-23 | 长沙理工大学 | Rapid image sparse decomposition method based on partial polyatomic matching pursuit |
| CN108550108A (en) * | 2017-09-28 | 2018-09-18 | 武汉大学 | A kind of Fourier's lamination image method for reconstructing minimized based on phase iteration |
| CN107942523A (en) * | 2017-12-05 | 2018-04-20 | 安徽大学 | Phase recovery system based on light intensity transmission measurement calculation |
| CN108537862A (en) * | 2018-04-11 | 2018-09-14 | 北京理工大学 | A kind of Fourier's Diffraction scans microscope imaging method of adaptive noise reduction |
Non-Patent Citations (2)
| Title |
|---|
| PAN X C: "Ptychographical imaging with partially saturated diffraction patterns[J].Journal of Modern Optics", 《JOURNAL OF MODERN OPTICS》, vol. 62, no. 15, 31 December 2015 (2015-12-31), pages 1270 - 1277 * |
| 王新华等: "压缩编码孔径红外成像超分辨重建", 《吉林大学学报(工学版)》, no. 04, 9 October 2015 (2015-10-09), pages 228 - 234 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110047048A (en) * | 2019-04-17 | 2019-07-23 | 清华大学深圳研究生院 | It is a kind of to select excellent phase recovery innovatory algorithm based on MSE |
| CN110231310A (en) * | 2019-06-28 | 2019-09-13 | 长沙理工大学 | A kind of imaging method for the target being hidden in opaque scattering medium |
| CN113030024A (en) * | 2021-03-24 | 2021-06-25 | 浙江大学 | Coherent diffraction imaging device and method based on low dynamic range spectrogram |
| CN113203485A (en) * | 2021-04-27 | 2021-08-03 | 浙江大学 | Device and method for realizing axial phase difference wavefront reconstruction through single exposure |
| CN113203485B (en) * | 2021-04-27 | 2022-08-05 | 浙江大学 | Device and method for realizing axial phase difference wavefront reconstruction through single exposure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109472842A (en) | A kind of phase recovery imaging reconstruction method of no lens imaging | |
| KR101226210B1 (en) | High resolution imaging | |
| EP2356487B1 (en) | Provision of image data | |
| JP5917507B2 (en) | Calibration method of probe by typography method | |
| Huang et al. | Continuous-wave terahertz multi-plane in-line digital holography | |
| CN106680240B (en) | A kind of double object distance lamination imaging methods of continuous THz wave | |
| CN104161536B (en) | Cone beam CT scatter correction method and device based on complementary gratings | |
| CN105548080B (en) | A kind of continuous THz wave spacescan coherent diffraction imaging system and method | |
| Berujon et al. | X-ray pulse wavefront metrology using speckle tracking | |
| BR112020007609B1 (en) | IMAGE RECONSTRUCTION METHOD AND DEVICE AND MICROSCOPIC IMAGING DEVICE | |
| CN109188881B (en) | Large-field-of-view high-resolution terahertz wave digital holographic imaging method and system | |
| Savransky et al. | Focal plane wavefront sensing and control for ground-based imaging | |
| CN110360924A (en) | A kind of laminated micro imaging method of double wave and system | |
| CN105675151A (en) | Application device based on light intensity transmission equation phase retrieval | |
| CN105973897A (en) | Measuring device and method for geometric size distribution of needle damage loci of KDP crystal | |
| Hefferan et al. | Tests of microstructure reconstruction by forward modeling of high energy X-ray diffraction microscopy data | |
| CN106679586B (en) | Modulation ePIE Phase Retrieve Algorithm is scanned based on entrance pupil | |
| CN110231310A (en) | A kind of imaging method for the target being hidden in opaque scattering medium | |
| Misell et al. | Phase determination from image intensity measurements in bright-field optics | |
| JPWO2017158695A1 (en) | Point spread function measuring device, measuring method, image acquiring device, and image acquiring method | |
| Whitehead et al. | Fresnel diffractive imaging: Experimental study of coherence and curvature | |
| JP5677751B2 (en) | X-ray CT apparatus and image processing program | |
| CA2836269A1 (en) | Improvements in and relating to ophthalmoscopes | |
| WO2019011692A1 (en) | Phase contrast imaging method | |
| CN112525935A (en) | Laminated diffraction imaging device and method |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190315 |