張緒國 尚志鳴 張躍東 曹桂麗

?

大氣色散對(duì)航空雙譜段高分辨率斜視成像影響

張緒國1,2尚志鳴1,2張躍東1,2曹桂麗1

（1 北京空間機(jī)電研究所，北京 100094） （2 北京市航空智能遙感裝備工程技術(shù)研究中心，北京 100094）

航空雙譜段高分辨率斜視成像載荷作為重要的成像手段，具有焦距長、分辨率高等特點(diǎn)，由于大氣層上疏下密的分層特性，光線在大氣中斜視傳播時(shí)，導(dǎo)致光譜展寬并產(chǎn)生色散和畸變，嚴(yán)重影響系統(tǒng)成像分辨率和目標(biāo)定位精度。文章分析了斜視成像的幾何特性；利用MODTRAN軟件仿真分析了斜視成像時(shí)大氣透過特性和光譜散射特性；通過光線追蹤法定量分析了大氣色散對(duì)可見光和中波紅外譜段的光線彎曲和光譜色散的影響，并提出了基于多光譜ZnS楔形窗口補(bǔ)償方法。結(jié)果表明，遠(yuǎn)距離斜視成像可見光近紅外0.50~0.95μm透過率高，散射??；中波3.70~4.80μm散射和路徑輻射效應(yīng)較小，光線斜視傳播時(shí)，光線彎曲角度基本一致；大氣色散對(duì)可見光譜段影響較大，航高20km斜視120km時(shí)，成像分辨率由0.375m退化至4.2m，大氣色散對(duì)中波紅外譜段影響較小；色散補(bǔ)償后，可見光譜段色散角度下降為原來的1/2，中波紅外譜段不受影響，在提升可見光成像品質(zhì)的同時(shí)，保持中波紅外譜段色散特性不變。

大氣色散 斜視成像 雙譜段 高分辨率 航空遙感

0 引言

航空高分辨率成像是高分對(duì)地觀測(cè)系統(tǒng)的重要組成部分，具有時(shí)效性強(qiáng)、機(jī)動(dòng)靈活等特點(diǎn)，在信息化中發(fā)揮著越來越大的作用。航空斜視高分辨率成像載荷因其焦距長、分辨率高、成像范圍廣等特點(diǎn)，在大縱深信息獲取等方面，體現(xiàn)出極其重要的意義[1-4]。國外航空斜視高分辨率成像始于20世紀(jì)60年代，經(jīng)歷了膠片相機(jī)向CCD相機(jī)、TDICCD相機(jī)、雙/多譜段相機(jī)的發(fā)展歷程[5-6]。比較著名的航空高分辨率斜視相機(jī)有DB-110、F-9120、LOROP、474/494HD等，隨著光學(xué)、成像器件、穩(wěn)定控制和圖像處理技術(shù)的發(fā)展，載荷焦距逐漸增加到5 500mm，100km斜視成像距離理論分辨率最高可達(dá)0.1m[7]。

航空斜視高分辨率成像工作在大氣環(huán)境下，由于大氣層的分層特性，不可避免受到大氣色散和湍流的影響，焦距越長，影響越明顯；大氣色散使光線展寬，大氣湍流效應(yīng)使光波相位發(fā)生畸變降低了光學(xué)成像系統(tǒng)分辨率，影響成像品質(zhì)。因此，分辨率的提升除對(duì)載荷穩(wěn)定控制精度的要求提高外，大氣色散、湍流等也成為制約航空斜視高分辨率像質(zhì)的重要因素[8-11]。早在膠片相機(jī)年代，國外就研究了大氣對(duì)航空高分辨率成像品質(zhì)的影響，20世紀(jì)80年代，美國研究了飛機(jī)氣流對(duì)成像的影響[11]，通過理論分析和建模仿真研究了大氣湍流對(duì)航空高分辨率相機(jī)成像品質(zhì)的影響。文獻(xiàn)[12-14]從大氣湍流引起的系統(tǒng)MTF退化理論、圖像復(fù)原和成像體制等方面，分析了大氣湍流引起的MTF退化及圖像復(fù)原方法，但是，針對(duì)大氣色散對(duì)成像品質(zhì)的研究較少。本文以雙譜段航空高分辨率斜視成像系統(tǒng)為例，利用光線追蹤法建立了不同波長與折射率的關(guān)系模型，重點(diǎn)研究了色散對(duì)高分辨率斜視成像角分辨率的影響，并提出了雙譜段色散補(bǔ)償方法。

