对具有条形结构的图像进行滤波的方法技术

一种借助伽柏滤波器对具有条形结构的图像进行滤波的方法，这些伽柏滤波器在空间域内是由沿主方向叠加了余弦函数的二维高斯钟形曲线形成的，其特征在于，将图像分成多个块，对于各个块确定条形结构的主导方向，并且以这样一种方式进行滤波：一次旋转一个块，直到主导方向与伽柏滤波器的主方向垂直，在主方向上进行一次滤波，而在与主方向垂直的方向上进行另一次滤波，并且将经过滤波的块再次旋转回来。（*该技术在2023年保护过期，可自由使用*）

【技术实现步骤摘要】
【国外来华专利技术】
本专利技术涉及一种借助伽柏(Gabor)滤波器，这些伽柏滤波器是由沿主方向叠加了余弦函数的二维高斯钟形曲线在空间域中形成的。本专利技术的重要应用领域涉及皮肤印痕图像的处理，尤其是指纹图像。从指纹的灰度图像(由各种各样可采用商业手段得到的传感器提供)开始，将要实现所包含的凹槽结构的强化，同时实现任何所包括的缺陷的修复。在进一步的步骤中，将这一经过处理的图像二值化，然后从有蚀痕的二值图像中提取出决定性特征。伽柏滤波器已经证明了它们在(局部)凹槽结构的强化和修复方面的价值，因为它们的原理是依赖于方向的强化和依赖于频率的强化之一。局部地需要与主方向和主频率相关的信息，与它们的局部变化一起，作为局部伽柏滤波的输入值。预先假定对这些的计算达到了足够的质量，并且所述计算不是本专利技术的组成部分。图像段(例如，32×32块)的伽柏滤波通常是这样进行的借助变换到频域(FFT＝快速傅立叶变换)的转换，然后逐点乘以由局部参数预先计算出来的块尺寸的伽柏滤波器掩码，然后反变换到时间和空间域(IFFT＝反快速傅立叶变换)。由于所发生的变换(FFT和IFFT)的伪像，各个得到的块的组合因此通常不够精确，如果这是在不与相邻块重叠的情况下进行的话。此类的窗口和重叠技术是现有技术部分。按照已知的方法，通过频域的转换得到了接受，因为二维卷积在时域内具有高度的复杂性，如果不采用通过频域的转换，则就必须完成时间域内的二维卷积。本专利技术的目的是实现图像的伽柏滤波，尤其是实现皮肤印痕图像的伽柏滤波，就程序实现而言其可能的复杂性最小，其中滤波要在很大程度上适用于局部。这个目的按照本专利技术是这样实现的，其中将图像分为多个块，对于各个块确定条形结构的主导方向，并且以这样一种方式进行滤波一次旋转一个块，直到主导方向与伽柏滤波器的主方向垂直；在主方向上进行一次滤波，而在与主方向垂直的方向上进行另一次滤波；并且将经过滤波的块再次旋转回来。本专利技术能够直接在时域中实现局部和自适应伽柏滤波，而不必通过频域进行转换，无需计算“真正”的二维卷积。这是借助将单独的二维伽柏滤波器分解为两个一维滤波器来实现的，这两个一维滤波器的向量积得出所求的伽柏滤波器，其中时间顺序可以自由选择。不过，仅仅在余弦振荡的波前和所叠加的高斯钟形曲线的正交主轴沿着轴向行进的情况下(根据本专利技术这是通过滤波器的固定排列伴同块的相应旋转来实现的)，才可能精确地分解为两个一维滤波器。与伽柏滤波在频域中的实现不同，时域中的伽柏滤波的特征在于具有下述优点-可以自由选择伽柏滤波器的大小(不必与所选择的块大小相当)，-采用给定的参数，一维滤波器的结构实现起来极其简单，-可以自由选择块的大小，不会造成所需时间和程序存储方面的任何额外的成本升高(可以采用非二次的块，不仅仅是并矢量，这优选用于FFT/IFFT)，-程序简单易行，-需要相对较小的程序存储，-能够很容易地实现合格的或所需的精度(在浮点算法和定点算法中都是如此)，-不需要对目标平台执行二维FFT/IFFT。滤波器对特定块的自适应优选是这样进行的逐块地，针对滤波之一，得出频率与垂直于主导方向的结构的频率相等的余弦振荡，并且由高斯钟形曲线对该余弦振荡进行调制，并且，如果可行，则逐块地，对于另一个滤波，高斯钟形曲线的宽度取决于块上结构的方向变化。如果块上结构的方向仅发生轻微变化，则可以并且优选采用较宽的高斯钟形曲线。于是可以将更大的区域包含在沿着这些结构的方向的滤波中，从而抑制了不属于这些结构的不规则事物，例如，小点。不过，在块发生显著的轮廓变化的情况下，即，急剧弯曲的结构，就需要沿着这个方向进行窄带滤波。还有可能将余弦振荡的方向上的高斯钟形曲线的宽度设置为取决于块上的频率变化，不过，如果块大小得到了适当匹配，则在指纹图像的情况下，通常不需要这样设置高斯钟形曲线的宽度。为了进一步减小执行和计算输入，可以在按照本专利技术的方法中对为旋转而定义的选定角度(在特定程序中实现)进行规定，然后将所定义的角度中与旋转本身所需要的最接近一致的一个角度用于滤波应用。这样，例如，可以允许±5°、±10°、±15°、…、±90°的旋转，在它们的复杂程度方面，使它们得到了最小化，尤其是对特定的角度(±45°，±90°)。依靠这一背景，可以实现在旋转的各个可选择角度的时间方面最优的程序。由于特定块的各个图像元素在旋转期间无论如何也要重新计算，因此按照本专利技术的方法还可以进行这样的简化在旋转期间，低通滤波是通过内插进行的，和/或在反向旋转期间同时进行二值化。在要进行滤波的块按照本专利技术的旋转期间，不可避免地会有更多的图像元素包含在滤波过程中，从而最终对比结果块大的块进行了滤波。这个块的大小取决于各种不同的影响，例如，包括滤波器的大小。为了对不与相邻块重叠的大小限定的块进行滤波，本专利技术在这方面具有优势的一种实施方式包括形成与相邻块重叠的更大的块，并且该更大的块具有至少根号二倍(double the root)的尺寸；并且，在旋转之后，在边长相当于该更大的块边长的至少根号二倍的正方形内对该更大的块进行滤波。在按照本专利技术方法的情况下，可以进行计算时间的进一步最小化在滤波期间，不予考虑位于阈值以下以及处于一维滤波器的边缘上的项(值)。将参照附图中所示的实施方式的实例对本专利技术加以进一步介绍，不过，本专利技术并不局限于这些实例。附附图说明图1示意性地表示按照本专利技术的方法的各个步骤。附图2表示在按照本专利技术的方法的不同阶段下要被滤波的块的大小。附图3表示指纹图像的一个实例。附图4表示伽柏滤波器的立体图。附图1表示指纹图像块1的伽柏滤波，其中，由于专利附图的规定所限，滤波本身的效果没有示出。块1是附图3所示的指纹的一部分，并且可以具有，例如，32×32个图像元素。不过，也可以采用其它大小和非方形、正方形的块。采用已知的伽柏滤波，在附图1中所示的例子中必须要进行二维滤波，其中伽柏滤波器是依照倾斜地且以曲线行进的条纹2排列的。不过，在按照本专利技术的方法中，这些条纹的平均方向是通过本身公知的算法确定的，根据这些算法，按照这样一种方式旋转块1条纹2的平均方向处于直角处。由滤波器4对这样旋转过的块3进行滤波，该滤波器4包括由高斯钟形曲线调制的余弦函数。该余弦函数的频率预先与条纹2的空间频率进行了匹配。高斯钟形曲线的宽度取决于条纹的曲率，并且如果可应用的话，取决于条纹空间频率的变化。在进一步的步骤中，由高斯钟形曲线5与之前的滤波过程成直角地对经过初始一维滤波的、经过旋转的块进行滤波。这一高斯钟形曲线的宽度取决于条纹2的平均曲率。随后，将经过二维滤波的块3’旋转回到其起始位置。附图2表示，为了对经过旋转的块进行滤波，必须在一开始时得到更大的块，以避免边缘处的伪像。由此起始点是具有32×32个图像元素的结果块。由于块的旋转，需要考虑更大数量的图像元素。这一变量是通过将边长的根号二倍作为对角线给出的。如果在旋转期间要进行双线性内插，则被滤波的大小为32×32的块必须通过具有进行缓冲的最小尺寸的、50×50的块的反向旋转来得出。如果将这一缓冲扩展到52×52个图像元素，则可以完全避免伪像。为了获得由15×15个图像元素的滤波器进行滤波的大小为52×52的有效区域，最终要产生大小为66×66的要进行缓冲的块。附图2表示在执行按照本发本文档来自技高网...

【技术保护点】

【技术特征摘要】
【国外来华专利技术】

【专利技术属性】
技术研发人员：S·肖尔茨，
申请(专利权)人：皇家飞利浦电子股份有限公司，
类型：发明
国别省市：