CN112051247B - Lens-free imaging device based on laminated imaging and phase recovery method thereof - Google Patents

Lens-free imaging device based on laminated imaging and phase recovery method thereof Download PDF

Info

Publication number
CN112051247B
CN112051247B CN202010847763.5A CN202010847763A CN112051247B CN 112051247 B CN112051247 B CN 112051247B CN 202010847763 A CN202010847763 A CN 202010847763A CN 112051247 B CN112051247 B CN 112051247B
Authority
CN
China
Prior art keywords
sample
scattering layer
image sensor
image
solution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010847763.5A
Other languages
Chinese (zh)
Other versions
CN112051247A (en
Inventor
颜成钢
朱嘉凯
郑国安
陈子阳
陈安琪
孙垚棋
张继勇
张勇东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202010847763.5A priority Critical patent/CN112051247B/en
Publication of CN112051247A publication Critical patent/CN112051247A/en
Application granted granted Critical
Publication of CN112051247B publication Critical patent/CN112051247B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6456Spatial resolved fluorescence measurements; Imaging
    • G01N21/6458Fluorescence microscopy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6486Measuring fluorescence of biological material, e.g. DNA, RNA, cells
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/50Constructional details
    • H04N23/54Mounting of pick-up tubes, electronic image sensors, deviation or focusing coils

Abstract

The invention provides a lens-free imaging device based on laminated imaging and a phase recovery method thereof. The device comprises a laser, a sample clamping device, a motor, a scattering layer and an image sensor. In the imaging process, a motor is continuously used for moving the position of the sample in the horizontal direction, and the scattering layer is used for modulating emergent light generated after laser irradiates the sample. And simultaneously keeping the relative distance between the sample and the image sensor unchanged, recording images of the sample at different positions by using the image sensor, and recovering the profile of the scattering layer and the profile of the sample by using a two-dimensional cross-correlation algorithm and an rPIE algorithm after 200 images are acquired. The two-dimensional cross-correlation algorithm and the rPIE algorithm were implemented using MATLAB2018a software. According to the invention, a scattering layer is added between the sample and the image sensor, so that the displacement of the sample at different positions can be calculated, and the limit of reaching the resolution is improved from 9-1 to 9-4.

Description

Lens-free imaging device based on laminated imaging and phase recovery method thereof
Technical Field
The invention relates to the field of computational imaging, in particular to a lens-free imaging device based on a laminated imaging structure modulation principle and a phase recovery method thereof.
Technical Field
Lens-free microscopy is a technique that directly brings a sample into close proximity with a photosensor, such as a charge-coupled device (CCD) or a complementary metal-oxide semiconductor Chip (CMOS), and directly images the sample without the need for optical elements. In a classical lensless microscope system, the main components are the light source, the sample and the photosensor, and therefore the resolution of the photosensor will directly limit the resolution of the image. Because no lens is added in the lensless microscope system, the problem that the traditional optical microscope imaging needs to be compromised between a field of view and resolution is well solved, and the problem of aberration caused by the lens can be well avoided. Compared with the traditional analysis and detection instrument which is expensive in price, large in size and needs to be operated by professional technical personnel, the lens-free microtechnical principle is simple, the price is low, the carrying is convenient, and the real-time observation and analysis of the sample by the detection personnel are facilitated.
In 1969, hoppe first proposed the concept of stacked-scan imaging (ptychodography), with the purpose of being able to deterministically solve the phase information of complex-amplitude objects. The laminated scanning coherent diffraction imaging system obtains a plurality of intensity diffraction patterns corresponding to different parts of a sample through a two-dimensional moving aperture diaphragm or the sample. Stack imaging improves the convergence of the phase recovery algorithm and reduces the requirements on experimental conditions by recording a high redundancy of data. With the development of the reconstruction algorithm, the imaging quality of the laminated scanning is remarkably improved, and the application range is continuously expanded. Rodenburg et al in 2004 proposed a robust stacked diffraction imaging algorithm (PIE), using multiple intensity patterns to replace the traditional single intensity pattern, greatly improving the convergence rate and noise immunity of the iterative phase recovery algorithm, and the recovered result has the characteristics of high contrast, high resolution, large field of view, quantitative phase, etc.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a lens-free imaging device based on laminated imaging and a phase recovery method thereof. The method can be used to quickly acquire raw images and calculate amplitude and phase information of an object from the raw images.
A lens-free imaging device based on laminated imaging comprises a laser, a sample clamping device, a motor, a scattering layer and an image sensor.
Wherein, the scattering layer is made by uniformly coating 1-5 μm fluorescent beads on the protective glass on the surface of the image sensor, and one side of the surface of the image sensor is provided with a certain space without coating the fluorescent beads, namely a blank area without the scattering layer. The laser is fixed right above the image sensor and the sample, the vertical distance between the laser and the sample is 10cm, the sample is fixed above the image sensor through a sample clamping device, the distance d1 between the sample and the scattering layer is 400 mu m, the motor is connected with the sample clamping device to control the movement of the sample, and the distance d2 between the scattering layer and the image sensor is fixed at 160 mu m because the scattering layer is coated on the protective glass on the surface of the image sensor.
The laser is green light with the wavelength of 532nm, and the power is 10mW;
further, the image sensor is a CMOS sensor manufactured by IMAGINGSOURSE corporation, the model is DMM37UX226-ML, and the pixel size of the sensor is 1.85 μm.
Further, the resolution of the image sensor is 4000 × 3000, and the size of the region where the fluorescent beads are not coated is 200 × 3000.
In the imaging process, a motor is continuously used for moving the position of the sample in the horizontal direction, and the scattering layer is used for modulating emergent light generated after laser irradiates the sample. And simultaneously keeping the relative distance between the sample and the image sensor unchanged, recording images of the sample at different positions by using the image sensor, and recovering the profile of the scattering layer and the profile of the sample by using a two-dimensional cross-correlation algorithm and an rPIE algorithm after 200 images are acquired. The two-dimensional cross-correlation algorithm and the rPIE algorithm are realized by using MATLAB2018a software
A lens-free imaging phase recovery method based on laminated imaging comprises the following steps:
step 1: using image transferThe sensor collects an image. In the process of acquiring images, two motors are used for moving the position of a sample in the horizontal direction, the images formed by the sample at different positions are recorded by using the image sensor, the number of the images to be acquired is 200, and the number is marked as I j (j =1,2.., 200), the size of a single image is 4000 × 3000.
And 2, step: the relative displacement between the images is calculated.
And calculating the initial moving position of the sample after the motor moves in the horizontal direction each time by adopting a two-dimensional cross-correlation algorithm. The method comprises the following specific steps:
(1) Selecting an image I j (j =1,2,.. 200) corresponds to the area without fluorescent beads. When the fluorescent ball is coated, a certain blank area is left on the surface of the image sensor, so that the information that the small ball does not exist in part of the area on the acquired image is also available, and the information is selected and recorded as I _ blank j (j=1,2,...,200);
(2) Will I _ blank j (j =1,2,.., 200) a first image is registered as a target image and subjected to a two-dimensional fourier transform operation;
(3) Sequentially reading 200 acquired images, respectively carrying out Fourier transformation, and calculating the relative displacement between each image and a target image in the Fourier plane by using a dftreistration () function of MATLAB (matrix laboratory), namely calculating the relative position between the images, wherein the relative displacement of the jth image is recorded as (x) j ,y j )。
And step 3: selecting an image I j (j =1,2,.., 200) the area corresponding to the coated fluorescent beads, i.e., the area where the scattering layer exists, is marked as I _ recovery j (j=1,2,...,200)。
And 4, step 4: the sample is initialized and the solution of the scattering layer to the wavefront is initialized. Two all-1 matrix variables are used as initial solutions, and the initial solution of the sample wavefront and the initial solution of the scattering layer wavefront are denoted as O (x, y) and D (x, y), respectively.
And 5: using I _ recovery j (j =1,2.., 200) with stacked imaging structure modulation;
the description is given through the jth image, and the specific steps are as follows:
(1) Calculating the relative displacement (x, y) of the initial solution O (x, y) of the sample wave front according to the jth image j ,y j ) Translating to obtain solution O of sample wave front j (x, y), the corresponding equation is:
O j (x,y)=O(x-x j ,y-y j );
(2) The solution of the sample wavefront is propagated to the scattering layer. The propagation function used is a point spread function of light propagating in the self-using space, wherein the propagation distance is d1, and the solution after the sample wavefront propagation is recorded as O D (x, y), the corresponding equation is:
O D (x,y)=PSF(d 1 )*O j (x,y);
(3) Multiplying the solution after sample wave front propagation with the initial solution D (x, y) of the scattering layer to obtain a multiplied solution
Figure BDA0002643660530000041
The corresponding formula is:
Figure BDA0002643660530000042
(4) Propagating the multiplied solution to the image sensor plane using a propagation function that is a point spread function of light propagating in the free space with a propagation distance d2, and noting the solution propagated to the image sensor as psi j (x, y), the corresponding equation is:
Figure BDA0002643660530000043
(5) By I _ recovery j Replacement psi j The amplitude of (x, y) corresponds to the formula:
Figure BDA0002643660530000044
(6) The result after replacement is psi' j (x,y) propagates back to the scattering layer, using a propagation function that is the point spread function of light propagating in the dead space, corresponding to the formula:
Figure BDA0002643660530000045
(7) Updating O Using rPIE Algorithm D (x, y) and D (x, y), obtaining the updated solution after the sample wavefront propagation
Figure BDA0002643660530000046
And updated initial solution D of the scattering layer update (x, y) corresponding to the formulas:
Figure BDA0002643660530000047
Figure BDA0002643660530000048
(8) Will be provided with
Figure BDA0002643660530000051
Backward propagating back to object layer to obtain updated solution of sample wavefront
Figure BDA0002643660530000052
The used propagation function is a point spread function of light propagating in the self-use space, and the corresponding formulas are respectively as follows:
Figure BDA0002643660530000053
(9) Will be provided with
Figure BDA0002643660530000054
Translating back to the initial position to obtain the updated initial solution O of the sample wavefront update (x, y) corresponding to the formulas:
Figure BDA0002643660530000055
step 6: and (5) iteratively updating the 200 collected images out of order, wherein the number of iterations is 2.
The invention has the following beneficial effects:
compared with the traditional lens-free imaging device, the lens-free imaging device has the advantages that the scattering layer is added between the sample and the image sensor, wherein the scattering layer is made by coating the fluorescent beads on the image sensor, and meanwhile, a part of the area, on the image sensor, where the scattering layer is not coated is reserved, so that the displacement of the sample at different positions can be calculated. The USAF1951 resolution plate was used for testing, and compared with the previous method, the resolution can reach the limit of 9-1 to 9-4, meanwhile, unstained rat kidney tissue section and stained esophagus cancer section are used for verifying the imaging effect of the system.
Drawings
FIG. 1 is a schematic view of a lensless imaging apparatus;
FIG. 2 is a schematic view of a scattering layer application;
FIG. 3 is pseudo code of a reconstruction algorithm flow;
FIG. 4 is a graph of the results of reconstructing a USAF1951 resolution plate;
FIG. 5 is a graph of the results of phase reconstruction of unstained mouse kidney tissue sections;
FIG. 6 is a graph of the results of reconstructing a large field of view of stained esophageal cancer sections.
Detailed Description
The invention is further described with reference to the following figures and examples;
a lens-free imaging device based on laminated imaging comprises a laser, a sample clamping device, a motor, a scattering layer and an image sensor.
Wherein, the scattering layer is made by uniformly coating 1-5 μm fluorescent beads on the protective glass on the surface of the image sensor, and one side of the surface of the image sensor is provided with a certain space without coating the fluorescent beads, namely a blank area without the scattering layer. The laser is fixed right above the image sensor and the sample, the vertical distance between the laser and the sample is 10cm, the sample is fixed above the image sensor through the sample clamping device, the distance d1 between the sample and the scattering layer is 400 micrometers, the motor is connected with the sample clamping device and controls the movement of the sample, and the distance d2 between the scattering layer and the image sensor is fixed at 160 micrometers as the scattering layer is coated on the protective glass on the surface of the image sensor.
The laser is green light with the wavelength of 532nm, and the power is 10mW;
further, the image sensor is a CMOS sensor manufactured by IMAGINGSOURSE corporation, the model is DMM37UX226-ML, and the pixel size of the sensor is 1.85 μm.
Further, the resolution of the image sensor is 4000 × 3000, and the size of the region where the fluorescent beads are not coated is 200 × 3000.
In the imaging process, a motor is continuously used for moving the position of the sample in the horizontal direction, and the scattering layer is used for modulating emergent light generated after laser irradiates the sample. And simultaneously keeping the relative distance between the sample and the image sensor unchanged, recording images of the sample at different positions by using the image sensor, and recovering the profile of the scattering layer and the profile of the sample by using a two-dimensional cross-correlation algorithm and an rPIE algorithm after 200 images are acquired. The two-dimensional cross-correlation algorithm and the rPIE algorithm are realized by using MATLAB2018a software
A lens-free imaging phase recovery method based on laminated imaging comprises the following steps:
step 1: an image is acquired using an image sensor. In the process of acquiring images, two motors are used for moving the position of a sample in the horizontal direction, the images formed by the sample at different positions are recorded by using the image sensor, the number of the images to be acquired is 200, and the number is marked as I j (j =1,2.., 200), the size of a single image is 4000 × 3000.
And 2, step: the relative displacement between the images is calculated.
And calculating the initial moving position of the sample after the motor moves in the horizontal direction each time by adopting a two-dimensional cross-correlation algorithm. The method comprises the following specific steps:
(1) Selecting an image I j (j =1,2.., 200) corresponds to the area without fluorescent beads. When the fluorescent ball is coated, a certain blank area is left on the surface of the image sensor, so that the information that the small ball does not exist in part of the area on the acquired image is also available, and the information is selected and recorded as I _ blank j (j=1,2,...,200);
(2) Will I _ blank j (j =1,2.., 200) the first image is registered as the target image and subjected to a two-dimensional fourier transform operation;
(3) Sequentially reading 200 acquired images, respectively carrying out Fourier transformation, and calculating the relative displacement between each image and a target image in the Fourier plane by using a dftreistration () function of MATLAB (matrix laboratory), namely calculating the relative position between the images, wherein the relative displacement of the jth image is recorded as (x) j ,y j )。
And 3, step 3: selecting an image I j (j =1,2,.., 200) the area corresponding to the coated fluorescent beads, i.e., the area where the scattering layer exists, is marked as I _ recovery j (j=1,2,...,200)。
And 4, step 4: the sample is initialized and the solution of the scattering layer to the wavefront is initialized. Two all-1 matrix variables are used as initial solutions, and the initial solution of the sample wavefront and the initial solution of the scattering layer wavefront are denoted as O (x, y) and D (x, y), respectively.
And 5: using I _ recovery j (j =1,2.., 200) with stacked imaging structure modulation;
the description is given through the jth image, and the specific steps are as follows:
(1) Calculating the relative displacement (x, y) of the initial solution O (x, y) of the sample wave front according to the jth image j ,y j ) Translating to obtain solution O of sample wave front j (x, y), the corresponding equation is:
O j (x,y)=O(x-x j ,y-y j );
(2) The solution of the sample wavefront is propagated to the scattering layer. The propagation function used is a point spread function of light propagating in the dead space, where the propagation distance is d1, after propagating the sample wavefrontIs solved as O D (x, y), the corresponding equation:
O D (x,y)=PSF(d 1 )*O j (x,y);
(3) Multiplying the solution after sample wave front propagation with the initial solution D (x, y) of the scattering layer to obtain a multiplied solution
Figure BDA0002643660530000081
The corresponding formula is:
Figure BDA0002643660530000082
(4) Propagating the multiplied solution to the image sensor plane using a propagation function that is a point spread function of light propagating in the free space with a propagation distance d2, and noting the solution propagated to the image sensor as psi j (x, y), the corresponding equation is:
Figure BDA0002643660530000083
(5) By I _ recovery j Replacement of psi j The amplitude of (x, y) corresponds to the formula:
Figure BDA0002643660530000084
(6) The result after replacement is psi' j (x, y) propagates back to the scattering layer using a propagation function that is the point spread function of light propagating in the dead space, corresponding to the formula:
Figure BDA0002643660530000085
(7) Updating O Using the rPIE Algorithm D (x, y) and D (x, y), obtaining the updated solution after the sample wavefront propagation
Figure BDA0002643660530000086
And updated initial solution D of the scattering layer update (x, y) corresponding to the formulas:
Figure BDA0002643660530000087
Figure BDA0002643660530000091
(8) Will be provided with
Figure BDA0002643660530000092
Back propagating back to the object layer to obtain an updated solution of the sample wavefront
Figure BDA0002643660530000093
The used propagation function is a point spread function of light propagating in the self-use space, and the corresponding formulas are respectively as follows:
Figure BDA0002643660530000094
(9) Will be provided with
Figure BDA0002643660530000095
Translating back to the initial position to obtain the updated initial solution O of the sample wavefront update (x, y) the corresponding equations are:
Figure BDA0002643660530000096
and 6: and (5) iteratively updating the 200 collected images out of order, wherein the number of iterations is 2.
FIG. 1 is a schematic view of a lensless imaging apparatus. The left figure is a device structure schematic diagram, and the right figure is a device 3D structure schematic diagram. The device mainly comprises four parts, namely laser, a sample, a scattering layer and an image sensor, wherein in the process of actually building the device, the used laser is green light with the wavelength of 532nm, and the power is 10mW; the distance from the laser to the sample is 10cm, the distance d1 from the sample to the scattering layer is about 400 μm, and the distance d2 from the scattering layer to the image sensor is fixed at 160 μm because the scattering layer coats the surface of the image sensor.
Fig. 2 is a schematic view of scattering layer application. The left side is the image sensor corresponding to the small fluorescent ball which is not coated, and the right side is the image sensor coated with the small fluorescent ball. The scattering layer is made by uniformly coating fluorescent small balls with the diameter of 1-5 mu m on the surface of the image sensor, and a certain space is reserved on one side of the image sensor to prevent the small balls from being coated in the process of coating the fluorescent small balls. The image sensor used in this patent is 4000 x 3000 in size, and the area without smeared beads is approximately 200 x 3000 in size.
FIG. 3 is pseudo code of a reconstruction algorithm flow;
FIG. 4 is a graph of the results of reconstructing a USAF1951 resolution plate;
unstained mouse kidney tissue sections and stained esophageal cancer sections were used to verify the imaging effect of the system. FIG. 5 is a graph of the results of phase reconstruction of unstained mouse kidney tissue sections; FIG. 6 is a graph of the results of reconstructing a large field of view of stained esophageal cancer sections.

