MTF的倾斜边缘法计算方法.docx

MTF的倾斜边缘法计算方法 简介 光学系统性能的衡量方法有很多，常见的有点扩散函数法、瑞利判断法、点 列图法、光学传递函数 (MTF)法等，其中 MTF法在光学系统和镜头加工制造中 使用 最为广泛。 MTF曲线真实的反映了成像系统将物方信息传递到像方的能力。 MTF曲线的横坐标一般是 [1][11] ，纵坐标是反映对比 cycle/mm或者 linepair/mm 度传递特性的像 /物方调制度的比值。 MTF的计算方法有很多，比如使用条形目标的方法、随机目标方法、点扩散函数法、带宽受限激光散斑法 [2]、倾斜边缘法 (SlantedEdgeMethod)等，其中 点扩散函数法和倾斜边缘法是工业界普遍使用的方法。点扩散函数法的的缺点在 于需要使用小于光学系统或镜头分辨极限的点光源，这样使得成像器件很难获得 足够强的信号。 图 1 倾斜的边缘 (SlantedEdge) 倾斜边缘法的算法流程如图 2：首先获取倾斜边缘的边缘扩散函数 (ESF)，然后求导得到对应的线扩散函数 (LSF)，最后傅里叶变换得到 MTF。 图 2 算法流程图 公式推导 [3] 倾斜边缘的响应函数可以由一个冲激函数表示： 当边缘响应函数由完善的 (没有像差 )的光学系统成像时，系统的成像质量 不会被劣化。因此边缘函数被线性不变的光学系统成像时，系统的输出 O(x)等于线传递函数 LSF与系统的响应函数 S(x)的卷积： 当 x-α 0时，阶跃函数 S(x)=0，其他情况下 S(x)=1，所以 ESF(x)可以表示 为： 因此， ESF(x)的导数可以写为： 所以可以将 MTF写作 LSF的如下函数： 通常， MTF会对零频率幅值归一化，同时由卷积定义及傅里叶变换理论可以推导得出级联系统的 MTF： 影响因素 倾斜角度 需要注意的是，边缘倾斜的角度对计算的结果有一定的影响。 一般来说，倾斜角需要控制在 2°-10之°间，当角度超出这个范围时，计算的结果会出现比较大的偏差 [4][5][6]。 图 3 倾斜角度对 MTF计算结果的影响 图 4 倾斜角度对 MTF50计算结果的影响 当倾斜角度比较大时，可以取样相邻的多行 /列的数据进行计算，这种方法 其实就是过采样 (oversampling)校正。 ISO12233建议 4倍的过采样 (获取邻近的 4 行/列线条的 ESF)[4]。 对比度 在没有噪声，并且边缘对比度足够高 (20%)的情况下，对比度对计算几乎没有影响 [5]。当边缘对比度低于 20%时，计算结果将出现大的偏差 [6]。 数据限幅 (dataclippig) 如图 5 所示，当相机的曝光时间或者增益偏大时，会导致边缘扩散函数的上升沿呈 直线状 或者顶部变得十分平坦。这种情况下，计算的 MTF都会出现极大的偏差。 使用算法、软件对图像去噪或者增强对比度都可能使 数据限幅 ，因此必须谨慎的使用图像增强处理。 图 5 数据限幅 (或图像增强过度 )时的边缘扩散函数 ImatestVSsfrmat3VSImagej(SE_MTF) Imatestmaster是一款用于检测数字相机、摄像机、工业相机拍摄的图像质 量及交流检测结果的商业软件。 Imatest包含多个模块，如 ① SFRplus，可以与 SFRplus测试卡配合使用，检测 MTF，横向色差，畸变，阶调响应和颜色还原。 ② SFR，配合测试卡，可检测相机和镜头的锐度，噪音，色差。 ③MTFCompare，使用 SFR的结果比较不同相机、镜头和影像系统的 MTF性能。 ④ Stepchart，配合灰阶测试卡检测相机的阶调响应，噪音，动态范围，还可用于检测曝光的准确性，镜头眩光。 ⑤ Distortion，检测镜头的畸变，使用可打印的正方形和矩形网格计算矫正系数。 ⑥TestCharts，可以生成图像文件用打于打印测试图。生成的测试卡包括 SFR斜边图像，星形测试图，各种频率和反差的图样。可选项 包括：反差，高光色彩，正弦线和条纹模式。可输出位图或 可绽放矢量图 (SVG)[7]。 sfrmat3是 PeterD.Burns编写的用于评估数字相机与扫描仪的成像空间频率响应 (SpatialFrequencyResponse)的Matlab程序。 sfrmat3遵循 ISO12233 的标准，使用倾斜边缘计算图像的空间频率响应 [8]。 ImageJ是由 NationalInstitutesofHealth开发的，基于 java的公共的图像 处理软件。 ImageJ被设计为一个开放的架构，可以通过 java插件与可读的宏来扩展功能。因为本身是插 件架构、并内建了开发环境，因此 ImageJ已经成为一个十分受欢迎的图像处理教学平台 [9]。ImageJ有一个使用倾斜边缘方法计算 M