[go: up one dir, main page]

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
light modulator
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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Lenses (AREA)
  • Liquid Crystal (AREA)

Abstract

本发明提出了一种基于数字复用透镜的多物面同时成像方法。由液晶空间光调制器编程调制的数字复用透镜,结合固定焦距的无限筒长显微物镜前组构成的轴向多物面同时成像系统,能够在样品的一定深度范围内,使任意截面成像清晰、或者多个截面同时成像清晰,在同一成像平面上形成沿轴向和横向同时偏移的具有衍射限像质的像面。系统的成像面个数、轴向间隔及像面间的横向偏移距离通过编程可实时控制调节。本发明可以实现无需轴向扫描、同时具有时间和空间分辨、灵活的三维成像系统,在生物医学成像、药理学、材料科学等领域具有重要的实际应用价值。

Figure 202011326887

The invention proposes a simultaneous imaging method of multiple object planes based on a digital multiplexing lens. The digital multiplexing lens programmed and modulated by the liquid crystal spatial light modulator, combined with the axial multi-object plane simultaneous imaging system composed of the front group of the fixed focal length infinite cylindrical microscope objective lens, can image any cross-section within a certain depth range of the sample Clear, or multiple cross-sections are imaged clearly at the same time, forming an image plane with diffraction-limited image quality that is simultaneously shifted in the axial and lateral directions on the same imaging plane. The number of imaging planes of the system, the axial spacing and the lateral offset distance between the imaging planes can be controlled and adjusted in real time through programming. The invention can realize a flexible three-dimensional imaging system that does not need axial scanning, has time and space resolution at the same time, and has important practical application value in the fields of biomedical imaging, pharmacology, material science and the like.

Figure 202011326887

Description

一种基于数字复用透镜的多物面同时成像方法A Multi-object Simultaneous Imaging Method Based on Digital Multiplex Lens

技术领域technical field

本发明属于光学成像领域,尤其涉及一种基于数字复用透镜的多物面同时成像方法。The invention belongs to the field of optical imaging, and in particular relates to a multi-object plane simultaneous imaging method based on a digital multiplexing lens.

背景技术Background technique

光学成像技术作为研究活细胞生命现象的直接可视化工具,实现灵活的时间和空间同时分辨成像一直以来是其面临的主要挑战。现有的成像技术,包括激光共聚焦扫描显微镜、双光子激光扫描显微镜、多焦点多光子显微成像技术,光学相干层析成像技术等,需要对样品实现横向和轴向扫描成像,限制了成像速度,丢失了时间分辨特性。Optical imaging technology, as a direct visualization tool for studying living phenomena in living cells, has always been the main challenge to achieve flexible time- and space-resolved imaging. Existing imaging techniques, including confocal laser scanning microscopy, two-photon laser scanning microscopy, multifocal multiphoton microscopy imaging, optical coherence tomography, etc., require lateral and axial scanning imaging of the sample, which limits imaging. speed, the time-resolved properties are lost.

基于变形光栅的多焦面并行成像方法,虽然能够实现多截面同时成像,但其截面成像能力主要集中在较低的衍射级次0级和±1级,并且各个级的衍射能量也不均匀。基于扭曲达曼光栅及多物面同时成像系统解决了截面成像能量分布不均匀性问题,并且实现了更多的截面成像,但涉及到复杂的加工工艺和检测技术手段,且灵活性、实时性不够,当改变截面个数和截面间隔时,需要重新加工扭曲达曼光栅。Although the multi-focal plane parallel imaging method based on deformed grating can realize simultaneous multi-section imaging, its cross-section imaging capability is mainly concentrated in the lower diffraction orders 0 and ±1, and the diffraction energy of each order is not uniform. Based on twisted Damman grating and multi-object simultaneous imaging system, the problem of energy distribution inhomogeneity in cross-sectional imaging is solved, and more cross-sectional imaging is realized, but it involves complex processing technology and detection technology, and is flexible and real-time. Not enough, when the number of sections and the section interval are changed, the twisted Damman grating needs to be reprocessed.

发明内容SUMMARY OF THE INVENTION

本发明将不同焦距的理想成像透镜和不同衍射角的闪耀光栅镶嵌在一起,由液晶空间光调制器(PLUTO)编程生成数字复用透镜。提出一种基于数字复用透镜的多物面同时成像方法。该成像方法,能够在样品的一定深度范围内,使任意截面成像清晰,或者多个截面同时成像清晰,具有相同衍射限成像质量,并且截面成像个数、截面成像间隔,截面横向偏移距离通过液晶空间光调制器(PLUTO)编程实时可控。本发明通过液晶空间光调制器实现数字复用透镜功能,从而分别实现所述多个理想成像透镜的焦距和多个闪耀光栅的衍射角。In the present invention, ideal imaging lenses with different focal lengths and blazed gratings with different diffraction angles are inlaid together, and a digital multiplexing lens is generated by programming a liquid crystal spatial light modulator (PLUTO). A multi-object simultaneous imaging method based on digital multiplex lens is proposed. The imaging method can make the imaging of any section clear within a certain depth range of the sample, or the imaging of multiple sections at the same time is clear, with the same diffraction limit imaging quality, and the number of section imaging, section imaging interval, section lateral offset distance passes through Liquid crystal spatial light modulator (PLUTO) programming is controllable in real time. The present invention realizes the digital multiplexing lens function through the liquid crystal spatial light modulator, thereby realizing the focal lengths of the plurality of ideal imaging lenses and the diffraction angles of the plurality of blazed gratings respectively.

