表面纹理分析中的样条滤波器

        样条滤波器之所以得名，是因为通过滤波获得的轮廓是样条曲线。样条函数涵盖了不同的函数，最著名的是三阶多项式样条，也称为三次样条。通过三阶多项式中把分段过渡组合且足够平滑。对于这样的样条曲线，例如，作为均匀各向弹簧的弹性线。

        三次样条曲线是平滑连接一系列点的几何曲线。样条曲线用于在CAD程序中绘制曲线形状，也用于插值数据点。

        应用三次样条滤波器后，样条曲线不通过数据点，而是通过悬挂点。这些点通过弹簧连接到数据点，弹簧的特征是弹簧常数β，称为张力参数。设置此弹簧常数（设置为介于0和1之间的值）可以控制过滤器的行为。

        基于三次样条的样条滤波器的传递函数比根据ISO 11562的高斯滤波器更陡峭，这意味着样条滤波器，可以实现比高斯滤波器更尖锐的分离然而，两个滤波器的传递函数之间的差异非常小，以至于应用样条滤波器代替高斯滤波器不会导致任何特性值的显著变化。

图1：样条滤波器（实线）和高斯滤波器的传递函数

        样条滤波器是为了避免高斯滤波器的缺点而开发的与后者相比，样条滤波器是一个纯数字滤波器。这意味着不再需要像高斯滤波器所要求的那样对加权函数进行采样。因此，样条滤波器的实现涉及更少的不确定性，相比于高斯滤波器使得样条滤波结果适用于不同的实现更容易达成一致的结果。

        与高斯滤波器不同，样条滤波器由滤波器方程定义，而不是由加权函数。这些滤波方程是简单结构的线性矩阵方程。他们的算法可以用结构最简单的程序实现，运行速度比高斯滤波器已知的所有算法快得多。，

        样条滤波器有两种不同的变体：非周期样条滤波器应用于开放轮廓，如平面。以及应用于闭合轮廓的周期样条滤波器例如，当测量圆柱形表面时。非周期性样条滤波器的设计避免了高斯滤波器的末端效应。从图2中的示例可以看出这一点。为了便于比较图中还显示了使用高斯滤波器的滤波。两个滤波器的截止波长是相等的。

        如图1所示，样条滤波器不会出现类似于高斯滤波器的末端效应。这是一个很大的优势，因为使用样条滤波器，整个测量长度的信息可以用于进一步评估，而高斯过滤后，必须省略剖面两端的，以避免末端效应带来的数据篡改，例如粗糙度的计算参数等。与高斯滤波器相比样条滤波器将测量长度减少约30%。这个可以在各种情况下非常重要，特别是对于小零件，测量时间可以显著减少。

 图2：用样条曲线（实体）计算的平均线比较线）和高斯滤波器（虚线）用于粗糙度轮廓。两个滤波器的截止波长相等。

        基于三次样条的样条滤波器在滤波过程中保持三次多项式不变-而高斯滤波器只保留线性函数不变。样条滤波器相比于高斯滤波器可以更容易获取轮廓表面的大的形状成份。这一特性可以在图3中被看出。这里，为了便于比较，还显示了高斯滤波器的滤波结果（两个滤波器的截止波长相等）。很明显，与高斯滤波器相比样条滤波器更适合作为形状滤波器被使用。

 图3：用样条曲线（实线） 和高斯滤波器（虚线）相比获得更大的形状成份。两个滤波器的截止波长相等。

        尽管闭合轮廓不会出现末端效应问题，但问题可能尽管如此，当过滤闭合轮廓时。如果使用非周期性过滤器轮廓起点和轮廓终点连接处出现不连续并且被认为是扭结。这在图中的右手侧可以清楚地看到图4与此相反，对于周期性滤波器，连接是平滑的。因此需要两种不同类型的过滤器，一种用于打开，另一种用于关闭，充分证明。

 图4：通过周期（粗实线）和非周期（粗虚线）样条曲线滤器平均拟合圆显示为虚线。

        请注意，在图4的示例中，一个非常小的截止波数（大截止波长）的滤波器被选择用来提取形状成份，但它不能精确的表达圆（与显示的平均拟合圆比较）。

        在圆度测量中，较大的截止波数（较小的截止波长）将被选择，以便基本上抑制粗糙度，保留形状和波纹度误差。

下面将再次总结样条滤波器的优点：

1.纯数字滤波器，即与高斯滤波器相比不需要加权函数；因此，不会发生偏置效应；

2.快速算法（是高斯滤波器最佳算法的两倍）；它节省内存空间，并能容忍舍入错误，这意味着适用于大数据量的有效处理；

3.比高斯滤波器更陡的传递函数，即更好的判别作用；

