CN116124153A - 一种对空天目标的双星共视定位方法及设备 - Google Patents
一种对空天目标的双星共视定位方法及设备 Download PDFInfo
- Publication number
- CN116124153A CN116124153A CN202310411789.9A CN202310411789A CN116124153A CN 116124153 A CN116124153 A CN 116124153A CN 202310411789 A CN202310411789 A CN 202310411789A CN 116124153 A CN116124153 A CN 116124153A
- Authority
- CN
- China
- Prior art keywords
- coordinates
- coordinate
- determining
- target
- optical
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000003287 optical effect Effects 0.000 claims abstract description 203
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000000630 rising effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 2
- 230000004807 localization Effects 0.000 claims 4
- 230000009977 dual effect Effects 0.000 claims 2
- 230000001174 ascending effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000002068 genetic effect Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Astronomy & Astrophysics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了对空天目标的双星共视定位方法及设备,包括:通过第一卫星设备和第二卫星设备在预设时刻各自获取目标的光学图像、设备在预设时刻的姿态参数、轨道运行参数和位置坐标;以及获取光学图像对应的光学传感器的参数;根据光学图像和光学传感器的参数,分别确定目标光学传感器的坐标系下的两个光学坐标;根据光学坐标、姿态参数和轨道运行参数,分别确定目标在J2000坐标系下的两个目标坐标;根据两个目标坐标和两个位置坐标,分别确定第一直线和第二直线;根据第一直线和第二直线,确定第一直线和第二直线之间的最短线段的端点坐标,根据端点坐标、第一位置坐标和第二位置坐标,确定目标在预设时刻的位置坐标。
Description
技术领域
本发明属于航天测量与控制技术领域,具体涉及一种对空天目标的双星共视定位方法及设备。
背景技术
当前,随着遥感卫星的作用越来越广泛,其飞行过程中观测获取的目标种类也越来越多样,其中,空天动目标(以下称为空天目标)就是观测的主要对象。这类目标涵盖范围广,包括人造卫星、飞机、恒星和再入空间目标等。随着空天目标的逐渐增多,其位置和运动速度也成为遥感领域越来越关注的信息之一。
通常,可以采用光学遥感卫星对空天目标进行定位,或者,采用地面的雷达和光学设备对空天目标进行定位。但是,光学遥感卫星对空天目标的探测属于被动探测,单颗卫星观测仅能得到目标的测角信息,无法得到距离信息,需要通过激光测距手段实现空天目标的定位。然而,激光受卫星的功率和频率的限制,使得对空天目标的定位精度较差。地面的雷达和光学设备跟踪快速目标的弧度较短,也难以得到空天目标的实时高精度的位置信息。
发明内容
为了解决相关技术中存在的上述问题,本发明提供了一种对空天目标的双星共视定位方法及设备。本发明要解决的技术问题通过以下技术方案实现:
本发明提供一种对空天目标的双星共视定位方法,包括:
分别获取预设时刻空天目标的第一光学图像与第二光学图像,以及所述第一光学图像对应的第一卫星设备在所述预设时刻的第一姿态参数、第一轨道运行参数和第一位置坐标,所述第二光学图像对应的第二卫星设备在所述预设时刻的第二姿态参数、第二轨道运行参数和第二位置坐标;
获取所述第一光学图像对应的第一光学传感器的参数,以及所述第二光学图像对应的第二光学传感器的参数;
根据所述第一光学图像和所述第一光学传感器的参数,确定所述空天目标在所述第一光学传感器的坐标系下的第一光学坐标,以及根据所述第二光学图像和所述第二光学传感器的参数,确定所述空天目标在所述第二光学传感器的坐标系下的第二光学坐标;
根据所述第一光学坐标、所述第一姿态参数和所述第一轨道运行参数,确定所述空天目标在J2000坐标系下的第一目标坐标,以及根据所述第二光学坐标、所述第二姿态参数和所述第二轨道运行参数,确定所述空天目标在J2000坐标系下的第二目标坐标;
根据所述第一目标坐标和所述第一位置坐标,确定第一直线,以及根据所述第二目标坐标和所述第二位置坐标,确定第二直线;
根据所述第一直线和所述第二直线,确定所述第一直线和所述第二直线之间的最短线段的端点坐标;
根据所述端点坐标、所述第一位置坐标和所述第二位置坐标,确定所述空天目标在所述预设时刻的位置坐标。
本发明具有如下有益技术效果:
通过使用第一卫星设备和第二卫星设备共同观测空天目标,并分别获得空天目标在第一卫星设备和第二卫星设备的光学传感器像平面的位置,通过坐标转换将第一卫星设备和第二卫星设备分别观测的空天目标的位置转换到J2000坐标系下,之后,结合第一卫星设备和第二卫星设备在观测到空天目标时的位置,确定出第一卫星设备观测空天目标的视向量所在的直线,以及第二卫星设备观测空天目标的视向量所在的直线,通过寻找出这两条异面直线的公垂线段,并根据公垂线段确定空天目标的位置;这样,通过双卫星设备进行一些参数的获取与坐标计算,便可获得空天目标的实时位置,不仅简化了定位过程,提高了定位效率,而且,由于位置是通过一系列精确计算获得,所以,还提高了定位精度,以及,由于卫星设备可以保持对空天目标的长时间的随动观测,所以,还提高了对空天目标的跟踪时长。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1为本发明实施例提供的对空天目标的双星共视定位方法的一个方法流程图;
图2为本发明实施例提供的对空天目标的双星共视定位方法的另一个流程图;
图3为本发明实施例提供的示例性的确定预设时刻的空天目标的位置的一个原理示意图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
在本发明的描述中,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。此外,本领域的技术人员可以将本说明书中描述的不同实施例或示例进行接合和组合。
尽管在此结合各实施例对本发明进行了描述,然而,在实施所要求保护的本发明过程中,本领域技术人员通过查看所述附图、公开内容、以及所附权利要求书,可理解并实现所述公开实施例的其他变化。在权利要求中,“包括”(comprising)一词不排除其他组成部分或步骤,“一”或“一个”不排除多个的情况。单个处理器或其他单元可以实现权利要求中列举的若干项功能。相互不同的从属权利要求中记载了某些措施,但这并不表示这些措施不能组合起来产生良好的效果。
图1是本发明实施例提供的对空天目标的双星共视定位方法的一个方法流程图,如图1所示,所述方法包括以下步骤:
S101、分别获取预设时刻空天目标的第一光学图像与第二光学图像,以及第一光学图像对应的第一卫星设备在预设时刻的第一姿态参数、第一轨道运行参数和第一位置坐标,第二光学图像对应的第二卫星设备在预设时刻的第二姿态参数、第二轨道运行参数和第二位置坐标。
本发明实施例中,空天目标可以是人造卫星、飞机、恒星、再入空间目标等,本发明实施例对此不作限定。
本发明实施例中,预设时刻可以是每个单位时刻,例如,每1分钟、每30秒等,也可以是预设的某个时间点,例如,10点30分等等,本发明实施例对此不作限定。
本发明实施例中,第一卫星设备具有第一光学传感器,第二卫星设备具有第二光学传感器,第一卫星设备和第二卫星设备在预设时刻可以同时通过自身的光学传感器采集同一空天目标的光学图像,分别得到第一光学图像和第二光学图像;同时,第一卫星设备在预设时刻还获取自身的第一姿态参数、第一轨道运行参数和第一位置坐标,第二卫星设备在预设时刻还获取自身的第二姿态参数、第二轨道运行参数和第二位置坐标。
这里,第一卫星设备和第二卫星设备处于不同位置。
这里,第一位置坐标是第一卫星设备在J2000坐标系下的位置坐标,同样的,第二位置坐标是第二卫星设备在J2000坐标系下的位置坐标。
这里,J2000坐标系又被称为J2000平赤道地心坐标系。该坐标系的原点是在地球质心,xy平面为J2000时刻的地球平赤道面,x轴指向J2000时刻的平春分点(即J2000时刻平赤道面与平黄道面的一个交点)。
S102、获取第一光学图像对应的第一光学传感器的参数,以及第二光学图像对应的第二光学传感器的参数。
本发明实施例中,第一卫星设备可以在获取到第一光学图像之后,或者在获取第一光学图像的同时,获取自身的第一光学传感器的参数,同样的,第二卫星设备可以在获取到第二光学图像之后,或者在获取第二光学图像的同时,获取自身的第二光学传感器的参数。
这里,第一光学传感器固定设置于第一卫星设备上,并且,第一光学传感器的视角的改变是通过调整第一卫星设备的飞行姿态而改变的,从而使得第一光学传感器的坐标系与第一卫星设备自身的坐标系一致。第二光学传感器同样固定设置于第二卫星设备上,并且,第二光学传感器的视角的改变也是通过调整第二卫星设备的飞行姿态而改变的,从而使得第二光学传感器的坐标系与第二卫星设备自身的坐标系一致。
这里,第一光学传感器的参数可以包括:第一光学传感器的焦距、第一光学传感器的像元个数,以及每个像元的尺寸;同样的,第二光学传感器的参数可以包括:第二光学传感器的焦距、第二光学传感器的像元个数,以及每个像元的尺寸。
这里,第一光学传感器的像元个数可以包括:第一光学传感器的坐标系中轴的像元个数和轴的像元个数。同样的,第二光学传感器的像元个数可以包括:第二光学传感器的坐标系中轴的像元个数和轴的像元个数。
S103、根据第一光学图像和第一光学传感器的参数,确定空天目标在第一光学传感器的坐标系下的第一光学坐标,以及根据第二光学图像和第二光学传感器的参数,确定空天目标在第二光学传感器的坐标系下的第二光学坐标。
本发明实施例中,第一光学图像中包含了空天目标的图像部分,还包含了除空天目标之外(背景)的图像部分,可以根据第一光学图像中包含的空天目标的图像部分,以及第一光学传感器的参数,确定出空天目标在第一光学传感器的坐标系下的第一光学坐标。同样的,第二光学图像中也包含了空天目标的图像部分,以及除空天目标之外(背景)的图像部分,可以根据第二光学图像中包含的空天目标的图像部分,以及第二光学传感器的参数,确定出空天目标在第二光学传感器的坐标系下的第二光学坐标。
在一些实施例中,上述S103中的根据第一光学图像和第一光学传感器的参数,确定空天目标在第一光学传感器的坐标系下的第一光学坐标,可以通过S1031~S1032实现:
S1031、根据第一光学图像,确定空天目标的第一质心位置坐标。
本发明实施例中,可以根据第一光学图像中包含的空天目标的图像部分,通过质心法计算出该空天目标的第一质心位置坐标。
示例性的,当第一光学图像通过公式(1)表示时,得到的第一质心位置坐标中的横坐标可以表示为公式(2),纵坐标可以表示为公式(3):
(1);
(2);
(3);
上式中,为第一光学图像,为第一光学图像中空天目标的图像部分,为第一光学图像中除空天目标之外的图像部分,为第一光学图像中每个像素点的坐标,1或0表示像素点对应的像素值。为第一质心位置坐标。
S1032、根据第一光学传感器的焦距、第一光学传感器的像元个数和每个像元的尺寸,将第一质心位置坐标转化为第一光学传感器的坐标系下的第一光学坐标。
本发明实施例中,可以根据第一光学传感器的焦距、第一光学传感器的坐标系的轴的像元个数和轴的像元个数,以及每个像元的尺寸,将第一质心位置坐标,转换为第一光学传感器的坐标系下的第一光学坐标。
示例性的,第一光学坐标可以通过公式(4)表示:
(4);
上式中,为第一光学坐标,为第一光学感器的坐标系的轴的像元个数,为第一光学传感器的坐标系的轴的像元个数,为第一光学传感器的每个像元的尺寸,为第一光学传感器的焦距。
在一些实施例中,上述S103中的根据第二光学图像和第二光学传感器的参数,确定空天目标在第二光学传感器的坐标系下的第二光学坐标,可以通过与S1031~S1032相同的原理实现。
示例性的,当第二光学图像通过公式(5)表示时,得到的第二质心位置坐标中的横坐标可以表示为公式(6),纵坐标可以表示为公式(7),以及得到的第二光学坐标可以表示为公式(8):
(5);
(6);
(7);
(8);
上式中,为第二光学图像,为第二光学图像中空天目标的图像部分,为第二光学图像中除空天目标之外的图像部分,为第二光学图像中每个像素点的坐标,1或0表示像素点对应的像素值;为第二光学感器的坐标系的轴的像元个数,为第二光学传感器的坐标系的轴的像元个数,为第二光学传感器的每个像元的尺寸,为第二光学传感器的焦距。为第二质心位置坐标。
S104、根据第一光学坐标、第一姿态参数和第一轨道运行参数,确定空天目标在J2000坐标系下的第一目标坐标,以及根据第二光学坐标、第二姿态参数和第二轨道运行参数,确定空天目标在J2000坐标系下的第二目标坐标。
本发明实施例中,第一轨道运行参数包括:第一卫星轨道参数和第一卫星地心向径;第二轨道运行参数包括:第二卫星轨道参数和第二卫星地心向径。
这里,第一卫星轨道参数可以是第一卫星设备所在的卫星轨道的基本参数,例如,可以是第一卫星设备的卫星轨道倾角、升交点赤经和升交点角距。第二卫星轨道参数同样可以是第二卫星设备所在的卫星轨道的基本参数,例如,可以是第二卫星设备的卫星轨道倾角、升交点赤经和升交点角距。
这里,第一卫星地心向径和第二卫星地心向径可以相同,例如,可以为地球半径。
在一些实施例中,如图2所示,上述S104中的根据第一光学坐标、第一姿态参数和第一轨道运行参数,确定空天目标在J2000坐标系下的第一目标坐标,可以通过S1041~S1045实现:
S1041、根据第一姿态参数,确定第一旋转矩阵。
示例性的,当第一姿态参数为第一卫星设备的欧拉角时,第一旋转矩阵可以表示为公式(9)~(11):
(9);
(10);
(11);
上式中,、和为第一欧拉角,、和为第一旋转矩阵。
S1042、根据第一旋转矩阵,将第一光学坐标转换为第一卫星设备的卫星轨道坐标系下的第一中间坐标。
这里,第一卫星设备的卫星轨道坐标系可以是UNW坐标系,该坐标系的原点为卫星质心,U轴沿卫星轨道的切线方向并以卫星运动方向为正,N轴沿卫星轨道的法线方向并以内法向为正,W轴沿轨道面的法线方向并以指向北极的方向为正。
示例性的,第一中间坐标可以通过公式(12)表示:
(12);
上式中,为第一中间坐标,为第一光学坐标,为的转置。
S1043、根据第一卫星地心向径和第二旋转矩阵,将第一中间坐标转换为地心轨道坐标系下的第一转换坐标。
本发明实施例中,第二旋转矩阵可以为和。
示例性的,第二旋转矩阵可以表示为公式(13)~(14),第一转换坐标可以表示为公式(15):
(13);
(14);
(15);
上式中,为第一转换坐标,为第一中间坐标,为的转置,为第一卫星地心向径。
S1044、根据第一卫星轨道参数,确定第三旋转矩阵。
本发明实施例中,第一卫星轨道参数可以包括:第一卫星轨道倾角、第一升交点赤经和第一升交点角距;第三旋转矩阵包括基于第一卫星轨道倾角确定的矩阵、基于第一升交点赤经确定的矩阵和基于第一升交点角距确定的矩阵。
示例性的,第三旋转矩阵可以通过公式(16)~(18)表示:
(16);
(17);
(18);
上式中,为第一升交点角距,为第一卫星轨道倾角,为第一升交点赤经,、和为第三旋转矩阵。
S1045、根据第三旋转矩阵,将第一转换坐标转换为J2000坐标系下的第一目标坐标。
示例性的,第一目标坐标可以通过公式(19)表示:
(19);
上式中,为第一目标坐标,为的转置。
在一些实施例中,上述S104中的根据第二光学坐标、第二姿态参数和第二轨道运行参数,确定空天目标在J2000坐标系下的第二目标坐标,同样可以通过与上述S1041~S1045相同的原理实现;其中,第二姿态参数为第二卫星设备的第二欧拉角、和;、和为第二卫星设备对应的第一旋转矩阵;第二卫星设备对应的第二旋转矩阵同样为和;第二卫星轨道参数可以包括:第二卫星轨道倾角、第二升交点赤经和第二升交点角距;第二卫星设备对应的第三旋转矩阵包括:基于第二卫星轨道倾角确定的矩阵、基于第二升交点赤经确定的矩阵和基于第二升交点角距确定的矩阵。
示例性的,第二卫星设备对应的第一旋转矩阵可以通过公式(20)~(22)表示,第二卫星设备对应的第二中间坐标可以通过公式(23)表示,第二卫星设备对应的第二转换坐标可以通过公式(24)表示,第二卫星设备对应的第三旋转矩阵可以通过公式(25)~(27)表示,以及第二卫星设备对应的第二目标坐标可以通过公式(28)表示:
(20);
(21);
(22);
(23);
(24);
(25);
(26);
(27);
(28);
上式中,、和为第二欧拉角,、和为第二卫星设备对应的第一旋转矩阵。为第二中间坐标,为第二光学坐标。为的转置。为第二中间坐标,为的转置。为第二卫星设备对应的第二卫星地心向径。为第二升交点角距,为第二卫星轨道倾角,为第二升交点赤经,、和为第二卫星设备对应的第三旋转矩阵。为第二转换坐标,为的转置。为第二目标坐标。
S105、根据第一目标坐标和第一位置坐标,确定第一直线,以及根据第二目标坐标和第二位置坐标,确定第二直线。
本发明实施例中,可以根据第一目标坐标和第一位置坐标,确定出第一目标坐标和第一位置坐标分别对应的坐标点所在的第一直线,以及根据第二目标坐标和第二位置坐标,确定出第二目标坐标和第二位置坐标分别对应的坐标点所在的第二直线。
示例性的,在第一目标坐标为,第一位置坐标为时,第一直线可以通过公式(29)表示:
(29);
上式中,为第一直线,、和为变量。
示例性的,在第二目标坐标为,第二位置坐标为时,第二直线可以通过公式(30)表示:
(30);
上式中,为第二直线,、和同样为变量。
S106、根据第一直线和第二直线,确定第一直线和第二直线之间的最短线段的端点坐标。
在一些实施例中,可以根据第一直线和第二直线上的点的坐标,通过最优化算法在第一直线上的点中寻找出第一目标点,在第二直线上的点中寻找出第二目标点,即第一目标点位于第一直线上,第二目标点位于第二直线上;其中,第一目标点与第二目标点构成的线段为最短线段;将第一目标点的坐标和第二目标点的坐标,作为最短线段的端点坐标。
这里,最优化算法例如可以为遗传算法、贪心算法等等,本发明实施例对此不作限定。
这里,第一直线和第二直线之间的最短线段即为第一直线和第二直线的公垂线段,端点坐标即为该公垂线段的垂足。
示例性的,当最优化算法为遗传算法,第一卫星设备确定的空天目标的位置为,第二卫星设备确定的空天目标的位置为,以及第一卫星设备的位置为(即第一卫星设备的第一光学传感器的坐标系的原点),第二卫星设备的位置为(即第二卫星设备的第二光学传感器的坐标系的原点)时,上述S106可以通过以下步骤实现:
S1、初始种群选取:将第一直线中由指向的射线L1分为n段,并在每段选取一个点的坐标,得到n个点,以及,将第二直线中由指向的射线L2也分为n段,并在每段选取一个点的坐标,同样得到n个点;并将L1与L2之间的最短线段作为目标函数;采用二进制编码对每个点的坐标进行编码,L1和L2分别对应得到n个二进制串,将得到的二进制串作为遗传算法的运算对象;
S2、选择当前射线的当前种群:对于L1对应的n个点中的任意一个点,根据点的坐标计算该点到L2对应的n个点之间的距离,将得到的n个距离值中的最小距离值的倒数,作为该点的适应度,如此,得到L1对应的n个点中每个点的适应度;根据适应度,对L1对应的n个点进行由高到低的排序,并将排序位于前预设百分比(例如,前20%)的点作为选择对象,将排序位于后预设百分比(例如,后20%)的点淘汰;之后,对保留的排序位于前预设百分比的点,继续采用S2步骤的原理进行一次选择,得到L1当前对应的m个点;
S3、对当前射线的当前种群进行交叉:将L1当前对应的m个点的二进制串进行两两随机组合,并按预设交叉概率(,为预设交叉概率,为这m个点中的最大适应度,为这m个点的平均适应度,为进行交叉的每个二进制串的组合中的最大适应度,和均为预设常数,例如,均为小于0.1的预设常数)进行交叉,其中,交叉位随机选择,从而重新形成m个点。计算重新形成的这m个点与L2上的n个点之间的最短距离,并选择出得到的最短距离中的最小值;将交叉后得到的最短距离的最小值与交叉前的最短距离的最小值进行比较,如果交叉后得到的最短距离的最小值小于交叉前的最短距离的最小值,则交叉成功;否则,交叉失败,并重新对L1当前对应的m个点进行交叉,直到交叉成功;
S4、对交叉成功后的个体进行变异:以预设变异概率(例如,0.01)从L1的交叉成功后得到的m个点中确定要变异的点,将要变异的点的二进制串中的任意一位的值进行改变,得到变异后的点,如此,得到L1对应的更新的m个点;
S5、选择当前射线的当前种群:对于L2对应的n个点中的任意一个点,根据点的坐标计算该点到L1对应的更新的m个点之间的距离,将得到的m个距离值中的最小距离值的倒数,作为该点的适应度,如此,得到L2对应的n个点中每个点的适应度;根据适应度,对L2对应的n个点进行由高到低的排序,并将排序位于前预设百分比(例如,前20%)的点作为选择对象,将排序位于后预设百分比(例如,后20%)的点淘汰;之后,对保留的排序位于前预设百分比的点,继续采用S5步骤的原理进行一次选择,得到L2当前对应的m个点;
S6、对当前射线的当前种群进行交叉:将L2当前对应的m个点的二进制串进行两两随机组合,并按预设交叉概率进行交叉,其中,交叉位随机选择,从而重新形成m个点。计算这m个点各自与L1对应的更新的m个点之间的最短距离,并选择出得到的最短距离中的最小值;将交叉后得到的最短距离的最小值与交叉前的最短距离的最小值进行比较,如果交叉后得到的最短距离的最小值小于交叉前的最短距离的最小值,则交叉成功;否则,交叉失败,并重新对L2当前对应的m个点进行交叉,直到交叉成功;
S7、对交叉成功后的个体进行变异:以预设变异概率(例如,0.01)从L2的交叉成功后得到的m个点中确定要变异的点,将要变异的点的二进制串中的任意一位的值进行改变,得到变异后的点,如此,得到L2对应的更新的m个点;
S8、采用与上述S5~S7相同的原理,继续对L1对应的更新的m个点进行处理,如此循环迭代,直至连续多次(例如,连续两次或三次)得到的交叉后得到的最短距离的最小值相同时,将该最小值作为L1对应的第一直线与L2对应的第二直线之间的最短距离,并将该最小值对应的两个点的坐标,作为第一直线与第二直线之间的最短距离的端点坐标。
S107、根据端点坐标、第一位置坐标和第二位置坐标,确定空天目标在预设时刻的位置坐标。
本发明实施例中,最短线段由第一目标点和第二目标点构成,端点坐标包括:第一目标点的坐标和第二目标点的坐标。基于此,可以确定第二目标点的坐标与第一目标点的坐标之间的坐标差值,以及根据第一目标点的坐标、第二目标点的坐标、第一位置坐标和第二位置坐标,确定权重系数;之后,根据第一目标点的坐标、坐标差值和权重系数,确定出空天目标在预设时刻的位置坐标。
在一些实施例中,可以根据第一目标点的坐标、第二目标点的坐标和第一位置坐标,确定第一距离;根据第一目标点的坐标、第二目标点的坐标和第二位置坐标,确定第二距离;确定第一距离与第二距离之和,得到总距离;根据第一距离与总距离,确定权重系数。
示例性的,当第一目标点的坐标为、第二目标点的坐标为、第一位置坐标为、第二位置坐标为时,第一距离可以通过公式(31)表示,第二距离可以通过公式(32)表示,权重系数可以通过公式(33)表示,具体如下:
(31);
(32);
(33);
上式中,为第一距离,为第二距离,为权重系数。
这里,第一位置坐标对应第一位置、第二位置坐标对应第二位置。在一些实施例中,还可以根据第一目标点的坐标和第二目标点的坐标,确定最短线段的中点的坐标;根据第一位置坐标和第一目标点的坐标,确定由第一位置指向第一目标点的第一向量;根据第一位置坐标和中点的坐标,确定由第一位置指向中点的第二向量;确定第一向量与第二向量之间的第一夹角;根据第二位置坐标和第二目标点的坐标,确定由第二位置指向第二目标点的第三向量;根据第二位置坐标和中点的坐标,确定由第二位置指向中点的第四向量;确定第三向量与第四向量之间的第二夹角;根据第一夹角与第二夹角,确定权重系数。
示例性的,当第一目标点的坐标为、第二目标点的坐标为时,最短线段的中点的坐标为。例如,如图3所示,在J2000坐标系中,在预设时刻时,第一卫星设备确定的空天目标的位置为,第二卫星设备确定的空天目标的位置为,以及第一卫星设备的位置为(即第一卫星设备的第一光学传感器的坐标系的原点),第二卫星设备的位置为(即第二卫星设备的第二光学传感器的坐标系的原点),则由和确定的第一直线与由和确定的第二直线之间的最短线段为 ,最短线段 的中点为。
示例性的,当第一目标点为、第二目标点为、第一位置坐标对应的第一位置为、第二位置坐标对应的第二位置为、最短线段的中点为时,第一向量为、第二向量为、第三向量为、第四向量为,并且,第一向量可以通过公式(34)表示、第二向量可以通过公式(35)表示,第一夹角可以通过(36)表示,权重系数可以通过公式(37)表示,具体如下:
(34);
(35);
(36);
(37);
其中,为第一向量,为第二向量,为第一夹角,为第二夹角,为权重系数。
这里,确定第三向量的原理可以与上述公式(34)所示的原理相同,确定第四向量的原理可以与上述公式(35)所示的原理相同,确定第二夹角的原理可以与上述公式(36)所示的原理相同。
本发明实施例,通过双卫星设备进行一些参数的获取与坐标计算,便可获得空天目标的实时位置,不仅简化了定位过程,提高了定位效率,而且,由于位置是通过一系列精确计算获得,所以,还提高了定位精度,以及,由于卫星设备可以保持对空天目标的长时间的随动观测,所以,还提高了对空天目标的跟踪时长。
本发明实施例还提供一种对空天目标的双星共视定位设备,包括:第一卫星设备和第二卫星设备。第一卫星设备,用于获取预设时刻空天目标的第一光学图像,以及所述预设时刻时自身的第一姿态参数、自身的第一轨道运行参数和自身的第一位置坐标;获取所述第一光学图像对应的第一光学传感器的参数;第二卫星设备,用于获取预设时刻时所述空天目标的第二光学图像,以及所述预设时刻时自身的第二姿态参数、自身的第二轨道运行参数和自身的第二位置坐标;获取所述第二光学图像对应的第二光学传感器的参数;所述第一卫星设备或所述第二卫星设备,还用于根据所述第一光学图像和所述第一光学传感器的参数,确定所述空天目标在所述第一光学传感器的坐标系下的第一光学坐标,以及根据所述第二光学图像和所述第二光学传感器的参数,确定所述空天目标在所述第二光学传感器的坐标系下的第二光学坐标;根据所述第一光学坐标、所述第一姿态参数和所述第一轨道运行参数,确定所述空天目标在J2000坐标系下的第一目标坐标,以及根据所述第二光学坐标、所述第二姿态参数和所述第二轨道运行参数,确定所述空天目标在J2000坐标系下的第二目标坐标;根据所述第一目标坐标和所述第一位置坐标,确定第一直线,以及根据所述第二目标坐标和所述第二位置坐标,确定第二直线;根据所述第一直线和所述第二直线,确定所述第一直线和所述第二直线之间的最短线段的端点坐标;根据所述端点坐标、所述第一位置坐标和所述第二位置坐标,确定所述空天目标在所述预设时刻的位置坐标。
在一些实施例中,第一卫星设备和第二卫星设备中的一个卫星设备可以在获得预设时刻自身的姿态参数、轨道运行参数、位置坐标,以及自身的光学传感器的参数后,将这些数据发送至另一个卫星设备,由该另一个卫星设备根据自身的预设时刻自身的姿态参数、轨道运行参数、位置坐标,自身的光学传感器的参数,以及接收到的这些数据,计算出空天目标在预设时刻的位置坐标。
在一些实施例中,第一卫星设备和第二卫星设备也可以各自根据获取的参数计算出目标坐标和自身的位置坐标,之后,由第一卫星设备和第二卫星设备中的一个卫星设备将计算出的目标坐标和自身的位置坐标,发送至另一个卫星设备,并由该另一个设备根据自己计算出的目标坐标和自身的位置坐标,以及接收到的目标坐标和位置坐标,计算出空天目标在预设时刻的位置坐标。
这里,第一卫星设备或第二卫星设备计算出空天目标在预设时刻的位置坐标的原理具体如上述方法部分所述,此处不再说明。
本发明基于双星观测同一目标场景,解算空间目标的位置信息,对空间目标的定位和跟踪具有很强的应用价值,同时也对光学遥感卫星星座设计具有指导意义。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (10)
1.一种对空天目标的双星共视定位方法,其特征在于,包括:
分别获取预设时刻空天目标的第一光学图像与第二光学图像,以及所述第一光学图像对应的第一卫星设备在所述预设时刻的第一姿态参数、第一轨道运行参数和第一位置坐标,所述第二光学图像对应的第二卫星设备在所述预设时刻的第二姿态参数、第二轨道运行参数和第二位置坐标;
获取所述第一光学图像对应的第一光学传感器的参数,以及所述第二光学图像对应的第二光学传感器的参数;
根据所述第一光学图像和所述第一光学传感器的参数,确定所述空天目标在所述第一光学传感器的坐标系下的第一光学坐标,以及根据所述第二光学图像和所述第二光学传感器的参数,确定所述空天目标在所述第二光学传感器的坐标系下的第二光学坐标;
根据所述第一光学坐标、所述第一姿态参数和所述第一轨道运行参数,确定所述空天目标在J2000坐标系下的第一目标坐标,以及根据所述第二光学坐标、所述第二姿态参数和所述第二轨道运行参数,确定所述空天目标在J2000坐标系下的第二目标坐标;
根据所述第一目标坐标和所述第一位置坐标,确定第一直线,以及根据所述第二目标坐标和所述第二位置坐标,确定第二直线;
根据所述第一直线和所述第二直线,确定所述第一直线和所述第二直线之间的最短线段的端点坐标;
根据所述端点坐标、所述第一位置坐标和所述第二位置坐标,确定所述空天目标在所述预设时刻的位置坐标。
2.根据权利要求1所述的对空天目标的双星共视定位方法,其特征在于,所述最短线段由第一目标点和第二目标点构成,所述端点坐标包括:所述第一目标点的坐标和所述第二目标点的坐标;
所述根据所述端点坐标、所述第一位置坐标和所述第二位置坐标,确定所述空天目标在所述预设时刻的位置坐标,包括:
确定所述第二目标点的坐标与所述第一目标点的坐标之间的坐标差值;
根据所述第一目标点的坐标、所述第二目标点的坐标、所述第一位置坐标和所述第二位置坐标,确定权重系数;
根据所述第一目标点的坐标、所述坐标差值和所述权重系数,确定出所述空天目标在所述预设时刻的位置坐标。
3.根据权利要求2所述的对空天目标的双星共视定位方法,其特征在于,所述第一目标点位于所述第一直线上,所述第二目标点位于所述第二直线上;所述根据所述第一目标点的坐标、所述第二目标点的坐标、所述第一位置坐标和所述第二位置坐标,确定权重系数,包括:
根据所述第一目标点的坐标、所述第二目标点的坐标和所述第一位置坐标,确定第一距离;
根据所述第一目标点的坐标、所述第二目标点的坐标和所述第二位置坐标,确定第二距离;
确定所述第一距离与所述第二距离之和,得到总距离;
根据所述第一距离与所述总距离,确定所述权重系数。
4.根据权利要求2所述的对空天目标的双星共视定位方法,其特征在于,所述第一目标点位于所述第一直线上,所述第二目标点位于所述第二直线上;所述第一位置坐标对应第一位置,所述第二位置坐标对应第二位置;所述根据所述第一目标点的坐标、所述第二目标点的坐标、所述第一位置坐标和所述第二位置坐标,确定权重系数,包括:
根据所述第一目标点的坐标和所述第二目标点的坐标,确定所述最短线段的中点的坐标;
根据所述第一位置坐标和所述第一目标点的坐标,确定由所述第一位置指向所述第一目标点的第一向量;
根据所述第一位置坐标和所述中点的坐标,确定由所述第一位置指向所述中点的第二向量;
确定所述第一向量与所述第二向量之间的第一夹角;
根据所述第二位置坐标和所述第二目标点的坐标,确定由所述第二位置指向所述第二目标点的第三向量;
根据所述第二位置坐标和所述中点的坐标,确定由所述第二位置指向所述中点的第四向量;
确定所述第三向量与所述第四向量之间的第二夹角;
根据所述第一夹角与所述第二夹角,确定所述权重系数。
5.根据权利要求1所述的对空天目标的双星共视定位方法,其特征在于,所述第一轨道运行参数包括:第一卫星轨道参数和第一卫星地心向径;所述根据所述第一光学坐标、所述第一姿态参数和所述第一轨道运行参数,确定所述空天目标在J2000坐标系下的第一目标坐标,包括:
根据所述第一姿态参数,确定第一旋转矩阵;
根据所述第一旋转矩阵,将所述第一光学坐标转换为所述第一卫星设备的卫星轨道坐标系下的第一中间坐标;
根据所述第一卫星地心向径和第二旋转矩阵,将所述第一中间坐标转换为地心轨道坐标系下的第一转换坐标;
根据所述第一卫星轨道参数,确定第三旋转矩阵;
根据所述第三旋转矩阵,将所述第一转换坐标转换为J2000坐标系下的所述第一目标坐标。
6.根据权利要求5所述的对空天目标的双星共视定位方法,其特征在于,所述第一姿态参数为第一欧拉角;所述第一中间坐标表示为:
;
;
;
;
其中,、和为所述第一欧拉角,、和为所述第一旋转矩阵,为所述第一光学坐标,为的转置。
7.根据权利要求5所述的对空天目标的双星共视定位方法,其特征在于,所述第一转换坐标表示为:
;
;
;
其中,为所述第一中间坐标,为的转置,和为所述第二旋转矩阵,为所述第一卫星地心向径。
8.根据权利要求5所述的对空天目标的双星共视定位方法,其特征在于,所述第一卫星轨道参数包括:第一卫星轨道倾角、第一升交点赤经和第一升交点角距;所述第一目标坐标表示为:
;
;
;
;
其中,为所述第一转换坐标,为的转置,为所述第一升交点角距,为所述第一卫星轨道倾角,为所述第一升交点赤经,、和为所述第三旋转矩阵。
9.根据权利要求1所述的对空天目标的双星共视定位方法,其特征在于,所述第一光学传感器的参数包括:所述第一光学传感器的焦距、所述第一光学传感器的像元个数和每个像元的尺寸;所述根据所述第一光学图像和所述第一光学传感器的参数,确定所述空天目标在所述第一光学传感器的坐标系下的第一光学坐标,包括:
根据所述第一光学图像,确定所述空天目标的第一质心位置坐标;
根据所述第一光学传感器的焦距、所述第一光学传感器的像元个数和所述每个像元的尺寸,将所述第一质心位置坐标转化为所述第一光学传感器的坐标系下的所述第一光学坐标。
10.一种对空天目标的双星共视定位设备,其特征在于,包括:
第一卫星设备,用于获取预设时刻空天目标的第一光学图像,以及所述预设时刻时自身的第一姿态参数、自身的第一轨道运行参数和自身的第一位置坐标;获取所述第一光学图像对应的第一光学传感器的参数;
第二卫星设备,用于获取预设时刻时所述空天目标的第二光学图像,以及所述预设时刻时自身的第二姿态参数、自身的第二轨道运行参数和自身的第二位置坐标;获取所述第二光学图像对应的第二光学传感器的参数;
所述第一卫星设备或所述第二卫星设备,还用于根据所述第一光学图像和所述第一光学传感器的参数,确定所述空天目标在所述第一光学传感器的坐标系下的第一光学坐标,以及根据所述第二光学图像和所述第二光学传感器的参数,确定所述空天目标在所述第二光学传感器的坐标系下的第二光学坐标;根据所述第一光学坐标、所述第一姿态参数和所述第一轨道运行参数,确定所述空天目标在J2000坐标系下的第一目标坐标,以及根据所述第二光学坐标、所述第二姿态参数和所述第二轨道运行参数,确定所述空天目标在J2000坐标系下的第二目标坐标;根据所述第一目标坐标和所述第一位置坐标,确定第一直线,以及根据所述第二目标坐标和所述第二位置坐标,确定第二直线;根据所述第一直线和所述第二直线,确定所述第一直线和所述第二直线之间的最短线段的端点坐标;根据所述端点坐标、所述第一位置坐标和所述第二位置坐标,确定所述空天目标在所述预设时刻的位置坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310411789.9A CN116124153B (zh) | 2023-04-18 | 2023-04-18 | 一种对空天目标的双星共视定位方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310411789.9A CN116124153B (zh) | 2023-04-18 | 2023-04-18 | 一种对空天目标的双星共视定位方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116124153A true CN116124153A (zh) | 2023-05-16 |
CN116124153B CN116124153B (zh) | 2023-06-16 |
Family
ID=86308534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310411789.9A Active CN116124153B (zh) | 2023-04-18 | 2023-04-18 | 一种对空天目标的双星共视定位方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116124153B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1305772C (en) * | 1988-08-23 | 1992-07-28 | Communications Satellite Corporation | Method of orienting a synchronous satellite |
WO2000023817A1 (de) * | 1998-10-15 | 2000-04-27 | Astrium Gmbh | Verfahren und vorrichtung zur lagebestimmung von kommunikations-satelliten |
CN104573251A (zh) * | 2015-01-16 | 2015-04-29 | 航天东方红卫星有限公司 | 一种星载光学遥感器全视场表观光谱辐亮度确定方法 |
US20160259044A1 (en) * | 2013-01-04 | 2016-09-08 | National Central University | Three-dimensional positioning method |
CN111609794A (zh) * | 2020-04-24 | 2020-09-01 | 北京控制工程研究所 | 一种基于双机械臂抓捕的目标卫星星箭对接环抓捕点定位方法 |
CN114942431A (zh) * | 2022-04-29 | 2022-08-26 | 西安电子科技大学 | 光雷同步观测的空间目标瞬时姿态估计方法 |
CN115183669A (zh) * | 2022-06-01 | 2022-10-14 | 中国空间技术研究院 | 一种基于卫星图像的目标定位方法 |
CN115793002A (zh) * | 2022-11-22 | 2023-03-14 | 中国空间技术研究院 | 基于测向误差权重的双星联合无源定位方法 |
CN115880328A (zh) * | 2022-09-29 | 2023-03-31 | 北京空间机电研究所 | 一种双星协同观测的空中运动目标实时定位方法 |
-
2023
- 2023-04-18 CN CN202310411789.9A patent/CN116124153B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1305772C (en) * | 1988-08-23 | 1992-07-28 | Communications Satellite Corporation | Method of orienting a synchronous satellite |
WO2000023817A1 (de) * | 1998-10-15 | 2000-04-27 | Astrium Gmbh | Verfahren und vorrichtung zur lagebestimmung von kommunikations-satelliten |
US20160259044A1 (en) * | 2013-01-04 | 2016-09-08 | National Central University | Three-dimensional positioning method |
CN104573251A (zh) * | 2015-01-16 | 2015-04-29 | 航天东方红卫星有限公司 | 一种星载光学遥感器全视场表观光谱辐亮度确定方法 |
CN111609794A (zh) * | 2020-04-24 | 2020-09-01 | 北京控制工程研究所 | 一种基于双机械臂抓捕的目标卫星星箭对接环抓捕点定位方法 |
CN114942431A (zh) * | 2022-04-29 | 2022-08-26 | 西安电子科技大学 | 光雷同步观测的空间目标瞬时姿态估计方法 |
CN115183669A (zh) * | 2022-06-01 | 2022-10-14 | 中国空间技术研究院 | 一种基于卫星图像的目标定位方法 |
CN115880328A (zh) * | 2022-09-29 | 2023-03-31 | 北京空间机电研究所 | 一种双星协同观测的空中运动目标实时定位方法 |
CN115793002A (zh) * | 2022-11-22 | 2023-03-14 | 中国空间技术研究院 | 基于测向误差权重的双星联合无源定位方法 |
Non-Patent Citations (1)
Title |
---|
金烨;郭鑫;汪精华;: "星载光学图像***目标定位算法研究", 兵工学报, no. 2 * |
Also Published As
Publication number | Publication date |
---|---|
CN116124153B (zh) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3321888B1 (en) | Projected image generation method and device, and method for mapping image pixels and depth values | |
US10704863B1 (en) | System for tracking a presumed target using network-connected lead and follower scopes, and scope for configured for use in the system | |
CN106871932B (zh) | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 | |
CN109443359B (zh) | 一种地面全景影像的地理定位方法 | |
EP3287736B1 (en) | Dynamic, persistent tracking of multiple field elements | |
CN103837150A (zh) | 一种ccd天顶望远镜地面快速天文定位的方法 | |
JP2023539615A (ja) | 宇宙空間における空間飛行体の姿勢決定のための方法、装置及びコンピュータプログラム製品 | |
Qiao et al. | Ground target geolocation based on digital elevation model for airborne wide-area reconnaissance system | |
Zhu et al. | Vision navigation for aircrafts based on 3D reconstruction from real-time image sequences | |
US10317218B2 (en) | Planetary surveillance system | |
CN116124153B (zh) | 一种对空天目标的双星共视定位方法及设备 | |
Willner et al. | New astrometric observations of Phobos with the SRC on Mars Express | |
Sheshadri et al. | Position estimation by registration to planetary terrain | |
CN116309798A (zh) | 一种无人机成像定位方法 | |
Nadezhdina et al. | A libration model for Enceladus based on geodetic control point network analysis | |
CN104458653A (zh) | 一种测量大天顶距处的大气折射值的方法和*** | |
Rawashdeh et al. | Design of a stellar gyroscope for visual attitude propagation for small satellites | |
CN109059866B (zh) | 基于图像的行星近景摄影测量相机安装参数的测量方法 | |
Kaiser et al. | Position and orientation of an aerial vehicle through chained, vision-based pose reconstruction | |
Palmer et al. | Mercator—Independent rover localization using stereophotoclinometry and panoramic images | |
Li et al. | Rapid star identification algorithm for fish-eye camera based on PPP/INS assistance | |
CN113348765B (zh) | 基于蒙特卡洛法的星敏感器全天识别星库构建方法 | |
Li et al. | IMU-Assisted Uncertainty-Weighted Attitude Estimation Algorithm from Non-Corresponding Points | |
Bai et al. | Multitarget location capable of adapting to complex geomorphic environment for the airborne photoelectric reconnaissance system | |
Zhu et al. | A landmark-based navigation method for autonomous aircraft |
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 |