留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

AOTF后楔角对光谱图像像质的影响

徐泽夫 赵慧洁 贾国瑞

徐泽夫, 赵慧洁, 贾国瑞. AOTF后楔角对光谱图像像质的影响[J]. 红外与激光工程, 2022, 51(7): 20210590. doi: 10.3788/IRLA20210590
引用本文: 徐泽夫, 赵慧洁, 贾国瑞. AOTF后楔角对光谱图像像质的影响[J]. 红外与激光工程, 2022, 51(7): 20210590. doi: 10.3788/IRLA20210590
Xu Zefu, Zhao Huijie, Jia Guorui. Influence of the AOTF rear cut angle on spectral image quality[J]. Infrared and Laser Engineering, 2022, 51(7): 20210590. doi: 10.3788/IRLA20210590
Citation: Xu Zefu, Zhao Huijie, Jia Guorui. Influence of the AOTF rear cut angle on spectral image quality[J]. Infrared and Laser Engineering, 2022, 51(7): 20210590. doi: 10.3788/IRLA20210590

AOTF后楔角对光谱图像像质的影响

doi: 10.3788/IRLA20210590
基金项目: 国家重点研发计划(2016YFB0500505);国防基础科研计划(JCKY2019601C101);国家自然科学基金(61227806)
详细信息
    作者简介:

    徐泽夫,男,博士生,主要从事AOTF器件建模以及AOTF成像光谱仪系统集成方面的研究

    赵慧洁,女,教授,博士生导师,博士,主要从事高光谱遥感、三维测量等方面的研究

    通讯作者: 贾国瑞,男,副教授,博士生导师,博士,主要从事高光谱成像建模、仿真和数据预处理等方面的研究。
  • 中图分类号: O439;TP731

Influence of the AOTF rear cut angle on spectral image quality