本发明的技术方案为一种基于数字复用透镜的多物面同时成像方法,具体如下:The technical solution of the present invention is a multi-object simultaneous imaging method based on a digital multiplexing lens, and the details are as follows:

步骤1:构建数字复用透镜的多物面同时成像系统模型。Step 1: Build a multi-object simultaneous imaging system model of a digital multiplex lens.

步骤2:根据构建的多物面同时成像系统模型,按照理想光学系统的组合法则,通过组合成像系统的焦距、放大率、像方主平面位置、物方主平面位置、以及组合成像系统的物距变化量、组合成像系统物方主平面改变的距离随理想成像透镜的焦距变化,实现截面成像数量、截面成像间距的控制。Step 2: According to the constructed multi-object plane simultaneous imaging system model, according to the combination rule of the ideal optical system, by combining the focal length, magnification, image-side principal plane position, object-side principal plane position, and object of the combined imaging system The amount of distance change and the changed distance of the main plane of the combined imaging system on the object side change with the focal length of the ideal imaging lens, so as to realize the control of the number of cross-section imaging and the distance between cross-section imaging.

步骤3:建立理想成像透镜在液晶空间光调制器中的相位模型,建立闪耀光栅在液晶空间光调制器中的相位模型;实现每一成像截面具有不同的衍射角。Step 3: establish the phase model of the ideal imaging lens in the liquid crystal spatial light modulator, establish the phase model of the blazed grating in the liquid crystal spatial light modulator; realize that each imaging section has different diffraction angles.

步骤4:根据成像截面个数和成像间距的需要,设定理想成像透镜的数量、组合系统物方主平面改变的距离,通过设定的理想成像透镜的数量、组合系统物方主平面改变的距离结合步骤2计算理想成像透镜的焦距,通过步骤3计算闪耀光栅的衍射角,构建数字复用透镜在液晶空间光调制器的相位模型,通过液晶空间光调制器的编程调制实现数字复用透镜中理想成像透镜、闪耀光栅的设计功能。Step 4: According to the needs of the number of imaging sections and the imaging distance, set the number of ideal imaging lenses and the distance that the object-side principal plane of the combination system changes. The distance is combined with step 2 to calculate the focal length of the ideal imaging lens, the diffraction angle of the blazed grating is calculated through step 3, the phase model of the digital multiplexed lens in the liquid crystal spatial light modulator is constructed, and the digital multiplexed lens is realized through the programming modulation of the liquid crystal spatial light modulator Design functions of ideal imaging lenses and blazed gratings.

作为优选,步骤1所述数字复用透镜的多物面同时成像系统模型包括:Preferably, the multi-object simultaneous imaging system model of the digital multiplex lens in step 1 includes:

无限筒长显微物镜前组、数字复用透镜,所述数字复用透镜位于所述无限筒长显微物镜前组的后焦面;an infinite barrel length microscope objective lens front group and a digital multiplexing lens, the digital multiplexing lens is located on the back focal plane of the infinite barrel length microscope objective lens front group;

所述数字复用透镜由多个理想成像透镜、多个闪耀光栅构成;所述多个理想成像透镜为不同焦距的理想成像透镜,所述多个闪耀光栅为不同衍射角的闪耀光栅;所述理想成像透镜与所述闪耀光栅的数量相等,均为n,n≥1且n为整数;The digital multiplexing lens is composed of multiple ideal imaging lenses and multiple blazed gratings; the multiple ideal imaging lenses are ideal imaging lenses with different focal lengths, and the multiple blazed gratings are blazed gratings with different diffraction angles; the The number of the ideal imaging lens and the blazed grating is equal, and both are n, n≥1 and n is an integer;

所述理想成像透镜和所述闪耀光栅镶嵌在一起,由液晶空间光调制器编程调制生成;The ideal imaging lens and the blazed grating are inlaid together, and are generated by programming and modulation of a liquid crystal spatial light modulator;

作为优选,步骤2所述组合成像系统焦距为:Preferably, the focal length of the combined imaging system described in step 2 is:

f′zi=f'obj f' zi = f' obj

其中,f'obj为设无限筒长显微物镜前组的焦距,f′zi为组合成像系统的焦距,即在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距f′i如何变化,组合成像系统的焦距f′zi保持不变,n为理想成像透镜的数量;Among them, f' obj is the focal length of the front group of the infinite barrel length microscope objective, f' zi is the focal length of the combined imaging system, that is, the image-side focal plane of the front group of the infinite barrel length microscope objective with the focal length f' obj is superimposed For an ideal imaging lens with a focal length of f′ i i∈[1,n], no matter how the focal length f′ i of the ith ideal imaging lens changes, the focal length f′ zi of the combined imaging system remains unchanged, and n is the focal length of the ideal imaging lens. quantity;

步骤2所述组合成像系统放大率为:The magnification of the combined imaging system described in step 2 is:

βzi=ββ zi = β

其中,β为无限筒长显微物镜前组的放大倍率,βzi为组合成像系统的放大倍率,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距fi'如何变化,组合成像系统的放大倍率不随f′i的变化而发生变化,n为理想成像透镜的数量;Among them, β is the magnification of the front group of the infinite barrel length microscope objective, and β zi is the magnification of the combined imaging system, that is, the superimposed focal length of the image-side focal plane of the front group of the infinite barrel length microscope objective with the focal length f'obj The ideal imaging lens is f' i i∈[1,n], no matter how the focal length f i ' of the ith ideal imaging lens changes, the magnification of the combined imaging system does not change with the change of f' i , n is ideal the number of imaging lenses;

步骤2所述组合成像系统的像方主平面位置为:The position of the main plane of the image side of the combined imaging system described in step 2 is:

