留言板

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

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

星载斜视等距扫描成像

王浩 石志城 宫辉 王军 唐绍凡

王浩, 石志城, 宫辉, 王军, 唐绍凡. 星载斜视等距扫描成像[J]. 红外与激光工程, 2022, 51(3): 20210390. doi: 10.3788/IRLA20210390
引用本文: 王浩, 石志城, 宫辉, 王军, 唐绍凡. 星载斜视等距扫描成像[J]. 红外与激光工程, 2022, 51(3): 20210390. doi: 10.3788/IRLA20210390
Wang Hao, Shi Zhicheng, Gong Hui, Wang Jun, Tang Shaofan. Spaceborne squint isometric scanning imaging[J]. Infrared and Laser Engineering, 2022, 51(3): 20210390. doi: 10.3788/IRLA20210390
Citation: Wang Hao, Shi Zhicheng, Gong Hui, Wang Jun, Tang Shaofan. Spaceborne squint isometric scanning imaging[J]. Infrared and Laser Engineering, 2022, 51(3): 20210390. doi: 10.3788/IRLA20210390

星载斜视等距扫描成像

doi: 10.3788/IRLA20210390
基金项目: 科工局民用航天十四五预研项目
详细信息
    作者简介:

    王浩,男,工程师,硕士,主要从事空间光学遥感总体设计工作

  • 中图分类号: TP732.2

Spaceborne squint isometric scanning imaging

