星载的单像素运动成像影响分析与补偿_于跃.pdf
《星载的单像素运动成像影响分析与补偿_于跃.pdf》由会员分享,可在线阅读,更多相关《星载的单像素运动成像影响分析与补偿_于跃.pdf(11页珍藏版)》请在咨信网上搜索。
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
17、)可生成航天器的轨道和姿态数据,是卫星计算机仿真分析中常用的工具。利用 STK 仿真地球静止轨道卫星的飞行,所生成的卫星轨道如图 5所示。卫星轨道六要素分别为半长轴(Semi-major Axis)42 152.139 604 km,偏心率(eccentricity)0.000 396,轨道倾角(inclination 图 3滚动角和俯仰角指向抖动光场追踪补偿原理Fig.3Schematic diagram of light field tracking compensation for roll angle and pitch angle图 4偏航角指向抖动光场追踪补偿原理图Fig.4Sch
18、ematic diagram of light field tracking compensation for yaw angle光子学报02110026angle)0.464,右 旋 升 交 点 赤 经(the right ascension of ascending node)12.451,真 近 点 角(true anomaly)322.316,近地点辐角(argument of perigee)34.529。利用卫星工具包输出卫星滚动角、俯仰角和偏航角的数据。4.2基于 STK 卫星数据的滚动角与俯仰角运动补偿为了验证星载的单像素运动补偿方案的可行性,针对卫星滚动角和俯仰角的变化进行如
19、下仿真实验,运动目标为小飞机,运动场景大小为 128128,空间光调制器频率为 1 kHz。由于卫星滚动角和俯仰角不断改变,小飞机位置对应每次采样过程其相应的数学表达为xp=x+(i)yp=y+(i)(14)式中,(i)为每次采样过程中的滚动角,(i)为每次采样过程中的俯仰角。(i)和(i)均随采样数 i改变,变化特性如图 6所示。进行星载滚动角和俯仰角运动补偿单像素成像仿真,对散斑光场区域投影 16 384张空间及时间上变化的散斑光场,对应仿真过程是 Hadamard矩阵与运动物体所处区域的内积运算,重构结果如图 7所示,分别计算峰值信噪比来比较运动补偿前后图像还原质量,其中,图 7(c)为
20、运动补偿前的重构结果,峰值信噪比为9.68 dB,图 7(d)为运动补偿后重构结果,峰值信噪比为 38.48 dB。仿真结果表明,光场追踪补偿策略适用于基于 STK卫星数据的滚动角与俯仰角运动补偿,可以有效提升图像重构质量。图 5卫星轨道Fig.5Satellite orbit图 6滚动角和俯仰角随采样数的变化特性Fig.6Characteristics of roll angle and pitch angle with sampling number于跃,等:星载的单像素运动成像影响分析与补偿021100274.3基于 STK 卫星数据的偏航角运动补偿为了验证卫星偏航角运动补偿方案的可行性
21、,针对卫星偏航角的变化进行如下仿真实验,实验条件均与上述仿真相同。由于卫星偏航角不断改变,小飞机位置对应每次采样过程其相应的数学表达为r=+(i)(15)式中,(i)为每次采样过程中的偏航角,随采样数 i改变,变化特性如图 8所示。进行星载偏航角运动补偿单像素成像仿真,对散斑光场区域投影 16 384张空间及时间上变化的散斑光场,对应仿真过程是 Hadamard矩阵与运动物体所处区域的内积运算,重构结果如图 9所示,分别计算峰值信图 7基于 STK卫星数据的滚动角与俯仰角运动补偿成像结果Fig.7Reconstruction results of object with roll angle
22、and pitch angle motion based on STK satellite data图 8偏航角随采样数的变化特性Fig.8Characteristics of yaw angle with sampling number光子学报02110028噪比来比较运动补偿前后图像还原质量,其中,图 9(c)为未补偿重构结果,峰值信噪比为 31.35 dB,图 9(d)为补偿后重构结果,峰值信噪比为 38.60 dB。仿真结果表明,光场追踪补偿策略适用于基于 STK 卫星数据的偏航角运动补偿,可以提升图像重构质量。4.4典型场景的星载单像素运动成像补偿在上述仿真结果的基础上继续对卫星探测
23、的实际图像进行卫星指向抖动的运动补偿仿真分析。运动目标分别为卫星观测的大坝,海岸线,公路和树木,运动场景即散斑光场覆盖范围为 128128,综合考虑由于卫星指向抖动引起的滚动角,俯仰角和偏航角的变化,进行 16 384次采样,对应每次采样过程其运动相应的数学表达为 xp=x+(i)yp=y+(i)r=+(i)(16)对应仿真过程是 Hadamard矩阵与运动物体所处区域的内积运算,重构结果如图 10所示,可以看到运动补偿方案有效地避免了由于桶探测器所探测到的总光强与照明散斑之间的关联性降低而造成的图像质量下降。进一步分别计算峰值信噪比来比较运动补偿前后图像还原质量,其中,大坝场景运动补偿前后的
24、峰值信噪比由 13.30 dB提升至 20.52 dB;海岸线场景运动补偿前后的峰值信噪比由 12.89 dB提升至 17.61 dB;公路场景运动补偿前后的峰值信噪比由 13.39 dB 提升至 21.03 dB;树木图场景运动补偿前后的峰值信噪比由 18.25 dB提升至 23.23 dB。仿真结果表明,光场追踪补偿策略同样适用于基于 STK 卫星数据综合考虑滚动角,俯仰角和偏航角的运动补偿,可以有效提升图像重构质量。图 9基于 STK卫星数据的偏航角运动补偿成像结果Fig.9Reconstruction results of object with yaw angle motion ba
25、sed on STK satellite data于跃,等:星载的单像素运动成像影响分析与补偿021100295 结论本文研究应用于卫星对地观测的单像素成像方法,设计的光场追踪补偿方法可以有效解决由于卫星滚动角,俯仰角和偏航角变化所造成的参考光场与探测器信号之间关联性丢失的问题,使重构图像峰值信噪比提升明显。运用 STK 软件仿真地球静止卫星数据,根据卫星的滚动角,俯仰角和偏航角对物体进行单像素成像追踪补偿,分别对卫星观测的大坝,海岸线,公路和树木进行运动补偿仿真,与未进行补偿的成像结果作对比,峰值信噪比提升了至少 4.7 dB,有效提高了成像质量,为星载单像素成像的运动补偿提供了一条有效技术
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 像素 运动 成像 影响 分析 补偿
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。