CN106842540B - 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 - Google Patents
基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 Download PDFInfo
- Publication number
- CN106842540B CN106842540B CN201710182694.9A CN201710182694A CN106842540B CN 106842540 B CN106842540 B CN 106842540B CN 201710182694 A CN201710182694 A CN 201710182694A CN 106842540 B CN106842540 B CN 106842540B
- Authority
- CN
- China
- Prior art keywords
- light intensity
- phase
- field microscope
- annular
- annular diaphragm
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 42
- 238000005286 illumination Methods 0.000 title claims abstract description 40
- 230000003287 optical effect Effects 0.000 title claims abstract description 39
- 230000005540 biological transmission Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000001427 coherent effect Effects 0.000 claims abstract description 14
- 238000012546 transfer Methods 0.000 claims description 34
- 238000010587 phase diagram Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 2
- 239000007769 metal material Substances 0.000 claims description 2
- 210000001747 pupil Anatomy 0.000 claims description 2
- 238000007747 plating Methods 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 230000009466 transformation Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 7
- 238000011084 recovery Methods 0.000 description 7
- 230000004044 response Effects 0.000 description 5
- 210000004271 bone marrow stromal cell Anatomy 0.000 description 4
- 210000002919 epithelial cell Anatomy 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- TZSMWSKOPZEMAJ-UHFFFAOYSA-N bis[(2-methoxyphenyl)methyl] carbonate Chemical compound COC1=CC=CC=C1COC(=O)OCC1=CC=CC=C1OC TZSMWSKOPZEMAJ-UHFFFAOYSA-N 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000276498 Pollachius virens Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009713 electroplating Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001879 gelation Methods 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 238000005191 phase separation Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/06—Means for illuminating specimens
- G02B21/08—Condensers
- G02B21/14—Condensers affording illumination for phase-contrast observation
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/06—Means for illuminating specimens
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
- G02B21/367—Control or image processing arrangements for digital or video microscopes providing an output produced by processing a plurality of individual source images, e.g. image tiling, montage, composite images, depth sectioning, image comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B5/00—Optical elements other than lenses
- G02B5/005—Diaphragms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10141—Special mode during image acquisition
- G06T2207/10148—Varying focus
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Microscoopes, Condenser (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,首先设计一个环形光阑,然后基于部分相干成像方法,引入弱物体近似,利用环形光阑以及明场显微镜的参数计算出弱物体光学传递函数WOTF,最后利用相机采集三幅强度图像并通过反卷积求解光强传输方程得到定量相位图。本发明可有效解决云雾状低频噪声与高频模糊难以兼顾的矛盾,大大提高了重构相位的空间分辨率,使其达到明场显微镜两倍物镜数值孔径的分辨率,且对低频噪声具有较好的鲁棒性。而且无需对传统明场显微镜进行复杂的改造,可赋予明场显微镜高分辨率定量相位成像的能力。
Description
技术领域
本发明属于光学测量、成像技术,特别是一种基于光强传输方程的环形光照明高分辨率定量相位显微成像方法。
背景技术
相位恢复是光学测量与成像的一个重要技术,无论在生物医学还是工业检测领域,相位成像技术都在发挥着重要的作用。纵观光学测量技术近半个世纪的进展,最经典的相位测量方法为非干涉测量方法。然而,干涉测量法的缺点也十分明显:(1)干涉测量一般需要高度相干性的光源(如激光),从而需要较为复杂的干涉装置;(2)额外的参考光路的引入导致对于测量环境的要求变得十分苛刻;(3)高相干性的光源引入的相干散斑噪声限制了成像系统的空间分辨率与测量精度。
不同于干涉测量方法,另一类非常重要的相位测量技术并不需要借助干涉,它们统称为相位恢复。由于直接测量光波场的相位分布非常困难,而测量光波场的振幅/强度十分容易。因此,可以将利用强度分布来恢复(估算)相位这一过程视为一个数学上的“逆问题”,即相位恢复问题。相位恢复的方法还可细分为迭代法与直接法。基于光强传输方程的相位恢复是一种较为典型的直接相位恢复方法。光强传输方程是一个二阶椭圆偏微分方程,阐明了沿着光轴方向上光强度的变化量与垂直于光轴的平面上光波相位的定量关系。在已知光强轴向微分和聚焦面光强分布的情况下,通过数值求解光强传输方程可直接获取相位信息。相比与干涉法与迭代相位恢复法,其主要优点包括:(1)非干涉,仅仅通过测量物面光强直接求解相位信息,不需要引入额外参考光;(2)非迭代,通过直接求解微分方程获得相位;(3)可以很好的应用于白光照明,如传统明场显微镜中的科勒照明(illumination);(4)无需相位解包裹,直接获取相位的绝对分布,不存在一般干涉测量方法中的2π相位包裹问题;(5)无须复杂的光学系统,对于实验环境没有苛刻的要求,振动不敏感。
求解光强传输方程最主要的问题就是得到在聚焦面处光强的轴向微分,该轴向微分可通过光强的数值有限差分获得。而如何选择适当的离焦距离来求取光强的轴向微分,这其中存在一个噪声与分辨率的权衡问题。研究表明,当离焦距离过小或过大时,最后重构的相位图像中会出现较多的云雾状低频噪声或高频信息丢失非常严重。为了解决这个问题,许多研究者提出多平面(>2)的强度测量方法对轴向微分进行估计。为了更方便地分析基于光强传输方程的相位重构方法,研究者引入了相衬传递函数(CTF)(文献[1]E.D.Barone-Nugent,A.Barty,K.A.Nugent,“Quantitative phase-amplitude microscopyI:optical microscopy,”Journal of Microscopy 206,194-203(2002).文献[2]T.E.Gureyev,A.Pogany,D.M.Paganin,S.W.Wilkins,“Linear algorithms for phaseretrieval in the Fresnel region,”Optics Communications 231,53-70(2004).)。以上研究表明,在小离焦情况下,随着空间频率的减小,相衬传递函数迅速地减小到0,这也就解释了基于光强传输方程的相位重构方法对低频噪声的高敏感性。在大离焦的情况下,较高空间频率区域的光强传输方程传递函数响应与相干情况下的CTF差值变大,这就导致了最终的重构相位高频信息丢失严重,细小特征模糊。
下面对光强传输方程进行较为详细的背景介绍。
考虑一个沿着z轴传播的单色相干近轴光波场,其复振幅U(x)为其中j是虚数单位,φ(x)为所要恢复的相位分布。光强传输方程可以表示为
其中表示横向梯度算子,·表示点乘算子。表示光强的轴向微分,该轴向微分可通过有限差分近似获得。一般情况下通过引入辅助函数方程可以被转化成以下的泊松方程:
通过求解第一个泊松方程,可以解得辅助函数ψ(x),再通过求解第二个泊松方程即可求得待测物体的定量相位。对于观察样品为未染色的细胞和生物组织来说,其光强沿横向变化很小,光强传输方程可以被简化为一个泊松方程:
至此,以上两步泊松方程的解法可总结为傅里叶空间中的逆拉普拉斯算子(1/πλ|u|2)。该算子可以看作是光强传输方程的响应或传递函数,而逆拉普拉斯算子不管是低频还是高频响应都非常差,这就影响了重构相位的信噪比以及可达到的最大成像分辨率。因此,这其中就存在低频噪声与高频模糊两个主要问题。为了解决这个问题,许多研究者提出多平面(>=2)的强度测量来估计轴向微分(文献[3]M.Beleggia,M.A.Schofield,V.V.Volkov,Y.Zhu,“On the transport of intensity technique for phaseretrieval,”Ultramicroscopy 102,37-49(2004).文献[4]C.Zuo,Q.Chen,Y.Yu,A.Asundi,“Transport-of-intensity phase imaging using Savitzky-Golay differentiationfilter-theory and applications,”Optic Express 21,5346-5362(2013).文献[5]C.Zuo,Q.Chen,L.Huang,A.Asundi,“Phase discrepancy analysis and compensationfor fast Fourier transform based solution of the transport of intensityequation,”Optics Express 22(14),17172-17186(2014).)。这些方法都是通过对相位传递函数中的某些空间频率进行选择、组合或最小二乘拟合,以达到传递函数的最优化。然而随着测量平面的增加,这就延长了数据获取和处理的时间,限制了系统的成像速度,即成像系统的通量。
此外,光强传输方程的局限是仅仅给出了拍摄强度图像的相位,而不是真实物体的相位,这是因为在方程中并没有包含成像系统的参数。通过使相干CTF推广到部分相干下的弱物体光学传递函数(WOTF),光强传输方程的应用范围扩展到了部分相干领域。通过引入包含成像系统参数的弱物体光学传递函数,光强传输方程实现了在部分相干情况下物体的相位恢复。虽然部分相干照明有益于提高定量相位成像的分辨率,但是其相位重构的好坏主要还是依赖于WOTF的形式。在传统的部分相干成像中,可以通过提升圆形照明光的数值孔径来增大WOTF的截止频率,以达到提高成像分辨率的目的。如图2(a)至图2(d),传统圆形光照明时随着相干系数s(s为照明数值孔径与物镜数值孔径的比值)的变大,WOTF的截止频率也随之变大。但当照明数值孔径接近于物镜数值孔径时,WOTF在各空间频率的响应变差,导致强度图像对比度逐渐减,也就无法实现准确的相位重构。因此,传统方法利用圆形光源照明中存在低频云雾状噪声与高频模糊难以同时兼顾的矛盾,其相位重构的准确性差。
发明内容
本发明的目的在于提供一种基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,可有效解决云雾状低频噪声与高频模糊难以兼顾的矛盾,大大提高了重构相位的空间分辨率,且对噪声具有较好的鲁棒性。
实现本发明目的的技术解决方案为:一种基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,首先设计一个环形光阑,然后基于部分相干成像方法,引入弱物体近似,利用环形光阑以及明场显微镜的参数计算出弱物体光学传递函数WOTF,最后利用相机采集三幅强度图像并通过反卷积求解光强传输方程得到定量相位图。
本发明与现有技术相比,其显著优点:(1)能显著提高成像空间截止频率到相干成像系统衍射极限的2倍,从而提高相位重构的空间分辨率,而且额外的高频部分有利于相位细节的恢复。(2)能使相位对比度提高2倍以上,使低频和高频的成分显著提高,由此解决低频噪声与高频模糊的矛盾,且对低频噪声具有较好的鲁棒性。(3)无需对传统明场显微镜进行复杂的改造,即可赋予其高分辨率定量相位成像的能力。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是本发明基于光强传输方程的环形光照明高分辨率定量相位显微成像方法步骤流程示意图。
图2是传统圆形光照明时不同相干系数s和不同离焦距离下WOTF(虚部)的响应曲线。(a)s=0.1;(b)s=0.4;(c)s=0.75;(d)s=0.99。
图3(a)是所拍摄的微透镜阵列的全息图。
图3(b)是数字全息重构得到的相位图。
图3(c)是基于本发明方法恢复得到的相位图。
图3(d)是相位横截面对比,分别对应图3(b)和图3(c)中的标线区域。
图4是对颊上皮细胞样品进行了定量相位成像的实验结果示意图。图4(a)是本发明重构出的定量相位图;图4(b)是颊上皮细胞样品厚度的伪彩色三维显示图。图4(c)~(e)是图4(a)中所选三个区域的局部放大图。图4(f)是图4(e)中箭头所标出的截面相位曲线图。
图5是对BMSC人体骨髓基质细胞样品进行定量相位成像的实验结果示意图。图5(a1)是相干系数为1的环形光照明下离焦距离为0.5μm沿光轴方向采集3幅光强图像后恢复的相位图;图5(b1)~(e1)分别是相干系数为0.22、0.44、0.78、0.99的圆形光照明下离焦距离为0.5μm沿光轴方向采集3幅光强图像后恢复的相位图;图5(a2)~(e2)分别是图5(a1)中对应区域的局部放大图;图5(a3)~(e3)分别是图5(a2)~(e2)中箭头所标出的截面的相位曲线图。
图6是环形光阑的结构示意图。
具体实施方式
结合图1,本发明基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,具体实施步骤如下:
步骤一,设计一个环形光阑。该环形光阑可由3D打印机打印,也可以通过玻璃或金属材料黑色阳极电镀工艺制成(与传统相差板加工工艺相同),还可以由可编程LED或LCD编码实现。所设计的环形光阑主要有两个参数:一个是环形光阑的外径r,另一个是圆环宽度△s。该环形光阑的环形区域透光,中心为黑色不透光区域,该中心黑色不透光区域具有较高的遮光率并且整个环形光阑具有良好的热稳定性,环形光阑的外径r等于明场显微镜的物镜光瞳半径,如图6所示。考虑到圆环宽度△s对透光强弱的影响,一般取外径r的5%~15%,本发明举例时选取10%(即△s=10%r),将该环形光阑放置在明场显微镜的聚光镜孔径光阑处,并将伯特兰透镜放入明场显微镜的目镜观察筒的筒镜中,通过目镜观察物镜后焦面以确保环形光阑内切于物镜光曈。
步骤二,基于部分相干成像的(方法)原理,引入弱物体近似,利用环形光阑以及明场显微镜的参数计算出弱物体光学传递函数(WOTF),即弱物体光学传递函数的公式为:
其中,u和u'是傅里叶空间的二维坐标,k是波数,λ是明场显微镜的照明光波长,S(u)表示环形光阑平面光源的光强分布,|P(u)|表示明场显微镜的物镜光曈函数,Δz表示明场显微镜沿光轴的离焦距离,WOTF为求出的弱物体光学传递函数。其中该传递函数的虚部为弱物体的相位传递函数。
所设计的环形光阑的数学表达式为:
其中,为频率坐标,该频率坐标被归一化到相干衍射极限(NAobj/λ,NAobj表示物镜数值孔径,λ表示系统照明光波长),s1和s2分别表示环形光阑内圈和外圈的数值孔径(根据所设计的光阑s2=1,s1=90%),表示环形光阑平面光源的光强分布,表示明场显微镜的物镜光曈函数。根据明场显微镜的照明波长λ以及离焦距离Δz(这里可以对Δz选取不同的值,但该量不宜过大,一般取2μm以内)代入弱物体光学传递函数公式(1)中,计算出在特定离焦距离和照明波长下该环形光阑的弱物体光学传递函数。
步骤三,利用相机沿光轴拍摄三幅强度图像,并通过反卷积求解光强传输方程得到定量相位图,具体步骤如下:
第一步,通过手动调焦或电控离焦控制相机同步采集欠焦、聚焦、过焦三幅光强图像,分别记作I△z(x),I0(x)与I-Δz(x),它们的离焦距离分别是Δz,0,-Δz。
第二步,对三幅光强图像进行傅里叶变换,分别记作与将以及步骤二中计算得出的弱物体光学传递函数(WOTF)代入方程(4)中,通过反卷积求出相位的傅里叶变换
第三步,对进行逆傅里叶变换即可求得定量相位φ(x),从而实现了高分辨率定量相位的显微成像。
综合以上步骤可知,本发明通过设计一个环形光阑,并求得由环形光阑以及明场显微镜参数决定的弱物体光学传递函数,从而使得传递函数各频率的响应显著提高。由于最终重构相位中的低频和高频成分增多,其可有效的解决传统圆形照明时低频云雾状噪声与高频模糊难以兼顾的矛盾。与此同时,由于使用高数值孔径的环形照明,这就增大明场显微镜的截止频率,其可显著提高明场显微镜的成像分辨率。此外,本方法无需对传统明场显微镜进行复杂的改造,可赋予明场显微镜高分辨率定量相位成像的能力。
下面给出本发明在微光学元件上的实验结果。样品是一间距为300μm的平凸石英微透镜阵列(MLA300-7AR,Thorlabs)。实验中,光强轴向微分信号是通过两幅±1μm距离的过焦与欠焦光强分布采用中心有限差分得到的。图3(a)给出了所拍摄的微透镜阵列的全息图。图3(c)给出了本发明重构出的定量相位。注意尽管相位的范围已经超过了2π,光强传输方程仍可以直接恢复出绝对相位且不需要相位解包裹。为了验证本发明方法的准确性,将结果与数字全息显微镜成像结果进行对比(λ=632.8nm,60x,NA=0.85),如图3(b)所示。本发明还从两者结果中选取了一条截面曲线进行定量比较,如图3(d)所示。两条曲线符合的很好,误差有可能是由于噪声的影响以及数字全息重构中定义最佳聚焦面的不准确性。
本发明对颊上皮细胞样品进行了定量相位成像,使用的油浸显微物镜(UPLSAPO100X,Olympus)放大率为100倍,数值孔径为1.4,照明光的中心波长为550nm。如图4所示,图4(a)是本发明重构出的定量相位图。图4(b)是颊上皮细胞样品厚度的伪彩色三维显示图。图4(c)~(e)是图4(a)中所选三个区域的局部放大图。图4(f)是图4(e)中箭头所标出的截面相位曲线图。从图4(f)可以看出,最小可分辨的球形结构间距为208nm(相当于2.66的有效数值孔径),这意味着在定量相位成像中横向分辨率至少可达到208nm。从实验结果可以看出,由本发明恢复的相位在重建分辨率上得到了很大的提升(达到2倍物镜数值孔径的分辨率)。
此外,本发明还选取了BMSC人体骨髓基质细胞样品进行了定量相位成像,使用的显微物镜(UPLSAPO 40X,Olympus)放大率为40倍,数值孔径为0.9,照明光的中心波长为550nm。如图5所示,图5(a1)是相干系数为1的环形光照明下离焦距离为0.5μm沿光轴方向采集3幅光强图像后恢复的相位图。图5(b1)~(e1)分别是相干系数为0.22、0.44、0.78、0.99的圆形照明下离焦距离为0.5μm沿光轴方向采集3幅光强图像后恢复的相位图。图5(a2)~(e2)分别是图5(a1)中对应区域的局部放大图。图5(a3)~(e3)分别是图5(a2)~(e2)中箭头所指截面的相位曲线图。对比这五个不同照明模式下恢复的相位图结果可以看出使用环形光照明恢复出的相位比任一种圆形光照明的结果都好。本发明重构的相位图背景更均匀,在去除低频噪声的基础下,高频细节得到最大程度的保留。由此证明本发明能够有效的解决低频云雾噪声与高频模糊难以兼顾的矛盾,对低频噪声具有较好的鲁棒性。
Claims (4)
1.一种基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,其特征在于首先设计一个环形光阑,然后基于部分相干成像方法,引入弱物体近似,利用环形光阑以及明场显微镜的参数计算出弱物体光学传递函数WOTF,最后利用相机采集三幅强度图像并通过反卷积求解光强传输方程得到定量相位图;
环形光阑有两个参数:一个是环形光阑的外径r,另一个是圆环宽度Δs,Δs=(5%~15%)r;该环形光阑的环形区域透光,中心为黑色不透光区域,外径r等于明场显微镜的物镜光瞳半径;将该环形光阑放置在明场显微镜的聚光镜孔径光阑处,并将伯特兰透镜放入明场显微镜的目镜观察筒的筒镜中,通过目镜观察物镜后焦面以确保该环形光阑内切于物镜光曈。
2.根据权利要求1所述的基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,其特征在于设计的环形光阑由3D打印机打印,或者通过玻璃/金属材料黑色阳极电镀工艺制成,或者由可编程LED或LCD编码实现。
3.根据权利要求1所述的基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,其特征在于计算弱物体光学传递函数WOTF,即弱物体光学传递函数的公式为:
其中,u和u'是傅里叶空间的二维坐标,k是波数,λ是明场显微镜的照明光波长,S(u)表示环形光阑平面光源的光强分布,|P(u)|表示明场显微镜的物镜光曈函数,Δz表示明场显微镜沿光轴的离焦距离,WOTF为求出的弱物体光学传递函数,该传递函数的虚部为弱物体的相位传递函数;
所设计的环形光阑为:
其中,为频率坐标,该频率坐标被归一化到相干衍射极限,即NAobj/λ,NAobj表示物镜数值孔径,λ表示系统照明光波长;s1和s2分别表示环形光阑内圈和外圈的数值孔径,表示环形光阑平面光源的光强分布,表示明场显微镜的物镜光曈函数;
根据明场显微镜的照明波长λ以及离焦距离Δz代入弱物体光学传递函数公式(1)中,计算出在特定离焦距离和照明波长下该环形光阑的弱物体光学传递函数。
4.根据权利要求1所述的基于光强传输方程的环形光照明高分辨率定量相位显微成像方法,其特征在于利用相机沿光轴拍摄三幅强度图像,并通过反卷积求解光强传输方程得到定量相位图,步骤如下:
第一步,通过手动调焦或电控离焦控制相机同步采集欠焦、聚焦、过焦三幅光强图像,分别记作IΔz(x),I0(x)与I-Δz(x),它们的离焦距离分别是Δz,0,-Δz;
第二步,对三幅光强图像进行傅里叶变换,分别记作与将以及计算得出的弱物体光学传递函数WOTF代入方程(4)中,通过反卷积求出相位的傅里叶变换
第三步,对进行逆傅里叶变换即可求得定量相位φ(x),从而实现了高分辨率定量相位的显微成像。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710182694.9A CN106842540B (zh) | 2017-03-24 | 2017-03-24 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
US16/496,548 US11106029B2 (en) | 2017-03-24 | 2018-02-26 | Annular-irradiation high-resolution quantitative phase microimaging method using an annular illumination aperture |
PCT/CN2018/077214 WO2018171383A1 (zh) | 2017-03-24 | 2018-02-26 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710182694.9A CN106842540B (zh) | 2017-03-24 | 2017-03-24 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842540A CN106842540A (zh) | 2017-06-13 |
CN106842540B true CN106842540B (zh) | 2018-12-25 |
Family
ID=59130786
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710182694.9A Active CN106842540B (zh) | 2017-03-24 | 2017-03-24 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11106029B2 (zh) |
CN (1) | CN106842540B (zh) |
WO (1) | WO2018171383A1 (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106842540B (zh) | 2017-03-24 | 2018-12-25 | 南京理工大学 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
CN107228846B (zh) * | 2017-05-10 | 2019-09-17 | 暨南大学 | 基于离轴光束焦面共轭的荧光成像光切片方法和装置 |
CN107272178B (zh) * | 2017-08-04 | 2019-06-21 | 南京理工大学 | 基于半环状led照明的差分相衬显微成像方法 |
CN107290846B (zh) | 2017-08-04 | 2019-06-21 | 南京理工大学 | 基于环状可编程led照明的定量相位显微成像方法 |
CN107966801A (zh) * | 2017-12-29 | 2018-04-27 | 南京理工大学 | 一种基于环形照明的高速傅立叶叠层成像装置及重构方法 |
CN109580457B (zh) | 2018-11-01 | 2021-05-04 | 南京理工大学 | 基于led阵列编码照明的三维衍射层析显微成像方法 |
CN109581645B (zh) * | 2018-11-22 | 2020-07-24 | 南京理工大学 | 基于光强传输方程的相衬与微分干涉相衬的显微成像方法 |
CN109685745B (zh) * | 2019-01-02 | 2022-07-15 | 西北工业大学 | 一种基于深度学习的相位显微成像方法 |
CN110132847A (zh) * | 2019-05-29 | 2019-08-16 | 东北大学 | 一种便携式共振光声池 |
CN110441983B (zh) * | 2019-07-24 | 2022-09-09 | 成都仲伯科技有限公司 | 基于光学传递函数的x光高分辨率成像方法 |
CN111610623B (zh) * | 2020-04-29 | 2022-05-27 | 南京理工大学 | 基于傅里叶叠层的显微成像方法 |
CN111929881A (zh) * | 2020-08-06 | 2020-11-13 | 东南大学江北新区创新研究院 | 一种基于色散的相位物体成像设备及方法 |
CN113759535B (zh) * | 2021-07-26 | 2024-05-31 | 南京理工大学智能计算成像研究院有限公司 | 一种基于多角度照明反卷积的高分辨率显微成像方法 |
CN114296228B (zh) * | 2021-12-31 | 2024-01-23 | 南京理工大学智能计算成像研究院有限公司 | 一种基于混合照明模式的高分辨率显微成像方法 |
CN114509869B (zh) * | 2021-12-31 | 2024-04-26 | 南京理工大学智能计算成像研究院有限公司 | 一种基于多半径环形照明的显微成像装置及反卷积方法 |
WO2024113316A1 (zh) * | 2022-12-01 | 2024-06-06 | 骆远 | 可携式量化差分相位对比显微模块以及其深度学习重建方法 |
CN116612051B (zh) * | 2023-05-05 | 2024-03-19 | 之江实验室 | 一种基于频谱域滤波的干涉散射图像弱信号的实时恢复方法和系统 |
CN117111283B (zh) * | 2023-08-25 | 2024-07-09 | 南京理工大学 | 基于电控变焦透镜的高速多模态景深延拓显微成像系统与方法 |
CN118587094A (zh) * | 2024-07-31 | 2024-09-03 | 南京理工大学 | 一种基于信息欠采样的无透镜定量相位成像及超分辨方法 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5774222A (en) * | 1994-10-07 | 1998-06-30 | Hitachi, Ltd. | Manufacturing method of semiconductor substrative and method and apparatus for inspecting defects of patterns on an object to be inspected |
US20060052321A1 (en) * | 2002-04-05 | 2006-03-09 | Raitano Arthur B | Nucleic acid and corresponding protein entitled 98P4B6 useful in treatment and detection of cancer |
JP2005294801A (ja) * | 2004-03-11 | 2005-10-20 | Advanced Lcd Technologies Development Center Co Ltd | レーザー結晶化装置及びレーザー結晶化方法 |
US7787588B1 (en) * | 2008-07-21 | 2010-08-31 | Xradia, Inc. | System and method for quantitative reconstruction of Zernike phase-contrast images |
EP2553513B1 (en) * | 2010-03-28 | 2020-10-28 | Ecole Polytechnique Federale de Lausanne (EPFL) | Complex index refraction tomography with improved resolution |
JP2015052663A (ja) * | 2013-09-06 | 2015-03-19 | キヤノン株式会社 | 画像処理方法、画像処理装置、撮像装置およびプログラム |
US10033996B2 (en) * | 2014-02-06 | 2018-07-24 | Bar Ilan University | System and method for imaging with pinhole arrays |
EP3123449B1 (en) * | 2014-03-25 | 2020-01-01 | Apple Inc. | Method and system for representing a virtual object in a view of a real environment |
EP3146501B1 (en) * | 2014-05-19 | 2020-09-09 | The Regents of the University of California | Fourier ptychographic microscopy with multiplexed illumination |
US10473916B2 (en) * | 2014-09-30 | 2019-11-12 | Washington University | Multiple-view compressed-sensing ultrafast photography (MV-CUP) |
US9967541B2 (en) * | 2014-11-05 | 2018-05-08 | The Regents Of The University Of Colorado, A Body Corporate | 3D imaging, ranging, and/or tracking using active illumination and point spread function engineering |
US10393655B2 (en) * | 2014-11-24 | 2019-08-27 | Everfine Photo-E-Information Co., Ltd. | Apparatus for retro-reflection measurement using a holed mirror and circular aperture for annular band detection |
AU2016211635A1 (en) * | 2015-01-26 | 2017-06-29 | California Institute Of Technology | Multi-well fourier ptychographic and fluorescence imaging |
CN105222901A (zh) * | 2015-09-24 | 2016-01-06 | 程灏波 | 基于光强传输方程的实时定量相位检测装置 |
CN105158894A (zh) * | 2015-09-29 | 2015-12-16 | 南京理工大学 | 基于彩色led阵列照明的无透镜相位显微层析装置及图像重构方法 |
CN105403508B (zh) * | 2015-09-29 | 2019-05-07 | 南京理工大学 | 基于合成相位传递函数的非干涉相位成像方法 |
CN105675151B (zh) * | 2016-01-06 | 2018-10-23 | 上海大学 | 一种基于光强传输方程相位恢复应用装置 |
US11628610B2 (en) * | 2017-02-08 | 2023-04-18 | Cryovac, Llc | Process for in-line inspection of functional film layer containing detectable component |
CN106842540B (zh) * | 2017-03-24 | 2018-12-25 | 南京理工大学 | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 |
CN107290846B (zh) | 2017-08-04 | 2019-06-21 | 南京理工大学 | 基于环状可编程led照明的定量相位显微成像方法 |
-
2017
- 2017-03-24 CN CN201710182694.9A patent/CN106842540B/zh active Active
-
2018
- 2018-02-26 US US16/496,548 patent/US11106029B2/en active Active
- 2018-02-26 WO PCT/CN2018/077214 patent/WO2018171383A1/zh active Application Filing
Also Published As
Publication number | Publication date |
---|---|
CN106842540A (zh) | 2017-06-13 |
US11106029B2 (en) | 2021-08-31 |
WO2018171383A1 (zh) | 2018-09-27 |
US20210103135A1 (en) | 2021-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842540B (zh) | 基于光强传输方程的环形光照明高分辨率定量相位显微成像方法 | |
Fan et al. | Smart computational light microscopes (SCLMs) of smart computational imaging laboratory (SCILab) | |
CN109580457B (zh) | 基于led阵列编码照明的三维衍射层析显微成像方法 | |
Zuo et al. | High-speed transport-of-intensity phase microscopy with an electrically tunable lens | |
US20150100278A1 (en) | Systems and methods for quantitative phase imaging with partially coherent illumination | |
WO2021097916A1 (zh) | 高保真图像重构方法、系统、计算机设备和存储介质 | |
Kostencka et al. | Autofocusing method for tilted image plane detection in digital holographic microscopy | |
CN111121675B (zh) | 一种用于微球表面显微干涉测量的视场扩展方法 | |
CN110058392A (zh) | 一种基于光强传输方程的散斑定量相位成像系统及其方法 | |
CN112327473A (zh) | 无透镜显微成像系统及基于平均投影迭代的图像重构方法 | |
CN110619680A (zh) | 一种基于图变分的三维断层相位显微镜重建方法 | |
Malik et al. | A practical criterion for focusing of unstained cell samples using a digital holographic microscope | |
Strassberg et al. | Single-shot quantitative phase imaging with polarization differential interference contrast | |
Xue et al. | Quantitative interferometric microscopy cytometer based on regularized optical flow algorithm | |
Kemper et al. | Techniques and applications of digital holographic microscopy for life cell imaging | |
CN108072614A (zh) | 一种基于非均匀傅里叶变换的干涉合成孔径显微方法 | |
Fan et al. | Accurate dynamic quantitative phase imaging using multi-wavelength multiplexing | |
Garud et al. | Volume visualization approach for depth-of-field extension in digital pathology | |
WO2003034010A1 (en) | Phase determination of a radiation wavefield | |
WO2021155378A1 (en) | Systems and methods for performing multiple-wavelength quantitative phase imaging (qpi) | |
Kume et al. | Theoretical analysis on coherent optical super-resolution method for inspection of functional micro-structured surfaces with complex-amplitude-response distribution | |
Sha et al. | Lensless coherent diffraction imaging based on spatial light modulator with unknown modulation curve | |
Kwon et al. | Numerical phase-detection autofocusing method for digital holography reconstruction processing | |
Hu et al. | Transport of intensity phase microscopy combined with accelerated iteration for quantitative phase imaging | |
Yılmaz | Quantitative determination of surface morphology of red blood cell |
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 |