Funds: The 14th Five-Year Civil Aerospace Technology Advance Pre-research Project
  • 摘要: 针对传统穿轨扫描成像方式存在的随扫描角度增大分辨率退化严重、大气程辐射差异大影响定量化应用等问题,提出了一种基于斜视等距扫描的成像方法。首先,介绍了该成像方法的原理,根据几何光学理论建立了几何成像模型。然后,基于几何成像模型,给出了大幅宽斜视等距扫描成像过程中扫描方向与垂直扫描方向空间分辨率随扫描角度的关系式,以及幅宽与扫描角度的关系式;进一步给出了该成像体制,通过卫星平台俯仰机动,实现多角度观测时分辨率、幅宽等关系式。最后,结合某预研星载相机进行了仿真分析。结果表明:分辨率退化大幅减小,在扫描角度60°时边缘分辨率退化相比穿轨扫描成像体制的9.3倍降低为2.5倍。为大幅宽成像边缘视场分辨率退化严重的问题提供了一种新的解决途径,对推动超大幅宽、高分辨率、多角度星载遥感的发展具有一定参考意义。
  • 图  1  穿轨扫描分辨率

    Figure  1.  GSD of cross-track scanning imaging

    图  2  大气程辐射在相机入瞳总能量的占比

    Figure  2.  Proportion of atmospheric radiation in the total energy of imager entrance pupil

    图  3  成像原理示意图

    Figure  3.  Schematic diagram of imaging principle

    图  4  分辨率计算示意图

    Figure  4.  Schematic diagram of resolution calculation

    图  5  幅宽计算示意图

    Figure  5.  Schematic diagram of width calculation

    图  6  分辨率随扫描角度变化曲线

    Figure  6.  Curve of GSD changes with scanning angle

    图  7  归一化分辨率随扫描角度变化曲线

    Figure  7.  Curve of normalized GSD changes with scanning angle

    图  8  幅宽随扫描角度变化曲线

    Figure  8.  Curve of width changes with scanning angle

    图  9  卫星机动和不机动分辨率对比

    Figure  9.  Comparison of GSD between satellite maneuver and non maneuver

    图  10  不同俯仰角度时地面轨迹

    Figure  10.  Ground trajectory at different pitch angles

    图  11  幅宽随俯仰角度变化曲线

    Figure  11.  Curve of width changes with pitch angle

  • [1] Liu Xiaoyong, Cao Kaiqin, Sun Dexin, et al. Scanning control of spaceborne infrared imaging with super-swath and low distortion [J]. Optics and Precision Engineering, 2018, 26(1): 208-217. (in Chinese) doi:  10.3788/OPE.20182601.0208
    [2] Zhao Yanhua, Zhao Limin. Analysis of advantages of multi-angle infrared observation [J]. Spacecraft Recovery & Remote Sensing, 2017, 38(1): 30-37. (in Chinese)
    [3] Zhao Zongcun, Ding Xuezhuan, Yang Bo, et al. Design and analysis of super large width and low distortion imaging system [J]. Journal of Infrared and Millimeter Waves, 2017, 36(6): 732-738. (in Chinese)
    [4] Xiong X, Sun J, Xie X, et al. On-orbit calibration and performance of aqua MODIS reflective solar bands [J]. IEEE Transactions on Geoscience & Remote Sensing, 2009, 48(1): 535-546.
    [5] Yang Z, Lu N, Shi J, et al. Overview of FY-3 payload and ground application system [J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(12): 4846-4853. doi:  10.1109/TGRS.2012.2197826
    [6] Xiong X, Butler J, Chiang K, et al. VIIRS on-orbit calibration methodology and performance [J]. Journal of Geophysical Research Atmospheres, 2014, 119(9): 5065-5078. doi:  10.1002/2013JD020423
    [7] Xiong X, Butler J, Chiang K, et al. VIIRS on-orbit calibration and performance update [C]//Geoscience and Remote Sensing Symposium, IEEE, 2014: 1389-1392.
    [8] Zhang Keli, Zhang Yongchao, Ma Yue. Design of on orbit cross calibration method for HY-1C/D satellite [J]. Spacecraft Engineering, 2019, 28(2): 24-29. (in Chinese)
    [9] Li Xiao, Yu Tingting, Wang Chun. Design and implementation of a new half-angle mirror steering mechanism with high accuracy [J]. Spacecraft Recovery & Remote Sensing, 2019, 40(4): 76-85.
    [10] Cheng Shaoyuan, Zhang Li, Gao Weijun, et al. Geometric parameters analysis of large FOV space camera when rolling [J]. Infrared and Laser Engineering, 2015, 44(6): 198-203. (in Chinese)
    [11] Shi Zhicheng. Research on scan width variation for optical mechanical scanner [J]. Spacecraft Recovery & Remote Sensing, 2011, 32(1): 12-17. (in Chinese)
    [12] Zhang Shuqing, Zhang Yuan, Zhou Chenghao, et al. Image motion model of azimuthally photography for satellite borne TDICCD camera[J]. Infrared and Laser Engineering, 2014, 43(6): 1823-1829. (in Chinese)
    [13] Liu Yinnian. Analysis of the imaging characteristics and scanning traces of the 45° rotating scanning mirror [J]. Optics and Precision Engineering, 2002, 10(1): 110-115. (in Chinese)
  • [1] 党诗沛, 李润泽, 周美玲, 千佳, 但旦, 于湘华, 姚保利.  显微成像系统分辨率问题讨论(特邀) . 红外与激光工程, 2022, 51(11): 20220735-1-20220735-8. doi: 10.3788/IRLA20220735
    [2] 张爱武, 张希珍, 赵江华.  斜模超分辨率成像最佳角度 . 红外与激光工程, 2019, 48(8): 826001-0826001(7). doi: 10.3788/IRLA201948.0826001
    [3] 郜鹏, 温凯, 孙雪莹, 姚保利, 郑娟娟.  定量相位显微中分辨率增强技术综述 . 红外与激光工程, 2019, 48(6): 603007-0603007(13). doi: 10.3788/IRLA201948.0603007
    [4] 余乐文, 张达.  采空区三维激光扫描空间分辨率增强方法 . 红外与激光工程, 2017, 46(10): 1006002-1006002(6). doi: 10.3788/IRLA201760.1006002
    [5] 刘波, 耿林, 刘琳, 曹昌东, 眭晓林, 颜子恒.  双通道接收的无扫描激光成像技术 . 红外与激光工程, 2016, 45(12): 1206008-1206008(5). doi: 10.3788/IRLA201645.1206008
    [6] 赵明, 谢晨波, 钟志庆, 王邦新, 王珍珠, 尚震, 谭敏, 刘东, 王英俭.  高光谱分辨率激光雷达探测大气透过率 . 红外与激光工程, 2016, 45(S1): 76-80. doi: 10.3788/IRLA201645.S130002
    [7] 张雷, 丁亚林, 徐正平, 张洪文, 张健, 郭万存.  长条形扫描反射镜的柔性支撑 . 红外与激光工程, 2015, 44(12): 3678-3683.
    [8] 巩晋南, 侯晴宇, 邓盼, 张伟.  扫描非均匀性对点目标探测影响 . 红外与激光工程, 2015, 44(11): 3282-3287.
    [9] 沈同圣, 史浩然, 娄树理, 李召龙.  填充因子对微扫描系统成像质量影响分析 . 红外与激光工程, 2015, 44(8): 2292-2297.
    [10] 李孟麟, 朱精果, 孟柘, 任建峰, 肖芳, 张珂殊.  轻小型机载激光扫描仪设计 . 红外与激光工程, 2015, 44(5): 1426-1431.
    [11] 张明凯, 高思田, 卢荣胜, 李伟, 李琪, 钱晓力.  紫外扫描线宽测量系统的研究 . 红外与激光工程, 2015, 44(2): 625-631.
    [12] 孙鸣捷, 于康龙.  像素有效形状对微扫描红外超分辨成像的影响 . 红外与激光工程, 2015, 44(1): 48-52.
    [13] 田祥瑞, 徐立军, 徐腾, 李小路, 张勤拓.  车载LiDAR扫描系统安置误差角检校 . 红外与激光工程, 2014, 43(10): 3292-3297.
    [14] 张树青, 张媛, 周程灏, 王治乐.  星载TDICCD 相机方位扫描像移模型研究 . 红外与激光工程, 2014, 43(6): 1823-1829.
    [15] 李鹏, 张常麟, 王文学, 刘连庆.  扫描离子电导显微镜的研制与实现 . 红外与激光工程, 2014, 43(6): 1894-1898.
    [16] 孙鸣捷, 于康龙, 孟照魁.  像素感光面形状对基于微扫描的红外超分辨成像的影响 . 红外与激光工程, 2014, 43(10): 3183-3187.
    [17] 林建粦, 平西建, 马德宝.  采用DBT的漂移扫描星图小目标检测方法 . 红外与激光工程, 2013, 42(12): 3440-3446.
    [18] 李海星, 姚园, 史磊.  单向驱动往复式扫描反射镜模型 . 红外与激光工程, 2013, 42(8): 2097-2102.
    [19] 崔智高, 李艾华, 姜柯, 苏延召, 金广智.  双目协同多分辨率主动跟踪方法 . 红外与激光工程, 2013, 42(12): 3509-3516.
    [20] 李珂, 张康伟, 罗淼.  基于LCTF的大幅面高分辨率多光谱仪光学系统设计 . 红外与激光工程, 2013, 42(3): 675-679.
  • 加载中
图(11)
计量
  • 文章访问数:  67
  • HTML全文浏览量:  15
  • PDF下载量:  19
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-06-09
  • 修回日期:  2021-09-28
  • 刊出日期:  2022-04-07

星载斜视等距扫描成像

doi: 10.3788/IRLA20210390
    作者简介:

    王浩,男,工程师,硕士,主要从事空间光学遥感总体设计工作

基金项目:  科工局民用航天十四五预研项目
  • 中图分类号: TP732.2

摘要: 针对传统穿轨扫描成像方式存在的随扫描角度增大分辨率退化严重、大气程辐射差异大影响定量化应用等问题,提出了一种基于斜视等距扫描的成像方法。首先,介绍了该成像方法的原理,根据几何光学理论建立了几何成像模型。然后,基于几何成像模型,给出了大幅宽斜视等距扫描成像过程中扫描方向与垂直扫描方向空间分辨率随扫描角度的关系式,以及幅宽与扫描角度的关系式;进一步给出了该成像体制,通过卫星平台俯仰机动,实现多角度观测时分辨率、幅宽等关系式。最后,结合某预研星载相机进行了仿真分析。结果表明:分辨率退化大幅减小,在扫描角度60°时边缘分辨率退化相比穿轨扫描成像体制的9.3倍降低为2.5倍。为大幅宽成像边缘视场分辨率退化严重的问题提供了一种新的解决途径,对推动超大幅宽、高分辨率、多角度星载遥感的发展具有一定参考意义。

English Abstract

    • 大幅宽成像可以缩短卫星的重访周期,提高观测的时间分辨率;高分辨遥感图像可以提供更精细的几何特征信息,提高目标识别能力[1];多角度遥感通过地面目标多个方向的观察,使得对目标的观测信息更加丰富,为定量遥感提供了多元的信息,有利于提高遥感数据定量化水平[2];大幅宽、高分辨率、多角度三者结合,将会是未来对地遥感的重要发展方向。

      目前超大幅宽对地成像方式主要有多载荷视场拼接、高轨成像和光机扫描成像。其中,光机扫描成像方式具有机动、灵活、容易实现且成本低的特点,在机载、航天大幅宽成像观测中得到广泛应用[3],如:美国MODIS[4]、我国风云三号中分辨率成像光谱仪采用反射镜摆扫的光机扫描方案[5],美国VIIRS、我国海洋一号C/D星上定标光谱仪以及正在研制的新一代海洋水温水色仪采用主光学望远扫描的光机扫描方案[6-9]。但超大幅宽成像时,传统光机扫描成像方式受观测斜距、地球曲率的影响,使得地面空间分辨率随扫描角度增大而迅速降低,降低视场边缘目标可探性的;同时带来观测天顶角逐渐增大,成像幅宽范围内大气斜程路径差异大,影响定量化应用的问题。

      文中提出了一种基于斜视等距扫描的成像方法,旨在解决大幅宽光机扫描成像过程中视场边缘地面空间分辨率急剧降低、不同像元大气程辐射差异大而影响定量化应用的问题。推导了大幅宽斜视等距扫描成像过程中扫描方向与垂直扫描方向(线阵方向)空间分辨率随扫描角度的关系式,幅宽随扫描角度的关系式;进一步给出了该成像体制,通过卫星平台俯仰机动实现多角度观测时分辨率、幅宽的数学表达式,并进行了仿真分析。为大幅宽成像提供了一种新的解决途径,对推动超大幅宽、高分辨率、多角度遥感的发展具有一定参考意义。

    • 由于受观测斜距、地球曲率的双重影响,穿轨光机扫描成像方式空间分辨率随扫描角度增大而迅速降低,导致地面轨迹存在明显的蝴蝶结效应[10-11]。以卫星运行轨道高度705 km、相机焦距4250 mm、探测器像元尺寸10 μm、扫描角度范围$ \left[ { - {{60}^ \circ },{{60}^ \circ }} \right] $为例,图1给出了中心视场像元沿轨$ GS{D_x} $及穿轨$ GS{D_y} $两个方向分辨率的退化情况,可以看出沿轨与穿轨方向地面分辨率从星下点1.66 m分别退化到扫描边缘4.25、15.44 m,分辨率退化严重。图2给出了1976年美国标准大气、乡村可视距离5 km气象条件下,该相机从星下点到视场边缘成像,可见到长波红外全谱段的大气程辐射在相机入瞳总能量的占比随观测天顶角的变化趋势,星下点观测天顶角为0°时大气程辐射占比为42.2%,扫描边缘观测天顶角为60°,大气程辐射增加到56.32%,大气程辐射差异大。为此,文中提出了一种斜视等距扫描成像方案用以解决上述问题。

      图  1  穿轨扫描分辨率

      Figure 1.  GSD of cross-track scanning imaging

      图  2  大气程辐射在相机入瞳总能量的占比

      Figure 2.  Proportion of atmospheric radiation in the total energy of imager entrance pupil

    • 图3所示,建立坐标系,X方向为卫星平台飞行方向,Z方向指向地心,Y方向满足右手坐标系扫描镜与XOY平面夹角为$ \theta $,扫描镜旋转轴为P,卫星不姿态机动时,P轴与Z方向平行。光学系统光轴方向与Y轴平行,探测器位置固定,线阵方向平行于X轴,垂直线阵方向平行于Y轴。扫描成像过程中,探测器任一像元$ i $扫描的物方矢量为圆锥形,对应的探测斜距相同,地面投影扫描轨迹为环带,中心视场像元对应的观测天顶角为$ 2\theta $。因此,可以消除传统穿轨光机扫描成像时探测斜距对分辨率的影响。此外,采用斜视等距观测时大气路径相同,大幅宽成像时不同视场的大气程辐射接近,利于降低大气程辐射的影响及载荷的定量化应用。

      图  3  成像原理示意图

      Figure 3.  Schematic diagram of imaging principle

    • 图1所示,任一像元$ i $的像方矢量用矩形四个角点的单位化矢量表示:${{\boldsymbol{r}}_I} = \left[ {{{\boldsymbol{r}}_a}\left( i \right),{{\boldsymbol{r}}_b}\left( i \right),{{\boldsymbol{r}}_c}\left( i \right),{{\boldsymbol{r}}_d}\left( i \right)} \right]$。根据光学反射矢量基本理论,对应的单位化物方矢量${{\boldsymbol{r}}_E}\left( {i,\;\beta } \right) = \left[ {{{\boldsymbol{r}}_{a'}}\left( {i,\;\beta } \right),{{\boldsymbol{r}}_{b'}}\left( {i,\;\beta } \right),{{\boldsymbol{r}}_{c'}}\left( {i,\;\beta } \right),{{\boldsymbol{r}}_{d'}}\left( {i,\;\beta } \right)} \right]$可以由像方矢量$ {{\boldsymbol{r}}_I} $及法线向量${{\boldsymbol{N}}_{\boldsymbol{\;\beta }}} = {[{N_{\;\beta ,x}},{N_{\;\beta ,y}},{N_{\;\beta ,z}}]^{\rm{T}}}$决定的反射矩阵$ {R_N}\left( \;\beta \right) $线性表示[12]

      $$ \begin{split}& {{\boldsymbol{r}}_E}\left( {i,\;\beta } \right) = {R_N}\left( \;\beta \right){{\boldsymbol{r}}_I}\left( i \right) = \\& \left[ \begin{gathered} 1 - 2N_{\;\beta ,x}^2{\text{ }} - 2{N_{\;\beta ,x}}{N_{\;\beta ,y}}{\text{ }} - 2{N_{\;\beta ,x}}{N_{\;\beta ,{\textit{z}}}} \\ - 2{N_{\;\beta ,x}}{N_{\;\beta ,y}}{\text{ }}1 - 2N_{\;\beta ,y}^2{\text{ }} - 2{N_{\;\beta ,y}}{N_{\;\beta ,{\textit{z}}}} \hfill \\ - 2{N_{\;\beta ,x}}{N_{\;\beta ,{\textit{z}}}}{\text{ }} - 2{N_{\;\beta ,y}}{N_{\;\beta ,{\textit{z}}}}{\text{ }}1 - 2N_{\;\beta ,{\textit{z}}}^2 \hfill \\ \end{gathered} \right]{{\boldsymbol{r}}_I}\left( i \right) \end{split} $$ (1)

      当扫描角度$ \;\;\beta = {0^ \circ } $时,扫描镜法线与X轴、Z轴、探测器线阵在同一平面内,初始法线矢量为${{\boldsymbol{N}}_{0}} = {\left[ {\sin \theta ,0,\cos\theta } \right]^{\rm{T}}}$。当扫描镜绕P轴转动时,法线矢量为[13]

      $$ {{\boldsymbol{N}}_{\boldsymbol{\;\beta }}} = {\boldsymbol{S}_{P,\;\;\beta }}{\left[ {\sin \theta ,0,\cos\theta } \right]^{\rm{T}}} $$ (2)

      式中:${\boldsymbol{S}_{P,\;\;\beta }}$为转动矩阵。

      $$ {\boldsymbol{S}_{P,\;\;\beta }} = \left[ \begin{gathered} \cos \;\beta + 2{P_x}^2{\sin ^2}\;\beta /2{\text{ }} - {P_{\textit{z}}}\sin \;\beta + 2{P_x}{P_y}{\sin ^2}\;\beta /2{\text{ }}{P_y}\sin \;\beta + 2{P_x}{P_{\textit{z}}}{\sin ^2}\;\beta /2 \hfill \\ {P_{\textit{z}}}\sin \;\beta + 2{P_x}{P_y}{\sin ^2}\;\beta /2{\text{ }}\cos \;\beta + 2{P_y}^2{\sin ^2}\;\beta /2{\text{ }} - {P_x}\sin \;\beta + 2{P_y}{P_{\textit{z}}}{\sin ^2}\;\beta /2 \hfill \\ - {P_y}\sin \;\beta + 2{P_x}{P_{\textit{z}}}{\sin ^2}\;\beta /2{\text{ }}{P_x}\sin \;\beta + 2{P_y}{P_{\textit{z}}}{\sin ^2}\;\beta /2{\text{ }}\cos \;\beta + 2{P_{\textit{z}}}^2{\sin ^2}\;\beta /2 \hfill \\ \end{gathered} \right] $$ (3)

      矢量${\boldsymbol{P}}{\text{ = }}{\left[ {{P_x},{P_y},{P_{\textit{z}}}} \right]^{\rm{T}}}$P轴的单位方向矢量。卫星不姿态机动时,P轴与Z方向平行。

    • 图4所示,相机探测器任一像元$ i $的四个物方矢量${{\boldsymbol{r}}_E}\left( {i,\;\;\beta } \right)$在地面的投影点为$ a',b',c',d' $,则投影点指向地心的向量$ {{\boldsymbol{s}}_{a'}} $可表示为:

      图  4  分辨率计算示意图

      Figure 4.  Schematic diagram of resolution calculation

      $$ {{\boldsymbol{s}}_{a'}} = {L_{a'}}\cdot {{\boldsymbol{r}}_{a'}}\left( {i,\;\beta } \right) + (R + H){\left[ {0,0,1} \right]^{\rm{T}}} $$ (4)

      式中:$ {L_{a'}} $为视线斜距。

      $$ {L_{a'}} = (R + H)\cos {\;\beta _{_{a'}}} - \sqrt {{R^2} - {{(R + H)}^2}{{\sin }^2}{\;\beta _{a'}}} $$ (5)

      式中:R为地球半径;H为轨道高度;$ {\;\beta _{_{a'}}} $为物方视线矢量$ {{\boldsymbol{r}}_{{\boldsymbol{a}}'}} $与星下点Z方向的夹角。需要说明的是:只有中心视场对应像元的$ {\;\beta _{_{a'}}} $等于扫描角度$ \;\beta $,其他像元取值虽然与扫描角度$ \;\beta $相关,但并不等同。

      同理可得$ {{\boldsymbol{s}}_b} $$ {{\boldsymbol{s}}_c} $$ {{\boldsymbol{s}}_d} $,则相机探测器线阵方向分辨率及垂直线阵方向分辨率分别为:

      $$ \begin{gathered} GS{D_x} = R\cdot \arccos \left(\frac{{{{\boldsymbol{s}}_{a'}}\left( {i,\;\beta } \right)\cdot {{\boldsymbol{s}}_{b'}}\left( {i,\;\beta } \right)}}{{\left| {{{\boldsymbol{s}}_{a'}}\left( {i,\;\beta } \right)} \right|\cdot \left| {{{\boldsymbol{s}}_{b'}}\left( {i,\;\beta } \right)} \right|}}\right) \hfill \\ GS{D_y} = R\cdot \arccos \left(\frac{{{{\boldsymbol{s}}_{a'}}\left( {i,\;\beta } \right)\cdot {{\boldsymbol{s}}_{d'}}\left( {i,\;\beta } \right)}}{{\left| {{{\boldsymbol{s}}_{a'}}\left( {i,\;\beta } \right)} \right|\cdot \left| {{{\boldsymbol{s}}_{d'}}\left( {i,\;\beta } \right)} \right|}}\right) \hfill \\ \end{gathered} $$ (6)

      图5所示,斜视等距扫描成像体制的有效成像幅宽为扫描环带在地表的弦长。设${{\boldsymbol{r}}_{{{{w}}_s}}}$${{\boldsymbol{r}}_{{{{w}}_e}}}$分别为探测器线列中心视场像元(任取一个角点)扫描对应的物方起始向量与终止向量,可用公式(1)~(3)联合求得。${{\boldsymbol{r}}_{{{{w}}_s}}}$${{\boldsymbol{r}}_{{{{w}}_e}}}$对应的地面扫描环带起始点与终止点分别为$ {w_s} $$ {w_e} $,对应的指向地心向量为${{\boldsymbol{s}}_{{{{w}}_s}}}$${{\boldsymbol{s}}_{{{{w}}_e}}}$,其计算过程同公式(4),则幅宽公式如下:

      图  5  幅宽计算示意图

      Figure 5.  Schematic diagram of width calculation

      $$ W = R\cdot \arccos \left(\frac{{{{\boldsymbol{s}}_{{{{w}}_s}}}\cdot {{\boldsymbol{s}}_{{{{w}}_e}}}}}{{\left| {{{\boldsymbol{s}}_{{{{w}}_s}}}} \right|\cdot \left| {{{\boldsymbol{s}}_{{{{w}}_e}}}} \right|}}\right) $$ (7)

      下面进行卫星运行轨道高度705 km、相机焦距4250 mm、探测器像元尺寸10 μm、$ \theta $取27.5°进行相机斜视等距扫描成像时,分辨率和幅宽随扫描角度的变化分析。图6给出了线阵方向分辨率$ GS{D_x} $及垂直线阵方向分辨率$ GS{D_y} $的变化曲线,可见分辨率随扫描角度波动明显优于图2所示的穿轨扫描成像方式。为与传统穿轨扫描方式进行定量对比,采用各自扫描角度为0°时的分辨率进行归一化处理,结果如图7所示,可见线阵方向分辨率$ GS{D_x} $两种扫描方式差距不大,但穿轨扫描方式垂直线阵方向分辨率$ GS{D_y} $随扫描角度增大急剧增加,扫描角度为60°时的分辨率为星下点(扫描角度0°)时的9.3倍,而文中斜视等距扫描仅为2.5倍。

      图  6  分辨率随扫描角度变化曲线

      Figure 6.  Curve of GSD changes with scanning angle

      图  7  归一化分辨率随扫描角度变化曲线

      Figure 7.  Curve of normalized GSD changes with scanning angle

      图8给出了扫描幅宽随扫描角度的变化趋势,在扫描范围为120°($\; \;\beta = \pm {60^ \circ } $)时,可实现幅宽≥2000 km。

      图  8  幅宽随扫描角度变化曲线

      Figure 8.  Curve of width changes with scanning angle

    • 相机通过卫星俯仰机动实现多角度成像,俯仰机动为绕Y轴逆时针旋转$ \gamma $,则探测器任一像元像方矢量、初始法线矢量、扫描镜旋转轴向量分别转化为:

      $$ {{\boldsymbol{r}}_{I,\rho }} = {R_y}(\rho ){{\boldsymbol{r}}_I} $$ (8)
      $$ {N}_{0\text{,}\rho }={R}_{y}(\rho ){N}_{0} $$ (9)
      $$ {{\boldsymbol{P}}_\rho }{\text{ = }}\left[ {{P_{x,\;\rho }},{P_{y,\;\rho }},{P_{{\textit{z}},\;\rho }}} \right] = {R_y}(\rho ){\boldsymbol{P}} $$ (10)

      将公式(10)中$\left[ {{P_{x,\;\rho }},{P_{y,\;\rho }},{P_{{\textit{z}},\;\rho }}} \right]$代入公式(3)中得到转动矩阵为${{\boldsymbol{S}}_{P,\;\;\beta ,\;\rho }}$,进而得到卫星机动后,扫描过程中的法线矢量为${{\boldsymbol{N}}}_{\;\beta ,\;\rho }={{\boldsymbol{S}}}_{{P}_{\gamma },\;\;\beta ,\;\rho }{N}_{0,\;\rho }$,对应的反射矩阵为${{\boldsymbol{R}}_N}\left( {\;\beta ,\rho } \right)$,结合像方矢量得到物方矢量为:

      $$ {{\boldsymbol{r}}_E}\left( {i,\;\;\beta ,\;\rho } \right) = {{\boldsymbol{R}}_N}\left( {\;\beta ,\;\rho } \right){R_y}(\rho ){{\boldsymbol{r}}_I} $$ (11)

      将公式(11)的计算结果代入公式(4)~(6)中可以计算得到卫星机动后的分辨率;代入公式(4)、(5)、(7)中可以计算得到卫星机动后的幅宽。下面同样以卫星运行轨道高度705 km、相机焦距4250 mm、探测器像元尺寸10 μm、$ \theta = {27.5^ \circ } $进行卫星俯仰机动后的分辨率、幅宽等几何参数分析。图9所示为卫星逆时针俯仰55°和卫星不机动时两个方向分辨率的对比,可见由于卫星机动导致相机视线向星下点方向移动,斜距缩短,分辨率提高。图10图11所示分别为给出相同扫描角度$\; \;\beta = \left[ { - {{60}^ \circ },{{60}^ \circ }} \right] $时,卫星绕Y轴逆时针俯仰机动$\; \rho $角度后的地面轨迹及幅宽,可见随着卫星俯仰机动角度的增加,幅宽先减小后增大,在$\; \rho = - {36^ \circ } $时,幅宽最小为1510 km。

      图  9  卫星机动和不机动分辨率对比

      Figure 9.  Comparison of GSD between satellite maneuver and non maneuver

      图  10  不同俯仰角度时地面轨迹

      Figure 10.  Ground trajectory at different pitch angles

      图  11  幅宽随俯仰角度变化曲线

      Figure 11.  Curve of width changes with pitch angle

    • 大幅宽、高分辨率、多角度是未来对地遥感的重要发展方向,针对目前穿轨扫描成像体制存在的分辨率退化严重、大气程辐射差异大影响定量化应用等问题,提出了一种基于斜视等距扫描的成像方法,对大幅宽、高分辨率、多角度遥感探测的发展具有参考意义。给出了:(1) 大幅宽斜视等距扫描成像过程中,扫描方向与垂直扫描方向空间分辨率随扫描角度的关系式;(2) 幅宽随扫描角度的关系式;(3) 通过卫星平台俯仰机动,实现多角度观测时分辨率、幅宽关系式。将文中方法应用于某预研星载大视场相机,视场边缘分辨率退化相比穿轨扫描成像体制的9.3倍降低为2.5倍。

      同时,注意到斜视等距扫描方式在两个方向分辨率变化规律不同,这主要是由扫描镜扫描引入像旋导致的,像旋会影响遥感图像的解译应用,如何针对这种斜视等距扫描体制进行消旋是载荷设计及影像应用的重点,后续将持续深入研究。

参考文献 (13)

目录

    /

    返回文章
    返回