本发明属于图像处理领域，尤其涉及一种印刷凹版的网穴以及以柔性版和橡胶版为代表的凸版等印刷版面上面向立体结构网点三维图像重建的层面数据获取方法。
背景技术：
凹版上网穴(也称为网点)和以柔性版、橡胶版为代表的凸版上网点的立体结构是决定相应印刷版印刷适性的主要技术参数。在凹版及凸版的制版过程中需要对印刷版上的网点质量参数进行有效的测量。目前，在使用单目成像系统时对上述网点质量参数的测量与解算主要是使用计算机视觉技术实现对网点顶部面积、网点高度的测量和底部面积的测算。但对网点侧壁曲面上的几何属性参数的测量与解算，上述类型的系统是无能为力的，而是需要在立体结构网点三维图像重建的基础上进行相应参数的测量与解算。该三维图像重建的核心支撑技术之一就是高精度层面数据的获取。
技术实现要素：
本发明所要解决的技术问题是，提供一种面向立体结构网点三维图像重建的层面数据获取方法。对具有立体结构网点的单目断层图像数据进行分析，提取网点结构中清晰成像的图像区域并对其中的误检像素进行处理，本案由此产生。
本发明解决其技术问题所采用的技术方案是：提供一种面向立体结构网点三维图像重建的层面数据获取方法，包括以下步骤，步骤1，对采集得到的网点单目层面图像序列s中的各层面图像上的像素计算相关系数，得到相关系数图像序列sr；步骤2，对图像序列s中的各层面图像上的像素进行模糊系数的计算，得到模糊系数图像序列sf；步骤3，定义相关度和模糊度的判别阈限值；设计低通滤波器，分别对sr和sf进行滤波处理，得到对应的二值图像序列br和bf，将br和bf两个图像序列中对应的层面图像上诸像素进行二值乘法运算得到判别图像序列bd；步骤4，利用图像序列s和判别图像序列bd解算出清晰要素图像序列se；步骤5，输出图像序列se，即所获取的面向立体结构网点三维图像重建的层面数据。
按上述技术方案，对图像序列s进行8位灰度化处理：比较彩色像素的三通道色度信号值，提取其中的最大值作为灰度化值并赋给像素，实现图像的灰度化处理，处理完成之后需要以图像的像素灰度值对图像序列s进行重构。
按上述技术方案，对图像序列s进行假彩化处理：将像素灰度值扩展为三元组数据结构，其中，三元组的三个域分别表示假彩色rgb图像的三通道色度信号值；假彩化的三个通道中，红和绿通道为保留域，初始化为0，分别用于漏检和多检属性的标记；蓝通道的域为像素的灰度值，处理完成之后需要对原图像序列进行重构。
按上述技术方案，对灰度图像序列s进行增强处理：先使用高斯滤波器对灰度化图像序列s的诸层图像进行滤波处理得到ig；再使用数学形态学提取各层图像上的边缘得到二值边缘图像序列ie；将ig和ie进行图像融合得到增强图像序列ih，处理完成之后需要以增强图像序列ih对图像序列s进行重构。
按上述技术方案，对图像进行空间褶接处理，具体处理方法如下：分别计算待处理图像行和列上的绝对差分投影均值dr和dc，当dr不小于dc时，图像平面沿垂直方向卷接，再将得到的图像柱面沿水平方向对接；当dr<dc时，图像平面沿水平方向卷接，再将图像柱面沿垂直方向对接。
按上述技术方案，网点层面图像拍摄时出现的上层对下层的局部遮挡，导致下层面图像采集不完整的情况。需要在下层断面图像上进行判断与标记。具体处理方法如下：计算经灰度化和增强处理的输入图像序列s上诸层面图像中所有像素的聚焦系数得到聚焦系数图像序列if；采用平均先进法计算s序列的最大聚焦系数集合的平均先进聚焦系数faa，以faa为判据对se图像序列中的像素进行伪彩色标记，处理完毕后重构待输出图像序列se。
按上述技术方案，在获取层面图像时，景深区域相交引入的层面内的多检像素，需要对相应的像素数据进行判断与标记。具体处理方法如下：将判别图像序列bd中相邻层面图像进行乘法运算，得到标记图像序列bi，若bi为非零集合，首先判断se图像序列属性，如果其为灰度图像，则将其扩展为伪彩色图像得到si；如果已经扩展为伪彩色图像，则将其进行复制且将复制的数据集定义为si，再以序列的诸层面图像中像素值为1作为条件在bi中检索得到层面信息k-1和k，以及像素的位置信息，进而比较if图像序中相应的两层对应位置像素的聚焦系数值并进行标记，处理完毕后按下式重构待输出图像序列se。
本发明产生的有益效果是：本发明旨在对短景深单目光学成像系统在恒定的进给位移量下获取的由网点顶部至底部的层面图像序列进行处理，实现一种面向立体结构网点三维图像重建的层面数据获取的方法。本发明所获取的层面数据可以直接应用于具有立体结构网点的三维图像重建并为其几何结构的多参数直接测量提供支持。
附图说明
图1是面向立体结构网点三维图像重建的层面数据获取方法实施的流程图。
具体实施方式
为了使本发明的技术方案更加清楚，以下结合附图，以柔印版网点三维图像重建的层面数据获取为例对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
以拍摄得到的网点单目层面图像序列s为输入数据集，s＝{i1,i2,…,in}，n为序列中图像的幅数，亦即网点图像拍摄时层面的数量。判断输入s中图像为彩色还是灰度图，如果为彩色图像，则对其进行灰度化处理：设输入图像序列s的第k层图像上(i,j)位置为当前处理像素的红、绿、蓝三个通道的色彩分量分别定义为ik(r；i,j)、ik(g；i,j)、ik(b；i,j)，则按下式进行灰度化处理，如图1所述，
遍历s中各层图像上的诸像素并进行灰度化处理，处理完毕后重构输入图像序列s，s＝{ik←i'k|k＝1,2，…，n}。
因为拍摄网点外表面为非均匀漫反射面，且成像系统工作时引入的噪声，需要对灰度图像序列s进行增强处理。
先使用高斯滤波器对灰度化图像序列s的诸层图像进行噪声平滑处理得到ig；再使用数学形态学提取各层图像的边缘得到二值边缘图像序列ie；将ig和ie进行图像融合得到增强图像序列ih。任取第k层图像上(i,j)位置为当前处理像素(0<k<n+1)，
按式(2)遍历s上诸层图像进行上述增强处理，处理完毕后重构输入图像序列s，s＝{ik＝ihk|k＝1,2，…，n}。
建立相关系数图像序列sr，sr＝{ir1,ir2,…,irn}。以第k层图像上当前处理像素(i,j)为图像分块的角点像素，设置4个大小为k×k像素数的图像信息子块；对各个子块进行傅里叶变换，对变换结果取模的对数操作得到矩阵
式中，cov和σ分别为协方差和均方差算子。按下式像素(i,j)与在(2k-1)×(2k-1)域上的相关系数：
irk(i,j)＝(c1+c2)/2.(4)
遍历s中各层图像上的诸像素并计算相关系数得到相关系数图像序列sr。
建立模糊系数图像序列sf，sf＝{if1,if2,…,ifn}。以第k层图像上当前处理像素(i,j)为中心，设置m×m像素数的图像窗口；计算窗口内像素灰度值的均方差σl和当前层面图像的全局均方差σg；计算当前像素的模糊系数：
ifk(i,j)＝|σl-σg|/σg.(5)
遍历s中各层图像上的诸像素，计算模糊系数得到模糊系数图像序列sf。
计算相关系数图像序列sr的方差，按照西电准则(westernelectricrules)定义相关度判别阈限值为hr。类似，利用糊系数图像序列sf计算模糊度的判别阈限值hf；设计滤波器lpfr和lpff：
按式(6-7)，对sr和sf各层图像上的诸像素进行滤波处理，得到对应的二值图像序列br，br＝{br1,br2,…,brn}和bf，bf＝{bf1,bf2,…,bfn}。
将br和bf两个图像序列中对应的层面图像上诸像素进行二值乘法运算，得到判别图像序列bd，bd＝{bd1,bd2,…,bdn}。设bd上第k层图像上当前处理像素(i,j)，则
bdk(i,j)＝brk(i,j)×bfk(i,j).(8)
进一步，利用图像序列s和判别图像序列bd解算出清晰要素图像序列se，se＝{ie1,ie2,…,ien}，
按式(9)，图像序列se的层面图像中，清晰区域的像素灰度值为相应s中对应层面图像中像素的灰度值，模糊区域的像素灰度值置0。
网点层面图像拍摄时出现的上层对下层的局部遮挡，导致下层面图像采集不完整的情况。需要在下层断面图像上进行判断与标记。
定义像素的聚焦系数为该像素8邻域内所有像素的均方差σ。计算s上诸层面图像中所有像素的聚焦系数得到聚焦系数图像序列if，并求取if中各层面图像内的最大值并将其记为该层面图像的最大聚焦系数
式(10)中，max、average、mean函数分别为求取最大值、平均值和中位值的函数。
为了对像素的特殊属性，如漏检、多检等属性进行标记，将灰度图像扩展到三通道的假彩色rgb图像。设图像i在(i,j)位置处的像素灰度值为i(i,j)，扩展为三元组数据结构ix(r,g,b；i,j)。其中，三元组的域ix(r；i,j)、ix(g；i,j)、ix(b；i,j)分别表示假彩色rgb图像的三通道色度信号值。按下式将灰度图像扩展到假彩色rgb图像：
三个通道中，r和g为保留域，分别用于漏检和多检属性的标记。
将图像序列se的灰度层面图像按式(11)扩展为假彩色rgb图像，得到扩展的要素图像序列so，so＝{io1,io2,…,ion}，计算平均先进聚焦系数faa并按下式对层面间遮挡而引入的漏检像素进行标记：
式(12)中，iek∈se，iok∈so，ifk∈if，k＝1,2,…,n。按式(12)遍历se诸层面图像的所有像素，以faa为判据对se图像序列中的像素进行伪彩色标记。处理完毕后重构待输出图像序列se，se＝{iek←iok|k＝1,2，…，n}。
在获取层面图像时，考虑到景深区域相交可能会引入的层面内的多检像素，需要对相应的像素数据进行判断与标记。
将判别图像序列bd中相邻层面图像进行乘法运算，得到标记图像序列bi，bi＝{bi1,bi2,…,bin}：
若bi为非零集合，首先检查se图像序列属性，如果其为灰度图像，则按式(11)将其扩展为伪彩色图像得到si，si＝{ii1,ii2,…,iin}；如果已经扩展为伪彩色图像，则将其进行复制且将复制的数据集定义为si。再bi序列的诸层面图像中以像素值为1作为条件在bi中检索得到层面信息k-1和k，以及像素的位置信息(i,j)，进而比较if图像序中相应的两层对应位置像素的聚焦系数值并按下式进行标记。
式(14)中，iek∈se，iik∈si，ifk∈if，k＝2,…,n。按式(14)遍历bi中值为1的像素位置对应的if序列的相应层面图像上的所有对应像素，并进行伪彩色标记。处理完毕后按下式重构待输出图像序列se。
输出图像序列se，即所获取的面向立体结构网点三维图像重建的层面数据。
在对图像像素进行邻域上的变换处理时，图像边界像素存在邻域像素缺失的情况，需要对图像进行空间褶接处理。
设图像i的大小为m×n(像素)，分别计算行和列上的绝对差分投影均值dr和dc：
式(16-17)中，r和c分别表示当前层面图像的行和列的序号。按下述逻辑进行空间褶接处理：
(1)当dr>dc时，i(1,·)的后继行为i(m,·)，i(m,·)的前驱行为i(1,·)，其余像素的空间前驱-后继维持关系不变。进一步，将图像柱面沿水平方向对接，i(·,1)的后继列为i(·,n)，i(·,n)的前驱列为i(·,1)，其余像素的空间前驱-后继维持关系不变；当dr＝dc时，水平方向优先，执行本逻辑。
(2)当dr<dc时，图像平面沿水平方向卷接。进一步，将图像柱面沿垂直方向对接。
本案所发明的技术手段也适用于丝网印刷品上的网点，以及工业及医学等领域中涉及微观立体结构的多层面三维图像重建的层面数据获取方法。
应当理解的是，对本领域技术人员来说，可以根据上述说明加以变换或改进，而所有这些变换和改进都应属于本发明所附权利要求的保护范围。