一种基于多角度照明反卷积的高分辨率显微成像方法技术

发明专利技术公开了一种基于多角度照明反卷积的高分辨率显微成像方法。与传统明场显微镜采用的科勒照明方式不同，本系统采用LED阵列作为照明光源，通过依次采集不同照明角度下的物体图像，获取不同空间频率成分的物体信息，然后根据系统参数计算出不同照明角度对应的强度传递函数对图像进行最小二乘反卷积重构。相比传统明场显微镜同时获取所有空间频率信息的成像方法，本发明专利技术通过使用最小二乘反卷积的算法，对噪声具有更强的鲁棒性。对噪声具有更强的鲁棒性。对噪声具有更强的鲁棒性。

【技术实现步骤摘要】
一种基于多角度照明反卷积的高分辨率显微成像方法


[0001]本专利技术属于光学显微成像技术，特别是一种基于多角度照明反卷积的高分辨率显微成像方法。

技术介绍

[0002]自从Leeuwenhoek于三百多年前第一次使用手工自制显微镜观察细胞以来，人类在显微成像领域陆续开展了诸多探索。其中，以看的更清楚，更精细为目的的高分辨率显微成像是其中重要的一环。1873年，Abbe在德国显微镜学报上发表了显微镜衍射成像理论，对显微镜的成像过程做了科学的阐释，首次提出了显微镜的光学分辨率是有限的，受到物镜数值的限制(张祥翔.现代显微成像技术综述[J].光学仪器,2015(6):550
‑
560.)。
[0003]在传统科勒照明下的明场显微成像中，由于衍射光线无法全部通过物镜，所以物空间的一个点经显微镜成像后会映射成一个衍射斑，导致分辨率的降低。在这里，系统的光学分辨率的计算方法常常采用瑞利判据，它规定当一点的衍射极大与另一点的第一衍射极小相重合时，则刚好能分辨两点。(Katsumasa Fujita,Follow
‑
up review:recent progress in the development of super
‑
resolution optical microscopy[J],Microscopy,Volume 65,Issue 4,August 2016,275
–
281)在瑞利判据下，显微镜的光学分辨率f
R
为：/>[0004][0005]综上所述，在照明光一定的情况下，显微镜的光学分辨率由物镜的数值孔径决定，所以要想对物体看的越精细，就需要越大数值孔径的物镜。然而瑞利判据下的光学分辨率截止频率只是人眼直接观察的分辨率极限，并不是光学成像系统的理论分辨率极限。光学系统的理论分辨率极限f(Chao Zuo,Transport of intensity equation:a tutorial[J],Optics and Lasers in Engineering,Volume 135,2020,106
‑
187,)为：
[0006][0007]其中，NA
ill
为照明光源的数值孔径，当物镜数值孔径NA
obj
较小时，NA
ill
通常与NA
obj
相等，此时光学系统的理论分辨率极限f为：
[0008][0009]但是光学成像系统往往并不能达到其理论分辨率极限，原因是光学成像系统对信息的传递能力随着频率的升高而衰减，到截止频率附近衰减到0，会导致一定范围内的高频信息被淹没在噪声中。所以光学成像系统的实际分辨率难以达到理论分辨率极限的问题亟待解决。这里我们提出了一种基于多角度扫描的高分辨率显微成像方法，可以解决上述问题。

技术实现思路