Claims (4)

1. A lens-free imaging device based on laminated imaging is characterized by comprising a laser, a sample clamping device, a motor, a scattering layer and an image sensor;
wherein, the scattering layer is made by uniformly coating 1-5 μm fluorescent beads on the protective glass on the surface of the image sensor, and one side of the surface of the image sensor is provided with a certain space without coating the fluorescent beads, namely a blank area without the scattering layer; the laser is fixed right above the image sensor and the sample, the vertical distance between the laser and the sample is 10cm, the sample is fixed above the image sensor through a sample clamping device, the distance d1 between the sample and the scattering layer is 400 mu m, the motor is connected with the sample clamping device to control the movement of the sample, and the distance d2 between the scattering layer and the image sensor is fixed at 160 mu m as the scattering layer is coated on the protective glass on the surface of the image sensor;
the laser is green light with the wavelength of 532nm, and the power is 10mW.
2. The lens-free IMAGING device based on stacked IMAGING as claimed in claim 1, wherein the image sensor is a CMOS sensor manufactured by IMAGING SOURSE corporation, model number DMM37UX226-ML, and the pixel size of the sensor is 1.85 μm.
3. The lens-free imaging device according to claim 2, wherein the image sensor has a resolution of 4000 × 3000, and the area not coated with the fluorescent beads has a resolution of 200 × 3000.
4. The phase recovery method of the lens-free imaging device based on the stacked imaging as claimed in claim 1, characterized by comprising the following steps:
step 1: acquiring an image using an image sensor; in the process of acquiring images, two motors are used for moving the position of a sample in the horizontal direction, the image sensors are used for recording images formed by the sample at different positions, the number of the images to be acquired is 200, and the number is marked as I j (j =1,2, …, 200), the size of a single image is 4000 × 3000;
step 2: calculating relative displacement between the images;
calculating the initial moving position of the sample after the motor moves in the horizontal direction each time by adopting a two-dimensional cross-correlation algorithm; the method comprises the following specific steps:
(1) Selecting an image I j (j =1,2, …, 200) corresponds to the area without fluorescent beads; when the fluorescent beads are coated, a certain blank area is left on the surface of the image sensor, so that information that the beads do not exist in partial area is also obtained on the acquired image, and the information is marked as I _ blank j (j=1,2,…,200);
(2) Will I _ blank j (j =1,2.., 200) the first image is registered as the target image and subjected to a two-dimensional fourier transform operation;
(3) Sequentially reading 200 acquired images, respectively carrying out Fourier transformation, and calculating the relative displacement between each image and a target image in the Fourier plane by using a dftreistration () function of MATLAB (matrix laboratory), namely calculating the relative position between the images, wherein the relative displacement of the jth image is recorded as (x) j ,y j );
And step 3: selecting an image I j (j =1,2,.., 200) the area corresponding to the coated fluorescent beads, i.e., the area where the scattering layer exists, is marked as I _ recovery j (j=1,2,...,200);
And 4, step 4: initializing a sample and a solution of a scattering layer corresponding to a wave front; taking two matrix variables of all 1 as initial solutions, and respectively recording the initial solution of the sample wave front and the initial solution of the scattering layer wave front as O (x, y) and D (x, y);
and 5: using I _ recovery j (j =1,2,.., 200) for stack imaging structure modulation;
the description is given through the jth image, and the specific steps are as follows:
(1) Calculating the relative displacement (x, y) of the initial solution O (x, y) of the sample wave front according to the jth image j ,y j ) Translating to obtain solution O of sample wave front j (x, y), the corresponding equation is:
O j (x,y)=O(x-x j ,y-y j );
(2) Propagating a solution of the sample wavefront to the scattering layer; the propagation function used is a point spread function of light propagating in the self-using space, wherein the propagation distance is d1, and the solution after the sample wavefront propagation is recorded as O D (x, y), the corresponding equation is:
O D (x,y)=PSF(d 1 )*O j (x,y);
(3) Multiplying the solution after sample wave front propagation with the initial solution D (x, y) of the scattering layer to obtain a multiplied solution
Figure FDA0003787824360000021
The corresponding formula is:
Figure FDA0003787824360000022
(4) Propagating the multiplied solution to the image sensor plane using a propagation function that is a point spread function of light propagating in the free space with a propagation distance d2, and noting the solution propagated to the image sensor as psi j (x, y), the corresponding equation is:
Figure FDA0003787824360000031
(5) By I _ recovery j Replacement of psi j The amplitude of (x, y) corresponds to the formula:
Figure FDA0003787824360000032
(6) The result after replacement is psi' j (x, y) propagates back to the scattering layer using a propagation function that is the point spread function of light propagating in the dead space, corresponding to the formula:
Figure FDA0003787824360000033
(7) Updating O Using the rPIE Algorithm D (x, y) and D (x, y), obtaining the updated solution after the sample wavefront propagation
Figure FDA0003787824360000034
And updated initial solution D of the scattering layer update (x, y) corresponding to the formulas:
Figure FDA0003787824360000035
Figure FDA0003787824360000036
(8) Will be provided with
Figure FDA0003787824360000037
Backward propagating back to object layer to obtain updated solution of sample wavefront
Figure FDA0003787824360000038
The used propagation function is a point spread function of light propagating in the self-use space, and the corresponding formulas are respectively as follows:
Figure FDA0003787824360000039
(9) Will be provided with
Figure FDA0003787824360000041
Translating back to the initial position to obtain the updated initial solution O of the sample wavefront update (x, y) corresponding to the formulas:
Figure FDA0003787824360000042
step 6: and (5) iteratively updating the 200 collected images out of order, wherein the number of iterations is 2.
CN202010847763.5A 2020-08-21 2020-08-21 Lens-free imaging device based on laminated imaging and phase recovery method thereof Active CN112051247B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010847763.5A CN112051247B (en) 2020-08-21 2020-08-21 Lens-free imaging device based on laminated imaging and phase recovery method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010847763.5A CN112051247B (en) 2020-08-21 2020-08-21 Lens-free imaging device based on laminated imaging and phase recovery method thereof

