CN112379516B - Multi-object-plane simultaneous imaging method based on digital multiplexing lens - Google Patents

Multi-object-plane simultaneous imaging method based on digital multiplexing lens Download PDF

Info

Publication number
CN112379516B
CN112379516B CN202011326887.5A CN202011326887A CN112379516B CN 112379516 B CN112379516 B CN 112379516B CN 202011326887 A CN202011326887 A CN 202011326887A CN 112379516 B CN112379516 B CN 112379516B
Authority
CN
China
Prior art keywords
imaging
lens
ideal
liquid crystal
spatial light
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
CN202011326887.5A
Other languages
Chinese (zh)
Other versions
CN112379516A (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.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
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 Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN202011326887.5A priority Critical patent/CN112379516B/en
Publication of CN112379516A publication Critical patent/CN112379516A/en
Application granted granted Critical
Publication of CN112379516B publication Critical patent/CN112379516B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B30/00Optical systems or apparatus for producing three-dimensional [3D] effects, e.g. stereoscopic images

Abstract

The invention provides a simultaneous imaging method for multiple object planes based on a digital multiplexing lens. The digital multiplexing lens modulated by the liquid crystal spatial light modulator in a programming mode is combined with an axial multi-object-plane simultaneous imaging system formed by an infinite cylinder-length microscope objective front group with a fixed focal length, so that any section can be imaged clearly or a plurality of sections can be imaged clearly simultaneously within a certain depth range of a sample, and an image plane which is offset simultaneously along the axial direction and the transverse direction and has diffraction limited image quality is formed on the same imaging plane. The number of imaging surfaces, the axial interval and the transverse offset distance between the imaging surfaces of the system can be controlled and adjusted in real time through programming. The invention can realize a flexible three-dimensional imaging system with time and space resolution without axial scanning, and has important practical application value in the fields of biomedical imaging, pharmacology, material science and the like.

Description

Multi-object-plane simultaneous imaging method based on digital multiplexing lens
Technical Field
The invention belongs to the field of optical imaging, and particularly relates to a simultaneous imaging method for multiple object planes based on a digital multiplexing lens.
Background
Optical imaging technology has been a major challenge to achieve flexible time and space simultaneous resolution imaging as a direct visualization tool for studying living cell life phenomena. The existing imaging technologies, including confocal laser scanning microscope, two-photon laser scanning microscope, multi-focus multi-photon microscopic imaging technology, optical coherence tomography technology, etc., need to realize horizontal and axial scanning imaging on samples, so that the imaging speed is limited, and the time resolution characteristic is lost.
Although the multi-focal plane parallel imaging method based on the deformed grating can realize multi-section simultaneous imaging, the section imaging capability is mainly concentrated on the lower diffraction order 0 and +/-1, and the diffraction energy of each order is not uniform. Based on the distorted Dammann grating and the multi-object-plane simultaneous imaging system, the problem of non-uniformity of cross-section imaging energy distribution is solved, more cross-section imaging is realized, but complex processing technology and detection technical means are involved, the flexibility and the real-time performance are not enough, and when the number of cross sections and the cross-section interval are changed, the distorted Dammann grating needs to be processed again.
Disclosure of Invention
The invention inlays ideal imaging lenses with different focal lengths and blazed gratings with different diffraction angles together, and generates a digital multiplexing lens by programming of a liquid crystal spatial light modulator (PLUTO). A method for simultaneously imaging multiple object planes based on a digital multiplexing lens is provided. The imaging method can enable any section to be imaged clearly within a certain depth range of a sample, or a plurality of sections to be imaged clearly at the same time, the imaging quality is the same, and the number of the section images, the section imaging interval and the section transverse offset distance are programmable and controllable in real time through a liquid crystal spatial light modulator (PLUTO). The invention realizes the function of the digital multiplexing lens through the liquid crystal spatial light modulator, thereby respectively realizing the focal lengths of the plurality of ideal imaging lenses and the diffraction angles of the plurality of blazed gratings.
The technical scheme of the invention is a simultaneous imaging method of multiple object planes based on a digital multiplexing lens, which comprises the following steps:
step 1: and constructing a multi-object-plane simultaneous imaging system model of the digital multiplexing lens.
Step 2: according to the constructed multi-object-plane simultaneous imaging system model and the combination rule of an ideal optical system, the control of the number of cross-section imaging and the cross-section imaging distance is realized by combining the focal length, the magnification, the image side main plane position, the object side main plane position of the imaging system, the object distance variation of the combined imaging system and the change distance of the object side main plane of the combined imaging system with the focal length of an ideal imaging lens.
And step 3: establishing a phase model of an ideal imaging lens in the liquid crystal spatial light modulator, and establishing a phase model of a blazed grating in the liquid crystal spatial light modulator; a different diffraction angle for each imaging section is achieved.
And 4, step 4: the method comprises the steps of setting the number of ideal imaging lenses and the changed distance of an object side main plane of a combined system according to the requirements of the number of imaging sections and the imaging interval, calculating the focal length of the ideal imaging lenses by combining the set number of the ideal imaging lenses and the changed distance of the object side main plane of the combined system with the step 2, calculating the diffraction angle of a blazed grating by the step 3, constructing a phase model of the digital multiplexing lens on a liquid crystal spatial light modulator, and realizing the design functions of the ideal imaging lenses and the blazed grating in the digital multiplexing lens by the programmed modulation of the liquid crystal spatial light modulator.
Preferably, the multiple object plane simultaneous imaging system model of the digital multiplexing lens in step 1 includes:
the digital multiplexing lens is positioned on the back focal plane of the front group of the infinite cylinder length microobjective;
the digital multiplexing lens is composed of a plurality of ideal imaging lenses and a plurality of blazed gratings; the plurality of ideal imaging lenses are ideal imaging lenses with different focal lengths, and the plurality of blazed gratings are blazed gratings with different diffraction angles; the number of the ideal imaging lens is equal to that of the blazed grating, the ideal imaging lens and the blazed grating are both n, n is more than or equal to 1, and n is an integer;
the ideal imaging lens and the blazed grating are embedded together and are generated by programming modulation of a liquid crystal spatial light modulator;
preferably, the focal length of the combined imaging system in step 2 is:
f′zi=f'obj
wherein, f'objIs set as the focal length of the front group of the microscope objective lens with infinite tube length, f'ziOf combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, the focal length f 'of the imaging system is combined'ziKeeping the same, wherein n is the number of ideal imaging lenses;
and 2, the magnification of the combined imaging system is as follows:
βzi=β
wherein beta is the magnification of the front group of the microscope objective with infinite cylinder length, betaziMagnification for combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the ith ideal imaging lensFocal length f of mirrori'how to vary, magnification of combined imaging system does not follow f'iN is the number of ideal imaging lenses;
the image space main plane position of the combined imaging system in the step 2 is as follows:
H'H′zi=0
wherein H ' is the image side main plane, H ' of the front group of the infinite tube length microscope objective lens 'ziIs the image-side principal plane of the combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, H'ziAlways coincides with the image-side principal plane H 'of the objective lens front group, i.e. the image plane of the combined imaging system remains unchanged and does not follow f'iN is the number of ideal imaging lenses;
the main plane position of the object space of the combined imaging system in the step 2 is as follows:
Figure BDA0002794031220000031
h is the principal plane of object space of the front group of the microscope objective with infinite tube length, HziIs the object-side principal plane of the combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Of the combined imaging system, object-side principal plane with f'iN is the number of ideal imaging lenses;
step 2 object distance position with f 'of combined imaging system'iThe change is as follows:
Figure BDA0002794031220000032
wherein lziFor combining the object distances of the imaging system, l' is the image distance of the front group of the microscope objective with infinite tube length, and the formula shows that the combination is combinedThe distance of the change of the object space of the image system is equal to the distance of the change of the main plane of the object space of the combined system, and n is the number of ideal imaging lenses;
the changed distance of the object space main plane of the combined system in the step 2 is as follows:
Figure BDA0002794031220000033
wherein, Δ ziThe distance changed by the main plane of the object space of the combined system after the ith imaging lens is loaded is n, and the n is the number of ideal imaging lenses;
calculating the focal length of the ideal imaging lens in the step 2 according to the changed distance of the object main plane of the combined system in the step 2;
by varying the number of ideal imaging lenses, i.e. n and fi' can be used for controlling the number of section imaging and the section imaging interval, i belongs to [1, n ]]
By changing the corresponding diffraction angle theta of each imaging section in the x and y directionsx-i,θy-iThe transverse offset distance of the cross section and the corresponding diffraction angle theta in the x and y directions can be controlledx-i,θy-iThe method is realized by blazed gratings;
preferably, the step 3 of establishing a phase model of the ideal imaging lens in the liquid crystal spatial light modulator is as follows:
the digital multiplexing lens is a phase diagram generated by a liquid crystal spatial light modulator, and according to the Fourier optics theory, the phase modulation amount of the ith ideal imaging lens on the imaging wavefront is as follows:
Figure BDA0002794031220000041
where λ is the incident wavelength, x, y are coordinates with the center of the lens as the origin, and fi' is the focal length of the ith ideal imaging lens, and n is the number of the ideal imaging lenses;
setting the pixel resolution of the liquid crystal spatial light modulator as M x N, the pixel center distance as a, M being the number of rows of the liquid crystal spatial light modulator pixels, and N being the number of columns of the liquid crystal spatial light modulator pixels;
the center of the liquid crystal spatial light modulator is taken as the origin of coordinates, and the focal length is fiThe phase model of the ideal imaging lens of' in the liquid crystal spatial light modulator can be expressed as:
Figure BDA0002794031220000042
Figure BDA0002794031220000043
in the formula (mod)Representing a 2 pi operation, k being a first coefficient, l being a second coefficient, a being a pixel center-to-center distance, fi' is the focal length of the ith ideal imaging lens, M is the number of rows of the liquid crystal spatial light modulator pixels, N is the number of columns of the liquid crystal spatial light modulator pixels, and N is the number of the ideal imaging lenses;
and 3, establishing a phase model of the blazed grating in the liquid crystal spatial light modulator:
in order to realize that different axial section images are not overlapped in different areas of an image plane, emergent light beams of ideal imaging lenses with different focal lengths have different diffraction angles;
is realized by using a liquid crystal spatial light modulator to generate a blazed grating with a focal length fi' the outgoing beam of the i-th ideal imaging lens, which is loaded with the phase model distribution of the corresponding i-th blazed grating, is represented as:
Figure BDA0002794031220000051
Figure BDA0002794031220000052
in the formula, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe method comprises the steps that a grating period corresponding to the ith ideal imaging lens in the y direction by taking a pixel as a unit is defined, k is a first coefficient, l is a second coefficient, M is the number of rows of pixels of the liquid crystal spatial light modulator, N is the number of columns of the pixels of the liquid crystal spatial light modulator, and N is the number of the ideal imaging lenses;
step 3, each imaging section has different diffraction angles:
after loading the ith blazed grating, corresponding to fiThe diffraction angles of the imaging cross-section light beam in the x and y directions are respectively as follows:
Figure BDA0002794031220000053
wherein a is the pixel center-to-center distance, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens in the y direction by taking the pixel as a unit, and n is the number of the ideal imaging lenses;
the magnitude of the diffraction angle is respectively equal to Tx-iAnd Ty-iIn relation, the direction of the diffraction angle depends on the sign of k, l;
preferably, the phase model of the digital multiplexing lens in the liquid crystal spatial light modulator in step 4 is:
Figure BDA0002794031220000054
Figure BDA0002794031220000055
where a is the pixel center-to-center distance, fi' is the focal length of the ith ideal imaging lens, M is the number of rows of liquid crystal spatial light modulator pixels, N is the number of columns of liquid crystal spatial light modulator pixels, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens in the unit of pixel in the y directionK is the first coefficient and l is the second coefficient.
The invention inlays ideal imaging lenses with different focal lengths and blazed gratings with different diffraction angles together, and generates a digital multiplexing lens by programming of a liquid crystal spatial light modulator (PLUTO). A method for simultaneously imaging multiple object planes based on a digital multiplexing lens is provided. The imaging method can enable any section to be imaged clearly within a certain depth range of a sample, or a plurality of sections to be imaged clearly at the same time, the imaging quality is the same, and the number of the section images, the section imaging interval and the section transverse offset distance are programmable and controllable in real time through a liquid crystal spatial light modulator (PLUTO). The invention realizes the function of the digital multiplexing lens through the liquid crystal spatial light modulator, thereby respectively realizing the focal lengths of the plurality of ideal imaging lenses and the diffraction angles of the plurality of blazed gratings.
The invention can realize a flexible three-dimensional imaging system with time and space resolution without axial scanning, and has important practical application value in the fields of biomedical imaging, pharmacology, material science and the like.
Drawings
FIG. 1: a multi-object plane simultaneous imaging system of an digital multiplexing lens.
FIG. 2: focal length schematic of the combined imaging system.
FIG. 3: schematic diagram of two-dimensional matrix.
FIG. 4: and (5) calibrating a phase gray scale image.
FIG. 5: the method of the invention is a flow chart.
Detailed Description
In order to facilitate the understanding and implementation of the present invention for those of ordinary skill in the art, the present invention is further described in detail with reference to the accompanying drawings and examples, it is to be understood that the embodiments described herein are merely illustrative and explanatory of the present invention and are not restrictive thereof.
The liquid crystal spatial light modulator adopts a PLUTO-VIS type liquid crystal spatial light modulator produced by HOLOEYE company in Germany, the image plane size is 15.36mm multiplied by 8.64mm, the resolution MXN is 1920 multiplied by 1080, the pixel size is 8 mu M, the interval between pixels is 1 mu M, and the size a of a single pixel is 9 mu M
The following describes an embodiment of the present invention with reference to fig. 1 to 5 as follows:
a method for simultaneously imaging multiple object planes based on a digital multiplexing lens comprises the following specific steps:
step 1: a model of a multi-object plane simultaneous imaging system of an digitally multiplexed lens is constructed as shown in fig. 1.
The multiple object plane simultaneous imaging system model of the digital multiplexing lens in the step 1 comprises:
the digital multiplexing lens is positioned on the back focal plane of the front group of the infinite cylinder length microobjective;
the multi-object-plane simultaneous imaging system of the digital multiplexing lens adopts helium neon laser as an illumination light source, the wavelength lambda is 0.6328 mu, and the focal length f 'of a front group of an infinite cylinder length objective lens is selected'obj25mm, numerical aperture NA 0.25, multiplying factor beta 10×
The digital multiplexing lens is composed of a plurality of ideal imaging lenses and a plurality of blazed gratings; the plurality of ideal imaging lenses are ideal imaging lenses with different focal lengths, and the plurality of blazed gratings are blazed gratings with different diffraction angles; the number of the ideal imaging lenses is equal to that of the blazed grating, n is 9, n is larger than or equal to 1, and n is an integer;
the ideal imaging lens and the blazed grating are embedded together and are generated by programming modulation of a liquid crystal spatial light modulator;
step 2: according to the constructed multi-object-plane simultaneous imaging system model and the combination rule of an ideal optical system, the control of the number of cross-section imaging and the cross-section imaging distance is realized by combining the focal length, the magnification, the image side main plane position, the object side main plane position of the imaging system, the object distance variation of the combined imaging system and the change distance of the object side main plane of the combined imaging system with the focal length of an ideal imaging lens.
As schematically shown in fig. 2;
step 2, the focal length of the combined imaging system is as follows:
f′zi=f'obj
wherein, f'objIs set as the focal length of the front group of the microscope objective lens with infinite tube length, f'ziOf combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, the focal length f 'of the imaging system is combined'ziKeeping the same, wherein n is the number of ideal imaging lenses;
focal length f 'of the ith ideal imaging lens'iAccording to the requirements of the cross-section imaging distance and the changed distance of the main plane of the object space of the combined system in the step 2 in the invention content;
wherein, for P0Cross section,. DELTA.z0=0,f0' -0, and P1To P8Cross section, set requirement Δ z1=Δz2=…=Δz80.1mm as shown in figure 1;
f′1=6250mm,f′2=3125mm,f′3=2083.33mm,f′4=1562.5mm
f′5=1250mm,f′6=1041.67mm,f′7=892.85mm,f′8=781.25mm
and 2, the magnification of the combined imaging system is as follows:
βzi=β
wherein beta is the magnification of the front group of the microscope objective with infinite cylinder length, betaziMagnification for combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow varied, the magnification of the combined imaging system is not f'iN is the number of ideal imaging lenses;
the image space main plane position of the combined imaging system in the step 2 is as follows:
H'H′zi=0
wherein H' is an infinite cylinder length micro-objectImage-side principal plane, H 'of the front mirror group'ziIs the image-side principal plane of the combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, H'ziAlways coincides with the image-side principal plane H 'of the objective lens front group, i.e. the image plane of the combined imaging system remains unchanged and does not follow f'iN is the number of ideal imaging lenses;
the main plane position of the object space of the combined imaging system in the step 2 is as follows:
Figure BDA0002794031220000081
h is the principal plane of object space of the front group of the microscope objective with infinite tube length, HziIs the object-side principal plane of the combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'i i∈[1,n]Of the combined imaging system, object-side principal plane with f'iN is the number of ideal imaging lenses;
step 2 object distance position with f 'of combined imaging system'iThe change is as follows:
Figure BDA0002794031220000082
wherein is 1ziThe object distance of the combined imaging system, l' is the image distance of the front group of the infinite tube-length microscope objective, the formula shows that the distance of the object space change of the combined imaging system is equal to the distance of the object space main plane change of the combined imaging system, and n is the number of ideal imaging lenses;
the changed distance of the object space main plane of the combined system in the step 2 is as follows:
Figure BDA0002794031220000091
wherein, Δ ziThe distance changed by the main plane of the object space of the combined system after the ith imaging lens is loaded is n, and the n is the number of ideal imaging lenses;
calculating the focal length of the ideal imaging lens in the step 2 according to the changed distance of the object main plane of the combined system in the step 2;
by varying the number of ideal imaging lenses, i.e. n and fi' can be used for controlling the number of section imaging and the section imaging interval, i belongs to [1, n ]]
By changing the corresponding diffraction angle theta of each imaging section in the x and y directionsx-i,θy-iThe transverse offset distance of the cross section and the corresponding diffraction angle theta in the x and y directions can be controlledx-i,θy-iThe method is realized by blazed gratings;
and step 3: establishing a phase model of an ideal imaging lens in the liquid crystal spatial light modulator, and establishing a phase model of a blazed grating in the liquid crystal spatial light modulator; a different diffraction angle for each imaging section is achieved.
Step 3, establishing a phase model of the ideal imaging lens in the liquid crystal spatial light modulator is as follows:
the digital multiplexing lens is a phase diagram generated by a liquid crystal spatial light modulator, and according to the Fourier optics theory, the phase modulation amount of the ith ideal imaging lens on the imaging wavefront is as follows:
Figure BDA0002794031220000092
where λ is the incident wavelength, x, y are coordinates with the center of the lens as the origin, and fi' is the focal length of the ith ideal imaging lens, and n is the number of the ideal imaging lenses;
setting the pixel resolution of the liquid crystal spatial light modulator as M x N, the pixel center distance as a, M being the number of rows of the liquid crystal spatial light modulator pixels, and N being the number of columns of the liquid crystal spatial light modulator pixels;
the center of the liquid crystal spatial light modulator is taken as the origin of coordinates, and the focal length is fiThe phase model of the ideal imaging lens of' in the liquid crystal spatial light modulator can be expressed as:
Figure BDA0002794031220000093
Figure BDA0002794031220000101
in the formula (mod)Representing a 2 pi operation, k being a first coefficient, l being a second coefficient, a being a pixel center-to-center distance, fi' is the focal length of the ith ideal imaging lens, M is the number of rows of the liquid crystal spatial light modulator pixels, N is the number of columns of the liquid crystal spatial light modulator pixels, and N is the number of the ideal imaging lenses;
and 3, establishing a phase model of the blazed grating in the liquid crystal spatial light modulator:
in order to realize that different axial section images are not overlapped in different areas of an image plane, emergent light beams of ideal imaging lenses with different focal lengths have different diffraction angles;
is realized by using a liquid crystal spatial light modulator to generate a blazed grating with a focal length fi' the outgoing beam of the i-th ideal imaging lens, which is loaded with the phase model distribution of the corresponding i-th blazed grating, is represented as:
Figure BDA0002794031220000102
Figure BDA0002794031220000103
in the formula, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iIs the grating period corresponding to the ith ideal imaging lens in the y direction by taking the pixel as a unit, k is a first coefficient, l is a second coefficient, and M is the pixel of the liquid crystal spatial light modulatorThe number of rows, N the number of columns of liquid crystal spatial light modulator pixels, and N the number of ideal imaging lenses;
step 3, each imaging section has different diffraction angles:
after loading the ith blazed grating, corresponding to fiThe diffraction angles of the imaging cross-section light beam in the x and y directions are respectively as follows:
Figure BDA0002794031220000104
wherein a is the pixel center-to-center distance, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens in the y direction by taking the pixel as a unit, and n is the number of the ideal imaging lenses;
the magnitude of the diffraction angle is respectively equal to Tx-iAnd Ty-iIn relation, the direction of the diffraction angle depends on the sign of k, l;
and 4, step 4: the method comprises the steps of setting the number of ideal imaging lenses and the changed distance of an object side main plane of a combined system according to the requirements of the number of imaging sections and the imaging interval, calculating the focal length of the ideal imaging lenses by combining the set number of the ideal imaging lenses and the changed distance of the object side main plane of the combined system with the step 2, calculating the diffraction angle of a blazed grating by the step 3, constructing a phase model of the digital multiplexing lens on a liquid crystal spatial light modulator, and realizing the design functions of the ideal imaging lenses and the blazed grating in the digital multiplexing lens by the programmed modulation of the liquid crystal spatial light modulator.
Step 4, the phase model of the digital multiplexing lens in the liquid crystal spatial light modulator is as follows:
Figure BDA0002794031220000111
Figure BDA0002794031220000112
where a is the pixel center-to-center distance, fiIs the focal length of the ith ideal imaging lens, M is the number of rows of liquid crystal spatial light modulator pixels, N is the number of columns of liquid crystal spatial light modulator pixels, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens is in the y direction by taking a pixel as a unit, k is a first coefficient, and l is a second coefficient.
The phase distribution of the digital multiplex lens is shown in table 1:
TABLE 1 phase distribution of digitally multiplexed lenses
Figure BDA0002794031220000121
And 5: matlab is called to generate a two-dimensional matrix, and 9 phase distributions are randomly assigned with initial values to the two-dimensional matrix, as shown in FIG. 3.
Step 6: as shown in fig. 4, which is a scaled phase grayscale diagram, the phase modulation amount corresponding to each pixel point in the matrix is calculated according to the phase expression in table 1, and a grayscale matrix is generated.
The PLUTO works by changing the addressing voltage to control the phase modulation amount of each pixel, the manufacturer of PLUTO already maps the driving voltage to the gray value displayed by the computer when the PLUTO leaves the factory, the manufacturer gives a (gray-phase) customer lookup table, but specifically, each device needs to be recalibrated to obtain the gray-phase calibration curve shown in fig. 4, and the gray corresponding to the curve is loaded on the driving software of the PLUTO to realize the required phase modulation.
And 7: a gray scale map is generated using the Imwrite function and loaded onto the spatial light modulator.
It should be understood that parts of the specification not set forth in detail are well within the prior art.
It should be understood that the above description of the preferred embodiments is given for clarity and not for any purpose of limitation, and that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (3)

1. A multi-object-plane simultaneous imaging method based on a digital multiplexing lens is characterized by comprising the following steps:
step 1: constructing a multi-object-plane simultaneous imaging system model of the digital multiplexing lens;
step 2: according to the constructed multi-object-plane simultaneous imaging system model and the combination rule of an ideal optical system, the control of the number of cross-section imaging and the cross-section imaging distance is realized by combining the focal length, the magnification, the image side main plane position, the object side main plane position of the imaging system, the object distance variation of the combined imaging system and the change distance of the object side main plane of the combined imaging system with the focal length change of an ideal imaging lens;
and step 3: establishing a phase model of an ideal imaging lens in the liquid crystal spatial light modulator, and establishing a phase model of a blazed grating in the liquid crystal spatial light modulator; different diffraction angles of each imaging section are realized;
and 4, step 4: setting the number of ideal imaging lenses and the changed distance of the object side main plane of the combined system according to the requirements of the number of imaging sections and the imaging interval, calculating the focal length of the ideal imaging lenses by combining the set number of the ideal imaging lenses and the changed distance of the object side main plane of the combined system with the step 2, calculating the diffraction angle of the blazed grating by the step 3, constructing a phase model of the digital multiplexing lens on the liquid crystal spatial light modulator, and realizing the design functions of the ideal imaging lenses and the blazed grating in the digital multiplexing lens by the programmed modulation of the liquid crystal spatial light modulator;
the multiple object plane simultaneous imaging system model of the digital multiplexing lens in the step 1 comprises:
the digital multiplexing lens is positioned on the back focal plane of the front group of the infinite cylinder length microobjective;
the digital multiplexing lens is composed of a plurality of ideal imaging lenses and a plurality of blazed gratings; the plurality of ideal imaging lenses are ideal imaging lenses with different focal lengths, and the plurality of blazed gratings are blazed gratings with different diffraction angles; the number of the ideal imaging lens is equal to that of the blazed grating, the ideal imaging lens and the blazed grating are both n, n is more than or equal to 1, and n is an integer;
the ideal imaging lens and the blazed grating are embedded together and are generated by programming modulation of a liquid crystal spatial light modulator;
step 2, the focal length of the combined imaging system is as follows:
f′zi=f′obj
wherein, f'objIs set as the focal length of the front group of the microscope objective lens with infinite tube length, f'ziOf combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'ii∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, the focal length f 'of the imaging system is combined'ziKeeping the same, wherein n is the number of ideal imaging lenses;
and 2, the magnification of the combined imaging system is as follows:
βzi=β
wherein beta is the magnification of the front group of the microscope objective with infinite cylinder length, betaziMagnification for combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'ii∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow varied, the magnification of the combined imaging system is not f'iN is the number of ideal imaging lenses;
the image space main plane position of the combined imaging system in the step 2 is as follows:
H'H′zi=0
wherein H ' is the image side main plane, H ' of the front group of the infinite tube length microscope objective lens 'ziIs the image-side principal plane of the combined imaging system, i.e. at focal length f'objImage space focal plane superposition focal of front group of infinite tube length microscope objectiveDistance is f'ii∈[1,n]Irrespective of the focal length f 'of the i-th ideal imaging lens'iHow to vary, H'ziAlways coincides with the image-side principal plane H 'of the objective lens front group, i.e. the image plane of the combined imaging system remains unchanged and does not follow f'iN is the number of ideal imaging lenses;
the main plane position of the object space of the combined imaging system in the step 2 is as follows:
Figure FDA0003555344610000021
i∈[1,n]
h is the principal plane of object space of the front group of the microscope objective with infinite tube length, HziIs the object-side principal plane of the combined imaging system, i.e. at focal length f'objF 'is superposed on the image focal plane of the front group of the infinite tube length microscope objective lens'ii∈[1,n]Of the combined imaging system, object-side principal plane with f'iN is the number of ideal imaging lenses;
step 2 object distance position with f 'of combined imaging system'iThe change is as follows:
Figure FDA0003555344610000022
i∈[1,n]
wherein lziThe distance of the object space of the combined imaging system is l' is the image distance of the front group of the infinite tube-length microscope objective, the formula shows that the distance of the object space change of the combined imaging system is equal to the distance of the change of the main plane of the object space of the combined system, and n is the number of ideal imaging lenses;
the changed distance of the object space main plane of the combined system in the step 2 is as follows:
Figure FDA0003555344610000031
i∈[1,n]
wherein, Δ ziThe distance changed by the main plane of the object space of the combined system after the ith imaging lens is loaded is n, and the n is the number of ideal imaging lenses;
calculating the focal length of the ideal imaging lens in the step 2 according to the changed distance of the object main plane of the combined system in the step 2;
by varying the number of ideal imaging lenses, i.e. n and f'iCan be used for controlling the number of cross-section imaging and the cross-section imaging interval, i belongs to [1, n ]]
By changing the corresponding diffraction angle theta of each imaging section in the x and y directionsx-i,θy-iThe transverse offset distance of the cross section and the corresponding diffraction angle theta in the x and y directions can be controlledx-i,θy-iBy means of blazed gratings.
2. The method of claim 1, wherein the method comprises:
step 3, establishing a phase model of the ideal imaging lens in the liquid crystal spatial light modulator is as follows:
the digital multiplexing lens is a phase diagram generated by a liquid crystal spatial light modulator, and according to the Fourier optics theory, the phase modulation amount of the ith ideal imaging lens on the imaging wavefront is as follows:
Figure FDA0003555344610000032
i∈[1,n]
wherein λ is an incident wavelength, x, y are coordinates with the center of the lens as an origin, and f'iIs the focal length of the ith ideal imaging lens, and n is the number of the ideal imaging lenses;
setting the pixel resolution of the liquid crystal spatial light modulator as M x N, the pixel center distance as a, M being the number of rows of the liquid crystal spatial light modulator pixels, and N being the number of columns of the liquid crystal spatial light modulator pixels;
the focal length is f 'by taking the center of the liquid crystal spatial light modulator as the origin of coordinates'iThe phase model of the ideal imaging lens in the liquid crystal spatial light modulator can be expressed as:
Figure FDA0003555344610000033
Figure FDA0003555344610000041
i∈[1,n]
in the formula (mod)Representing 2 pi complementation operation, k being a first coefficient, l being a second coefficient, a being a pixel center distance, f'iThe focal length of the ith ideal imaging lens is, M is the number of rows of the pixels of the liquid crystal spatial light modulator, N is the number of columns of the pixels of the liquid crystal spatial light modulator, and N is the number of the ideal imaging lenses;
and 3, establishing a phase model of the blazed grating in the liquid crystal spatial light modulator:
in order to realize that different axial section images are not overlapped in different areas of an image plane, emergent light beams of ideal imaging lenses with different focal lengths have different diffraction angles;
realized by generating a blazed grating with a liquid crystal spatial light modulator, for a focal length of f'iThe loaded phase model distribution of the corresponding ith blazed grating of the outgoing beam of the ith ideal imaging lens is expressed as:
Figure FDA0003555344610000042
Figure FDA0003555344610000043
i∈[1,n]
in the formula, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iIs a y directionIn the grating period corresponding to the ith ideal imaging lens taking a pixel as a unit, k is a first coefficient, l is a second coefficient, M is the number of rows of pixels of the liquid crystal spatial light modulator, N is the number of columns of the pixels of the liquid crystal spatial light modulator, and N is the number of the ideal imaging lenses;
step 3, each imaging section has different diffraction angles:
f 'after loading the ith blazed grating'iThe diffraction angles of the imaging cross-section light beams in the x and y directions are respectively as follows:
Figure FDA0003555344610000044
i∈[1,n]
wherein a is the pixel center-to-center distance, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens is in the y direction by taking a pixel as a unit, and n is the number of the ideal imaging lenses;
the magnitude of the diffraction angle is respectively equal to Tx-iAnd Ty-iIn this regard, the direction of the diffraction angle depends on the sign of k, l.
3. The method of claim 1, wherein the method comprises:
step 4, the phase model of the digital multiplexing lens in the liquid crystal spatial light modulator is as follows:
Figure FDA0003555344610000051
Figure FDA0003555344610000052
i∈[1,n]
wherein a is pixel center distance, f'iIs the focal length of the ith ideal imaging lens,m is the number of rows of liquid crystal spatial light modulator pixels, N is the number of columns of liquid crystal spatial light modulator pixels, Tx-iIs the grating period, T, corresponding to the ith ideal imaging lens in the x direction by taking the pixel as a unity-iThe grating period corresponding to the ith ideal imaging lens is in the y direction by taking a pixel as a unit, k is a first coefficient, and l is a second coefficient.
CN202011326887.5A 2020-11-23 2020-11-23 Multi-object-plane simultaneous imaging method based on digital multiplexing lens Active CN112379516B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011326887.5A CN112379516B (en) 2020-11-23 2020-11-23 Multi-object-plane simultaneous imaging method based on digital multiplexing lens

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011326887.5A CN112379516B (en) 2020-11-23 2020-11-23 Multi-object-plane simultaneous imaging method based on digital multiplexing lens

Publications (2)

Publication Number Publication Date
CN112379516A CN112379516A (en) 2021-02-19
CN112379516B true CN112379516B (en) 2022-04-26

Family

ID=74587336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011326887.5A Active CN112379516B (en) 2020-11-23 2020-11-23 Multi-object-plane simultaneous imaging method based on digital multiplexing lens

Country Status (1)

Country Link
CN (1) CN112379516B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114371549B (en) * 2021-12-27 2022-10-28 华中科技大学 Quantitative phase imaging method and system based on multi-focus multiplexing lens

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60186801A (en) * 1984-03-07 1985-09-24 Olympus Optical Co Ltd Multiplex lens structure
EP0523033A1 (en) * 1991-07-10 1993-01-13 IMS Ionen Mikrofabrikations Systeme Gesellschaft m.b.H. Ion optical imaging system
CN102628970A (en) * 2012-04-12 2012-08-08 中国科学院上海光学精密机械研究所 Distorted Dammann grating and system for simultaneously imaging multiple object planes
CN104122609A (en) * 2014-07-29 2014-10-29 深圳大学 Method for computing actual focal length of variable-focal-length lens based on liquid crystal spatial light modulator
CN110441311A (en) * 2019-07-22 2019-11-12 中国科学院上海光学精密机械研究所 The multifocal camera lens of multiaxis for the imaging of more object planes
CN111474611A (en) * 2020-06-05 2020-07-31 浙江舜宇光学有限公司 Super lens group and super imaging device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8356900B2 (en) * 2006-01-20 2013-01-22 Clarity Medical Systems, Inc. Large diopter range real time sequential wavefront sensor

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60186801A (en) * 1984-03-07 1985-09-24 Olympus Optical Co Ltd Multiplex lens structure
EP0523033A1 (en) * 1991-07-10 1993-01-13 IMS Ionen Mikrofabrikations Systeme Gesellschaft m.b.H. Ion optical imaging system
CN102628970A (en) * 2012-04-12 2012-08-08 中国科学院上海光学精密机械研究所 Distorted Dammann grating and system for simultaneously imaging multiple object planes
CN104122609A (en) * 2014-07-29 2014-10-29 深圳大学 Method for computing actual focal length of variable-focal-length lens based on liquid crystal spatial light modulator
CN110441311A (en) * 2019-07-22 2019-11-12 中国科学院上海光学精密机械研究所 The multifocal camera lens of multiaxis for the imaging of more object planes
CN111474611A (en) * 2020-06-05 2020-07-31 浙江舜宇光学有限公司 Super lens group and super imaging device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于ZEMAX 的一款显微镜附加镜头的光学设计;曾爱云等;《蚌埠学院学报》;20190430;105-107页 *

Also Published As

Publication number Publication date
CN112379516A (en) 2021-02-19

Similar Documents

Publication Publication Date Title
Tao et al. Live imaging using adaptive optics with fluorescent protein guide-stars
US7742213B2 (en) Adaptive-scanning optical microscope
Baddeley et al. Three-dimensional sub-100 nm super-resolution imaging of biological samples using a phase ramp in the objective pupil
Girkin et al. Adaptive optics for deeper imaging of biological samples
Ji et al. Adaptive optics via pupil segmentation for high-resolution imaging in biological tissues
WO2006081031A2 (en) Adaptive scanning optical microscope
CN110632726B (en) Method and microscope for calibrating a phase mask
JP2016507078A (en) Depth of field 3D imaging SLM microscope
CN110378473B (en) Phase chromatography method and device based on deep learning and random pattern
JP5661221B2 (en) Method and optical device for super-resolution localization of particles
US11327288B2 (en) Method for generating an overview image using a large aperture objective
EP3918402B1 (en) Calibration of a light-field imaging system
CN105026916A (en) Fluorescence detection device and fluorescence detection method
JP2006031017A (en) Method to grasp at least one sample region using optical raster microscope
WO2018172161A1 (en) Microscope and method for imaging an object
Philipp et al. Diffraction-limited axial scanning in thick biological tissue with an aberration-correcting adaptive lens
CN112379516B (en) Multi-object-plane simultaneous imaging method based on digital multiplexing lens
Vladimirov et al. Dual-view light-sheet imaging through a tilted glass interface using a deformable mirror
US9930241B2 (en) Image processing apparatus, image processing program, and image processing method
Luo et al. Area scanning method for 3D surface profilometry based on an adaptive confocal microscope
CN210922541U (en) Double-camera parallel confocal differential microscopic 3D morphology measuring device
Ren et al. Aberration-free large-area stitch-free 3D nano-printing based on binary holography
Hung et al. SOLEIL: single-objective lens inclined light sheet localization microscopy
Rai et al. Quantitative analysis of illumination and detection corrections in adaptive light sheet fluorescence microscopy
Malivert et al. Active image optimization for lattice light sheet microscopy in thick samples

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