H'H′zi=0H'H' zi = 0

其中,H'为无限筒长显微物镜前组的像方主平面,H′zi为组合成像系统的像方主平面,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距f′i如何变化,H′zi始终与物镜前组的像方主平面H'重合,即组合成像系统的像平面保持不变,不随f′i的变化而发生变化,n为理想成像透镜的数量;Among them, H' is the principal plane of the image side of the front group of the infinite barrel length microscope objective, H' zi is the principal plane of the image side of the combined imaging system, that is, the image side principal plane of the front group of the infinite barrel length microscope objective with the focal length f' obj The image-side focal plane superimposes an ideal imaging lens with a focal length of f′ i i∈[1,n]. No matter how the focal length f′ i of the i-th ideal imaging lens changes, H′ zi is always connected to the image-side principal plane of the front group of the objective lens. H' coincides, that is, the image plane of the combined imaging system remains unchanged and does not change with the change of f' i , and n is the number of ideal imaging lenses;

步骤2所述组合成像系统的物方主平面位置为:The position of the object-side principal plane of the combined imaging system described in step 2 is:

Figure BDA0002794031220000031
Figure BDA0002794031220000031

H为无限筒长显微物镜前组的物方主平面,Hzi为组合成像系统的物方主平面,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,组合成像系统的物方主平面随f′i的变化而发生变化,n为理想成像透镜的数量;H is the object-side principal plane of the front group of the infinite barrel length microscope objective, H zi is the object side principal plane of the combined imaging system, that is, the image side focal plane of the front group of the infinite barrel length microscope objective with the focal length f'obj Superimpose an ideal imaging lens with a focal length of f′ i i∈[1,n], the object-side principal plane of the combined imaging system changes with the change of f′ i , and n is the number of ideal imaging lenses;

步骤2所述组合成像系统的物距位置随f′i变化为:The object distance position of the combined imaging system described in step 2 changes with f′ i as follows:

Figure BDA0002794031220000032
Figure BDA0002794031220000032

其中,lzi为组合成像系统的物距,l'为无限筒长显微物镜前组的像距,该公式表明组合成像系统物方变化的距离等于组合系统物方主平面改变的距离,n为理想成像透镜的数量;Among them, l zi is the object distance of the combined imaging system, l' is the image distance of the front group of the infinite tube long microscope objective, this formula shows that the distance of the combined imaging system object side change is equal to the combined system object side principal plane change distance, n is the number of ideal imaging lenses;

步骤2所述组合系统物方主平面改变的距离为:The distance of the change of the main plane of the composite system in the step 2 is:

Figure BDA0002794031220000033
Figure BDA0002794031220000033

其中,Δzi为加载第i个成像透镜后组合系统物方主平面改变的距离,n为理想成像透镜的数量;Among them, Δzi is the distance that the object-side principal plane of the combined system changes after the ith imaging lens is loaded, and n is the number of ideal imaging lenses;

通过步骤2所述组合系统物方主平面改变的距离计算步骤2所述理想成像透镜的焦距;Calculate the focal length of the ideal imaging lens described in step 2 by the distance changed by the object-side principal plane of the combination system described in step 2;

通过改变理想成像透镜的数量即n和fi',可以用于控制截面成像数量及截面成像间距,i∈[1,n]By changing the number of ideal imaging lenses, namely n and f i ', it can be used to control the number of cross-section imaging and the distance between cross-section imaging, i∈[1,n]

通过改变每一成像截面在x,y方向对应的衍射角θx-i,θy-i,可以控制截面的横向偏移距离,在x,y方向对应的衍射角θx-i,θy-i通过闪耀光栅实现;By changing the diffraction angles θ xi , θ yi corresponding to each imaging section in the x, y directions, the lateral offset distance of the section can be controlled, and the corresponding diffraction angles θ xi , θ yi in the x, y directions are realized by blazed gratings;

作为优选,步骤3所述建立理想成像透镜在液晶空间光调制器中的相位模型为:Preferably, the phase model of the ideal imaging lens in the liquid crystal spatial light modulator described in step 3 is:

数字复用透镜是一种由液晶空间光调制器生成的相位图,根据傅里叶光学理论,第i个理想成像透镜对成像波前的相位调制量为:The digital multiplexing lens is a phase map generated by a liquid crystal spatial light modulator. According to Fourier optics theory, the phase modulation amount of the i-th ideal imaging lens to the imaging wavefront is:

Figure BDA0002794031220000041
Figure BDA0002794031220000041

式中,λ为入射波长,x,y为以透镜中心为原点的坐标,fi'为第i个理想成像透镜的焦距,n为理想成像透镜的数量;where λ is the incident wavelength, x and y are the coordinates with the center of the lens as the origin, f i ' is the focal length of the ith ideal imaging lens, and n is the number of ideal imaging lenses;

设液晶空间光调制器的像素分辨率为M*N,像素中心间距为a,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量;Let the pixel resolution of the liquid crystal spatial light modulator be M*N, the pixel center spacing be a, M is the number of rows of liquid crystal spatial light modulator pixels, and N is the number of columns of liquid crystal spatial light modulator pixels;

以液晶空间光调制器中心为坐标原点,则焦距为fi'的理想成像透镜在液晶空间光调制器中的相位模型可表示为:Taking the center of the liquid crystal spatial light modulator as the coordinate origin, the phase model of the ideal imaging lens with focal length f i ' in the liquid crystal spatial light modulator can be expressed as:

Figure BDA0002794031220000042
Figure BDA0002794031220000042

Figure BDA0002794031220000043
Figure BDA0002794031220000043

