1. 引言
无人机凭借其灵活机动、小型化和智能化的特点,已广泛应用于农业管理[1]、环境监测[2]、军事侦察[3] [4]等领域。然而,传统无人机载相机存在视场小[5]-[7]、体积大等问题,获取的仅是三维场景的二维投影信息,这限制了其在测速与定位、深度估计等任务中的应用[8]。为突破这些局限,部分研究者提出了多相机系统[9]。然而,该系统不仅体积庞大、数据处理复杂,而且在实际应用中多相机的同步问题仍面临较大挑战。
曲面复眼与平面复眼系统相比,曲面复眼系统通过子眼阵列构建非平面的成像结构,能够有效减少像差并扩大视场范围。虽然鱼眼镜头[10]也能实现大视场成像,但同时伴随着严重的畸变和轴外像差问题。基于这些问题,诸多科研人员开始广泛地对曲面复眼系统进行研究。
习啸天[11]等人设计了一款复眼相机系统,其最大视场为122˚,在对地距离1000 m的高空拍摄时,地面分辨率可达0.8 m,系统总体积为123 mm × 123 mm × 125 mm。许黄蓉[12]等人设计了一个采用六边形环带排布的曲面大孔径复眼系统,总视场达到98˚ × 98˚,在对地距离500 m的高空拍摄时,地面分辨率达到1 m,系统体积为123 mm × 123 mm × 195 mm。尽管上述系统实现了大视场,但在地面分辨率和体积控制方面仍存在一定不足。
针对这一问题,本文聚焦于无人机载相机对大视场、高分辨率及小型化的需求,研究在小型化设计的前提下如何兼顾大视场与高分辨率的性能优化。
2. 系统设计
2.1. 子眼阵列的设计
曲面复眼相较于平面复眼更贴近自然界复眼的结构特性,其高灵敏度的优势源于子眼间的视场重叠。相较于平面排列,曲面排列的子眼系统能够实现更大的视场,因此,子眼的排列方式成为研究的关键。本文采用六边形排布方式,以进一步提高光能利用率[13],为大视场、高分辨率成像提供优化设计方案。
图1中R是子眼曲面排列基底曲率半径,d是相邻子眼之间的间距,D是子眼孔径,φ是相邻子眼之间的夹角,2ω是子眼的视场。
Figure 1. Schematic diagram of the surface arrangement of the sub-eye
图1. 子眼曲面排列示意图
相邻子眼间的光轴夹角为:
(1)
相邻子眼间的视场重叠率为:
(2)
式中σ为视场重叠率。当
时,
;当
时,在无穷远处,三个相邻子眼间中间子眼的视场会被相邻两个子眼完全覆盖[14],所以取
。
复眼光学系统总视场从下式可得:
(3)
其中,
为复眼光学系统总视场,
为六边形环带层数。
子眼个数
与六边形环带层数关系为[13]:
(4)
对于二次成像系统,系统焦距与子系统焦距的关系为:
(5)
系统采用1/3" ICX445作为接收器,其感光表面对角线为6 mm,像素尺寸大小为3.75 um。地面分辨率公式为:
(6)
其中
为系统焦距,
为像元尺寸,
为对地高度,
为地面分辨率。考虑到小型化的需求,令子眼焦距为
,中继系统焦距为
,作为初始结构。
根据设计指标要求,参数结合公式(5)和公式(6),系统总体参数如表1所示。
Table 1. Design indicators for compound ophthalmic systems
表1. 复眼光学系统设计指标
系统参数 |
设计指标 |
总视场θ |
112˚ |
波长范围 |
可见光 |
F数 |
3.5 |
子眼视场ω |
6˚ |
子眼孔径D |
2 mm |
光轴夹角φ |
10˚ |
六边形环带层数S |
6 |
子眼个数n |
91 |
子眼焦距f1 |
12 mm |
材料库 |
成都光明 |
根据奈奎斯特计算公式可得复眼系统的截止频率为:
(7)
子眼由三块镜片构成,图2和图3分别展示了优化后子眼的结构;点列图及调制传递函数(Modulation Transfer Function, MTF)图。
Figure 2. Sub-eye structure diagram
图2. 子眼结构图
(a) MTF图 (b) 点列图
Figure 3. Simulation analysis diagram of the optical system for sub-eye imaging
图3. 子眼成像光学系统仿真分析图
经过仿真分析,从图3(a)可看出MTF截止频率处的值均大于0.3,且子眼各视场的RMS半径均小于单个像元尺寸,能够满足光学成像系统与接收器的匹配要求。
2.2. 中继转像系统的设计
受到工艺影响,曲面接收器目前难以制造,且成本较高,所以选择中继转像系统将曲面微透镜阵列形成的曲面像转为平面像,便于平面接收器接收。
通过Zemax全局优化算法,在子眼阵列MTF值最大化和畸变最小化的约束条件下,最终确定了在曲率半径为R = 30 mm的曲面基底上,以六边形结构排列子眼阵列。子眼阵列在焦曲面上形成了曲面像,但由于技术限制,该曲面像无法直接被接收器接收。因此,需要设计一个中继转像系统,将子眼阵列形成的曲面像转换为平面像,以便接收器进行有效成像。为了确保成像质量,中继转像系统的物面曲率半径必须与子眼阵列的焦曲面曲率半径一致,同时其物面口径应大于子眼阵列焦曲面的口径,从而确保从子眼阵列发出的光线能够完全进入中继转像系统。
Figure 4. Schematic diagram of the structure of the relay image conversion system
图4. 中继转像系统结构示意图
(a) MTF图 (b) 点列图
Figure 5. Simulation analysis diagram of the relay image system
图5. 中继转像系统仿真分析图
中继转像系统通常可视为近物距的鱼眼镜头,其前1至2片镜片为负弯月镜片,主要作用是压缩视场。当光束以较大入射角照射到前光组的光学面时,经过光学系统成像后,子午和平面内的聚焦位置与波阵面参数可能完全不一致。因此,传统的赛德尔像差理论主要适用于轴对称光学系统的像差分析,而不适用于此类光学系统的像差计算,需要运用平面对称光学系统的波像差理论来计算系统的初始结构[15]。
基于平面对称光学系统的波像差理论,本文中继系统采用负弯月透镜与双高斯转像结构的复合架构作为初始光学拓扑。在优化过程中,通过引入非球面并策略性地将其置于光阑远端实现了系统轴向尺寸压缩与视场边缘像差的协同优化。负弯月透镜不仅能有效压缩视场,还能通过反向曲率设计补偿双高斯结构产生的场曲。此外,配合双高斯透镜的对称性对系统彗差进行了有效控制。
经过仿真优化分析,图4展示了中继系统的结构图,图5则展示了中继系统的成像质量。
从图5(a)看出中继转像系统的MTF在截止频率处的值均大于0.3,且各视场RMS半径值均小于像素尺寸,满足光学成像系统与接收器的匹配要求。
2.3. 复眼光学系统组合
在完成子眼阵列和中继转像系统的设计后,需要将其拼接成完整的复眼光学系统。在拼接过程中,必须遵循光瞳衔接原则,即子眼阵列的出瞳与中继转像系统的入瞳重合。由于系统采用离轴结构,拼接时不仅需要调整子眼的倾斜参数,还需对其偏心参数进行调整,确保每个子眼发出的主光线在进入中继转像系统前表面时,能够与中继转像系统中物面上各视场发出的主光线进入前表面的位置一致。
考虑到复眼光学系统中边缘视场与中心视场的差异较大,在分析系统像质时,需要对每个子眼的像质进行独立评估。在优化软件中,将每个子眼设为独立的多重结构,并确保在优化过程中,每个结构都能单独优化,同时保持除倾斜、偏心及与中继转像系统前表面距离不一致外,其余参数的一致性。最终,通过构建评价函数,对各个结构在MTF截止频率处的值和畸变进行约束,从而完成优化过程。
拼接后的复眼光学系统主阵列如图6所示:
复眼光学系统是二次成像系统,整个系统是由多个子眼拼接而成,需要对各位置子眼单独进行像质评估。由于系统是对称结构,且边缘视场像差较为敏感,本文仅展示一侧的中心视场与边缘视场的像质图。其中图7为系统MTF图,图8展示了系统的RMS半径值,图9展示了系统的场曲/畸变值。
Figure 6. Schematic diagram of the structure of the compound vision system
图6. 复眼光学系统结构示意图
复眼光学系统中心与边缘子眼MTF截止频率处的值如图7所示。从图7中可得,各子眼MTF在截止频率处的值均大于0.3。0˚和50˚视场子眼点列图如图8所示。各子眼全视场RMS半径值均小于接收器像元尺寸。场曲/畸变如图9所示,各结构子眼畸变值均小于4%,完全符合光学成像系统与接收器的匹配要求。
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 7. MTF simulation analysis diagram of compound vision system
图7. 复眼光学系统MTF仿真分析图
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 8. Diagram of the simulation analysis points of the double-vision optical system
图8. 复眼光学系统仿真分析点列图
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 9. Simulation and analysis of field curvature/distortion diagram of a compound ophthalmic system
图9. 复眼光学系统仿真分析场曲/畸变图
3. 公差分析
Table 2. Tolerance analysis results
表2. 公差分析结果
蒙特卡洛样本百分比/% |
几何MTF平均值 |
90 |
0.193 013 08 |
80 |
0.211 302 30 |
50 |
0.235 151 84 |
20 |
0.256 399 06 |
10 |
0.261 649 66 |
在光学系统的设计过程中,不仅要满足光学成像系统与接收器的匹配要求,还要满足产品生产加工的需求,公差分析的目的在于了解系统在多大公差范围内,仍然能够满足成像要求。本文选用“几何MTF平均”的公差评价方法,用“灵敏度”模式,“蒙特卡洛运行数”选择“200”,根据光学成像系统公差表Q5等级对系统边缘结构进行公差分析。该系统边缘结构名义值为0.27541188,具体结果如表2所示。分析表2可得由加工或装配等因素造成误差引起MTF截止频率处的值下降后,仍然有90%的产品其MTF截止频率处的值高于0.19301308,从结果可以看出大多数样品的MTF截止频率处的值大于0.15,满足设计指标。
公差分配选用光学成像系统公差表中Q5等级,具体参数如表3:
Table 3. Tolerance distribution range
表3. 公差分配范围
名称 |
参数 |
折射率 |
±0.001 |
不规则度 |
±0.5 |
阿贝误差(%) |
±0.1 |
厚度(mm) |
±0.0375 |
表面偏心(mm) |
±0.01 |
元件偏心(mm) |
±0.0375 |
表面倾斜(`) |
±1 |
元件倾斜(`) |
±1 |
4. 无热化分析
在完成公差分析并确定系统对制造装配误差的容忍度后,我们发现光学性能的稳定性不仅取决于静态的几何误差控制,动态环境因素特别是温度变化也会影响光学系统的成像质量。为构建具有环境适应性的光学系统,需要将热力学变量纳入综合优化框架,通过光学被动消热差的方法保证光学系统在满足温度范围内仍然满足成像要求。根据无人机产品的工作环境温度标准,无人机应能够在−20℃至50℃的温度范围内正常工作。为了满足无人机载光学系统在宽温度范围内的稳定性要求,本文通过Zemax热–光学耦合仿真对系统进行无热化分析,仿真温度范围覆盖−20℃至50℃。当温度为基准温度20℃时,系统的无热化分析仿真结果如图7和图8所示。
当温度为−20℃时,系统的无热化分析仿真结果如图10和图11所示。
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 10. Simulated analysis of MTF diagram of compound vision system at −20˚C
图10. 复眼光学系统−20˚C仿真分析MTF图
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 11. Diagram of the −20˚C simulation analysis point of the compound vision system
图11. 复眼光学系统−20˚C仿真分析点列图
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 12. Simulated analysis of MTF diagram of compound vision system at 50˚C
图12. 复眼光学系统50˚C仿真分析MTF图
(a) 中心0˚视场 (b) 边缘50˚视场
Figure 13. Diagram of the 50˚C simulation analysis point of the compound vision system
图13. 复眼光学系统50˚C仿真分析点列图
当温度为50℃时,系统的无热化分析仿真结果如图12和图13所示。
经过Zemax无热化仿真分析验证,该系统在温度范围为−20℃至50℃时,成像质量稳定。
5. 结论
本文设计了一款适用于无人机的小型化、大视场、高分辨率曲面复眼相机系统,通过适配转像系统将曲面子眼阵列形成的焦曲面像转化为平面像,便于接收器接收。经拼接与系统优化后,系统总体积为60 mm × 60 mm × 72 mm,焦距为3.75 mm,满足对地距离500 m实现0.5 m地面分辨率的要求。系统全视场在奈奎斯特截止频率处的MTF值高于0.3,子眼全视场RMS半径小于像元尺寸,畸变值小于4%。公差分析结果表明,系统经过合理的公差分配后,90%的样品在截止频率处的MTF值仍高于0.15,成像质量良好。无热化分析结果表明,系统在−20℃至50℃的温度范围内成像质量稳定。该系统结构简单、加工难度低,可广泛应用于大视场成像、军事目标探测与识别、导航定位等多个领域。
基金项目
国家自然基金青年基金(62303441);中国科学院青年创新促进会(2021233);上海市优秀学术带头人(22XD1424500)。
NOTES
*通讯作者。