Publications (2)

Publication Number Publication Date
CN112051247A CN112051247A (en) 2020-12-08
CN112051247B true CN112051247B (en) 2022-10-18

Family

ID=73599648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010847763.5A Active CN112051247B (en) 2020-08-21 2020-08-21 Lens-free imaging device based on laminated imaging and phase recovery method thereof

Country Status (1)

Country Link
CN (1) CN112051247B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102660457A (en) * 2012-04-17 2012-09-12 南昌航空大学 Device and method for analyzing and counting blood cells by lensless holographic diffraction imaging
JP2018120025A (en) * 2017-01-23 2018-08-02 セイコーエプソン株式会社 Lighting system and projector
CN110927115A (en) * 2019-12-09 2020-03-27 杭州电子科技大学 Lens-free dual-type fusion target detection device and method based on deep learning

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1889111A2 (en) * 2005-05-25 2008-02-20 Massachusetts Institute of Technology Multifocal imaging systems and methods

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102660457A (en) * 2012-04-17 2012-09-12 南昌航空大学 Device and method for analyzing and counting blood cells by lensless holographic diffraction imaging
JP2018120025A (en) * 2017-01-23 2018-08-02 セイコーエプソン株式会社 Lighting system and projector
CN110927115A (en) * 2019-12-09 2020-03-27 杭州电子科技大学 Lens-free dual-type fusion target detection device and method based on deep learning

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Machine learning based single-frame super-resolution processing for lensless blood cell counting;Xiwei huang et al.;《Sensors》;20161102;第1-16页 *
Wide-field, high-resolution lensless on-chip;Jiang S et al.;《Lab on a chip》;20200211;第1058-1065页 *
小型化无透镜微流控片上生物成像检测;陈津 等;《物联网学报》;20190630;第9-19页 *
无透镜片上显微成像技术:理论、发展与应用;张佳琳 等;《红外与激光工程 》;20190625;第121-153 *