1 斜視成像分析

1.1 幾何模型

航空遠(yuǎn)距離斜視成像，通常飛行高度在8~20km，獲取縱深達(dá)50~150km信息，斜視成像幾何模型如圖1所示。由圖1可知，地面收容寬度（即成像時(shí)刈幅圖像對(duì)應(yīng)的地面區(qū)域的寬度）與航高、相機(jī)拍攝俯角、地球半徑和張角之間的關(guān)系為

式中=6 400km，以航高10km和20km為例，不同俯角對(duì)應(yīng)的收容寬度如表1所示。

由表1可知，為了實(shí)現(xiàn)地面收容寬度50~150km，在航高10km時(shí)，俯角在5°~10°；航高20km時(shí)，俯角在10°~20°。由此可見，斜視成像通常在5°~20°俯角范圍內(nèi)，基本上相當(dāng)于平視通過大氣層，由于大氣層的分層結(jié)構(gòu)和密度的非均勻性，導(dǎo)致光線彎曲和色散，圖像品質(zhì)退化，造成的影響無法忽略。

圖1 航空遠(yuǎn)距離斜視成像示意

表1 不同俯角對(duì)應(yīng)的收容寬度

Tab.1 Coverage range with different depression angle

1.2 成像譜段分析

航空雙譜段高分辨率斜視成像載荷作為重要的裝備，需要覆蓋可見光近紅外和紅外譜段，以滿足全天時(shí)晝夜成像需求，因此，需要考慮大氣散射、吸收等光線傳輸?shù)挠绊?。?0km成像高度，斜視120km為例，通過MODTRAN對(duì)可見光近紅外和中長波紅外譜段進(jìn)行了透過特性和散射特性的分析。

如圖2所示，斜視成像模式下，可見光近紅外譜段透過率在0.50μm以下透過率低于0.1，在0.40μm時(shí)僅有0.03，在0.95~1.00μm有明顯的吸收峰。0.40~0.50μm譜段由于波長較短，太陽光散射和單次散射較大，散射的存在會(huì)導(dǎo)致圖像噪聲增加，圖像品質(zhì)退化。綜合考慮斜視成像的大氣透過率特性、散射、吸收特性等，可見光近紅外斜視成像譜段范圍選擇0.50~0.95μm。

圖2 可見光近紅外譜段吸收和散射特性

由圖3和圖4可知，中波3~5μm和長波8~12μm均有明顯的大氣窗口，透過率均值相當(dāng)，但是長波紅外路徑熱輻射和大氣散射影響比中波高一個(gè)數(shù)量級(jí)，路徑越長，長波受路徑輻射和太陽散射影響越嚴(yán)重。因此，綜合考慮中波紅外探測(cè)器的光譜響應(yīng)特性，對(duì)于斜視成像，中波3.70~4.80μm具有更好的效果。

圖3 中波紅外譜段透過和輻射特性

圖4 長波紅外譜段透過和輻射特性

2 色散影響分析

2.1 折射率變化分析

環(huán)繞地球的整個(gè)大氣層呈現(xiàn)上疏下密的狀態(tài)，且不同波長的光線折射率也有所不同，空氣折射率的大小可由Edlen公式求出[15]

式中為波數(shù)（cm–1）；為總的大氣壓（hPa）；H2O為水蒸氣對(duì)應(yīng)的大氣壓（hPa）；0為標(biāo)準(zhǔn)大氣壓，0=1 013.25hPa；為溫度（K）；0、1、2、1、2、1、2為常數(shù)，0=83.43，1=185.08，2=4.11，1=1.14×105，2=6.24×104，0=43.49，1=1.70×104。