式中,mod表示以2π取余操作,k为第一系数,l为第二系数,a为像素中心间距,fi'为第i个理想成像透镜的焦距,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,n为理想成像透镜的数量;In the formula, mod represents the operation of taking the remainder of 2π, k is the first coefficient, l is the second coefficient, a is the pixel center distance, f i ' is the focal length of the ith ideal imaging lens, and M is the liquid crystal spatial light modulator. The number of rows of pixels, N is the number of columns of liquid crystal spatial light modulator pixels, and n is the number of ideal imaging lenses;

步骤3所述建立闪耀光栅在液晶空间光调制器中的相位模型为:The phase model of the blazed grating in the liquid crystal spatial light modulator established in step 3 is:

为了实现不同的轴向截面成像在像平面的不同区域而互不重叠,不同焦距的理想成像透镜其出射光束应具有不同的衍射角;In order to realize the imaging of different axial sections in different regions of the image plane without overlapping each other, the outgoing beams of ideal imaging lenses with different focal lengths should have different diffraction angles;

利用液晶空间光调制器产生闪耀光栅来实现,针对焦距为fi'的第i个理想成像透镜的出射光束,其加载的对应的第i个闪耀光栅的相位模型分布表示为:It is realized by using the liquid crystal spatial light modulator to generate a blazed grating. For the outgoing beam of the i-th ideal imaging lens with a focal length f i ', the phase model distribution of the corresponding i-th blazed grating loaded is expressed as:

Figure BDA0002794031220000051
Figure BDA0002794031220000051

Figure BDA0002794031220000052
Figure BDA0002794031220000052

式中,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,k为第一系数,l为第二系数,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,n为理想成像透镜的数量;In the formula, T xi is the grating period corresponding to the i-th ideal imaging lens in the x-direction in pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y-direction in pixels, and k is the first coefficient, l is the second coefficient, 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, and n is the number of ideal imaging lenses;

步骤3所述每一成像截面具有不同的衍射角为:Each imaging section described in step 3 has a different diffraction angle:

加载第i个闪耀光栅后,对应fi'的成像截面光束在x,y方向的衍射角分别为:After loading the i-th blazed grating, the diffraction angles of the imaging section beam corresponding to f i ' in the x and y directions are:

Figure BDA0002794031220000053
Figure BDA0002794031220000053

其中,a为像素中心间距,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,n为理想成像透镜的数量;Among them, a is the pixel center spacing, T xi is the grating period corresponding to the i-th ideal imaging lens in the x-direction in pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y-direction in pixels, n is the number of ideal imaging lenses;

其衍射角大小分别与Tx-i和Ty-i有关,衍射角的方向取决于k,l的符号;The size of the diffraction angle is related to T xi and T yi respectively, and the direction of the diffraction angle depends on the sign of k, l;

作为优选,步骤4所述数字复用透镜在液晶空间光调制器的相位模型为:Preferably, the phase model of the digital multiplexing lens in the liquid crystal spatial light modulator in step 4 is:

Figure BDA0002794031220000054
Figure BDA0002794031220000054

Figure BDA0002794031220000055
Figure BDA0002794031220000055

其中,a为像素中心间距,fi'为第i个理想成像透镜的焦距,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,k为第一系数,l为第二系数。where a is the pixel center spacing, f i ' 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, and T xi is The x-direction is the grating period corresponding to the i-th ideal imaging lens in pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y-direction in pixels, k is the first coefficient, and l is the second coefficient.

本发明将不同焦距的理想成像透镜和不同衍射角的闪耀光栅镶嵌在一起,由液晶空间光调制器(PLUTO)编程生成数字复用透镜。提出一种基于数字复用透镜的多物面同时成像方法。该成像方法,能够在样品的一定深度范围内,使任意截面成像清晰,或者多个截面同时成像清晰,具有相同衍射限成像质量,并且截面成像个数、截面成像间隔,截面横向偏移距离通过液晶空间光调制器(PLUTO)编程实时可控。本发明通过液晶空间光调制器实现数字复用透镜功能,从而分别实现所述多个理想成像透镜的焦距和多个闪耀光栅的衍射角。In the present invention, ideal imaging lenses with different focal lengths and blazed gratings with different diffraction angles are inlaid together, and a digital multiplexing lens is generated by programming a liquid crystal spatial light modulator (PLUTO). A multi-object simultaneous imaging method based on digital multiplex lens is proposed. The imaging method can make the imaging of any section clear within a certain depth range of the sample, or the imaging of multiple sections at the same time is clear, with the same diffraction limit imaging quality, and the number of section imaging, section imaging interval, section lateral offset distance passes through Liquid crystal spatial light modulator (PLUTO) programming is controllable in real time. The present invention realizes the digital multiplexing lens function through the liquid crystal spatial light modulator, thereby realizing the focal lengths of the plurality of ideal imaging lenses and the diffraction angles of the plurality of blazed gratings respectively.

本发明可以实现无需轴向扫描、同时具有时间和空间分辨、灵活的三维成像系统,在生物医学成像、药理学、材料科学等领域具有重要的实际应用价值。The invention can realize a flexible three-dimensional imaging system that does not need axial scanning, has time and space resolution at the same time, and has important practical application value in the fields of biomedical imaging, pharmacology, material science and the like.

附图说明Description of drawings

图1:数字复用透镜的多物面同时成像系统。Figure 1: Multi-object simultaneous imaging system with digital multiplex lens.

图2:组合成像系统焦距示意图。Figure 2: Schematic diagram of the focal length of the combined imaging system.

图3:二维矩阵示意图。Figure 3: Schematic diagram of a two-dimensional matrix.

