CN110187627B - Complex amplitude sensing imaging device and method - Google Patents
Complex amplitude sensing imaging device and method Download PDFInfo
- Publication number
- CN110187627B CN110187627B CN201910571757.9A CN201910571757A CN110187627B CN 110187627 B CN110187627 B CN 110187627B CN 201910571757 A CN201910571757 A CN 201910571757A CN 110187627 B CN110187627 B CN 110187627B
- Authority
- CN
- China
- Prior art keywords
- complex amplitude
- image sensor
- plane
- distribution diagram
- amplitude distribution
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims description 27
- 238000005070 sampling Methods 0.000 claims abstract description 47
- 238000010586 diagram Methods 0.000 claims description 75
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000000149 penetrating effect Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000003287 optical effect Effects 0.000 description 10
- 238000011084 recovery Methods 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000001093 holography Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
- G03H2001/0454—Arrangement for recovering hologram complex amplitude
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Holo Graphy (AREA)
Abstract
The invention discloses a complex amplitude sensing imaging device and a complex amplitude sensing imaging method. The device comprises: the system comprises a lens-free digital camera, a sampling plate, a connecting shell and a processor; the lens-free digital camera comprises an image sensor; the sampling plate is arranged in parallel with the image sensor; the lens-free digital camera is connected with the sampling plate through the connecting shell; a lighting area with the same size and shape as the image sensor is arranged on the sampling plate; a plurality of lighting holes are arranged on the lighting area; the central position of the light collecting area and the central position of the image sensor are positioned on the same horizontal straight line. When light penetrates through an object, the light irradiates the sampling plate, after the light penetrates through the lighting hole, the light is diffracted and spread to the image sensor, the image sensor acquires a light intensity distribution graph of the light wave, and then the processor calculates to obtain a complex amplitude distribution graph of the object by adopting the imaging method disclosed by the invention, so that holographic imaging of the object is realized. The complex amplitude sensing imaging device and the complex amplitude sensing imaging method improve the accuracy of object imaging under the condition of no reference light.
Description
Technical Field
The invention relates to the technical field of holographic imaging, in particular to a complex amplitude sensing imaging device and method.
Background
Because the frequency of the optical wave is very high, the existing photoelectric detector (including the image sensor) can only measure the intensity information of the optical wave, but can not measure the phase information of the optical wave, so the existing image sensor can only obtain the intensity distribution information of the optical wave. At present, to obtain the phase distribution information of the optical wave, a reference beam is usually introduced, which is the digital holography technology. The digital holographic technology for wavefront detection or holographic imaging comprises the following two steps: the first step is to record a hologram formed by the mutual interference of object light and reference light by using a digital image sensor (such as a CCD or CMOS device). And the second step is to input the hologram into the computer and to simulate the light diffraction propagation process to obtain the reconstructed image of the object. The digital holography technology can obtain the intensity information and the phase information of the object light at the same time. However, because reference light is needed, the optical path is complicated, and errors are easily generated due to disturbance of the reference light, and the holographic image of the object cannot be accurately obtained, so that the method is greatly limited in practical application.
For example, in document [1] Horisaki R, Ogura Y, Aino M, et al, single-shot phase imaging with a coded image [ J ]. optlet 2014,39(22):6466-9, the transmission function of the sampling plate is a known binary function, the position of the light transmission part thereof is randomly distributed, and after the object light wave irradiates the sampling plate, the light of the transmission part thereof propagates by diffraction to reach the image sensor surface. During imaging, the image sensor records the diffraction light intensity. Then, the complex amplitude of the light-transmitting part on the sampling plate is reconstructed from the diffracted light intensity recorded by the image sensor by utilizing a phase recovery algorithm. And finally, reconstructing the complex amplitude in the object plane by using the complex amplitude of the light transmission part on the sampling plate by using a compressed sensing algorithm, namely realizing the holographic imaging of the object. The method only is suitable for imaging sparse objects due to the adoption of a compressive sensing technology, and when the phases of the objects are continuously distributed and the fluctuation quantity of the phases is more than 2 pi, an accurate imaging result cannot be obtained.
In documents [2] Cheng ZJ, Wang BY, Xie YY, et al, phase regenerative and differential imaging based on Babinet's principle and complementary random sampling [ J ]. Opt express.2015,23(22):28874-82, a spatial light modulator is used as a sampling plate, and the position of the light-transmitting part thereof may be controlled BY a computer. During imaging, all pixels on the spatial light modulator are randomly divided into four parts. The computer controls one quarter of randomly distributed pixels in the spatial light modulator to be transparent, the rest pixels are opaque, the image sensor records the next diffraction light intensity, and the complex amplitude of the transparent pixels on the sampling plate is reconstructed from the diffraction light intensity recorded by the image sensor by utilizing a phase recovery algorithm. Then, the other quarter of pixels of the spatial light modulator are made to transmit light, the image sensor records a second piece of diffraction light intensity, and the complex amplitude of the second part of light-transmitting pixels on the sampling plate is reconstructed from the second piece of diffraction light intensity recorded by the image sensor by utilizing a phase recovery algorithm. Then, the other quarter of pixels of the spatial light modulator are made to transmit light, the image sensor records a third diffracted light intensity, and the complex amplitude of the third part of light-transmitting pixels on the sampling plate is reconstructed from the third diffracted light intensity recorded by the image sensor by utilizing a phase recovery algorithm. Then, the remaining quarter of pixels of the spatial light modulator are made to transmit light, the image sensor records a fourth diffracted light intensity, and the complex amplitude of the fourth part of pixels of the spatial light modulator on the sampling plate is reconstructed from the fourth diffracted light intensity recorded by the image sensor by utilizing a phase recovery algorithm. Therefore, the complex amplitude of all pixel points in the sampling plate (spatial light modulator) can be reconstructed from the recorded four diffracted light intensities, and then the complex amplitude is reversely propagated into the object plane, so that the complex amplitude in the object plane can be obtained, and the holographic imaging of the object is realized. However, this method uses the spatial light modulator as a sampling plate, which is expensive in equipment cost and requires multiple recordings, and thus accurate imaging of a dynamic object cannot be achieved.
The imaging method disclosed in the document [3] Wang BY, Han L, Yang Y, et al, wave front sensing based on an optical light modulator and an interferometric combining method [ J ]. Opt Lett.2017,42(3):603-6 is similar to the method disclosed in the document [2], except that the number of pixels through which the spatial light modulator is controlled BY a computer is gradually increased. During imaging, when the diffraction light intensity is recorded for the first time, the spatial light modulator is provided with about one-fourth light-transmitting pixels; when the diffraction light intensity is recorded for the second time, the spatial light modulator is provided with about two quarters of light-transmitting pixels; when the diffraction light intensity is recorded for the third time, the spatial light modulator is provided with about three quarters of light-transmitting pixels; when the diffracted intensity is recorded for the fourth time, all pixels of the spatial light modulator are fully transparent. When the image is reconstructed, firstly, the complex amplitude of the first light-transmitting part pixel is reconstructed by the first diffracted light intensity by utilizing a phase recovery algorithm; then, reconstructing the complex amplitude of the second-time light-transmitting two-quarter light-transmitting part pixel by using a phase recovery algorithm according to the diffraction light intensity recorded for the second time and the complex amplitude of the first-time light-transmitting part pixel; then, reconstructing the complex amplitude of the three-quarter light-transmitting part pixel of the third light transmission by utilizing a phase recovery algorithm according to the diffraction light intensity recorded for the third time and the complex amplitude of the second light-transmitting part pixel; and finally, reconstructing the complex amplitude of all the spatial light modulation pixels by the fourth recorded diffraction light intensity and the complex amplitude of the third light transmission part pixels, and reversely transmitting the complex amplitude to the object plane to obtain the complex amplitude in the object plane so as to realize holographic imaging of the object. In this method, a spatial light modulator is used as a sampling plate, which is expensive in equipment cost, and a plurality of recordings are required, and accurate imaging of a dynamic object cannot be realized.
Disclosure of Invention
The invention aims to provide a complex amplitude sensing imaging device and a complex amplitude sensing imaging method, which can improve the accuracy of object imaging under the condition of no reference light.
In order to achieve the purpose, the invention provides the following scheme:
a complex amplitude sensing imaging device, comprising: the system comprises a lens-free digital camera, a sampling plate, a connecting shell and a processor;
the lens-free digital camera comprises an image sensor; the sampling plate is arranged in parallel with the image sensor; the lens-free digital camera and the sampling plate are connected through the connecting shell;
a lighting area with the same size and shape as the image sensor is arranged on the sampling plate; a plurality of lighting holes are formed in the lighting area;
the center position of the light collecting area and the center position of the image sensor are positioned on the same horizontal straight line;
the image sensor is used for acquiring a light intensity distribution diagram of diffracted light penetrating through the light collecting area;
the processor is used for converting the light intensity distribution diagram obtained by the image sensor into a complex amplitude distribution diagram in an image sensor plane and correcting the complex amplitude distribution diagram in the image sensor plane; the processor is further configured to convert the corrected complex amplitude profile in the image sensor plane to a complex amplitude profile of the sample object.
Optionally, the size of the lighting hole is the same as the size of a pixel on the image sensor;
the shape of the lighting hole is the same as that of the pixel on the image sensor.
Optionally, the distance between the lighting holes is equal to the size of a pixel of the image sensor.
Optionally, the connecting shell is a cylindrical tube made of a light-tight material.
A method of complex amplitude sensing imaging, comprising:
step A, acquiring a light intensity distribution diagram of diffracted light of a sample object by adopting the complex amplitude sensing imaging device;
b, generating a first complex amplitude distribution graph in the plane of the image sensor from the light intensity distribution graph;
step C, converting the first complex amplitude distribution diagram in the plane of the image sensor into a first complex amplitude distribution diagram in the plane of a lighting area;
step D, setting the complex amplitude value of the opaque part corresponding to the lighting area in the first complex amplitude distribution map in the lighting area plane to be 0, and obtaining a second complex amplitude distribution map in the lighting area plane;
step E, converting the second complex amplitude distribution diagram in the plane of the lighting area into a second complex amplitude distribution diagram in the plane of the image sensor;
step F, judging whether the difference value between the second complex amplitude distribution diagram in the image sensor plane obtained by the nth calculation and the second complex amplitude distribution diagram in the image sensor plane obtained by the (n-1) th calculation is smaller than a set value or not; wherein n is an integer greater than or equal to 2;
step G, if the difference value is larger than or equal to the set value, extracting the phase distribution of a second complex amplitude distribution diagram in the plane of the image sensor; converting the second complex amplitude distribution map in the image sensor plane into a third complex amplitude distribution map in the image sensor plane according to the phase distribution, replacing the first complex amplitude distribution map in the image sensor plane in the step C with the third complex amplitude distribution map in the image sensor plane, and then repeating the steps C-F;
step H, if the difference value is smaller than the set value, extracting the phase distribution of a second complex amplitude distribution diagram in the plane of the image sensor; converting the second complex amplitude distribution map in the image sensor plane into a third complex amplitude distribution map in the image sensor plane according to the phase distribution;
step I, converting the third complex amplitude distribution diagram in the plane of the image sensor into a third complex amplitude distribution diagram in the plane of a lighting area;
step J, filling a part with a complex amplitude value of 0 in the third complex amplitude distribution diagram in the plane of the lighting area to obtain a fourth complex amplitude distribution diagram in the plane of the lighting area; the fourth complex amplitude distribution diagram in the lighting area plane is the complex amplitude distribution diagram in the lighting area plane;
and K, converting the complex amplitude distribution diagram in the plane of the lighting area into the complex amplitude distribution diagram of the sample object.
Optionally, the converting the second complex amplitude distribution map in the image sensor plane into a third complex amplitude distribution map in the image sensor plane according to the phase distribution includes: by the formula
Calculating to obtain a third complex amplitude distribution diagram in the plane of the image sensor; where exp is an exponential function with a natural constant e as the base, I is a light intensity distribution diagram, j is an imaginary symbol, phi(n)For the n-th extracted phase distribution, n represents the number of iterations and is an integer equal to or greater than 2.
Optionally, an angular spectrum propagation method is used to perform mutual transformation between the complex amplitude distribution diagram in the plane of the image sensor and the complex amplitude distribution diagram in the plane of the lighting area.
Optionally, the filling a portion, with a complex amplitude value of 0, of the third complex amplitude distribution map in the lighting zone plane includes: and carrying out interpolation filling on the opaque part with the complex amplitude value of 0 in the third complex amplitude distribution diagram in the plane of the lighting area by adopting the complex amplitude value corresponding to each lighting hole in the third complex amplitude distribution diagram in the plane of the lighting area.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: in the complex amplitude sensing imaging device provided by the invention, the sampling plate is provided with the lighting area with the same size and shape as the image sensor, and the lighting area is also provided with a plurality of lighting holes, so that light irradiated on an object can be diffracted and transmitted to the image sensor after penetrating through the sampling plate, and the image sensor acquires the light intensity distribution diagram of the light wave. The processor calculates a complex amplitude distribution diagram of the object based on the imaging method disclosed by the invention, so as to realize holographic imaging of the object. In addition, in the whole imaging method process, the object imaging accuracy is further improved by correcting the complex amplitude distribution diagram in the image sensor plane through iterative operation of simulating the repeated propagation of the optical wave complex amplitude between the sampling plate plane and the image sensor plane.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
Fig. 1 is a schematic structural diagram of a complex amplitude sensing imaging apparatus according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a sampling plate in a complex amplitude sensing imaging device according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a complex amplitude sensing imaging apparatus according to an embodiment of the present invention;
fig. 4 is a flowchart of the operation of a complex amplitude sensing imaging method according to an embodiment of the present invention.
Reference numerals: 1-lens-free digital camera, 11-image sensor, 2-sampling plate, 21-light collecting area, 211-light collecting hole and 3-connecting shell.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a complex amplitude sensing imaging device and a complex amplitude sensing imaging method, which can improve the accuracy of object imaging under the condition of no reference light.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Fig. 1 is a schematic structural diagram of a complex amplitude sensing imaging apparatus according to an embodiment of the present invention, and as shown in fig. 1, the complex amplitude sensing imaging apparatus includes: the system comprises a lens-free digital camera 1, a sampling plate 2, a connecting shell 3 and a processor;
the lens-less digital camera 1 includes an image sensor 11; the sampling plate 2 is arranged in parallel with the image sensor 11; the lens-less digital camera 1 and the sampling plate 2 are connected through the connecting shell 3;
a lighting area 21 with the same size and shape as the image sensor 11 is arranged on the sampling plate 2; a plurality of lighting holes 211 are formed in the lighting area 21;
taking a square pixel arrangement as an example, as shown in fig. 2, the size and shape of the lighting hole 211 are the same as those of the pixels on the image sensor 11. The distance between the light collecting holes 211 is equal to the size of the pixel of the image sensor 11. And if the number of pixels of the image sensor 11 is M × N, the number of the light collecting holes 211 on the light collecting area is M × N
When the size of the pixel of the image sensor 11 is Δ, the lighting hole 211 is a square lighting hole having a side of Δ.
In order to further improve the imaging accuracy, the arrangement of the lighting holes 211 in the lighting area 21 is set to be completely consistent with the arrangement of the pixels in the image sensor 11, and each pixel is provided with one lighting hole 211 in the lighting area 21.
The central position of the light collecting area 21 and the central position of the image sensor 11 are located on the same horizontal straight line, so that the light collecting area 21 can be vertically projected onto the image sensor 11, and the projection is just overlapped with the image sensor 11.
As shown in fig. 3, during imaging, after a parallel light passes through an object, the parallel light is irradiated onto the sampling plate 2, wherein a part of the light passes through the light collecting hole 211 and is then diffracted and transmitted onto the image sensor 11, the image sensor 11 records a diffracted light intensity distribution graph (I), the light intensity distribution graph (I) acquired by the image sensor 11 is converted into a complex amplitude distribution graph in the plane of the image sensor by the processor, the complex amplitude distribution graph in the plane of the image sensor 11 is corrected, and finally the corrected complex amplitude distribution graph in the plane of the image sensor 11 is converted into a complex amplitude distribution graph of the sample object.
Wherein, the connecting shell 3 is a cylinder made of opaque material. So that the connection housing 3 can shield surrounding light.
In addition, the complex amplitude sensing imaging method disclosed by the invention is similar to the traditional phase recovery algorithm, iterative operation is realized by simulating the repeated propagation of the light wave complex amplitude between the plane of the sampling plate and the plane of the image sensor, and the complex amplitude of the daylight hole is gradually recovered by using the transmission distribution characteristic of the sampling plate as the constraint condition in the plane of the sampling plate. Fig. 4 is a flowchart of the operation of a complex amplitude sensing imaging method according to an embodiment of the present invention, and as shown in fig. 4, the method includes:
step A, acquiring a light intensity distribution diagram of diffracted light of a sample object by adopting the complex amplitude sensing imaging device;
step B, generating a first complex amplitude distribution graph in the plane of the image sensor 11 from the light intensity distribution graph; the first complex amplitude distribution diagram is
Where j represents an imaginary symbol and the initial value of phi is set to 0 or a random value between 0 and pi.
Step C, adopting an angular spectrum propagation method to distribute a first complex amplitude distribution diagram O in the plane of the image sensor 11r(1)Into a first complex amplitude profile O in the plane of the lighting zone 21s(n)；
Step D, adopting an angular spectrum propagation method to enable the light in the plane of the light collecting area 21First complex amplitude distribution diagram Os(n)The complex amplitude value of the opaque part corresponding to the lighting area is set to be 0, and a second complex amplitude distribution diagram O 'in the plane of the lighting area 21 is obtained's(n)；
Step E, adopting an angular spectrum propagation method to obtain a second complex amplitude distribution diagram O 'in the plane of the light collecting area 21's(n)Conversion to a second complex amplitude profile O 'in the plane of the image sensor 11'r(n)；
Step F, judging a second complex amplitude distribution diagram O 'in the plane of the image sensor 11 obtained by the n-th calculation'r(n)And a second complex amplitude distribution map o 'in the plane of the image sensor 11 obtained by n-1 time calculation'r(n-1)Whether the difference value of (a) is less than a set value; and the set value is a threshold value which is manually set according to actual needs. Wherein n is an integer greater than or equal to 2;
g, if the difference value is larger than or equal to the set value, extracting a second complex amplitude distribution diagram O 'in the plane of the image sensor 11'r(n)The phase distribution of (a); according to the phase distribution, adopting an angular spectrum propagation method to carry out secondary complex amplitude distribution O 'in the plane of the image sensor 11'r(n)To a third complex amplitude profile in the plane of the image sensor 11; a third complex amplitude profile in the plane of the image sensor 11 of
And the first complex amplitude distribution diagram O in the image sensor 11 plane in the step C is obtainedr(1)Is replaced by a third complex amplitude profile O in the plane of the image sensor 11r(n)Then, repeating the steps C-F;
step H, if the difference value is smaller than the set value, extracting a second complex amplitude distribution diagram O 'in the plane of the image sensor 11'r(n)The phase distribution of (a); according to the phase distribution, adopting an angular spectrum propagation method to carry out secondary complex amplitude distribution O 'in the plane of the image sensor 11'r(n)To a third complex amplitude profile in the plane of the image sensor 11; a third complex amplitude profile in the plane of the image sensor 11 of
Where exp is an exponential function with a natural constant e as the base, I is a light intensity distribution diagram, j is an imaginary symbol, phi(n)For the n-th extracted phase distribution, n represents the number of iterations and is an integer equal to or greater than 2.
Step I, converting the third complex amplitude distribution diagram or (n) in the plane of the image sensor 11 into a third complex amplitude distribution diagram O' in the plane of the light collecting area 21 by adopting an angular spectrum propagation methods(n)；
Step J, a third complex amplitude distribution diagram O' in the plane of the light collecting area 21s(n)Filling the part with the medium complex amplitude value of 0 to obtain a fourth complex amplitude distribution diagram O 'in the plane of the lighting area 21's(n)(ii) a A fourth complex amplitude distribution diagram O 'in the plane of the light collecting area 21's(n)Namely a complex amplitude distribution diagram Oo in the plane of the light collecting area 21; wherein the third complex amplitude distribution diagram O' in the plane of the light collecting area 21s(n)Interpolation is used to fill in the part with complex amplitude value 0. During the filling process, the third complex amplitude distribution diagram O' in the plane of the light collecting area 21 is obtaineds(n)The complex amplitude value corresponding to each of the light collection holes 211 is interpolated and filled in the opaque portion having a complex amplitude value of 0 in the third complex amplitude distribution diagram in the plane of the light collection area 21, for example, the complex amplitude values of two light collection holes adjacent to a certain opaque portion are averaged and then filled as the complex amplitude value of the opaque portion.
K, adopting an angular spectrum propagation method to distribute a complex amplitude distribution diagram O in the plane of the light collecting area 21oInto a complex amplitude profile of said sample object.
In the complex amplitude sensing imaging device provided by the invention, the sampling plate is provided with the lighting area with the same size and shape as the image sensor, and the lighting area is also provided with a plurality of lighting holes, so that light irradiated on an object can be diffracted and transmitted to the image sensor after penetrating through the sampling plate, and the image sensor acquires the light intensity distribution diagram of the light wave. The processor calculates a complex amplitude distribution diagram of the object based on the imaging method disclosed by the invention, so as to realize holographic imaging of the object. In addition, in the whole imaging method process, the object imaging accuracy is further improved by correcting the complex amplitude distribution diagram in the image sensor plane through iterative operation of simulating the repeated propagation of the optical wave complex amplitude between the sampling plate plane and the image sensor plane.
In addition, the invention also has the following technical effects:
1. based on the existing phase recovery algorithm, the transmission distribution characteristic of the sampling plate is used as a constraint condition in the plane of the sampling plate, the complex amplitude distribution map of the light wave on each light collecting hole is gradually recovered through iterative operation, the complex amplitude value of the light-tight part of the sampling plate is filled through an interpolation method, the complete light wave complex amplitude distribution map in the plane of the sampling plate is obtained, and holographic imaging of an object to be detected with a complex shape or a phase fluctuation larger than 2 pi can be realized under the conditions of no need of reference light and single exposure.
2. The distance between the centers of the adjacent lighting holes in the device disclosed by the invention is equal to 2 times of the pixel size of the image sensor, namely the sampling interval is equal to 2 times of the pixel size of the image sensor. Thus, the theoretical minimum resolvable distance of the imaging system is equal to 2 times the pixel size of the image sensor.
3. The complex amplitude sensing imaging device and the method provided by the invention can also calculate a complex amplitude distribution diagram O in a lighting area planeoThe phase of the optical wave is detected accurately.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (7)
1. A complex amplitude sensing imaging apparatus, comprising: the system comprises a lens-free digital camera, a sampling plate, a connecting shell and a processor;
the lens-free digital camera comprises an image sensor; the sampling plate is arranged in parallel with the image sensor; the lens-free digital camera and the sampling plate are connected through the connecting shell;
a lighting area with the same size and shape as the image sensor is arranged on the sampling plate; a plurality of lighting holes are formed in the lighting area;
the size of the lighting hole is the same as that of the pixel on the image sensor; the shape of the lighting hole is the same as that of a pixel on the image sensor;
the center position of the light collecting area and the center position of the image sensor are positioned on the same horizontal straight line;
the image sensor is used for acquiring a light intensity distribution diagram of diffracted light penetrating through the light collecting area;
the processor is used for converting the light intensity distribution diagram obtained by the image sensor into a complex amplitude distribution diagram in an image sensor plane and correcting the complex amplitude distribution diagram in the image sensor plane; the processor is further configured to convert the corrected complex amplitude profile in the image sensor plane to a complex amplitude profile of the sample object.
2. A complex amplitude sensing imaging device according to claim 1, wherein the distance between the lighting holes is equal to the size of the pixels of the image sensor.
3. A complex amplitude sensing imaging device according to claim 1, wherein said connecting housing is a cylindrical tube made of a light-impermeable material.
4. A complex amplitude sensing imaging method, characterized by being applied to the complex amplitude sensing imaging apparatus according to any one of claims 1 to 3; the method comprises the following steps:
step A, acquiring a light intensity distribution diagram of diffracted light of a sample object by adopting the complex amplitude sensing imaging device;
b, generating a first complex amplitude distribution graph in the plane of the image sensor from the light intensity distribution graph;
step C, converting the first complex amplitude distribution diagram in the plane of the image sensor into a first complex amplitude distribution diagram in the plane of a lighting area;
step D, setting the complex amplitude value of the opaque part corresponding to the lighting area in the first complex amplitude distribution map in the lighting area plane to be 0, and obtaining a second complex amplitude distribution map in the lighting area plane;
step E, converting the second complex amplitude distribution diagram in the plane of the lighting area into a second complex amplitude distribution diagram in the plane of the image sensor;
step F, judging whether the difference value between the second complex amplitude distribution diagram in the image sensor plane obtained by the nth calculation and the second complex amplitude distribution diagram in the image sensor plane obtained by the (n-1) th calculation is smaller than a set value or not; wherein n is an integer greater than or equal to 2;
step G, if the difference value is larger than or equal to the set value, extracting the phase distribution of a second complex amplitude distribution diagram in the plane of the image sensor; converting the second complex amplitude distribution map in the image sensor plane into a third complex amplitude distribution map in the image sensor plane according to the phase distribution, replacing the first complex amplitude distribution map in the image sensor plane in the step C with the third complex amplitude distribution map in the image sensor plane, and then repeating the steps C-F;
step H, if the difference value is smaller than the set value, extracting the phase distribution of a second complex amplitude distribution diagram in the plane of the image sensor; converting the second complex amplitude distribution map in the image sensor plane into a third complex amplitude distribution map in the image sensor plane according to the phase distribution;
step I, converting the third complex amplitude distribution diagram in the plane of the image sensor into a third complex amplitude distribution diagram in the plane of a lighting area;
step J, filling a part with a complex amplitude value of 0 in the third complex amplitude distribution diagram in the plane of the lighting area to obtain a fourth complex amplitude distribution diagram in the plane of the lighting area; the fourth complex amplitude distribution diagram in the lighting area plane is the complex amplitude distribution diagram in the lighting area plane;
and K, converting the complex amplitude distribution diagram in the plane of the lighting area into the complex amplitude distribution diagram of the sample object.
5. A method of complex amplitude sensing imaging according to claim 4, wherein said converting the second complex amplitude profile in the image sensor plane to a third complex amplitude profile in the image sensor plane according to the phase distribution comprises: by the formula
Calculating to obtain a third complex amplitude distribution diagram in the plane of the image sensor; where exp is an exponential function with a natural constant e as the base, I is a light intensity distribution diagram, j is an imaginary symbol, phi(n)For the n-th extracted phase distribution, n represents the number of iterations and is an integer equal to or greater than 2.
6. The method of claim 4, wherein the inter-conversion between the complex amplitude profile in the plane of the image sensor and the complex amplitude profile in the plane of the illuminated area is performed using angular spectrum propagation.
7. The method according to claim 4, wherein said filling the portion of the third complex amplitude profile in the plane of the illuminated area with complex amplitude values of 0 comprises: and carrying out interpolation filling on the opaque part with the complex amplitude value of 0 in the third complex amplitude distribution diagram in the plane of the lighting area by adopting the complex amplitude value corresponding to each lighting hole in the third complex amplitude distribution diagram in the plane of the lighting area.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910571757.9A CN110187627B (en) | 2019-06-28 | 2019-06-28 | Complex amplitude sensing imaging device and method |
| PCT/CN2020/086411 WO2020259042A1 (en) | 2019-06-28 | 2020-04-23 | Complex amplitude sensing imaging device and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910571757.9A CN110187627B (en) | 2019-06-28 | 2019-06-28 | Complex amplitude sensing imaging device and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110187627A CN110187627A (en) | 2019-08-30 |
| CN110187627B true CN110187627B (en) | 2021-04-02 |
Family
ID=67724067
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910571757.9A Expired - Fee Related CN110187627B (en) | 2019-06-28 | 2019-06-28 | Complex amplitude sensing imaging device and method |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110187627B (en) |
| WO (1) | WO2020259042A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110187627B (en) * | 2019-06-28 | 2021-04-02 | 赣南师范大学 | Complex amplitude sensing imaging device and method |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3670046B2 (en) * | 1995-02-10 | 2005-07-13 | 大日本印刷株式会社 | Diffraction grating recording medium and method and apparatus for producing the same |
| CN101046882A (en) * | 2006-03-27 | 2007-10-03 | 谷京陆 | Mixed resolution image making technology capable of direct output |
| WO2010072065A1 (en) * | 2008-12-25 | 2010-07-01 | 深圳市泛彩溢实业有限公司 | Hologram three-dimensional image information collecting device and method, reproduction device and method |
| CN101726366B (en) * | 2009-12-02 | 2011-11-23 | 山东师范大学 | Wavefront measurement method and device based on multiple-pinhole plate |
| WO2012030357A1 (en) * | 2010-09-03 | 2012-03-08 | Arges Imaging, Inc. | Three-dimensional imaging system |
| US10359513B2 (en) * | 2017-05-03 | 2019-07-23 | Elwha Llc | Dynamic-metamaterial coded-aperture imaging |
| CN106950811B (en) * | 2017-05-08 | 2019-01-25 | 赣南师范大学 | A kind of digital composite holographic imaging method and device |
| CN107861360B (en) * | 2017-12-20 | 2020-01-10 | 清华大学 | Single-exposure lens-free imaging system and method based on multi-angle illumination multiplexing |
| CN108645336B (en) * | 2018-05-11 | 2020-06-09 | 赣南师范大学 | Reference-light-free digital holographic camera and calibration method |
| CN109859127A (en) * | 2019-01-17 | 2019-06-07 | 哈尔滨工业大学 | Object phase recovery technology based on code aperture |
| CN110187627B (en) * | 2019-06-28 | 2021-04-02 | 赣南师范大学 | Complex amplitude sensing imaging device and method |
-
2019
- 2019-06-28 CN CN201910571757.9A patent/CN110187627B/en not_active Expired - Fee Related
-
2020
- 2020-04-23 WO PCT/CN2020/086411 patent/WO2020259042A1/en active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| CN110187627A (en) | 2019-08-30 |
| WO2020259042A1 (en) | 2020-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6549240B2 (en) | Lensless imaging system using an image sensor with one or more attenuation layers | |
| Salvi et al. | A state of the art in structured light patterns for surface profilometry | |
| CN107306333B (en) | High-speed single-pixel imaging method | |
| CN104407506A (en) | Compressive sensing theory-based digital holographic imaging device and imaging method | |
| US3669673A (en) | Recording of a continuous tone focused image on a diffraction grating | |
| CN103913229B (en) | Coding template multi-target active imaging spectrum system and method | |
| KR102404416B1 (en) | Apparatus and method of image-quality measurement for hologram images | |
| CN105404128B (en) | Multiframe phase-shifted digital holographic method and device | |
| CN113114882B (en) | Fourier single-pixel imaging method with high sampling efficiency | |
| CN110187627B (en) | Complex amplitude sensing imaging device and method | |
| JP2011226871A (en) | Shape measurement method, shape measurement device, distortion measurement method, and distortion measurement device | |
| CN104765263A (en) | Compressing, coding and transmitting method for electronic holographic three-dimensional information | |
| US20220113674A1 (en) | Differential holography | |
| CN113504202B (en) | Coherent modulation imaging method based on axial translation binary amplitude mask | |
| CN204360096U (en) | Based on the digital hologram imaging device of compressed sensing theory | |
| CN102095503B (en) | Wavefront detection and reconstruction method based on differential sensor | |
| CN104090476A (en) | Three-dimensional scene information acquisition method for holographic display | |
| CN108957999B (en) | Phase shift holographic device based on phase type vortex lens and imaging method | |
| CN116337010A (en) | Associated imaging reconstruction recovery method for intelligent optimization of speckle | |
| CN115307577A (en) | Target three-dimensional information measuring method and system | |
| EP4160319A1 (en) | Electromagnetic wave phase/amplitude generation device, electromagnetic wave phase/amplitude generation method, and electromagnetic wave phase/amplitude generation program | |
| Vetal et al. | Reconstruction of apertured Fourier Transform Hologram using compressed sensing | |
| CN112702486A (en) | High-speed coherent imaging camera based on electrically adjustable lens and phase recovery method | |
| CN114764220B (en) | Method for improving speckle autocorrelation reconstruction effect based on off-axis digital holography | |
| CN103969659B (en) | Celestial body spectrum image imaging method based on compressed sensing |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210402 Termination date: 20210628 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |