-
地球静止轨道卫星成像光谱仪在气象观测,资源调查,环境监测等领域具有重大的应用价值。而覆盖可见光到长波红外的光谱仪器具有昼夜观测的技术优势[1]。地球静止轨道凝视型面阵成像光谱仪对其时间分辨率、空间分辨率和成像幅宽,遥感成像动态范围提出了更高的要求,以期获得更加丰富的目标层次及细节信息,以便对探测目标信息进行有效地获取及识别,由于卫星资源有限,大数据量会产生存储、传输及处理的困难。压缩感知技术以其先压缩再采样的特点,在现有信号传输水平下可以有效提高获取图像的效率,降低数据采集的压力。在此背景下,利用地球静止轨道平台可以长期驻留固定区域上空的特点[2-5],结合压缩感知光谱成像技术,对全谱段成像光学系统进行分析,并对其光学系统进行设计。
-
将压缩感知理论应用于光谱成像技术中,可以实现对光谱图像信息的压缩采集,解决了卫星数据存储、传输、处理的问题,其原理如图1所示。待测光经过物镜成像于离散编码模板上,离散编码模板将得到的二维信息调制准直后,由色散元件色散为光谱带,会聚后,探测器获取经过压缩后的二维离散光谱混叠信息,通过优化解混后即可高概率获得原始光谱信息,压缩感知成像光谱技术获得的采样数据量远小于传统成像光谱技术的数据量,可将将卫星信号数据量降低两个数量级[6-8]。
根据此模型提出本文的光学系统结构形式,根据光学系统全谱段的要求,需将系统分为可见光、中波红外、长波红外3个谱段光路分别进行设计[7],采用分色片分谱段,选择DMD作为离散编码模板,可见光所选DMD像元大小为7.6 μm,采用一个DMD像元与4个探测器像元大小进行匹配,中波红外光路所选DMD大小为7.6 μm,采用一个DMD像元与一个探测器像元大小匹配,长波红外光路所选DMD像元大小为10.8 μm,采用一个DMD像元与1个探测器像元大小进行匹配,DMD的开合角度为±12°,通过中继镜组的放大倍率匹配编码模板与探测器像元面积的不同,目前常用DMD的像元数受芯片工艺限制,无法与探测器实现完美的匹配,当探测器分辨率高于编码DMD的分辨率时,可通过对DMD像元进行灰度分级的办法,保证分级后的编码像元与探测器像元一一对应,即可保证系统的整体分辨率。由于DMD带来的额外像差可由像面倾斜进行补偿,基于DMD的目标模拟器应用于长波红外波段时,照明光束需要选取合适的入射角以减少DMD的衍射效应对系统成像对比度造成的影响[9-11]。整体光学系统的模型如图2所示,其由物镜,分光片,一次成像系统,准直系统,色散元件,DMD,以及光谱信息采集处理系统七部分组成,待测光经过物镜入射到分光系统进行波段选择,经成像系统成像于编码模板上,由准直系统准直为平行光,在色散元件表面色散为光谱带,经会聚系统会聚在探测器上。光学系统覆盖可见光、中波红外、长波红外,谱段较宽。根据压缩感知成像光谱系统模型,需对系统进行模块化设计、可降低整体设计的难度,在保证各模块的成像质量后,将各模块整合最终对整体进行设计[12-20]。
-
本文光学系统分为3个谱段,各谱段设计要求如表3.1所示
表 3.1 光学系统设计要求
可见光 中波红外 长波红外 空间分辨率/m 40 400 625 幅宽/km 400×400 400×400 320×320 光谱分辨率/nm 20 50 150 根据空间分辨率以及幅宽的要求,所选探测器参数如表3.2所示。
表 3.2 光学系统探测器参数
可见光 中波红外 长波红外 像元数 10000×10000 1024×1024 512×512 像元尺寸/ μm 6.4×6.4 15×15 24×24 根据公式1对系统,焦距及视场角进行计算
$$B=2H{\rm{tan}}\omega =\frac{LH}{f}$$ (1) 其中H为卫星的轨道高度,B为相机的幅宽,L为像面的长度,a为像元尺寸大小,ω为半视场角,
$f$ 为系统焦距。相机镜头孔径D的大小影响相机的分辨本领,调制传递函数(MTF),曝光量,信噪比和相机尺寸。相机孔径D的大小需考虑探测器的信噪比、噪声等效温差、噪声等效反射率差的要求。综合考虑后,光学系统孔径为700 mm。
表 3.3 光学系统参数
可见光 中波红外 长波红外 系统孔径/mm 700 视场角2ω/° 0.66×0.66 0.66×0.66 0.52×0.52 系统焦距/mm 4500 1350 1375 物镜设计为同轴反射无焦系统,可以实现无色差设计以及卫星长时间凝视成像的轨道的震动补偿,其放大倍率为5x,要求系统遮拦比小于1/3。同时根据DMD像元大小与探测器像元大小计算不同光路准直光路的焦距,准直光路焦距与会聚光路焦距保持一致。选择光栅作为系统的色散元件,光栅的刻痕数直接决定系统的光谱分辨率。根据公式(2)(3)与各谱段会聚光路焦距计算光栅刻线数,如表3.4
$$\delta \lambda =\frac{d\lambda }{d\theta }\Delta \theta =\frac{d{\rm{cos}}\theta }{m}\frac{\lambda }{Nd{\rm{cos}}\theta }=\frac{\lambda }{mN}$$ (2) $$\theta =2{\rm{arctan}}\frac{\delta }{2f}$$ (3) 其中
$\delta $ 为像元的大小,$\delta \lambda $ 为系统的光谱分辨率,各光学子系统参数如表3.3所示。表 3.4 子系统参数
可见光一次会聚光路/准直光路 中波红外一次会聚光路/准直光路 长波红外一次会聚光路/准直光路 系统焦距/mm 900/534 1350/137 1375/124 光栅线对数lp/mm 170 50 45 设计结果要求一次会聚系统,准直系统,会聚系统各部分像质良好,对于高分辨率的成像光谱仪系统,二级光谱的存在将严重影响最终的成像质量,二级光谱的色差大小与系统焦距成正比,因此对于长焦距的可见光路系统,需对其二级光谱进行优化及评价。由于系统空间维很大,必然存在谱线弯曲和色畸变,谱线弯曲是指不同波长弯曲图像与直线的偏离程度,色畸变是由光谱成像系统对不同波长的像放大率不同造成的,严重时,会影响成像质量,谱线弯曲过大还会降低系统的空间分辨率,为获得准确的多光谱图像,谱线弯曲以及色畸变需控制在1/2个像元以内。在中波、长波光路中选择制冷型红外探测器,因此需对进行冷光阑匹配。
-
基于上述参数分别对3个谱段进行光学系统设计,其整体结构如图3所示。
物镜选择以卡塞格林为基础的同轴三反无焦系统,在二镜,三镜中加入反射镜以折转光路,使系统紧凑,同时在主镜与三镜中引入高阶非球面在保证系统放大率和遮拦比的情况下,像质达到最佳。其设计结果如图4所示,其放大倍率为5x,遮拦比为0.185。
对可见光、中波红外、长波红外系统整体要求结构紧凑,对各光路不同系统进行模块化设计[15-20],可以降低系统的设计难度。
可见光路的优化中,约束一次会聚光路、准直光路、会聚光路焦距分别为900 mm、534 mm、534 mm,约束各镜片厚度以及系统整体长度,对于长焦距系统,对二级光谱的优化选择阿贝数相差较大、部分色散相同或相近的玻璃对其进行优化。采用棱镜-光栅结构作为可见光路的分光元件,补偿系统的谱线弯曲,通过真实光线的追迹控制系统的色畸变大小,利用像面倾斜的方式校正轴向色差并达到平谱面的目的,在对各模块优化后,对整体结构进行优化,优化结果可见光会聚光路各视场MTF在波长500 nm、700 nm、900 nm时,在78.125 lp/mm处分别高于0.674、0.562、0.455,点列斑83%以上位于艾里斑内,弥散斑半径最大分别为5.231 μm、4.694 μm、4.113 μm,小于一个探测器像元尺寸。在不同视场下光谱变化下20 nm可以清晰分辨,满足20 nm的光谱分辨率,其中不同颜色光斑分别代表900 nm、880 nm、720 nm、700 nm、520 nm、500 nm波长在光学系统的点列图大小,可以看出光学系统在各视场均满足20 nm的光谱分辨率,可见光系统的二级光谱小于0.05 mm,谱线弯曲很小,色畸变最大值为2.753 μm,小于1/2个像元大小,如图5~图15所示。
中波红外的优化中,约束一次会聚光路、准直光路、会聚光路焦距分别为270 mm、137 mm、137 mm,约束各镜片厚度以及系统整体长度,并对部分表面引入二次曲面以及高阶非球面,针对中波红外系统的谱线弯曲以及色畸变,采用离轴透镜的方式对其补偿,先优化得到一次成像未离轴时的光学系统,然后根据系统存在的谱线弯曲量,确定透镜的离轴量,最后适当改变透镜的半径完成谱线弯曲的校正,通过真实光线的追迹控制系统的色畸变大小,采用二次成像的方式对光学系统的冷光阑进行匹配,控制系统的焦距以及各视场的边缘光线出瞳像高对系统进行优化。优化后各模块成像质量良好、最终得到的中波红外光学系统在波长3500 nm、3800 nm、4100 nm处的MTF在33.3 lp/mm时分别高于0.607、0.566、0.518。弥散斑半径最大分别为5.419 μm、5.922 μm、6.678 μm。在不同视场下光谱变化50 nm时各波长点列图光斑可以清晰的分辨,其中不同颜色光斑分别代表4100 nm、4050 nm、3800 nm、3750 nm、3550 nm、3500 nm波长光学在光学系统的点列图大小,可以看出光学系统满足50 nm的光谱分辨率,谱线弯曲很小,色畸变最大值为5.48 μm,小于1/2个像元,如图16~图25所示。
图 14 不同视场下可见光路光谱分辨率
Figure 14. Spectral resolution of the visible spectrum in different fields of view
长波红外的优化中,约束一次会聚光路、准直光路、会聚光路焦距分别为275 mm、126 mm、126 mm,约束各镜片厚度以及系统整体长度并对部分表面引入二次曲面以及高阶非球面,针对长波红外系统的谱线弯曲以及色畸变,与中波红外系统相似,采用离轴透镜的方式对其优化,通过真实光线的追迹控制系统的色畸变的大小,采用二次成像的方式对光学系统的冷光阑进行匹配,约束系统焦距以及各视场光线出瞳像高对系统进行优化。优化后各模块成像质量良好、最终得到的长波红外光学系统在波长7700 nm、8600 nm、9500 nm处的MTF在20.8 lp/mm时分别高于0.568、0.536、0.498。弥散斑半径最大分别为10.16 μm、4.007 μm、9.163 μm。在不同视场下光谱变化150 nm时各波长点列图光斑可以清晰的分辨,其中不同颜色光斑分别代表9500 nm、9350 nm、8600 nm、8450 nm、7850 nm、7700 nm波长光学在光学系统各视场的点列图大小,可以看出光学系统满足150 nm的光谱分辨率,谱线弯曲很小,色畸变最大值为9.35 μm,小于1/2个像元大小,如图26~图35所示。
-
为了获取地球静止轨道的高动态光谱数据,设计了一种基于压缩感知的全谱段成像光学系统。该系统通过物镜无焦化设计以满足卫星轨道震动补偿的需求并使用分色片将光学系统谱段细化,以实现全谱段的探测,经过优化设计,系统结构紧凑,点列图能量集中于艾里斑,各部分MTF均接近衍射极限,像质良好,光谱分辨率分别为可见光20 nm、中波红外50 nm、长波红外150 nm。满足设计要求。
Optical design of visual and infrared large-aperture imaging system based on a space-based platform
-
摘要: 针对用于地球静止轨道卫星的遥感面阵快照式成像光谱仪传输数据量过大引起的数据传输困难,信号采集处理时间长的问题,利用地球静止轨道平台可以长期驻留固定区域上空的特点,提出采用压缩感知的大口径全谱段快照式光谱仪的方案。对其光学系统结构进行设计,并对相关参数进行了计算。物镜采用同轴三反式无焦系统,用分色片对系统分光,经过对各系统优化处理,最终获得了幅宽为400 km◊400 km,可见光分地面像元分辨率为40 m,中波红外地面像元分辨率为400 m,长波红外地面像元为625 m的光学系统,设计结果MTF可见光路在78.125 lp/mm高于0.455,中波红外光路在33.3 lp/mm高于0.518,长波红外光路在20.8 lp/mm高于0.498;系统光谱分辨率分别为可见光20 nm、中波红外50 nm、长波红外光150 nm;可见光路二级光谱小于0.05 mm,设计结果具有良好的成像质量,各部分光学系统成像质量接近衍射极限,设计结果满足应用和指标需求。Abstract: To solve the data transmission difficulties and long-signal acquisition and processing time problem caused by excessive data transmission of the Geostationary Orbit array staring spectrometer, a scheme for a large-aperture visual and infrared snapshot spectrometer based on compressive sensing is proposed, which takes advantage of the fact that a geostationary orbit platform can stay over the fixed area for a long time. This paper analyzes the physical model of compressive sensing spectral imaging. The structure of the optical system is designed, and the relevant parameters are calculated. The objective lens uses a coaxial three-mirror afocal optical system, and dichroic films are used to split the spectrum. The relevant parameters were calculated according to requirements and the system was designed using optical software. After optimization, the optical system was shown to have a width of 400 km◊400 km, a visible area of 40 m GSD, an MWIR area of 400 m GSD, and an LWIR area of 625 m GSD. The results show that the MTF of the visible area is higher than 0.455 at 78.125 lp / mm. In mid-wave infrared, the MTF is higher than 0.518 at 33 lp / mm, and the MTF is higher than 0.498 at 20.8 lp / mm in long-wave infrared. The spectral resolution of the visible light is 20 nm, the mediumwave infrared region’s resolution is 50 nm, and the long-wave infrared spectrum resolution is 150 nm. The second-order spectrum of the visual area is less than 0.05 mm. The optical system has good imaging performance and the imaging quality of each part of the optical system is close to the diffraction limit, which meets the needs of applications and indicators.
-
表 3.1 光学系统设计要求
可见光 中波红外 长波红外 空间分辨率/m 40 400 625 幅宽/km 400×400 400×400 320×320 光谱分辨率/nm 20 50 150 表 3.2 光学系统探测器参数
可见光 中波红外 长波红外 像元数 10000×10000 1024×1024 512×512 像元尺寸/ μm 6.4×6.4 15×15 24×24 表 3.3 光学系统参数
可见光 中波红外 长波红外 系统孔径/mm 700 视场角2ω/° 0.66×0.66 0.66×0.66 0.52×0.52 系统焦距/mm 4500 1350 1375 表 3.4 子系统参数
可见光一次会聚光路/准直光路 中波红外一次会聚光路/准直光路 长波红外一次会聚光路/准直光路 系统焦距/mm 900/534 1350/137 1375/124 光栅线对数lp/mm 170 50 45 -
[1] 戴立群, 唐绍凡, 徐丽娜, 等. 从可见光到热红外全谱段探测的星载多光谱成像仪器技术发展概述[J]. 红外技术,2019,41(2):107-117. DAI L Q, TANG SH F, XU L N, et al. Development overview of space-borne multi-spectral imager with band range from visible to thermal infrared[J]. Infrared Technology, 2019, 41(2): 107-117. (in Chinese [2] 马文坡, 练敏隆. “高分四号”卫星凝视相机的技术特点[J]. 航天返回与遥感,2016,37(4):26-31. doi: 10.3969/j.issn.1009-8518.2016.04.004 MA W P, LIAN M L. Technical characteristics of the staring camera on board GF-4 satellite[J]. Spacecraft Recovery &Remote Sensing, 2016, 37(4): 26-31. (in Chinese doi: 10.3969/j.issn.1009-8518.2016.04.004 [3] 陶家生, 孙治国, 孙英华, 等. 静止轨道高分辨率光学遥感探索[J]. 光电工程,2012,39(6):1-6. doi: 10.3969/j.issn.1003-501X.2012.06.001 TAO J SH, SUN ZH G, SUN Y H, et al. Exploration of high resolution optical remote sensing of the geostationary orbit[J]. Opto-Electronic Engineering, 2012, 39(6): 1-6. (in Chinese doi: 10.3969/j.issn.1003-501X.2012.06.001 [4] 罗秀娟, 刘辉, 张羽, 等. 地球同步轨道暗弱目标地基光学成像技术综述[J]. 中国光学,2019,12(4):753-766. doi: 10.3788/co.20191204.0753 LUO X J, LIU H, ZHANG Y, et al. Review of ground-based optical imaging techniques for dim GEO objects[J]. Chinese Optics, 2019, 12(4): 753-766. (in Chinese doi: 10.3788/co.20191204.0753 [5] 黄思婕. 地球静止轨道大动态范围信息获取技术研究[D]. 上海: 中国科学院研究生院(上海技术物理研究所), 2015. HUANG S J. Research on the technology of geosynchonous orbit high dynamic range information acquisition[D]. Shanghai: Shanghai Institute of Technical Physics of the Chinese Academy of Sciences, 2015. (in Chinese). [6] 刘铭鑫. 基于压缩感知的编码孔径光谱成像技术研究[D]. 长春: 中国科学院大学(中国科学院长春光学精密机械与物理研究所), 2019. LIU M X. Research on coded aperture spectral imaging technology based on compressed sensing[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, 2019. (in Chinese). [7] 钱路路. 计算光谱成像技术研究[D]. 合肥: 中国科学技术大学, 2013. QIAN L L. Research on computational imaging spectroscopy[D]. Hefei: University of Science and Technology of China, 2013. (in Chinese). [8] 闫歌, 许廷发, 马旭, 等. 动态测量的高光谱图像压缩感知[J]. 中国光学,2018,11(4):550-559. doi: 10.3788/co.20181104.0550 YAN G, XU T F, MA X, et al. Hyperspectral image compression sensing based on dynamic measurement[J]. Chinese Optics, 2018, 11(4): 550-559. (in Chinese doi: 10.3788/co.20181104.0550 [9] 韩庆, 王健, 熊峥, 等. 用于长波红外目标模拟器的DMD衍射特性分析[J]. 红外与激光工程,2017,46(5):93-99. HAN Q, WANG J, XIONG ZH, et al. Diffraction characteristics analysis for DMD-based scene projectors in the long-wave infrared[J]. Infrared and Laser Engineering, 2017, 46(5): 93-99. (in Chinese [10] 吕伟振, 刘伟奇, 魏忠伦, 等. 基于DMD的高动态范围成像光学系统设计[J]. 红外与激光工程,2014,43(4):1167-1171. doi: 10.3969/j.issn.1007-2276.2014.04.026 LV W ZH, LIU W Q, WEI ZH L, et al. Design of high dynamic range imaging optical system based on DMD[J]. Infrared and Laser Engineering, 2014, 43(4): 1167-1171. (in Chinese doi: 10.3969/j.issn.1007-2276.2014.04.026 [11] 孙永强, 胡源, 王月旗, 等. 数字微镜器件在会聚成像光路中的像差分析[J]. 光学学报,2019,39(3):0311001. doi: 10.3788/AOS201939.0311001 SUN Y Q, HU Y, WANG Y Q, et al. Analysis on aberration of digital micromirror device in convergent imaging path[J]. Acta Optica Sinica, 2019, 39(3): 0311001. (in Chinese doi: 10.3788/AOS201939.0311001 [12] 刑振冲. 灵巧型长焦多波段共口径光学系统的研究[D]. 长春: 中国科学院大学(中国科学院长春光学精密机械与物理研究所), 2018. XING ZH CH. Research on miniature telefocal multiband common aperture optical system[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, 2018. (in Chinese). [13] 张天一, 朱永田, 侯永辉, 等. LAMOST高分辨率光谱仪研制[J]. 中国光学,2019,12(1):148-155. doi: 10.3788/co.20191201.0148 ZHANG T Y, ZHU Y T, HOU Y H, et al. Construction of a LAMOST high resolution spectrograph[J]. Chinese Optics, 2019, 12(1): 148-155. (in Chinese doi: 10.3788/co.20191201.0148 [14] 曹佃生, 林冠宇, 杨小虎, 等. 紫外双光栅光谱仪结构设计与波长精度分析[J]. 中国光学,2018,11(2):219-230. doi: 10.3788/co.20181102.0219 CAO D SH, LIN G Y, YANG X H, et al. Structure design and wavelength accuracy analysis of ultraviolet double grating spectrometer[J]. Chinese Optics, 2018, 11(2): 219-230. (in Chinese doi: 10.3788/co.20181102.0219 [15] 潘君骅. 光学非球面的设计、加工与检验[M]. 苏州: 苏州大学出版社, 2004. PAN J H. The Design, Manufacture and Test of the Aspherical Optical Surfaces[M]. Suzhou: Suzhou University Press, 2004. (in Chinese). [16] 孙武, 韩诚山, 吕恒毅, 等. 推扫式多光谱遥感相机动态范围拓展方法[J]. 中国光学,2019,12(4):905-912. doi: 10.3788/co.20191204.0905 SUN W, HAN CH SH, LV H Y, et al. Dynamic range extending method for push-broom multispectral remote sensing cameras[J]. Chinese Optics, 2019, 12(4): 905-912. (in Chinese doi: 10.3788/co.20191204.0905 [17] 袁立银, 谢佳楠, 侯佳, 等. 紧凑型红外成像光谱仪光学设计[J]. 红外与激光工程,2018,47(4):186-191. YUAN L Y, XIE J N, HOU J, et al. Optical design of compact infrared imaging spectrometer[J]. Infrared and Laser Engineering, 2018, 47(4): 186-191. (in Chinese [18] 胡斌, 黄颖, 马永利, 等. 高分辨率红外成像仪五反无焦主系统设计[J]. 红外与激光工程,2016,45(5):0518001. doi: 10.3788/irla201645.0518001 HU B, HUANG Y, MA Y L, et al. Design of five-mirror afocal principal system for high spatial resolution infrared imager[J]. Infrared and Laser Engineering, 2016, 45(5): 0518001. (in Chinese doi: 10.3788/irla201645.0518001 [19] 张营. 长波红外高光谱成像仪光学技术研究[D]. 上海: 中国科学院研究生院(上海技术物理研究所), 2016. ZHANG Y. Optical technology of long-wave infrared hyperspectral imaging[D]. Shanghai: Shanghai Institute of Technical Physics of the Chinese Academy of Sciences, 2016. (in Chinese). [20] 韩军, 李殉, 吴玲玲, 等. 一种光栅型成像光谱仪光学系统设计[J]. 应用光学,2012,33(2):233-239. HAN J, LI X, WU L L, et al. Optical system design of grating-based imaging spectrometer[J]. Journal of applied Optics, 2012, 33(2): 233-239. (in Chinese -