技术领域
本发明属于光学成像相关领域,涉及一种用于结构光照明显微成像中结构光初相位的提取方法。
背景技术
光学显微成像已被广泛用于生物学和生命科学等领域的研究,但成像分辨率受到光学衍射极限的限制,无法满足对细微生物结构观察的需求。为了打破衍射极限的限制,近年来发展了多种超分辨光学显微成像技术。其中,2000年Gustafsson提出的结构光照明显微术(Structured Illumination Microscopy,SIM),基于“莫尔条纹”原理,通过余弦型结构光条纹照明的方式,将高频信息编码到低频区域被探测器接收,实现2倍于普通宽场成像的分辨率。SIM仅改变了成像系统的照明部分,其它部分无需改变,能与多种成像系统兼容,成为当下使用最广泛的超分辨成像技术。
SIM并不是一种直接的超分辨成像技术,通常需要使用相移法,采集三幅不同初相位的结构光照明图像,共同重建一幅超分辨图像。SIM超分辨重建中如何精确地提取三个不同照明结构光场的初相位,是获得高质量超分辨重建成像的最关键问题。目前常用的结构光场初相位提取方法有峰值相位法、自相关法和图像重组变换法,但均存在一定的局限性。其中,峰值相位法与自相关法均是在满足前提近似条件下取得了结构光初相位的近似解。但由于近似中限定了结构光场调制度足够高,及成像中不含有未加载结构光条纹的离焦背景信息等条件,因此峰值相位法与自相关法不适用于结构光场调制度较低及厚样品成像的情况。而图像重组变换法是一种精确的初相位提取方法,不受前提近似条件的限制,其基本提取过程如下:
在结构光场间相移步长相等且为特殊值π/2的情况下,采集三幅不同初相位结构光照明样品的图像,数学表达可写为:
式中:r表示空间域的二维坐标;m、p与I
对上述的三幅图像D
式中,j为虚数单位。对上式重组图像做傅里叶变换,获得图像的频谱图可表示为:
式中,k表示频域的二维坐标;~表示傅里叶变换。
在余弦结构光的空间频率p处,仅含有唯一复指数项
从上述提取初相位的过程,可以发现图像重组变换法通过图像间相减的方式去除了离焦背景的影响,并且整个过程未作任何近似,不受结构光调制度的影响,能够提取余弦结构光场初相位的精确值。但该方法仅适用于步长相等且为特殊值π/2的情况,而实际相移过程中,往往会由于空气扰动、环境震动及相移器不可避免的误差等原因,造成余弦结构条纹间的相移值并非理想的相移值π/2,产生的相移误差会严重降低传统的图像重组变换法获取初相位的精度。
发明内容
为了解决现有技术中存在的不足,本发明提出了一种用于结构光照明显微成像中结构光初相位的提取方法。在图像重组变换法的基础上,对图像重组方式进行了改进,将图像重组变换法推广到了更一般的形式,实现在任意未知相移值下,对结构光场初相位的提取。
本发明的技术路线为:结构光照明样品成像,会引入两个额外的高频频谱分量与低频频谱分量相互混叠在一起。记录三幅不同初相位结构光照明的图像,并结合图像间的相移值,对混叠的频谱分量分离。如果准确获知图像间的实际相移值,那么就可以准确的分离出两个高频频谱分量,且此时分离的两个高频频谱分量在结构光空间频率p处的强度差值最大,但若未能准确获知图像间的实际相移值,那么在空间频率p处的强度差值偏小。基于这种特性,可通过迭代优化估计相移值的方式确定图像间的实际相移值。利用获取的图像间的实际相移值对三幅记录图像重组,在重组图像频谱图的p位置提取照明结构光场的初相位。
具体步骤如下:
步骤1记录三幅不同初相位的结构光照明图像:
通过激光照明,经空间光调制器调制在样品面生成余弦结构光场可表示为:
式中:I(r)为余弦结构光场;r表示空间域的二维坐标;m、p与I
控制空间光调制器改变结构光场的相移值,依次产生3个不同初相位的结构光场;结构光照明样品经显微系统成像,在成像系统的像面由CCD依次记录三幅不同初相位的结构光照明图像1、2、3,图像的表达式为:
式中:D
步骤2获取图像间的实际相移值:
对采集的三幅图像做傅里叶变换,获得三幅频谱图像4、5、6,频域图像表达式为:
式中:k表示频域的二维坐标;~表示傅里叶变换;j为虚数单位;
设定初始估计相移值
构造下式计算分离出的频谱分量ψ
SD=|ψ
基于频谱分量间强度差的特性,不断迭代优化估计相移值
步骤3:提取照明结构光初相位:
利用步骤2中已获取的图像间实际相移值
对重组获得的图像7做傅里叶变换(FFT)得到其频谱图8,并计算频谱图8在空间频率p处的值,可表示为:
通过MATLAB取辐角函数arg计算上式的辐角值,即可计算获得初相位
利用通过获取的实际相移值
本发明的有益效果是:
(1)本发明避免了峰值相位法与自相关法中需要限定结构光场调制度足够高,及忽视离焦背景等近似条件。本发明中未作近似处理,不会受结构光场调制度的限制,且通过图像重组的方式去除了离焦背景的影响。
(2)本发明摆脱了图像重组变换法中对特定相移值的要求。本发明中,首先通过分离频谱强度差之间的特性,获取了图像间的实际相移值,再用测得的实际相移值对图像进行重组,进而提取初相位,拓宽了图像重组变换法的使用范围,大大降低了对成像系统性能的要求。
附图说明
图1为本发明方法流程图;
图2为基于激光干涉式结构光照明荧光显微成像系统;
图3为分离频谱分量强度分布随估计相移值与准确相移值间差值的变化;
图4(a)-图4(d)为不同结构光调制度下,本发明方法与传统方法提取结构光初相位结果对比。
附图2中标记为:1—激光光源,2—第一透镜,3—第二透镜,8—第三透镜,10—第四透镜,4—第一半波片,6—第二半波片,5—反射式空间光调制器SLM,7—反射镜,9—空间滤波器,11—二向色镜,12—显微物镜,13—载物台,14—滤光片,15—筒镜,16—数字相机CCD。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步说明。
本发明的初相位提取流程如图1所示。结构光照明样品成像,记录三幅不同初相位的结构光照明图像。对记录的图像做傅里叶变换获得三幅包含未知相移值的频谱图像。为获取未知相移值的实际值,通过迭代设定初始估计相移值分离因余弦结构光照明引入频谱图像中的两个高频频谱分量,并计算分离的两个高频频谱分量在频率p处的强度差值,当强度差绝对值的负值最小时输出对应的估计相移值,即为结构光场间的实际相移值。再利用获得的实际相移值重新组合三幅记录的原始图像,得到一幅新的重组图像。对重组图像傅里叶变换进频域,在频谱图中的频率p处取辐角计算初相位的值。
基于激光干涉式结构光照明荧光显微成像系统如图2所示:包括激光光源1、设置在激光光源1后的光束扩束准直透镜组(透镜2与透镜3)、设置在光束扩束准直透镜组后的半波片4、设置在半波片4后面的反射式空间光调制器SLM5、设置在反射式空间光调制器SLM5后面的半波片6、设置在半波片6后面的反射镜7、设置在反射镜7后面的透镜8、设置在透镜8后面的空间滤波器9、设置在空间滤波器9后的透镜10、设置在透镜10后面的二向色镜11、设置在二向色镜11后面的显微物镜12与载物台13、设置在荧光光路中的滤波片14、筒镜15及数字相机CCD16。
基于上述流程,下面实施例中以荧光微球作为样品对本发明进一步解释,但本发明并不仅仅限于本实施例。
步骤1记录三幅不同初相位的结构光照明图像:
如图2中基于激光干涉式结构光照明荧光显微成像系统,以457nm的激光经半波片4调制偏振态后照射到加载余弦光栅的反射式空间光调制器SLM5上,反射产生多级衍射光束,通过空间滤波器9滤除±1级以外的衍射光束,经二向色镜反射进入显微物镜相互干渉产生二维结构光照明样品。
将荧光微球样品置于载物台上,经结构光激发照明,荧光微球发光荧光成像,所成图像被数字相机CCD记录为图像1。
通过空间光调制器对结构光完成π/2步长的相移,用CCD记录成像,记为图像2。再次通过空间光调制器对结构光完成π/2步长的相移,通过CCD记录图像,记为图像3。图像1、2、3可用数学式表示为:
式中:r表示空间域的二维坐标;m、p与I
步骤2获取图像间的实际相移值:
在计算机Matlab的编程环境下编码完成后续流程。
对记录的图像1、2、3做傅里叶变换(FFT),获得三幅频谱图像4、5、6:
式中:k表示频域的二维坐标;~表示傅里叶变换;j为虚数单位。在频谱图
设定初始估计相移值
构造下式计算分离出的频谱分量ψ
SD=|ψ
图3中给出了分离频谱分量ψ
基于图3中频谱分量的特性,迭代寻找实际的相移值
步骤3:提取余弦结构光场初相位:
在步骤2中获取了图像间实际相移值
对重组获得的图像7做傅里叶变换(FFT)得到其频谱图8,并计算频谱图8在空间频率p处的值,可表示为:
通过MATLAB取辐角函数arg计算上式的辐角值,即可获得初相位
利用步骤2中获取的实际相移值
下面进行计算机模拟,设理想相移值为π/2、π,实际相移值为1.8rad和3.5rad,初相位
本发明提出一种用于结构光照明显微成像中结构光初相位的提取方法。该方法不受任何前提条件限制,且无需相等且精确的相移,大大降低了对成像系统性能的要求,为结构光照明成像提供了一种更普适、更方便的初相位提取方法。
机译: 显微镜照明装置具有用于将源光耦合到光透射元件(例如光导元件)中的装置。一种微结构玻璃纤维,可影响耦合到元件中的光的相位和幅度
机译: 用于在各种光亮状态下对样品进行成像的光学显微镜的测试样品装置,具有纳米结构,使得光学远场在表面的变化位置处解耦并产生不同的光照条件
机译: 用于控制高分辨率nmr的成像质量的测试结构-医学和生物学应用以及材料或材料研究中的成像位置末端(nmr-显微镜)及其制备方法