CN115200573A - 空间目标的测量装备定位方法、***和存储介质 - Google Patents

空间目标的测量装备定位方法、***和存储介质 Download PDF

Info

Publication number
CN115200573A
CN115200573A CN202211092032.XA CN202211092032A CN115200573A CN 115200573 A CN115200573 A CN 115200573A CN 202211092032 A CN202211092032 A CN 202211092032A CN 115200573 A CN115200573 A CN 115200573A
Authority
CN
China
Prior art keywords
coordinate system
axis
dimensional position
data
space target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211092032.XA
Other languages
English (en)
Other versions
CN115200573B (zh
Inventor
宗永红
高昕
李希宇
孙亮亮
胡蕾
雷呈强
郑东昊
师恒
孙锐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
63921 Troops of PLA
Original Assignee
63921 Troops of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 63921 Troops of PLA filed Critical 63921 Troops of PLA
Priority to CN202211092032.XA priority Critical patent/CN115200573B/zh
Publication of CN115200573A publication Critical patent/CN115200573A/zh
Application granted granted Critical
Publication of CN115200573B publication Critical patent/CN115200573B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • G01C21/025Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/66Analysis of geometric attributes of image moments or centre of gravity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Astronomy & Astrophysics (AREA)
  • Geometry (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及一种空间目标的测量装备定位方法、***和存储介质,获取空间目标相对于测量装备的位置测量数据;将空间目标相对于测量装备的位置测量数据从以测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到第一三维位置数据;将测量装备的第一三维位置数据从第二天球坐标系转换至地球坐标系,得到第二三维位置数据;根据空间目标的精密星历,得到预设观测时刻在地球坐标系下的空间目标的第三三维位置数据;基于第二三维位置数据和第三三维位置数据,确定地球坐标系下的测量装备的位置定位数据;由此,本发明能够不依赖于GNSS和大地测量的新方法,以完成天基测量装备的高精度定位。

Description

空间目标的测量装备定位方法、***和存储介质
技术领域
本发明涉及空间目标测量技术领域,具体涉及一种空间目标的测量装备定位方法、***和存储介质。
背景技术
空间目标的测量装备是太空态势感知体系的重要设备,为完成对空间目标的高精度测量与定轨,需要知道测量装备自身的精确位置。根据测量装备自身的精确位置和获取的测量数据,才能完成对目标的定位。
目前空间目标测量装备主要分为地基和天基设备,现在的地基设备大多通过高精度大地测量或者通过全球导航卫星***(GNSS)测量获得自身精确位置坐标,或者通过拍摄恒星方式解算站址自身的天文经纬度,结合高程地图获得精确站址。而对于天基(主要为机载或星载)测量装备而言,在无GNSS信号或GNSS信号受到干扰时,测量装备将无法完成自身定位,进而导致空间目标测量精度下降或者失效,为此,需要研究不依赖于GNSS和大地测量的新方法,以完成天基测量装备的高精度定位。
发明内容
本发明主要解决的技术问题是如何对天基测量装备进行高精度定位。
根据第一方面,一种实施例中提供一种空间目标的测量装备定位方法,包括:
获取空间目标相对于测量装备的位置测量数据,所述位置测量数据为以所述测量装备为中心的第一天球坐标系下的位置测量数据;
将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据;
将所述测量装备的第一三维位置数据从所述第二天球坐标系转换至所述地球坐标系,得到在地球坐标系下的所述测量装备相对于所述空间目标的第二三维位置数据;
获取所述空间目标的精密星历,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据;
基于所述第二三维位置数据和所述第三三维位置数据,确定所述地球坐标系下的所述测量装备的位置定位数据。
根据第二方面,一种实施例中提供一种空间目标的测量装备定位***,包括:
空间目标,所述空间目标为已知精密星历的空间目标;
测量装备,用于观测并跟踪所述空间目标,以获取所述空间目标的图像数据;
搭载平台,用于搭载所述测量装备;
数据处理服务器,用于获取所述空间目标的图像数据,并根据所述空间目标的图像数据,确定所述空间目标相对于测量装备的位置测量数据;
所述数据处理服务器还用于将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据;将所述测量装备的第一三维位置数据从所述第二天球坐标系转换至所述地球坐标系,得到在地球坐标系下的所述测量装备相对于所述空间目标的第二三维位置数据;获取所述空间目标的精密星历,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据;基于所述第二三维位置数据和所述第三三维位置数据,确定所述地球坐标系下的所述测量装备的位置定位数据。
根据第三方面,一种实施例中提供一种计算机可读存储介质,所述介质上存储有程序,所述程序能够被处理器执行以实现如上述实施例所述的方法。
依据上述实施例的空间目标的测量装备定位方法/***,获取空间目标相对于测量装备的位置测量数据;将空间目标相对于测量装备的位置测量数据从以测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到第一三维位置数据;将测量装备的第一三维位置数据从第二天球坐标系转换至地球坐标系,得到第二三维位置数据;根据空间目标的精密星历,得到预设观测时刻在地球坐标系下的空间目标的第三三维位置数据;基于第二三维位置数据和第三三维位置数据,确定地球坐标系下的测量装备的位置定位数据;由此,本发明能够不依赖于GNSS和大地测量方法,完成天基测量装备的高精度定位。
附图说明
图1为一种实施例的空间目标的测量装备定位***的结构示意图;
图2为一种实施例的空间目标的测量装备定位方法流程图;
图3为第一天球坐标系和过渡坐标系之间的转换关系示意图;
图4为过渡坐标系和第二天球坐标系之间的转换关系示意图;
图5为第二天球坐标系和地球坐标系之间的转换关系示意图。
具体实施方式
下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中,很多细节描述是为了使得本申请能被更好的理解。然而,本领域技术人员可以毫不费力的认识到,其中部分特征在不同情况下是可以省略的,或者可以由其他元件、材料、方法所替代。在某些情况下,本申请相关的一些操作并没有在说明书中显示或者描述,这是为了避免本申请的核心部分被过多的描述所淹没,而对于本领域技术人员而言,详细描述这些相关操作并不是必要的,他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。
另外,说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时,方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此,说明书和附图中的各种顺序只是为了清楚描述某一个实施例,并不意味着是必须的顺序,除非另有说明其中某个顺序是必须遵循的。
本文中为部件所编序号本身,例如“第一”、“第二”等,仅用于区分所描述的对象,不具有任何顺序或技术含义。而本申请所说“连接”、“联接”,如无特别说明,均包括直接和间接连接(联接)。
天基测量装备是指搭载在空中平台(例如飞机、卫星)的测量装备,针对天基的空间目标测量装备在GNSS信号受到干扰或者测量装备高于GPS卫星轨道时,无法完成自身定位的问题,本发明实施例通过测量装备对已知精密星历的空间目标的拍摄测量实现,无需额外设备保障,操作简单,运算量小,实现天基测量装备的高精度定位,从而实现对其他空间目标的定轨测量。
请参考图1,图1为一种实施例的空间目标的测量装备定位***的结构示意图,以下简称测量装备定位***,该测量装备定位***包括:空间目标101、测量装备102、搭载平台103和数据处理服务器104。
空间目标101为已知精密星历的空间目标,例如人造卫星,以北斗卫星为例,其精密星历的定位误差小于0.2m。
测量设备102稳定固定在搭载平台103上,测量设备102用于观测并跟踪空间目标101,常见的搭载平台103包括:飞机、卫星等,在天气允许的情况下,测量装备102先进行观测前的星校,测量装备102接收空间目标101的观测计划,并加载空间目标101的预报信息并对其进行解算,计算出空间目标101的过境具体时间,以及测量装备102的方位角度、俯仰角度等信息,最后生成测量装备102自身的观测任务,并在指定的天区进行等待。一旦空间目标101进入测量装备102的视场,测量装备102将自动跟踪观测空间目标101,获得多幅空间目标101和恒星的同帧画幅图像数据,即获取到空间目标101的图像数据。
数据处理服务器104与测量装备102通信连接,数据处理服务器104获取测量装备102输出的空间目标101的图像数据,通过天文定位方式获取空间目标101相对于测量装备102的视赤经和视赤纬数据,即获取到测量装备102的位置测量数据。
数据处理服务器104获取测量装备102的位置测量数据具体方法如下:
(1)根据获取的空间目标101的图像数据、观测时间等信息,快速检索星表,找出测量装备102的视场中的恒星(定标星),根据找出的恒星位置,计算恒星对应的理想坐标数据,并通过修正矩分析法计算恒星的实测坐标数据。
(2)根据恒星的个数,采用不同的探测器处理模型,计算出探测器处理模型的系数。
(3)对测量装备102输出的空间目标101的图像数据进行背景非均匀性消除、星图配准、差帧去除星图、图像多帧叠加、二值化处理等处理,采用基于各向异性高斯曲面拟合的星点质心提取方法提取空间目标101的星象,并通过修正矩分析法计算空间目标101的实测坐标数据。
(4)根据探测器处理模型的系数和空间目标101的实测坐标数据,计算空间目标101的理想坐标数据,进而得到空间目标101在以测量装备102为中心的第一天球坐标系中的位置测量数据,即空间目标101相对于测量装备102的视赤经和视赤纬。
本实施例通过以上方法实现对不少于2个空间目标多个弧段的观测,获取多组空间目标相对于测量设备102的位置测量数据。
基于以上方法获取的空间目标101相对于测量设备102的位置测量数据,数据处理服务器104对测量设备102进行定位计算。请参考图2,图2为一种实施例的空间目标的测量装备定位方法流程图,该测量装备定位方法应用于数据处理服务器104,其包括步骤201至步骤205。
步骤201:获取空间目标101相对于测量装备102的位置测量数据,该位置测量数据为以测量装备102为中心的第一天球坐标系下的位置测量数据。
步骤202:将空间目标101相对于测量装备102的位置测量数据从以测量装备102为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到测量装备102在第二天球坐标系下的第一三维位置数据。本实施例基于构建的过渡坐标系,将位置测量数据从第一天球坐标系转换到第二天球坐标系。
步骤203:将测量装备102的第一三维位置数据从第二天球坐标系转换至地球坐标系,得到在地球坐标系下的测量装备相对于空间目标的第二三维位置数据。
步骤204:获取空间目标101的精密星历,根据空间目标101的精密星历,得到预设观测时刻在地球坐标系下的空间目标101的第三三维位置数据。
步骤205:基于第二三维位置数据和第三三维位置数据,确定地球坐标系下的测量装备102的位置定位数据。
下面对本实施例涉及的坐标系进行说明。
第一天球坐标系(O1-X1Y1Z1)以测量装备102的质心为原点O1,以原点O1指向春分点方向为X1轴的正方向,Z1轴平行于地球自转轴,由地心指向北天极方向为Z1轴的正方向,Y1轴位于过原点的平面内,且Y1轴与X1轴 、Z1轴构成右手直角坐标系。
过渡坐标系(O2-X2Y2Z2)为以空间目标101的质心为原点,以原点O2指向春分点方向为X2轴的正方向,Z2轴平行于地球自转轴,由地心指向北天极方向为Z2轴的正方向,Y2轴位于过原点的平面内,且Y2轴与X2轴、Z2轴构成右手直角坐标系。
第二天球坐标系(O3-X3Y3Z3)以地心为原点O3,Z3轴平行于地球自转轴,由地心指向北天极方向为Z3轴的正方向,X3轴在赤道面上,由地心指向春分点方向为X3轴的正方向,Y3轴在赤道面上,且Y3轴与X3轴、Z3轴构成右手直角坐标系。
地球坐标系(O4-X4Y4Z4)以地心为原点O4,由原点O4指向给定的协议地球极CIP的方向为Z4轴的正方向,X4轴在赤道平面内,由原点O4指向给定的经度零点为正方向,Y4轴位于赤道平面内,Y4轴与X4轴、Z4轴构成右手直角坐标系。
在一实施例中,步骤202中,将空间目标101相对于测量装备102的位置测量数据从以测量装备102为中心的第一天球坐标系(O1-X1Y1Z1)转换至以地心为中心的第二天球坐标系(O3-X3Y3Z3),得到测量装备102在第二天球坐标系(O3-X3Y3Z3)下的第一三维位置数据,包括:
步骤2021:构建过渡坐标系(O2-X2Y2Z2)。为了实现测量装备的定位,本发明提出建立过渡坐标系(O2-X2Y2Z2),过渡坐标系(O2-X2Y2Z2)与以测量装备102为中心的第一天球坐标系相比,除坐标原点不同之外,各坐标轴的指向方向相同。
步骤2022:将空间目标101相对于测量装备102的位置测量数据从第一天球坐标系(O1-X1Y1Z1)转换到过渡坐标系(O2-X2Y2Z2),得到测量装备102在过渡坐标系(O2-X2Y2Z2)下的第四三维位置数据。
如图3所示,图3示出了第一天球坐标系(O1-X1Y1Z1)和过渡坐标系(O2-X2Y2Z2)之间的转换关系,具体为:
空间目标
Figure DEST_PATH_IMAGE001
在第
Figure DEST_PATH_IMAGE002
时刻相对于测量装备102的位置测量数据为(
Figure DEST_PATH_IMAGE003
Figure DEST_PATH_IMAGE004
),位置测量数据(
Figure 442135DEST_PATH_IMAGE003
Figure 202018DEST_PATH_IMAGE004
)是在第一天球坐标系(O1-X1Y1Z1)下。假定此时空间目标
Figure DEST_PATH_IMAGE005
在第
Figure 247335DEST_PATH_IMAGE002
时刻到测量装备102的距离为
Figure DEST_PATH_IMAGE006
,则在过渡坐标系(O2-X2Y2Z2)中,测量装备102在第
Figure 455593DEST_PATH_IMAGE002
时刻的位置数据(
Figure DEST_PATH_IMAGE007
Figure DEST_PATH_IMAGE008
)及测量装备102与空间目标之间在第
Figure 177955DEST_PATH_IMAGE002
时刻的距离
Figure DEST_PATH_IMAGE009
分别为:
Figure DEST_PATH_IMAGE011
其中,
Figure 531576DEST_PATH_IMAGE003
为空间目标
Figure 142817DEST_PATH_IMAGE001
在第
Figure 810559DEST_PATH_IMAGE002
时刻的视赤经,
Figure 338492DEST_PATH_IMAGE004
为空间目标
Figure 800697DEST_PATH_IMAGE001
在第
Figure 397770DEST_PATH_IMAGE002
时刻的视赤纬,
Figure 931519DEST_PATH_IMAGE007
为测量装备102在第
Figure 923746DEST_PATH_IMAGE002
时刻的视赤经,
Figure 963377DEST_PATH_IMAGE008
为测量装备102在第
Figure 408265DEST_PATH_IMAGE002
时刻的视赤纬,
Figure 542443DEST_PATH_IMAGE006
为在第一天球坐标系(O1-X1Y1Z1)下空间目标
Figure 654756DEST_PATH_IMAGE001
在第
Figure 823876DEST_PATH_IMAGE002
时刻到测量装备102的距离,
Figure 21639DEST_PATH_IMAGE009
为过渡坐标系(O2-X2Y2Z2)下测量装备102与空间目标
Figure 693929DEST_PATH_IMAGE001
之间在第
Figure 660748DEST_PATH_IMAGE002
时刻的距离。
将测量装备102在第
Figure 776603DEST_PATH_IMAGE002
时刻的位置数据(
Figure 196083DEST_PATH_IMAGE007
Figure 875326DEST_PATH_IMAGE008
)及测量装备102与空间目标
Figure 493389DEST_PATH_IMAGE001
之间在第
Figure 481942DEST_PATH_IMAGE002
时刻的距离
Figure 513352DEST_PATH_IMAGE009
转换为三维空间直角坐标,可得到第四三维位置数据
Figure DEST_PATH_IMAGE012
,具体公式为:
Figure DEST_PATH_IMAGE013
其中,
Figure DEST_PATH_IMAGE014
为测量装备102在过渡坐标系(O2-X2Y2Z2)下第
Figure DEST_PATH_IMAGE015
时刻的第四三维位置数据。
步骤2023:将测量装备102在过渡坐标系(O2-X2Y2Z2)下的第四三维位置数据转换到第二天球坐标系(O3-X3Y3Z3)中,得到测量装备102在第二天球坐标系(O3-X3Y3Z3)下的第一三维位置数据。
除坐标原点不同外,过渡坐标系(O2-X2Y2Z2)和第二天球坐标系(O3-X3Y3Z3)的坐标轴指向天***置的表示方法相同,过渡坐标系(O2-X2Y2Z2)可看作第二天球坐标系(O3-X3Y3Z3)的原点平移到以空间目标101为原点的坐标系。如图4所示,图4示出了过渡坐标系(O2-X2Y2Z2)和第二天球坐标系(O3-X3Y3Z3)的转换关系,则测量装备102在第二天球坐标系(O3-X3Y3Z3)的第一三维位置数据为:
Figure DEST_PATH_IMAGE016
其中,
Figure DEST_PATH_IMAGE017
为测量装备102在第二天球坐标系(O3-X3Y3Z3)的第一三维位置数据;
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE019
Figure DEST_PATH_IMAGE020
为空间目标
Figure DEST_PATH_IMAGE021
在第二天球坐标系(O3-X3Y3Z3)的三维位置数据。
在一实施例中,步骤203中,通过以下方式将测量装备102的第一三维位置数据从第二天球坐标系(O3-X3Y3Z3)转换至地球坐标系系(O4-X4Y4Z4),得到在地球坐标系(O4-X4Y4Z4)下的测量装备102相对于空间目标101的第二三维位置数据。
第二天球坐标系(O3-X3Y3Z3)与地球坐标系系(O4-X4Y4Z4)之间的转换关系具体为:
定义旋转矩阵为:
Figure DEST_PATH_IMAGE023
其中,
Figure DEST_PATH_IMAGE024
为绕X轴旋转
Figure DEST_PATH_IMAGE025
的旋转矩阵;
Figure DEST_PATH_IMAGE026
为绕Y轴旋转
Figure DEST_PATH_IMAGE027
的旋转矩阵;
Figure DEST_PATH_IMAGE028
为绕Z轴旋转
Figure DEST_PATH_IMAGE029
的旋转矩阵。
假设通过精密星历获取的空间目标
Figure 248483DEST_PATH_IMAGE021
Figure DEST_PATH_IMAGE030
时刻的三维坐标为
Figure DEST_PATH_IMAGE031
。空间目标
Figure 517791DEST_PATH_IMAGE021
Figure DEST_PATH_IMAGE032
时刻转换至第二天球坐标系(O3-X3Y3Z3)的转换矩阵为
Figure DEST_PATH_IMAGE033
,则空间目标
Figure 539230DEST_PATH_IMAGE021
Figure 933302DEST_PATH_IMAGE032
时刻在第二天球坐标系(O3-X3Y3Z3)的三维坐标
Figure DEST_PATH_IMAGE034
为:
Figure DEST_PATH_IMAGE035
其中,
Figure DEST_PATH_IMAGE036
为转换矩阵。
在一实施例中,步骤204中,根据空间目标101的精密星历,得到预设观测时刻在地球坐标系下的空间目标101的第三三维位置数据,包括:
在一插值区间内,对空间目标101的精密星历内插任意时刻的空间目标,以获取预设观测时刻在地球坐标系下的空间目标101的第三三维位置数据。
本实施例根据以下公式对空间目标
Figure 610140DEST_PATH_IMAGE021
在插值区间中,内插任意时刻的空间目标:
Figure DEST_PATH_IMAGE037
Figure DEST_PATH_IMAGE038
Figure DEST_PATH_IMAGE039
其中,
Figure DEST_PATH_IMAGE040
表示预设观测时刻t空间目标
Figure 281424DEST_PATH_IMAGE021
的插值位置数据,
Figure DEST_PATH_IMAGE041
为插值样本点历元时刻,
Figure DEST_PATH_IMAGE042
Figure 473764DEST_PATH_IMAGE041
时刻空间目标
Figure 355132DEST_PATH_IMAGE021
的位置数据。通过不断变换插值区间,待插值点t位于插值区间的中心位置,计算预设观测时刻t的空间目标
Figure 976607DEST_PATH_IMAGE021
的位置数据,即得到地球坐标系下的空间目标
Figure 627031DEST_PATH_IMAGE021
的第三三维位置数据。
在本实施例中,为了提高测量装备102的定位精度,空间目标101的数量不少于2个,即,测量装备102对不少于2颗人造卫星进行观测测量,设根据三角公式可求得测量装备102在地球坐标系(O4-X4Y4Z4)中的三维坐标
Figure DEST_PATH_IMAGE043
,通过上述的坐标转换,可得:
Figure DEST_PATH_IMAGE044
Figure DEST_PATH_IMAGE045
式中,
Figure DEST_PATH_IMAGE046
分别表示在
Figure DEST_PATH_IMAGE047
时刻第1个、第2个…第
Figure DEST_PATH_IMAGE048
个空间目标101与测量装备102之间的距离;
Figure DEST_PATH_IMAGE049
Figure DEST_PATH_IMAGE050
表示在
Figure 613441DEST_PATH_IMAGE047
时刻通过精密星历获取的第
Figure 106740DEST_PATH_IMAGE048
个空间目标101的第三三维位置数据;
Figure DEST_PATH_IMAGE051
表示在
Figure 407271DEST_PATH_IMAGE047
时刻测量装备102给出的第
Figure 769593DEST_PATH_IMAGE048
个空间目标101的视赤经;
Figure DEST_PATH_IMAGE052
表示在
Figure 989222DEST_PATH_IMAGE047
时刻测量装备102给出第
Figure 923811DEST_PATH_IMAGE048
个空间目标101的视赤纬。
由以上公式联立方程组可求得未知解
Figure 762454DEST_PATH_IMAGE046
,从而得到测量装备在地球坐标系中的位置为:
Figure DEST_PATH_IMAGE053
其中,
Figure DEST_PATH_IMAGE054
表示中间参数;
Figure DEST_PATH_IMAGE055
Figure DEST_PATH_IMAGE056
Figure DEST_PATH_IMAGE057
表示在
Figure 433475DEST_PATH_IMAGE047
时刻通过精密星历获取的第
Figure DEST_PATH_IMAGE058
个空间目标101在地球坐标系下在X4轴和Y4轴方向的位置数据,
Figure DEST_PATH_IMAGE059
表示在
Figure 558426DEST_PATH_IMAGE047
时刻测量装备102给出第
Figure 216197DEST_PATH_IMAGE058
个空间目标101的视赤经,
Figure DEST_PATH_IMAGE060
表示在
Figure 452006DEST_PATH_IMAGE047
时刻测量装备102给出第
Figure 931529DEST_PATH_IMAGE048
个空间目标101的视赤经。
在本发明实施例中,通过测量装备102自身对空间目标101的拍摄测量实现,无需额外设备保障,操作简单,运算量小,实现任意超GEO轨道的测量装备等天基观测点装备的高精度定位,从而实现对其他空间目标的定轨测量。
本领域技术人员可以理解,上述实施方式中各种方法的全部或部分功能可以通过硬件的方式实现,也可以通过计算机程序的方式实现。当上述实施方式中全部或部分功能通过计算机程序的方式实现时,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器、随机存储器、磁盘、光盘、硬盘等,通过计算机执行该程序以实现上述功能。例如,将程序存储在设备的存储器中,当通过处理器执行存储器中程序,即可实现上述全部或部分功能。另外,当上述实施方式中全部或部分功能通过计算机程序的方式实现时,该程序也可以存储在服务器、另一计算机、磁盘、光盘、闪存盘或移动硬盘等存储介质中,通过下载或复制保存到本地设备的存储器中,或对本地设备的***进行版本更新,当通过处理器执行存储器中的程序时,即可实现上述实施方式中全部或部分功能。
以上应用了具体个例对本发明进行阐述,只是用于帮助理解本发明,并不用以限制本发明。对于本发明所属技术领域的技术人员,依据本发明的思想,还可以做出若干简单推演、变形或替换。

Claims (10)

1.一种空间目标的测量装备定位方法,其特征在于,包括:
获取空间目标相对于测量装备的位置测量数据,所述位置测量数据为以所述测量装备为中心的第一天球坐标系下的位置测量数据;
将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据;
将所述测量装备的第一三维位置数据从所述第二天球坐标系转换至地球坐标系,得到在地球坐标系下的所述测量装备相对于所述空间目标的第二三维位置数据;
获取所述空间目标的精密星历,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据;
基于所述第二三维位置数据和所述第三三维位置数据,确定所述地球坐标系下的所述测量装备的位置定位数据。
2.如权利要求1所述的测量装备定位方法,其特征在于,所述第一天球坐标系以测量装备的质心为原点O1,以原点O1指向春分点方向为X1轴的正方向;Z1轴平行于地球自转轴,由地心指向北天极方向为Z1轴的正方向;Y1轴位于过原点的平面内,且Y1轴与X1轴、Z1轴构成右手直角坐标系。
3.如权利要求1所述的测量装备定位方法,其特征在于,所述第二天球坐标系以地心为原点O3,Z3轴平行于地球自转轴,由地心指向北天极方向为Z3轴的正方向;X3轴在赤道面上,由地心指向春分点方向为X3轴的正方向;Y3轴在赤道面上,且Y3轴与X3轴、Z3轴构成右手直角坐标系。
4.如权利要求1所述的测量装备定位方法,其特征在于,所述地球坐标系以地心为原点O4,由原点O4指向给定的协议地球极CIP的方向为Z4轴的正方向;X4轴在赤道平面内,由原点O4指向给定的经度零点为正方向;Y4轴位于赤道平面内,且Y4轴与X4轴、Z4轴构成右手直角坐标系。
5.如权利要求1所述的测量装备定位方法,其特征在于,将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据,包括:
构建过渡坐标系,所述过渡坐标系为以所述空间目标的质心为原点O2,以原点O2指向春分点方向为X2轴的正方向, Z2轴平行于地球自转轴,由地心指向北天极方向为Z2轴的正方向,Y2轴位于过原点的平面内,且Y2轴与X2轴、Z2轴构成右手直角坐标系;
将所述空间目标相对于测量装备的位置测量数据从第一天球坐标系转换到所述过渡坐标系,得到所述测量装备在所述过渡坐标系下的第四三维位置数据;
将所述测量装备在所述过渡坐标系下的第四三维位置数据转换到所述第二天球坐标系中,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据。
6.如权利要求1所述的测量装备定位方法,其特征在于,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据,包括:
在一插值区间内,对所述空间目标的精密星历内插任意时刻的空间目标,以获取预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据。
7.一种空间目标的测量装备定位***,其特征在于,包括:
空间目标,所述空间目标为已知精密星历的空间目标;
测量装备,用于观测并跟踪所述空间目标,以获取所述空间目标的图像数据;
搭载平台,用于搭载所述测量装备;
数据处理服务器,用于获取所述空间目标的图像数据,并根据所述空间目标的图像数据,确定所述空间目标相对于测量装备的位置测量数据;
所述数据处理服务器还用于将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据;将所述测量装备的第一三维位置数据从所述第二天球坐标系转换至地球坐标系,得到在地球坐标系下的所述测量装备相对于所述空间目标的第二三维位置数据;获取所述空间目标的精密星历,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据;基于所述第二三维位置数据和所述第三三维位置数据,确定所述地球坐标系下的所述测量装备的位置定位数据。
8.如权利要求7所述的测量装备定位***,其特征在于,将所述空间目标相对于测量装备的位置测量数据从以所述测量装备为中心的第一天球坐标系转换至以地心为中心的第二天球坐标系,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据,包括:
构建过渡坐标系,所述过渡坐标系为以所述空间目标的质心为原点,以原点O2指向春分点方向为X2轴的正方向,Z2轴平行于地球自转轴,由地心指向北天极方向为Z2轴的正方向,Y2轴位于过原点的平面内,且Y2轴与X2轴、Z2轴构成右手直角坐标系;
将所述空间目标相对于测量装备的位置测量数据从第一天球坐标系转换到所述过渡坐标系,得到所述测量装备在所述过渡坐标系下的第四三维位置数据;
将所述测量装备在所述过渡坐标系下的第四三维位置数据转换到所述第二天球坐标系中,得到所述测量装备在所述第二天球坐标系下的第一三维位置数据。
9.如权利要求7所述的测量装备定位***,其特征在于,根据所述空间目标的精密星历,得到预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据,包括:
在一插值区间内,对所述空间目标的精密星历内插任意时刻的空间目标,以获取预设观测时刻在地球坐标系下的所述空间目标的第三三维位置数据。
10.一种计算机可读存储介质,其特征在于,所述介质上存储有程序,所述程序能够被处理器执行以实现如权利要求1-6中任一项所述的方法。
CN202211092032.XA 2022-09-08 2022-09-08 空间目标的测量装备定位方法、***和存储介质 Active CN115200573B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211092032.XA CN115200573B (zh) 2022-09-08 2022-09-08 空间目标的测量装备定位方法、***和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211092032.XA CN115200573B (zh) 2022-09-08 2022-09-08 空间目标的测量装备定位方法、***和存储介质

Publications (2)

Publication Number Publication Date
CN115200573A true CN115200573A (zh) 2022-10-18
CN115200573B CN115200573B (zh) 2022-12-27

Family

ID=83573665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211092032.XA Active CN115200573B (zh) 2022-09-08 2022-09-08 空间目标的测量装备定位方法、***和存储介质

Country Status (1)

Country Link
CN (1) CN115200573B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115639849A (zh) * 2022-11-28 2023-01-24 中国人民解放军63921部队 一种机电复合的目标过顶跟踪方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004224171A (ja) * 2003-01-22 2004-08-12 Mitsubishi Electric Corp 人工衛星の姿勢決定装置
CN109655064A (zh) * 2018-11-27 2019-04-19 孙秀聪 一种基于相对转动的地球固连系-惯性系快速高精度转换方法
CN109813526A (zh) * 2018-12-29 2019-05-28 中国科学院紫金山天文台 一种基于天文定位的光学望远镜外场跟踪精度检测方法
CN113204917A (zh) * 2021-04-25 2021-08-03 中国科学院国家空间科学中心 针对geo目标的天基光学测角弧段初定轨及关联方法
CN114545374A (zh) * 2021-12-31 2022-05-27 中国空间技术研究院 一种星载激光雷达在轨定标方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004224171A (ja) * 2003-01-22 2004-08-12 Mitsubishi Electric Corp 人工衛星の姿勢決定装置
CN109655064A (zh) * 2018-11-27 2019-04-19 孙秀聪 一种基于相对转动的地球固连系-惯性系快速高精度转换方法
CN109813526A (zh) * 2018-12-29 2019-05-28 中国科学院紫金山天文台 一种基于天文定位的光学望远镜外场跟踪精度检测方法
CN113204917A (zh) * 2021-04-25 2021-08-03 中国科学院国家空间科学中心 针对geo目标的天基光学测角弧段初定轨及关联方法
CN114545374A (zh) * 2021-12-31 2022-05-27 中国空间技术研究院 一种星载激光雷达在轨定标方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨宏 等,: ""空基跟踪测量导弹的方法研究"", 《导弹与航天运载技术》 *
陆长明 等,: ""基于北斗卫星精密星历的中高轨望远镜测量精度标定方法"", 《光子学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115639849A (zh) * 2022-11-28 2023-01-24 中国人民解放军63921部队 一种机电复合的目标过顶跟踪方法及装置
CN115639849B (zh) * 2022-11-28 2023-05-02 中国人民解放军63921部队 一种机电复合的目标过顶跟踪方法及装置

Also Published As

Publication number Publication date
CN115200573B (zh) 2022-12-27

Similar Documents

Publication Publication Date Title
US9194954B2 (en) Method for geo-referencing an imaged area
CN105910607B (zh) 基于地面控制的卫星长周期姿态误差修正方法
JP6271047B2 (ja) 衛星信号受信特性推定装置、その方法及びプログラム
US8630804B2 (en) Method for geolocating an object by multitelemetry
CA2925850C (en) Target direction determination method and system
Mulawa On-orbit geometric calibration of the OrbView-3 high resolution imaging satellite
CN101743453A (zh) 任务后高精确度定位和定向***
CN103837150B (zh) 一种ccd天顶望远镜地面快速天文定位的方法
CN111398980B (zh) 一种机载LiDAR数据处理的方法及装置
Tang et al. Geometric accuracy analysis model of the Ziyuan-3 satellite without GCPs
CN115469308B (zh) 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质
CN115200573B (zh) 空间目标的测量装备定位方法、***和存储介质
CN114001756A (zh) 一种小视场星敏感器外场地面寻星方法
Pelc-Mieczkowska et al. GNSS obstacle mapping as a data preprocessing tool for positioning in a multipath environment
JP4423378B2 (ja) 静止衛星の位置座標表示方法およびそれを用いた座標表示装置
Seo et al. Kompsat-2 direct sensor modeling and geometric calibration/validation
Dawidowicz Antenna phase center variations corrections in processing of GPS observations with use of commercial software
Grm et al. Testing the functionality and applicability of smart devices for a handheld celestial navigation system
Ma et al. One-step solution to local tie vector determination at co-located GNSS/VLBI sites
El-Hallaq et al. Positional accuracy of the Google Earth imagery in the Gaza strip
Gupta et al. Long strip modelling for cartosat-1 with minimum control
Grimm GNSS antenna orientation based on modification of received signal strengths
RU2711834C1 (ru) Способ определения орбиты космического аппарата с аппаратурой для съёмки подстилающей поверхности
Yoo et al. Determination of the invariant point of the Korean VLBI network radio telescopes: First results at the Ulsan and Tamna observatories
Samaan et al. Compass star tracker for GPS-like applications

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant