太阳直射光/行星反射光到达时间差分估测方法技术

本发明专利技术公开了一种太阳直射光/行星反射光到达时间差分估测方法，其特征在于：在航天器上搭载两个光电设备，分别用于获取太阳视频和观测行星反射光强轮廓；根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓；将构建的行星反射光强轮廓与观测的行星反射光强轮廓进行对比，得到到达时间差。本发明专利技术构建的行星反射光轮廓精度高，实时性好，对设备性能要求不高，能够支持高精度深空探测导航。

A Differential Estimation Method for the Time of Arrival of Direct Solar Light/Planetary Reflected Light

The invention discloses a method for estimating the arrival time difference of solar direct light/planetary reflected light, which is characterized by: two photoelectric devices are mounted on the spacecraft to obtain solar video and observe planetary reflected light intensity profiles respectively; planetary reflected light intensity profiles are constructed by using solar video according to the formation mechanism of planetary reflected light; planetary reflected light intensity profiles are constructed and observed. The time difference of arrival is obtained by comparing the reflected light intensity profiles of planets. The planetary reflected light contour constructed by the invention has high precision, good real-time performance, low requirement for equipment performance, and can support high precision deep space exploration and navigation.

【技术实现步骤摘要】
太阳直射光/行星反射光到达时间差分估测方法
本专利技术属于航天器天文自主导航领域，特别涉及一种太阳直射光与行星反射光到达时间之间差分值的估测方法。
技术介绍
导航信息对于深空探测的成败至关重要。受远距离和长时延的影响，地面站无法提供实时高精度的导航信息，特别是在捕获段。而航天器天文自主导航系统可以做到这一点。因此，对于捕获段而言，天文自主导航信息是极其重要的。面向捕获段的天文自主导航最终目的是提供相对于目标天体的位置和速度信息。近年来，能直接提供这类导航信息的天文导航方法已称为研究热点。目前，申请人已专利技术了《一种太阳耀斑到达时间差分量测及组合导航方法、系统》(2018年授权)，并发表SCI论文《SolarflareTDOAnavigationmethodusingdirectandreflectedlightforMarsexploration》IEEETransactionsonAerospaceandElectronicSystems。该方法将太阳直射光强波形与行星反射光强波形进行比较，得到二者的到达时间差，直接提供了航天器相对于目标天体的距离信息，并作为导航观测量。可在已有设备(超高速相机)、低精度的目标天体星历、无太阳耀斑计时模型的条件下正常工作。但经过进一步研究发现，该专利技术仅仅构建了太阳耀斑到达时间差分导航的基本框架，尚未涉及误差源分析。误差源分析是高精度量测值估计方法的前提条件。文献《GeometryerroranalysisinsolarDopplerdifferencenavigationforthecapturephase(面向捕获段的太阳多普勒差分导航中的几何误差分析)》，IEEETransactionsonAerospaceandElectronicSystems发现太阳这一巨大球形会对到太阳直射光与行星反射光达差分时间造成扩展效应。在转移段末端，扩展时间达到0.001s～0.01s量级，引起的测距误差约为300km～3000km。因此，将太阳直射光强波形与行星反射光强波形直接进行比较，会引起较大的误差。
技术实现思路
为了克服以上技术问题，本专利技术针对这一误差源的特性，提出了一种太阳直射光/行星反射光到达时间差分估测方法。本专利技术技术方案提供一种太阳直射光/行星反射光到达时间差分估测方法，在航天器上搭载两个光电设备，分别用于获取太阳视频和观测行星反射光强轮廓；根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓；将构建的行星反射光强轮廓与观测的行星反射光强轮廓进行对比，得到到达时间差。而且，所述根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓，实现如下，设有拍摄速率f、拍摄时长T、图像帧数S，在观测时间间隔[0,T]内，超高速相机获得的太阳直射光序列图像为Is,s＝1,2,…,S；设Is(m,n)为太阳直射光序列图像Is的第m行第n列像素，第s幅太阳直射光序列图像Is对应的行星反射光强子轮廓为ps，图像Is的第n列像素值之和等于ps的第n个元素值ps(n)；一幅太阳图像对应的时间扩展宽度为ΔT，设太阳图像直径为D个像素，一个像素宽度对应的时间展宽为ΔT/D；将观测期间内的行星反射光强子轮廓移位叠加，构建成行星反射光强轮廓。而且，将构建的行星反射光强轮廓与观测的行星反射光强轮廓进行对比，得到到达时间差，实现如下，设接收的行星反射光强轮廓矢量为PR，进行互相关计算如下，其中，J为互相关函数取最大值时对应的位移；i,j为变量，P(i)为构建的行星反射光强轮廓矢量P的第i个元素，PR(i+j)为接收的行星反射光强轮廓矢量PR的第i+j个元素；将该位移值转换为到达时间差δT，而且，光电设备采用超高速相机。而且，以火星为观测行星。而且，用于深空探测导航。本专利技术与现有技术相比的优点在于：(1)构建的行星反射光轮廓精度高。本专利技术考虑了太阳这一巨大球形的扩展效应，将太阳图像转化为一个光强子轮廓。因此，本专利技术构造的轮廓精度高。(2)实时性好。太阳图像转化为行星反射光强子轮廓的计算量小。且该计算过程可与太阳图像采集过程同步进行。因而，该专利技术的实时性好，能够用于高精度深空探测导航。(3)本专利技术对设备性能要求不高。市场上已有的超高速相机即可满足本专利技术的要求，无需研制新仪器，因此成本较低，具有重要的市场价值。附图说明图1为本专利技术实施例的流程示意图。图2为本专利技术实施例的直角坐标系示意图。图3为本专利技术实施例的太阳图像转行星反射光强子轮廓原理示意图。图4为本专利技术实施例的行星反射光强轮廓移位叠加的原理示意图。具体实施方式为了便于本领域普通技术人员理解和实施本专利技术，下面结合附图及实施例对本专利技术作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本专利技术，并不用于限定本专利技术。实施例中，以火星为观测行星。航天器搭载两个光电设备：第一超高速相机，用于获取太阳视频；第二超高速相机，用于观测行星反射光强轮廓。根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓，将构建的轮廓与观测的轮廓进行对比，即可得到到达时间差。将太阳视频转换为行星反射光强轮廓的过程如下：太阳视频可视为图像序列。设S为图像帧数，将第s幅太阳图像转换为行星反射光强第s个子轮廓，其中，s为编号，s＝1,2,…,S。。这样，太阳视频可转换为一系列的行星反射光强子轮廓。将这些子轮廓移位叠加，即可构建出行星反射光强轮廓。参见图1，本专利技术实施例提供了一种太阳直射光与行星反射光到达时间之间差分值的估测方法，包括以下步骤：步骤1：建立坐标系。本专利技术建立航天器位置的直角坐标系和超高速相机的敏感面阵坐标系，以实现图像配准的效果。建立航天器位置的直角坐标系为(x,y,z)直角坐标系，如图2所示。原点O为太阳质心。太阳质心、火星质心、航天器质心三点确定YZ平面，Z轴指向火星，X轴的正轴方向是根据右手定则确定的。火星质心和航天器质心坐标分别为(0,0,z1)、(0,y2,z2)。建立超高速相机的敏感面阵的坐标系。以太阳质心在敏感面阵上的投影为原点，以X’轴和Y’轴建立敏感面阵的坐标系。其中，X’轴与X轴平行，而Y’轴与X’轴垂直。拍摄图像的每一行与X’轴平行，每一列与Y’轴平行。步骤2：利用太阳视频(太阳序列图像)构建行星反射光强轮廓。利用超高速相机获取太阳序列图像。以Kirana-01M超高速摄像机为例，其像素数为924(W)×768(H)，最高拍摄速率为107帧/秒。在本专利技术中，拍摄速率f为105帧/秒。图像中心对准太阳质心。在观测时间间隔[0,T]内，超高速相机获得的太阳直射光序列图像为Is,s＝1,2,…,S，其中，S为图像帧数。拍摄速率f、拍摄时长T、图像帧数S之间的关系如下：S＝f·T(1)构建行星反射光强轮廓包括以下子步骤：步骤21：构建单幅图像的行星反射光强子轮廓。太阳图像转行星反射光强子轮廓的基本原理如图3所示。其过程如下：Is(m,n)为Is的第m行、第n列像素。第s幅图像对应的行星反射光强子轮廓为ps。图像Is的第n列像素值之和等于ps的第n个元素值ps(n)：一幅太阳图像对应的时间扩展宽度为ΔT。其中，火星位置矢量为rM，卫星位置矢量为r。R为太阳直径，c为光速。设太阳图像直径为D个像素，那么一个像素宽度对应的时间展宽为Δ本文档来自技高网...

