一种基于马尔科夫随机场约束的显微三维重建方法技术

本发明专利技术公开了一种基于马尔科夫随机场约束的显微三维重建方法：用显微镜对物体扫描拍摄得到多张不同聚焦程度的图像，组成多聚焦图像序列；对图像序列中每张图像进行非下采样轮廓波变换，并构建金字塔形分布的低通滤波器组对非下采样轮廓波变换的结果进行处理；通过比较不同图像处理结果得到隐含在图像序列中的深度信息和全焦信息；利用马尔可夫随机场的先验约束对深度信息进行修正；利用获取的多聚焦图像序列的深度信息，构建显微图像的三维形态。该算法实现原理简单，可以有效地应用于精细结构检测、超精密加工以及医疗手术等领域。该方法解决了现有技术中聚焦评价函数难以平衡检测灵敏性和对噪声鲁棒性的问题，以及提高了弱纹理样本的重建能力。了弱纹理样本的重建能力。了弱纹理样本的重建能力。

【技术实现步骤摘要】
一种基于马尔科夫随机场约束的显微三维重建方法


[0001]本专利技术涉及数字图像处理、计算机视觉等
，特别涉及一种基于马尔科夫随机场约束的显微三维重建方法。

技术介绍

[0002]精密加工技术的发展对显微三维重建技术提出了更高的要求。三维重建技术的关键在于准确的提取显微样本的深度信息和全焦信息。在这种情况下，由于聚焦法测量范围大、对待测样本表面光滑程度的要求较低，成为了在精密加工领域内进行质量检测的可靠方法，尤其是对切削刀具等具有大倾角斜面样本的检测。
[0003]尽管聚焦法具有快速获取样本表面深度信息和全焦信息的能力，但其关键技术——聚焦评价函数仍存在一些固有问题，从而严重限制了显微三维重建算法应用范围的扩展。第一，为了达到精确的效果，聚焦评价函数需要对高频信息敏感，但同时还需要对噪声保持鲁棒性，这是很难平衡的。第二，聚焦评价函数表现的上限仍然受样本表面纹理信息丰富性的影响，对表面光滑且纹理弱的样本很难进行准确的聚焦检测，进而限制了三维重建的精度。目前现有技术中多以提高聚焦评价函数的敏锐性和鲁棒性为目的，例如ydin T,Akgul Y S.A New Adaptive Focus Measure for Shape From Focus[C]//BMVC.2008:1
‑
10.；以及公开号为CN108036739B、CN110443882B、CN111624756A和CN111932677A等中国专利公开文本。
[0004]综上所述，可以发现，尽管与其他三维重建算法相比，聚焦法在显微三维重建领域存在一定优势，但其重建算法仍有待提高。

技术实现思路

