李軼庭,王靈杰,張玉慧,劉銘鑫
(1.長春理工大學 光電測控與光信息傳輸技術教育部重點實驗室,吉林 長春 130022;2.中國科學院 長春光學精密機械與物理研究所 中國科學院光學系統(tǒng)先進制造技術重點實驗室,吉林 長春 130033)
地球靜止軌道衛(wèi)星成像光譜儀在氣象觀測、資源調查、環(huán)境監(jiān)測等領域具有重大的應用價值。而覆蓋可見光到長波紅外的光譜儀器具有晝夜觀測的技術優(yōu)勢[1]。地球靜止軌道凝視型面陣成像光譜儀對其時間分辨率、空間分辨率和成像幅寬、遙感成像動態(tài)范圍提出了更高的要求,以期獲得更加豐富的目標層次及細節(jié)信息,以便對探測目標信息進行有效地獲取及識別,由于衛(wèi)星資源有限,大數(shù)據(jù)量會產生存儲、傳輸及處理的困難。壓縮感知技術以其先壓縮再采樣的特點,在現(xiàn)有信號傳輸水平下可以有效提高獲取圖像的效率,降低數(shù)據(jù)采集的壓力。在此背景下,利用地球靜止軌道平臺可以長期駐留固定區(qū)域上空的特點[2-5],結合壓縮感知光譜成像技術,對寬譜段成像光學系統(tǒng)進行分析,并對其光學系統(tǒng)進行設計。
將壓縮感知理論應用于光譜成像技術中,可以實現(xiàn)對光譜圖像信息的壓縮采集,解決了衛(wèi)星數(shù)據(jù)存儲、傳輸、處理的問題,其原理如圖1所示。待測光經(jīng)過物鏡成像于離散編碼模板上,離散編碼模板將得到的二維信息調制準直后,由色散元件色散為光譜帶,會聚后,探測器獲取經(jīng)過壓縮后的二維離散光譜混疊信息,通過優(yōu)化解混后即可高概率獲得原始光譜信息,壓縮感知成像光譜技術獲得的采樣數(shù)據(jù)量遠小于傳統(tǒng)成像光譜技術的數(shù)據(jù)量,可將將衛(wèi)星信號數(shù)據(jù)量降低兩個數(shù)量級[6-8]。
圖1 壓縮感知光譜成像系統(tǒng)原理圖Fig.1 Schematic diagram of compressive sensing spectral imaging system
根據(jù)此模型提出本文的光學系統(tǒng)結構形式,根據(jù)光學系統(tǒng)全譜段的要求,需將系統(tǒng)分為可見光、中波紅外、長波紅外3個譜段光路分別進行設計[7],采用分色片分譜段,選擇DMD作為離散編碼模板,可見光所選DMD像元大小為7.6 μm,采用一個DMD像元與4個探測器像元大小進行匹配,中波紅外光路所選DMD大小為7.6 μm,采用一個DMD像元與一個探測器像元大小匹配,長波紅外光路所選DMD像元大小為10.8 μm,采用一個DMD像元與1個探測器像元大小進行匹配,DMD的開合角度為±12°,通過中繼鏡組的放大倍率匹配編碼模板與探測器像元面積的不同。目前常用DMD的像元數(shù)受芯片工藝限制,無法與探測器實現(xiàn)完美的匹配,當探測器分辨率高于編碼DMD的分辨率時,可通過對DMD像元進行灰度分級的辦法,保證分級后的編碼像元與探測器像元一一對應,即可保證系統(tǒng)的整體分辨率。由于DMD帶來的額外像差可由像面傾斜進行補償,基于DMD的目標模擬器應用于長波紅外波段時,照明光束需要選取合適的入射角以減少DMD的衍射效應對系統(tǒng)成像對比度造成的影響[9-11]。整體光學系統(tǒng)的模型如圖2(彩圖見期刊電子版)所示,其由物鏡、分光片、一次成像系統(tǒng)、準直系統(tǒng)、色散元件、DMD,以及光譜信息采集處理系統(tǒng)7部分組成,待測光經(jīng)過物鏡入射到分光系統(tǒng)進行波段選擇,經(jīng)成像系統(tǒng)成像于編碼模板上,由準直系統(tǒng)準直為平行光,在色散元件表面色散為光譜帶,經(jīng)會聚系統(tǒng)會聚在探測器上。光學系統(tǒng)覆蓋可見光、中波紅外、長波紅外,譜段較寬。根據(jù)壓縮感知成像光譜系統(tǒng)模型,需對系統(tǒng)進行模塊化設計、可降低整體設計的難度,在保證各模塊的成像質量后,將各模塊整合最終對整體進行設計[12-20]。
圖2 基于壓縮感知的大口徑多譜段光譜儀物理模型Fig.2 Physical model of visual and infrared large aperture multispectral sensor based on compressive sensing
本文光學系統(tǒng)包括3個譜段,各譜段設計要求如表1所示。
表1 光學系統(tǒng)設計要求Tab.1 Requirements for optical system design
根據(jù)空間分辨率以及幅寬的要求,所選探測器參數(shù)如表2所示。
表2 光學系統(tǒng)探測器參數(shù)Tab.2 Parameters of the optical system’s detector
根據(jù)式(1)對系統(tǒng)、焦距及視場角進行計算,
其中H為衛(wèi)星的軌道高度,B為相機的幅寬,L為像面的長度,a為像元尺寸,ω為半視場角,f為系統(tǒng)焦距。
相機鏡頭孔徑D的大小影響相機的分辨率、調制傳遞函數(shù)(MTF)、曝光量、信噪比和相機尺寸。D的確定需考慮探測器的信噪比、噪聲等效溫差、噪聲等效反射率差的要求。綜合考慮后,光學系統(tǒng)孔徑設定為700 mm,見表3。
表3 光學系統(tǒng)最終設計參數(shù)Tab.3 Parameters of the designed optical system
物鏡設計為同軸反射無焦系統(tǒng),可以實現(xiàn)無色差設計以及衛(wèi)星長時間凝視成像的軌道的震動補償,其放大倍率為5×,要求系統(tǒng)遮攔比小于1/3。同時根據(jù)DMD像元大小與探測器像元大小計算不同光路準直光路的焦距,準直光路焦距與會聚光路焦距保持一致。選擇光柵作為系統(tǒng)的色散元件,光柵的刻痕數(shù)直接決定系統(tǒng)的光譜分辨率。根據(jù)式(2)、式(3)與各譜段會聚光路焦距計算光柵刻線數(shù),如表4所示。
表4 子系統(tǒng)參數(shù)Tab.4 Parameters of the sub-optical system
其中δ 為像元大小,δ λ為系統(tǒng)的光譜分辨率。
設計結果要求一次會聚系統(tǒng)、準直系統(tǒng)、會聚系統(tǒng)各部分像質良好,對于高分辨率的成像光譜儀系統(tǒng),二級光譜的存在將嚴重影響最終的成像質量,二級光譜的色差大小與系統(tǒng)焦距成正比,因此對于長焦距的可見光路系統(tǒng),需對其二級光譜進行優(yōu)化及評價。由于系統(tǒng)空間維很大,必然存在譜線彎曲和色畸變,譜線彎曲是指不同波長彎曲圖像與直線的偏離程度,色畸變是由光譜成像系統(tǒng)對不同波長的像放大率不同造成的,嚴重時,會影響成像質量,譜線彎曲過大還會降低系統(tǒng)的空間分辨率,為獲得準確的多光譜圖像,譜線彎曲以及色畸變需控制在1/2個像元以內。在中波、長波光路中選擇制冷型紅外探測器,因此需對冷光闌進行匹配。
基于上述參數(shù)分別對3個譜段進行光學系統(tǒng)設計,其整體結構如圖3所示。
圖3 光學系統(tǒng)結構圖Fig.3 Structure of the optical system
物鏡選擇以卡塞格林為基礎的同軸三反無焦系統(tǒng),在二鏡、三鏡中加入反射鏡以折轉光路,使系統(tǒng)緊湊,同時在主鏡與三鏡中引入高階非球面在保證系統(tǒng)放大率和遮攔比的情況下,像質達到最佳。其放大倍率為5×,遮攔比為0.185。
對可見光、中波紅外、長波紅外系統(tǒng)整體要求結構緊湊,對各光路不同系統(tǒng)進行模塊化設計[15-20],可以降低系統(tǒng)的設計難度。
可見光路優(yōu)化中,約束一次會聚光路、準直光路、會聚光路焦距分別為900、534、534 mm,約束各鏡片厚度以及系統(tǒng)整體長度,對于長焦距系統(tǒng),對二級光譜的優(yōu)化選擇阿貝數(shù)相差較大、部分色散相同或相近的玻璃對其進行優(yōu)化。采用棱鏡-光柵結構作為可見光路的分光元件,補償系統(tǒng)的譜線彎曲,通過真實光線的追跡控制系統(tǒng)的色畸變大小,利用像面傾斜的方式校正軸向色差并達到平譜面的目的,在對各模塊優(yōu)化后,對整體結構進行優(yōu)化,優(yōu)化結果可見光會聚光路在波長分別為500、700、900 nm時,各視場MTF在78.125 lp/mm處分別高于0.674、0.562、0.455,點列斑83%以上位于艾里斑內,彌散斑半徑最大分別為5.231、4.694、4.113 μm,小于一個探測器像元尺寸。在不同視場下光譜變化下20 nm可以清晰分辨,滿足20 nm的光譜分辨率,其中不同顏色光斑分別代表900、880、720、700、520、500 nm波長在光學系統(tǒng)的點列圖大小,可以看出光學系統(tǒng)在各視場均滿足20nm的光譜分辨率,可見光系統(tǒng)的二級光譜小于0.05mm,譜線彎曲很小,色畸變最大值為2.753 μm,小于1/2個像元大小,如圖4~圖8(彩圖見期刊電子版)所示。
圖4 可見光一次會聚光路Fig.4 First converging lens in visual part
圖5 可見光準直及會聚光路Fig.5 Collimating and converging lens in visual part
圖8 可見光路二級光譜Fig.8 Secondary spectrum in visual part
圖6 波長為 (a) 500 nm; (b) 700 nm; (c) 900 nm時的MTF值Fig.6 MTF at (a) 500 nm; (b) 700 nm; (c) 900 nm
圖7 波長為(a) 500 nm; (b) 700 nm及(c) 900 nm時點列圖Fig.7 Spot diagram at (a) 500 nm; (b) 700 nm and (c) 900 nm
中波紅外的優(yōu)化中,約束一次會聚光路、準直光路、會聚光路焦距分別為270、137、137 mm,約束各鏡片厚度以及系統(tǒng)整體長度,并對部分表面引入二次曲面以及高階非球面,針對中波紅外系統(tǒng)的譜線彎曲以及色畸變,采用離軸透鏡的方式對其進行補償,先優(yōu)化得到一次成像未離軸時的光學系統(tǒng),然后根據(jù)系統(tǒng)存在的譜線彎曲量,確定透鏡的離軸量,最后適當改變透鏡的半徑完成譜線彎曲的校正,通過真實光線的追跡控制系統(tǒng)的色畸變大小,采用二次成像的方式對光學系統(tǒng)的冷光闌進行匹配,控制系統(tǒng)的焦距以及各視場的邊緣光線出瞳像高對系統(tǒng)進行優(yōu)化。優(yōu)化后各模塊成像質量良好、最終得到的中波紅外光學系統(tǒng)在波長3500 nm、3800 nm、4100 nm處的MTF在33.3 lp/mm時分別高于0.607、0.566、0.518。彌散斑半徑最大分別為5.419、5.922、6.678 μm。在不同視場下光譜變化50 nm時各波長點列圖光斑可以清晰的分辨,其中不同顏色光斑分別代表4100、4050、3800、3750、3550、3500 nm波長光學在光學系統(tǒng)的點列圖大小,可以看出光學系統(tǒng)滿足50 nm的光譜分辨率,譜線彎曲很小,色畸變最大值為5.48 μm,小于1/2個像元,如圖9~圖14(彩圖見期刊電子版)所示。
圖9 不同視場下光譜分辨率Fig.9 Spectral resolution in different fields of view
圖11 中波紅外一次會聚光路Fig.11 First converging lens of MWIR part
圖12 中波紅外準直及會聚光路Fig.12 Collimating and converging lens of MWIR part
圖13 波長為(a) 3500 nm;(b) 3800 nm及(c) 4100 nm時的MTFFig.13 MTF at (a) 3500 nm; (b) 3800 nm; (c) 4100 nm
圖14 波長為(a) 3500 nm; (b) 3800 nm;及(c) 4100 nm點列圖Fig.14 Spot diagram at (a) 3500 nm; (b) 3800 nm; (c) 4100 nm
長波紅外的優(yōu)化中,約束一次會聚光路、準直光路、會聚光路焦距分別為275、126、126 mm,約束各鏡片厚度以及系統(tǒng)整體長度并對部分表面引入二次曲面以及高階非球面,針對長波紅外系統(tǒng)的譜線彎曲以及色畸變,與中波紅外系統(tǒng)相似,采用離軸透鏡的方式對其優(yōu)化,通過真實光線的追跡控制系統(tǒng)的色畸變的大小,采用二次成像的方式對光學系統(tǒng)的冷光闌進行匹配,約束系統(tǒng)焦距以及各視場光線出瞳像高對系統(tǒng)進行優(yōu)化。優(yōu)化后各模塊成像質量良好、最終得到的長波紅外光學系統(tǒng)在波長7700、8600、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、9350、8600、8450、7850、7700nm波長光學在光學系統(tǒng)各視場的點列圖大小,可以看出光學系統(tǒng)滿足150 nm的光譜分辨率,譜線彎曲很小,色畸變最大值為9.35 μm,小于1/2個像元大小,如圖9~圖10、圖15~圖18(彩圖見期刊電子版)所示。
圖10 系統(tǒng)譜線彎曲示意圖Fig.10 Schematic diagram of spectral smile
圖15 長波紅外一次會聚光路Fig.15 First converging lens of LWIR part
圖16 長波紅外準直及會聚光路Fig.16 Collimating and converging lens of LWIR part
圖17 (a) 7700 nm;(b) 8600 nm; (c) 9500 nm處MTFFig.17 MTF at (a) 7700 nm; (b) 8600 nm; (c) 9500 nm
圖18 波長為(a) 7700 nm; (b) 8600 nm; (c) 9500 nm時的點列圖Fig.18 Spot diagram at (a) 7700 nm; (b) 8600 nm; (c) 9500 nm
為了獲取地球靜止軌道的高動態(tài)光譜數(shù)據(jù),設計了一種基于壓縮感知的全譜段成像光學系統(tǒng)。該系統(tǒng)通過物鏡無焦化設計以滿足衛(wèi)星軌道震動補償?shù)男枨蟛⑹褂梅稚瑢⒐鈱W系統(tǒng)譜段細化,以實現(xiàn)全譜段的探測,經(jīng)過優(yōu)化設計,系統(tǒng)結構緊湊,點列圖能量集中于艾里斑,各部分MTF均接近衍射極限,像質良好,光譜分辨率分別為可見光20 nm、中波紅外50 nm、長波紅外150 nm,滿足設計要求。