图4:标定的相位灰度图。Figure 4: The calibrated phase grayscale map.

图5:本发明方法流程图。Figure 5: Flow chart of the method of the present invention.

具体实施方式Detailed ways

为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。In order to facilitate the understanding and implementation of the present invention by those of ordinary skill in the art, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the embodiments described herein are only used to illustrate and explain the present invention, but not to limit it. this invention.

本发明实施方式汇总,液晶空间光调制器采用德国HOLOEYE公司生产的PLUTO-VIS型液晶空间光调制器,像面尺寸15.36mm×8.64mm,分辨率M×N为1920×1080,像元尺寸为8μm,,像素之间的间隔为1μm,则单个像素的大小a为9μmSummary of the embodiments of the present invention, the liquid crystal spatial light modulator adopts the PLUTO-VIS liquid crystal spatial light modulator produced by HOLOEYE company in Germany, the image size is 15.36mm×8.64mm, the resolution M×N is 1920×1080, and the pixel size is 8μm, and the interval between pixels is 1μm, the size a of a single pixel is 9μm

下面结合图1至图5描述本发明的具体实施方式为:Below in conjunction with Fig. 1 to Fig. 5, describe the specific embodiment of the present invention as follows:

一种基于数字复用透镜的多物面同时成像方法,具体如下:A multi-object simultaneous imaging method based on a digital multiplexing lens, the details are as follows:

步骤1:构建数字复用透镜的多物面同时成像系统模型,如图1所示。Step 1: Build a multi-object simultaneous imaging system model of a digital multiplex lens, as shown in Figure 1.

步骤1所述数字复用透镜的多物面同时成像系统模型包括:The multi-object simultaneous imaging system model of the digital multiplex lens in step 1 includes:

无限筒长显微物镜前组、数字复用透镜,所述数字复用透镜位于所述无限筒长显微物镜前组的后焦面;an infinite barrel length microscope objective lens front group and a digital multiplexing lens, the digital multiplexing lens is located on the back focal plane of the infinite barrel length microscope objective lens front group;

所述数字复用透镜的多物面同时成像系统采用氦氖激光作照明光源,波长λ=0.6328μ,选择无限筒长物镜前组焦距f′obj=25mm,数值孔径NA=0.25,倍率β=10×The multi-object plane simultaneous imaging system of the digital multiplexing lens adopts helium-neon laser as the illumination light source, the wavelength λ=0.6328μ, the focal length of the front group of the infinite barrel long objective lens is selected f' obj =25mm, the numerical aperture NA=0.25, the magnification β= 10 × ;

所述数字复用透镜由多个理想成像透镜、多个闪耀光栅构成;所述多个理想成像透镜为不同焦距的理想成像透镜,所述多个闪耀光栅为不同衍射角的闪耀光栅;所述理想成像透镜与所述闪耀光栅的数量相等,均为n=9,n≥1且n为整数;The digital multiplexing lens is composed of multiple ideal imaging lenses and multiple blazed gratings; the multiple ideal imaging lenses are ideal imaging lenses with different focal lengths, and the multiple blazed gratings are blazed gratings with different diffraction angles; the The number of the ideal imaging lens and the blazed grating is equal, and both are n=9, n≥1 and n is an integer;

所述理想成像透镜和所述闪耀光栅镶嵌在一起,由液晶空间光调制器编程调制生成;The ideal imaging lens and the blazed grating are inlaid together, and are generated by programming and modulation of a liquid crystal spatial light modulator;

步骤2:根据构建的多物面同时成像系统模型,按照理想光学系统的组合法则,通过组合成像系统的焦距、放大率、像方主平面位置、物方主平面位置、以及组合成像系统的物距变化量、组合成像系统物方主平面改变的距离随理想成像透镜的焦距变化,实现截面成像数量、截面成像间距的控制。Step 2: According to the constructed multi-object plane simultaneous imaging system model, according to the combination rule of the ideal optical system, by combining the focal length, magnification, image-side principal plane position, object-side principal plane position, and object of the combined imaging system The amount of distance change and the changed distance of the main plane of the combined imaging system on the object side change with the focal length of the ideal imaging lens, so as to realize the control of the number of cross-section imaging and the distance between cross-section imaging.

如图2示意所示;As shown schematically in Figure 2;

步骤2所述组合成像系统焦距为:The focal length of the combined imaging system described in step 2 is:

f′zi=f'obj f' zi = f' obj

其中,f'obj为设无限筒长显微物镜前组的焦距,f′zi为组合成像系统的焦距,即在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距f′i如何变化,组合成像系统的焦距f′zi保持不变,n为理想成像透镜的数量;Among them, f' obj is the focal length of the front group of the infinite barrel length microscope objective, f' zi is the focal length of the combined imaging system, that is, the image-side focal plane of the front group of the infinite barrel length microscope objective with the focal length f' obj is superimposed For an ideal imaging lens with a focal length of f′ i i∈[1,n], no matter how the focal length f′ i of the ith ideal imaging lens changes, the focal length f′ zi of the combined imaging system remains unchanged, and n is the focal length of the ideal imaging lens. quantity;

所述第i个理想成像透镜的焦距f′i根据截面成像间距的要求和发明内容中步骤2所述组合系统物方主平面改变的距离确定;The focal length f′ i of the ith ideal imaging lens is determined according to the requirements of the cross-sectional imaging distance and the distance changed by the object-side principal plane of the combined system in step 2 of the Summary of the Invention;

其中,对P0截面,Δz0=0,f0'=0,而P1至P8截面,设定要求Δz1=Δz2=…=Δz8=0.1mm,如图1所示;Among them, for the P 0 section, Δz 0 =0, f 0 '=0, and for the P 1 to P 8 sections, the setting requirements are Δz 1 =Δz 2 =...=Δz 8 =0.1mm, as shown in Figure 1;

