留言板

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

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

地球临边场景红外遥感成像仿真方法

陈雪琦 万敏杰 许运凯 钱惟贤 陈钱 顾国华

陈雪琦, 万敏杰, 许运凯, 钱惟贤, 陈钱, 顾国华. 地球临边场景红外遥感成像仿真方法[J]. 红外与激光工程, 2022, 51(2): 20210896. doi: 10.3788/IRLA20210896
引用本文: 陈雪琦, 万敏杰, 许运凯, 钱惟贤, 陈钱, 顾国华. 地球临边场景红外遥感成像仿真方法[J]. 红外与激光工程, 2022, 51(2): 20210896. doi: 10.3788/IRLA20210896
Chen Xueqi, Wan Minjie, Xu Yunkai, Qian Weixian, Chen Qian, Gu Guohua. Infrared remote sensing imaging simulation method for earth’s limb scene[J]. Infrared and Laser Engineering, 2022, 51(2): 20210896. doi: 10.3788/IRLA20210896
Citation: Chen Xueqi, Wan Minjie, Xu Yunkai, Qian Weixian, Chen Qian, Gu Guohua. Infrared remote sensing imaging simulation method for earth’s limb scene[J]. Infrared and Laser Engineering, 2022, 51(2): 20210896. doi: 10.3788/IRLA20210896

地球临边场景红外遥感成像仿真方法

doi: 10.3788/IRLA20210896
基金项目: 国家自然科学基金青年科学基金(62001234);江苏省自然科学基金青年项目( BK20200487);中国博士后科学基金面上项目(2020M681597);江苏省博士后科研资助(2020Z051);上海航天科技创新基金(SAST2020-071);中央高校基本科研业务费专项(JSGP202102)
详细信息
    作者简介:

    陈雪琦,女,博士生,主要从事红外图像处理与目标检测方面的研究

    通讯作者: 万敏杰,男,副教授,博士,主要从事图像处理、计算机视觉和计算成像方面的研究。; 顾国华,男,教授,博士,主要从事研究领域为光学设计、计算机视觉与机器学习方面的研究。
  • 中图分类号: TP722.5

Infrared remote sensing imaging simulation method for earth’s limb scene

