一种基于傅里叶叠层成像技术的显微镜错位校正方法与流程

本发明涉及计算领域，尤其是涉及了一种基于傅里叶叠层成像技术的显微镜错位校正方法。

背景技术

傅立叶叠层成像显微镜(fpm)是一种新开发的计算成像技术，可以提供高分辨率(hr)和宽视场(fov)的十亿像素的图像。傅里叶叠层成像技术可以应用于光学显微领域，获得高分辨率的显微图像，也可以用于生物医学领域，对病毒等微生物以及人体器官进行精确的研究治疗；在生命科学领域，利用傅里叶叠层成像技术可以更加方便地对微小的生物组织进行深入的科学研究。然而，以往的傅里叶叠层成像位置修正方法会增加算法复杂度、计算机负载和计算时间。

本发明中提出了一种基于傅里叶叠层成像技术的显微镜错位校正方法，先初始化重建的高分辨率物体和光瞳函数，然后依次点亮所有的led，分别生成低分辨率图像，并施加强度约束，接着更新物体函数和光瞳函数，重复以上步骤直到算法收敛，最后直接校正led阵列的全局偏移直到搜索全局偏移的成本函数达到最小值。本发明能够快速有效地校正显微镜成像的错位，同时使用的算法简单实用。

技术实现要素：

针对以往的傅里叶叠层成像位置修正方法会增加算法复杂度、计算机负载和计算时间的问题，本发明的目的在于提供一种基于傅里叶叠层成像技术的显微镜错位校正方法，先初始化重建的高分辨率物体和光瞳函数，然后依次点亮所有的led，分别生成低分辨率图像，并施加强度约束，接着更新物体函数和光瞳函数，重复以上步骤直到算法收敛，最后直接校正led阵列的全局偏移直到搜索全局偏移的成本函数达到最小值。

为解决上述问题，本发明提供一种基于傅里叶叠层成像技术的显微镜错位校正方法，其主要内容包括：

(一)全局偏移模型；

(二)高分辨率图像的重建；

(三)错位校正。

其中，所述的全局偏移模型，是一个只具有两个因子的简化位置模型，x和y是确定每个发光二极管(led)的位置的两个因子，(xm,n,ym,n)代表第m行第n列的led的位置。

其中，所述的高分辨率图像的重建，一共分为六个步骤：

(1)初始化重建的高分辨率(hr)物体oj(kx,ky)和光瞳函数pj(kx,ky)；

(2)生成低分辨率(lr)图像；

(3)施加强度约束；

(4)更新物体和光瞳函数；

(5)对所有的led重复步骤(2)到(4)；

(6)重复步骤(2)到(5)直到算法收敛。

进一步地，所述的生成低分辨率图像，是指用以下公式生成第m行第n列led对应的低分辨率图像：

其中表示用第m行第n列led照射获得的lr图像的傅里叶频谱。

进一步地，所述的施加强度约束，是用以下公式对捕获的图像施加强度约束：

其中和分别是lr图像带强度约束和不带强度约束的复数场。

进一步地，所述的更新物体和光瞳函数，是指用以下公式更新物体和光瞳函数：

其中，δ1和δ2是两个正则化常数，用于保证数量的稳定性。

其中，所述的错位校正，是指直接校正led阵列的全局偏移，直到全局偏移的成本函数达到最小，并且为了提高效率，只计算在目标的照明领域内的入射向量。

进一步地，所述的led阵列，每个led的入射波矢量为：

其中考虑到两个相邻led的间距是d，全局偏移的范围δx和δy设为[-d,d]，δx和δy初始化的值都是0mm。

进一步地，所述的成本函数，是指判断搜索全局偏移δx和δy的成本函数，每次重建高分辨率图像后判断成本函数是否达到最小值，如果没有则再次搜索。

进一步地，所述的搜索全局偏移，是指使用模拟退火(sa)算法来在一个循环中搜索δx和δy，这个循环一直持续到成本函数达到最小值为止，然后在校正led阵列的全局偏移后，通过传统的傅立叶叠层成像算法，利用所有的lr图像重建样本的hr图像；在某些情况下，校正led阵列的全局偏移后，需要使用常规sa算法的一次迭代来修正每个led固有的局部细微位置误差，最后，可以消除由led错位引起的重构的hr振幅和相位的退化。

附图说明

图1是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的系统流程图。

图2是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的成像原理图。

图3是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的全局偏移模型图。

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互结合，下面结合附图和具体实施例对本发明作进一步详细说明。

图1是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的系统流程图。主要包括全局偏移模型、高分辨率图像的重建和错位校正。

其中，高分辨率图像的重建一共分为六个步骤：

(1)初始化重建的高分辨率(hr)物体oj(kx,ky)和光瞳函数pj(kx,ky)；

(2)生成低分辨率(lr)图像；

(3)施加强度约束；

(4)更新物体和光瞳函数；

(5)对所有的led重复步骤(2)到(4)；

(6)重复步骤(2)到(5)直到算法收敛。

进一步地，所述的生成低分辨率图像，是指用以下公式生成第m行第n列led对应的低分辨率图像：

其中表示用第m行第n列led照射获得的lr图像的傅里叶频谱。

进一步地，所述的施加强度约束，是用以下公式对捕获的图像施加强度约束：

其中和分别是lr图像带强度约束和不带强度约束的复数场。

进一步地，所述的更新物体和光瞳函数，是指用以下公式更新物体和光瞳函数：

其中，δ1和δ2是两个正则化常数，用于保证数量的稳定性。

其中，错位校正是指直接校正led阵列的全局偏移，直到全局偏移的成本函数达到最小，并且为了提高效率，只计算在目标的照明领域内的入射向量。

其中每个led的入射波矢量为：

其中考虑到两个相邻led的间距是d，全局偏移的范围δx和δy设为[-d,d]，δx和δy初始化的值都是0mm。

进一步地，所述的成本函数，是指判断搜索全局偏移δx和δy的成本函数，每次重建高分辨率图像后判断成本函数是否达到最小值，如果没有则再次搜索。

进一步地，所述的搜索全局偏移，是指使用模拟退火(sa)算法来在一个循环中搜索δx和δy，这个循环一直持续到成本函数达到最小值为止，然后在校正led阵列的全局偏移后，通过传统的傅立叶叠层成像算法，利用所有的lr图像重建样本的hr图像；在某些情况下，校正led阵列的全局偏移后，需要使用常规sa算法的一次迭代来修正每个led固有的局部细微位置误差，最后，可以消除由led错位引起的重构的hr振幅和相位的退化。

图2是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的成像原理图。

傅里叶叠层成像技术的成像原理是按顺序开启led，捕获物体在不同照明角度下的低分辨率图像，重构物体的高分辨率图像。其中在空间域中，利用低分辨率图像的强度测量值作为约束条件来保证求解的收敛，在傅里叶域中，将目标透镜的约束相干传递函数作为支撑约束。经过多次迭代后，可以得到物体的高分辨率复数域和光瞳函数。

图3是本发明一种基于傅里叶叠层成像技术的显微镜错位校正方法的全局偏移模型图。

全局偏移模型是一个只具有两个因子的简化位置模型，x和y是确定每个发光二极管(led)的位置的两个因子，δx和δy是全局偏移的大小，圆点代表没有偏移时led所在的位置，菱形表示偏移后led所在的位置。

对于本领域技术人员，本发明不限制于上述实施例的细节，在不背离本发明的精神和范围的情况下，能够以其他具体形式实现本发明。此外，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围，这些改进和变型也应视为本发明的保护范围。因此，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。