f′1=6250mm,f′2=3125mm,f′3=2083.33mm,f′4=1562.5mmf' 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.25mmf' 5 =1250mm, f' 6 =1041.67mm, f' 7 =892.85mm, f' 8 =781.25mm

步骤2所述组合成像系统放大率为:The magnification of the combined imaging system described in step 2 is:

βzi=ββ zi = β

其中,β为无限筒长显微物镜前组的放大倍率,βzi为组合成像系统的放大倍率,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距f′i如何变化,组合成像系统的放大倍率不随f′i的变化而发生变化,n为理想成像透镜的数量;Among them, β is the magnification of the front group of the infinite barrel length microscope objective, and β zi is the magnification of the combined imaging system, that is, the superimposed focal length of the image-side focal plane of the front group of the infinite barrel length microscope objective with the focal length f'obj For an ideal imaging lens of f′ i i∈[1,n], no matter how the focal length f′ i of the ith ideal imaging lens changes, the magnification of the combined imaging system does not change with the change of f′ i , and n is ideal the number of imaging lenses;

步骤2所述组合成像系统的像方主平面位置为:The position of the main plane of the image side of the combined imaging system described in step 2 is:

H'H′zi=0H'H' zi = 0

其中,H'为无限筒长显微物镜前组的像方主平面,H′zi为组合成像系统的像方主平面,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,不管第i个理想成像透镜的焦距f′i如何变化,H′zi始终与物镜前组的像方主平面H'重合,即组合成像系统的像平面保持不变,不随f′i的变化而发生变化,n为理想成像透镜的数量;Among them, H' is the principal plane of the image side of the front group of the infinite barrel length microscope objective, H' zi is the principal plane of the image side of the combined imaging system, that is, the image side principal plane of the front group of the infinite barrel length microscope objective with the focal length f' obj The image-side focal plane superimposes an ideal imaging lens with a focal length of f′ i i∈[1,n]. No matter how the focal length f′ i of the i-th ideal imaging lens changes, H′ zi is always connected to the image-side principal plane of the front group of the objective lens. H' coincides, that is, the image plane of the combined imaging system remains unchanged and does not change with the change of f' i , and n is the number of ideal imaging lenses;

步骤2所述组合成像系统的物方主平面位置为:The position of the object-side principal plane of the combined imaging system described in step 2 is:

Figure BDA0002794031220000081
Figure BDA0002794031220000081

H为无限筒长显微物镜前组的物方主平面,Hzi为组合成像系统的物方主平面,即,在焦距为f'obj的无限筒长显微物镜前组的像方焦面叠加焦距为f′i i∈[1,n]的理想成像透镜,组合成像系统的物方主平面随f′i的变化而发生变化,n为理想成像透镜的数量;H is the object-side principal plane of the front group of the infinite barrel length microscope objective, H zi is the object side principal plane of the combined imaging system, that is, the image side focal plane of the front group of the infinite barrel length microscope objective with the focal length f'obj Superimpose an ideal imaging lens with a focal length of f′ i i∈[1,n], the object-side principal plane of the combined imaging system changes with the change of f′ i , and n is the number of ideal imaging lenses;

步骤2所述组合成像系统的物距位置随f′i变化为:The object distance position of the combined imaging system described in step 2 changes with f′ i as follows:

Figure BDA0002794031220000082
Figure BDA0002794031220000082

其中,为lzi个组合成像系统的物距,l'为无限筒长显微物镜前组的像距,该公式表明组合成像系统物方变化的距离等于组合系统物方主平面改变的距离,n为理想成像透镜的数量;Among them, is the object distance of l zi combined imaging systems, and l' is the image distance of the front group of the infinite tube long microscope objective. n is the number of ideal imaging lenses;

步骤2所述组合系统物方主平面改变的距离为:The distance of the change of the main plane of the composite system in the step 2 is:

Figure BDA0002794031220000091
Figure BDA0002794031220000091

其中,Δzi为加载第i个成像透镜后组合系统物方主平面改变的距离,n为理想成像透镜的数量;Among them, Δzi is the distance that the object-side principal plane of the combined system changes after the ith imaging lens is loaded, and n is the number of ideal imaging lenses;

通过步骤2所述组合系统物方主平面改变的距离计算步骤2所述理想成像透镜的焦距;Calculate the focal length of the ideal imaging lens described in step 2 by the distance changed by the object-side principal plane of the combination system described in step 2;

通过改变理想成像透镜的数量即n和fi',可以用于控制截面成像数量及截面成像间距,i∈[1,n]By changing the number of ideal imaging lenses, namely n and f i ', it can be used to control the number of cross-section imaging and the distance between cross-section imaging, i∈[1,n]

通过改变每一成像截面在x,y方向对应的衍射角θx-i,θy-i,可以控制截面的横向偏移距离,在x,y方向对应的衍射角θx-i,θy-i通过闪耀光栅实现;By changing the diffraction angles θ xi , θ yi corresponding to each imaging section in the x, y directions, the lateral offset distance of the section can be controlled, and the corresponding diffraction angles θ xi , θ yi in the x, y directions are realized by blazed gratings;

步骤3:建立理想成像透镜在液晶空间光调制器中的相位模型,建立闪耀光栅在液晶空间光调制器中的相位模型;实现每一成像截面具有不同的衍射角。Step 3: establish the phase model of the ideal imaging lens in the liquid crystal spatial light modulator, establish the phase model of the blazed grating in the liquid crystal spatial light modulator; realize that each imaging section has different diffraction angles.