Funds: National Key Research and Development Program of China(2016YFB0500505);National Defense Basic Scientific Research Program of China(JCKY2019601C101);National Natural Science Foundation of China (61227806)
  • 摘要: 针对AOTF成像光谱仪的光谱图像退化程度难以定量化评价的问题,提出了包含后楔角影响的AOTF器件调制传递函数(MTF)定量计算方法。该方法通过建立AOTF的光谱-空间维响应模型获取AOTF器件的线扩散函数(LSF),然后通过线扩散函数傅里叶变换得到理论MTF。在验证实验中,截止频率内的MTF理论值与实测值相对偏差小于15%,可以为AOTF的像质评估提供理论依据。进而,通过该定量计算方法讨论不同AOTF后楔角对光谱图像像质的影响。通过仿真分析得出,AOTF的后楔角难以同时满足横向色差校正与图像退化抑制,需要根据实际应用对于空间分辨率以及横向色差的要求进行取舍。所以,针对AOTF不同后楔角进行光谱图像像质的定量分析对于评估AOTF光谱图像像质具有重要意义,为AOTF器件设计提供了理论依据。
  • 图  1  动量匹配的波矢量布局图

    Figure  1.  Wave vector layout of momentum matching

    图  2  非理想声束的声光互作用

    Figure  2.  Acousto-optic interaction of non ideal sound beam

    图  3  AOTF声光作用以及后楔角的色散

    Figure  3.  Acousto-optic interaction of AOTF and dispersion of rear cut angle

    图  4  (a) 59.6 MHz (中心波长750 nm)的光谱-空间维响应(红线)、光谱维响应(蓝线)、空间维响应(绿线);(b)MTF曲线,蓝线为59.6 MHz空间维响应的傅里叶变换,红线为斜刀口靶实测MTF;(c) MTF实测值与理论值的相对偏差@59.6 MHz;(d) 59.6 MHz的光谱图像; (e) 99.4 MHz (中心波长500 nm)的光谱-空间维响应(红线)、光谱维响应(蓝线)、空间维响应(绿线);(f) MTF曲线,蓝线为99.4 MHz空间维响应的傅里叶变换,红线为斜刀口靶实测MTF; (g) MTF实测值与理论值的相对偏差@99.4 MHz;(h) 99.4 MHz的光谱图像

    Figure  4.  (a) Response of spectrum-space dimensional (red line), spectral response (blue line), spatial response (green line) at 59.6 MHz (central wavelength 750 nm); (b) MTF curve, blue line is Fourier transform of 59.6 MHz spatial response, red line is measured MTF by oblique knife edge target; (c) Deviation between the measured MTF and the theoretical MTF @59.6 MHz; (d) Spectral image of 59.6 MHz; (e) Response of spectrum-space dimensional (red line), spectral response (blue line), spatial response (green line) at 99.4 MHz (central wavelength 500 nm); (f) MTF curve, blue line is Fourier transform of 99.4 MHz spatial response, red line is measured MTF by oblique knife edge target; (g) Deviation between the measured MTF and the theoretical MTF @99.4 MHz; (h) Spectral image of 99.4 MHz

    图  5  (a) 59.6 MHz (中心波长750 nm)不同后楔角的MTF; (b) 99.4 MHz(中心波长500 nm)不同后楔角的MTF

    Figure  5.  (a) MTF of 59.6 MHz (central wavelength 750 nm) with different rear cut angles; (b) MTF of 99.4 MHz (central wavelength 500 nm) with different rear cut angles

    表  1  99.4 MHz (中心波长500 nm)和59.6 MHz (中心波长750 nm)不同后楔角的MTF30空间频率

    Table  1.   MTF30 with different rear cut angles at 99.4 MHz (center wavelength 500 nm) and 59.6 MHz (center wavelength 750 nm)

    99.4 MHz59.6 MHz
    Rear cut angle/(°)Spatial frequency/cycles·mm−1 @MTF30Rear cut angle/(°)Spatial frequency/cycles·mm−1 @ MTF30
    8 75.9 8 55.7
    6 45.6 6 50.9
    4 32.6 4 47.1
    2 25.6 2 43.4
    0 21.0 0 40.2
    −2 17.7 −2 37.4
    −4 15.3 −4 34.7
    −4.8 14.5 −4.8 33.7
    −6 11.7 −6 32.2
    −8 7.2 −8 29.9
    下载: 导出CSV
  • [1] 唐绍凡, 鲁之君, 王伟刚, 李欢. 航天高光谱成像仪简述(特邀)[J]. 红外与激光工程, 48(3): 0303003-0303003(9), (2019) doi:  10.3788/IRLA201948.0303003

    Tang Shaofan, Lu Zhijun, Wang Weigang, et al. Brief description of space hyperspectral imager (invited) [J]. Infrared and Laser Engineering, 2019, 48(3): 0303003. (in Chinese) doi:  10.3788/IRLA201948.0303003
    [2] 张瑞, 解琨阳, 景宁, 王志斌, 陈媛媛, 李克武, 乔亚. 基于超消色差1/4波片和AOTF的高光谱全偏振成像技术[J]. 红外与激光工程, 48(8): 814004-0814004(8), (2019) doi:  10.3788/IRLA201948.0814004

    Zhang Rui, Xie Kunyang, Jing Ning, et al. Technology of hyperspectral full polarization imaging based on superachromatic 1/4 wave plate and acousto-optic tunable filter [J]. Infrared and Laser Engineering, 2019, 48(8): 0814004. (in Chinese) doi:  10.3788/IRLA201948.0814004
    [3] Batshev V, Machikhin A, Martynov G. et al. Polarizer-free AOTF-based SWIR hyperspectral imaging for biomedical applications [J]. Sensors, 2020, 20(16): 4439. doi:  10.3390/s20164439
    [4] Korablev O, Ivanov A, Fedorova A, et al. Development of a mast or robotic arm-mounted infrared AOTF spectrometer for surface Moon and Mars probes[C]//Proceedings of SPIE, 2015, 9608: 960807.
    [5] Wu Changkun, Zhang Wei, Hao Yazhe. Design of a control system for a visible/near-infrared real-time imaging spectrometer [J]. Chinese Optics, 2022, 15(2): 348-354. (in Chinese) doi:  10.37188/CO.2021-0119
    [6] Chang I C. Noncollinear acousto-optic filter with large angular aperture [J]. Applied Physics Letters, 1974, 25(7): 370-372. doi:  10.1063/1.1655512
    [7] Yano T, Watanabe A. Acousto-optic TeO2 tunable filter using far-off-axis anisotropic bragg diffraction [J]. Applied Optics, 1976, 15(9): 2250-2258. doi:  10.1364/AO.15.002250
    [8] Gass P A, Sambles J R. Accurate design of a noncollinear acousto-optic tunable filter [J]. Optics Letters, 1991, 16(6): 429-431. doi:  10.1364/OL.16.000429
    [9] Voloshinov V B, Molchanov V Y, Mosquera J C. Spectral and polarization analysis of optical images by means of acousto-optics [J]. Optics & Laser Technology, 1996, 28(2): 119-127.
    [10] Voloshinov V B, Mosquera J C. Wide-aperture acousto-optic interaction in birefringent crystals [J]. Optics & Spectroscopy, 2006, 101(4): 635-641.
    [11] Krauz L, Páta P, Bednář J, et al. Quasi-collinear IR AOTF based on mercurous halide single crystals for spatio-spectral hyperspectral imaging [J]. Optics Express, 2021, 29(9): 12813-12832. doi:  10.1364/OE.420571
    [12] Gorevoy A V, Machikhin A S, Martynov G N, et al. Spatio-spectral transformation of noncollimated light beams diffracted by ultrasound in birefringent crystals [J]. Photonics Research, 2021, 9(5): 687-693. doi:  10.1364/PRJ.417992
    [13] Mantsevich S N, Kostyleva E I. Examination of the temperature influence on phase matching frequency in tunable acousto-optic filters [J]. Ultrasonics, 2018, 91: 45-51.
    [14] 常凌颖, 宋晶晶, 靳梦竹. AOTF成像光谱仪声光晶体光谱传递函数的研究[J]. 应用光学, 040(002): 202-209. (2019)

    Chang Lingying, Song Jingjing, Jin Mengzhu. Study on spectral transfer function based on dual AOTF acousto-optic crystal [J]. Journal of Applied Optics, 2019, 40(2): 202-209. (in Chinese)
    [15] Suhre D R, Gupta N. Acousto-optic tunable filter sidelobe analysis and reduction using telecentric confocal optics [J]. Applied Optics, 2005, 44(27): 5797-5801. doi:  10.1364/AO.44.005797
    [16] Wachman E S, Niu W, Farkas D L. Imaging acousto-optic tunable filter with 0.35-micrometer spatial resolution [J]. Applied Optics, 1996, 35(25): 5220-5226. doi:  10.1364/AO.35.005220
    [17] Denes L J, Gottlieb M S, Kaminsky B. Acousto-optic tunable filters in imaging applications [J]. Optical Engineering, 1998, 37(4): 1262-1267. doi:  10.1117/1.601962
    [18] Zhang C, Wang H, Qiu Y. Sidelobe reduction of Acousto-optic tunable filter by double-filtering[C]//Proceedings of SPIE, 2010, 7656: 76566F.
    [19] Anikin S P, Esipov V F, Molchanov V Y, et al. An acousto-optical imaging spectrometer for astrophysical measurements [J]. Optics & Spectroscopy, 2016, 121(1): 115-122.
    [20] Machikhin A S, Pozhar V E. Image transformation caused by wide-angle acousto-optic interaction [J]. Quantum Electronics, 2010, 40(9): 837-841. doi:  10.1070/QE2010v040n09ABEH014371
    [21] Voloshinov V B, Yushkov K B, Linde B. Improvement in performance of a TeO2 acousto-optic imaging spectrometer [J]. Journal of Optics A Pure and Applied Optics, 2007, 9(4): 341-347. doi:  10.1088/1464-4258/9/4/006
    [22] 俞宽新, 丁晓红, 庞兆广. 声光原理与声光器件[M]. 北京: 科学出版社, (2011).

    Yu Kuanxin, Ding Xiaohong, Pang Zhaoguagn. Acousto-optic Principle and Acousto-optic Devices[M]. Beijing: Science Press, 2011. (in Chinese)
  • [1] 李文璐, 裴丽, 白冰, 左晓燕, 王建帅, 郑晶晶, 李晶, 宁提纲.  集成光子储备池的可调谐光滤波器 . 红外与激光工程, 2023, 52(9): 20220915-1-20220915-11. doi: 10.3788/IRLA20220915
    [2] 余凯, 郭琦, 李娜, 成驰, 赵慧洁.  中红外AOTF成像仪光谱校准方法 . 红外与激光工程, 2023, 52(12): 20230291-1-20230291-9. doi: 10.3788/IRLA20230291
    [3] 赵雨时, 贺文俊, 刘智颖, 付跃刚.  编码孔径光谱成像仪中凸面闪耀光栅的研制 . 红外与激光工程, 2022, 51(3): 20220007-1-20220007-10. doi: 10.3788/IRLA20220007
    [4] 李居尚, 战荫泽, 张立东.  组合沃拉斯顿棱镜组的目标偏振检测方法 . 红外与激光工程, 2021, 50(11): 20210339-1-20210339-5. doi: 10.3788/IRLA20210339
    [5] 黄见, 时东锋, 孟文文, 查林彬, 孙宇松, 苑克娥, 胡顺星, 王英俭.  光谱编码计算关联成像技术研究 . 红外与激光工程, 2021, 50(1): 20200120-1-20200120-8. doi: 10.3788/IRLA20200120
    [6] 赵雨时, 贺文俊, 刘智颖, 付跃刚.  光谱维编码中红外光谱成像系统的光学设计 . 红外与激光工程, 2021, 50(12): 20210700-1-20210700-9. doi: 10.3788/IRLA20210700
    [7] 李美萱, 张斯淇, 李宏, 李楠, 任雨轩, 田嘉龙.  基于单次曝光多光谱关联成像系统带通滤波器的研究 . 红外与激光工程, 2020, 49(9): 20200169-1-20200169-7. doi: 10.3788/IRLA20200169
    [8] 唐钊, 张钧翔, 付士杰, 白晓磊, 盛泉, 史伟, 姚建铨.  基于MMI滤波器的可调谐连续光全光纤OPO . 红外与激光工程, 2019, 48(5): 520002-0520002(6). doi: 10.3788/IRLA201948.0520002
    [9] 张晨, 刘舒扬, 赵安娜, 王天鹤, 贾晓东.  基于小型化便携式光谱成像技术的分类识别系统研究 . 红外与激光工程, 2019, 48(10): 1023001-1023001(6). doi: 10.3788/IRLA201948.1023001
    [10] 刘智颖, 高柳絮, 黄蕴涵.  offner型连续变焦中波红外光谱成像系统设计 . 红外与激光工程, 2019, 48(7): 718003-0718003(9). doi: 10.3788/IRLA201948.0718003
    [11] 张瑞, 解琨阳, 景宁, 王志斌, 陈媛媛, 李克武, 乔亚.  基于超消色差1/4波片和AOTF的高光谱全偏振成像技术 . 红外与激光工程, 2019, 48(8): 814004-0814004(8). doi: 10.3788/IRLA201948.0814004
    [12] 陈娇, 童峥嵘, 张卫华, 薛力芳.  采用复合滤波器的温度可调谐多波长光纤激光器 . 红外与激光工程, 2018, 47(1): 105001-0105001(7). doi: 10.3788/IRLA201847.0105001
    [13] 贺超, 廖同庆, 吴昇, 魏小龙.  可调谐窄带宽的负系数微波光子滤波器 . 红外与激光工程, 2017, 46(9): 920003-0920003(6). doi: 10.3788/IRLA201746.0920003
    [14] 毛红敏, 马锡英, 王晓丹, 徐国定.  多通道可调谐1.55μm光子晶体滤波器 . 红外与激光工程, 2017, 46(6): 620002-0620002(5). doi: 10.3788/IRLA201746.0620002
    [15] 何志平, 秦侠格, 徐映宇.  面阵凝视声光光谱成像技术及其轻小型无人机载应用探讨 . 红外与激光工程, 2016, 45(2): 203001-0203001(7). doi: 10.3788/IRLA201645.0203001
    [16] 米仁杰, 万助军, 汪涵.  MEMS可调谐平顶窄带光学滤波器 . 红外与激光工程, 2016, 45(7): 720001-0720001(5). doi: 10.3788/IRLA201645.0720001
    [17] 侯晴宇, 张树青.  单探测器共孔径多光谱成像系统设计 . 红外与激光工程, 2015, 44(5): 1638-1642.
    [18] 刘书勤, 王斌永, 徐睿, 陈凯, 何志平.  基于多谱段集成检测的宽光谱AOTF性能测试系统 . 红外与激光工程, 2015, 44(4): 1343-1348.
    [19] 付安邦, 张怀东, 张新宇, 桑红石, 季安, 谢长生.  中红外波段可调谐液晶法布里珀罗高光谱成像探测仪 . 红外与激光工程, 2013, 42(7): 1853-1857.
    [20] 姜庆辉, 邱跃洪, 文延, 汪欢, 许维星.  AOTF偏振光谱成像数据采集系统设计 . 红外与激光工程, 2012, 41(1): 218-222.
  • 加载中