Also Published As

Publication number Publication date
CN112051247A (en) 2020-12-08

Similar Documents

Publication Publication Date Title
US11946854B2 (en) Systems and methods for two-dimensional fluorescence wave propagation onto surfaces using deep learning
CN106990694B (en) Non-iterative phase recovery device and method under partially-dry-light illumination
CN109884018B (en) Submicron lens-free microscopic imaging method and system based on neural network
JP2018537709A (en) Computer microscope and method for generating images under different illumination conditions
CN106324853B (en) A kind of double object distance lamination imaging methods of visible domain
CN111366557B (en) Phase imaging method based on thin scattering medium
CN109685745B (en) Phase microscopic imaging method based on deep learning
US11482021B2 (en) Adaptive sensing based on depth
JP2013545138A (en) On-chip 4D light field microscope
CN102057397A (en) System and method for producing an optically sectioned image using both structured and uniform illumination
CN106600687B (en) Multi-direction flame emission chromatography system
CN108594418A (en) A kind of light field micro imaging system and its method based on array single pixel detector
CN110675451B (en) Digital self-adaptive correction method and system based on phase space optics
JP2013531268A (en) Measuring distance using coded aperture
JP2022529366A (en) Systems and methods for deep learning-based color holographic microscopes
CN108364342B (en) Light field microscopic system and three-dimensional information reconstruction method and device thereof
CN113039493A (en) System and method for converting holographic microscopy images into various modality microscopy images
CN115144371A (en) Lensless Fourier laminated diffraction tomography microscopic imaging method based on wavelength scanning
CN108537862B (en) Fourier diffraction scanning microscope imaging method with self-adaptive noise reduction function
CN112051247B (en) Lens-free imaging device based on laminated imaging and phase recovery method thereof
CN113504202B (en) Coherent modulation imaging method based on axial translation binary amplitude mask
CN114972284A (en) Lens-free microscopic imaging system and method based on self-supervision deep learning
CN113218914B (en) Non-invasive scattering medium point spread function acquisition device and method
CN110989155B (en) Lens-free microscopic imaging device based on optical filter array and reconstruction method
CN112747822A (en) Three-dimensional imaging system 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
GR01 Patent grant
GR01 Patent grant