CN112270741B - 一种基于偏振的光栅相位衬度成像系统及方法 - Google Patents
一种基于偏振的光栅相位衬度成像系统及方法 Download PDFInfo
- Publication number
- CN112270741B CN112270741B CN202011092226.0A CN202011092226A CN112270741B CN 112270741 B CN112270741 B CN 112270741B CN 202011092226 A CN202011092226 A CN 202011092226A CN 112270741 B CN112270741 B CN 112270741B
- Authority
- CN
- China
- Prior art keywords
- grating
- light
- polarization
- angle
- phase contrast
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 80
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000003287 optical effect Effects 0.000 claims abstract description 22
- 230000001427 coherent effect Effects 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims description 32
- 238000005286 illumination Methods 0.000 claims description 17
- 238000007781 pre-processing Methods 0.000 claims description 9
- 239000007769 metal material Substances 0.000 claims description 6
- 238000013500 data storage Methods 0.000 claims description 3
- 230000001678 irradiating effect Effects 0.000 claims 1
- 238000011084 recovery Methods 0.000 abstract description 13
- 230000010354 integration Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000001493 electron microscopy Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/42—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
- G02B27/4205—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect having a diffractive optical element [DOE] contributing to image formation, e.g. whereby modulation transfer function MTF or optical aberrations are relevant
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/42—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
- G02B27/4272—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect having plural diffractive elements positioned sequentially along the optical path
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B5/00—Optical elements other than lenses
- G02B5/18—Diffraction gratings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
- G06T15/506—Illumination models
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Optics & Photonics (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种基于偏振的光栅相位衬度成像系统及方法,涉及领域,包括:光源装置,发射部分相干光;准直装置,将部分相干光转换为准直光;线偏振片,将准直光转换为线偏振光;光栅组件,上方放置有一样品,不同偏振角度的线偏振光分别以一预设照射角度依次照射在样品上形成折射光线,经由光栅组件处理形成强度变化的光信号;图像探测器,接收光信号以形成样品在每个偏振角度对应的相位衬度图像;三维重建系统,用于根据各相位衬度图像及对应的偏振角度进行三维重建得到样品的三维图像。有益效果是具备了无透镜成像兼具大视场和高分辨率的特点,同时降低了系统的复杂度,进一步提高了设备的集成度,可以更好地推动无透镜相位恢复系统的实用化。
Description
技术领域
本发明涉及相位恢复技术领域,尤其涉及一种基于偏振的光栅相位衬度成像系统及方法。
背景技术
相位恢复是光学测量与成像技术领域的一个重要课题。在某些特定领域,如光学测量、材料物理学、自适应光学、X射线衍射光学、电子显微学、生物医学成像领域,大部分样本都属于相位物体。这类物体的振幅透射率分布均匀,但折射率或厚度的空间分布不均匀,因此相位物体的光波振幅改变甚小,相位改变却非常大。人眼或其他光学探测器都只能判断物体的振幅变化而无法判断其相位的变化,因此也就不能“看见”相位物体,即不能区分相位物体内厚度或折射率不同的各个部分。所以对于这些领域,获取相位信息显得尤为重要。
由于直接测量光波场的相位分布非常困难,而测量光波场的振幅或强度十分容易。因此,可以将由强度分布来恢复相位这一过程作为数学上的一个逆问题,即相位恢复问题。基于非干涉的相位恢复方法可以大体分为两大类,即迭代法与直接法。
迭代法:基于迭代运算的相位复原法于1972年由Gerchberg等在研究电子显微镜的相位恢复问题时首次提出,该算法称为GS(Gerchberg-Saxton)算法。该方法指出,当待测光波场在像平面和远场衍射平面的光强分布已知时,可以通过衍射计算迭代的方式求解光场波前相位。GS算法的提出具有开创性的意义,但也存在很多问题。这些问题一方面是求解逆问题时所固有的,譬如解的存在性与唯一性问题。另一方面则是由迭代算法本身引起的,如算法经过最初几次迭代后收敛速度减慢甚至陷入停滞、陷于局部(非全局)极小值。
直接法:光强传输方程于1983年由Teague利用亥姆霍兹方程在傍轴近似条件下首次推导得出。TIE是一个二阶椭圆型偏微分方程,阐明了平行于光轴方向上光强度的变化量与垂直于光轴平面上光波的相位的定量关系。该方法区别于迭代相位恢复算法的一大特点是没有利用传统的衍射计算公式来迭代恢复相位,而是在待求平面上的光强分布(直接测量)以及光强轴向微分(通过采集离焦光强进行数值差分估计得到)已知的情况下,通过数值求解光强传输方程直接获取相位信息,从而不需要任何迭代求解过程。
目前,国内外大多数的相位恢复研究都是在显微镜系统上开展的。然而显微镜系统在大视场的情况下难以获得高分辨率的固有矛盾限制了它在实际中的进一步应用。无透镜成像技术具有大视场与高分辨率兼得的特点,但目前针对无透镜的相位恢复技术操作繁琐,算法复杂度高,使其难以在实际中获得广泛的应用。
发明内容
针对现有技术中存在的问题,本发明提供一种基于偏振的光栅相位衬度成像系统,包括:
一光源装置,用于发射部分相干光;
一准直装置,设置于所述光源装置的光线出射口,用于将所述部分相干光转换为准直光;
一线偏振片,设置于所述准直装置的下方,用于将所述准直光转换为线偏振光;
一光栅组件,设置于所述线偏振片的下方,所述光栅组件的上方放置有一样品,不同偏振角度的所述线偏振光分别以一预设照射角度依次照射在所述样品上形成折射光线,所述折射光线经由所述光栅组件处理形成强度变化的光信号;
一图像探测器,设置于所述光栅组件的下方,用于接收所述光信号以形成所述样品在每个所述偏振角度对应的相位衬度图像;
一三维重建系统,连接所述图像探测器,用于根据各所述相位衬度图像及对应的所述偏振角度进行三维重建得到所述样品的三维图像。
优选的,所述三维重建系统包括:
一第一处理单元,用于分别提取每个所述相位衬度图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第一比值;
一第二处理单元,连接所述第一处理单元,用于根据各所述第一比值及对应的所述偏振角度处理得到对应的所述折射光线入射至所述光栅组件的的入射角度;
一重建单元,连接所述第二处理单元,用于根据各所述入射角度进行三维重建得到所述样品的三维图像。
优选的,所述三维重建系统还包括一预处理单元,连接所述第二处理单元,所述预处理单元包括:
一图像获取模块,用于获取所述光栅组件上方未放置所述样品时,不同偏振角度的所述线偏振光分别由预定的复数个基准照射角度直接照射在所述光栅组件上,所述图像探测器采集得到的每个所述基准照射角度对应的基准图像;
一图像处理模块,连接所述图像获取模块,用于分别提取每个所述基准图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第二比值;
一数据存储模块,分别连接所述图像获取模块和所述图像处理模块,用于建立所述偏振角度、所述基准照射角度与对应的所述第二比值之间的对应关系曲线并保存;
所述第二处理单元根据各所述相位衬度图像中对应的所述第一比值和所述偏振角度于所述对应关系曲线中匹配得到对应的所述基准照射角度作为所述入射角度。
优选的,所述光栅组件集成于所述图像探测器的上方。
优选的,所述图像探测器为包含多个像素的面阵探测器。
优选的,所述光栅组件包括一第一光栅和一第二光栅,所述第一光栅与所述第二光栅相互垂直。
优选的,所述第一光栅和所述第二光栅之间的间距的取值范围为[400纳米,600纳米]。
优选的,所述第一光栅和所述第二光栅由金属材料制成,所述金属材料的厚度的取值范围为[100纳米,300纳米]。
优选的,所述第一光栅的周期的取值范围为[400纳米,600纳米],所述第二光栅的周期的取值范围为[800纳米,1000纳米]。
优选的,所述线偏振片沿所述图像探测器的X轴或Y轴旋转不同的旋转角度形成所述基准照射角度。
优选的,所述基准照射角度的取值范围为[-90度,90度]。
一种基于偏振的光栅相位衬度成像方法,应用于上述的基于偏振的光栅相位衬度成像系统,包括:
步骤S1、发射部分相干光;
步骤S2、将所述部分相干光转换为准直光;
步骤S3、将所述准直光转换为线偏振光;
步骤S4、不同偏振角度的所述线偏振光分别以一预设照射角度照射在放置于光栅组件上的样品上形成折射光线,所述折射光线经由所述光栅组件处理形成强度变化的光信号;
步骤S5、接收所述光信号以形成所述样品的在每个所述偏振角度对应的相位衬度图像;
步骤S6、根据各所述相位衬度图像及对应的所述偏振角度进行三维重建得到所述样品的三维图像。
优选的,所述步骤S6包括:
步骤S61、分别提取每个所述相位衬度图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第一比值;
步骤S62、根据各所述第一比值和对应的所述偏振角度处理得到对应的所述折射光线入射至所述光栅组件的的入射角度;
步骤S63、根据各所述入射角度进行三维重建得到所述样品的三维图像。
优选的,执行所述步骤S62之前,还包括一预处理过程,包括:
步骤A1、获取所述光栅组件上方未放置所述样品时,不同偏振角度的所述线偏振光分别由预定的复数个基准照射角度直接照射在所述光栅组件上,所述图像探测器采集得到的每个所述基准照射角度对应的基准图像;
步骤A2、分别提取每个所述基准图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第二比值;
步骤A3、建立所述偏振角度、所述基准照射角度与对应的所述第二比值之间的对应关系曲线并保存;
所述步骤S62中,根据各所述第一比值和对应的所述偏振角度于所述对应关系曲线中匹配得到对应的所述基准照射角度作为所述入射角度。
上述技术方案具有如下优点或有益效果:
1)将光栅组件与图像探测器相结合,使成像系统具备了无透镜成像兼具大视场和高分辨率的特点,同时降低了系统的复杂度,降低了操作的复杂度与成像设备的体积,进一步提高了设备的集成度,可以更好地推动无透镜相位恢复系统的实用化;
2)由于物体对线偏振光不敏感,通过设置线偏振片,能够有效提升获取的相位衬度图像的准确度,进而提升三维重建的精准度,尤其是对于厚度相对较大的物体。
附图说明
图1为本发明的较佳的实施方式中,一种基于偏振的光栅相位衬度成像系统的结构示意图;
图2为本发明的较佳的实施方式中,一种基于偏振的光栅相位衬度成像系统的原理示意图;
图3为本发明的较佳的实施方式中,对应关系曲线的示意图;
图4为本发明的较佳的实施方式中,未放置样品时,图像探测器记录的基准图像;
图5为本发明的较佳的实施方式中,放置样品后,图像探测器记录的相位衬度图像;
图6为本发明的较佳的实施方式中,三维重建结果的示意图;
图7为本发明的较佳的实施方式中,光栅组件的结构示意图;
图8为本发明的较佳的实施方式中,一种基于偏振的光栅相位衬度成像方法的流程示意图;
图9为本发明的较佳的实施方式中,根据相位衬度图像进行三维重建的过程示意图;
图10为本发明的较佳的实施方式中,预处理过程的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本发明并不限定于该实施方式,只要符合本发明的主旨,则其他实施方式也可以属于本发明的范畴。
本发明的较佳的实施例中,基于现有技术中存在的上述问题,现提供一种基于偏振的光栅相位衬度成像系统,如图1和图2所示,包括:
一光源装置1,用于发射部分相干光;
一准直装置2,设置于光源装置1的光线出射口,用于将部分相干光转换为准直光;
一线偏振片3,设置于准直装置2的下方,用于将准直光转换为线偏振光;
一光栅组件4,设置于线偏振片3的下方,光栅组件4的上方放置有一样品5,不同偏振角度的线偏振光分别以一预设照射角度依次照射在样品5上分别形成折射光线,折射光线经由光栅组件4处理形成强度变化的光信号;
一图像探测器6,设置于光栅组件4的下方,用于接收光信号以形成样品5的在每个偏振角度对应的相位衬度图像;
一三维重建系统7,连接图像探测器6,用于根据各相位衬度图像及对应的偏振角度进行三维重建得到样品的三维图像。
具体地,在成像领域中,为了在高分辨率的前提下获得大的视场,常采用的技术手段是图像拼接法,但是这样做需要采用精密的机械移动与复杂的算法计算,不仅耗时而且精度也较低。为了进一步解决该问题,无透镜成像是一个趋势。本实施方式中,通过将光栅组件4和图像探测器6相结合,使得在相衬成像领域中应用无透镜成像技术成为可能。进一步地,由于准直光由不同的照射角度直接照射在样品上时,会出现成像点的偏移,特别是对于厚度相对较大的物体,为提升成像准确度,还需要计算成像点的偏移部分,计算过程较为复杂,而物体对于线偏振光不敏感,通过设置线偏振片3,使得准直光转换为线偏振光,进而通过线偏振光照射样品,消除成像点偏移对成像准确度的影响,操作简单,无需复杂的计算过程。
作为优选的实施方式,上述光源装置1可以是LED光源,上述准直装置2和光源装置1可以一体化设置,在调整线偏振光的照射角度时,可以通过沿图像探测器的X轴或Y轴方向旋转线偏振片实现。光源装置1产生的部分相干光经由准直器成为准直光,该准直光为平行光,该平行光进一步由线偏振片转换为线偏振光,该线偏振光照射到样品5上并穿透样品时,会发生吸收、散射、折射等物理现象,折射后的出射光线携带有特定空间位置中的材料信息,通过在样品5和图像探测器6之间设置光栅组件4,上述出射光线经过光栅组件4会产生相位变化信息,该相位变化信息的外在表现为光强的变化,即该出射光线经由光栅组件4后会形成强度变化的光信号,该强度变化的光信号中包含被折射光线的折射角信息,最后由图像探测器6接收携带有相位变化信息的光信号并记录形成相位衬度图像。由于上述相位衬度图像包含了准直光照射在样品5上产生的相位变化信息,三维重建系统7可以根据不同偏振角度的线偏振片形成的线偏振光在预设照射角度的相位衬度图像对样品4内部的相位分布进行三维重建,进而获取样品5的三维图像。该预设照射角度可以是垂直照射。整个成像系统具备了无透镜成像兼具大视场和高分辨率的特点,同时降低了系统的复杂度,只需对偏振角度进行调整,省去了繁琐的操作,可以更好地推动无透镜相位恢复系统的实用化。
作为优选的实施方式,三维重建系统7包括:
一第一处理单元71,用于分别提取每个相位衬度图像中的各像素点的光强值,并计算位于偶数行的像素点的光强值与位于奇数行的像素点的光强值之间的第一比值;
一第二处理单元72,连接第一处理单元71,用于根据各第一比值和对应的偏振角度处理得到对应的折射光线入射至光栅组件3的的入射角度;
一重建单元73,连接第二处理单元72,用于根据各入射角度进行三维重建得到样品的三维图像。
作为优选的实施方式,以垂直照射角度为例,线偏振光在以垂直于图像探测器6表面的照射角度照射样品5时,此时由于样品内部折射率不均匀或者厚度不均匀,将会使原本垂直于线偏振光的波面发生改变,而在样品5内部,线偏振光的传播方向与波面的法线重合。上述第二处理单元72优选利用已经建立好的偏振角度、入射角度与上述比值的对应关系曲线,恢复出经过样品后不同点折射光的角度,进而恢复出经过样品后的波前,由此波前则可以恢复出样品内部折射率分布以及三维结构。由波前恢复出样品内部折射率分布及三维结构的详细过程请参考此文献:“Boundary-artifact-free phase retrieval with thetransport of intensity equation II:applications to microlenscharacterization”(Optics Express,第22卷,第15期,第18310-18324页,2014年)。
作为优选的实施方式,三维重建系统7还包括一预处理单元74,连接第二处理单元72,预处理单元74包括:
一图像获取模块741,用于获取光栅组件上方未放置样品时,准直光由预定的复数个基准照射角度直接照射在光栅组件上,图像探测器采集得到的每个基准照射角度对应的基准图像;
一图像处理模块742,连接图像获取模块741,用于分别提取每个基准图像中的各像素点的光强值,并计算位于偶数行的像素点的光强值与位于奇数行的像素点的光强值之间的第二比值;
一数据存储模块743,分别连接图像获取模块741和图像处理模块742,用于建立偏振角度、基准照射角度与对应的第二比值之间的对应关系曲线并保存;
第二处理单元72根据各相位衬度图像中对应的第一比值和偏振角度于对应关系曲线中匹配得到对应的基准照射角度作为入射角度。
作为优选的实施方式,如图3所示,为上述对应关系曲线,每条对应关系曲线对应一个偏振角度,每条对应关系曲线上的各标记点分别代表每个照射角度对应的第二比值,上述预定的复数个照射角度可以是沿图像探测器6的X轴或Y轴的[-90度,90度]范围内每间隔5度进行一次照射。由于光栅组件上方未放置样品,线偏振光直接入射至光栅组件,并形成一个特定光强变化的基准图像,如图4所示,该预定的复数个基准照射角度即为入射至光栅组件的入射角度。在放置样品后,为了能提取出定量的折射角信号,最简单的途径是在光强和样品折射角之间建立线性关系。为此,在放入样品前,在不同偏振角度下,将线偏振光以不同基准照射角度照射光栅组件,从而获得不同偏振角度下的不同基准照射角度与光强变化的线性关系,在放入样品后,线偏振光照射在样品上,样品会在水平方向和垂直方向产生正的或者负的折射角,不同折射角的光线会在图像探测器上对应不同的强度,从而获得样品沿水平方向的折射衬度像和沿垂直方向的折射衬度像形成样品的相位衬度图像,如图5所示。基于无样品时所测得的不同偏振角度下的入射角度与图像探测器的偶数行像素所记录的光强与奇数行像素所记录的光强的比值的对应关系曲线和有样品时不同偏振角度的相位衬度图像,能够对样品的结构进行三维重建,如图6所示,为基于如图5所示的相位衬度图像和如图3所示的对应关系曲线进行三维重建的结果示意。
进一步地,如图3中可以看出,入射角度与对应比值之间的对应关系曲线并非单调的,在实际的测量过程中,需要改变线偏振片的偏振角度,通过在不同线偏振角度下,同一折射角度光线的比值不同,从而唯一确定该折射光线的入射角度。在实际测量中,可以参考图3,选取线偏振角度为40(P-40),60(P-60),80(P-80),100(P-100)五个线偏振角度来唯一的确定某一折射光线角度。
作为优选的实施方式,光栅组件4集成于图像探测器6的上方。
作为优选的实施方式,图像探测器6为包含多个像素的面阵探测器。
作为优选的实施方式,如图7所示,光栅组件4包括一第一光栅41和一第二光栅42,第一光栅41与第二光栅42相互垂直。
具体地,本实施方式中,面阵探测器具有多个相互垂直布置的像素,为保证图像探测器6能够准确成像,光栅组件4的结构需要与图像探测器6像素分布一致,因此需要设置第一光栅41与第二光栅42相互垂直,该相互垂直关系为空间上相互垂直,即第一光栅41和第二光栅42分布于相互平行的两个平面上,作为优选的实施方式,第一光栅41和第二光栅42之间的间距的取值范围为[400纳米,600纳米]。
作为优选的实施方式,第一光栅41和第二光栅42由金属材料制成,金属材料的厚度的取值范围为[100纳米,300纳米],该金属材料可以是铜。
作为优选的实施方式,第一光栅41的周期的取值范围为[400纳米,600纳米],第二光栅42的周期的取值范围为[800纳米,1000纳米]。
作为优选的实施方式,线偏振片沿图像探测器6的X轴或Y轴转不同的旋转角度形成基准照射角度,从而便于获得样品准确的组织结构信息。
作为优选的实施方式,基准照射角度的取值范围为[-85度,85度]。
一种基于偏振的光栅相位衬度成像方法,应用于上述的基于偏振的光栅相位衬度成像系统,如图8所示,包括:
步骤S1、发射部分相干光;
步骤S2、将部分相干光转换为准直光;
步骤S3、将准直光转换为线偏振光;
步骤S4、不同偏振角度的线偏振光分别以一预设照射角度照射在放置于光栅组件上的样品上分别形成折射光线,折射光线经由光栅组件处理形成强度变化的光信号;
步骤S5、接收光信号以形成样品的在每个偏振角度对应的相位衬度图像;
步骤S6、根据各相位衬度图像和对应的偏振角度进行三维重建得到样品的三维图像。
作为优选的实施方式,如图9所示,步骤S6包括:
步骤S61、分别提取每个相位衬度图像中的各像素点的光强值,并计算位于偶数行的像素点的光强值与位于奇数行的像素点的光强值之间的第一比值;
步骤S62、根据各第一比值和对应的偏振角度处理得到对应的折射光线的入射角度;
步骤S63、根据各入射角度进行三维重建得到样品的三维图像。
作为优选的实施方式,执行步骤S62之前,还包括一预处理过程,如图10所示,包括:
步骤A1、获取光栅组件上方未放置样品时,不同偏振角度的线偏振光分别由预定的复数个基准照射角度直接照射在光栅组件上,图像探测器采集得到的每个偏振角度对应的基准图像;
步骤A2、分别提取每个基准图像中的各像素点的光强值,并计算位于偶数行的像素点的光强值与位于奇数行的像素点的光强值之间的第二比值;
步骤A3、建立偏振角度、基准照射角度与对应的第二比值之间的对应关系曲线并保存;
步骤S62中,根据各第一比值和对应的偏振角度于对应关系曲线中匹配得到对应的照射角度作为入射角度。
以上所述仅为本发明较佳的实施例,并非因此限制本发明的实施方式及保护范围,对于本领域技术人员而言,应当能够意识到凡运用本说明书及图示内容所作出的等同替换和显而易见的变化所得到的方案,均应当包含在本发明的保护范围内。
Claims (10)
1.一种基于偏振的光栅相位衬度成像系统,其特征在于,包括:
一光源装置,用于发射部分相干光;
一准直装置,设置于所述光源装置的光线出射口,用于将所述部分相干光转换为准直光;
一线偏振片,设置于所述准直装置的下方,用于将所述准直光转换为线偏振光;
一光栅组件,设置于所述线偏振片的下方,所述光栅组件的上方放置有一样品,不同偏振角度的所述线偏振光分别以一预设照射角度依次照射在所述样品上分别形成折射光线,所述折射光线经由所述光栅组件处理形成强度变化的光信号;
一图像探测器,设置于所述光栅组件的下方,用于接收所述光信号以形成所述样品的在每个所述偏振角度对应的相位衬度图像;
一三维重建系统,连接所述图像探测器,用于根据各所述相位衬度图像及对应的所述偏振角度进行三维重建得到所述样品的三维图像;
所述三维重建系统包括:
一第一处理单元,用于分别提取每个所述相位衬度图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第一比值;
一第二处理单元,连接所述第一处理单元,用于根据各所述第一比值和对应的所述偏振角度处理得到对应的所述折射光线入射至所述光栅组件的入射角度;
一重建单元,连接所述第二处理单元,用于根据各所述入射角度进行三维重建得到所述样品的三维图像;
所述三维重建系统还包括一预处理单元,连接所述第二处理单元,所述预处理单元包括:
一图像获取模块,用于获取所述光栅组件上方未放置所述样品时,不同偏振角度的所述线偏振光分别由预定的复数个基准照射角度直接照射在所述光栅组件上,所述图像探测器采集得到的每个所述基准照射角度对应的基准图像;
一图像处理模块,连接所述图像获取模块,用于分别提取每个所述基准图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第二比值;
一数据存储模块,分别连接所述图像获取模块和所述图像处理模块,用于建立所述偏振角度、所述基准照射角度与对应的所述第二比值之间的对应关系曲线并保存;
所述第二处理单元根据各所述第一比值和对应的所述偏振角度于所述对应关系曲线中匹配得到对应的所述基准照射角度作为所述入射角度。
2.根据权利要求1所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述光栅组件集成于所述图像探测器的上方。
3.根据权利要求1所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述图像探测器为包含多个像素的面阵探测器。
4.根据权利要求1所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述光栅组件包括一第一光栅和一第二光栅,所述第一光栅与所述第二光栅相互垂直。
5.根据权利要求4所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述第一光栅和所述第二光栅之间的间距的取值范围为[400纳米,600纳米]。
6.根据权利要求4所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述第一光栅和所述第二光栅由金属材料制成,所述金属材料的厚度的取值范围为[100纳米,300纳米]。
7.根据权利要求4所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述第一光栅的周期的取值范围为[400纳米,600纳米],所述第二光栅的周期的取值范围为[800纳米,1000纳米]。
8.根据权利要求1所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述线偏振片沿所述图像探测器的X轴或Y轴旋转不同的旋转角度形成所述基准照射角度。
9.根据权利要求8所述的基于偏振的光栅相位衬度成像系统,其特征在于,所述基准照射角度的取值范围为[-85度,85度]。
10.一种基于偏振的光栅相位衬度成像方法,其特征在于,应用于如权利要求1-9中任意一项所述的基于偏振的光栅相位衬度成像系统,包括:
步骤S1、发射部分相干光;
步骤S2、将所述部分相干光转换为准直光;
步骤S3、将所述准直光转换为线偏振光;
步骤S4、不同偏振角度的所述线偏振光分别以一预设照射角度照射在放置于光栅组件上的样品上分别形成折射光线,所述折射光线经由所述光栅组件处理形成强度变化的光信号;
步骤S5、接收所述光信号以形成所述样品的在每个所述偏振角度对应的相位衬度图像;
步骤S6、根据各所述相位衬度图像及对应的所述偏振角度进行三维重建得到所述样品的三维图像;
所述步骤S6包括:
步骤S61、分别提取每个所述相位衬度图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第一比值;
步骤S62、根据各所述第一比值和对应的所述偏振角度处理得到对应的所述折射光线入射至所述光栅组件的的入射角度;
步骤S63、根据各所述入射角度进行三维重建得到所述样品的三维图像;
执行所述步骤S62之前,还包括一预处理过程,包括:
步骤A1、获取所述光栅组件上方未放置所述样品时,不同偏振角度的所述线偏振光分别由预定的复数个基准照射角度直接照射在所述光栅组件上,所述图像探测器采集得到的每个所述基准照射角度对应的基准图像;
步骤A2、分别提取每个所述基准图像中的各像素点的光强值,并计算位于偶数行的所述像素点的所述光强值与位于奇数行的所述像素点的所述光强值之间的第二比值;
步骤A3、建立所述偏振角度、所述基准照射角度与对应的所述第二比值之间的对应关系曲线并保存;
所述步骤S62中,根据各所述第一比值和对应的所述偏振角度于所述对应关系曲线中匹配得到对应的所述基准照射角度作为所述入射角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092226.0A CN112270741B (zh) | 2020-10-13 | 2020-10-13 | 一种基于偏振的光栅相位衬度成像系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092226.0A CN112270741B (zh) | 2020-10-13 | 2020-10-13 | 一种基于偏振的光栅相位衬度成像系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112270741A CN112270741A (zh) | 2021-01-26 |
CN112270741B true CN112270741B (zh) | 2024-04-09 |
Family
ID=74338020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011092226.0A Active CN112270741B (zh) | 2020-10-13 | 2020-10-13 | 一种基于偏振的光栅相位衬度成像系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112270741B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
CN108169257A (zh) * | 2017-12-08 | 2018-06-15 | 中国科学院高能物理研究所 | 具有三种定量成像机制的纳米ct三维成像方法 |
CN108680589A (zh) * | 2018-05-31 | 2018-10-19 | 北京航空航天大学 | 基于横向错位光栅的x射线光栅差分相位衬度三维锥束计算机层析成像方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102781327B (zh) * | 2009-12-10 | 2015-06-17 | 皇家飞利浦电子股份有限公司 | 相衬成像 |
US10359375B2 (en) * | 2013-10-23 | 2019-07-23 | Nanovision Technology (Beijing) Co., Ltd. | Photon count-based radiation imaging system, method and device thereof |
-
2020
- 2020-10-13 CN CN202011092226.0A patent/CN112270741B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
CN108169257A (zh) * | 2017-12-08 | 2018-06-15 | 中国科学院高能物理研究所 | 具有三种定量成像机制的纳米ct三维成像方法 |
CN108680589A (zh) * | 2018-05-31 | 2018-10-19 | 北京航空航天大学 | 基于横向错位光栅的x射线光栅差分相位衬度三维锥束计算机层析成像方法及装置 |
Non-Patent Citations (6)
Title |
---|
"A high numerical aperture, polarization-insensitive metalens for long-wavelength infrared imaging";Fan, Qingbin, et al.;《Applied Physics Letters》;全文 * |
"Anisotropy imaging using polarization and angular multiplexing";Aparna, V., et al.;《Quantitative Phase Imaging III》;全文 * |
"Digital holography and quantitative phase contrast imaging using computational shear interferometry";Falldorf, Claas, Mostafa Agour, and Ralf B. Bergmann.;《Optical Engineering 》;全文 * |
"Polarization grating based on diffraction phase microscopy for quantitative phase imaging of paramecia";Zhang, Meiling, et al.;《Optics Express》;全文 * |
"基于迭代重建算法的 X 射线光栅相位 CT 成像";戚俊成, et al.;《物理学报》;全文 * |
"X射线微结构阵列光源光栅相衬成像系统研究";昝贵彬.;《中国优秀博士学位论文全文数据库 基础科学辑》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112270741A (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kuś et al. | Holographic tomography: hardware and software solutions for 3D quantitative biomedical imaging | |
EP2736403B1 (en) | Apparatus for quantitative phase tomography through linear scanning with coherent and non-coherent detection | |
US11781966B2 (en) | 3D diffraction tomography microscopy imaging method based on LED array coded illumination | |
Kaufmann et al. | An electro-optical read-out system for analysis of speckle photographs | |
US7039553B2 (en) | Phase determination of a radiation wave field | |
US8129703B2 (en) | Intraoral imaging system and method based on conoscopic holography | |
WO2013040345A1 (en) | Systems and methods of dual-plane digital holograghic microscopy | |
US20110032586A1 (en) | Light microscope with novel digital method to achieve super-resolution | |
WO1996042066A1 (en) | Methods and apparatus for image reconstruction | |
US20130147919A1 (en) | Multi-View Difraction Grating Imaging With Two-Dimensional Displacement Measurement For Three-Dimensional Deformation Or Profile Output | |
Lu et al. | Quantitative phase imaging camera with a weak diffuser | |
Dong et al. | Automatic filtering for zero-order and twin-image elimination in off-axis digital holography | |
CN105403508B (zh) | 基于合成相位传递函数的非干涉相位成像方法 | |
Lefebvre et al. | Fully automated dual-resolution serial optical coherence tomography aimed at diffusion MRI validation in whole mouse brains | |
Chaumet et al. | Quantitative phase microscopies: accuracy comparison | |
CN112270741B (zh) | 一种基于偏振的光栅相位衬度成像系统及方法 | |
CN112258628B (zh) | 一种基于光栅相位衬度的成像系统及方法 | |
US11385164B2 (en) | Method for calibrating an analysis device, and associated device | |
Nakajima | Noniterative Phase Retrieval from a Single Diffraction Intensity Pattern<? format?> by Use of an Aperture Array | |
CN105739279A (zh) | 基于全息图强度信息的迭代重构成像方法 | |
Wang et al. | Comparison of the refocus criteria for the phase, amplitude, and mixed objects in digital holography | |
Morgan et al. | High-fidelity direct coherent diffractive imaging of condensed matter | |
Pan et al. | Digital holographic microtomography for geometric parameter measurement of optical fiber | |
US20230073901A1 (en) | Systems and methods for performing multiple-wavelength quantitative phase imaging (qpi) | |
CN109724690A (zh) | 基于光谱角度填图的纹影法远场测量焦斑自动重构方法 |
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 |