地球表面大氣折射率變化曲線如圖5所示，由圖可知可見光近紅外0.50~0.95μm折射率變化相對(duì)差異遠(yuǎn)大于中波紅外3.70~4.80μm變化差異，即大氣色散對(duì)可見光近紅外譜段影響較嚴(yán)重。

圖5 地球表面大氣折射率曲線

在0~20km范圍內(nèi)，可見光和中波紅外譜段折射率隨高度變化，如圖6所示。由圖6放大部分可知，3.70μm和4.80μm折射率變化基本保持一致，而0.50μm和0.95μm折射率差別較大。

圖6 折射率隨高度變化曲線

2.2 折射率變化引起的光線彎曲

由于折射率隨高度的變化，大氣呈現(xiàn)分層結(jié)構(gòu)，斜視成像光線傳輸路徑示意如圖7所示。當(dāng)觀測(cè)目標(biāo)偏離天頂位置時(shí)，光線在斜視傳輸時(shí)會(huì)發(fā)生偏折，假設(shè)觀測(cè)點(diǎn)為，觀測(cè)角為′，理想情況下，光線沿傳播，由于大氣的分層結(jié)構(gòu)，光線傳播時(shí)會(huì)發(fā)生彎曲，實(shí)際傳播光線為。

式中為地球半徑；為地球半徑與航高的和（為地心）；′為入射到鏡頭的光線角度；為目標(biāo)光線的出射角。

對(duì)于10~30km航高，斜視成像120km，大氣折射率變化引起的光線彎曲角度變化如圖8所示。

圖7 斜視成像光線傳輸路徑示意

圖8 不同高度斜視成像大氣色散引起的光線彎曲

圖9 不同高度大氣色散角度變化

由圖8可知，由于大氣折射變化的影響，可見光光線彎曲角度大于中波紅外彎曲角度，且對(duì)于相同的成像距離，航高越低，大氣折射率變化對(duì)光線彎曲角度影響越嚴(yán)重；對(duì)于10km航高，彎曲角度最大為2.2mrad，對(duì)于120km的成像距離，理想情況下位置偏差達(dá)到260m以上，嚴(yán)重影響目標(biāo)實(shí)時(shí)定位精度。

2.3 折射率變化引起的色散

無論是可見光還是紅外，均具有一定的光譜范圍，不同波長光線大氣折射率不同，光線在大氣中傳播時(shí)，不同波長彎曲的角度不同，大氣色散使像面上的點(diǎn)像由一個(gè)點(diǎn)變成一條離散的光譜，能量分散，使本來清晰的圖像變得模糊，產(chǎn)生畸變和彌散，降低圖像品質(zhì)。對(duì)于可見光和中波紅外，在航高10~30km，成像距離120km，大氣色散角度變化如圖9所示。

由圖9可知，對(duì)于相同的成像距離，航高越低，色散角度越大，航高10km，斜視120km，可見光色散角度達(dá)到35μrad，以DB-110相機(jī)為例，可見光焦距2 794mm，斜視120km時(shí)，理論分辨率0.375m，由于色散影響，分辨率彌散為4.2m，分辨率退化為原來的1/11，因此必須采取色散補(bǔ)償措施，才能保證成像分辨率。

由圖8和圖9可知，中波紅外僅存在光線彎曲，色散角度很小，基本上不受大氣色散的影響；航高越高，大氣色散影響越小。

3 色散補(bǔ)償方法分析

通過對(duì)大氣折射率變化導(dǎo)致的光線彎曲和色散引起的譜段展寬分析可知，對(duì)于航空遠(yuǎn)距離斜視成像：1）大氣影響可等效成光楔；2）可見光和中波紅外譜段均存在光線彎曲，且彎曲角度接近；3）大氣色散對(duì)可見光影響較大，對(duì)中波紅外譜段影響較小。為了提高目標(biāo)成像分辨率，必須采取補(bǔ)償措施，色散補(bǔ)償需要同時(shí)滿足能夠減小可見光光譜色散，但不能對(duì)中波紅外色散產(chǎn)生影響的條件。