Funds: National Natural Science Foundation of China Youth Science Foundation(62001234);Youth Project of Jiangsu Natural Science Foundation( BK20200487);General Project of China Postdoctoral Science Foundation(2020M681597);Postdoctoral Research Fund of Jiangsu Province(2020Z051);Shanghai Aerospace Science and Technology Innovation Fund(SAST2020-071);Basic Research Funds for Central universities(JSGP202102)
  • 摘要: 地球临边场景仿真是卫星红外探测领域的关键组成部分,是空中高速目标远距离探测场景模拟的重要基础。临边观测下的地球表面近似于球面,传统的基于海洋三维形态并计算表面辐射特性的海洋红外图像仿真方法不适用。云层的厚度和高度对红外辐射传输特性的计算有重要影响,视云层为粒子团的处理方法会大大降低仿真的计算速度。因此,研究了海洋和云的红外辐射模型、地球-空间坐标系与红外相机坐标系的转换关系和大气传输模型,提出地球临边场景红外遥感成像仿真方法。根据场景组分的差异,分别建立海洋分布模型、多层云分布模型,并根据海洋和云层的红外辐射与反射特性,构建地球临边场景红外辐射模型。通过地球-空间坐标系与相机坐标系的转换关系,利用大气传输理论和传感器效应仿真,计算各观测角度的地球临边场景卫星遥感红外仿真图像。实验结果表明:仿真得到的红外图像画质清晰,符合地球临边场景红外辐射特性,其平均拉普拉斯算子和可达0.15,平均灰度梯度可达0.70。
  • 图  1  中点偏移算法示意图

    Figure  1.  Schematic diagram of midpoint displacement

    图  2  地球临边场景红外图像仿真示意图

    Figure  2.  Schematic diagram of infrared image simulation of earth limb scene

    图  3  坐标系转换关系示意图

    Figure  3.  Schematic diagram of transformation between coordinate systems

    图  4  地球临边场景红外仿真效果示意图(无噪声)

    Figure  4.  Schematic diagram of infrared simulation effect of earth limb scene (Without noise)

    图  5  地球临边场景红外仿真效果示意图(有噪声)

    Figure  5.  Schematic diagram of infrared simulation effect of earth limb scene (Including noise)

    表  1  仿真参数

    Table  1.   Simulation parameters

    Parameters(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)
    Focal length/mm202020204040404060606080
    Obliquity of satellite orbit/rad000020202020020020
    Altitude of satellite orbit/km2 0002 0002 0002 00010 00010 00010 00010 00020 00020 00020 00020 000
    Longitude of virtual target/rad5π/125π/127π/127π/125π/125π/127π/127π/125π/125π/127π/127π/12
    Latitude of virtual target/radπ/18π/6π/18π/6π/18π/6π/18π/6π/18π/6π/18π/6
    Altitude of virtual target/km1002002002005005003 0003 0002002005 0005 000
    Longitude of direct sunlight/rad02π/34π/302π/34π/302π/34π/302π/34π/3
    Latitude of direct sunlight/rad−π/9−π/9−π/9000π/9π/9π/9−π/90π/9
    下载: 导出CSV

    表  2  图像评价指标

    Table  2.   Image evaluation metrics

    Evaluation index
    image in figures
    LSGMGRuntime/s
    Fig.4Fig.5Fig.4Fig.5
    (a)0.0850.0780.340.3398.57
    (b)0.0900.0800.440.4098.71
    (c)0.110.100.500.50164.22
    (d)0.130.110.710.6187.30
    (e)0.130.120.580.53105.50
    (f)0.140.120.630.57111.58
    (g)0.170.150.810.76154.30
    (h)0.160.150.850.77108.92
    (i)0.220.230.981.16156.75
    (j)0.190.170.940.85139.99
    (k)0.190.190.930.99144.22
    (l)0.140.130.680.68111.33
    Average0.150.140.700.68123.45
    下载: 导出CSV
  • [1] Zhang Jing, Wang Shurong, Huang Yu, et al. Status and development of limb imaging spectrometers [J]. Chinese Optics, 2013, 6(5): 692-700. (in Chinese)
    [2] Wang Yan, Xie Xiaofang, Yang Jian, et al. The research on dynamic simulation of IR sea background scene [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2015, 35(5): 199-202+207. (in Chinese)
    [3] Song Bo, Cui Wenyu, Du Lili. et al. High-resolution satellite ocean background imaging simulation method [J]. Infrared and Laser Engineering, 2021, 50(9): 20200514. (in Chinese)
    [4] Shi Bo, Li Rui, Liu Chunsheng, et al. Research on earth′s atmosphere limb background infrared radiation characteristics [J]. Laser & Optoelectronics Progress, 2015, 52(12): 32-37. (in Chinese)
    [5] Wang Feixiang, Guo Jie, Xu Fangyu, et al. Calculation and measurement of infrared atmosphere transmittance at difference altitudes [J]. Chinese Optics, 2019, 12(4): 844-853. (in Chinese)
    [6] Nicodemus F E. Directional reflectance and emissivity of an opaque surface [J]. Applied Optics, 1965, 4(7): 767-773. doi:  10.1364/AO.4.000767
    [7] Cook R L. A reflectance models for computer graphics [J]. Acm Transactions on Graphics, 1982, 1(1): 7-14.
    [8] Wang Chunlei, Li Fuxing, Ma Junjun. Estimation of downwelling surface longwave radiation under cirrus cloud with remotely sensed data [J]. Journal of Henan University (Nature Science), 2021, 51(2): 183-192. (in Chinese)
    [9] Wang Feng, Niu Shibo, Yue Chengfei, et al. Design of attitude control system for ASRTU microsatellite [J]. Optics and Precision Engineering, 2020, 28(10): 2192-2202. (in Chinese) doi:  10.37188/OPE.20202810.2192
    [10] Liu Lianwei, Dong Shikui, Chen Qianrong, et al. Infrared radiation imaging calculation of aerial target based on CUDA parallel computing [J]. Infrared and Laser Engineering, 2020, 49(4): 20200514. (in Chinese)
    [11] Yu Kun, Guo Biao, Cong Mingyu. Infrared imaging modelling and image simulation of limb background for space target detection [J]. Infrared and Laser Engineering, 2019, 48(9): 0904005. (in Chinese)
    [12] Liu Dong, Dai Congming, Wei Heli. Comparison of LTE and non-LTE model for limb infrared radiation simulation of middle and upper atmosphere [J]. Journal of Atmospheric and Environmental Optics, 2019, 14(5): 337-344. (in Chinese)
    [13] Xue Qingsheng, Wang Shurong, Li Futian, et al. Limb imaging spectrometer for atmospheric remote sensing [J]. Optics and Precision Engineering, 2010, 18(4): 823-830. (in Chinese)
    [14] Cui Guangmang, Zhang Keqi, Mao Lei, et al. Micro-image definition evaluation using multi-scale decomposition and gradient absolute value [J]. Opto-Electronic Engineering, 2019, 46(6): 59-69. (in Chinese)
  • [1] 魏合理, 戴聪明, 武鹏飞, 李建玉, 孙凤萤, 唐超礼, 黄宏华, 李学彬, 朱文越, 饶瑞中, 王英俭.  通用大气辐射传输软件CART在场景计算中的应用(特邀) . 红外与激光工程, 2022, 51(5): 20210916-1-20210916-10. doi: 10.3788/IRLA20210916
    [2] 吴双, 李超, 高传卫, 佟岐.  时序平滑多尺度叠加动态红外云场景仿真 . 红外与激光工程, 2022, 51(8): 20220656-1-20220656-5. doi: 10.3788/IRLA20220656
    [3] 于鲲, 段雨晗, 丛明煜, 戴文聪.  飞行器红外物理成像仿真优化计算方法 . 红外与激光工程, 2021, 50(4): 20200241-1-20200241-14. doi: 10.3788/IRLA20200241
    [4] 宋波, 方薇, 杜丽丽, 崔文煜, 王涛, 易维宁.  海面目标高分辨率卫星成像仿真方法 . 红外与激光工程, 2021, 50(12): 20210127-1-20210127-14. doi: 10.3788/IRLA20210127
    [5] 宋波, 崔文煜, 杜丽丽, 周川杰, 易维宁.  高分辨率卫星海洋背景成像仿真方法 . 红外与激光工程, 2021, 50(9): 20200514-1-20200514-12. doi: 10.3788/IRLA20200514
    [6] 吴俊政, 倪维平, 严卫东, 张晗.  圆周阵列太赫兹干涉成像中目标场景仿真 . 红外与激光工程, 2019, 48(1): 125004-0125004(8). doi: 10.3788/IRLA201948.0125004
    [7] 于鲲, 郭彪, 丛明煜.  空间目标临边探测背景红外成像建模与图像仿真 . 红外与激光工程, 2019, 48(9): 904005-0904005(10). doi: 10.3788/IRLA201948.0904005
    [8] 李霞, 刘建国, 董雁冰, 刘兴润.  基于遥感数据的地球背景中红外场景仿真 . 红外与激光工程, 2018, 47(11): 1104004-1104004(7). doi: 10.3788/IRLA201847.1104004
    [9] 鲁佳颖, 谷小婧, 顾幸生.  面向微光/红外融合彩色夜视的场景解析方法 . 红外与激光工程, 2017, 46(8): 804002-0804002(6). doi: 10.3788/IRLA201746.0804002
    [10] 谢锋, 刘成玉, 邵红兰, 张长兴, 杨贵, 王建宇.  基于场景的热红外高光谱数据光谱定标 . 红外与激光工程, 2017, 46(1): 138001-0138001(6). doi: 10.3788/IRLA201746.0138001
    [11] 张作宇, 廖守亿, 张金城, 苏德伦, 闫循良.  基于物理模型的战场烟幕实时红外仿真 . 红外与激光工程, 2016, 45(4): 404004-0404004(8). doi: 10.3788/IRLA201645.0404004
    [12] 薛莲, 张力, 刘佳琪, 李志峰, 张陆萍.  基于遥感反演的空间下视系统红外场景仿真 . 红外与激光工程, 2016, 45(7): 704004-0704004(5). doi: 10.3788/IRLA201645.0704004
    [13] 许家林, 王晓东, 李丙玉, 王鹤, 孙强.  基于DMD的红外场景模拟器图像数据传输和分割存储方法 . 红外与激光工程, 2015, 44(9): 2622-2626.
    [14] 贺明, 王亚弟, 王新赛, 路建方, 李志军.  场景自适应的红外焦平面成像系统灰度超分辨技术 . 红外与激光工程, 2014, 43(7): 2138-2142.
    [15] 郭凯, 孙红胜, 张虎, 张玉国, 魏树弟, 孙广尉.  仿真用红外动态场景模拟器现场校准技术 . 红外与激光工程, 2014, 43(5): 1380-1384.
    [16] 彭晨, 陈钱, 钱惟贤, 徐富元.  复杂地面场景下的红外运动目标跟踪 . 红外与激光工程, 2013, 42(6): 1410-1414.
    [17] 何永强, 唐德帅, 胡文刚.  基于DMD的红外场景仿真系统投影光路消热差设计 . 红外与激光工程, 2013, 42(9): 2319-2323.
    [18] 方义强, 樊祥, 程正东, 朱斌, 张发强, 施展.  晴空背景下地基红外告警系统作用距离 . 红外与激光工程, 2013, 42(7): 1654-1659.
    [19] 李卓, 钱丽勋, 欧文.  动态红外场景生成新技术 . 红外与激光工程, 2012, 41(1): 1-6.
    [20] 于洋, 汤心溢, 陈元林.  实时闭环红外场景仿真技术 . 红外与激光工程, 2012, 41(1): 15-19.
  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  189
  • HTML全文浏览量:  47
  • PDF下载量:  60
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-20
  • 修回日期:  2021-12-25
  • 网络出版日期:  2022-03-04
  • 刊出日期:  2022-02-28

地球临边场景红外遥感成像仿真方法

doi: 10.3788/IRLA20210896
    作者简介:

    陈雪琦,女,博士生,主要从事红外图像处理与目标检测方面的研究

    通讯作者: 万敏杰,男,副教授,博士,主要从事图像处理、计算机视觉和计算成像方面的研究。; 顾国华,男,教授,博士,主要从事研究领域为光学设计、计算机视觉与机器学习方面的研究。
基金项目:  国家自然科学基金青年科学基金(62001234);江苏省自然科学基金青年项目( BK20200487);中国博士后科学基金面上项目(2020M681597);江苏省博士后科研资助(2020Z051);上海航天科技创新基金(SAST2020-071);中央高校基本科研业务费专项(JSGP202102)
  • 中图分类号: TP722.5

摘要: 地球临边场景仿真是卫星红外探测领域的关键组成部分,是空中高速目标远距离探测场景模拟的重要基础。临边观测下的地球表面近似于球面,传统的基于海洋三维形态并计算表面辐射特性的海洋红外图像仿真方法不适用。云层的厚度和高度对红外辐射传输特性的计算有重要影响,视云层为粒子团的处理方法会大大降低仿真的计算速度。因此,研究了海洋和云的红外辐射模型、地球-空间坐标系与红外相机坐标系的转换关系和大气传输模型,提出地球临边场景红外遥感成像仿真方法。根据场景组分的差异,分别建立海洋分布模型、多层云分布模型,并根据海洋和云层的红外辐射与反射特性,构建地球临边场景红外辐射模型。通过地球-空间坐标系与相机坐标系的转换关系,利用大气传输理论和传感器效应仿真,计算各观测角度的地球临边场景卫星遥感红外仿真图像。实验结果表明:仿真得到的红外图像画质清晰,符合地球临边场景红外辐射特性,其平均拉普拉斯算子和可达0.15,平均灰度梯度可达0.70。

English Abstract

    • 近年来,利用红外卫星临边观测获取的地球临边场景遥感图像对空中运动目标进行探测和识别是遥感技术的重要应用方向之一。临边观测模式得到的图像具有高空间覆盖范围、高垂直分辨率的特点[1]。然而,公开可获取的图像较少,且难以覆盖不同观测角度、气候条件、云区分布、观测区域的情况。地球临边场景红外仿真可以为空中运动目标检测识别算法的研究提供不同条件下的数据样本。

      目前,地球临边场景仿真的相关研究包括海洋场景仿真,大气场景红外辐射特性研究两个方面。针对海洋场景仿真,王彦等[2]根据海面微面法向量计算海洋辐射特性,渲染海面纹理,计算出红外图像;宋波等[3]提出了一种高分辨率海洋背景成像仿真方法,通过亚米级海面零视距辐射场,模拟不同海况下的卫星遥感图像。这些研究集中在如何增强海洋红外图像的波纹细节信息,没有考虑到大气环境对海洋辐射模型的复杂影响,不适用于远距离、大范围的卫星遥感领域。针对大气场景红外辐射特性,史波等[4]建立了临边大气红外辐射计算模型,研究了影响大气红外辐射特性的主要因素;王飞翔等[5]研究了不同海拔下的大气中波红外透过率。针对大气场景红外辐射的研究主要集中于地面观测角度,不同于卫星遥感领域需要考虑地表出射的红外辐射对大气、云层的影响。

      地表包括陆地和海洋,不同于在陆地上可以使用陆地-卫星联合探测,当飞行器目标飞过海洋上空时,对空中目标的观测主要依赖于卫星。深空背景是地球大气层外附近的宇宙空间,其中的红外辐射主要来源于大气层反射的太阳辐射和大气层逸散的辐射。宇宙中除太阳以外的恒星等红外辐射源距离较远,其影响可以忽略不计。因此,文中的研究对象集中在海洋、大气和云。文中研究了海洋和云的红外辐射模型、地球-空间坐标系与相机坐标系的转换关系、大气传输模型,提出了地球临边场景红外遥感成像仿真方法,可以为空中高速目标探测提供有力的数据支撑。

      文中主要贡献有:

      (1)提出了一种根据云覆盖域的变化计算海洋红外辐射和大气辐射的图像仿真模型,能够反映云区对海洋接收的红外辐射造成的削弱效果;

      (2)提出了一种根据空中目标实时位置进行场景定位的图像仿真模型,能够仿真卫星对不同地理位置、不同飞行高度的空中目标进行遥感观测的情况;

      (3)根据太阳、地球、卫星探测器的位置变化,实施图像仿真实验,能够快速提供卫星在不同探测条件下的地球临边场景图像。

    • 地球临边场景的红外辐射包含海洋和云的红外辐射。海洋出射的红外辐射包括对太阳辐射的反射,对大气辐射的反射和海洋自发辐射。大气层对红外图像的贡献体现在云的成像、大气传输效应和大气自身热辐射。

    • 地球临边场景图像中每个像元对应的地表面积较大,海洋的纹理细节不明显。和传统的海洋图像仿真算法重视三维形态的仿真不同,地球临边场景中的海洋红外辐射模型主要考察海面在区域范围的辐射和反射特性,其中海洋自发辐射在海洋总体红外辐射中的比重最大。

    • 为了计算海面自身热辐射,可把海面近似作为灰体处理。首先计算同温度下的黑体辐射亮度,再将黑体辐射亮度与海面自发辐射率相乘,得到海面自发辐射亮度:

      $$ {L}_{h}={\varepsilon }_{sea}\cdot {L}_{BB}\left(\lambda ,T\right) $$ (1)

      式中:${\varepsilon }_{sea}=0.98\left[1-{\left(1-\mathrm{cos}\;\theta \right)}^{5}\right]$为海面自发辐射率,$ \theta $表示海平面法向量与探测方向的夹角;$ {L}_{BB}\left(\lambda ,T\right) $为黑体红外波段的自发辐射亮度。

      海面同温黑体的光谱辐射出射度为:

      $$ {M}_{\lambda bb}=\frac{{c}_{1}}{{\lambda }^{5}}\cdot \frac{1}{{\rm{exp}}\left({c}_{2}/\lambda T\right)-1} $$ (2)

      式中:${c}_{1}=\left(3.741\;5\pm 0.000\;3\right)\times {10}^{8}\;\mathrm{W}\cdot {\text{μ}\mathrm{m}}^{4}\cdot {\mathrm{m}}^{-2}$为第一辐射常数;${c}_{2}=\left(1.438\;79\pm 0.000\;19\right)\times {10}^{4}\;\text{μ}\mathrm{m}\cdot \mathrm{K}$为第二辐射常数;$\lambda$是辐射波长,单位:μm;T 代表海面的绝对温度,单位:K。在红外波段,黑体自发辐射亮度为:

      $$ {L}_{BB}\left(\lambda ,T\right)=\frac{{M}_{\lambda bb}}{\pi }={\int }_{{\lambda }_{1}}^{{\lambda }_{2}}\frac{1}{\pi }\frac{{c}_{1}}{{\lambda }^{5}}\cdot \frac{{\rm{d}}\lambda }{{\rm{exp}}\left({c}_{2}/\lambda T\right)-1} $$ (3)

      式中:${\lambda }_{1}$${\lambda }_{2}$分别为红外波段范围的下限和上限,单位:μm。海洋在$8\sim14\;\text{μ}\mathrm{m}$的自发热辐射强度比在$3\sim5\;\text{μ}\mathrm{m}$波段大一个数量级左右。

    • 计算海面的反射模型,首先计算海洋接收的红外辐照度,再将接收的红外辐照度与海面在探测方向的反射率相乘,得到海面反射辐射亮度。

      利用辐射传输软件MODTRAN计算在不同入射角度、不同地理条件下的太阳光谱辐照度$ {E}_{sun,\lambda } $,从而计算红外波段海面接收的太阳辐照度:

      $$ {E}_{sun}={\int }_{{\lambda }_{1}}^{{\lambda }_{2}}{E}_{sun,\lambda }{\rm{d}}\lambda $$ (4)

      平静的海面对光源的反射包括漫反射和镜面反射。使用双向反射分布函数(Bidirectional Reflectance Distribution Function, BRDF)可以用来描述海面的反射特性。BRDF描述了出射方向上的反射光强相对于入射辐照度的比例$ f $

      $$ f\left({\theta }_{i},{\varphi }_{i},{\theta }_{r},{\varphi }_{r},\lambda \right)=\frac{{\rm{d}}{L}_{r}\left({\theta }_{r},{\varphi }_{r},\lambda \right)}{{\rm{d}}{E}_{i}\left({\theta }_{i},{\varphi }_{i},\lambda \right)} $$ (5)

      式中:$ {\theta }_{i} $为入射天顶角;$ {\theta }_{r} $为出射天顶角;${\varphi }_{i}$为入射方位角;${\varphi }_{r}$为出射方位角;$ \lambda $为入射光波长;${L}_{r}$为出射方向的辐照度,单位:$\mathrm{W}/\left({\mathrm{m}}^{2}\cdot {\rm{sr}}\cdot {\text{μ}}\mathrm{m}\right)$$ {E}_{i} $为入射方向的辐照度,单位:$\mathrm{W}/\left({\mathrm{m}}^{2}\cdot {\text{μ}}\mathrm{m}\right)$。BRDF函数表示为:

      $$ f\left({\theta }_{i},{\varphi }_{i},{\theta }_{r},{\varphi }_{r},\lambda \right)=\frac{D\cdot G}{4\mathit{{\rm{cos}}}\;{\theta }_{i}\mathit{{\rm{cos}}}\;{\theta }_{r}} $$ (6)

      式中:D为坡度分布系数;G为几何衰减因子。DG的数值与海面的粗糙度、光线的入射方向与反射方向有关[6-7]

      根据海面单位面积接收的太阳辐照度与海面对太阳光的反射率,求得海面对太阳辐射的反射光强为:

      $$ {L}_{r,sun}=f\cdot {E}_{sun} $$ (7)

      同理可计算海面对大气背景辐射的反射模型。大气背景光来自各个方向,为了简化计算,使用4个方位上的45°天顶角入射的大气辐射强度代替整个大气背景光。根据双向反射分布原理,计算海面对大气背景光的反射光强,结合海面对太阳辐射的反射,计算海洋对红外辐射的反射光强。大气辐射在$8\sim14\;{\text{μ}}\mathrm{m}$的自发热辐射强度比在$3\sim5\;{\text{μ}}\mathrm{m}$波段大一个数量级。此外,还需要考虑云对太阳辐射的遮挡与对大气辐射的衰减[8]

    • 云的红外辐射包括云表层反射的太阳辐射、穿过云层出射的背景辐射、以及云的自发辐射。根据云的分布规律生成多层云图,结合云的辐射基本原理模型,计算红外辐射强度。

    • 云的自发辐射集中在$8\sim14\;{\text{μ}}\mathrm{m}$波段,采用灰体近似方法来计算。云的辐亮度公式为:

      $$ {L}_{cloud}=\left(1-{t}_{c}\right){L}_{bb}\left({T}_{cloud}\right) $$ (8)

      式中:$ {t}_{c} $为云的透过率;$ {L}_{bb}\left({T}_{cloud}\right) $是与云温度相同的黑体辐射亮度;$ {T}_{cloud} $为云的温度。黑体辐射计算公式同公式(3)。灰体的温度稳定时,吸收率等于发射率,吸收率与透过率之和为1,因此$ 1-{t}_{c} $可以看作是云的发射率。

      假设云的密度均匀且各向同性,视线上云的透过率可以通过下式来计算:

      $$ {t}_{c}={{\rm{e}}}^{-{\sigma }_{e}d} $$ (9)

      式中:$ {\sigma }_{e} $为云的容积消光系数;d为视线方向上的云层厚度。

      云对太阳辐射的反射遵循朗伯体公式:

      $$ {L}_{reflection}=\left(1-{t}_{c}\right){\rho }_{c}\frac{E}{\pi }\mathit{{\rm{cos}}}\;\theta $$ (10)

      式中:$ {\rho }_{c} $为云的反射系数;E为辐射源照射到云上的辐照度;$ \theta $为太阳辐射传输方向与云表面法向量的夹角。与自发辐射相比,云的反射在零视距总体红外辐射中所占比重较小。

      穿过云层出射的背景辐射包括海洋背景辐射和近地大气背景辐射。二者在云层中传输时经过吸收、散射等作用后大幅衰减,能够穿出云层出射的辐射强度非常小,可以忽略不计。

    • 采用中点偏移算法(Midpoint Displacement, MPD)来生成云的外观模型数据场。数据场中,数据分布范围代表云的分布范围,各点的数值代表云层厚度。算法基本过程如下:

      (1)给四边形的4个角点赋初始高度值,如图1 (a)所示。

      图  1  中点偏移算法示意图

      Figure 1.  Schematic diagram of midpoint displacement

      (2)在正方形的中点生成一个随机值,中点值为4个角点高度平均值叠加随机量,如图1 (b)所示。

      (3)取正方形中点和任意一条边的两个端点,计算高度平均值,叠加随机偏移量,分别得到四条边上的中点值,如图1 (c)所示。

      上述过程为0级递归,能将2×2区域扩充为3×3区域。重复上述过程,通过多级递归求得二维数据场d

      将数据场d转换为云的厚度场H,点$ \left(x,y\right) $处的云厚度值H为:

      $$ \begin{split} H\left(x,y\right)=\left\{\begin{array}{cc}{H}_{min}+\Delta H\dfrac{{{d}}{\text{'}}\left(x,y\right)-{{d}}_{avg}}{{{d}}_{max}-{{d}}_{avg}}& \left({{d}}{\text{'}}\left(x,y\right)\ne 0\right)\\ 0& \left({{d}}{\text{'}}\left(x,y\right)=0\right)\end{array}\right. \\[-10pt] \end{split}$$ (11)

      式中:$ {d}_{max} $$ {d}_{avg} $分别为$ d $的最大值和平均值;$ d\text{'} $是由$ d $计算得到的数据场,$ d\text{'} $中的0值像素是$ d $中小于平均值$ {d}_{avg} $的像素,代表无云点;${H}_{min}\sim{H}_{max}$代表云厚度的分布范围,$ \Delta H={H}_{max}-{H}_{min} $

      由此得到单层云的三维分布信息。地球上的云分布在不同的海拔高度。在高云、中云、低云的高度范围内分别选取代表值作为云层高度,仿真单层云的外观模型,结合云的红外辐射模型计算多个云层叠加后的仿真图像。

    • 地球临边场景红外图像仿真的实质是计算探测器各像素的灰度值。首先根据探测器的地理位置、探测方向,判断探测器各像元的探测对象是海洋或云层,或无对象(深空场景);其次,根据第1节所述的海洋场景和云场景红外辐射模型,结合红外辐射的大气传输理论计算探测对象的红外辐射在传输路径上的衰减,从而得到探测器各像元接收的红外辐射强度;最后将红外辐射强度信号转换为图像灰度值,得到地球临边场景红外仿真图像。仿真过程如图2所示。

      图  2  地球临边场景红外图像仿真示意图

      Figure 2.  Schematic diagram of infrared image simulation of earth limb scene

    • 地球-空间坐标系以地心为原点,可以直观地描述海洋与云层的分布情况,卫星的经纬度、高度和太阳直射点的经纬度。卫星坐标系以卫星为原点,光轴方向指向地球中心。相机坐标系以卫星搭载的红外传感器为原点,光轴方向指向空中虚拟目标。通过卫星姿态控制系统[9],可以调整传感器方向,将其对准虚拟目标,从而将卫星坐标系转换为相机坐标系。虚拟目标不参与红外成像的过程。

      通过计算地球-空间坐标系、卫星坐标系和相机坐标系之间的位置和方向,得到地球-空间坐标系与相机坐标系的转换关系,从而确定地球-空间坐标系中的任意点在相机坐标系中的位置。结合成像投影公式,得到该点在像面的投影位置,即可确定红外相机单个像元的探测对象。在计算中,将面元细分为更小面元,再计算面元成像位置,能够提高仿真图像的准确性[10]。坐标系之间的转换关系如图3所示,以地球中心(Center of earth and O-earth)为原点的地球-空间坐标系;以卫星(O-satellite)为原点,Z轴指向地球中心的卫星坐标系;以相机(O-camera)为原点,Z轴指向虚拟目标(Virtual target)的相机坐标系。卫星坐标系的中心与相机坐标系的中心可近似视为同一个点。

      图  3  坐标系转换关系示意图

      Figure 3.  Schematic diagram of transformation between coordinate systems

      相机坐标系中的点$\left(x,y,{\textit{z}}\right)$在像素坐标系中的位置$ \left(u,v\right) $遵循相机投影公式:

      $$ s\left[\begin{array}{c}u\\ v\\ 1\end{array}\right]=H\times RT\times \left[\begin{array}{c}\begin{array}{c}X\\ Y\end{array}\\ \begin{array}{c}Z\\ 1\end{array}\end{array}\right] $$ (12)

      式中:$ H=\left[\begin{array}{ccc}{f}_{x}& 0& {u}_{0}\\ 0& {f}_{y}& {v}_{0}\\ 0& 0& 1\end{array}\right] $表示相机内部参数;fxfy分别表示x轴、y轴上的焦距;(u0v0)表示像面的中心坐标;RT是一个3×4矩阵,描述了地球-空间坐标系与相机坐标系的转换关系;s为归一化常数。

    • 根据地球-空间坐标系、卫星坐标系和相机坐标系之间的转换关系,结合相机投影公式,判断卫星红外传感器单个像元的探测对象。根据探测对象的不同,计算探测器接收的红外辐射强度。当红外卫星传感器的探测对象是海洋时,传感器接收的辐射能量包括海洋表面零视距辐射出射度穿过整个大气层到达大气层外的部分,以及视线路径上各点的大气背景辐射到达大气层外的部分[11]

      $$ {L}_{s}={\tau }_{sea}\left({L}_{h}+{L}_{r,sun}+{L}_{r,sky}\right)+\int {\tau }_{sky}\cdot {E}_{sky} $$ (13)

      式中:$ {\tau }_{sea} $为海洋辐射在大气传输路径上的衰减率;$ {L}_{h} $为海面自发辐射;$ {L}_{r,sun} $为海面对太阳辐射的反射;$ {L}_{r,sky} $为海面对大气背景辐射的反射;$ {\tau }_{sky} $表示大气中某点在传输路径上的辐射衰减率;$ {E}_{sky} $为该点处的大气背景辐照度。

      当探测对象是云时,传感器接收的辐射能量为云层表面零视距辐射出射度和视线路径上各点的大气背景辐射到达大气层外的部分:

      $$ {L}_{c}={\tau }_{cloud}\left({L}_{cloud}+{L}_{reflection}\right)+\int {\tau }_{sky}\cdot {E}_{sky} $$ (14)

      式中:$ {\tau }_{cloud} $为云辐射在大气传输路径上的衰减率;$ {L}_{cloud} $表示云的自发辐射;$ {L}_{reflection} $为云对太阳辐射的反射;$ {\tau }_{sky} $$ {E}_{sky} $与公式(13)中的同名参数意义相同,而数值不同,这是因为不同层流的大气辐射强度有差异 [12]$\int {\tau }_{sky}\cdot {E}_{sky}$中包含单次和多次散射的太阳光谱辐射[13]

    • 传感器单像元接收的辐射强度光学信号经过光电转换和灰度映射转化为灰度值:

      $$ {V}_{s}={L}_{d\left(i,j\right)}\times {R}_{d}\times {A}_{d} $$ (15)

      式中:$ {V}_{s} $表示单像元的电压值;$ {L}_{d\left(i,j\right)} $表示像面上$ \left(i,j\right) $处的像元接收的辐射强度,根据探测对象的不同分别为$ {L}_{s} $(海洋)、$ {L}_{c} $(云)或0 (深空);$ {R}_{d} $为传感器的光谱响应率;$ {A}_{d} $为像元接收面积。电压经过灰度映射转化为灰度值:

      $$ F\left(i,j\right)=\frac{{V}_{s}+{V}_{n}-{V}_{min}}{{V}_{max}-{V}_{min}}\cdot M $$ (16)

      式中:$ F\left(i,j\right) $表示$ \left(i,j\right) $处的像元灰度值;$ {V}_{n} $为随机噪声;$ {V}_{max} $$ {V}_{min} $分别为附加随机噪声后的传感器最大电压和最小电压;M表示图像的灰度等级。随机噪声包括散粒噪声、热噪声,以及传感器的随机振动造成的像点模糊等。其中散粒噪声和热噪声都具有白噪声的特性,采用均值为0的高斯分布来模拟。

      像点模糊效果采用点扩散函数(Point Spread Function, PSF)来计算:

      $$ PSF\left({i}{\text{'}},{j}{\text{'}}\right)=\pi \cdot {p}^{2}\cdot {{\rm{e}}}^{-{\left(\pi \cdot p\right)}^{2}\cdot \left({i}^{2}+{j}^{2}\right)} $$ (17)

      式中:$ p $为光学响应指数;$\left({i}{\mathrm{\text{'}}},{j}{\mathrm{\text{'}}}\right)$表示像面上的点$ \left(i,j\right) $经过点扩散映射的位置。利用点扩散函数对图像进行卷积,可以模拟传感器随机振动的影响。

      根据海洋场景和云场景的红外辐射模型计算探测对象的零视距辐射出射度,经过坐标系转换、大气传输效应计算、传感器效应仿真,得到地球临边场景仿真图像。

    • 为了验证仿真模型的真实性和准确性,选取几组不同拍摄时间、探测方向、虚拟目标位置的仿真参数,模拟不同日照条件下的仿真图像。算法运行环境为Win10 64位操作系统,Intel i7-9700CPU,主存32 GB,软件平台为MATLAB 2012b。模型中的探测器像素尺寸为320×256,像元间距为20 μm,探测波段为8~14 μm。卫星处于东经60°上空。仿真参数设置见表1,仿真效果如图4图5所示。

      表 1  仿真参数

      Table 1.  Simulation parameters

      Parameters(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)
      Focal length/mm202020204040404060606080
      Obliquity of satellite orbit/rad000020202020020020
      Altitude of satellite orbit/km2 0002 0002 0002 00010 00010 00010 00010 00020 00020 00020 00020 000
      Longitude of virtual target/rad5π/125π/127π/127π/125π/125π/127π/127π/125π/125π/127π/127π/12
      Latitude of virtual target/radπ/18π/6π/18π/6π/18π/6π/18π/6π/18π/6π/18π/6
      Altitude of virtual target/km1002002002005005003 0003 0002002005 0005 000
      Longitude of direct sunlight/rad02π/34π/302π/34π/302π/34π/302π/34π/3
      Latitude of direct sunlight/rad−π/9−π/9−π/9000π/9π/9π/9−π/90π/9

      图  4  地球临边场景红外仿真效果示意图(无噪声)

      Figure 4.  Schematic diagram of infrared simulation effect of earth limb scene (Without noise)

      图  5  地球临边场景红外仿真效果示意图(有噪声)

      Figure 5.  Schematic diagram of infrared simulation effect of earth limb scene (Including noise)

      图4所示为不考虑噪声的仿真图像,图5所示为考虑传感器传输噪声与镜头振动的仿真图像。图中可以明显区分海洋、云和深空因红外辐射特性的不同而导致的图像灰度差异。云的图像存在单方向的条纹,离地球边缘越近,条纹现象越明显。在倾斜角度拍摄真实的云图时,存在单方向的条纹特性,而在地球临边观测中,同样采用倾斜角度来观测云,同样出现了条纹,符合观测经验。

      使用拉普拉斯算子和(Laplacian Sum, LS)、灰度梯度(Grayscale Mean Gradient, GMG)评价仿真图像的清晰度[14]。拉普拉斯算子和比较了3×3 pixel邻域内的像素偏差程度。首先使用拉普拉斯算子对图像求卷积,该算子定义如下:

      $$ L=\left[\begin{array}{ccc}-1& -1& -1\\ -1& 8& -1\\ -1& -1& -1\end{array}\right] $$ (18)

      基于拉普拉斯算子和的图像清晰度定义为:

      $$ {A}_{LS}=\frac{1}{\left(M-2\right)\left(N-2\right)}\sum\nolimits _{i=2}^{M-1}\sum\nolimits _{j=2}^{N-1}\left|G\left(i,j\right)\right| $$ (19)

      式中:$ \left|G\left(i,j\right)\right| $为像素点$ \left(i,j\right) $处的拉普拉斯算子卷积;MN为图像的行数和列数;LS值越大,图像信息越丰富,边缘越锐利,图像质量越好。

      灰度梯度函数反映了图像的细节和纹理变化特征,定义为:

      $$ \begin{split} &{A}_{GMG}=\frac{1}{\left(M-1\right)\left(N-1\right)}\times\\ &\sum\nolimits _{j}\sum\nolimits _{i} \sqrt{\frac{{\left|f\left(i,j\right)-f\left(i,j-1\right)\right|}^{2}+{\left|f\left(i,j\right)-f\left(i+1,j\right)\right|}^{2}}{2}} \end{split} $$ (20)

      式中:$ f\left(i,j\right) $为像素点$ \left(i,j\right) $处的图像灰度;GMG值越大,高频信息越丰富,图像越清晰。

      对仿真实验结果进行图像清晰度评价,分别计算图像的拉普拉斯算子和、灰度梯度,其结果见表2。不含噪声的仿真结果图像的平均拉普拉斯算子和(LS值)为0.15,平均灰度梯度(GMG值)为0.70。含噪声图像的LS值为0.14,GMG值为0.68。仿真图像边缘易检测,高频信息丰富,画质清晰。传感器噪声对图像质量的影响较小。图像仿真时间见表2第6列,平均单帧图像仿真时间为123.45 s,运行速度较快。

      表 2  图像评价指标

      Table 2.  Image evaluation metrics

      Evaluation index
      image in figures
      LSGMGRuntime/s
      Fig.4Fig.5Fig.4Fig.5
      (a)0.0850.0780.340.3398.57
      (b)0.0900.0800.440.4098.71
      (c)0.110.100.500.50164.22
      (d)0.130.110.710.6187.30
      (e)0.130.120.580.53105.50
      (f)0.140.120.630.57111.58
      (g)0.170.150.810.76154.30
      (h)0.160.150.850.77108.92
      (i)0.220.230.981.16156.75
      (j)0.190.170.940.85139.99
      (k)0.190.190.930.99144.22
      (l)0.140.130.680.68111.33
      Average0.150.140.700.68123.45
    • 文中研究了地球临边场景的红外卫星遥感成像仿真方法。首先,根据场景组分的差异,分别建立海洋分布模型、多层云分布模型,并根据海洋和云层的红外辐射和反射特性,建立地球临边场景红外辐射模型。其次,建立了地球-空间坐标系与红外相机坐标系的转换关系,结合大气传输理论和传感器效应,计算各观测角度的地球临边场景红外遥感仿真图像。文中提出了根据云覆盖域的变化计算海洋红外辐射和大气辐射的图像仿真模型,能够反映云区对海洋接收红外辐射造成的削弱效果;根据虚拟目标进行场景定位,能够仿真卫星对不同地理位置、不同飞行高度的空中目标进行遥感观测的情况。仿真实验结果表明,该方法能够仿真不同时间下卫星在各轨道高度探测的地球临边场景,图像平均拉普拉斯算子和可达0.15,平均灰度梯度可达0.70,图像画质清晰;单帧图像平均仿真时间为123 s,可以快速模拟地球临边场景。

参考文献 (14)

目录

    /

    返回文章
    返回