【技术保护点】
1.一种太阳直射光/行星反射光到达时间差分估测方法，其特征在于：在航天器上搭载两个光电设备，分别用于获取太阳视频和观测行星反射光强轮廓；根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓；将构建的行星反射光强轮廓与观测的行星反射光强轮廓进行对比，得到到达时间差。

【技术特征摘要】
1.一种太阳直射光/行星反射光到达时间差分估测方法，其特征在于：在航天器上搭载两个光电设备，分别用于获取太阳视频和观测行星反射光强轮廓；根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓；将构建的行星反射光强轮廓与观测的行星反射光强轮廓进行对比，得到到达时间差。2.根据权利要求1所述太阳直射光/行星反射光到达时间差分估测方法，其特征在于：所述根据行星反射光形成的机理，利用太阳视频构建行星反射光强轮廓，实现如下，设有拍摄速率f、拍摄时长T、图像帧数S，在观测时间间隔[0,T]内，超高速相机获得的太阳直射光序列图像为Is,s＝1,2,…,S；设Is(m,n)为太阳直射光序列图像Is的第m行第n列像素，第s幅太阳直射光序列图像Is对应的行星反射光强子轮廓为ps，图像Is的第n列像素值之和等于ps的第n个元素值ps(n)；一幅太阳图像对应的时间扩展宽度为ΔT，设太阳图像直径为D个像素，一个像素宽度对应的时...

【专利技术属性】
技术研发人员：刘劲，李阳阳，潘超，刘兰兰，何宇飞，康志伟，
申请(专利权)人：武汉科技大学，
类型：发明
国别省市：湖北,42