[0005]本专利技术的目的在于提供一种基于马尔科夫随机场约束的显微三维重建方法，该方法解决了现有技术中聚焦评价函数难以平衡检测灵敏性和对噪声鲁棒性的问题，以及提高了弱纹理样本的重建能力。
[0006]为实现上述专利技术目的，本专利技术提供的技术方案为：
[0007]一种基于马尔可夫随机场约束的显微三维重建算法，包括以下步骤：
[0008](1.1)利用显微镜和扫描平台搭建显微三维重建系统对待测样本进行纵向扫描成像，得到反映样本表面不同深度区域聚焦情况的多聚焦图像序列；
[0009](1.2)、借助基于非下采样轮廓波变换的聚焦评价函数对多聚焦图像序列进行处理，得到与多聚焦图像序列对应的反映图像聚焦区域的掩模版序列；
[0010](1.3)借助掩模版序列得到样本表面初始深度信息和全焦信息；
[0011](1.4)借助马尔科夫随机场和初始深度信息构建能量函数，通过对能量函数进行求解，得到深度信息的优化结果；
[0012](1.5)利用提取的全焦信息和优化的深度信息进行三维重建。
[0013]其中，在(1.1)中，借助显微镜小景深的成像特点，通过等间距移动显微镜载物台Z
轴对待测样本进行纵向扫描成像。
[0014]优选的，所述借助基于非下采样轮廓波变换的聚焦评价函数对多聚焦图像序列进行处理，得到与多聚焦图像序列对应的反映图像聚焦区域的掩模版序列的具体方法包括：
[0015]1，将多聚焦图像转换为单通道图像，并借助非下采样轮廓波变换对聚焦图像序列中的每张图像进行变换；
[0016]2，对于每张图像变换后的结果，利用金字塔形低通滤波器组进行处理，得到一个相同大小的矩阵；
[0017]3，为多聚焦图像序列中每张图像构建相同大小的、初始值为0的掩模版图像，构成掩模版序列；
[0018]4，对图像之间处理结果(2中得到的矩阵)逐像素进行比较，记录比较结果中最值所在的像素位置和通道数目，并在掩模版序列中修改相应的数值为1，以更新掩模版序列。
[0019]优选的，所述将多聚焦图像转换为单通道图像，并借助非下采样轮廓波变换对图像序列中的每张图像进行变换的具体方法包括：
[0020]对多聚焦图像序列中的每张图像进行非下采样轮廓波变换，提取3个频段的高频信息，在低频段到高频段依次构建2、3、4级方向滤波器组进行方向滤波，以获得28个矩阵表示不同频段、不同方向的纹理分布。
[0021]优选的，所述对于每张图像变换后的结果，利用金字塔形低通滤波器组进行处理，得到一个相同大小的矩阵的具体方法包括：
[0022]1，按照均值μ＝0，方差σ＝1的参数设置，从低频段到高频段分别构建尺寸为21、13、5的高斯滤波器，形成金字塔形低通滤波器组；
[0023]2，对不同尺度下的方向系数矩阵进行绝对值求和处理，每次加上一个系数矩阵的绝对值后，用对应尺度的低通滤波器组对结果进行处理，直到该尺度下所有方向系数矩阵相加完成，以得到不同尺度的处理结果；
[0024]3，对不同尺度的处理结果求和，以得到一个反映图像中每个像素聚焦程度的矩阵。
[0025]优选的，所述借助掩模版序列得到样本表面初始深度信息和全焦信息包括：
[0026]1，将不同图像的掩模版与对应的通道相乘，并将相乘的结果相加，以得到基于聚焦评价算法估计到的初始深度信息；
[0027]2，将不同的图像与对应的掩模版相乘，并将相乘的结果相加，以得到反映样本表面纹理信息的全焦图像。
[0028]优选的，所述借助马尔科夫随机场和初始深度信息构建能量函数，通过对能量函数进行求解，得到深度信息的优化结果包括：
[0029]1，基于最大后验概率的能量函数框架，代入初始深度信息作为能量函数的保真项，代入样本表面的马尔可夫随机场作为能量函数中的惩罚项；
[0030]2，将能量函数表示为无向图的形式，借助以α拓展图割算法为核心的迭代优化算法对无向图进行求解，以得到修正后的深度信息。
[0031]优选的，所述将能量函数表示为无向图的形式，借助以α拓展图割算法为核心的迭代优化算法对无向图进行求解，以得到修正后的深度信息包括：
[0032]1，以像素与八邻接像素之间差的绝对值作为无向图中的n
‑
links，以当前深度信
息与当前通道之间差的绝对值作为无向图中的t
‑
links，并将无向图的能量作为初始的最小割能量；
[0033]2，基于α拓展图割对当前无向图进行求解，以得到无向图的最小割；
[0034]3，计算最小割能量是否降低，如果是，以当前最小割重新构建无向图的t
‑
links，并返回步骤2，否则将上次迭代的最小割作为优化结果输出。
[0035]优选的，所述用提取的全焦信息和优化的深度信息进行三维重建，是以全焦信息作为三维图形的底面和颜色，以深度信息作为三维图形的高度，构建三维图形。
[0036]与现有技术相比，本专利技术具有以下技术效果：
[0037]利用非下采样轮廓波变换与金字塔形低通滤波器结合的方法进行聚焦评价，不仅发挥了非下采样轮廓波变换多尺度、多方向提取纹理信息的优势，而且针对不同尺度的变换结果，进行不同程度的低通处理，很好地平衡了聚焦评价函数的敏本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于马尔科夫随机场约束的显微三维重建方法，其特征在于，包括以下步骤：(1.1)利用显微镜和扫描平台搭建显微三维重建系统对待测样本进行纵向扫描成像，得到反映样本表面不同深度区域聚焦情况的多聚焦图像序列；(1.2)借助基于非下采样轮廓波变换的聚焦评价函数对多聚焦图像序列进行处理，得到与多聚焦图像序列对应的反映图像聚焦区域的掩模版序列；(1.3)借助掩模版序列得到样本表面初始深度信息和全焦信息；(1.4)借助马尔科夫随机场和初始深度信息构建能量函数，通过对能量函数进行求解，得到深度信息的优化结果；(1.5)利用提取的全焦信息和优化的深度信息进行三维重建。2.根据权利要求1所述的基于马尔科夫随机场约束的显微三维重建方法，其特征在于，所述(1.2)通过以下子步骤来实现：(2.1)将多聚焦图像转换为单通道图像，并借助非下采样轮廓波变换对多聚焦图像序列中的每张图像进行变换；(2.2)对于每张图像变换后的结果，利用金字塔形低通滤波器组进行处理，得到一个相同大小的矩阵；(2.3)为多聚焦图像序列中每张图像构建相同大小的、初始值为0的掩模版图像，构成掩模版序列；(2.4)对图像之间处理结果逐像素进行比较，记录比较结果中最值所在的像素位置和通道数目，并在掩模版序列中修改相应的数值为1，以更新掩模版序列。3.根据权利要求2所述的基于马尔科夫随机场约束的显微三维重建方法，其特征在于，所述步骤(2.1)通过如下方式实现：对多聚焦图像序列中的每张图像进行非下采样轮廓波变换，提取3个频段的高频信息，在低频段到高频段依次构建2、3、4级方向滤波器组进行方向滤波，以获得28个矩阵表示不同频段、不同方向的纹理分布情况。4.根据权利要求2所述的基于马尔科夫随机场约束的显微三维重建方法，其特征在于，所述步骤(2.2)通过以下子步骤来实现：(2.2.1)按照均值μ＝0，方差σ＝1的参数设置，从低频段到高频段分别构建尺寸为21、13、5的高斯滤波器，形成金字塔形低通滤波器组；(2.2.2)对不同尺度下的方向系数矩阵进行绝对值求...

【专利技术属性】
技术研发人员：尚明皓，周海洋，余飞鸿，
申请(专利权)人：杭州图谱光电科技有限公司，
类型：发明
国别省市：