由于斜視成像大氣效應(yīng)類似于“光楔”，如圖10所示，因此，可以通過在成像系統(tǒng)前增加一個(gè)與大氣“光楔”相反的楔形窗口玻璃，補(bǔ)償大氣色散對(duì)可見光和中波紅外譜段的影響。

既透可見光，又透中波紅外的玻璃，有人造藍(lán)寶石、BGG玻璃、多光譜ZnS等。其中人造藍(lán)寶石和多光譜ZnS不同波長的折射率隨波長的變化曲線如圖11所示。

圖10 大氣色散及補(bǔ)償?shù)刃Ч麍D

圖11 人造藍(lán)寶石和多光譜ZnS變化曲線

由圖11可知，人造藍(lán)寶石在可見光和中波紅外譜段折射率變化基本上成線性變化，可以滿足可見光譜段的色散補(bǔ)償，但是，會(huì)導(dǎo)致中波紅外譜段產(chǎn)生色散。多光譜ZnS在可見光譜段折射率變化較大，在中波紅外譜段變化較小，可以滿足補(bǔ)償可見光色散，同時(shí)不會(huì)對(duì)中波紅外譜段色散產(chǎn)生影響，因此，多光譜ZnS可以滿足斜視成像色散補(bǔ)償要求。

由于采用固定楔形角度無法滿足所有傾斜角度，為了滿足更大覆蓋范圍的使用要求，由圖9可知，航高10km，可見光0.50~0.90μm色散角度約為35μrad；航高20km，可見光色散角度約為26μrad。由 式（3），可求出需要多光譜ZnS窗口的楔形角度分別為41″和23″，經(jīng)過色散補(bǔ)償后，可見光和中波紅外譜段色散角度對(duì)比如圖12所示。

由圖12可知，采用楔形多光譜ZnS窗口進(jìn)行色散補(bǔ)償對(duì)可見光具有較好的效果，同時(shí)，對(duì)中波紅外基本上沒有影響。航高10km時(shí)，經(jīng)過色散補(bǔ)償，色散角由35μrad，下降至12μrad。航高20km時(shí)，經(jīng)過色散補(bǔ)償，色散角由26μrad，下降至13μrad。

圖12 色散補(bǔ)償前后結(jié)果對(duì)比

4 結(jié)束語

隨著分辨率的提升，大氣逐漸成為影響航空高分辨率斜視成像品質(zhì)的主要因素。大氣折射率隨高度的變化，導(dǎo)致光線在斜視傳輸時(shí)產(chǎn)生彎曲和色散，進(jìn)而引起圖像品質(zhì)退化和畸變。針對(duì)大氣色散的影響，本文利用MODTRAN軟件和光線追蹤法，仿真分析了航高10km、20km，斜視120km，折射率變化引起的光線彎曲和色散角度，結(jié)果表明，1）大氣色散對(duì)航空遠(yuǎn)距離斜視成像系統(tǒng)的影響不可忽略，折射率變化引起的光線彎曲角度達(dá)到2.2mrad，直接定位誤差260m；2）色散導(dǎo)致的可見光譜段展寬角度最大約為35μrad，是理想角分辨率的1/11；3）色散對(duì)中波紅外基本沒有影響；4）對(duì)于相同的成像距離，航高越低，色散影響越嚴(yán)重；5）采用楔形多光譜ZnS進(jìn)行補(bǔ)償，在保證中波紅外不受影響的情況下，可見光色散角度降低為原來的1/2。

[1] Davis L, Mrinal I, Larry M, et al. The Goodrich 3rd Generation DB-110 System: Successful Flight Test on the F-16 Aircraft[C]//SPIE Vol.6546, Airborne Intelligence, Surveillance, Reconnaissance (ISR) Systems and Applications IV. SPIE, 2007.

[2] Davis L, William A, Mrinal I, et al. The Goodrich DB-110 System: Multi-band Operation Today and Tomorrow[C]//SPIE, Vol. 5109, Airborne Reconnaissance XXVII. SPIE, 2003.