4.无末端效应的影响，即可以使用测量的全部信息和/或测量长度可减少约30%；

5.滤波前不需要调整轮廓；

6.适合用作为形状滤波器；

7.可对所有边界条件进行调整，即可用开放和封闭轮廓。

        可以在没有任何问题的情况下扩展样条滤波器，以便可以作为形状滤波器使用。它们也可以被考虑作为3D曲线或非等距数据的滤波器。这意味着样条曲线过滤器也可以用于坐标测量机。

2 强健的样条滤波器

        滤波器的强健性是指在测量时，经常有的特性必须在受干扰的环境中确定数据。异常值必须考虑在其中。与非强健线性滤波器一样，对于异常值数据的平滑是敏感的，必须使用非线性滤波器。为此，各种方法被开发出来。基于已知标准方法的Lp范数强健样条滤波器，已知p=1…1,5。

图5：稳健样条滤波器（实线）和符合DIN EN ISO 11562 高斯滤波器（虚线）在多孔表面的过滤效果，两个滤波器的截至长度都相等（lamda c = 0.8）

        特别开发了强健的样条滤波器，以取代根据ISO 13545-1基于高斯滤波器的样条滤波器，与后者相比，它类似于普通样条滤波器，是一种纯数字滤波器。

        与样条滤波器类似，稳健样条滤波器有两种变体：非周期和周期强健样条滤波器。图5显示了将非周期强健样条滤波器应用于多孔表面的轮廓。正如可能的那样可见，孔隙被很好地抑制。

图6：用强健样条滤波器计算的轮廓线（实线）和符合DIN EN ISO 13565-1的特殊滤波器（虚线）过滤具有异常值的轮廓。两个滤波器截止波长相等（lamda c=0.8 mm）。

        示例中非常好的孔隙抑制的原因是该滤波器是强健的，如果滤波器能够抵抗异常值，则它们是强健的。正如如图6所示，根据DIN EN ISO 13565-1特殊滤波器，尽管负向异常值（指向下方）不会干扰过滤结果，但是正向值（指向上）会导致不希望得到的轮廓线。强健样条滤波器就不会出现这一状况。

图7：用强健滤波器计算的轮廓线（实线）和符合DIN EN ISO 11562的高斯滤波器（虚线）的轮廓的线。两个滤波器截止波长相等（lamda c＝8mm）。

        另一个典型强健样条滤波器的示例如图7，这里，滤波器用于获得在汽车工业中使用的钣金表面轮廓。这些钣金通常显示较大形状偏差。图7中轮廓的总测量长度为100 mm。然而，仅显示了45 mm和60 mm之间的部分，因为其它的细节将不可见。为了消除轮廓的形状偏差，选用了截止值为8mm的滤波器。如图7所示。强健样条滤波器给出了期望的结果，即接近曲面的平均线而根据ISO 11562的高斯滤波器给出了平均值线错误地向下移动。

3 结论

        已经表明，样条滤波器和强健样条滤波器具有显著的优势-与标准化高斯滤波器相比。样条滤波器是纯数字的滤波器，即与高斯滤波器相反，没有必要对加权函数进行离散化处理-。它们可以通过快速（两倍的高斯滤波器算法速度）节省内存算法得出结果，即它非常适合于大数据量的数据处理。样条滤波器的传递函数比高斯滤波的传递函数更陡峭，即更好的鉴别效果。与高斯滤波器相比没有末端效应。即，可以使用测量的全部信息，从而测量长度可减少约30%。当使用样条曲线滤波器时，任何过滤前的轮廓调整是不必要的。因为样条曲线可以很容易地跟踪任何形状偏差，样条滤波器非常适合应用于形状的过滤。调整样条滤波器可以适用于任何边界条件，例如可用于开放式和封闭式外形。因此样条滤波器可用于直线度测量和圆度测量没有任何困难。最后，样条滤波器也可以设计为不是等间隔数据或从数据点到数据点精度不相等的数据。样条滤波器易于扩展到区域过滤。因此，它们非常适合于表面纹理分析。对于这个应用领域，已经开发了快速算法。甚至可以过滤3D曲线。这使得样条滤波器适合在扫描模式坐标测量机中的应用。此类算法也已经实施。由于其优点，样条滤波器可用于需要平均线滤波器的任何情况。强烈建议以在标准化高斯滤波器不适用的所有情况下使用它们。本文给出了这些情况的一些示例。因为细节ISO文件中给出了样条滤波器和强健样条滤波器的可以预期，这两个滤波器将用于大多数表面纹理测量仪器的下一代分析软件、粗糙度测量、圆度测量和形状测量，以及在坐标测量机的预处理软件中。