1、第 52 卷第 2 期2023 年 2 月Vol.52 No.2February 2023光子学报ACTA PHOTONICA SINICA02110021星载的单像素运动成像影响分析与补偿于跃1，杨照华1，余远金2（1 北京航空航天大学 仪器科学与光电工程学院，北京 100191）（2 北京理工大学 自动化学院，北京 100081）摘要：单像素成像具有高灵敏度，抗干扰等特点，但得到一张高质量图像需要多次采样，因而在卫星对地观测应用时，由于卫星姿态变化使得指向抖动从而造成参考光场与探测器信号之间关联性丢失，使得重构出的图像质量严重退化。针对这一问题，本文设计一种光场追踪补偿方案以补偿卫星姿态对
2、成像带来的影响，结合卫星工具包生成卫星数据分别进行滚动角，俯仰角的运动补偿验证和偏航角的运动补偿验证。结果表明进行滚动角和俯仰角运动补偿后成像的峰值信噪比提升 28.8 dB，进行偏航角运动补偿后成像的峰值信噪比提升 7.25 dB，同时考虑滚动角，俯仰角和偏航角对典型场景进行运动补偿后的成像结果峰值信噪比提升至少 4.7 dB，仿真结果证明了补偿方案有效提高了成像质量，为星载单像素成像提供了一条有效技术途径。关键词：单像素成像；图像重构；差分关联成像；星载成像；追踪补偿中图分类号：O436 文献标识码：A doi：10.3788/gzxb20235202.02110020 引言单像素成像1又
3、称鬼成像，具有非视域、超分辨、抗干扰等特点，自二十世纪九十年代初被提出以来，引起了广泛的关注。随着工程化和实用化研究的不断深入，单像素成像技术在光学遥感、医学显微、视觉传感等多领域发挥其独特优势，同时，在国防军事及生命科学领域具有极高的潜在应用价值。单像素成像技术最早可追溯到利用光场强度的二阶关联特性完成恒星角距和恒星半径的测量。在单像素成像研究的初始阶段，认为纠缠光子对2是单像素成像的必要条件，后来有学者质疑这一条件，并通过实验验证最终证实赝热光甚至热光3也可以实现单像素成像。纠缠光单像素成像中存在光源制备困难的问题，热光、赝热光双光臂单像素成像实验中必须满足物臂和参考臂等距，在实际应用中很
4、难满足这一要求，因此有学者提出计算关联成像4-5，实现单光臂关联成像。此后，国内外在单像素成像领域的研究发展迅速，在纠缠光成像、赝热光成像、计算量子关联成像等方面取得了极大进展6-13，这些研究成果将推动和促进单像素成像技术的应用。单像素成像包括主动成像（或称基于结构化照明）和被动成像（或称基于结构化探测）两种实现形式，其中利用激光作为光源的主动成像方式在远距离成像探测中更有价值14。光场强度调制的质量和速度是影响成像的关键。而且需要调制后照射到物体的光场和收集到的光强保持精准的同步，否则物体与调制光的时空关联性丢失，将难以重构出清晰的像，甚至导致成像失败。研究运动对单像素成像的影响对于远距离
5、遥感、目标探测与识别、生物医疗诊断等具有极大的现实意义和潜在的应用价值，基于此本文研究星载应用下的单像素成像方案，考虑卫星姿态角变化，补偿卫星指向抖动对单像素成像的影响，实现相对运动下的目引用格式：YU Yue，YANG Zhaohua，YU Yuanjin.Motion Analysis and Compensation of Space-borne Single Pixel ImagingJ.Acta Photonica Sinica，2023，52（2）：0211002于跃，杨照华，余远金.星载的单像素运动成像影响分析与补偿 J.光子学报，2023，52（2）：0211002基金项目：国
6、 家 自 然 科 学 基 金(Nos.61973018,61975229,62173039)，国 家 民 用 航 天 计 划(No.D30401)，国 防 基 础 科 研 计 划(No.JCKY2021602B036)第一作者：于跃，yuyue_通讯作者：杨照华，；余远金，收稿日期：2022 09 09；录用日期：2022 11 03http:/光子学报02110022标图像高质量重构。1 单像素成像理论基础1.1单像素成像光路如图 1所示，计算关联成像方案采用可人为设计参考光场的空间光调制器件进行光场调制，代替传统关联成像的激光照射毛玻璃方案产生参考散斑光场，从而可只利用物臂来实现对物体图像
7、的关联重构。其最大特点在于单光臂测量操作简便、约束较少，从成像光路结构上，更具有实用性。综合以上优点，实际工程广泛采用计算关联成像方案。1.2单像素成像基本原理在计算关联成像方案中，空间光调制器件对二维光场进行时空强度调制，t时刻其光场强度分布表示为It(x，y)，同时It(x，y)也是 t时刻空间光调制器件中所被加载的调制矩阵，其中(x，y)表示光场的二维空间位置坐标。空间光调制器、成像透镜和物体三者之间满足薄透镜成像高斯公式1L1+1L2=1f（1）式中，L1为空间光调制器到成像透镜之间的距离，L2为成像透镜到物体之间的距离，f为成像透镜的焦距。假设物体的空间透过率函数用T(x，y)表示，
8、则 t时刻单像素桶探测器的光强响应信号为Bt=It(x，y)T(x，y)dxdy（2）根据二阶关联成像理论，经过N次测量后的强度关联函数为G(x，y)=1Nr=1N()Br-BIr(x，y)=BrIr(x，y)-BI(x，y)（3）式中，B=1Nr=1NBr代表数据的系综平均。整个测量过程需要参考臂记录散斑光场分布用以重构图像，而重构出的图像中每个像素点的强度值是进行系综平均运算后的计算结果，也就是光场分布对应光强信号相乘相加最后求平均的过程。差分关联成像（Differential Ghost Imaging，DGI）多用面阵单元对参考光场进行记录，其后端往往要 伴 随计算单元进行累加求和操作
9、，以获取差分信号，相较于直接关联算法成像可以一定程度上消除图像的噪声。物体的瞬时透过率函数为T=AsI(x)T(x)d2xAsI(x)d2x（4）图 1计算关联成像光路图Fig.1Computational correlation imaging scheme于跃，等：星载的单像素运动成像影响分析与补偿02110023此时，经过差分计算后的光强信号S-应为S-=AsIB(xB)T(xB)d2x（5）式中，T(x)=T(x)-T利用S-代替原有的桶探测器信号SB，可得 DGI的数学表示为G2(x)DGI=S-IR(x)=SBIR(x)-SBSRSBIR(x)（6）1.3评价指标峰值信噪比（Pea
10、k Signal-to-Noise Ratio，PSNR）是一个信号的最大可能功率与影响它的表示精度的破坏性噪声功率的比值，是经常用于图像压缩等领域中信号重建质量的测量方法。峰值信噪比定义为PSNR=10 lg()MAX2i1MNi=1Mj=1NT0(i，j)-T(i，j)2（7）通过计算原图和关联重构后的成像结果间的误差关系来评价图像还原程度，不仅适合传统图像，而且完全符合单像素成像的质量评价标准，本文后续内容均以此为基础进行图像重构质量评估。2 星载的单像素运动成像影响分析静态目标单像素成像可以保证桶探测器所探测到的总光强与照明散斑之间的关联性，由此进行关联重构可以重构出高质量图像。然而进
11、行星载单像素成像时，由于卫星本身存在指向抖动，使卫星的滚动角，俯仰角和偏航角时刻发生改变，因此在星载单像素成像过程中散斑光场和成像物体之间发生偏移，当单像素成像过程中卫星姿态角变化大于单像素成像的角分辨率时，桶探测器所探测到的总光强与照明散斑之间的关联性丢失，直接进行关联重构所重构出的图像质量会大大降低。如图 2 所示，令星载单像素成像的视场角为，成像数字分辨率为N N，地心惯性坐标系表示为oxiyizi，卫星轨道坐标系表示为oxoyozo，此时卫星的姿态角表示为滚动角，俯仰角和偏航角。星载单像素成像可直接重构的条件为图 2卫星探测地球示意图Fig.2Diagram of satellite
12、exploration of the earth光子学报02110024 tNtNtN（8）式中，t是卫星滚动角，t是卫星俯仰角，t是卫星偏航角。当卫星姿态角变化不符合直接单像素成像条件时，可以通过恢复桶探测器所探测到的探测信号与散斑光场之间的关联性以恢复物体高质量图像，根据卫星姿态运动设计光场追踪补偿的单像素成像方法。假设卫星位于地球静止轨道，且不存在指向抖动，则探测器收集到的光强值应为Si=Ii(x，y)T(x，y)dxdy（9）式中，Si是第 i次采样的桶探测器值，Ii(x，y)为由哈达玛矩阵第 i行或第 i列变形成为的第 i个光场序列，T(x，y)代表物体的透过率函数。实际星载探测过程
13、中，卫星的滚动角、俯仰角和偏航角均会因指向抖动而改变，这导致每次采样时物体的有效信息不断改变，当仅考虑滚动角和俯仰角的改变时，桶探测器的探测值可以表示为Sip=Ii(x，y)T()x+(i)，y+(i)dxdy（10）式中，(i)是第 i次采样过程中的卫星滚动角，(i)是第 i次采样过程中的卫星俯仰角。当仅考虑偏航角的改变时，桶探测器的探测值可以表示为Sir=020rIi(，)T()，+(i)dd（11）式中，(i)是第 i次采样过程中的卫星偏航角。由于卫星滚动角、俯仰角和偏航角的变化，导致桶探测器的探测值和经空间光调制器调制的参考光场失去关联性，这是造成重构图像质量下降的根本原因。因此，星载
14、的单像素运动成像的关键是对卫星的姿态运动进行补偿，恢复探测值与参考光场的关联性。3 基于卫星姿态信息的星载单像素运动成像补偿针对星载的单像素运动成像采用一种光场追踪补偿单像素成像方案，主要思路是：桶探测器探测的光强值是一组数值，当卫星指向发生抖动时，若将可控的散斑光场按照卫星的指向运动方式进行对应的运动补偿，可以使得每次对目标物体采样时桶探测器收集到的光强值与对静止物体采样获得的探测结果大致相同，再进行二阶关联运算就可以有效提升星载成像质量。基于上述思想，在计算关联成像重构图像的过程中，将人为设定的散斑光场Ii(x，y)按照卫星的指向运动方式进行对应补偿，就可以得到同样的强度信号序列Sip=I
15、i()x+(i)，y+(i)T()x+(i)，y+(i)d()x+(i)d()y+(i)=Ii()xp，ypT()xp，ypdxpdyp （12）Sir=020rIi()，+(i)T()，+(i)dd()+(i)=020rIi(，r)T(，r)ddr（13）式中，xp=x+(i)，yp=y+(i)，r=+(i)。式（12）和式（10），式（13）和式（11）在实施光场追踪补偿方案后数学表达形式均相同。滚动角和俯仰角运动的光场追踪补偿原理图如图 3所示，图 3（a）为静态物体，图 3（b）为静态物体对应参考光场，图 3（c）为卫星姿态变化过程中某一时刻物体所处位置，图 3（d）为与姿态变化对应的
16、补偿参考光场。目标物体相对于卫星发生(，)的指向抖动，使得物体在视场中的位置由图 3（a）位置变化到图 3（c）位置，将可控的散斑光场根据卫星姿态运动进行对应的运动补偿，可以获取同静止状态相同的桶探测器光强信号。图 4 用来描述偏航角运动的光场追踪补偿方案，其中图 4（a）为静态物体，图 4（b）为静态物体对应参考光场，图 4（c）为卫星姿态变化过程中某一时刻物体所处位置，图 4（d）为与姿态变化对应的补偿参考光场。于跃，等：星载的单像素运动成像影响分析与补偿021100254 基于 STK 的星载单像素运动成像仿真4.1卫星数据生成卫星工具软件包（Satellite Tool Kit，STK）可生成航天器的轨道和姿态数据，是卫星计算机仿真分析中常用的工具。利用 STK 仿真地球静止轨道卫星的飞行，所生成的卫星轨道如图 5所示。卫星轨道六要素分别为半长轴（Semi-major Axis）42 152.139 604 km，偏心率（eccentricity）0.000 396，轨道倾角（inclination 图 3滚动角和俯仰角指向抖动光场追踪补偿原理Fig.3Schematic di