[0010]本专利技术的目的在于解决了传统科勒照明下，因截止频率处信噪比过低而无法达到衍射极限分辨率的问题。
[0011]本专利技术的技术方案如下：
[0012]一种基于多角度照明反卷积的高分辨率显微成像方法，步骤如下：
[0013]步骤一.原始图像采集，使用可编程LED阵列作为成像系统的照明光源，依次点亮LED阵列的LED单元，由相机在同步触发信号下采集一系列不同照明角度下的图像序列I
m,n
(x,y)，其中，I是光强图像，(x,y)是空域坐标，m,n是光强图像所对应的LED单元在LED阵列中的位置坐标；
[0014]步骤二.计算出每个LED单元的所对应的照明光空间频率矢量(u
m
,u
n
)；
[0015]步骤三.计算LED阵列坐标为(m,n)的LED单元所对应的相干强度传递函数ATF
m,n
；
[0016]步骤四.对步骤一中拍到的图像序列进行最小二乘反卷积，首先对拍到的强度图序列I
m,n
(x,y)做傅里叶变换得到强度图频谱序列(u,v)是频域坐标，然后进行最小二乘反卷积获得重构图像I
deconv
(x,y)。
[0017]优选的，步骤二中，每个LED单元所对应的照明光空间频率矢量(u
m
,u
n
)的计算公式为：
[0018][0019][0020]其中，m为LED阵列X方向的第m个LED单元，n为LED阵列Y方向的第n个LED单元，h为可编程LED阵列距离待测样品的距离，d为LED单元的间距，λ为照明波长。
[0021]优选的，步骤三中，相干强度传递函数ATF
m,n
的计算公式为：
[0022][0023]其中(u,v)是频域坐标，为物镜相干照明截止频率,NA
obj
是物镜数值孔径。
[0024]优选的，步骤四中，最小二乘反卷积公式为：
[0025][0026]其中，I
deconv
(x,y)为反卷积后图像的强度图，F
‑1傅里叶逆变换。
[0027]本专利技术的有益效果是：(1)使用可编程LED阵列使用角度扫描的方法拍摄多角度照明下待测物体的光强图序列，可以对非相干照明中的不同空间频率成分进行高信噪比，高
动态范围的信息采集，具有成像分辨率更高，成像效果更好的特点。(2)基于系统实际参数，对不同照明角度下的光强传递函数进行了计算，然后基于最小二乘法对所拍到的光强图序列进行反卷积重构，所得结果更精确，对噪声的鲁棒性更强。
附图说明
[0028]图1是基于多角度照明反卷积的高分辨率显微成像方法的流程图
[0029]图2是不同角度照明下的强度图像、强度图像频谱和与之对应的相干强度传递函数ATF
m,n
。
[0030]图3是使用该方法的合成强度传递函数和传统非相干照明强度传递函数对比。
[0031]图4是使用该方法的成像结果、使用传统非相干成像结果的对比。
具体实施方式
[0032]为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本申请，并不用于限定本申请。
[0033]结合图1，本实施例是一种实现基于多角度照明反卷积的高分辨率显微成像方法，包括以下步骤：
[0034]步骤一：多角度照明下的强度图像采集。使用可编程LED阵列作为系统的照明光源，依次点亮LED单元，LED阵列的中心单元须对准成像系统光轴，LED单元每次发出不同角度的单色或白色近似平面波照射在样品上。同时单色或彩色相机接收来自可编程LED阵列的同步触发信号后开始曝光，对多角度照明下的待测物体所成像进行采集，形成单色或彩色光强图像序列I
m,n
(x,y)。其中，I(x,y)是光强图像，m,n是光强图像所对应的LED单元在LED阵列中的位置坐标。...

【技术保护点】

【技术特征摘要】
1.一种基于多角度照明反卷积的高分辨率显微成像方法，其特征在于步骤如下：步骤一.原始图像采集，使用可编程LED阵列作为成像系统的照明光源，依次点亮LED阵列的LED单元，由相机在同步触发信号下采集一系列不同照明角度下的图像序列I
m,n
(x,y)，其中，I是光强图像，(x,y)是空域坐标，m,n是光强图像所对应的LED单元在LED阵列中的位置坐标；步骤二.计算出每个LED单元的所对应的照明光空间频率矢量(u
m
,u
n
)；步骤三.计算LED阵列坐标为(m,n)的LED单元所对应的相干强度传递函数ATF
m,n
；步骤四.对步骤一中拍到的图像序列进行最小二乘反卷积，首先对拍到的强度图序列I
m,n
(x,y)做傅里叶变换得到强度图频谱序列(u,v)是频域坐标，然后进行最小二乘反卷积获得重构图像I
deconv
(x,y)。2.根据权利...

【专利技术属性】
技术研发人员：张晓磊，左超，孙佳嵩，胡岩，沈德同，尹维，
申请(专利权)人：南京理工大学智能计算成像研究院有限公司，
类型：发明
国别省市：