CN201421432Y - A coherent diffraction imaging processing device - Google Patents
A coherent diffraction imaging processing device Download PDFInfo
- Publication number
- CN201421432Y CN201421432Y CN2009200194783U CN200920019478U CN201421432Y CN 201421432 Y CN201421432 Y CN 201421432Y CN 2009200194783 U CN2009200194783 U CN 2009200194783U CN 200920019478 U CN200920019478 U CN 200920019478U CN 201421432 Y CN201421432 Y CN 201421432Y
- Authority
- CN
- China
- Prior art keywords
- pinhole
- plate
- pinholes
- processing device
- pinhole plate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 50
- 238000012545 processing Methods 0.000 title claims abstract description 40
- 230000001427 coherent effect Effects 0.000 title claims abstract description 37
- 238000009826 distribution Methods 0.000 claims description 25
- 238000005259 measurement Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- 238000000034 method Methods 0.000 description 23
- 238000005516 engineering process Methods 0.000 description 8
- 238000005314 correlation function Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000007493 shaping process Methods 0.000 description 5
- 238000011084 recovery Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012804 iterative process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000005311 autocorrelation function Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本实用新型公开了一种相干衍射成像处理装置。它包括光源、放置样品的样品台和图像记录及处理装置,其特征在于,它还包括扩束及波面整形器、多针孔板和步进旋转支架;多针孔板安装在步进旋转支架上,多针孔板上的参考针孔中心与步进旋转支架的旋转中心重合;沿光束前进的方向依次排列光源、扩束及波面整形器、样品、多针孔板和图像记录及处理装置;图像记录及处理装置包括图像传感器,图像传感器置于步进旋转支架之后,并与计算机连接;同时计算机还与步进旋转支架连接。它结构简单,调节方便,成本较低,适用于多种不同光源,可实现复数物体或三维物体的成像,不需要成像透镜,尤其适用于像X射线这类难以制备出高质量成像透镜的场合。
The utility model discloses a coherent diffraction imaging processing device. It includes a light source, a sample stage for placing samples, and an image recording and processing device. It is characterized in that it also includes a beam expander and wave front shaper, a multi-pinhole plate and a stepping rotation bracket; the multi-pinhole plate is installed on the stepping rotation bracket Above, the center of the reference pinhole on the multi-pinhole plate coincides with the rotation center of the stepping rotating bracket; the light source, beam expander and wave surface shaper, sample, multi-pinhole plate and image recording and processing device are arranged in sequence along the direction of beam advancement The image recording and processing device includes an image sensor, which is placed behind the stepping and rotating bracket and connected with the computer; meanwhile, the computer is also connected with the stepping and rotating bracket. It is simple in structure, easy to adjust, and low in cost. It is suitable for a variety of different light sources, and can realize the imaging of multiple objects or three-dimensional objects without imaging lenses. It is especially suitable for occasions such as X-rays where it is difficult to prepare high-quality imaging lenses. .
Description
技术领域 technical field
本实用新型涉及一种相干衍射成像技术,特别是一种利用多针孔板实现一般复振幅物体记录和重现的相干衍射成像处理装置。The utility model relates to a coherent diffraction imaging technology, in particular to a coherent diffraction imaging processing device which uses a multi-pinhole plate to realize the recording and reproduction of general complex amplitude objects.
背景技术 Background technique
相干衍射成像是一种利用物波的远场或夫琅和费衍射的强度分布实现复数或三维物体成像的技术。由于可以避免成像透镜孔径和像差对分辨率的限制,这种技术特别适用于像X射线、电子束这类缺乏或很难制备出高质量成像透镜的领域。相干衍射成像技术的关键问题是怎样从一幅或多幅衍射强度图样中准确快速地将被测物波的振幅和相位信息恢复出来。目前,解决这一问题的一个主要途径是采用迭代算法。传统的迭代方法一般都需要较长的迭代时间,迭代结果存在不确定性;对被测物体也存在一些苛刻的限制条件,如物体要特殊固定、要求被测物体是纯振幅物体或纯相位物体,等等。Coherent diffraction imaging is a technology that utilizes the far field of object waves or the intensity distribution of Fraunhofer diffraction to realize imaging of complex or three-dimensional objects. Because it can avoid the limitations of imaging lens aperture and aberration on resolution, this technology is especially suitable for fields such as X-rays and electron beams that lack or are difficult to manufacture high-quality imaging lenses. The key problem of coherent diffraction imaging technology is how to accurately and quickly recover the amplitude and phase information of the measured object wave from one or more diffraction intensity patterns. At present, one of the main ways to solve this problem is to use iterative algorithm. Traditional iterative methods generally require a long iteration time, and there are uncertainties in the iterative results; there are also some harsh restrictions on the measured object, such as the object must be specially fixed, and the measured object is required to be a pure amplitude object or a pure phase object ,etc.
为了提高相位恢复的准确性和成像速度,近年来人们已经发展了一些改进的衍射成像技术,如超抽样衍射成像技术、分区扫描衍射成像技术等(参见Phys.Rev.Lett.93,023903(2004);Phys.Rev.Lett.98,034801(2007);Phys.Rev.Lett.100,155503(2008);Science 321,379(2008))。但上述技术都不能从根本上避免迭代算法和迭代结果的不确定性。其相位恢复过程都还需要涉及大量的迭代过程;这就意味着用这些方法实现动态或实时成像仍然是比较困难的。In order to improve the accuracy and imaging speed of phase recovery, some improved diffraction imaging techniques have been developed in recent years, such as super-sampling diffraction imaging technology, partition scanning diffraction imaging technology, etc. (see Phys.Rev.Lett.93, 023903 (2004 ); Phys. Rev. Lett. 98, 034801 (2007); Phys. Rev. Lett. 100, 155503 (2008); Science 321, 379 (2008)). But none of the above techniques can fundamentally avoid the uncertainty of iterative algorithms and iterative results. The phase recovery process still needs to involve a large number of iterative processes; this means that it is still relatively difficult to achieve dynamic or real-time imaging with these methods.
实用新型内容 Utility model content
本实用新型的目的就是为了克服现有技术所存在的问题和不足,提供一种相干衍射成像处理装置,它是让被成像物体所发出的物波先通过一个可旋转的多针孔板,然后再用图像传感器记录物波经过多针孔板后的夫琅和费衍射强度分布,最后借助计算机图像处理技术提取出被测物波的振幅和相位分布信息。The purpose of the utility model is to overcome the existing problems and deficiencies in the prior art, and provide a coherent diffraction imaging processing device, which allows the object wave emitted by the imaged object to pass through a rotatable multi-pinhole plate first, and then The image sensor is used to record the Fraunhofer diffraction intensity distribution of the object wave passing through the multi-pinhole plate, and finally the amplitude and phase distribution information of the measured object wave is extracted by means of computer image processing technology.
为实现上述目的,本实用新型采用如下技术方案:In order to achieve the above object, the utility model adopts the following technical solutions:
一种相干衍射成像处理装置,它包括光源、放置样品的样品台和图像记录及处理装置,其特征在于,它还包括扩束及波面整形器、多针孔板和步进旋转支架;多针孔板安装在步进旋转支架上,多针孔板上的参考针孔中心与步进旋转支架的旋转中心重合;沿光束前进的方向依次排列光源、扩束及波面整形器、样品、多针孔板和图像记录及处理装置;图像记录及处理装置包括图像传感器,图像传感器置于步进旋转支架之后,并与计算机连接;同时计算机还与步进旋转支架连接。A coherent diffraction imaging processing device, which includes a light source, a sample stage for placing samples, and an image recording and processing device, and is characterized in that it also includes a beam expander and wavefront shaper, a multi-pinhole plate and a stepping rotary support; a multi-pin The orifice plate is installed on the stepping rotating bracket, and the center of the reference pinhole on the multi-pinhole plate coincides with the rotation center of the stepping rotating bracket; the light source, beam expander and wave surface shaper, sample, multi-pin are arranged in sequence along the direction of beam advance The orifice plate and the image recording and processing device; the image recording and processing device includes an image sensor, the image sensor is placed behind the stepping and rotating bracket, and is connected with a computer; meanwhile, the computer is also connected with the stepping and rotating bracket.
所述样品与多针孔板之间的距离Z1保证照射到多针孔板上的物波在每个针孔通光孔径上的相位变化不大于π/4弧度;所述图像传感器与多针孔板之间的距离Z2满足夫琅和费衍射条件,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。The distance Z1 between the sample and the multi-pinhole plate ensures that the phase change of the object wave irradiated on the multi-pinhole plate on the clear aperture of each pinhole is not greater than π/4 radians; the image sensor and the multi-pinhole plate The distance Z2 between them satisfies the Fraunhofer diffraction condition, that is, the light field distribution on the image sensor recording plane is proportional to the Fourier transform of the object wave passing through the multi-pinhole plate.
一种相干衍射成像处理装置,它包括光源、放置样品的样品台和图像记录及处理装置,其特征在于,它还包括扩束及波面整形器、多针孔板和步进旋转支架以及傅里叶变换透镜,多针孔板安装在步进旋转支架上,多针孔板上的参考针孔中心与步进旋转支架的旋转中心重合;沿光束前进的方向依次排列光源、扩束及波面整形器、样品、多针孔板、傅里叶变换透镜和图像记录及处理装置;图像记录及处理装置包括图像传感器,图像传感器置于步进旋转支架之后,并与计算机连接;同时计算机还与步进旋转支架连接。A coherent diffraction imaging processing device, which includes a light source, a sample stage for placing samples, and an image recording and processing device, and is characterized in that it also includes a beam expander and wavefront shaper, a multi-pinhole plate, a stepping and rotating bracket, and a Fourier Leaf transformation lens, the multi-pinhole plate is installed on the stepping rotation bracket, the center of the reference pinhole on the multi-pinhole plate coincides with the rotation center of the stepping rotation bracket; the light source, beam expansion and wave surface shaping are arranged in sequence along the direction of the beam advance device, sample, multi-pinhole plate, Fourier transform lens, and image recording and processing device; the image recording and processing device includes an image sensor, which is placed behind the stepping and rotating bracket and connected to the computer; at the same time, the computer is also connected to the stepper into the swivel bracket connection.
所述样品与多针孔板之间的距离Z1保证照射到多针孔板上的物波在每个针孔通光孔径上的相位变化不大于π/4弧度;所述图像传感器位于傅里叶变换透镜的后焦面上,即图像传感器记录平面上的光场分布正比于透过多针孔板的物波的傅立叶变换。The distance Z1 between the sample and the multi-pinhole plate ensures that the phase change of the object wave irradiated on the multi-pinhole plate on the clear aperture of each pinhole is not greater than π/4 radians; the image sensor is located in the Fourier transform The light field distribution on the back focal plane of the lens, that is, the recording plane of the image sensor, is proportional to the Fourier transform of the object wave passing through the multi-pinhole plate.
所述多针孔板是一块不透明薄板上设置一个参考针孔和若干测量针孔组成;各针孔在多针孔板上的分布满足以下要求:两块完全相同的多针孔板重叠放置,当平移其中一块使其参考针孔和另一块的任一测量针孔重叠时,其它针孔都不会发生重叠;针孔的大小满足一个针孔范围内物波的相位变化不大于四分之一波长。The multi-pinhole plate is composed of a reference pinhole and several measuring pinholes arranged on an opaque thin plate; the distribution of each pinhole on the multi-pinhole plate meets the following requirements: two identical multi-pinhole plates are placed overlappingly, When one of the pinholes is translated so that the reference pinhole overlaps with any measurement pinhole of the other piece, the other pinholes will not overlap; the size of the pinhole satisfies that the phase change of the object wave within the range of a pinhole is not more than 1/4 one wavelength.
所述多针孔板上的针孔均匀分布在一个半圆环上,针孔中心的间隔相等,不小于针孔直径的1.5倍,位于针孔板中心半圆环一端处的针孔为参考针孔,其他针孔均为测量针孔。The pinholes on the multi-pinhole plate are evenly distributed on a semicircular ring, the center of the pinholes are equally spaced, not less than 1.5 times the diameter of the pinhole, and the pinholes at one end of the semicircular ring in the center of the pinhole plate are used as a reference pinhole, other pinholes are measuring pinholes.
所述多针孔板上的针孔分布在三条互成120度的直线上,其中位于针孔板中心三直线相交处的针孔为参考针孔,其他针孔均为测量针孔;参考针孔中心与相邻的测量针孔中心的距离D0不小于针孔直径的1.5倍,其他测量针孔中心的间隔相等,均为D0的两倍。The pinholes on the multi-pinhole plate are distributed on three straight lines at 120 degrees to each other, wherein the pinhole located at the intersection of the three straight lines in the center of the pinhole plate is a reference pinhole, and the other pinholes are all measuring pinholes; the reference needle The distance D0 between the center of the hole and the center of the adjacent measuring pinhole is not less than 1.5 times the diameter of the pinhole, and the intervals between the centers of other measuring pinholes are equal, which are twice as long as D0.
本方法完全不需要迭代过程。相位恢复过程也只在一幅衍射强度图像中进行,不涉及多幅图像的交互处理,因此相位恢复过程可以和衍射图样的扫描记录过程同步进行。这样不仅可以大大减小对扫描记录过程及定位精度的要求,而且也可以极大的提高衍射成像速度,从而为真正实现实时相干衍射成像提供了可能。This method does not require an iterative process at all. The phase recovery process is only performed in one diffraction intensity image, and does not involve the interactive processing of multiple images, so the phase recovery process can be performed simultaneously with the scanning and recording process of the diffraction pattern. This can not only greatly reduce the requirements on the scanning recording process and positioning accuracy, but also greatly increase the diffraction imaging speed, thus providing the possibility for real-time coherent diffraction imaging.
为了实现本实用新型,必须采用专门设计的衍射成像处理装置。本实用新型的衍射成像处理装置包括光源、扩束及波面整形器、样品台、被成像物体、多针孔板、步进旋转架、图像传感器和计算机。沿光束前进的方向依次共轴(光轴)排列着光源、扩束及波面整形器、被成像物体、多针孔板和图像传感器。被成像物体放在样品台上。多针孔板固定在可精密步进旋转的步进旋转架上。图像传感器和步进旋转架与计算机相连,通过计算机控制旋转架的定量转动和图像传感器的同步记录。相邻部件之间的距离可调。In order to realize the utility model, a specially designed diffraction imaging processing device must be used. The diffraction imaging processing device of the utility model comprises a light source, a beam expander and a wave surface shaper, a sample stage, an object to be imaged, a multi-pinhole plate, a stepping rotating frame, an image sensor and a computer. A light source, a beam expander and a wavefront shaper, an imaged object, a multi-pinhole plate and an image sensor are arranged coaxially (optical axis) in sequence along the light beam advancing direction. The object to be imaged is placed on the sample stage. The multi-pinhole plate is fixed on a step-rotation frame that can rotate in precise steps. The image sensor and the stepping rotating frame are connected with the computer, and the quantitative rotation of the rotating frame and the synchronous recording of the image sensor are controlled by the computer. The distance between adjacent parts is adjustable.
本实用新型的方法为:The method of the present utility model is:
(1)由光源发出一单色光,经扩束及波面整形处理后照明被成像物体。透过物体的光波(即物波)传播距离Z1后照射到一个特别设计的多针孔板上。该多针孔板上有一个参考针孔和多个测量针孔。(1) A monochromatic light is emitted from the light source, and the object to be imaged is illuminated after beam expansion and wavefront shaping. The light wave (that is, the object wave) passing through the object travels a distance Z1 and then irradiates a specially designed multi-pinhole plate. This multi-pinhole board has a reference pinhole and multiple measurement pinholes.
(2)透过多针孔板的物波继续传播距离Z2到达记录平面,并在该记录平面上形成它的夫琅和费衍射光场。在该记录平面处用一个二维图像传感器(如CCD)记录该夫琅和费衍射光场的强度分布图样。(2) The object wave passing through the multi-pinhole plate continues to travel distance Z2 to reach the recording plane, and forms its Fraunhofer diffraction light field on the recording plane. A two-dimensional image sensor (such as a CCD) is used to record the intensity distribution pattern of the Fraunhofer diffraction light field at the recording plane.
(3)在计算机中对该夫琅和费衍射强度图样做逆傅里叶变换得到透过多针孔板的物波复振幅的相关函数图样(一般为复变函数)。将该相关函数图样中与多针孔板上各测量针孔的中心位置所对应的点的函数值提取出来,就得到被测物波在各测量针孔处的相对振幅和相位值。(3) Perform an inverse Fourier transform on the Fraunhofer diffraction intensity pattern in the computer to obtain a correlation function pattern (generally a complex variable function) of the complex amplitude of the object wave passing through the multi-pinhole plate. By extracting the function value of the point corresponding to the center position of each measuring pinhole on the multi-pinhole plate in the correlation function pattern, the relative amplitude and phase value of the measured object wave at each measuring pinhole can be obtained.
(4)以参考针孔的中心为轴旋转多针孔板,使测量针孔在多针孔板所在平面上扫描。记录不同旋转角度下透过多针孔板的物波的夫琅和费衍射强度图样,并对记录的每幅衍射强度图样重复步骤(3),就得到物波在多针孔板所在平面上的复振幅分布的二维抽样阵列。利用该抽样阵列就可以在计算机中重现被成像物体。(4) Rotate the multi-pinhole plate with the center of the reference pinhole as the axis, so that the measuring pinhole scans on the plane where the multi-pinhole plate is located. Record the Fraunhofer diffraction intensity patterns of the object wave passing through the multi-pinhole plate at different rotation angles, and repeat step (3) for each recorded diffraction intensity pattern, to obtain the object wave on the plane where the multi-pinhole plate is located A 2D sampled array of the complex amplitude distribution of . Using this sampling array, the imaged object can be reproduced in a computer.
上述方法中,步骤(1)中的光源可以是可见光源,也可以是其他相干波源,如紫外光、X射线、电子束等。扩束及波面整形处理单元主要包括扩束器、限制光阑、准直器或会聚器,它们沿光束传播方向依次共轴排列,目的是对从光源发出的光波进行波面整形,以产生后续光路所需要的相干平面波或球面波。步骤(1)中所用多针孔板是通过在一块不透明薄板上制备若干微小针孔形成的。该多针孔板上的针孔包含一个参考针孔和若干个测量针孔。它们在多针孔板上的排列方式满足以下要求:若将两块完全相同的多针孔板重叠放置,当平移其中一块使其参考针孔和另一块的任一测量针孔重叠时,其它针孔都不会发生重叠。多针孔板上的针孔形状可以是圆孔也可以是其它形状的通光孔。针孔的大小取决于传播到多针孔板处的物波的复振幅分布的空间变化情况,一般要求在一个针孔范围内物波的相位变化不大于π/4弧度。当针孔大小一定时,可通过适当选择多针孔板与被成像物体的距离Z1使满足上述条件。In the above method, the light source in step (1) can be a visible light source, or other coherent wave sources, such as ultraviolet light, X-rays, electron beams, etc. The beam expansion and wave surface shaping processing unit mainly includes a beam expander, a limiting diaphragm, a collimator or a converging device, which are arranged coaxially in sequence along the beam propagation direction, and the purpose is to perform wave front shaping on the light wave emitted from the light source to generate a subsequent optical path Coherent plane or spherical waves required. The multi-pinhole plate used in step (1) is formed by preparing several tiny pinholes on an opaque thin plate. The pinholes on the multi-pinhole board include a reference pinhole and several measurement pinholes. Their arrangement on the multi-pinhole board meets the following requirements: If two identical multi-pinhole boards are placed on top of each other, when one of them is translated so that the reference pinhole overlaps with any measurement pinhole of the other, the other None of the pinholes overlap. The shape of the pinholes on the multi-pinhole board can be a round hole or a light-through hole of other shapes. The size of the pinhole depends on the spatial variation of the complex amplitude distribution of the object wave propagating to the multi-pinhole plate, and it is generally required that the phase change of the object wave within a pinhole range is not greater than π/4 radians. When the size of the pinholes is constant, the above conditions can be met by properly selecting the distance Z1 between the multi-pinhole plate and the object to be imaged.
步骤(2)中的记录平面与多针孔板之间的距离Z2满足夫琅和费衍射条件,即记录平面上的衍射光场正比于多针孔板平面上透过多针孔板的物波的傅立叶变换。设透过多针孔板上第m个针孔的物波的复振幅值为其中Am为第m个针孔处的物波的振幅,为第m个针孔处的物波的相位,为针孔板所在平面上的坐标变量(以参考针孔的中心为坐标原点),为第m个针孔中心的坐标,j为虚数,为第m个针孔的孔径函数;在记录平面上得到的衍射光场的强度分布可表示为:The distance Z2 between the recording plane and the multi-pinhole plate in step (2) satisfies the Fraunhofer diffraction condition, that is, the diffracted light field on the recording plane is proportional to the object passing through the multi-pinhole plate on the plane of the multi-pinhole plate Fourier transform of a wave. Let the complex amplitude of the object wave passing through the mth pinhole on the multi-pinhole plate be where A m is the amplitude of the object wave at the mth pinhole, is the phase of the object wave at the mth pinhole, is the coordinate variable on the plane where the pinhole plate is located (with the center of the reference pinhole as the coordinate origin), is the coordinate of the center of the mth pinhole, j is an imaginary number, is the aperture function of the mth pinhole; the intensity distribution of the diffracted light field obtained on the recording plane can be expressed as:
其中,I0为积分常数,为傅立叶变换算符,N为针孔板上测量针孔的数目,m=0对应参考针孔,为记录平面上的坐标变量。Among them, I 0 is the integral constant, is the Fourier transform operator, N is the number of measuring pinholes on the pinhole plate, m=0 corresponds to the reference pinhole, is the coordinate variable on the recording plane.
步骤(3)中对衍射强度图样进行逆傅立叶变换可以通过将该衍射强度图样读入计算中用计算机程序实现,也可以通过专用DSP芯片来完成。由于(2)式描述的物波的夫琅和费衍射的强度分布正比于物波的傅立叶变换函数的模方,对其逆傅立叶变换的结果刚好是该物波的相关函数。一般情况下,从该相关函数中恢复物波的复振幅函数是一项非常困难的任务,要用迭代算法才行。本方法中,由于所记录的衍射强度分布是通过前述的多针孔板后的物波的夫琅和费衍射强度,该相关函数在与任意测量针孔中心位置相对应的点的值刚好正比于透过该针孔的物波复振幅与透过参考针孔的物波的复共轭振幅值的乘积(以参考针孔中心为坐标中心)。因此,本方法不需要任何迭代算法。透过多针孔板的物波复振幅可以直接从衍射强度图样的逆傅立叶变换图样中提取。具体计算过程如下:对(1)式作逆傅立叶变换,得到的相关函数为The inverse Fourier transform of the diffraction intensity pattern in step (3) can be realized by reading the diffraction intensity pattern into the calculation with a computer program, or by a dedicated DSP chip. Since the intensity distribution of the Fraunhofer diffraction of the object wave described by formula (2) is proportional to the modulus of the Fourier transform function of the object wave, the result of its inverse Fourier transform is just the correlation function of the object wave. In general, it is a very difficult task to recover the complex amplitude function of the object wave from this correlation function, and an iterative algorithm is required. In this method, since the recorded diffraction intensity distribution is the Fraunhofer diffraction intensity of the object wave passing through the aforementioned multi-pinhole plate, the value of the correlation function at the point corresponding to the center position of any measuring pinhole is just proportional to The product of the complex amplitude of the object wave passing through the pinhole and the complex conjugate amplitude of the object wave passing through the reference pinhole (with the center of the reference pinhole as the coordinate center). Therefore, the method does not require any iterative algorithm. The complex amplitude of the object wave transmitted through the multi-pinhole plate can be directly extracted from the inverse Fourier transform pattern of the diffraction intensity pattern. The specific calculation process is as follows: Inverse Fourier transform is performed on formula (1), and the obtained correlation function for
其中,An为第n个针孔处的物波的振幅,为第n个针孔处的物波的相位,Among them, A n is the amplitude of the object wave at the nth pinhole, is the phase of the object wave at the nth pinhole,
其中,为积分变量,*号为共轭符号。根据(3)式所示的相关运算的几何意义,并考虑到前述针孔板上的针孔分布特点,可知(2)式在坐标矢量
其中P0为正比于针孔面积的常数。(4)式表明,由(2)式给出的逆傅立叶变换图样中与各测量针孔中心点相对应的位置上的值刚好就是我们想要得到的透过该针孔的物波的相对复振幅值。where P 0 is a constant proportional to the area of the pinhole. Equation (4) shows that the value at the position corresponding to the center point of each measuring pinhole in the inverse Fourier transform pattern given by Equation (2) is just the relative value of the object wave passing through the pinhole we want to obtain. Complex amplitude value.
步骤(4)中的可旋转多针孔板是以垂直通过参考针孔中心的直线为旋转轴。这样,在多针孔板的旋转过程中,参考针孔的位置始终保持不变,从而保证了在不同旋转角度下测得的物波复振幅分布的参考点保持一致。The rotatable multi-pinhole plate in step (4) takes a straight line passing through the center of the reference pinhole vertically as the axis of rotation. In this way, during the rotation process of the multi-pinhole plate, the position of the reference pinhole remains unchanged, thereby ensuring that the reference point of the object wave complex amplitude distribution measured at different rotation angles remains consistent.
本实用新型的原理是:被成像物体发出的物波先经过一个可旋转的多针孔板抽样,然后由一个图像传感器记录下透过多针孔板的物波的夫琅和费衍射的强度分布图样,接着对该衍射强度分布做逆傅立叶变换得到被记录物波的自相关函数,从该自相关函数中与各测量针孔中心位置坐标对应的点就可直接提取出被测物波的振幅和相位信息从而实现复振幅物体的衍射成像。The principle of the utility model is: the object wave emitted by the imaged object first passes through a rotatable multi-pinhole plate for sampling, and then an image sensor records the intensity of the Fraunhofer diffraction of the object wave passing through the multi-pinhole plate Then do inverse Fourier transform on the diffraction intensity distribution to obtain the autocorrelation function of the recorded object wave, and from the autocorrelation function, the points corresponding to the center position coordinates of each measuring pinhole can directly extract the measured object wave Amplitude and phase information to achieve diffraction imaging of complex amplitude objects.
本实用新型具有以下优点和有益效果:The utility model has the following advantages and beneficial effects:
(1)不需要成像透镜。物波的振幅和相位信息是从物波的夫琅和费衍射强度图样中通过简单的算法恢复出来的。(1) No imaging lens is required. The amplitude and phase information of the object wave is recovered from the Fraunhofer diffraction intensity pattern of the object wave by a simple algorithm.
(2)从物波的夫琅和费衍射强度图样中提取物波的振幅和相位信息的过程中不需要任何迭代过程,并且该振幅和相位信息是从该夫琅和费衍射强度的逆傅里叶变换图样中直接提取出来的。(2) The process of extracting the amplitude and phase information of the object wave from the Fraunhofer diffraction intensity pattern of the object wave does not require any iterative process, and the amplitude and phase information are obtained from the inverse Fu of the Fraunhofer diffraction intensity It is directly extracted from the leaf transform pattern.
(3)振幅和相位信息的提取过程可以和衍射图样的记录过程同步进行,效率高,速度快,为动态相干衍射成像的实现提供了一条可行的途径。(3) The extraction process of amplitude and phase information can be carried out synchronously with the recording process of the diffraction pattern, with high efficiency and fast speed, which provides a feasible way for the realization of dynamic coherent diffraction imaging.
(4)振幅和相位信息的提取不涉及多幅图样之间的交互处理,这样可以大大降低对记录系统的稳定性和扫描精度的要求。(4) The extraction of amplitude and phase information does not involve interactive processing between multiple patterns, which can greatly reduce the requirements for the stability and scanning accuracy of the recording system.
(5)本实用新型所采用的多针孔板容易制备,并且特别适用于像X射线这样的难以使用传统成像技术进行成像的波段。(5) The multi-pinhole plate used in the present invention is easy to prepare, and is especially suitable for the wave bands that are difficult to be imaged by traditional imaging techniques such as X-rays.
附图说明 Description of drawings
图1是本实用新型的相干衍射成像方法的过程方框图;Fig. 1 is the process block diagram of coherent diffraction imaging method of the present utility model;
图2是本实用新型的相干衍射成像方法中的多针孔板的第一种设计实例的针孔分布示意图;Fig. 2 is the pinhole distribution schematic diagram of the first design example of the multi-pinhole plate in the coherent diffraction imaging method of the present invention;
图3是本实用新型的多针孔板的第二种设计实例的针孔分布示意图;Fig. 3 is a schematic diagram of pinhole distribution of the second design example of the multi-pinhole plate of the present invention;
图4是本实用新型的相干衍射成像处理装置的第一种实施方式的结构示意图;Fig. 4 is a schematic structural view of the first embodiment of the coherent diffraction imaging processing device of the present invention;
图5是本实用新型的相干衍射成像处理装置的第二种实施方式的结构示意图;5 is a schematic structural view of a second embodiment of the coherent diffraction imaging processing device of the present invention;
图6a是采用本实用新型的相干衍射成像方法和装置进行相干衍射成像的一个实验结果实例;Figure 6a is an example of an experimental result of coherent diffraction imaging using the coherent diffraction imaging method and device of the present invention;
图6b是采用本实用新型的相干衍射成像方法和装置进行相干衍射成像的一个实验结果实例;Figure 6b is an example of an experimental result of coherent diffraction imaging using the coherent diffraction imaging method and device of the present invention;
图6c是采用本实用新型的相干衍射成像方法和装置进行相干衍射成像的一个实验结果实例;Figure 6c is an example of an experimental result of coherent diffraction imaging using the coherent diffraction imaging method and device of the present invention;
图6d是采用本实用新型的相干衍射成像方法和装置进行相干衍射成像的一个实验结果实例。Fig. 6d is an example of an experimental result of coherent diffraction imaging using the coherent diffraction imaging method and device of the present invention.
其中,光源1、扩束及波面整形器2、样品台3、样品4、步进旋转支架5、多针孔板6、图像传感器7、计算机8、傅里叶变换透镜9。Among them,
具体实施方式 Detailed ways
为了更好地理解本实用新型,下面结合附图与实施例对本实用新型做进一步的描述。In order to better understand the utility model, the utility model will be further described below in conjunction with the accompanying drawings and embodiments.
如图1所示,本实用新型方法包括如下步骤:(1)由光源发出的一束单色光,经扩束及波面整形处理单元处理形成一束相干平面波(或汇聚球面光);用该相干平面波(或球面波)照明置于样品台上的被成像物体;(2)透过物体的物波传播一定距离Z1后照射到一个固定在旋转支架上的具有特殊针孔分布的多针孔板上,该多针孔板上的针孔由一个参考针孔和若干个测量针孔构成;(3)透过多针孔板的物波继续传播,在远场(或通过一个傅里叶变换透镜后)形成夫琅和费衍射光场;(4)用一个二维图像传感器(如CCD)记录透过多针孔板的物波的夫琅和费衍射强度图样;(5)利用图像处理技术对该夫琅和费衍射强度图样作逆傅里叶变换得到透过多针孔板的物波复振幅的相关函数图样;将该相关函数图样中与多针孔板上各测量针孔的中心位置所对应的点的函数值提取出来,就得到透过多针孔板的物波复振幅在各测量针孔处的相对振幅和位相值;(6)以参考针孔的中心为轴旋转多针孔板,使测量针孔在多针孔板所在平面上扫描。记录不同旋转角度下透过多针孔板的物波的夫琅和费衍射强度图样,并对记录的每幅衍射强度图样重复步骤(5),就得到物波在多针孔板所在平面上的复振幅分布的二维抽样阵列。利用该抽样阵列就可以在计算机中重现被成像物体。As shown in Figure 1, the utility model method comprises the following steps: (1) a beam of monochromatic light emitted by the light source is processed by a beam expanding and wavefront shaping processing unit to form a beam of coherent plane waves (or converging spherical light); Coherent plane wave (or spherical wave) illuminates the imaged object placed on the sample stage; (2) The object wave that passes through the object propagates a certain distance Z1 and then irradiates a multi-pinhole fixed on a rotating support with a special pinhole distribution The pinholes on the multi-pinhole board consist of a reference pinhole and several measuring pinholes; (3) the object wave passing through the multi-pinhole board continues to propagate, and in the far field (or through a Fourier transform after changing the lens) to form the Fraunhofer diffraction light field; (4) use a two-dimensional image sensor (such as CCD) to record the Fraunhofer diffraction intensity pattern of the object wave passing through the multi-pinhole plate; (5) use the image The processing technology performs inverse Fourier transform on the Fraunhofer diffraction intensity pattern to obtain the correlation function pattern of the complex amplitude of the object wave passing through the multi-pinhole plate; Extract the function value of the point corresponding to the central position of , and then get the relative amplitude and phase value of the complex amplitude of the object wave passing through the multi-pinhole plate at each measuring pinhole; (6) take the center of the reference pinhole as the axis Rotate the multi-pinhole plate to make the measurement pinholes scan on the plane where the multi-pinhole plate is located. Record the Fraunhofer diffraction intensity patterns of the object wave passing through the multi-pinhole plate at different rotation angles, and repeat step (5) for each recorded diffraction intensity pattern, to obtain the object wave on the plane where the multi-pinhole plate is located. A 2D sampled array of the complex amplitude distribution of . Using this sampling array, the imaged object can be reproduced in a computer.
图2是本实用新型的多针孔板的一种典型实施方式。该针孔板上的针孔均匀分布在一个半圆环上,针孔中心的间隔不小于针孔直径的1.5倍。位于针孔板中心半圆环一端处的针孔为参考针孔,其他针孔均为测量针孔。将该多针孔板放置在旋转支架上时,要让参考针孔的中心与旋转轴重合。Fig. 2 is a typical implementation of the multi-pinhole plate of the present invention. The pinholes on the pinhole plate are evenly distributed on a semicircular ring, and the distance between the centers of the pinholes is not less than 1.5 times the diameter of the pinholes. The pinhole located at one end of the semicircular ring in the center of the pinhole plate is a reference pinhole, and the other pinholes are measuring pinholes. Place the multi-pinhole plate on the rotation stand so that the center of the reference pinhole coincides with the axis of rotation.
图3是本实用新型的方法中所述的多针孔板的第二种实施方式。该多针孔板上的针孔分布在三条互成120度的直线上。位于针孔板中心三直线相交处的针孔为参考针孔,其他针孔均为测量针孔。参考针孔中心与相邻的测量针孔中心间的距离D0不小于针孔直径的1.5倍。其他测量针孔中心的间隔相等,均为D0的两倍。将该多针孔板放置在旋转支架上时,要让参考针孔的中心与旋转轴重合。Fig. 3 is the second embodiment of the multi-pinhole plate described in the method of the present invention. The pinholes on the multi-pinhole plate are distributed on three straight lines at 120 degrees to each other. The pinhole located at the intersection of the three straight lines in the center of the pinhole plate is the reference pinhole, and the other pinholes are measurement pinholes. The distance D0 between the center of the reference pinhole and the center of the adjacent measuring pinhole is not less than 1.5 times the diameter of the pinhole. The intervals between the centers of other measuring pinholes are equal, which are twice of D0. Place the multi-pinhole plate on the rotation stand so that the center of the reference pinhole coincides with the axis of rotation.
图4是本实用新型的相干衍射成像处理装置的一种典型实施方式。该装置包括光源1、扩束及波面整形器2、样品台3、样品4、步进旋转支架5、多针孔板6、图像传感器7、计算机8。沿光束前进的方向依次共轴(光轴)排列着光源1、扩束及波面整形器2、样品台3、样品4、步进旋转支架5、多针孔板6、图像传感器7。光源1、扩束及波面整形器2之间的相互位置保证由光源1发出的单色光经扩束及波面整形器3后形成一束平面波或球面波。样品(被成像物体)4和多针孔板6之间的相对距离Z1保证照射到多针孔板上的物波在每个针孔通光孔径上的相位变化不大于π/4弧度。而多针孔板6到图像传感器7的感光平面之间的距离Z2保证满足夫琅和费衍射条件,使透过多针孔板6的物波在图像传感器感光平面上的衍射光场正比于该物波的傅里叶变换。步进旋转支架5由计算机8控制的步进驱动器控制。步进旋转支架7和图像传感器8由同一计算机8控制实现多针孔板的旋转与衍射强度图像的同步记录。从所记录的衍射强度图样中提取物波的复振幅信息并成像由根据本实用新型方法设计的计算机程序完成。Fig. 4 is a typical implementation of the coherent diffraction imaging processing device of the present invention. The device includes a
图5是本实用新型的装置的第二种实施方式,比图4多了一个傅里叶变换透镜9,即在多针孔板6和图像传感器7之间放置傅里叶变换透镜9,图像传感器记录平面位于傅里叶变换透镜9的后焦面上。这样可以使装置紧凑,并且可以通过改变透镜9的焦距方便地放大或缩小传播到传感器记录平面上的夫琅和费衍射强度图样的大小。Fig. 5 is the second kind of embodiment of the device of the present utility model, has one more Fourier transform lens 9 than Fig. 4, promptly places Fourier transform lens 9 between
图6a-图6d是采用本实用新型的相干衍射成像方法和装置进行相干衍射成像的一个实验结果实例。实验中,光源采用He-Ne激光,输出光波的波长为0.6328微米;被成像物体为字符“DM”的缩微图片;多针孔板采用图2所示的针孔排列方式。针孔直径为28微米,针孔间隔为56微米。实验装置采用图3所示的第二种实施方式,其中傅里叶变换透镜的焦距为240毫米,图像传感器为CCD数码相机,像素数为1300×1030,像素大小为6.7微米。图6a为用CCD记录的被测物波的夫琅和费衍射强度图样的其中一幅。图6b和图6c分别为利用本实用新型的方法从所记录的夫琅和费衍射强度图样中恢复出来的物波的振幅和相位分布的二维抽样阵列。图6d为利用图6b和图6c所示的测量数据在计算机中通过数字衍射得到的被测物体的一个重现像。Figures 6a-6d are examples of experimental results of coherent diffraction imaging using the coherent diffraction imaging method and device of the present invention. In the experiment, the light source is He-Ne laser, and the wavelength of the output light wave is 0.6328 microns; the object to be imaged is a miniature picture of the character "DM"; the multi-pinhole plate adopts the pinhole arrangement shown in Figure 2. The pinhole diameter is 28 microns, and the pinhole spacing is 56 microns. The experimental device adopts the second embodiment shown in Figure 3, in which the focal length of the Fourier transform lens is 240mm, the image sensor is a CCD digital camera, the number of pixels is 1300×1030, and the pixel size is 6.7 microns. Figure 6a is one of the Fraunhofer diffraction intensity patterns of the measured object waves recorded by the CCD. Fig. 6b and Fig. 6c are the two-dimensional sampling arrays of the amplitude and phase distribution of the object wave recovered from the recorded Fraunhofer diffraction intensity pattern by using the method of the present utility model respectively. Fig. 6d is a reconstructed image of the measured object obtained by digital diffraction in a computer using the measurement data shown in Fig. 6b and Fig. 6c.
上述方法及实施例都是通过记录被测物波通过一个多针孔板后的夫琅和费衍射强度图样和通过对该图样的图像处理来达到重现物波的振幅和相位分布的目的。本实用新型的实施不局限于上述具体实施方案。只要是通过记录被测物波经过上述方法所述的多针孔板后的夫琅和费衍射强度图样来实现重现物波的振幅和相位分布目的的方法、装置与系统,均属于本实用新型的保护范围。The above methods and embodiments achieve the purpose of reproducing the amplitude and phase distribution of the object wave by recording the Fraunhofer diffraction intensity pattern after the object wave passes through a multi-pinhole plate and by image processing the pattern. The implementation of the present utility model is not limited to the specific embodiments described above. As long as the method, device and system for reproducing the amplitude and phase distribution of the object wave by recording the Fraunhofer diffraction intensity pattern after the object wave passes through the multi-pinhole plate described in the above method, all belong to this utility model. A new type of protection.
Claims (7)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009200194783U CN201421432Y (en) | 2009-02-20 | 2009-02-20 | A coherent diffraction imaging processing device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009200194783U CN201421432Y (en) | 2009-02-20 | 2009-02-20 | A coherent diffraction imaging processing device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN201421432Y true CN201421432Y (en) | 2010-03-10 |
Family
ID=41808092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009200194783U Expired - Fee Related CN201421432Y (en) | 2009-02-20 | 2009-02-20 | A coherent diffraction imaging processing device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN201421432Y (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103502203A (en) * | 2010-10-06 | 2014-01-08 | 欧塞拉治疗有限公司 | Methods of making L-ornithine phenyl acetate |
| CN107320120A (en) * | 2017-08-01 | 2017-11-07 | 中国工程物理研究院激光聚变研究中心 | A kind of X-ray frame imaging device |
-
2009
- 2009-02-20 CN CN2009200194783U patent/CN201421432Y/en not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103502203A (en) * | 2010-10-06 | 2014-01-08 | 欧塞拉治疗有限公司 | Methods of making L-ornithine phenyl acetate |
| CN107320120A (en) * | 2017-08-01 | 2017-11-07 | 中国工程物理研究院激光聚变研究中心 | A kind of X-ray frame imaging device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101482503B (en) | A coherent diffraction imaging method and processing device thereof | |
| Shi et al. | Parametric study on light field volumetric particle image velocimetry | |
| JP7161777B2 (en) | Holographic imaging device and data processing method used in same | |
| Hinsch | Three-dimensional particle velocimetry | |
| Hinsch | Holographic particle image velocimetry | |
| CN102749834B (en) | Optical microscopic imaging system and imaging method | |
| Soria et al. | Towards 3C-3D digital holographic fluid velocity vector field measurement—tomographic digital holographic PIV (Tomo-HPIV) | |
| CN103838124A (en) | Imaging-view-field-increase-oriented lamination scanning digital holography | |
| CN109297874A (en) | A holographic real-time measurement method and device for measuring particle size of moving particles | |
| CN103322940B (en) | A kind of method obtaining microscopic image in three-dimensional shape | |
| CN106054570A (en) | Method for realizing large-phase reconstruction of single digital hologram by adopting intensity transmission equation | |
| CN107101974B (en) | A Novel Three-Step Lensless Coherent Diffraction Imaging Method | |
| CN105204310A (en) | Complex amplitude object reconstruction device based on optical scanning holographic technique | |
| CN110568080A (en) | Digital holographic detection system and method for wafer laser ultrasonic field | |
| CN103322941A (en) | Method for accurately acquiring three-dimensional microscopic image | |
| CN110942423B (en) | Far-field super-resolution reconstruction method based on Fourier laminated imaging | |
| CN109151419A (en) | A kind of single pixel imaging method based on Radon transform | |
| Huang et al. | Multiplane digital holography based on extrapolation iterations | |
| CN201421432Y (en) | A coherent diffraction imaging processing device | |
| CN103323938A (en) | Method for acquiring stereo microscopic image | |
| CN109597291B (en) | An Optical Scanning Holographic Image Recognition Method Based on Convolutional Neural Network | |
| CN102494614A (en) | High-precision digital speckle related measurement method | |
| CN112114501A (en) | Method for calculating Langqi shearing interference image in photoetching projection objective | |
| Lozano et al. | Use of holography in particle image velocimetry measurements of a swirling flow | |
| CN108180833A (en) | Reflective synchronous phase-shifted digital holographic apparatus and method based on light splitting pupil |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100310 Termination date: 20110220 |