CN106990694B - A non-iterative phase recovery device and method under partially coherent light illumination - Google Patents
A non-iterative phase recovery device and method under partially coherent light illumination Download PDFInfo
- Publication number
- CN106990694B CN106990694B CN201710197745.5A CN201710197745A CN106990694B CN 106990694 B CN106990694 B CN 106990694B CN 201710197745 A CN201710197745 A CN 201710197745A CN 106990694 B CN106990694 B CN 106990694B
- Authority
- CN
- China
- Prior art keywords
- phase
- light
- coherent light
- measured
- partially coherent
- 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
Links
- 230000001427 coherent effect Effects 0.000 title claims abstract description 90
- 238000011084 recovery Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005286 illumination Methods 0.000 title claims abstract description 27
- 239000005337 ground glass Substances 0.000 claims abstract description 12
- 238000012216 screening Methods 0.000 claims description 37
- 238000012545 processing Methods 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 claims description 6
- 238000007493 shaping process Methods 0.000 claims description 4
- 239000000758 substrate Substances 0.000 claims description 4
- 238000010329 laser etching Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims 7
- 238000005530 etching Methods 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 15
- 230000008878 coupling Effects 0.000 abstract 1
- 238000010168 coupling process Methods 0.000 abstract 1
- 238000005859 coupling reaction Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 19
- 230000008569 process Effects 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 5
- 239000005338 frosted glass Substances 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000001093 holography Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003574 free electron Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002086 nanomaterial Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000002207 retinal effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/10—Processes or apparatus for producing holograms using modulated reference beam
- G03H1/12—Spatial modulation, e.g. ghost imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
技术领域technical field
本发明涉及光学技术领域,尤其涉及一种部分相干光照明下的非迭代相位恢复装置及方法。The present invention relates to the field of optical technology, and in particular, to a non-iterative phase recovery device and method under partially coherent light illumination.
背景技术Background technique
对于一个包含振幅和相位信息的未知物体,可以通过电荷耦合元件对物体振幅信息进行直接观测,而相位信息并不能被直接获得,因此如何从强度信息中获取未知物体的相位信息成为了人们研究的一个重要课题。从强度信息中获取物体相位信息的技术叫做波前探测技术或者相位恢复,通过衍射或者干涉的方法获取物体相位信息以实现二维、三维成像的方法叫做相干衍射成像,该技术广泛运用于图像信息处理、微纳光学、自适应光学、材料科学以及量子层析等领域。随着相干衍射成像技术的飞速发展,分辨率已达到纳米量级,恢复装置也愈加智能、实时化。For an unknown object containing amplitude and phase information, the amplitude information of the object can be directly observed through the charge-coupled element, but the phase information cannot be directly obtained. Therefore, how to obtain the phase information of the unknown object from the intensity information has become a research topic. an important subject. The technology of obtaining the phase information of the object from the intensity information is called wavefront detection technology or phase recovery. The method of obtaining the phase information of the object through diffraction or interference to realize two-dimensional and three-dimensional imaging is called coherent diffraction imaging, which is widely used in image information. processing, micro-nano optics, adaptive optics, materials science, and quantum tomography. With the rapid development of coherent diffraction imaging technology, the resolution has reached the nanometer level, and the recovery device has become more intelligent and real-time.
在研究波前探测技术和相位恢复时,大多都假设照明光源为完全相干光,但实际应用中,例如高分辨波前探测,往往使用x射线或者电子束作为光源进行相位恢复,这些都不是完全相干光。此外,当一束完全相干光经过介质传输以后,其空间相干性也会降低,在这些情况下,仍然将其作为完全相干光进行处理,会存在一些问题。When studying wavefront detection technology and phase recovery, most of them assume that the illumination light source is completely coherent light, but in practical applications, such as high-resolution wavefront detection, X-rays or electron beams are often used as light sources for phase recovery, which are not completely coherent light. In addition, when a beam of fully coherent light is transmitted through the medium, its spatial coherence will also be reduced. In these cases, it is still treated as fully coherent light, there will be some problems.
实现相位恢复的方法有很多种,最早是在1952年,David Sayre提出利用香农定理通过测量更高密度的光强来实现相位恢复。至今为止,人们已经研究出了一系列相位恢复的方法,例如哈特曼波前传感技术、全息干涉技术、计算相位恢复技术和叠层技术等。There are many ways to achieve phase recovery. As early as 1952, David Sayre proposed to use Shannon's theorem to achieve phase recovery by measuring higher-density light intensity. So far, a series of phase recovery methods have been developed, such as Hartmann wavefront sensing technology, holographic interference technology, computational phase recovery technology and stacking technology.
哈特曼波前传感技术主要是通过测量波前斜率来恢复相位信息(Platt B C,Shack R.History and principles of Shack-Hartmann wavefront sensing[J].Journalof Refractive Surgery,2001,17(5):S573-S577.),夏克哈特曼传感器基于此技术,由微透镜阵列和电荷耦合元件构成,通过测量经过微透镜聚焦在电荷耦合元件上的光斑坐标,利用泽尼克多项式模式法算出泽尼克系数来重建波前。该技术恢复速度快、灵敏度高,已被广泛运用于天文望远镜的高分辨成像、人眼视网膜细胞分辨成像等领域。Hartmann wavefront sensing technology mainly recovers phase information by measuring the wavefront slope (Platt B C, Shack R. History and principles of Shack-Hartmann wavefront sensing[J]. Journal of Refractive Surgery, 2001, 17(5): S573-S577.), based on this technology, the Shack Hartmann sensor is composed of a microlens array and a charge-coupled element. By measuring the coordinates of the light spot focused on the charge-coupled element through the microlens, the Zernike polynomial mode method is used to calculate the Zernike coefficients to reconstruct the wavefront. The technology has fast recovery speed and high sensitivity, and has been widely used in high-resolution imaging of astronomical telescopes and retinal cell-resolved imaging of human eyes.
全息干涉技术是利用干涉原理恢复物波光场的技术(Eisebitt S,Lüning J,Schlotter W F,et al.Lensless imaging of magnetic nanostructures by X-rayspectro-holography[J].Nature,2004,432(7019):885-888.),分为拍摄和恢复两个过程,当一束参考光与物光进行干涉,干涉图中将记录该物体的相位和振幅信息,再用参考光照射全息干涉图,即可以恢复出物光的光场信息,从而提取物体的相位和振幅信息。随着全息干涉技术的飞速发展,使用电荷耦合元件记录干涉图样,而由计算机进行相位恢复的过程,该技术成为数字全息技术,被广泛运用在三维图像重构、数字显微成像、材料无损探伤、医学诊断等方面。Holographic interference technology is a technology that uses the principle of interference to restore the optical field of matter waves (Eisebitt S, Lüning J, Schlotter WF, et al. Lensless imaging of magnetic nanostructures by X-rayspectro-holography [J]. Nature, 2004, 432(7019): 885-888.), which is divided into two processes: shooting and recovery. When a beam of reference light interferes with the object light, the phase and amplitude information of the object will be recorded in the interferogram, and then the holographic interferogram can be illuminated with the reference light. The light field information of the object light is recovered, so as to extract the phase and amplitude information of the object. With the rapid development of holographic interference technology, charge-coupled elements are used to record the interference pattern, and the computer performs the phase recovery process. This technology has become a digital holography technology, which is widely used in three-dimensional image reconstruction, digital microscopic imaging, and non-destructive testing of materials. , medical diagnosis, etc.
计算相位恢复方法相比于前两种方法,应用范围更广,可见光以及极紫外波段皆可适用。该恢复方法是1972年提出的(Gerchberg R W.A practical algorithm for thedetermination of phase from image and diffraction plane pictures[J].Optik,1972,35:237.),利用迭代算法从拍摄到的衍射光斑的强度信息中恢复相位信息,该恢复方法可以用于以X射线或者自由电子束为光源的成像系统中,实现无透镜衍射成像和相位恢复,从而减小成像系统误差,简化系统结构,因此具有极大的前景。另一种利用迭代算法的叠层成像技术(Rodenburg J M,Faulkner H M L.A phase retrieval algorithm forshifting illumination[J].Applied physics letters,2004,85(20):4795-4797.),是近十年来新兴的技术,通过寻找样本进行重叠扫描的模式下,满足多幅远场衍射强度图像约束的唯一复数解,因不受光学聚焦原件的限制,从而可以突破衍射极限实现超分辨成像。Compared with the previous two methods, the computational phase recovery method has a wider application range, and is applicable to both visible light and extreme ultraviolet wavelengths. This recovery method was proposed in 1972 (Gerchberg R W.A practical algorithm for the determination of phase from image and diffraction plane pictures[J]. Optik, 1972, 35:237.), which uses an iterative algorithm from the intensity information of the diffracted light spots captured. This method can be used in imaging systems using X-rays or free electron beams as light sources to achieve lensless diffraction imaging and phase recovery, thereby reducing imaging system errors and simplifying system structure. prospect. Another stack imaging technique using an iterative algorithm (Rodenburg J M, Faulkner H M L. A phase retrieval algorithm for shifting illumination[J]. Applied physics letters, 2004, 85(20): 4795-4797.) is an emerging technology in the past decade. By searching for a unique complex solution that satisfies the constraints of multiple far-field diffraction intensity images in the mode of overlapping scanning of the sample, it is not limited by the optical focusing element, so it can break through the diffraction limit to achieve super-resolution imaging.
然而这些技术都存在一定的缺点和弊端,计算相位恢复方法的迭代算法需要大量的迭代次数和迭代时间,对于复杂相位物体,无法实现信息的快速、实时恢复,甚至会出现得不到唯一解的情况。且现有的迭代算法处理部分相干光时,会假设光源的互相关函数为高斯谢尔模,利用模式展开进行处理,然而当相干度很低的情况下,需要很多的模式数才可以正确恢复相位信息,且对于关联函数未知的光源,该模式展开的方法将不再适用。However, these technologies all have certain shortcomings and drawbacks. The iterative algorithm for calculating the phase recovery method requires a large number of iterations and iteration time. For complex phase objects, the fast and real-time recovery of information cannot be achieved, and there may even be cases where a unique solution cannot be obtained. Happening. Moreover, when the existing iterative algorithm processes partially coherent light, it assumes that the cross-correlation function of the light source is a Gaussian Shell mode, and uses mode expansion for processing. However, when the coherence degree is very low, a large number of modes are required to correctly recover the phase information. , and for light sources whose correlation function is unknown, this mode expansion method will no longer apply.
基于以上部分相干光照明下的物体信息恢复方面存在的缺陷,本专利创新性的提出一种新型的部分相干光照明下的相位信息恢复非迭代方法,本方法具有应用范围广、恢复速度快、装置简单等优点,具有重要的科研及实际应用价值。Based on the above defects in object information recovery under partially coherent light illumination, this patent innovatively proposes a new non-iterative method for phase information recovery under partially coherent light illumination, which has the advantages of wide application range, fast recovery speed, The device is simple and has important scientific research and practical application value.
发明内容SUMMARY OF THE INVENTION
为解决上述技术问题,本发明的目的是提供一种部分相干光照明下的非迭代相位恢复装置及方法,避开迭代算法的冗长和复杂,弥补模式展开法的弊端,可以实现在传统关联或关联结构复杂甚至未知的光源照明下,相位物体信息的正确、实时恢复。In order to solve the above-mentioned technical problems, the purpose of the present invention is to provide a non-iterative phase recovery device and method under partial coherent light illumination, which avoids the tediousness and complexity of the iterative algorithm, makes up for the drawbacks of the mode expansion method, and can realize in traditional correlation or Correct and real-time recovery of phase object information under the illumination of light sources with complex or even unknown correlation structures.
本发明的部分相干光照明下的非迭代相位恢复装置,包括部分相干光产生单元和物体相位测量单元,所述物体相位测量单元包括The non-iterative phase recovery device under partially coherent light illumination of the present invention includes a partially coherent light generating unit and an object phase measuring unit, wherein the object phase measuring unit includes
-分束镜,用于透射所述部分相干光产生单元产生的部分相干光,并反射经由空间光调制器调制后的光束;- a beam splitter for transmitting the partially coherent light generated by the partially coherent light generating unit and reflecting the light beam modulated by the spatial light modulator;
-空间光调制器,垂直于所述部分相干光产生单元的光轴放置,用于加载待测相位物体和对待测相位物体进行相位扰动的扰动点,所述空间光调制器反射所述分束镜透射的光,并让调制后的光经过所述分束镜再次发生反射;- a spatial light modulator, placed perpendicular to the optical axis of the partially coherent light generating unit, for loading the phase object to be measured and the perturbation point for phase perturbation of the phase object to be measured, the spatial light modulator reflects the split beam the light transmitted by the mirror, and let the modulated light reflect again through the beam splitter;
-多孔阵列板,供所述分束镜反射的光束穿过,所述多孔阵列板上设有周期排列的二维小孔阵列并在阵列中心附近设有一个参考小孔,所述多孔阵列板上的参考小孔对准所述分束镜反射的光束,所述多孔阵列板与空间光调制器之间的距离满足z≥d*L/λ,其中,d是多孔阵列板上小孔间的间隔,L是待测相位物体最宽处的尺寸,λ是部分相干光产生单元中激光光源的波长;- a perforated array plate through which the beams reflected by the beam splitter pass, the perforated array plate is provided with a periodically arranged two-dimensional array of small holes and a reference hole near the center of the array, the perforated array plate The reference hole on the beam splitter is aligned with the beam reflected by the beam splitter, and the distance between the multi-hole array plate and the spatial light modulator satisfies z ≥ d*L/λ, where d is the distance between the holes on the multi-hole array plate interval, L is the size of the widest point of the phase object to be measured, λ is the wavelength of the laser light source in the partially coherent light generating unit;
-傅里叶透镜,紧挨所述多孔阵列板之后放置,或能够使所述多孔阵列板位于傅里叶透镜的前焦面上,用于对穿过所述多孔阵列板的光束进行傅里叶变换;- a Fourier lens, placed immediately after the multiwell array plate, or enabling the multiwell array plate to be located on the front focal plane of the Fourier lens, for Fourier analysis of the light beam passing through the multiwell array plate leaf transform;
-电荷耦合元件,放置在傅里叶平面处拍摄光强信息;- Charge-coupled element, placed at the Fourier plane to capture light intensity information;
-计算机,与所述空间光调制器和电荷耦合元件连接,控制所述空间光调制器上的相位加载,并对拍摄得到的光强进行实时反傅里叶变换、筛选以及反传输处理,获得物体的相位信息。- a computer, connected with the spatial light modulator and the charge-coupled element, controls the phase loading on the spatial light modulator, and performs real-time inverse Fourier transform, screening and inverse transmission processing on the light intensity obtained by shooting to obtain Phase information of the object.
进一步的,当目的为产生传统高斯关联的部分相干光时,所述部分相干光产生单元包括依次设置的激光器、对激光器发出的激光束进行扩束的扩束镜、对光束进行准直的准直透镜和对光束进行整形的高斯滤波片,由所述高斯滤波片出来的光透射过所述分束镜,到达所述空间光调制器。Further, when the purpose is to generate the partially coherent light associated with the traditional Gaussian, the partially coherent light generating unit includes a laser arranged in sequence, a beam expander that expands the laser beam emitted by the laser, and a collimator that collimates the beam. A straight lens and a Gaussian filter for shaping the beam, the light from the Gaussian filter is transmitted through the beam splitter and reaches the spatial light modulator.
进一步的,当目的为产生拉盖尔高斯关联的部分相干光时,所述部分相干光产生单元包括依次设置的所述激光器、对激光器发出的激光束进行扩束的扩束镜、对扩束后的光进行相位改变的螺旋位相板、对光束进行准直的准直透镜和对光束进行整形的高斯滤波片,由所述高斯滤波片出来的光透射过所述分束镜,到达所述空间光调制器。Further, when the purpose is to generate the partially coherent light of the Laguerre-Gaussian correlation, the partially coherent light generating unit comprises the laser arranged in sequence, the beam expander for expanding the laser beam emitted by the laser, and the beam expander. The helical phase plate for changing the phase of the rear light, the collimating lens for collimating the beam, and the Gaussian filter for shaping the beam, the light from the Gaussian filter is transmitted through the beam splitter and reaches the Spatial Light Modulator.
进一步的,所述部分相干光产生单元还包括相干度调节组件,所述相干度调节组件包括对经所述扩束镜扩束后的光束或经所述螺旋位相板改变相位的光束进行聚焦的透镜,以及对聚焦后的光束进行散射的旋转毛玻璃,经所述旋转毛玻璃散射出来的光由所述准直透镜进行准直。Further, the partially coherent light generating unit further includes a coherence degree adjustment component, and the coherence degree adjustment component includes a focusing device for the beam expanded by the beam expander or the beam whose phase is changed by the helical phase plate. a lens, and a rotating frosted glass that scatters the focused light beam, and the light scattered by the rotating frosted glass is collimated by the collimating lens.
进一步的,所述分束镜反射的光束还可通过反射镜反射到所述多孔阵列板上。Further, the light beam reflected by the beam splitter can also be reflected to the porous array plate by a mirror.
进一步的,所述扰动点位于待测相位物体的任意位置,其尺寸远小于待测物体尺寸,其相位赋值有异于原待测相位物体于该位置的相位即可。Further, the disturbance point is located at any position of the object to be measured, its size is much smaller than the size of the object to be measured, and its phase assignment is different from the phase of the original object to be measured at this position.
进一步的,通过激光在不透光的基板上刻蚀一个参考小孔和以参考小孔为圆心向x和y方向各偏移一定距离再对称激光刻蚀二维小孔阵列,形成所述多孔阵列板,各二维小孔间的间隔需满足d≤z*λ/L,其中z为物平面到所述多孔阵列板的距离,L是待测相位物体最宽处的尺寸,λ是部分相干光产生单元中激光光源的波长;偏移量的大小满足a/2≤Δx=Δy≤d/2-a/2,其中a为二维小孔的尺寸;参考小孔和各二维小孔的尺寸一致,需满足远小于待测物体尺寸,且小于二维小孔间隔的三分之一;所述多孔阵列板的结构也可以通过所述空间光调制器模拟。Further, a reference small hole is etched on the opaque substrate by laser, and the reference small hole is used as the center of the circle to offset a certain distance in the x and y directions, and then the two-dimensional small hole array is symmetrically laser etched to form the porous structure. For the array plate, the interval between the two-dimensional holes should satisfy d≤z*λ/L, where z is the distance from the object plane to the porous array plate, L is the size of the widest point of the phase object to be measured, and λ is the part The wavelength of the laser light source in the coherent light generating unit; the size of the offset satisfies a/2≤Δx=Δy≤d/2-a/2, where a is the size of the two-dimensional small hole; the reference small hole and each two-dimensional small hole are The size of the holes is the same, which needs to be much smaller than the size of the object to be measured, and less than one-third of the interval between the two-dimensional small holes; the structure of the porous array plate can also be simulated by the spatial light modulator.
进一步的,所述空间光调制器为反射式纯相位空间光调制器。Further, the spatial light modulator is a reflective pure-phase spatial light modulator.
进一步的,所述分束镜是光强50:50的半透半反镜。Further, the beam splitter is a half mirror with a light intensity of 50:50.
本发明的利用上述部分相干光照明下的非迭代相位恢复装置进行非迭代相位恢复的方法,包括步骤:The method for non-iterative phase recovery using the above-mentioned non-iterative phase recovery device under partially coherent light illumination of the present invention includes the steps of:
(1)通过部分相干光产生单元获取高斯关联的部分相干光或拉盖尔高斯关联的部分相干光,通过调整用来聚焦的所述透镜在光轴上的前后位置,来改变聚焦在所述旋转毛玻璃上的光斑大小以调整光源的空间相干性;(1) Obtaining partially coherent light associated with Gaussian or partially coherent light associated with Laguerre Gaussian through the partially coherent light generating unit, and changing the focus on the Rotate the spot size on the frosted glass to adjust the spatial coherence of the light source;
(2)将部分相干光传输至所述物体相位测量单元,使部分相干光经过待测相位物体并记录光强信息:(2) The partial coherent light is transmitted to the object phase measurement unit, so that the partial coherent light passes through the object to be measured and records the light intensity information:
(21)当部分相干光为高斯关联的部分相干光时,记录两次光强,第一次拍摄——加载待测相位物体:将准备好的图片设置为灰度图模式,当该图片加载到纯相位的所述空间光调制器上时,灰度值将转化为对应的相位;第二次拍摄——在待测相位物体任意位置放置一个扰动点,其尺寸需远小于待测物体尺寸,其相位赋值有异于原待测相位物体于该位置的相位即可;将两次拍摄得到的光强信息,输送到计算机进行处理:首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出物体相位信息;所述筛选阵列的筛选小孔为周期排列的二维小孔,由所述计算机产生,各筛选小孔的边长及间隔与所述多孔阵列板一致;(21) When the partially coherent light is partially coherent light with Gaussian correlation, record the light intensity twice, the first shot - load the phase object to be measured: set the prepared image to grayscale image mode, when the image is loaded When it reaches the pure-phase spatial light modulator, the gray value will be converted into the corresponding phase; the second shot - place a disturbance point at any position of the object to be measured, and its size needs to be much smaller than the size of the object to be measured. , the phase assignment is different from the phase of the original phase object to be measured at this position; the light intensity information obtained by the two shots is sent to the computer for processing: First, the two groups of light intensity information are respectively Fourier transform, Then, it is screened by the screening array, and then the screened results are subtracted, and the subtracted results are reversely transmitted to extract the phase information of the object; the screening holes of the screening array are periodically arranged two-dimensional holes , generated by the computer, and the side length and interval of each screening hole are consistent with the multi-hole array plate;
(22)当部分相干光为拉盖尔高斯关联的部分相干光时,记录四次光强,第一次拍摄——加载待测相位物体:将准备好的图片设置为灰度图模式,当该图片加载到纯相位的所述空间光调制器上时,灰度值将转化为对应的相位;第二次拍摄——在待测相位物体任意位置放置一个扰动点,其尺寸远小于待测物体尺寸,其相位赋值有异于原待测相位物体于该位置的相位即可;第三次拍摄——移走待测相位物体但不移走扰动点,即将原待测相位物体所在区域的相位设置为0;第四次拍摄——移走扰动相位;将第三次和第四次拍摄得到的两次光强信息,输送到计算机进行处理:首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出光源的交叉谱密度信息及其相位信息;将第一次和第二次拍摄得到的两次光强信息,输送到计算机进行处理:首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出受光源相位影响的物体相位信息;最后用受光源相位影响的物体相位除去光源的相位,即可获得正确的物体相位信息;所述筛选阵列的筛选小孔为周期排列的二维小孔,由所述计算机产生,各筛选小孔的边长及间隔与所述多孔阵列板一致。(22) When the partially coherent light is the partially coherent light of the Laguerre-Gaussian correlation, record the light intensity four times, and take the first shot - load the object to be measured: set the prepared image to grayscale mode, when When the picture is loaded on the pure-phase spatial light modulator, the gray value will be converted into the corresponding phase; the second shot - place a disturbance point at any position of the phase object to be measured, and its size is much smaller than that of the phase to be measured. The size of the object, its phase assignment is different from the phase of the original phase object to be measured at this position; the third shot - remove the object to be measured but do not remove the disturbance point, that is, the area where the original phase object to be measured is located. The phase is set to 0; the fourth shot—removes the disturbed phase; the two light intensity information obtained from the third and fourth shots are sent to the computer for processing: first, Fourier is performed on the two sets of light intensity information respectively. Leaf transformation, and then screened by the screening array, secondly, the filtered results are subtracted, and the subtracted results are inversely transmitted to extract the cross-spectral density information and phase information of the light source; The two light intensity information obtained by the second shooting is sent to the computer for processing: first, the two sets of light intensity information are subjected to Fourier transform respectively, and then screened by the screening array, and then the screened results are subtracted, and the subtracted After inverse transmission of the result, the phase information of the object affected by the phase of the light source can be extracted; finally, the phase of the light source can be removed by the phase of the object affected by the phase of the light source, and the correct phase information of the object can be obtained; the screening of the screening array is small. The holes are periodically arranged two-dimensional small holes, which are generated by the computer, and the side length and interval of each screening small hole are consistent with the multi-hole array plate.
借由上述方案,本发明至少具有以下优点:By means of the above scheme, the present invention has at least the following advantages:
针对部分相干光照明下物体相位信息的恢复,与迭代算法相比,本发明更为快捷、实时化,较模式展开这种处理部分相干光的迭代算法,应用范围更广,对于处理关联结构复杂甚至未知的部分相干光照明情况下,物体相位的恢复有着独特的优势;本发明的无透镜衍射成像的相位恢复装置,可以拓展到X射线成像系统中,因而在实际应用中有着极为重要的意义。Compared with the iterative algorithm, the present invention is more rapid and real-time for the recovery of the phase information of the object under partial coherent light illumination. Compared with the mode expansion, the iterative algorithm for processing partial coherent light has a wider range of applications, and has a wider range of applications for processing complex correlation structures. Even in the case of unknown partially coherent light illumination, the phase recovery of the object has unique advantages; the lensless diffraction imaging phase recovery device of the present invention can be extended to the X-ray imaging system, so it has extremely important significance in practical applications .
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。The above description is only an overview of the technical solution of the present invention. In order to understand the technical means of the present invention more clearly, and implement it according to the content of the description, the preferred embodiments of the present invention are described in detail below with the accompanying drawings.
附图说明Description of drawings
图1是本发明的部分相干光照明下的非迭代相位恢复装置的结构示意图;1 is a schematic structural diagram of a non-iterative phase recovery device under partially coherent light illumination of the present invention;
图2是在空间光调制器上加载待测相位物体和扰动的一个示例,其中图2(a)为只加载物体信息,图2(b)在物体信息上加载了一个扰动;Fig. 2 is an example of loading the object to be measured and disturbance on the spatial light modulator, wherein Fig. 2(a) only loads the object information, and Fig. 2(b) loads a disturbance on the object information;
图3(a)是本发明所用的多孔阵列板的中心部分细节图;图3(b)是计算机恢复时使用的筛选阵列的中心部分细节图。Fig. 3(a) is a detail view of the center portion of the multiwell array plate used in the present invention; Fig. 3(b) is a detail view of the center portion of the screening array used in computer recovery.
具体实施方式Detailed ways
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。The specific embodiments of the present invention will be described in further detail below with reference to the accompanying drawings and embodiments. The following examples are intended to illustrate the present invention, but not to limit the scope of the present invention.
本发明提出的相位恢复方法包含三个过程:产生高斯关联或者特殊关联的部分相干光源、经过携带相位信息的物体并记录光强和计算机处理以恢复相位。其对应的结构装置如附图1所示。The phase recovery method proposed by the present invention includes three processes: generating a partially coherent light source with Gaussian correlation or special correlation, passing through an object carrying phase information and recording the light intensity, and computer processing to recover the phase. The corresponding structural device is shown in FIG. 1 .
首先产生部分相干光源,照射到待测物体后经过多孔阵列板,利用电荷耦合元件记录傅里叶平面的光斑光强。一般情况下,恢复待测物体完整的相位信息只需要进行两次拍摄。第一次拍摄,是让光源照射待测物体后,传输一段距离到达经过特别设计的多孔阵列板,再由放置在傅里叶平面的电荷耦合元件记录光强信息;第二次拍摄,在待测物体的中间添加扰动点,再经过同样的传输过程,由电荷耦合元件记录光强信息。两次拍摄后,利用计算机程序处理,恢复出物光信息。但是若光源不是高斯关联,恢复得到的物光信息会受到光源关联结构的影响,在这样的情况下,需要移走物体,其他不变,按照上述的步骤,再进行两次拍摄,并由计算机处理恢复,得到光源的交叉谱密度方程,最后,将两次恢复结果相除即可得到正确的物体振幅和相位信息。First, a partially coherent light source is generated, which irradiates the object to be measured and passes through the multi-hole array plate, and uses the charge-coupled element to record the intensity of the spot light in the Fourier plane. Generally, only two shots are required to recover the complete phase information of the object to be measured. The first shot is to let the light source illuminate the object to be tested, and then transmit a distance to the specially designed multi-hole array plate, and then record the light intensity information by the charge-coupled element placed on the Fourier plane. The disturbance point is added in the middle of the measuring object, and then the light intensity information is recorded by the charge-coupled element after the same transmission process. After two shots, the computer program is used to recover the object light information. However, if the light source is not Gaussian correlation, the recovered object light information will be affected by the light source correlation structure. In this case, the object needs to be removed, and the rest remain unchanged. According to the above steps, two more shots are taken, and the computer After processing the recovery, the cross spectral density equation of the light source is obtained. Finally, the correct object amplitude and phase information can be obtained by dividing the two recovery results.
本发明中高斯关联的部分相干光产生单元的结构包括:激光器1发出的激光束,由扩束镜2进行扩束,再由聚焦透镜4聚焦到旋转毛玻璃片5上,出射光由准直透镜6进行准直,并由高斯滤波片7进行整形。聚焦透镜4和旋转毛玻璃片5构成相干度调节系统,通过改变聚焦透镜4和旋转毛玻璃片5的位置,可以改变出射光相干度的大小,这是因为聚焦到旋转毛玻璃片5上的光斑大小直接影响着出射光的相干性,聚焦光斑越大,相干性越低,聚焦光斑越小,相干性越高。经过高斯滤波片7出来以后的光束就是高斯关联的部分相干光,即所需的光源。若在扩束镜2与聚焦透镜4间放置一个螺旋相位板3,而其他不变,则高斯滤波片7后产生的就是拉盖尔高斯关联的部分相干光,即所需的特殊关联的部分相干光源。The structure of the partially coherent light generating unit of Gaussian correlation in the present invention includes: the laser beam emitted by the
本发明中物体相位测量单元的结构具体包括:产生部分相干光源后,经过一个分束镜8,垂直入射到空间光调制器9上,该分束镜8是光强50:50的半透半反镜,空间光调制器9上加载了携带相位信息的物体,经过空间光调制器9后,由其反射出来的光经过分束镜8,再反射到多孔阵列板13上,空间光调制器9到多孔阵列板13具有一定距离,该距离满足公式z≥d*L/λ,其中d是多孔阵列板13上小孔间的间隔,L是待测物体最宽处的尺寸,λ是激光光源的波长。多孔阵列板13由周期排列的二维小孔阵列组成,并在阵列中间有一个参考小孔。到达多孔阵列板13的光束需要对准参考孔,傅里叶透镜12紧挨着多孔阵列板13之后放置,或能够使多孔阵列板位于傅里叶透镜的前焦面上,最后由电荷耦合元件11在傅里叶平面处拍摄光强信息。电荷耦合元件这里拍摄得到的光强需要输送到计算机10进行实时反傅里叶变换、筛选阵列筛选以及反传输处理,以获得物体的相位信息。The structure of the object phase measurement unit in the present invention specifically includes: after generating a partially coherent light source, it is vertically incident on the spatial light modulator 9 through a
本发明的依据和原理如下:The basis and principle of the present invention are as follows:
将光源交叉谱密度表示为W0(ρ1,ρ2),将待测物体表示为O(ρ),那光源照射到待测物体后,传输到多孔阵列板的交叉谱密度方程可以表示为:The cross spectral density of the light source is expressed as W 0 (ρ 1 ,ρ 2 ), and the object to be measured is expressed as O(ρ). After the light source illuminates the object to be measured, the cross spectral density equation transmitted to the multi-hole array plate can be expressed as :
W(r1,r2)=∫∫W0(ρ1,ρ2)O(ρ1)O(ρ2)*P(ρ1,r1)P(ρ2,r2)*dρ1dρ2 (1)W(r 1 ,r 2 )=∫∫W 0 (ρ 1,ρ2 )O(ρ 1 )O(ρ 2 ) * P(ρ 1 ,r 1 )P(ρ 2 ,r 2 ) * dρ 1 dρ 2 (1)
其中P(ρ,r)是从物平面到多孔阵列板的传输项。多孔阵列板可以用δ函数表示:where P(ρ,r) is the transfer term from the object plane to the multiwell array plate. Well array plates can be represented by the delta function:
M(r)=δ(r)+∑mnδ(r-rmn) (2)M(r)=δ(r)+∑ mn δ(rr mn ) (2)
其中rmn=(md+Δx,nd+Δy)是多孔阵列板上周期小孔的坐标,m和n是整数,d是阵列孔之间的间距,Δx和Δy表示的是参考孔附近的阵列孔在x和y方向由中心点产生的偏移,δ(r)表示中心的参考孔。经过多孔阵列板到达傅里叶平面的光强I(κ),再进行反傅里叶变换,对应的则是光经过多孔阵列板后的光场的交叉谱密度方程:where r mn =(md+Δx,nd+Δy) are the coordinates of the periodic holes on the multiwell array plate, m and n are integers, d is the spacing between the array holes, and Δx and Δy represent the array near the reference hole The offset of the hole in the x and y directions from the center point, δ(r) is the reference hole at the center. The light intensity I(κ) reaches the Fourier plane after passing through the porous array plate, and then performs inverse Fourier transform, which corresponds to the cross-spectral density equation of the light field after the light passes through the porous array plate:
其中rmn-rpq=[(m-p)d,(n-q)d],这里的p和q是整数。紧接着,用计算机程序模拟让反傅里叶变换后的交叉谱密度经过一个筛选阵列,筛选阵列和多孔阵列板的分布类似,但是缺少一个参考小孔。经过该筛选阵列,可以滤出W(-rmn,0)或者W(0,rmn)*,因为多孔阵列并不是严格以参考小孔为中心,因此,只能筛选出一个交叉谱密度方程。由公式(3)得到:where r mn -r pq =[(mp)d,(nq)d], where p and q are integers. Next, a computer program was used to simulate the inverse Fourier-transformed cross-spectral density through a screening array that had a similar distribution to the multiwell array plate, but lacked a reference well. After this screening array, W(-r mn ,0) or W(0,r mn ) * can be filtered out, because the porous array is not strictly centered on the reference hole, therefore, only one cross spectral density equation can be screened out . It is obtained by formula (3):
其中TO(ρ1,ρ2)=∫W0(ρ1,ρ2)O(ρ2)*P(ρ2,0)*dρ2,此时如果对W(rmn,0)进行反传输,得到的是TO(ρ1,ρ2)O(ρ1),这并不是所要的正确信息,因此,还需要进行第二次拍摄。在待测物体上加载一个扰动,数学上表示为Cδ(ρ-ρ0),其中,C是一个复数常数,ρ0是扰动的坐标。第二次拍摄得到的W′(rmn,0)表示为:where T O (ρ 1 ,ρ 2 )=∫W 0 (ρ 1 ,ρ 2 )O(ρ 2 ) * P(ρ 2 ,0) * dρ 2 , then if W(r mn ,0) is performed Inverse transmission, what is obtained is T O (ρ 1 ,ρ 2 )O(ρ 1 ), which is not the correct information required, therefore, a second shot is required. A disturbance is loaded on the object to be measured, which is mathematically expressed as Cδ(ρ-ρ 0 ), where C is a complex constant and ρ 0 is the coordinate of the disturbance. The W′(r mn ,0) obtained from the second shot is expressed as:
前后结果进行相减得到:Subtract the before and after results to get:
此时再进行反向传输,就只剩下W0(ρ1,ρ0)O(ρ1),对于高斯关联的部分相干光,其相位为1,那恢复出来的相位就是物体O(ρ1)携带的相位,当照明光源不是高斯关联而是特殊关联的部分相干光时,需要在光路中移去物体,用同样的实验和处理方法得到W0(ρ1,ρ0),再将前后结果相除,即可得到正确的相位信息。At this time, the reverse transmission is performed, and only W 0 (ρ 1 ,ρ 0 )O(ρ 1 ) is left. For the partially coherent light associated with Gaussian, its phase is 1, and the recovered phase is the object O(ρ 1 ) The phase carried, when the illumination light source is not Gaussian correlation but partially coherent light with special correlation, the object needs to be removed in the optical path, and W 0 (ρ 1 , ρ 0 ) is obtained by the same experiment and processing method, and then the The correct phase information can be obtained by dividing the before and after results.
下面结合附图及具体实施方式对本发明作进一步说明。The present invention will be further described below with reference to the accompanying drawings and specific embodiments.
实施例一:高斯关联部分相干光照明下物体相位信息恢复。Embodiment 1: Recovery of phase information of an object under Gaussian correlation partially coherent light illumination.
1、高斯关联的部分相干光源的产生:其结构包括功率可调的半导体泵浦固体激光器1,其发射出的激光波长为532nm,该激光器发出的激光束由扩束镜2进行扩束,扩束后的光束由聚焦透镜4聚焦到旋转毛玻璃片5上,散射出来的光由准直透镜6进行准直,再由高斯滤波片7进行整形。由高斯滤波片出来的光即为高斯关联的部分相干光。这里,聚焦透镜4的焦距为80mm,准直透镜6的焦距为150mm。1. Generation of partially coherent light source with Gaussian correlation: its structure includes a power-adjustable semiconductor-pumped solid-
1.1、调整光源的空间相干性:光源的空间相干性和聚焦在旋转毛玻璃片上的光斑大小有着直接的关系,因此,通过调整聚焦透镜4在光轴上的前后位置,来改变聚焦在旋转毛玻璃片5上的光斑大小,当聚焦光斑越小时,出射光的相干性越高,反之聚焦光斑越大,出射光的相干性越低。1.1. Adjust the spatial coherence of the light source: the spatial coherence of the light source is directly related to the spot size focused on the rotating ground glass sheet. Therefore, by adjusting the front and rear positions of the focusing
2、光源经过物体并记录光强的单元:上一步骤产生的部分相干光源透射过分束镜8,到达反射式纯相位空间光调制器9。该分束镜是光强50:50的半透半反镜,反射式空间光调制器垂直于部分相干光产生单元的光轴放置。空间光调制器9连接着计算机10,由计算机10控制空间光调制器上的相位加载。经相位调制的出射光再次经过分束镜,发生反射,为节省装置所占的空间,该反射光经过一个反射镜14再次发生反射,传输1170mm后到达多孔阵列板13。光束由空间光调制器到多孔阵列板的传输距离需满足条件z≥d*L/λ,其中d是多孔阵列板小孔的间隔,L是待测物体最宽处的尺寸,λ是激光光源的波长。紧挨着多孔阵列板后,放置一个傅里叶透镜12,其焦距为150mm,连接着计算机10的电荷耦合元件11放置在傅里叶平面处,记录光强。2. The unit where the light source passes through the object and records the light intensity: the partially coherent light source generated in the previous step is transmitted through the
整个过程共需要两次记录光强,两次拍摄过程中的唯一区别在于空间光调制器9上的相位加载不同:第一次拍摄,加载待测相位物体,如图2(a);第二次拍摄,在待测相位物体中间放置一个扰动点,如图2(b),该扰动点是边长为240μm的正方形,其相位赋值等于在原待测物体于该位置的相位再减去0.8π,以达到扰动效果。The whole process needs to record the light intensity twice. The only difference between the two shooting processes is that the phase loading on the spatial light modulator 9 is different: the first shooting, the object to be measured is loaded, as shown in Figure 2(a); the second For the second shooting, place a disturbance point in the middle of the object to be measured, as shown in Figure 2(b). The disturbance point is a square with a side length of 240 μm, and its phase assignment is equal to the phase of the original object to be measured at this position minus 0.8π , in order to achieve the perturbation effect.
2.1、在反射式空间光调制器上加载相位信息:首先,将准备好的图片设置为灰度图模式。当该图片加载到纯相位的空间光调制器上(型号Holoeye-Pluto,像素尺寸1920×1080,像素大小8μm),灰度值将转化为对应的相位。2.1. Load phase information on the reflective spatial light modulator: First, set the prepared image to grayscale image mode. When the image is loaded on a phase-only spatial light modulator (model Holoeye-Pluto, pixel size 1920×1080, pixel size 8μm), the gray value will be converted to the corresponding phase.
2.2、设计实验中使用的多孔阵列板:利用激光刻蚀,制作一个18mm×14mm的多孔阵列板,基片整体不透光,再进行打孔,其中原点位置放置一个54μm边长的正方形小孔,再从原点向x和y方向各偏移Δx=117μm和Δy=117μm,以此为圆心对称布置小孔阵列,如图3(a),其中x方向66个小孔,y方向48个小孔,每个小孔边长均为54μm,小孔间的间隔d=270μm。值得注意的是Δx=Δy≠d/2。(多孔阵列板也可以用透射式空间光调制器代替)。实验中,到达多孔阵列板的光束要对准参考孔。2.2. The multi-hole array plate used in the design experiment: use laser etching to make a 18mm×14mm multi-hole array plate, the substrate as a whole is opaque to light, and then drill holes, where a square hole with a side length of 54 μm is placed at the origin. , and then offset Δx=117 μm and Δy=117 μm from the origin to the x and y directions respectively, and arrange the small hole array symmetrically with this as the center of the circle, as shown in Figure 3(a), in which there are 66 small holes in the x direction and 48 small holes in the y direction. The side length of each small hole is 54 μm, and the interval between the small holes is d=270 μm. It is worth noting that Δx=Δy≠d/2. (The multiwell array plate can also be replaced with a transmissive spatial light modulator). In the experiment, the beam reaching the multiwell array plate is aimed at the reference well.
3、计算恢复相位:两次拍摄得到的光强信息,输送到计算机10,进行处理。首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出物体相位信息。3. Calculate the recovery phase: the light intensity information obtained by the two shots is sent to the
3.1、设计计算机恢复时使用的筛选阵列:筛选阵列不是实际的物体,而是在程序处理时需要用到的筛选信息用的阵列。其分布如图3(b),与实验用的多孔阵列板的唯一区别在于缺少中间的参考小孔,其他参数一致。3.1. Design the screening array used in computer recovery: the screening array is not an actual object, but an array for screening information that needs to be used in program processing. Its distribution is shown in Figure 3(b), the only difference from the experimental multi-well array plate is the lack of the reference hole in the middle, and other parameters are consistent.
实施例二:特殊关联部分相干光照明下物体相位信息恢复。Embodiment 2: Recovering the phase information of an object under the illumination of special associated partially coherent light.
1、特殊关联的部分相干光源的产生:其结构包括功率可调的半导体泵浦固体激光器1,其波长为532nm,该激光器发出的激光束由扩束镜2进行扩束,扩束后的光束经过拓扑荷数为2的螺旋位相板3,再由聚焦透镜4聚焦到旋转毛玻璃片5上,散射出来的光由准直透镜6进行准直,再由高斯滤波片7进行整形。这里,聚焦透镜4的焦距为80mm,准直透镜6的焦距为150mm。从高斯滤波片出来的光即为拉盖尔高斯关联的部分相干光。用拉盖尔高斯关联的部分相干光作为照明光源,来阐述关联结构复杂情况下的相位恢复过程。1. The generation of a special associated partially coherent light source: its structure includes a power-adjustable semiconductor-pumped solid-
1.1、调整光源的空间相干性:光源的空间相干性和聚焦在毛玻璃片上的光斑大小有着直接的关系,因此,通过调整聚焦透镜4在光轴上的前后位置,来改变聚焦在旋转毛玻璃片5上的光斑大小,当聚焦光斑越小时,出射光的相干性越高,反之聚焦光斑越大,出射光的相干性越低。1.1. Adjust the spatial coherence of the light source: The spatial coherence of the light source is directly related to the spot size focused on the frosted glass sheet. Therefore, by adjusting the front and rear positions of the focusing
2、光源经过物体并记录光强的单元:上一步骤产生的拉盖尔高斯关联的部分相干光源透射过分束镜8,到达反射式纯相位空间光调制器9。该分束镜是光强50:50的半透半反镜,反射式空间光调制器垂直于部分相干光产生单元的光轴放置。空间光调制器9连接着计算机10,由计算机10控制空间光调制上的相位加载。经相位调制的出射光再次经过分束镜,发生反射,为节省装置所占的空间,该反射光经过一个反射镜14再次发生反射,传输1170mm后到达多孔阵列板13。光束由空间光调制器到多孔阵列板的传输距离需满足条件z≥d*L/λ,其中d是多孔阵列板小孔的间隔,L是待测物体最宽处的尺寸,λ是激光光源的波长。紧挨着多孔阵列板后,放置一个傅里叶透镜12,其焦距为150mm,连接着计算机10的电荷耦合元件11放置在傅里叶平面处,记录光强。2. The unit where the light source passes through the object and records the light intensity: the partially coherent light source of the Laguerre-Gaussian correlation generated in the previous step is transmitted through the
整个过程共需要4次光强记录,4次拍摄过程中的唯一区别在于空间光调制器9上的相位加载不同:第1次拍摄,加载待测相位物体,如图2(a);第2次拍摄,在待测相位物体中间放置一个扰动点,如图2(b),该扰动点是边长为240μm的正方形,其相位赋值等于在原待测物体于该位置的相位再减去0.8π,以达到扰动效果;第3次拍摄,移走相位物体而不移走扰动加载(移走相位物体指的并不是移走空间光调制器,而是将原物体所在区域的相位设置为0);第4次拍摄,移走扰动加载。The whole process requires a total of 4 light intensity recordings. The only difference in the 4 shooting processes is that the phase loading on the spatial light modulator 9 is different: the first shooting, loading the phase object to be measured, as shown in Figure 2(a); the second For the second shooting, place a disturbance point in the middle of the object to be measured, as shown in Figure 2(b). The disturbance point is a square with a side length of 240 μm, and its phase assignment is equal to the phase of the original object to be measured at this position minus 0.8π , to achieve the perturbation effect; for the third shot, remove the phase object without removing the perturbation loading (removing the phase object does not mean removing the spatial light modulator, but setting the phase of the area where the original object is located to 0) ; 4th shot, remove perturbation loading.
2.1、在反射式空间光调制器上加载相位信息:首先,将准备好的图片设置为灰度图模式。当该图片加载到纯相位的空间光调制器上(型号Holoeye-Pluto,像素尺寸1920×1080,像素大小8μm),灰度值将转化为对应的相位。2.1. Load phase information on the reflective spatial light modulator: First, set the prepared image to grayscale image mode. When the image is loaded on a phase-only spatial light modulator (model Holoeye-Pluto, pixel size 1920×1080, pixel size 8μm), the gray value will be converted to the corresponding phase.
2.2、设计实验中使用的多孔阵列板:利用激光刻蚀,制作一个18mm×14mm的多孔阵列板,基片整体不透光,再进行打孔,其中原点位置放置一个54μm边长的正方形小孔,再从原点向x和y方向各偏移Δx=117μm和Δy=117μm,以此为圆心对称布置小孔阵列,如图3(a),其中x方向66个小孔,y方向48个小孔,每个小孔边长均为54μm,小孔间的间隔d=270μm。值得注意的是Δx=Δy≠d/2。(多孔阵列板也可以用透射式空间光调制器代替)。实验中,到达多孔阵列板的光束要对准参考孔。2.2. The multi-hole array plate used in the design experiment: use laser etching to make a 18mm×14mm multi-hole array plate, the substrate as a whole is opaque to light, and then drill holes, where a square hole with a side length of 54 μm is placed at the origin. , and then offset Δx=117 μm and Δy=117 μm from the origin to the x and y directions respectively, and arrange the small hole array symmetrically with this as the center of the circle, as shown in Figure 3(a), in which there are 66 small holes in the x direction and 48 small holes in the y direction. The side length of each small hole is 54 μm, and the interval between the small holes is d=270 μm. It is worth noting that Δx=Δy≠d/2. (The multiwell array plate can also be replaced with a transmissive spatial light modulator). In the experiment, the beam reaching the multiwell array plate is aimed at the reference well.
3、计算恢复相位:针对上述第3和第4次移走相位物体的拍摄结果:两次拍摄得到的光强信息,输送到计算机10,进行处理。首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出光源的交叉谱密度信息及其相位信息。针对上述第1和第2次带有相位物体的拍摄结果:两次拍摄得到的光强信息,输送到计算机10,进行处理。首先对两组光强信息分别进行傅里叶变换,再由筛选阵列筛选,其次将筛选后的结果相减,并将相减后的结果进行反传输,即可提取出受光源相位影响的物体相位信息。最后用受光源相位影响的物体相位除去光源的相位,即可获得正确的物体相位信息。3. Calculating the restored phase: for the shooting results of the above-mentioned 3rd and 4th time-removed phase objects: the light intensity information obtained by the two shots is sent to the
3.1、设计计算机恢复时使用的筛选阵列:筛选阵列不是实际的物体,而是在程序处理时需要用到的筛选信息用的阵列。其分布如图3(b),与实验用的多孔阵列板的唯一区别在于缺少中间的参考小孔,其他参数一致。3.1. Design the screening array used in computer recovery: the screening array is not an actual object, but an array for screening information that needs to be used in program processing. Its distribution is shown in Figure 3(b), the only difference from the experimental multi-well array plate is the lack of the reference hole in the middle, and other parameters are consistent.
整个过程包含四次光强记录和数据处理,处理过程简单,因此整个过程耗时极短,几乎可以实现实时恢复。The whole process includes four times of light intensity recording and data processing. The processing process is simple, so the whole process takes very little time, and almost real-time recovery can be achieved.
以上所述仅是本发明的优选实施方式,并不用于限制本发明,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. It should be pointed out that for those skilled in the art, some improvements can be made without departing from the technical principles of the present invention. These improvements and modifications should also be regarded as the protection scope of the present invention.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710197745.5A CN106990694B (en) | 2017-03-29 | 2017-03-29 | A non-iterative phase recovery device and method under partially coherent light illumination |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710197745.5A CN106990694B (en) | 2017-03-29 | 2017-03-29 | A non-iterative phase recovery device and method under partially coherent light illumination |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106990694A CN106990694A (en) | 2017-07-28 |
CN106990694B true CN106990694B (en) | 2022-07-12 |
Family
ID=59412401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710197745.5A Active CN106990694B (en) | 2017-03-29 | 2017-03-29 | A non-iterative phase recovery device and method under partially coherent light illumination |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106990694B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107329275B (en) * | 2017-08-30 | 2020-05-22 | 苏州大学 | Method and system for generating high-quality quasi-Bessel array beam |
CN107741691B (en) * | 2017-11-21 | 2019-07-26 | 清华大学 | Quasi-common optical path multi-angle illumination digital holographic imaging system and method |
CN108923880B (en) * | 2018-04-23 | 2019-06-07 | 中山大学 | A kind of photon trajectory angular momentum mode measurement method and system based on Screw transformation |
CN109242753B (en) * | 2018-08-22 | 2020-09-29 | 中国科学院西安光学精密机械研究所 | Method and device for stochastic resonance reconstruction of low-light images based on cross-modulation instability |
CN109782432A (en) * | 2019-01-17 | 2019-05-21 | 哈尔滨工业大学 | Single-shot stacked phase recovery technology of converging lens based on spatial light modulator |
CN109886882A (en) * | 2019-01-17 | 2019-06-14 | 哈尔滨工业大学 | Color single-shot stacked phase recovery technology based on spatial light modulator |
CN109785266A (en) * | 2019-01-17 | 2019-05-21 | 哈尔滨工业大学 | Single-shot stacked phase recovery technology based on spatial light modulator |
CN110361098B (en) * | 2019-08-26 | 2020-05-29 | 苏州大学 | A Method for Realizing Complex Coherence Measurement of Random Light Field Using Light Intensity Correlation |
CN110702383B (en) * | 2019-10-10 | 2020-11-20 | 浙江大学 | A wavefront reconstruction device and method based on extended rotationally symmetric structured light illumination |
CN111028974B (en) * | 2019-12-06 | 2023-04-07 | 上海理工大学 | Novel method and system for classifying and screening micro particles by array vortex light beams |
CN114895542B (en) * | 2022-04-27 | 2023-11-21 | 安徽大学 | Non-iterative rapid generation method of three-dimensional calculation hologram |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105677A (en) * | 2013-02-27 | 2013-05-15 | 苏州大学 | System and method used for producing Laguerre-Gaussian correlated partial coherence gauss beam |
CN103592768A (en) * | 2013-11-27 | 2014-02-19 | 苏州大学 | Cosine-gauss correlative light beam generation system and method and measuring equipment of cosine-gauss correlative light beam |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6525875B1 (en) * | 1998-04-15 | 2003-02-25 | Vincent Lauer | Microscope generating a three-dimensional representation of an object and images generated by such a microscope |
US7649160B2 (en) * | 2005-02-23 | 2010-01-19 | Lyncee Tec S.A. | Wave front sensing method and apparatus |
CN100573065C (en) * | 2006-08-07 | 2009-12-23 | 中国科学院光电技术研究所 | A Wavefront Measurement Method Based on Linear Phase Inversion |
FR2915281B1 (en) * | 2007-04-20 | 2010-05-14 | Commissariat Energie Atomique | METHOD FOR DETERMINING A DISTURBANCE OF AN OPTICAL WAVE |
US8040595B2 (en) * | 2007-11-02 | 2011-10-18 | Wavefront Analysis, Inc. | Light microscope with novel digital method to achieve super-resolution |
CN201434340Y (en) * | 2008-11-10 | 2010-03-31 | 长春理工大学 | Near-Field Uniform Illumination Device with Diffraction Hole Array |
CN101482503B (en) * | 2009-02-20 | 2011-06-29 | 山东师范大学 | A coherent diffraction imaging method and processing device thereof |
GB2501112B (en) * | 2012-04-12 | 2014-04-16 | Two Trees Photonics Ltd | Phase retrieval |
GB201215558D0 (en) * | 2012-08-31 | 2012-10-17 | Phase Focus Ltd | Improvements in phase retrieval |
CN102878930B (en) * | 2012-10-19 | 2015-08-19 | 华南师范大学 | A kind of method for quantitative measuring of phase object phase mehtod and device and application thereof |
CN202975600U (en) * | 2012-11-23 | 2013-06-05 | 苏州大学 | Measuring device of partially coherent vortex light beams |
CN102981391A (en) * | 2012-12-18 | 2013-03-20 | 苏州大学 | Method for generating random electromagnetic gaussian schell model vortex beam |
CN103063162B (en) * | 2013-01-11 | 2015-09-23 | 苏州大学 | A kind of method of measure portion coherent Gaussian beam Wave-front phase radius |
CN203163700U (en) * | 2013-01-11 | 2013-08-28 | 苏州大学 | Apparatus for measuring partially coherent Gaussian light beam wavefront phase radius |
CN103941405A (en) * | 2014-04-11 | 2014-07-23 | 苏州大学 | Method and device for generating perfect Laguerre-Gaussian beam |
CN103941407B (en) * | 2014-05-12 | 2016-08-24 | 苏州大学 | Generation system, production method and the measurement apparatus of partially coherent multimode Gaussian beam |
CN104331857A (en) * | 2014-10-28 | 2015-02-04 | 南京理工大学 | Phase position difference iteration compensation method in light intensity transmission equation phase retrieval |
CN207281476U (en) * | 2017-03-29 | 2018-04-27 | 苏州大学 | A kind of non-iterative phase recovery device under partially coherent optical illumination |
-
2017
- 2017-03-29 CN CN201710197745.5A patent/CN106990694B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105677A (en) * | 2013-02-27 | 2013-05-15 | 苏州大学 | System and method used for producing Laguerre-Gaussian correlated partial coherence gauss beam |
CN103592768A (en) * | 2013-11-27 | 2014-02-19 | 苏州大学 | Cosine-gauss correlative light beam generation system and method and measuring equipment of cosine-gauss correlative light beam |
Also Published As
Publication number | Publication date |
---|---|
CN106990694A (en) | 2017-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106990694B (en) | A non-iterative phase recovery device and method under partially coherent light illumination | |
CN101726366B (en) | A wavefront measurement method and device based on a multi-pinhole plate | |
Palacios et al. | 3D image reconstruction of transparent microscopic objects using digital holography | |
US12211241B2 (en) | Method and apparatus for super-resolution optical imaging | |
CN110455834B (en) | X-ray single exposure imaging device and method based on light intensity transmission equation | |
JP2011508896A (en) | Optical microscope with a novel digital method for achieving super-resolution | |
Flores-Moreno et al. | DHI contemporary methodologies: a review and frontiers | |
WO2023109428A1 (en) | All-light-field imaging camera and imaging method thereof, and all-light-field imaging device | |
CN109343321B (en) | X-ray single exposure phase shift radial shear digital holographic imaging method | |
CN109283821B (en) | Phase-shift digital holographic single-exposure imaging device and method based on vortex lens | |
CN207281476U (en) | A kind of non-iterative phase recovery device under partially coherent optical illumination | |
CN104793475A (en) | Non-coherent telescoping digital holographic imaging method and assorted device | |
Almoro et al. | Object wave reconstruction by speckle illumination and phase retrieval | |
CN112702486A (en) | A high-speed coherent imaging camera and phase recovery method based on electrically tunable lens | |
CN111553971B (en) | Fast full-field complex-number domain imaging method and device based on triple binarization modulation | |
Qiao et al. | Snapshot interferometric 3D imaging by compressive sensing and deep learning | |
He et al. | On-line beam diagnostics based on single-shot beam splitting phase retrieval | |
CN111982014A (en) | Micro-interference-based microsphere surface morphology large-field-of-view measurement method | |
Yao et al. | Phase retrieval based on coded splitting modulation | |
CN110132901A (en) | Systems and methods for synthetic aperture imaging through scattering media | |
CN102506755A (en) | Digital holographic tomography image recording device based on directional frequency spectrum separation | |
Shevkunov et al. | CNN-assisted quantitative phase microscopy for biological cell imaging. | |
Xu et al. | Sparse scanning Hartmann wavefront sensor | |
Baker | Curvature wave-front sensors for electron density characterization in plasmas | |
Stępień et al. | Block-matching-based filtration in holographic tomography reconstruction |
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 |