MTF和空间频率计算
在使用分辨率测试卡检测系统高分辨率时,常常设计到MTF的计算问题。但是大多数时候这种计算是由测试软件完成,但是了解这个计算过程对我们熟悉测试卡使用还是有好处的。
可以通过首先找到条的空间频率,然后基于傅立叶系数导出新的系数集,从楔的横截面更精确地计算MTF。空间频率以每像素的周期(C/P)为单位计算,可以使用简单方程将其转换为每个图像高度的线宽(对于任何指定的高度)
f(LW/PH)=2f(C/P)x相高 (5)
1.在y的每个值找到楔子边界(注意忽略干扰图案,如tic标记或数字)。
2.对于在分辨率测试卡楔上的每个扫描线处测量的信号振幅I,在边界内找到平均值(I)以用作找到I -平均值(I)从正到负或从负到正的交叉位置的阈值。这些交叉用于确定检测到的条的数量在低频值以下,对应于f = falias。
3.对于其中检测到条的数量等于低频值(f <falias)的每个扫描线(y位置),正到负和负到正交叉间隔的平均值是条的周期像素,等于以周期/像素(C/P)为单位的空间频率f的倒数。
4.对于双曲楔形,空间频率随距离线性变化。在每个扫描线使用线性回归拟合f <falias以找到用于所有计算和绘图中的空间频率f(y)的一阶方程。 [注:我们使用f <0.9 falias为回归拟合获得更好的结果。]
因为《ISO 16505 MTF和空间频率计算》等式(1)仅对正弦波严格正确,并且因为Imax-Imin受到来自分辨率测试卡上相邻条的噪声、采样相位和干扰的产生误差,所以用等式(1)中的Imax-Imin代替I内的傅里叶系数 其中边界X1和X2表示离散数量的检测周期。
使用√(ccos2 + csin2)而不是Imax-Imin是更准确的,因为它表示从积分(一种同步检测)导出的正弦系数,减少了除f(y)以外频率处的信号噪声影响。