Shack-Hartman Detector Real-time Wavefront Processor Based on FPGA
-
摘要: 夏克-哈特曼波前传感器是目前自适应光学系统中应用最为广泛的实时波前探测器。本文针对具有高分辨、高帧速、大规模子孔径数的夏克-哈特曼传感器，根据其波前处理计算量和实时性的要求，提出了一种基于现场可编程门阵列（FPGA，field-programmable gate array）的实时波前处理机结构及波前斜率计算方法。该方法利用核心处理模块重复利用的方式完成子孔径内光斑质心的计算，并通过USB3.0将处理后的质心数据实时传输给PC机。处理机以一片XILINX公司Kintex7-XC7K325T的FPGA作为处理芯片进行了设计，结果表明：该算法可对560帧/s的1020×1020图像（580 MB/s数据量），56×56子孔径哈特曼传感器，进行低延时实时光斑质心计算，提高了系统的波前处理速度和整个自适应光学系统的控制速度。Abstract: The Shack-Hartman wavefront sensor is the most widely used real-time wavefront detector in adaptive optics systems. In this study, a Shack-Hartmann sensor with high resolution, high frame rate, and a large-scale sub-aperture number is proposed. Based on the requirements of wavefront processing calculations and real-time performance, a field-programmable gate array (FPGA) is also proposed. The real-time wavefront processor structure and wavefront slope calculation method are investigated. The system employed the core processing module to reuse the method to calculate the centroid of the spot in the sub-aperture and transmitted the processed centroid data to the PC in real time through USB 3.0. The processor was designed with a XILINX Kintex7-325T FPGA processing chip. The results demonstrate that the algorithm can perform low-latency, real-time operations on 1020×1020 images and 56×56 sub-aperture Hartmann sensors at 560 frames per second. The spot centroid calculation increased the wavefront processing speed of the system and the control speed of the entire adaptive optics system.
-
Keywords:
- Shack-Hartmann sensor /
- real-time wavefront detector /
- FPGA
-
0. 引言
现代导弹的制导模式多种多样，具有单模、双模和多模3种制导模式，每种制导模式主要使用雷达、毫米波、主动/半主动激光、电视、红外等制导技术，各种制导技术根据自身的技术特点具有特定的使用范围。红外制导技术具有精度高、隐蔽性好、抗干扰能力强、能昼夜作战，在精确制导武器中备受青睐[1]。红外导引头分为热点式导引头和红外成像导引头[2]。红外成像导引头是红外制导导弹的关键部件之一，而红外导引头光学系统是红外成像导引头的“眼睛”[3]。红外成像导引头技术从单像元、线阵发展到面阵成像，其空间分辨率和抗干扰能力得到不断提高[4-5]。
制冷型红外探测器由于具有较高的灵敏度和较短的热响应时间，在早期的红外成像导引头中得到广泛的应用。随着非制冷型红外探测器制造技术的快速发展，热响应时间和灵敏度性能的不断提升，非制冷型红外探测器在红外成像导引头中逐渐得到应用。非制冷型红外成像导引头由于摒弃了体积大、重量重且价格昂贵的制冷机，与制冷型红外成像导引头相比具有体积小、重量轻、价格低的优点，可大大地减小导引头的尺寸和重量，实现了导引头的小型化，并且已成为红外成像导引头的重要成员之一。经过多年的发展，非制冷型红外成像导引头已广泛用于反坦克导弹、精确攻击导弹、小直径炸弹、反舰导弹等[6]。
基于以上分析，本文设计了一款适用于152mm中口径弹径的非制冷型、红外导引头的光学系统，并以仿真及试验结果说明此款红外成像导引头的性能。
1. 红外导引头系统组成
红外成像导引头主要由位标器及电子舱所组成，位标器主要包含红外成像组及伺服系统；电子舱主要包含伺服控制电路及图像处理电路。根据红外成像组与弹体耦合方式的不同，位标器可分为捷联式、万向支架式及陀螺式。捷联式位标器的红外成像组与弹体直接固定，不具备视场扫描能力；万向支架式位标器的红外成像组固定在万向支架上，可随万向支架转动形成扫描视场；陀螺式位标器的红外成像组全部或部分随陀螺转动形成扫描视场。捷联式红外成像导引头不具备视场扫描能力，只能跟踪固定大型目标或者运动速度较慢的目标，系统结构简单、制导精度低；万向支架式红外成像导引头具备视场扫描能力，可跟踪中速运动的目标，系统结构较复杂、制导精度较高；陀螺式红外成像导引头也具备视场扫描能力，可跟踪速度高的运动目标，系统结构更复杂、制导精度高。
成都鼎屹信息技术有限公司已经开发出采用384×288非制冷型长波红外机芯的万向支架式红外成像导引头，此款导引头在保证探测距离的条件下，瞬时视场比较小，需要辅助一定的框架角以提高搜索、跟踪的视场范围。红外成像导引头为了得到更远的探测距离、更大的瞬时视场，需要采用更长的光学焦距、更高分辨率的红外机芯，例如常用的640×512非制冷型长波红外机芯。
本文所设计的用于152mm中口径弹径红外导引头光学系统要求跟踪移动的中型坦克目标，其速度中等，可选用万向支架式的结构型式，此结构型式复杂度中等，精度也能满足要求，同时可降低制造的成本。
2. 红外导引头光学系统设计
2.1 设计指标
为了提高红外成像导引头的搜索、跟踪距离及瞬时视场，本文采用640×512非制冷型长波红外机芯；为了提高红外成像导引头的搜索、跟踪能力，采取了万向支架式的结构型式。红外成像导引头着重于搜索、跟踪远距离的目标，此目标所成图像为一个具有一定尺寸、一定对比度的热点图像。通过增大红外导引头光学系统的有效口径（减小光学系统的F/#），可增加接收的目标辐射能量，增大目标成像对比度及探测概率。640×512面阵、非制冷型红外导引头光学系统设计参数如表 1所示。
表 1 红外导引头光学系统设计参数Table 1. Design parameters of optical system for infrared seekerFocal length 50mm Field 12.4°×9.94° F/# 0.9 Wavelength 8-12μm Detector type UFPA 640×512, 17μm Frame angle ≥-30°~+30°（Pitch）
≥-20°~+20°（Level）Maximum blind area (target size 2.3m×2.3m) ≤50m Stable tracking range ≥2.5km 2.2 红外导引头光学系统性能优化分析
红外导引头光学系统由整流罩及红外成像组所组成。整流罩的面型大小及外形尺寸受红外制导导弹的气动特性及框架角的大小所决定。整流罩材料选用透过率高、工艺性比较好的热压ZnS材料。红外制导导弹在飞行过程中需要搜索、跟踪运动的目标，红外成像组在俯仰及偏航方向具有一定的转动角量。为了保证转动时光学成像特性不变，则整流罩需要做成同心、等厚的球冠，并且使红外成像组的转动中心与整流罩的球心相重合。本文所设计的整流罩内、外半径分别为72mm、76mm，厚度为4mm，外径为ϕ126mm。
红外制导导弹在发射前随季节及地理位置的变化环境温度变化较大，发射后飞行段弹体的温度也会发生比较大的变化，工作温度的变化使红外导引头光学系统的焦面发生一定量的位移并且不易人为调节补偿，造成成像模糊、不利于搜索、跟踪，因此红外导引头光学系统需要进行光学被动式消热差设计。
光学被动式消热差光学设计的基本原理在孙爱平[7]等相关论文中已有详细论述。红外成像组采取三片式设计型式，配合整流罩进行像差优化。在设计过程中加入衍射面，衍射光学元件其特殊的光学特性可单独作为一种特殊的光学材料参与像差优化。红外导引头光学系统选用线膨胀系数较小的Ge材料（α＝5.8×10－6/K）、线膨胀系数较大的IG6材料（α＝21.2×10－6/K）、整流罩热压ZnS材料（α＝6.6×10－6/K）及衍射元件配合镜筒Al材料（α＝23.6×10－6/K）进行光学被动式消热差设计。根据光学被动式消热差的基本原理，需要3种以上的红外材料进行消热差设计，本文采用3种光学材料及一种特殊的光学材料——衍射光学元件，足以满足光学被动式消热差设计的要求。红外导引头光学系统为大孔径光学系统（F/#＝0.9），相应的孔径像差比较难校正，它需要校正轴向球差、垂轴球差、轴向色差、垂轴色差、慧差、场曲、畸变7种像差及对应的高级像差。由于红外材料比较昂贵，在设计时需要遵循使用最少透镜的原则。本文4个透镜所提供的变量不足以校正7种像差及对应的高级像差，因此在设计时加入非球面以增加变量个数，实现比较好的像差校正效果。
红外导引头光学系统的衍射元件为基诺衍射光学元件，此衍射面理论上可放置在除球罩外的任何表面上，经过光学的优化设计结果比对，衍射面放置在第一透镜的第二表面上，像质效果最优。
红外导引头光学系统衍射面的基底面型为偶次非球面，对应衍射效率计算公式如下：
$${\eta _m} = {[\frac{{\sin [{\rm{ \mathsf{ π} }}(m - \frac{{n - 1}}{\lambda }d)]}}{{{\rm{ \mathsf{ π} }}(m - \frac{{n - 1}}{\lambda }d)}}]^2}_{}$$ (1) 红外导引头光学系统衍射面衍射级次为第一级（m＝1），中心波长为9.6μm，基底材料为IRG206（对应9.6μm波长的折射率为2.77908），对应刻蚀深度计算如下：
$$d = \frac{\lambda }{{n - 1}} = \frac{{9.6}}{{2.77908 - 1}}\,{\rm{ \mathsf{ μ} m = 5}}{\rm{.39605}}\;{\rm{ \mathsf{ μ} m}}$$ (2) 则中心波长处的衍射效率为100%。
基诺衍射面可由金刚石车床车削加工，基底可为平面、球面及非球面，其最小尺寸受刀具半径的限制。此衍射面的最小环带半径Δr间隔为1.055953mm，能满足加工厂家最小刀具的要求。基诺衍射面加工的理想面型如图 1所示。衍射面可使用轮廓仪进行检测，由峰-谷值和均方根值来判断衍射面的加工质量。基诺衍射面检验结果如图 2所示。
红外导引头光学系统的光学布局型式如图 3所示。使用ZEMAX软件进行像质优化，在常温（20℃）、低温（－40℃）、高温（60℃）特定工作条件下，中心频率点30lp/mm处的MTF曲线对比度大部分在0.5以上；通过判读点列图及能量包围曲线，约有85%的能量集中在一个像素内；通过场曲及畸变曲线可知畸变值均小于0.4%。通过以上分析，红外导引头光学系统消热差效果好，能够满足设计要求。图 4、图 5、图 6分别为常温（20℃）、低温（－40℃）、高温（60℃）条件下的MTF曲线、点列图、能量包围曲线、场曲与畸变曲线。
![]()
图 4 在20℃工作条件下红外导引头光学系统的MTF曲线、点列图、能量包围曲线、场曲与畸变曲线Figure 4. MTF curves, point diagram, energy enclosing curves, field curve and distortion curve of infrared seeker optical system at 20℃![]()
图 5 在－40℃工作条件下红外导引头光学系统的MTF曲线、点列图、能量包围曲线、场曲与畸变曲线Figure 5. MTF curves, point diagram, energy enclosing curves, field curve and distortion curve of infrared seeker optical system at -40℃![]()
图 6 在60℃工作条件下红外导引头光学系统的MTF曲线、点列图、能量包围曲线、场曲与畸变曲线Figure 6. MTF curves, point diagram, energy enclosing curves, field curve and distortion curve of infrared seeker optical system at 60℃红外导引头光学系统的焦深为±2λ(F/#)＝±2×9.6μm×0.92≈±17.7μm，所设计的光学系统工作在－40℃~60℃范围内时，离焦量在一个焦深以内，如表 2所示，则此光学被动消热差系统在整个工作温度内均能清晰成像。
表 2 红外导引头光学系统离焦量表Table 2. Infrared seeker optical system defocus gaugeTemperature Defocusing Depth of focus Remarks -40℃ 3.7μm 17.7μm Minimum temperature -20℃ 1.9μm - 0℃ 0.5μm - 20℃ 0μm Nominal temperature 40℃ -1.7μm - 60℃ -2.6μm Maximum temperature 2.3 框架角分析
随着红外制导导弹及目标的运动轨迹、运动速度的不同，在导弹飞行过程中出现目标移出瞬时视场的现象，造成目标丢失，此时需要红外导引头光学系统随目标移动方向转动形成一定的扫描视场，使目标始终处于红外导引头的瞬时视场内。红外导引头光学系统转动的角量对应着框架角值。
红外导引头光学系统转动角量为0°时，此时上边缘光线与球罩的交点到光轴的距离为h1，对应圆心角为α；当转动角度为β时，上边缘光线与球罩最外边缘相重合，并且与光轴距离为h2，此时对应圆心角为β；r为球罩外圆半径，如图 7所示。则框架角β计算如下：
$$\beta = \gamma - \alpha = \arcsin (\frac{{{h_2}}}{r}) - \arcsin (\frac{{{h_1}}}{r})$$ (3) 本文所设计红外导引头光学系统对应的h1＝29.5mm、h2＝63mm、r＝76mm，代入上式可得β＝33.1°，如图 8所示。由于红外导引头光学系统为球对称系统，故其框架角能够满足俯仰±30°、偏航±20°的要求。
![]()
图 8 导引头框架角在0°、33.1°工作时光学位置示意图Figure 8. Optical schematic diagram when seeker frame angle is 0° and 33.1°2.4 盲区分析
红外制导导弹是依据目标所成图像的轮廓特征实现搜索、跟踪。当目标充满成像机芯任意一个方向时，形成跟踪盲区，后续阶段依据导弹的惯性制导攻击目标。盲区距离越小，惯性制导段越短，攻击成功率就越高。
成像机芯分辨率为640×512，像元间距为17μm，对应的靶面尺寸为10.88mm×8.704mm，2.3m×2.3m的中型坦克目标充满水平、竖直向的条件下对应的盲区距离计算如下：
水平向盲区距离计算：
$$\frac{{10.88\,{\rm{mm}}}}{{50\,{\rm{mm}}}} = \frac{{2.3\,{\rm{m}}}}{{{L_1}}} \Rightarrow {L_1} = 10.6\,{\rm{m}}$$ 竖直向盲区距离计算：
$$\frac{{8.704\,{\rm{mm}}}}{{50\,{\rm{mm}}}} = \frac{{2.3\,{\rm{m}}}}{{{L_2}}} \Rightarrow {L_2} = 13.2\,{\rm{m}}$$ 综上所述，2.3m×2.3m的中型坦克目标充满成像机芯的短边时，盲区距离最远（13.2m），此距离值小于50m的要求。
2.5 公差分析
红外导引头光学系统按表 3分配零件公差，焦面位移作为补偿，以平均概率分布方式分配实际装配及加工时的公差值，并采用蒙特卡罗分析方法模拟加工装配后的虚拟镜头的MTF变化，依此判断实际镜头的成像效果。表 4的蒙特卡罗分析结果表明90%的镜头在奈奎斯特频率处具有不小于0.239的MTF值，焦面补偿在±0.5mm以内。此镜头加工及装配工艺比较成熟，且整个镜头的成像质量较好，即公差分配合理。
表 3 零件公差表Table 3. Part tolerance tableParameter Tolerance Parameter Tolerance Sphere error ±3 aperture Surface tilt ±50" Surface irregularity ±0.7 aperture Air thickness ±0.05 mm Aspheric error ±7×10-5 mm Element tilt 0.075mm Thickness ±0.03 mm Element eccentricity 3.6′ Focal plane
displacement compensation±0.5 mm - - 表 4 蒙特卡罗虚拟镜头分析结果Table 4. Analysis results of Monte Carlo virtual lensField Average 1 2 3 4 5 Nominal MTF 0.563 0.598 0.599 0.497 0.599 0.497 Optimum MTF 0.461 0.563 0.544 0.421 0.480 0.394 Worst MTF 0.217 0.224 0.222 0.190 0.251 0.198 Average MTF 0.343 0.358 0.355 0.314 0.366 0.300 Standard deviation 0.078 0.096 0.088 0.069 0.073 0.067 Compensation parameter statistics（Change of back focal plane） Maximum 0.074 Minimum -0.101 Nominal 0.009 Standard deviation 0.043 MTF analysis of Monte Carlo virtual lens（Nyquist frequency） MTF of Monte Carlo virtual lens（90%）≥0.239 MTF of Monte Carlo virtual lens（80%）≥0.272 MTF of Monte Carlo virtual lens（50%）≥0.345 3. 红外成像导引头实际效果
图 9是调试时截取观察到的高架桥上的私家车辆的图片，可清晰分辨出小尺寸的私家车辆目标；图 10是红外成像导引头外场打靶试验视频截图。通过以上试验可得出此红外成像导引头成像质量好，跟踪捕捉精度高，能够满足使用要求。
4. 结论
本文设计了一款适用于152mm弹径的非制冷型、红外导引头成像系统。详细介绍了红外导引头各种设计方案的优劣，提出使用万向支架式的设计方案。根据此方案对红外导引头光学系统进行像质优化及分析，得到比较好的成像效果，并对框架角及盲区开展分析计算。通过红外成像导引头调试试验截图、外场调试试验截图及挂飞试验截图的分析，得出成像质量良好，能够实现搜索、跟踪目标的功能。此红外导引头光学系统的设计，可以后类似系统的开发提供参考。
-
![]()
图 5 行缓冲将10像素80bit输入转成串行1像素8bit
Figure 5. Line buffer converts 10-pixel 80-bit parallel input into serial 1-pixel 8-bit
表 1 处理能力比较
Table 1 Comparison of processing power
Processor Number of sub-apertures The data volume of the sensor Processing delay 2002 Institute of Optics and Electronics, Chinese Academy of Sciences[1] 17DSP 61 7.8MB/s 340μs 2011 Institute of Optics and Electronics, Chinese Academy of Sciences[5] FPGA 177 No data 338μs 2014 Technische Universität Ilmenau[6] FPGA 120 43.3MB/s 1050μs 2015 Institute of Optics and Electronics, Chinese Academy of Sciences[2] Computer 949 88.2MB/s 163μs 2016 Changchun Institute of Optics and Mechanics[4] GK104 GPU+CPU 1600 115.2MB/s 18μs 2018 Changchun Institute of Optics and Mechanics[7] FPGA 349 113.5MB/s 198μs 2018 Indian Institute of Astrophysics[8] FPGA 2500 No data 24-26μs System and method of this paper FPGA 3136 580MB/s μs 
下载: 导出CSV
-
[1] 周璐春, 王春鸿, 李梅, 等. 基于FPGA技术的波前斜率处理方法[J]. 光电工程, 2002, 29(3): 28-31. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC200203007.htm ZHOU Luchun, WANG Chunhong, LI Mei, et al. A wavefront slope processing method based on FPGA technique[J]. Opto-Electronic Engineering, 2002, 29(3): 28-31 https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC200203007.htm
[2] 樊志华, 王春鸿, 姜文汉. 基于累加器的哈特曼-夏克波前斜率处理器[J]. 光学精密工程, 2011, 19(3): 501-507. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201103002.htm FAN Zhihua, WANG Chunhong, JIANG Wenhan. Accumulator -based wavefront slope processor for Shack-Hartmann sensors[J]. Optics and Precision Engineering, 2011, 19(3): 501-507. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201103002.htm
[3] 陈善球, 刘超, 许冰, 等. 通用自适应光学波前实时处理机的设计[J]. 中国激光, 2015, 42(12): 1212001 https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ201512036.htm CHEN Shanqiu, LIU Chao, XU Bing, et al. Design of generic adaptive optics wave-front real-time processor[J]. Chinese Journal of Lasers, 2015, 42(12): 1212001 https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ201512036.htm
[4] 李大禹. 基于多GPU的液晶自适应光学波前处理器[J]. 液晶与显示, 2016, 31(5): 491-496. https://www.cnki.com.cn/Article/CJFDTOTAL-YJYS201605011.htm LI Dayu. Liquid crystal adaptive optics wavefront processor based on multi-GPU[J]. Chinese Journal of Liquid Crystals and Displays, 2016, 31(5): 491-496. https://www.cnki.com.cn/Article/CJFDTOTAL-YJYS201605011.htm
[5] 贾建禄, 王建立, 赵金宇. 基于FPGA的自适应光学系统波前处理机[J]. 光学精密工程, 2011, 19(8): 1716-1722. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201108003.htm JIA Jianlu, WANG Jianli, ZHAO Jinyu, et al. Adaptive optical wave-front processor based on FPGA[J]. Optics and Precision Engineering, 2011, 19(8): 1716-1722. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201108003.htm
[6] Mauch S, Reger J. Real-time spot detection and ordering for a Shack-Hartmann wavefront sensor with a low-cost FPGA[J]. IEEE Transactions on Instrumentation and Measurement, 2014, 63(10): 2379-2386. DOI: 10.1109/TIM.2014.2310616
[7] 王建立. 349单元自适应光学波前处理器[J]. 光学精密工程, 2018, 26(5): 1007-1013. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201805001.htm WANG Jianli. Three hundred and forty-nine unit adaptive optical wave front processor[J]. Optics and Precision Engineering, 2018, 26(5): 1007-1013. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201805001.htm
[8] Avinash Surendran, Mahesh P Burse, Ramaprakash A N, et al. Scalable platform for adaptive optics real-time control (SPARC) part 2: field programmable gate array (FPGA) implementation and performance[J]. Journal of Astronomical Telescopes, Instruments, and Systems, 2018, 4(3): 039002. DOI: 10.1117/1.JATIS.4.3.039002.full
计量
- 文章访问数: 243
- HTML全文浏览量: 116
- PDF下载量: 60