步骤3所述建立理想成像透镜在液晶空间光调制器中的相位模型为:The phase model of the ideal imaging lens in the liquid crystal spatial light modulator described in step 3 is:

数字复用透镜是一种由液晶空间光调制器生成的相位图,根据傅里叶光学理论,第i个理想成像透镜对成像波前的相位调制量为:The digital multiplexing lens is a phase map generated by a liquid crystal spatial light modulator. According to Fourier optics theory, the phase modulation amount of the i-th ideal imaging lens to the imaging wavefront is:

Figure BDA0002794031220000092
Figure BDA0002794031220000092

式中,λ为入射波长,x,y为以透镜中心为原点的坐标,fi'为第i个理想成像透镜的焦距,n为理想成像透镜的数量;where λ is the incident wavelength, x and y are the coordinates with the center of the lens as the origin, f i ' is the focal length of the ith ideal imaging lens, and n is the number of ideal imaging lenses;

设液晶空间光调制器的像素分辨率为M*N,像素中心间距为a,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量;Let the pixel resolution of the liquid crystal spatial light modulator be M*N, the pixel center spacing be a, M is the number of rows of liquid crystal spatial light modulator pixels, and N is the number of columns of liquid crystal spatial light modulator pixels;

以液晶空间光调制器中心为坐标原点,则焦距为fi'的理想成像透镜在液晶空间光调制器中的相位模型可表示为:Taking the center of the liquid crystal spatial light modulator as the coordinate origin, the phase model of the ideal imaging lens with focal length f i ' in the liquid crystal spatial light modulator can be expressed as:

Figure BDA0002794031220000093
Figure BDA0002794031220000093

Figure BDA0002794031220000101
Figure BDA0002794031220000101

式中,mod表示以2π取余操作,k为第一系数,l为第二系数,a为像素中心间距,fi'为第i个理想成像透镜的焦距,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,n为理想成像透镜的数量;In the formula, mod represents the operation of taking the remainder of 2π, k is the first coefficient, l is the second coefficient, a is the pixel center distance, f i ' is the focal length of the ith ideal imaging lens, and M is the liquid crystal spatial light modulator. The number of rows of pixels, N is the number of columns of liquid crystal spatial light modulator pixels, and n is the number of ideal imaging lenses;

步骤3所述建立闪耀光栅在液晶空间光调制器中的相位模型为:The phase model of the blazed grating in the liquid crystal spatial light modulator established in step 3 is:

为了实现不同的轴向截面成像在像平面的不同区域而互不重叠,不同焦距的理想成像透镜其出射光束应具有不同的衍射角;In order to realize the imaging of different axial sections in different regions of the image plane without overlapping each other, the outgoing beams of ideal imaging lenses with different focal lengths should have different diffraction angles;

利用液晶空间光调制器产生闪耀光栅来实现,针对焦距为fi'的第i个理想成像透镜的出射光束,其加载的对应的第i个闪耀光栅的相位模型分布表示为:It is realized by using the liquid crystal spatial light modulator to generate a blazed grating. For the outgoing beam of the i-th ideal imaging lens with a focal length f i ', the phase model distribution of the corresponding i-th blazed grating loaded is expressed as:

Figure BDA0002794031220000102
Figure BDA0002794031220000102

Figure BDA0002794031220000103
Figure BDA0002794031220000103

式中,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,k为第一系数,l为第二系数,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,n为理想成像透镜的数量;In the formula, T xi is the grating period corresponding to the i-th ideal imaging lens in the x-direction in pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y-direction in pixels, and k is the first coefficient, l is the second coefficient, 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, and n is the number of ideal imaging lenses;

步骤3所述每一成像截面具有不同的衍射角为:Each imaging section described in step 3 has a different diffraction angle:

加载第i个闪耀光栅后,对应fi'的成像截面光束在x,y方向的衍射角分别为:After loading the i-th blazed grating, the diffraction angles of the imaging section beam corresponding to f i ' in the x and y directions are:

Figure BDA0002794031220000104
Figure BDA0002794031220000104

其中,a为像素中心间距,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,n为理想成像透镜的数量;Among them, a is the pixel center spacing, T xi is the grating period corresponding to the i-th ideal imaging lens in the x-direction in pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y-direction in pixels, n is the number of ideal imaging lenses;

其衍射角大小分别与Tx-i和Ty-i有关,衍射角的方向取决于k,l的符号;The size of the diffraction angle is related to T xi and T yi respectively, and the direction of the diffraction angle depends on the sign of k, l;

步骤4:根据成像截面个数和成像间距的需要,设定理想成像透镜的数量、组合系统物方主平面改变的距离,通过设定的理想成像透镜的数量、组合系统物方主平面改变的距离结合步骤2计算理想成像透镜的焦距,通过步骤3计算闪耀光栅的衍射角,构建数字复用透镜在液晶空间光调制器的相位模型,通过液晶空间光调制器的编程调制实现数字复用透镜中理想成像透镜、闪耀光栅的设计功能。Step 4: According to the needs of the number of imaging sections and the imaging distance, set the number of ideal imaging lenses and the distance that the object-side principal plane of the combination system changes. The distance is combined with step 2 to calculate the focal length of the ideal imaging lens, the diffraction angle of the blazed grating is calculated through step 3, the phase model of the digital multiplexed lens in the liquid crystal spatial light modulator is constructed, and the digital multiplexed lens is realized through the programming modulation of the liquid crystal spatial light modulator Design functions of ideal imaging lenses and blazed gratings.

