MTF Measurement Method For Infrared Dual-Mirror-Refracted Optical System
-
摘要: 为了压缩总长度，红外热像仪普遍采用双反射镜折布局（U型）及箱式壳体，而现有传函仪仅能测量非折转布局（一字型）的光学系统的调制传递函数（Modulation Transfer Function，MTF），导致传函仪无法接到U型光学系统的像面，因此无法测量其MTF。为了解决该难题，本文采用光瞳衔接的方式，设计和加工一个中继镜组把U型光学系统的像面引出，通过测量组合系统的MTF，再计算出U型光学系统的MTF，从而解决其MTF测量的难题。然后使用6个同一个型号的红外热像仪对测量方法进行验证，在排除影响测量的干扰因素以后，通过测量组合系统MTF，再计算出U型光学系统的MTF。同时，在测量过程中发现5#、6#红外热像仪的U型光学系统的MTF下降较多, 通过红外中心偏差测量仪的测量，发现第三、四镜片的中心偏差较大，故重新装调和控制，从而使U型光学系统的MTF有大幅提升。该方法为U型光学系统的MTF测量和性能提升提供了新的思路。Abstract: To compress their total length, infrared thermal imagers commonly have a double-reflector folded layout (U-shaped) and a box-shaped shell. Existing transmission instruments can only measure the modulation transfer function (MTF) of optical systems with a non-folded (straight) layout owing to the inability of the transmission instrument to receive the image plane of U-shaped optical systems. To solve this problem, this study used the pupil connection method and designed and processed a relay lens group to direct the image plane of the U-shaped optical system. The MTF of the U-shaped optical system was calculated by measuring the MTF of the combined system, thereby solving the MTF measurement problem. Six infrared thermal imagers of the same model were used to verify the measurement method. After eliminating the interference factors affecting the measurements, the MTF of the U-shaped optical system was calculated by measuring the MTF of the combined system. During measurement, the the U-shaped optical system MTF of infrared thermal imagers 5 and 6 decreased significantly. Using an infrared center-deviation measuring instrument, it was found that the center deviations of the third and fourth lenses were relatively large. Therefore, they were reinstalled and controlled, resulting in a significant improvement in the MTF of the U-shaped optical system. This method provides a new approach for MTF measurement and performance improvement in U-shaped optical systems.
-
Keywords:
- double reflective mirror folding /
- MTF measurement /
- pupil convergence
-
0. 引言
在红外热成像设备的设计和应用中，经常针对不同厂家或者不同型号的红外探测器定制非标的红外图像处理算法，以及标定不同的算法参数。由于算法模型或者算法参数经常变化，导致设计算法、仿真、在硬件上调试算法的时间很多，开发周期很长，开发效率低下等问题。而常规的红外图像处理算法，包括数据采集、时域滤波、非均匀性校正、空域滤波、锐化增强、坏点校正以及调光映射等。设计的初期，需要从算法原理设计、算法模型的建立、数学过程的计算，转化成可实现的编程语言，然后实现和显示。
随着小型化、低功耗、低成本的红外热成像设备应用的普及和民用。FPGA处理器非常适合于该种应用需求，FPGA是高速并行的处理器，适合于图像阵列的算法处理和加速，并且功耗低，可反复编程设计，所以常应用于红外热成像设备中。
1. 传统的红外图像仿真方法
基于FPGA的红外图像处理算法，仿真流程示意图如图 1所示。
从图 1中可知，利用FPGA处理器进行图像处理，需要两步的仿真结果。
第一步：根据算法需求，设计算法原理，在Matlab软件上编程实现算法过程，先用浮点数进行设计和计算，代码设计和编写完成后，导入测试图像，进行仿真和显示，在显示窗口中显示出来[1-4]。对显示出来的图像效果进行评估，看看图像效果有没有达到算法原理期望的结果，如果没有，调整算法参数和实现过程，反复迭代调整直至图像效果达到期望的结果。达到期望结果后，进行定点化设计，将算法实现过程中的浮点数，全部转化成定点化数，比如INT8、INT16、INT32等，并且转化完成后，保证算法实现的结果仍然达到期望值，不损失精度，不丢失细节。然后再更换测试图像，测试数据设计成依次递增或者递减的特征数据，再进行仿真。
第二步：文献[5]、[6]提出定点化之后的算法原理和实现过程，编写FPGA的硬件描述语言，设计FPGA的逻辑，将第一步过程中的算法原理实现出来，再进行ModelSim仿真，仿真测试的输入激励是第一步中同样依次递增或者递减的特征数据，得到仿真结果，与第一步中每个关键节点（关键步骤）的结果进行比对，保证每一步的结果都与Matlab仿真的结果相同，才能确保算法效果的达到。
传统的红外图像仿真方法，分成上述两个步骤，先要进行Matlab仿真，显示评判效果达到之后，再启动定点化，将定点化之后的算法过程转化成FPGA中可以实现的过程，并用FPGA的硬件描述语言实现相关的逻辑，再进行ModelSim仿真，由于ModelSim不能像Matlab工具一样方便导入一张图像（图片数据）进行仿真，仿真结果不能直观地可视化显示。所以在第一步、第二步过程中，输入激励的时候制作一些递增、递减的相同特征数据进行仿真。对每一步关键步骤和最后结果的仿真比对，保证数据正确，才能保证算法原理的实现和算法效果的达到。
2. 基于ModelSim可视化的红外图像仿真方法
从上述传统的红外图像仿真过程看，既要实现Matlab代码编写、仿真、还要做定点化设计。等FPGA代码编写完成后，在ModelSim仿真的时候还要进行每一步过程的计算结果比对，并且不能直接可视化的方式看到算法运行的结果。存在实现过程复杂、开发效率低下，开发周期长等问题。针对这些问题，本文提出了一种新的基于ModelSim可视化仿真的红外图像仿真方法。如图 2所示为本文的实现框图。
从图 2中可知，基于ModelSim可视化的红外图像仿真方法相比传统的红外图像仿真方法，仿真过程简单很多，开发效率也会提高很多。在算法原理设计完成后，直接编写FPGA的硬件描述语言实现算法计算过程，再进行ModelSim仿真，仿真的时候，不再用特定的递增数据或者递减的特征数据做输入激励，而是用之前原始采集到的红外图像，如果算法原理是针对16 bit的原始红外图像处理，可以是Y16数据，针对8 bit的调光之后红外图像处理可以是Y8数据，之前采集到的红外图像Y16或者Y8数据，存放在RAW文件中。而RAW文件中的图像数据是可以用ImageJ查看的，ImageJ工具可以直接显示16位或者8位的RAW图像。
对于红外图像算法所用到的数据文件——RAW文件，是不能直接导入ModelSim进行仿真的，需要经过binary工具将raw文件转换成txt文件，在ModelSim中通过系统调用＄readmemb/ readmemh函数读入TXT文件的方式读入到内存中，然后进行图像算法处理，算法处理完成，调用＄writememb/ ＄writememh函数将算法处理之后的图像写入txt文件中，通过txt2bin工具，将txt格式的数据转换成bin文件，再通过直接改后缀名，将.bin改成.raw文件。这样输入输出都是raw格式的数据文件，导入ImageJ工具中显示，以可视化的方式显示出了图像。
所以，这样红外图像算法处理的输入是图像数据，输出是处理之后的图像数据，可以直观地在ImageJ工具中显示出来对比，分析红外图像算法的处理过程是否正确。并且每一步的计算过程和关键步骤的结果也可以通过如图 3内存查看的方式、图 4内存读写过程仿真的方式查看数据在内存中存储、写入和读出是否异常。并通过如图 5红外图像算法处理仿真过程来查看数据计算过程正确与否。
![]()
图 4 红外图像数据内存读写仿真波形数据Figure 4. Infrared image data memory reading and writing simulation waveform data![]()
图 5 红外图像算法处理仿真波形数据Figure 5. Infrared image algorithm processing simulation waveform data这样数据输入，写入内存，然后从内存中读出，再进行红外图像算法处理，每一步数据的处理都是可以查询和追溯的，所以既方便又快捷。如图 6所示为本文进行红外图像滤波降噪算法处理过程中的输入源图和处理完成之后输出结果对比图Image显示，图 7为本文中基于ModelSim仿真过程中输入输出文件结构示意图。
![]()
图 6 红外图像滤波降噪输入输出对比图Figure 6. Infrared image filtering noise reduction input and output contrast diagram![]()
图 7 基于ModelSim仿真过程文件结构示意图Figure 7. Schematic diagram of file structure based on ModelSim在反复调整算法过程，以及算法参数的情况下，能够很快得到想要的算法结果，比如图 6中右边是降噪之后的效果图，降噪之后，图像噪声能消除了，但是图像变得模糊了，说明降噪算法参数设置不合适，通过调整滤波器的系数和阈值之后，降噪效果有明显的改善，如图 8所示。
![]()
图 8 红外图像滤波降噪后优化参数效果图Figure 8. Optimization parameter effect after infrared image filtering and noise reduction上述图 6与图 8中的红外图像滤波降噪算法为非局部平均滤波降噪算法，算法流程如图 9所示。这种降噪算法是一种利用图像自相似性的空域滤波算法，相比传统高斯平滑、双边滤波等方法，保护更多的图像细节。图像各点之间可以同时进行运算，非常适合于在FPGA并行处理器上运算。根据相似性计算权重的时候，如果相似性阈值选择太大，算法滤波出来的效果如图 6所示，过度平滑了，图像变模糊了，而选择合适的相似性阈值，得到滤波出来的效果如图 8所示，既能滤除噪声，又能保留图像细节。在该种红外图像仿真方法中能很快得到一组适合的算法参数。
基于该种ModelSim可视化的红外图像仿真方法，能够快速地修改算法参数和调整算法计算过程，运行之后得到图像算法的仿真结果，并且输入、算法计算、输出等环节里的每一步数据都是可以通过内存或者波形查看的方式追溯数据处理过程，保证数据向着期望的方向计算和处理。
3. 分析讨论
本文提出两种红外图像的仿真方法，传统的红外图像仿真方法，需要进行Matlab仿真，仿真结束后，再进行定点化设计，才能指导FPGA进行编程和实现，最后启动ModelSim的仿真，并且传统的ModelSim仿真并不能直接对图像数据的输入进行仿真并且可视化，而是制作相关的递增或者递减的特征数据作为输入，并且Matlab仿真同样以该特征数据作为输入，进行仿真，最后对2组仿真数据结果进行比对，比对结果一致代表该算法在FPGA中运行正常，能达到期望的算法效果。该过程复杂、繁琐、开发效率低下、开发周期长等问题不可忽视。所以本文在传统的红外图像仿真方法上进行了改进和优化，提出了一种新的基于ModelSim可视化的红外图像仿真方法与系统。该种方法在传统方法上省去了Matlab仿真的过程和环节，同时也省去了定点化的设计过程。在算法原理设计完成后，启动FPGA硬件描述语言的编写和仿真，只是在仿真的时候，输入图像可以用之前采集好的红外图像或者用PS工具生成一张图片，保存为.raw文件，在仿真输入的时候用binary工具将raw格式文件转换成txt文件，就可以启动仿真过程了，仿真结束后，生成的结果TXT，再通过txt2bin工具，将txt格式的数据转换成bin文件，再通过直接改后缀名的方式，将.bin改成.raw文件。用ImageJ工具打开source.raw和result.raw文件就可以直观地对比显示在图像窗口上了。并且对于算法参数、计算过程的调整，也是很快得到相应的输出结果。该过程相对传统的Matlab仿真和ModelSim联合仿真要方便很多，尤其是在算法调整和参数的调整的情况下，又要重复一遍Matlab仿真，再去ModelSim仿真的过程。但是该种方法对于红外图像算法原理的理解和FPGA逻辑实现图像算法的能力要求比较高。并且本文针对红外图像的仿真方法与传统的红外图像仿真方法，在适用范围上是一样的，对象都是红外图像的算法仿真。针对具体的算法模型，需要将算法分解成各计算步骤，然后针对各计算步骤进行数据仿真，所以对于其他文献[7]中不同的算法，比如弱小目标的背景的抑制，突出目标的增强算法等都是一样的仿真过程。
4. 结论
对于红外图像处理算法在FPGA的实现过程中，本文阐述了传统的仿真方法，先进行Matlab仿真，仿真效果达到的前提下，再进行定点化设计，指导FPGA进行逻辑设计，完成算法原理的实现，最后启动ModelSim仿真，而且仿真结果的正确评判标准是与Matlab仿真结果比对一致，不一致的话，要优化定点化过程，再进行仿真，重复直至结果一致。而本文的仿真方法基于ModelSim的可视化仿真，省去了上述繁琐的开发过程，直接将图像源文件转化后导入ModelSim仿真，算法处理完成后，转成图像结果文件进行比对，用可视化的方式评判图像经过算法处理之后是否达到算法原理期望达到的算法效果。同时如果算法未到达期望效果，通过调整参数、算法实现过程，以及追溯输入图像、写入、读出内存过程，判断图像算法计算过程中数据的正确与否。所以该方法与传统的红外图像仿真方法相比，评估算法效果以可视化显示输入、输出对比图差异性来达到。评估算法计算过程的正确性，同样可以追溯数据处理过程的每一步的正确性来达到。所以该仿真方法的仿真效果是与传统的方法一致的。但是开发过程缩减了和效率上大大提升了，开发周期短。尤其是在反复调参和优化过程中，显得更为明显和方便便捷。
-
![]()
图 2 红外热像仪光机机构示意图
Figure 2. Schematic diagram of opto mechanical structure of infrared thermal imager
![]()
图 13 红外热像仪和传函仪的安装基准对准
Figure 13. Align the installation benchmark of infrared thermal imager and transmitter
![]()
图 14 接口工装左、右对和存在夹角安装方式的MTF
Figure 14. MTF with multiple installation methods for interface tooling
![]()
图 16 U型红外光学系统的MTF理论值
Figure 16. Theoretical value of MTF for U-shaped infrared optical system
表 1 技术指标要求
Table 1 Technical indicators requirements
Relay lens group U-shaped infrared
optical systemWorking band/μm 3-5 3-5 F# or NA 0.2 3 Magnification 1× - Object height/mm 12 - 
下载: 导出CSV
表 2 中继镜组的MTFr理论值与实测值
Table 2 Theoretical and measured values of MTFr for relay lens group
Frequency/(lp/mm) MTF The measured value of MTF for Auxiliary lens The measured value of MTF for combination system Ⅰ Calculated value of MTF for relay lens group Theoretical value of MTF for relay lens group 0 1.0 1.0 1.0 1.0 5 0.96 0.89 0.92 0.93 10 0.93 0.82 0.86 0.88 15 0.88 0.72 0.80 0.82 20 0.85 0.66 0.73 0.77 25 0.82 0.57 0.66 0.70 30 0.79 0.52 0.60 0.67
下载: 导出CSV
表 3 6个组合系统的MTF测量值
Table 3 The measured value of MTFc for 6 Combination System
Order Frequency 5lp/mm 10lp/mm 15lp/mm 20lp/mm 25lp/mm 30lp/mm 1# 0.73 0.5 0.38 0.326 0.28 0.241 2# 0.72 0.48 0.35 0.28 0.25 0.22 3# 0.71 0.483 0.35 0.277 0.22 0.188 4# 0.71 0.46 0.35 0.26 0.22 0.18 5# 0.6 0.298 0.21 0.166 0.14 0.125 6# 0.65 0.395 0.25 0.155 0.1 0.07
下载: 导出CSV
表 4 红外热像仪U型红外光学系统的MTF测量值
Table 4 The measured value of MTF for U-shaped infrared optical system
Order Frequency/(lp/mm) 5 10 15 20 25 30 1# 0.777 0.581 0.475 0.453 0.424 0.402 2# 0.766 0.558 0.438 0.389 0.379 0.367 3# 0.755 0.562 0.438 0.385 0.333 0.313 4# 0.755 0.535 0.438 0.361 0.333 0.300 5# 0.638 0.347 0.263 0.231 0.212 0.208 6# 0.691 0.459 0.313 0.215 0.152 0.117
下载: 导出CSV
表 5 红外热像仪各镜组的中心偏差
Table 5 Center deviation of each lens for Infrared thermal imager
Lens Order Center deviation /μm 5# 6# First lens 1 5 3 2 8 5 Second lens 3 8 11 4 10 6 Third lens 5 33 22 6 56 38 Fourth lens 7 24 54 8 18 46 Fifth lens 9 9 7 10 6 9 Sixth lens 11 10 8 12 3 5
下载: 导出CSV
表 6 U型红外光学系统重新装调后的MTF测量值（5#、6#红外热像仪）
Table 6 The measured value of MTF for U-shaped infrared optical system After reinstallation and adjustment (5#、6# Infrared thermal imager)
Order Frequency/(lp/mm) 5 10 15 20 25 30 5# 0.761 0.575 0.442 0.386 0.342 0.315 6# 0.760 0.58 0.438 0.372 0.337 0.306
下载: 导出CSV
-
[1] 李林, 黄一帆, 王涌天. 现代光学设计方法[M]. 北京: 北京理工大学出版社, 2018: 152. LI lin, HUANG Yifan, WANG Yongtian. Modern Optical Design[M]. Beijing: Beijing Institute of Technology Press, 2018: 152.
[2] 张春仙, 李忠升, 张昦润. 红外热成像系统传递函数的测量研究[J]. 激光与红外, 2019, 49(4): 442-446. ZHANG Chunxian, LI Zhongsheng, ZHANG Haorun. Study on measurement of MTF in IR systems[J]. Laser & Infrared, 2019, 49(4): 442-446.
[3] 李航, 颜昌翔, 于平, 等. 红外成像系统的调制传递函数测试[J]. 光学精密工程, 2016, 24(4): 699-707. LI Hang, YAN Changxiang, YU Ping, et al. Measurement of modulation transfer function for IR imaging system[J]. Optics and Precision Engineering, 2016, 24(4): 699-707.
[4] 胡涛, 张冬冬. 高精度刃边法测量红外相机MTF的研究[J]. 红外, 2015, 36(7): 10-15. HU Tao, ZHANG Dongdong. Research on high-precision edge method for measuring MTF of infrared cameras[J]. Infrared, 2015, 36(7): 10-15.
[5] 黄海乐, 朱雷鸣, 彭宇, 等. 刃边法检测空间相机MTF的方法研究[J]. 光学与光电技术, 2013, 11(2): 71-73. HUANG Haile, ZHU Leiming, PENG Yu, et al. Research on edge method for detecting MTF of space cameras[J]. Optics and Optoelectronics Technology, 2013, 11(2): 71-73.
[6] 于谦. 光学检查镜头的MTF测量方法研究[D]. 杭州: 浙江大学, 2010: 15-18. YU Qian. Research on MTF Testing Methods of Optical Inspection Lenses[D]. Hangzhou: Zhejiang University, 2010: 15-18.
计量
- 文章访问数: 126
- HTML全文浏览量: 53
- PDF下载量: 38
