一种基于优化结构探测函数的显微成像方法与流程

本发明涉及结构探测显微成像技术领域，具体涉及一种优化结构探测函数的方法。

背景技术：

光学显微术是一种历史悠久且十分重要的无破坏性技术，被广泛应用于生物和材料科学等领域。

2009年，美国哈佛大学谢晓亮、jeffw.lichtman课题组在nanoletters杂志上发表论文，提出扫描图案探测显微技术，通过将探测调制和非解扫单元探测器完成的空间累积成像相结合来实现空间调制。仿真分析得出，在荧光非相干成像情况下，横向分辨率均可达到普通显微系统的2倍。

专利号：zl201510868029.6描述的一种超分辨结构探测共焦相干成像装置及其成像方法，其将结构探测成像方法与共焦相干显微系统相结合,提高了共焦相干成像系统的空间截止频率,拓宽空间频域带宽,从而改善了成像系统横向分辨力，但其不能充分提高成像分辨率，实现超分辨成像。

技术实现要素：

本发明的目的是：针对现有技术的缺陷。提供一种基于优化结构探测函数的显微成像方法，充分提高成像分辨力，实现超分辨成像。

本发明采用如下技术方案实现：一种基于优化结构探测函数的显微成像方法，包括以下步骤：

步骤一、在扫描显微成像系统中，将扫描标准样品采样点得到的探测光斑光强分布与随机生成的结构探测函数相乘，得到一个新的图像；

步骤二、将步骤一中得到的图像中的每个像素的灰度值求和后得到灰度和值；

步骤三、重复步骤一和步骤二求得标准样品所有采样点的灰度和值；

步骤四、将步骤三中标准样品所有采样点的数值归一化至0—255区间内，作为重构图像对应采样点的光强灰度值，得到重构图像；

步骤五、将重构图像与标准样品逐点比较获得每一像素点的灰度值误差；

步骤六、将步骤五中重构图像每一像素点的误差的平方作为基准误差，将基准误差求和作为在当前结构探测函数下的误差，将当前结构探测函数下的误差进行反向传递至结构探测函数部分，调整结构探测函数；

步骤七、通过迭代训练，得到最优结构探测函数；然后将最优结构探测函数带入步骤一中替换随机生成的结构探测函数，并经过步骤二、步骤三和步骤四得到重构图像，即超分辨图像。

本发明具有如下有益效果：本发明通过利用误差向结构探测函数的反向传播，不断地通过误差对当前结构探测函数的梯度调整结构探测函数，降低重构图像与标准样品的误差。从而求解结构探测函数最优解，通过图像重构得到超分辨图像，提高了基于结构探测的显微成像分辨率，实现超分辨成像。

附图说明

图1是本发明的原理框图。

具体实施方式

具体实施方式一：下面结合图1具体说明本实施方式。本实施方式一种基于优化结构探测函数的显微成像方法，包括以下步骤：

步骤一、在扫描显微成像系统中，将扫描标准样品采样点得到的探测光斑光强分布与随机生成的结构探测函数相乘，得到一个新的图像；

步骤二、将步骤一中得到的图像中的每个像素的灰度值求和后得到灰度和值；

步骤三、重复步骤一和步骤二求得标准样品所有采样点的灰度和值；

步骤四、将步骤三中标准样品所有采样点的数值归一化至0—255区间内，作为重构图像对应采样点的光强灰度值，得到重构图像；

步骤五、将重构图像与标准样品逐点比较获得每一像素点的灰度值误差；

步骤六、将步骤五中重构图像每一像素点的误差的平方作为基准误差，将基准误差求和作为在当前结构探测函数下的误差，将当前结构探测函数下的误差进行反向传递至结构探测函数部分，调整结构探测函数；

步骤七、通过迭代训练，得到最优结构探测函数；然后将最优结构探测函数带入步骤一中替换随机生成的结构探测函数，并经过步骤二、步骤三和步骤四得到重构图像，即超分辨图像。

通过将标准样品的重构图像与标准样品逐点比较获得的误差反向传播，实现对结构探测函数的调整，从而实现最优结构探测函数的求取。在激光扫描显微成像系统中，将扫描标准样品得到的探测光斑光强分布与结构探测函数相乘后再求和，将重构图像与标准样品逐点比较获得每一点的灰度值误差，将整幅图像每一点的误差的平方求和作为在当前结构探测函数下的误差，并将误差进行反向传递至结构探测函数部分，据此调整结构探测函数从而降低重构图像与标准样品的误差，通过迭代训练，最终实现最优结构探测函数的求取。进而能用最优结构探测函数实现更高分辨率的结构探测超分辨成像。

每个扫描标准样品采样点得到的探测光斑光强通过步骤一和步骤二得到其对应的数值，所有采样点的数值重新排列得到重构图像的排列顺序与对标准样品的扫描顺序是一致的。

具体实施方式二：本实施方式是对具体实施方式一的进一步说明，本实施方式中所述步骤一中，假设在扫描显微成像系统中，光斑的图像采集区域为m行n列的像素区域，标准样品在显微成像过程中，扫描显微成像系统实现对标准样品的u×v个点的采样；在采样点ij处获得的探测光斑图像的光强分布为s^ij，在m×n的像素区域中，第x行y列的像素的光强为该像素所对应的结构探测函数为mxy，探测光斑的光强分布与结构探测函数相乘后再求和可得重构图像在采样点ij处的光强imageij，因此图像的重构由如下公式实现；

imageij为重构图像在采样点ij处的光强灰度值，1≤x≤m,1≤y≤n。

具体实施方式三：本实施方式是对具体实施方式二的进一步说明，本实施方式中在所述步骤五中，重构图像与标准样品在采样点ij处的基准误差errorij可由如下公式求得，因为灰度值的差值有正有负，因此取平方避免抵消；

objectij为标准样品在采样点ij处的光强灰度值的真值，errorij是重构图像与标准样品在采样点ij处的基准误差。

具体实施方式四：本实施方式是对具体实施方式三的进一步说明，本实施方式中对于标准样品u×v的采样区域，当前结构探测函数下的误差error由如下公式求得：

1≤i≤u,1≤j≤v。

具体实施方式五：本实施方式是对具体实施方式四的进一步说明，本实施方式中所述步骤六中，将求得的总误差error反向传递至结构探测函数，据此调整结构探测函数总误差error对当前结构探测函数的梯度，即结构探测函数的调整量为：

具体实施方式六：本实施方式是对具体实施方式五的进一步说明，本实施方式中依据反向的误差，调整后的结构探测函数通过如下公式实现；

其中为第t次调整时坐标点x,y处的像素点采用的结构探测函数，为第t次调整后坐标点x,y处的像素点的结构探测函数，r是用来控制调整幅度大小的正整数比例因子。

虽然本发明所揭示的实施方式如上，但其内容只是为了便于理解本发明的技术方案而采用的实施方式，并非用于限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭示的核心技术方案的前提下，可以在实施的形式和细节上做任何修改与变化，但本发明所限定的保护范围，仍须以所附的权利要求书限定的范围为准。