步骤4所述数字复用透镜在液晶空间光调制器的相位模型为:The phase model of the digital multiplexing lens in the liquid crystal spatial light modulator in step 4 is:

Figure BDA0002794031220000111
Figure BDA0002794031220000111

Figure BDA0002794031220000112
Figure BDA0002794031220000112

其中,a为像素中心间距,fi为第i个理想成像透镜的焦距,M为液晶空间光调制器像素的行的数量,N为液晶空间光调制器像素的列的数量,Tx-i为x方向以像素为单位第i个理想成像透镜对应的光栅周期,Ty-i为y方向以像素为单位第i个理想成像透镜对应的的光栅周期,k为第一系数,l为第二系数。where a is the pixel center spacing, f i 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, and T xi is x The direction is the grating period corresponding to the i-th ideal imaging lens in the unit of pixels, T yi is the grating period corresponding to the i-th ideal imaging lens in the y direction in the unit of pixel, k is the first coefficient, and l is the second coefficient.

所述数字复用透镜的相位分布如表1所示:The phase distribution of the digital multiplexing lens is shown in Table 1:

表1数字复用透镜相位分布Table 1 Phase distribution of digital multiplexing lens

Figure BDA0002794031220000121
Figure BDA0002794031220000121

步骤5:调用Matlab生成二维矩阵,将9个相位分布随机给二维矩阵赋初值,如图3所示。Step 5: Call Matlab to generate a two-dimensional matrix, and assign the initial values of the nine phase distributions to the two-dimensional matrix randomly, as shown in Figure 3.

步骤6:如图4所示为标定的相位灰度图,根据表1的相位表达式计算矩阵中每一像素点对应的相位调制量,生成灰度矩阵。Step 6: Figure 4 shows the calibrated phase grayscale map, calculates the phase modulation amount corresponding to each pixel in the matrix according to the phase expression in Table 1, and generates a grayscale matrix.

PLUTO的工作过程是通过改变寻址电压控制每个像素的相位调制量,PLUTO在出厂时制造商已经将驱动电压映射为计算机显示的灰度值,厂家给出了(灰度-相位)客户查寻表,但具体到每一个器件需要对器件进行重新标定,得到如图4所示的灰度-相位标定曲线,将该曲线对应的灰度加载在PLUTO的驱动软件上,实现要求的相位调制。The working process of PLUTO is to control the phase modulation amount of each pixel by changing the addressing voltage. When PLUTO leaves the factory, the manufacturer has mapped the driving voltage to the gray value displayed by the computer. The manufacturer gives the (gray-phase) customer search However, specific to each device, the device needs to be re-calibrated to obtain the grayscale-phase calibration curve as shown in Figure 4, and the grayscale corresponding to the curve is loaded on the PLUTO driver software to achieve the required phase modulation.

步骤7:利用Imwrite函数生成灰度图,加载在空间光调制器上。Step 7: Use the Imwrite function to generate a grayscale image and load it on the spatial light modulator.

应当理解的是,本说明书未详细阐述的部分均属于现有技术。It should be understood that the parts not described in detail in this specification belong to the prior art.

应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。It should be understood that the above description of the preferred embodiments is relatively detailed, and therefore should not be considered as a limitation on the protection scope of the patent of the present invention. In the case of the protection scope, substitutions or deformations can also be made, which all fall within the protection scope of the present invention, and the claimed protection scope of the present invention shall be subject to 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 深圳大学 Calculation method of actual focal length of variable focus 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 深圳大学 Calculation method of actual focal length of variable focus 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
Zuo et al. Noninterferometric single-shot quantitative phase microscopy
Baddeley et al. Three-dimensional sub-100 nm super-resolution imaging of biological samples using a phase ramp in the objective pupil
Tao et al. Live imaging using adaptive optics with fluorescent protein guide-stars
Dean et al. Uniform and scalable light-sheets generated by extended focusing
Girkin et al. Adaptive optics for deeper imaging of biological samples
Schrader et al. 4Pi-confocal imaging in fixed biological specimens
Tsvirkun et al. Widefield lensless endoscopy with a multicore fiber
Pan et al. Single-camera microscopic stereo digital image correlation using a diffraction grating
Li et al. Conjugate adaptive optics in widefield microscopy with an extended-source wavefront sensor
Yang et al. Influence of the photonic nanojet of microspheres on microsphere imaging
Hubert et al. Adaptive optics light-sheet microscopy based on direct wavefront sensing without any guide star
Paudel et al. Axial range of conjugate adaptive optics in two-photon microscopy
Philipp et al. Diffraction-limited axial scanning in thick biological tissue with an aberration-correcting adaptive lens
Larsson et al. Microscopic 3-D displacement field measurements using digital speckle photography
WO2018172161A1 (en) Microscope and method for imaging an object
US20160202462A1 (en) System and method of edge-illumination microscopy
CN112379516B (en) Multi-object-plane simultaneous imaging method based on digital multiplexing lens
Luo et al. Area scanning method for 3D surface profilometry based on an adaptive confocal microscope
Ren et al. Aberration-free large-area stitch-free 3D nano-printing based on binary holography
US9930241B2 (en) Image processing apparatus, image processing program, and image processing method
Turcotte et al. Deconvolution for multimode fiber imaging: modeling of spatially variant PSF
Vladimirov et al. Dual-view light-sheet imaging through a tilted glass interface using a deformable mirror
Chen et al. Efficient light focusing through an MMF based on two-step phase shifting and parallel phase compensating
Shao et al. Addressable multiregional and multifocal multiphoton microscopy based on a spatial light modulator
Li et al. Design and fabrication of a freeform prism array for 3D microscopy

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