[3] Petrushevsky V. Calibration Method for IR Channel of Dual Band Long Range Airborne Camera[C]//SPIE, Vol. 5406, Infrared Technology and Applications XXX. SPIE, 2004.

[4] Russo M J. An All-Beryllium-Aluminum Optical System for Reconnaissance Applications[C]//SPIE, Vol. 7425, Optical Materials and Structures Technologies IV. SPIE, 2009.

[5] Augustyn T. The KS-146A Long Range Oblique Photography (LOROP) Camera System[C]//SPIE, Vol. 309, Airborne Reconnaissance V. SPIE, 1981.

[6] 李永昆, 林招榮, 張緒國. 國外遠(yuǎn)距斜視航空相機(jī)發(fā)展概況[J]. 航天返回與遙感, 2017, 38(6): 11-18. LI Yongkun, LIN Zhaorong, ZHANG Xuguo. Development Survey of Foreign Aerial Cameras for Distant Oblique Reconnaissance[J]. Spacecraft Recovery & Remote Sensing, 2017, 38(6):11-18. (in Chinese)

[7] 龔惠興. 中國電子信息工程科技發(fā)展研究—遙感技術(shù)及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2018: 36-51. GONG Huixing. Research on the Development of Electronic Information Engineering Technology: Remote Sensing Technology and Application[M]. Beijing: Science Press, 2018: 36-51.

[8] Eismann M T, LeMaster D A. Aerosol Modulation Transfer Function Model for Passive Long-range Imaging Over a Nonuniform Atmospheric Path[J]. Optical Engineering, 2013, 52(4): 1-12.

[9] 張緒國, 曹桂麗, 張躍東, 等. 大氣湍流對(duì)航空遠(yuǎn)距離斜視成像系統(tǒng)成像質(zhì)量的影響[C]//光學(xué)遙感載荷與信息處理技術(shù). 北京: 中國宇航學(xué)會(huì)空間遙感專業(yè)委員會(huì), 2013. ZHANG Xuguo, CAO Guili, ZHANG Yuedong, et al. Influence of Atmospheric Turbulence on Aerial Standoff Imaging System[C]//Remote Sensing Payload and Information Processing. Beijing: Chinese Society of Astronautics, 2013. (in Chinese)

[10] LeMaster D A, Eismann M T. Impact of Atmospheric Aerosols on Long Range Image Quality[C]//SPIE, Vol. 8355, Infrared Imaging Systems: Design, Analysis, Modeling, and Testing XXIII. SPIE, 2012.

[11] Kopeika N S. Imaging Through the Atmosphere for Airborne Reconnaissance[J]. Optical Engineering, 1987, 26(11): 1146-1154.

[12] SU C M, RAO R Z, HUANG H H, et al. Image Resolution Measurement for Horizontal Propagation in the Atmosphere[C]// SPIE, Vol. 5832, Optical Technologies for Atmospheric, Ocean, and Environmental Studies. SPIE, 2005.

[13] ARBEL D, MOLDOVAN O, JACOBSON R, et al. Imaging Vertically Through the Atmosphere: Restoration of Satellite Images Based on Atmospheric MTF Evaluation[C]//SPIE, Vol.3433, Propagation and Imaging through the Atmosphere II. SPIE, 1998.

[14] ZHAO S M, WANG B, GONG L Y, et al. Improving the Atmosphere Turbulence Tolerance in Holographic Ghost Imaging System by Channel Coding[J]. Journal of Light Wave Technology, 2013, 31(17): 2823-2828.

[15] Mrozek F, Bennett R A. Aerial Reconnaissance Camera System with Atmospheric Dispersion Correction[P]: US 20120038769. 2012.

[16] 吉書鵬, 李同海. 一種光機(jī)掃描型機(jī)載廣域偵察監(jiān)視系統(tǒng)設(shè)計(jì)[J]. 紅外技術(shù), 2018, 40(1): 20-26. JI Shupeng, LI Tonghai. A Design of Wide Area Airborne Surveillance System Based on Optical and Mechine Scanning[J]. Infrared Technology, 2018, 40(1): 20-26. (in Chinese)