图(5) / 表(1)
计量
  • 文章访问数:  107
  • HTML全文浏览量:  56
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-08-20
  • 修回日期:  2021-09-16
  • 网络出版日期:  2022-08-13
  • 刊出日期:  2022-08-05

AOTF后楔角对光谱图像像质的影响

doi: 10.3788/IRLA20210590
    作者简介:

    徐泽夫,男,博士生,主要从事AOTF器件建模以及AOTF成像光谱仪系统集成方面的研究

    赵慧洁,女,教授,博士生导师,博士,主要从事高光谱遥感、三维测量等方面的研究

    通讯作者: 贾国瑞,男,副教授,博士生导师,博士,主要从事高光谱成像建模、仿真和数据预处理等方面的研究。
基金项目:  国家重点研发计划(2016YFB0500505);国防基础科研计划(JCKY2019601C101);国家自然科学基金(61227806)
  • 中图分类号: O439;TP731

摘要: 针对AOTF成像光谱仪的光谱图像退化程度难以定量化评价的问题,提出了包含后楔角影响的AOTF器件调制传递函数(MTF)定量计算方法。该方法通过建立AOTF的光谱-空间维响应模型获取AOTF器件的线扩散函数(LSF),然后通过线扩散函数傅里叶变换得到理论MTF。在验证实验中,截止频率内的MTF理论值与实测值相对偏差小于15%,可以为AOTF的像质评估提供理论依据。进而,通过该定量计算方法讨论不同AOTF后楔角对光谱图像像质的影响。通过仿真分析得出,AOTF的后楔角难以同时满足横向色差校正与图像退化抑制,需要根据实际应用对于空间分辨率以及横向色差的要求进行取舍。所以,针对AOTF不同后楔角进行光谱图像像质的定量分析对于评估AOTF光谱图像像质具有重要意义,为AOTF器件设计提供了理论依据。

