首页 / 专利库 / 核能 / 仿星器 / 一种导航星座星间观测目标卫星图像模拟方法

一种导航星座星间观测目标卫星图像模拟方法

阅读:774发布:2020-05-14

专利汇可以提供一种导航星座星间观测目标卫星图像模拟方法专利检索,专利查询,专利分析的服务。并且一种导航 星座 星间观测目标卫星图像模拟方法,涉及星空图像模拟技术领域;主要包括如下步骤:步骤(一)、建立视线 坐标系 Olxlylzl、观测相机测量坐标系Ocxcyczc、观测相机像平面坐标系Opxpyp和观测卫星本体坐标系Obxbybzb;步骤(二)、计算目标卫星在星间观测相机像平面坐标系下实时 位置 坐标;步骤(三)、计算目标卫星的等效视星等和灰度值;步骤(四)、根据计算出的目标卫星的等效视星等m、观测相机像平面上各 像素 的灰度值gij以及目标卫星在星间观测相机像平面坐标系下位置坐标(x2d,y2d),进行仿真模拟;本 发明 提供一种星间观测目标卫星图像的模拟生成方法,该方法能够完成导航星座星间观测目标卫星图像的模拟生成。,下面是一种导航星座星间观测目标卫星图像模拟方法专利的具体信息内容。

1.一种导航星座星间观测目标卫星图像模拟方法,其特征在于:包括如下步骤:
步骤(一)、建立视线坐标系Olxlylzl、观测相机测量坐标系Ocxcyczc、观测相机像平面坐标系Opxpyp和观测卫星本体坐标系Obxbybzb;
步骤(二)、计算目标卫星在观测相机像平面坐标系下实时位置坐标;
计算观测卫星本体系相对J2000惯性系的初始时刻姿态四元数,即相机的初始惯性姿态Qbi0,具体计算方法如下:
计算J2000惯性系到视线坐标系Olxlylzl的转换矩阵Cli:
Cli=Ry(qβ)·Rz(-qε)                  (1)
式中,qε为目标卫星相对观测卫星的方位
qβ为目标卫星相对观测卫星高度角;
Ry(qβ)为视线坐标系绕着自身y轴旋转qβ角度;
Rz(qε)为视线坐标系绕着自身z轴旋转qε角度;
计算J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi:
Cbi=Cbl·Cli                      (2)
式中,Cbl为视线坐标系到观测卫星本体坐标系的转换矩阵,即卫星平台控制误差矩阵;
将Cbi转化为四元数,即得到初始时刻观测卫星本体坐标系相对于J2000惯性系观测相机的初始惯性姿态Qbi0;
将初始惯性姿态Qbi0作为初始值,代入姿态运动学方程,计算得到后续不同时刻的姿态四元数:
根据姿态运动学方程:
式中,Qbi为相机惯性姿态四元数;
为相机惯性姿态四元数的导数;
Ωbi为观测卫星本体坐标系相对于J2000惯性系的旋转角速度在观测卫星本体坐标系内的分量引入四元数表示;
表示四元数乘法;
计算目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d):
式中,(ex,ey,ez)为目标卫星相对矢量在观测相机像平面坐标系下的坐标;
f为观测相机的焦距;
步骤(三)、计算目标卫星的等效视星等和灰度值;
根据等效视星等的估算公式计算目标卫星的等效视星等m;
计算观测相机像平面上各像素的灰度值gij;
步骤(四)、根据计算出的目标卫星的等效视星等m、观测相机像平面上各像素的灰度值gij以及目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d),进行仿真模拟。
2.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(一)中,所述视线坐标系Olxlylzl的建立方法为:原点Ol位于观测卫星质心;xl轴与视线方向重合指向目标卫星;yl轴位于包含xl轴的纵向平面内并垂直于xl轴;zl轴由右手定则确定。
3.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(一)中,所述观测相机测量坐标系Ocxcyczc的建立方法为:原点Oc位于观测相机中心,xc轴与yc轴分别平行于相机探测器阵列两条互相垂直的边,zc轴沿着相机的光轴方向。
4.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(一)中,所述观测相机像平面坐标系Opxpyp的建立方法为:原点Op位于观测相机中心,xp轴与yp轴分别与相机测量坐标系的xc轴与yc轴重合。
5.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(一)中,所述观测卫星本体坐标系的建立方法为:原点Ob位于观测卫星质心;xb轴、yb轴和zb轴分别沿着观测卫星本体的惯量主轴方向且遵循右手定则。
6.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(二)中,将Cbi转化为四元数的具体方法为:
设根据公式(2)得到J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi为:
式中,c11、c12和c13为Cbi的第一行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的xb轴形成的方向余弦值;
c21、c22和c23为Cbi的第二行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的yb轴形成的方向余弦值;
c31、c32和c33为Cbi的第三行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体系的zb轴形成的方向余弦值;
将转换矩阵Cbi转化为四元数,得到观测相机的初始惯性姿态Qbi0=[q0 q1 q2 q3]T式中,q0为四元数Qbi0的标量部分,代表J2000惯性系转到观测卫星本体系所转过的欧拉角;
q1、q2和q3为四元数Qbi0的矢量部分,代表J2000惯性系转到观测卫星本体系所绕的欧拉轴的方向。
7.根据权利要求1所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(三)中,目标卫星的等效视星等m的计算方法为:
式中,m为目标卫星等效视星等;
S为目标卫星的截面积;
ρ为目标卫星的表面反射率;
f(σ)为相机、目标卫星和太阳的相位角函数;
R为观测卫星和目标卫星之间的距离。
8.根据权利要求7所述的一种导航星座星间观测目标卫星图像模拟方法,其特征在于:
所述步骤(三)中,观测相机像平面上各像素的灰度值gij的计算方法为:
式中,G为目标卫星在观测相机成像的灰度值;
gij为观测相机像平面上横坐标为i,纵坐标为j的像素的灰度值;
(xp0,yp0)为目标卫星像平面坐标(x2d,y2d)取整后得到的散焦像点的中心像素坐标;
Δx=x2d-xp0,为像素横坐标偏差;
Δy=y2d-yp0,为像素纵坐标偏差。