[17] 項(xiàng)劍鋒, 李保霖, 劉明, 等. 軍用航空相機(jī)發(fā)展分析與展望[J]. 學(xué)術(shù)研究, 2018, 288(12): 222-225. XIANG Jianfeng, LI Baolin, LIU Ming, et al. Analysis and Overview of the Development About Military Aerial Camera[J]. China Science & Technology Overview, 2018, 288(12): 222-225. (in Chinese)

[18] 張春曉, 林招榮, 姚毅剛, 等. 被動(dòng)成像廣域空中監(jiān)視系統(tǒng)綜述[J]. 航天返回與遙感, 2014, 35(1):37-45. ZHANG Chunxiao, LIN Zhaorong, YAO Yigang, et al. Overview of Passive Imaging Wide Area Airborne Surveillance[J]. Space Recovery & Remote Sensing, 2014, 35(1): 37-45. (in Chinese)

[19] 張緒國, 張躍東, 吳春楠, 等. 振動(dòng)對(duì)高分辨率航空相機(jī)成像質(zhì)量的影響[C]//中國宇航學(xué)會(huì), 第19屆中國遙感大會(huì). 西安, 2014. ZHANG Xuguo, ZHANG Yuedong, WU Chunnan, et al. Influences of Vibrations on High Resolution Aerial Camera[C]//19thConference on Remote Sensing of China, Xian, 2014. (in Chinese)

[20] 喬川, 丁亞林, 許永森, 等. 大角度傾斜成像航空相機(jī)對(duì)地目標(biāo)定位[J]. 光學(xué)精密工程, 2017, 25(7): 1714-1722. QIAO Chuan, DING Yalin, XU Yongsen, et al. Ground Target Geo-location Using Imaging Aerial Camera with Large Inclined Angles[J]. Optics and Precision Engineering, 2017, 25(7): 1714-1722. (in Chinese)

Influence of Atmospheric Chromatic Dispersion on Aerial Dual-band High Resolution Standoff Imaging

ZHANG Xuguo1,2SHANG Zhiming1,2ZHANG Yuedong1,2CAO Guili1

（1 Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China） （2 Beijing Engineering Technology Research Center of Aerial Intelligent Remote Sensing Equipments, Beijing 100094, China）

As an important product in imaging, the aerial dual-band high resolution standoff imaging payload has the characteristics of long focus length and high resolution. However, spectra are broadened with the layer effects of atmospheric density variation when light transmitted obliquely, which results in chromatic dispersion, distortion, resolution degradation and targeting error. In this paper, geometric characteristics of standoff imaging are analyzed. Characteristics of transmittance and scattering effects are simulated with MODTRAN. Ray tracing method is adopted to quantitatively analyze the influences of light ray deflection and chromatic dispersion. A dispersion compensation method is proposed based on multispectral ZnS window wedge. The results show as followed: 1) visible near infrared (VNIR) 0.5~0.95μm has good transmittance and low scattering, and middle wave infrared (MWIR) 3.7~4.8μm has lower scattering and path radiation effects; 2) visible and middle wave infrared spectral bands have similar deflection angle; 3) visible spectral band is influenced more seriously by dispersion than middle wave infrared, and the ground resolution degraded from 0.375m to 4.2m for altitude 20km standoff range 120km imaging; 4) after adopting compensation method, dispersion angle is greatly reduced in visible band, with about half of the original value, and has no effect in middle wave infrared band.

atmospheric chromatic dispersion; standoff imaging; dual-band; high resolution; aerial remote sensing

TP701

A

1009-8518(2019)03-0073-09

10.3969/j.issn.1009-8518.2019.03.010

張緒國，男，1982年生，2011年獲北京航空航天大學(xué)物理電子學(xué)專業(yè)博士學(xué)位，高級(jí)工程師。研究方向?yàn)楹娇崭叻直媛市币暢上駛刹毂O(jiān)視技術(shù)。E-mail：zxg_ltt@163.com。

（編輯：王麗霞）

2018-09-05

國家重大科技專項(xiàng)工程