English Abstract

    • 声光可调谐滤波器(AOTF)是一种电控可调谐的分光器件,在凝视型成像光谱探测领域已有较多成功的应用[1-5]。AOTF器件的性能参数与AOTF器件的设计参数相关,所以建立AOTF器件的响应模型可以在设计阶段进行性能指标估计,为实际应用需求提供保障。为了描述AOTF器件的响应特性,科研人员在AOTF调谐曲线关系、光谱传输函数等方面已开展了很多研究。例如,Chang[6]、Yano[7]以及Gass[8]等人最早开展了AOTF光谱传输函数与调谐曲线关系的研究,Voloshinov、Krauz[9-11]等人也对AOTF光谱传输函数以及不同晶体材料的光谱响应进行了分析。Gorevoy[12]讨论了不同入射角下AOTF光谱响应的变化,Mantsevich[13]讨论了温度对调谐曲线的影响,Chang[14]讨论了双AOTF级联对光谱响应的影响。所以,对光谱维的AOTF响应已经有了较为成熟的研究。但是与此同时,对于光谱成像应用,AOTF器件由于其声光作用还会产生图像退化,主要体现在声光互作用方向上,使得子午与弧矢方向的光学空间分辨率不一致。这种图像退化不仅影响了光谱图像的空间分辨率,也让临近像元的光谱产生了串扰,即相邻像元的光谱混叠。但目前这方面的定量计算还相对缺乏,主要的研究集中在基于定性分析的工程抑制方法上。例如,针对图像退化分别提出了远心共焦光路[15]、光谱图像复原方法[16]、增加换能器长度[17]、AOTF级联[18]等方式来获得更清晰的图像。但是以上几种方式都有一定的局限性。远心共焦光路尽管图像退化抑制较好,但是同时会带来光谱带宽展宽且轴向色差不易抑制的问题。光谱图像复原方法在图像退化程度较小时会有比较好的复原效果,但是当图像退化程度较严重时,高频分量严重丢失,该方法并不能很好的还原出截止频率以外的高频信息。增加换能器长度和AOTF级联的方式都是通过减小光谱带宽的方式,使图像变得更为锐利,但是这两种方法会抑制AOTF的辐通量,使得系统信噪比降低。为了寻求更好的图像退化抑制方法,目前有研究发现改变AOTF后楔角可以抑制准直光路中的图像退化现象[19],这说明AOTF后楔角的不同也会引起AOTF器件的调制传输函数(MTF)改变。但是,AOTF后楔角的改变对像质的影响还缺少定量计算方法。目前公开文献中,对于AOTF空间维响应的定量计算主要是Machikhin[20]在研究中指出AOTF器件由于声光作用存在畸变的影响,但该研究未讨论AOTF器件的图像退化问题。Voloshinov[21]讨论了AOTF晶体内部的带宽内衍射光衍射角漂移,但该计算进行了大量线性近似,并且并未计算晶体内衍射光经过后楔角后的的角漂移量。

      综上所述,现有的国内外研究主要针对AOTF器件的调谐曲线、光谱响应函数进行了分析,对于空间维响应的定量计算较少,对于考虑后楔角的AOTF像质评估更是缺乏。为了能够在AOTF设计阶段定量计算AOTF器件后楔角对像质的影响,文中提出了不同波长下AOTF器件的光学调制函数计算方法,用以评估不同后楔角下的像质变化。该方法为AOTF器件优化减小工程迭代成本提供了依据,也为提高AOTF成像系统总体设计的效率提供了可能性。

    • 声光效应主要有三种不同的分类方式[22]:据入射光与衍射光的偏振特性分,可以分为正常声光效应与反常声光效应两类。正常声光效应中衍射光的偏振方向与入射光相同,而反常声光效应中衍射光的偏振方向与入射光不同。根据声波传输的不同分为表面波以及体波声光效应两类。根据声光互作用长度分类,可以分为拉曼奈斯声光效应与布拉格声光效应。AOTF是一种基于声光效应的带通滤光片,属于反常声光效应中的布拉格衍射体波器件,具有入射光与出射光偏振态不同,声光互作用长度较长,主要产生±1级衍射光,声波在晶体内传输等特性。目前,在可见光到短波红外范围内主要采用的声光材料是二氧化碲(TeO2)晶体。

      AOTF声光作用需要满足晶体内的“声光栅”与入射光产生动量匹配[6],使得入射光波矢与声波矢满足几何约束关系,如图1所示。

      图  1  动量匹配的波矢量布局图

      Figure 1.  Wave vector layout of momentum matching

      图1中的矢量关系可以表述为公式(1)。当入射光为O光时,则满足如图1中红色矢量的约束,当入射光为E光时,则满足图1中蓝色矢量的约束。

      $$ {{\boldsymbol{k}}_{{i}}} \pm {{\boldsymbol{k}}_{{a}}}{\boldsymbol{ = }}{{\boldsymbol{k}}_{{d}}} $$ (1)

      式中:$\left| {{{\boldsymbol{k}}_{{i}}}} \right| = \dfrac{{2\pi {n_i}}}{\lambda }$$\left| {{{\boldsymbol{k}}_{{d}}}} \right| = \dfrac{{2\pi {n_d}}}{\lambda }$$\left| {{{\boldsymbol{k}}_{{a}}}} \right| = \dfrac{{2\pi {f_a}}}{{{V_a}}}$$ \lambda $为波长,$ {n_i} $$ {n_d} $分别为入射光以及衍射光波矢量方向的折射率,$ {f_a} $$ {V_a} $分别为声波频率与声波相声速。

    • AOTF成像光谱仪获取的光谱图像质量取决于两个因素。一是光学系统的口径限制了衍射极限,且光学系统本身存在几何像差,该部分与成像光谱仪光学设计参数相关,可根据光学系统设计参数用光学设计商业软件进行MTF计算,文中不予讨论;二是由于AOTF器件自身的声光作用导致光谱图像在子午方向和弧矢方向的锐利程度不一致。由于AOTF声波矢量的方向并不单一,当入射光与理想声矢量不满足动量匹配时,依旧会有声光作用产生,引起光谱图像像质退化,此时的波矢量布局图如图2所示。

      图  2  非理想声束的声光互作用

      Figure 2.  Acousto-optic interaction of non ideal sound beam

      图2波矢量布局图中的入射光矢量与波矢面产生交点,计算得到入射光波矢末端的坐标$\left( {{x_c}^\prime ,{y_c}^\prime ,{z_c}^\prime } \right)$,其中,${x_c}^\prime = 0$$\left( {{y_c}^\prime ,{z_c}^\prime } \right)$由公式(2)计算,沿着110轴为y轴,沿着001轴为z轴:

      $$ \left\{ \begin{array}{l} \dfrac{{{y_c}{'^2}}}{{n_{{e}}^2}} + \dfrac{{{z_c}{'^2}}}{{n_{{o}}^2}} = {\left( {\dfrac{{2\pi }}{\lambda }} \right)^2} \\ \dfrac{{{y_c}'}}{{{M_1}}}{\text{ = }}\dfrac{{{z_c}'}}{{{N_1}}} \end{array} \right. $$ (2)

      式中:${n_e},{n_o}$分别为双折射晶体110轴E光折射率与O光折射率;$ {M_1},{N_1} $分别为入射光向量的y轴和z轴方向余弦。在确定入射光矢量后,入射光矢量与非理想声矢量产生声光作用,从而确定晶体内的衍射光矢量方向。根据波矢量布局的几何关系约束,可得到晶体内衍射光出射角的隐式解析表达,如公式(3)所示:

      $$ \left\{ \begin{array}{l} \dfrac{{y{'^2}}}{{n_{{o}}^2}} + \dfrac{{z{'^2}}}{{n_{{o}}^2}} = {\left( {\dfrac{{2\pi }}{\lambda }} \right)^2} \\ \dfrac{{{{\left( {y' - {y_c}'} \right)}^2}}}{{{{\left( {\dfrac{{2\pi {f_a}}}{{{V_t}}}} \right)}^2}}} + \dfrac{{{{\left( {z' - {z_c}'} \right)}^2}}}{{{{\left( {\dfrac{{2\pi {f_a}}}{{{V_z}}}} \right)}^2}}} = 1 \end{array} \right. $$ (3)

      式中:$\left( {y',z'} \right)$为非理想声矢量末端坐标;${f_a}$为声频率;${V_t},{V_z}$分别为声波沿着晶体110轴以及001轴的相速度。以E光进O光出的AOTF为例,根据波矢面圆球的约束,在yoz平面上$ z' $可以表示为:

      $$ z' = \sqrt {{{\left( {\frac{{2\pi }}{\lambda }} \right)}^2}n_{{o}}^2 - y{'^2}} $$ (4)

      $ z' $代入等声频率椭球的公式,则求解等声频率椭圆与波矢面圆球的交点问题变为求解方程:

      $$ \dfrac{{{{\left( {y' - {y_c}'} \right)}^2}}}{{{{\left( {\dfrac{{2\pi {f_a}}}{{{V_t}}}} \right)}^2}}} + \frac{{{{\left( {\sqrt {{{\left( {\dfrac{{2\pi }}{\lambda }} \right)}^2}n_{{o}}^2 - y{'^2}} - {z_c}'} \right)}^2}}}{{{{\left( {\dfrac{{2\pi {f_a}}}{{{V_z}}}} \right)}^2}}} - 1 = 0 $$ (5)

      该方程可以通过最优化迭代求解快速得到衍射光的方向余弦$\left( {0,y',z'} \right)$数值解,再根据晶体内衍射光计算晶体外的衍射光方向。在获取动量失配的光线衍射方向后,只要获取该失配光线的相对衍射效率即可得到空间维的线扩散函数。由于非理想声矢量的强度较弱,所以该波长的衍射效率也会较低。为了计算失配光线的衍射效率需要先引入动量失配量$\Delta \boldsymbol k$,动量失配量可表示为:

      $$ \Delta \boldsymbol k{\text{ = }}\sqrt {{{\left( {y' - {y_c}^\prime } \right)}^2}{\text{ + }}{{\left( {z' - {z_c}^\prime } \right)}^2}} $$ (6)

      该失配量越大,则代表与入声光产生声光作用的声矢量越远离理想声矢量方向。当失配量较小时,即$ \left( {y',z'} \right) $$ \left( {{y_c}^\prime ,{z_c}^\prime } \right) $较近时,此时$\Delta \boldsymbol k$方向近似为等声频率椭球泰勒级数展开的一阶项,与声群速度方向相互垂直。在计算得到动量失配量后,即可根据公式(7)计算动量失配下的相对衍射效率:

      $$ \eta = {u^2}\dfrac{{{{\sin }^2}\left( {\sqrt {{u^2} + {{\left( {\dfrac{{\Delta kL}}{2}} \right)}^2}} } \right)}}{{{u^2} + {{\left( {\dfrac{{\Delta kL}}{2}} \right)}^2}}} $$ (7)

      式中:$ u \sim \sqrt {{\eta _0}} $${\eta _0}$为峰值衍射效率;$L$为声光互作用长度,声光互作用长度与换能器长度成正比。当衍射光的出射方向以及相对衍射效率确定后,即可确定AOTF在某个固定频率下的光谱-空间维耦合响应。不同的入射波长对应不同的衍射方向以及衍射效率。晶体内的衍射光方向通过波矢量布局的约束条件进行计算,而空气中的衍射光方向则会受到AOTF后楔角的影响,如图3所示。

      图  3  AOTF声光作用以及后楔角的色散

      Figure 3.  Acousto-optic interaction of AOTF and dispersion of rear cut angle

    • 可见光AOTF一般采用的晶体为TeO2,这是一种双折射晶体,双折射晶体本身存在色散,该色散与声光作用带宽内的色散耦合,导致AOTF衍射光在空气中的出射角同时受到声光作用和AOTF后楔角影响。根据公式(5)、(7)进行数值计算可以得到AOTF不同波长衍射光在空气中出射方向以及衍射效率,从而获取AOTF的光谱-空间维响应。文中以E光进O光出的AOTF为例,超声切角为6.5°,入射光与001轴夹角为15°,后楔角为−4.8°,80 MHz以下声光互作用长度为9 mm,80 MHz以上的声光互作用长度为2.6 mm,数值计算得到的光谱-空间维响应如图4(a),(e)红线所示。将光谱-空间维响应投影在光谱维即可得到AOTF的光谱响应函数,如图4(a),(e)蓝线所示。同理,将光谱-空间维响应投影在空间维即可得到AOTF的空间维响应函数,如图4(a),(e)绿线所示。根据空间维响应函数可以计算得到AOTF器件的线扩散函数(LSF),LSF通过傅里叶变换后可取得AOTF器件的MTF曲线,如图4(b),(f)蓝线所示。该MTF为AOTF器件由于声光作用产生的MTF,并不包含光学系统像差以及衍射极限的影响。

      图  4  (a) 59.6 MHz (中心波长750 nm)的光谱-空间维响应(红线)、光谱维响应(蓝线)、空间维响应(绿线);(b)MTF曲线,蓝线为59.6 MHz空间维响应的傅里叶变换,红线为斜刀口靶实测MTF;(c) MTF实测值与理论值的相对偏差@59.6 MHz;(d) 59.6 MHz的光谱图像; (e) 99.4 MHz (中心波长500 nm)的光谱-空间维响应(红线)、光谱维响应(蓝线)、空间维响应(绿线);(f) MTF曲线,蓝线为99.4 MHz空间维响应的傅里叶变换,红线为斜刀口靶实测MTF; (g) MTF实测值与理论值的相对偏差@99.4 MHz;(h) 99.4 MHz的光谱图像

      Figure 4.  (a) Response of spectrum-space dimensional (red line), spectral response (blue line), spatial response (green line) at 59.6 MHz (central wavelength 750 nm); (b) MTF curve, blue line is Fourier transform of 59.6 MHz spatial response, red line is measured MTF by oblique knife edge target; (c) Deviation between the measured MTF and the theoretical MTF @59.6 MHz; (d) Spectral image of 59.6 MHz; (e) Response of spectrum-space dimensional (red line), spectral response (blue line), spatial response (green line) at 99.4 MHz (central wavelength 500 nm); (f) MTF curve, blue line is Fourier transform of 99.4 MHz spatial response, red line is measured MTF by oblique knife edge target; (g) Deviation between the measured MTF and the theoretical MTF @99.4 MHz; (h) Spectral image of 99.4 MHz

      为了验证理论计算的MTF,通过斜刀口靶测试获取了AOTF器件的MTF,如图4(b),(f)红线所示,测试得到的MTF与理论计算的MTF偏差小于0.1,通过$ \left({\text{MTF}}_{实测}-{\text{MTF}}_{理论}\right)/{\text{MTF}}_{理论} $计算得到实测与理论值的相对偏差,如图4(c),(g)所示,截止频率内的相对偏差小于15%。此处,实测值与理论值的偏差会受到多方面因素影响,例如刀口靶边缘锐利程度、探测器噪声、图像非均匀性均会影响MTF的实测值。

      结果表明,对于文中所采用的AOTF设计参数,59.6 MHz的MTF比99.4 MHz的MTF更优,在光谱图像上则应该表现为750 nm波长的图像比500 nm的光谱图像更为清晰。而且,无论是59.6 MHz还是99.4 MHz的MTF均存在截止频率。如果采用图像复原的方式,只能够在截止频率的带宽内放大高频信号,对于超出截止频率的信号,则难以再通过数字图像处理的手段进行复原。

    • 由前述分析可知:AOTF后楔角会影响衍射光在空气中的出射角度,所以在可见光范围内可利用合理的后楔角产生色散对声光作用的带宽内色散进行补偿,从而抑制光谱图像的退化,但现有的AOTF后楔角普遍被设计为补偿AOTF的横向色差以减少不同波段光谱图像之间的漂移。为了评估不同后楔角下AOTF器件的像质,文中以超声切角为6.5°的AOTF器件为例,根据前述理论方法计算出不同后楔角下的MTF。该AOTF波段范围400~1000 nm,光谱带宽2~10 nm,分为两个换能器通道,AOTF切型与声光互作用长度与前述示例一致。文中将不同后楔角的AOTF像质进行横向比较,后楔角范围为−8°~8°,其中−4.8°为针对该AOTF横向色差抑制所优化的最优角度。同时,文中分别计算了两个驱动频率下AOTF的MTF,99.4 MHz对应的声光互作用长度为2.6 mm,59.6 MHz对应的声光互作用长度为9 mm。MTF计算与线扩散函数相关,而线扩散函数的计算过程中会有系统等效焦距的影响,文中以50 mm等效焦距为例。

      图5可知,无论是短波还是长波,随着后楔角增加,AOTF的MTF随之增大,但是从截止频率的变化范围来看,改变后楔角对短波的影响更大。为了表征不同后楔角下AOTF像质优劣,文中采用MTF30作为量化指标,该指标取MTF等于30%处的空间频率作为评价像质的标准,30% MTF处的空间频率越高代表像质越好。文中计算不同波段、不同后楔角对应的MTF30如表1所示。

      图  5  (a) 59.6 MHz (中心波长750 nm)不同后楔角的MTF; (b) 99.4 MHz(中心波长500 nm)不同后楔角的MTF

      Figure 5.  (a) MTF of 59.6 MHz (central wavelength 750 nm) with different rear cut angles; (b) MTF of 99.4 MHz (central wavelength 500 nm) with different rear cut angles

      表 1  99.4 MHz (中心波长500 nm)和59.6 MHz (中心波长750 nm)不同后楔角的MTF30空间频率

      Table 1.  MTF30 with different rear cut angles at 99.4 MHz (center wavelength 500 nm) and 59.6 MHz (center wavelength 750 nm)

      99.4 MHz59.6 MHz
      Rear cut angle/(°)Spatial frequency/cycles·mm−1 @MTF30Rear cut angle/(°)Spatial frequency/cycles·mm−1 @ MTF30
      8 75.9 8 55.7
      6 45.6 6 50.9
      4 32.6 4 47.1
      2 25.6 2 43.4
      0 21.0 0 40.2
      −2 17.7 −2 37.4
      −4 15.3 −4 34.7
      −4.8 14.5 −4.8 33.7
      −6 11.7 −6 32.2
      −8 7.2 −8 29.9

      在实际应用中,一般对成像系统的奈奎斯特频率处MTF更为关注,例如当应用要求奈奎斯特频率处MTF优于30%时,则可以通过查找表中的空间频率推导出满足应用需求的AOTF后楔角角度。同时,通过分析发现在现有AOTF −4.8°后楔角下的MTF30比0°后楔角的更低。所以,横向色差校正与更优的空间响应本身是一对矛盾,无法同时满足。当横向色差得到校正,则MTF会下降。如果针对MTF来优化后楔角,则横向色差会引起不同波长的光谱图像漂移。所以,设计阶段的MTF定量计算十分重要,可用以保证AOTF的后楔角设计能满足实际应用需求。

    • AOTF光谱图像退化是AOTF成像光谱仪的固有现象。文中聚焦于AOTF器件的空间维响应,率先提出了AOTF器件的空间响应定量计算方法。与现有的AOTF模型相比,文中的模型同时考虑了声光作用导致的带宽内色散以及后楔角的色散,根据AOTF器件的理论MTF判断AOTF光谱成像的像质是否改善。基于此方法,以超声切角为6.5°的AOTF为例,计算了该AOTF的理论MTF,理论结果与实际测试的MTF保持了一致,MTF偏差小于0.1,相对偏差小于15%。该方法对于AOTF器件的设计具有重要意义,基于文中方法,可在设计阶段就进行AOTF的像质评价,减少AOTF器件切型优化的工程迭代成本。最后,针对超声切角6.5°的AOTF分析了不同后楔角下的MTF,指出MTF与横向色差之间是一对矛盾的性能指标。如果针对MTF进行后楔角优化,则无法校正横向色差,但对于横向色差要求较低的应用,可采用定标的方式进行横向色差校正,然后针对光谱图像像质优化AOTF后楔角,减少邻近像元的光谱混叠。文中针对AOTF提出MTF理论模型,为AOTF器件的后楔角选择提供了重要参考,也为获取高空间分辨率的光谱图像奠定了基础。

参考文献 (22)

目录

    /

    返回文章
    返回