说明书全文

一种导航星座星间观测目标卫星图像模拟方法

技术领域

[0001] 本发明涉及一种星空图像模拟技术领域,特别是一种导航星座星间观测目标卫星图像模拟方法。

背景技术

[0002] 为保障在全球卫星导航系统地面站布设受限、突发故障、甚至战时被摧毁的情况下,卫星导航系统仍然能稳定连续地提供导航信息,导航卫星星座必须具备长时间、高精度的自主定轨能。目前,基于星间测距结合星间定向的星座自主定轨方法正得到广泛和深入地研究,而且有望得到全球卫星导航系统的工程应用。该方法通过导航卫星上安装的观测相机对同轨道相邻卫星及其背景恒星成像获取恒星-目标卫星距信息,联合星间测距信息一同输入至导航星座滤波器完成整星座的自主定轨。观测相机作为星座自主定轨系统的关键单机,在其研制、测试和验证过程中所必须的包含有目标卫星的星空图像至关重要。由于在轨实拍星空图像成本较大、技术要求较高、机会较少,而且难以获得,因此尽可能精确、尽可能逼近真实状况的目标卫星图像模拟需求迫切。

发明内容

[0003] 本发明的目的在于克服现有技术的上述不足,提供一种导航星座星间观测目标卫星图像模拟方法,该方法能够完成导航星座星间观测目标卫星图像的模拟生成。
[0004] 本发明的上述目的是通过如下技术方案予以实现的:
[0005] 一种导航星座星间观测目标卫星图像模拟方法,包括如下步骤:
[0006] 步骤(一)、建立视线坐标系Olxlylzl、观测相机测量坐标系Ocxcyczc、观测相机像平面坐标系Opxpyp和观测卫星本体坐标系Obxbybzb;
[0007] 步骤(二)、计算目标卫星在观测相机像平面坐标系下实时位置坐标;
[0008] 计算观测卫星本体系相对J2000惯性系的初始时刻姿态四元数,即相机的初始惯性姿态Qbi0,具体计算方法如下:
[0009] 计算J2000惯性系到视线坐标系Olxlylzl的转换矩阵Cli:
[0010] Cli=Ry(qβ)·Rz(-qε)                   (1)
[0011] 式中,qε为目标卫星相对观测卫星的方位角;
[0012] qβ为目标卫星相对观测卫星高度角;
[0013] Ry(qβ)为视线坐标系绕着自身y轴旋转qβ角度;
[0014] Rz(qε)为视线坐标系绕着自身z轴旋转qε角度;
[0015] 计算J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi:
[0016] Cbi=Cbl·Cli                      (2)
[0017] 式中,Cbl为视线坐标系到观测卫星本体坐标系的转换矩阵,即卫星平台控制误差矩阵;
[0018] 将Cbi转化为四元数,即得到初始时刻观测卫星本体坐标系相对于J2000惯性系观测相机的初始惯性姿态Qbi0;
[0019] 将初始惯性姿态Qbi0作为初始值,代入姿态运动学方程,计算得到后续不同时刻的姿态四元数:
[0020] 根据姿态运动学方程:
[0021]
[0022] 式中,Qbi为相机惯性姿态四元数;
[0023] 为相机惯性姿态四元数的导数;
[0024] Ωbi为观测卫星本体坐标系相对于J2000惯性系的旋转角速度在观测卫星本体坐标系内的分量引入四元数表示;
[0025] 表示四元数乘法;
[0026] 计算目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d):
[0027]
[0028] 式中,(ex,ey,ez)为目标卫星相对矢量在观测相机像平面坐标系下的坐标;
[0029] f为观测相机的焦距;
[0030] 步骤(三)、计算目标卫星的等效视星等和灰度值;
[0031] 根据等效视星等的估算公式计算目标卫星的等效视星等m;
[0032] 计算观测相机像平面上各像素的灰度值gij;
[0033] 步骤(四)、根据计算出的目标卫星的等效视星等m、观测相机像平面上各像素的灰度值gij以及目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d),进行仿真模拟。
[0034] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(一)中,所述视线坐标系Olxlylzl的建立方法为:原点Ol位于观测卫星质心;xl轴与视线方向重合指向目标卫星;yl轴位于包含xl轴的纵向平面内并垂直于xl轴;zl轴由右手定则确定。
[0035] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(一)中,所述观测相机测量坐标系Ocxcyczc的建立方法为:原点Oc位于观测相机中心,xc轴与yc轴分别平行于相机探测器阵列两条互相垂直的边,zc轴沿着相机的光轴方向。
[0036] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(一)中,所述观测相机像平面坐标系Opxpyp的建立方法为:原点Op位于观测相机中心,xp轴与yp轴分别与相机测量坐标系的xc轴与yc轴重合。
[0037] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(一)中,所述观测卫星本体坐标系的建立方法为:原点Ob位于观测卫星质心;xb轴、yb轴和zb轴分别沿着观测卫星本体的惯量主轴方向且遵循右手定则。
[0038] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(二)中,将Cbi转化为四元数的具体方法为:
[0039] 设根据公式(2)得到J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi为:
[0040]
[0041] 式中,c11、c12和c13为Cbi的第一行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的xb轴形成的方向余弦值;
[0042] c21、c22和c23为Cbi的第二行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的yb轴形成的方向余弦值;
[0043] c31、c32和c33为Cbi的第三行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体系的zb轴形成的方向余弦值;
[0044] 将转换矩阵Cbi转化为四元数,得到观测相机的初始惯性姿态Qbi0=[q0 q1 q2 q3]T[0045]
[0046] 式中,q0为四元数Qbi0的标量部分,代表J2000惯性系转到观测卫星本体系所转过的欧拉角;
[0047] q1、q2和q3为四元数Qbi0的矢量部分,代表J2000惯性系转到观测卫星本体系所绕的欧拉轴的方向。
[0048] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(三)中,目标卫星的等效视星等m的计算方法为:
[0049]
[0050] 式中,m为目标卫星等效视星等;
[0051] S为目标卫星的截面积;
[0052] ρ为目标卫星的表面反射率;
[0053] f(σ)为相机、目标卫星和太阳的相位角函数;
[0054] R为观测卫星和目标卫星之间的距离。
[0055] 在上述的一种导航星座星间观测目标卫星图像模拟方法,所述步骤(三)中,观测相机像平面上各像素的灰度值gij的计算方法为:
[0056]
[0057] 式中,G为目标卫星在观测相机成像的灰度值;
[0058] gij为观测相机像平面上横坐标为i,纵坐标为j的像素的灰度值;
[0059] (xp0,yp0)为目标卫星像平面坐标(x2d,y2d)取整后得到的散焦像点的中心像素坐标;
[0060] Δx=x2d-xp0,为像素横坐标偏差;
[0061] Δy=y2d-yp0,为像素纵坐标偏差。
[0062] 本发明与现有技术相比具有如下优点:
[0063] (1)本发明紧密结合导航星座观测相机的成像模式,全过程、高逼真地对观测相机所能观测的目标卫星图像进行模拟,针对性较强;
[0064] (2)本发明仅根据观测卫星J2000惯性系下位置、目标卫星J2000惯性系下位置和卫星平台的姿态特性信息,通过一系列相对简单的数学运算,即可实现对目标卫星图像的模拟,工程可行性较高。附图说明
[0065] 图1为本发明导航星座星间观测目标卫星图像模拟方法的流程图

具体实施方式

[0066] 下面结合附图和具体实施例对本发明作进一步详细的描述:
[0067] 如图1所示为导航星座星间观测目标卫星图像模拟方法的流程图,由图可知,一种导航星座星间观测目标卫星图像模拟方法,包括如下步骤:
[0068] 步骤(一)、导航星座观测相机成像模式为:初始时刻观测相机转台控制相机光轴指向目标卫星,然后转台死;相机和卫星平台固联,姿态跟随卫星平台变化;在设定的周期内,转台控制相机光轴重新指向目标卫星,重复上述过程。为模拟该成像过程,建立如下坐标系:视线坐标系Olxlylzl、观测相机测量坐标系Ocxcyczc、观测相机像平面坐标系Opxpyp和观测卫星本体坐标系Obxbybzb;
[0069] 其中,所述视线坐标系Olxlylzl的建立方法为:原点Ol位于观测卫星质心;xl轴与视线方向重合指向目标卫星;yl轴位于包含xl轴的纵向平面内并垂直于xl轴;zl轴由右手定则确定。
[0070] 所述观测相机测量坐标系Ocxcyczc的建立方法为:原点Oc位于观测相机中心,xc轴与yc轴分别平行于相机探测器阵列两条互相垂直的边,zc轴沿着相机的光轴方向。
[0071] 所述观测相机像平面坐标系Opxpyp的建立方法为:原点Op位于观测相机中心,xp轴与yp轴分别与相机测量坐标系的xc轴与yc轴重合。
[0072] 所述观测卫星本体坐标系的建立方法为:原点Ob位于观测卫星质心;xb轴、yb轴和zb轴分别沿着观测卫星本体的惯量主轴方向且遵循右手定则。
[0073] 步骤(二)、根据观测卫星和目标卫星在J2000惯性系下的实时位置信息和观测卫星的姿态运动特性,通过相机的初始惯性姿态确定、惯性姿态运动学解算、坐标系转换,完成目标卫星在观测相机像平面坐标系下实时位置坐标计算;
[0074] 计算观测卫星本体系相对J2000惯性系的初始时刻姿态四元数,即相机的初始惯性姿态Qbi0,具体计算方法如下:
[0075] 计算J2000惯性系到视线坐标系Olxlylzl的转换矩阵Cli:
[0076] Cli=Ry(qβ)·Rz(-qε)                  (1)
[0077] 式中,qε为目标卫星相对观测卫星的方位角;
[0078] qβ为目标卫星相对观测卫星高度角;
[0079] Ry(qβ)为视线坐标系绕着自身y轴旋转qβ角度;
[0080] Rz(qε)为视线坐标系绕着自身z轴旋转qε角度;
[0081] 计算J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi:
[0082] Cbi=Cbl·Cli                      (2)
[0083] 式中,Cbl为视线坐标系到观测卫星本体坐标系的转换矩阵,即卫星平台控制误差矩阵;
[0084] 将Cbi转化为四元数,即得到初始时刻观测卫星本体坐标系相对于J2000惯性系观测相机的初始惯性姿态Qbi0;
[0085] 将Cbi转化为四元数的具体方法为:
[0086] 设根据公式(2)得到J2000惯性系到观测卫星本体坐标系的转换矩阵Cbi为:
[0087]
[0088] 式中,c11、c12和c13为Cbi的第一行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的xb轴形成的方向余弦值;
[0089] c21、c22和c23为Cbi的第二行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体坐标系的yb轴形成的方向余弦值;
[0090] c31、c32和c33为Cbi的第三行元素,分别表示J2000惯性系的三个坐标轴与观测卫星本体系的zb轴形成的方向余弦值;
[0091] 将转换矩阵Cbi转化为四元数,得到观测相机的初始惯性姿态Qbi0=[q0 q1 q2 q3]T[0092]
[0093] 式中,q0为四元数Qbi0的标量部分,代表J2000惯性系转到观测卫星本体系所转过的欧拉角;
[0094] q1、q2和q3为四元数Qbi0的矢量部分,代表J2000惯性系转到观测卫星本体系所绕的欧拉轴的方向。
[0095] 将初始惯性姿态Qbi0作为初始值,代入姿态运动学方程,计算得到后续不同时刻的姿态四元数:
[0096] 根据姿态运动学方程:
[0097]
[0098] 式中,Qbi为相机惯性姿态四元数;
[0099] 为相机惯性姿态四元数的导数;
[0100] Ωbi为观测卫星本体坐标系相对于J2000惯性系的旋转角速度在观测卫星本体坐标系内的分量引入四元数表示;
[0101] 表示四元数乘法;
[0102] 最后计算目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d):
[0103]
[0104] 式中,(ex,ey,ez)为目标卫星相对矢量在观测相机像平面坐标系下的坐标;
[0105] f为观测相机的焦距;
[0106] 步骤(三)、计算目标卫星的等效视星等和灰度值,完成像点高斯灰度扩散,最终完成目标卫星图像模拟;
[0107] 根据等效视星等的估算公式计算目标卫星的等效视星等:
[0108]
[0109] 式中,m为目标卫星等效视星等;
[0110] S为目标卫星的截面积;
[0111] ρ为目标卫星的表面反射率;
[0112] f(σ)为相机、目标卫星和太阳的相位角函数;
[0113] R为观测卫星和目标卫星之间的距离;
[0114] 计算观测相机像平面上各像素的灰度值gij
[0115]
[0116] 式中,G为目标卫星在观测相机成像的灰度值;
[0117] gij为观测相机像平面上横坐标为i,纵坐标为j的像素的灰度值;
[0118] (xp0,yp0)为目标卫星像平面坐标(x2d,y2d)取整后得到的散焦像点的中心像素坐标;
[0119] Δx=x2d-xp0,为像素横坐标偏差;
[0120] Δy=y2d-yp0,为像素纵坐标偏差;
[0121] 步骤(四)、根据计算出的目标卫星的等效视星等m、观测相机像平面上各像素的灰度值gij以及目标卫星在观测相机像平面坐标系下位置坐标(x2d,y2d),进行仿真模拟。
[0122] 本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