CN114998423B - 一种望远镜自主定向方法 - Google Patents
一种望远镜自主定向方法 Download PDFInfo
- Publication number
- CN114998423B CN114998423B CN202210639586.0A CN202210639586A CN114998423B CN 114998423 B CN114998423 B CN 114998423B CN 202210639586 A CN202210639586 A CN 202210639586A CN 114998423 B CN114998423 B CN 114998423B
- Authority
- CN
- China
- Prior art keywords
- star
- calibration
- telescope
- stars
- matching
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000000007 visual effect Effects 0.000 claims abstract description 14
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 238000013316 zoning Methods 0.000 claims description 2
- 238000003491 array Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000005192 partition Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
- G01C15/002—Active optical surveying means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明属于空间目标观测技术领域,具体涉及一种望远镜自主定向方法。本发明通过获取天空图像、剔除非恒星目标、选取定标星、匹配标准星、望远镜视轴中心指向解算实现望远镜的自主定向。在剔除非恒星目标环节,通过长短曝光图像结合的方式,剔除了图像中的卫星、撞星和其他假目标情况,提高了匹配效率和准确性;采用“回”形区域分区的方式确保了定标星之间的角距尽量大,以减少匹配冗余;通过角距进行第一轮匹配,有匹配冗余时,继续通过星等匹配进一步消除冗余。在望远镜视轴中心指向解算环节,通过已知标准星的理想坐标和像面坐标建立望远镜测量模型,在望远镜站址精确天文坐标未知的情况下,依然可以解算出望远镜视轴中心指向。
Description
技术领域
本发明属于空间目标观测技术领域,具体涉及一种望远镜自主定向方法。
背景技术
望远镜在进行精密观测之前,必须进行方位方向基准的标定,确定方位零位,也就是定向。通常的做法是在望远镜周边架设多个方位标,通过测量已知大地坐标的方位标实现定向。在夜间利用北极星进行定向也是一种方案。但以上两种方案都具有一定的局限性,
对于地面机动望远镜,随着机动站布站位置的变化,跟随望远镜布设方位标的方案不能完全实现,给望远镜的使用带来了很多的不便,使观测无法顺利进行。在低纬度地区因不易观察到北极星也会导致定向失败。同时,对于望远镜精确站址不能提前得知的情况,以上两种方法均不能实现精确定向。
发明名称为《一种光电经纬仪快速自主定向的方法》提供了一种通过多颗恒星进行标定的方法,但该方法没有提供星图匹配的具体方法,同时该专利仅适用于观测站天文经纬度精确已知的情况,对于船载或运动测站等站址坐标不能精确得知的望远镜,无法使用该定向方法。
发明内容
本发明提供了一种可以全自动实现的望远镜自主定向方法,即使在机动望远镜天文经纬度未知的情况下也可实现。本方法通过获取天空图像、剔除非恒星目标、选取定标星、匹配标准星、望远镜视轴中心指向解算等步骤,可快速实现望远镜的自主定向。在剔除非恒星目标环节,通过长短曝光图像结合的方式,剔除了图像中的卫星、撞星和其他假目标情况,提高了匹配效率和准确性;在定标星选取环节,采用“回”形区域分区的方式确保了定标星之间的角距尽量大,以减少匹配冗余;在恒星匹配环节,通过角距进行第一轮匹配,有匹配冗余时,继续通过星等匹配进一步消除冗余。在望远镜视轴中心指向解算环节,在将定标星和标准星匹配后,通过已知定标星的理想坐标和像面坐标建立望远镜测量模型,在望远镜站址精确天文坐标未知的情况下,依然可以解算出望远镜视轴中心指向。
本发明的技术方案是:
一种望远镜自主定向方法,包括以下步骤:
S1、获取天空图像:望远镜指向一个天区静止不动,连续获取长曝光天空背景图像和短曝光天空背景图像,分别统计两张图像上所有目标的像元数、像面坐标、灰度值,并记录曝光中心时刻;
S2、剔除非恒星目标:
统计长曝光图像中像元数大于阈值的目标的像面坐标,在短曝光图像中扣除以上目标,可剔除图像中的卫星、撞星等目标,剩下恒星目标;
对比长曝光图像和短曝光图像,若短曝光图像中有目标的位置,但长曝光图像中该位置没有目标,在短曝光图像中扣除该目标,可剔除图像中的假目标;
S3、将短曝光图像分割成“回”形区域,按照星像的平均灰度从大到小在四个“一”字型区域中的三个区域各选择一颗星作为定标星,分别记为N1、N2、N3,三个定标星俯仰角从高到低依次为ξ1、ξ2、ξ3,根据定标星的像面坐标计算出三颗星两两之间的角距,为σ12、σ23、σ13;
S4、检索恒星星库中对应赤纬从高到低分别为ξ1′、ξ2′、ξ3′的恒星作为标准星,构建出三个顶点分别位于纬度ξ1′、ξ2′、ξ3′上的所有三角形,计算出三角形的角距,形成配对用标准数据库,设用于配对的恒星编号分别为M1、M2、M3,三颗星的理论角距依次排列为γ12、γ23、γ13;所述恒星星库的构建方式为:将设定纬度范围内的恒星,按照纬度区域增加的方式进行分区,在各纬度区域,按星等增加的方式排列,其中可以根据望远镜的视场调节星等,视场越大,星库中保留的恒星星等就越低,可以提高匹配效率;
对定标星和标准星进行按以下条件进行匹配,确定出三颗定标星对应的标准星:
条件1:
min(σ12、σ23、σ13)≥min(γ12、γ23、γ13)+∈1
max(σ12、σ23、σ13)≤max(γ12、γ23、γ13)+∈1
条件2:
其中,∈1是设定误差门限1,∈2是设定误差门限2;
当符合要求的标准星超过一组时,根据星等进行再次匹配,将每组标准星按照星等从亮到暗排序,对应的赤纬分别为ξMmax′、ξMmid′、ξMmin′,对其两两做差,记录差值的正负性;对应的定标星按照星象平均灰度从高到低排序,对应的俯仰角分别为ξMmax、ξMmid、ξMmin,对其两两做差,记录差值的正负性,两组记录正负性的数组完全一致时,认为匹配成功;匹配成功后,通过查询定标星对应标准星在星库中的数据,得出定标星的赤经、赤纬。
S5、对于已知站点天文经纬度的测站,直接根据第i颗定标星的赤经、赤纬和地方时t求出定标星相对测站的理论方位角,再根据理论方位角和实测方位角之间的偏差ΔAi对方位编码器进行置数,实现定向;对于站点精确天文经纬度未知的测站,先根据站址的地理经纬度粗略估算出一个零位,再根据望远镜的误差模型进行修正,得出望远镜的精确零位,具体为:
根据测站粗略地理经纬度站址和第i颗定标星的赤经、赤纬(αi,δi)和地方时t粗略估算出测站中心点相对正北的偏差,为了提高精度,对多颗定标星的偏差值求平均值:
ΔAi=Ai-Ai'
式中:Ai'、Ai为第i颗定标星相对测站的理论方位角和实测方位角;
计算得出定标星的理想坐标(Xi,Yi):
式中:(αi,δi)为第i颗定标星的赤经、赤纬,(α0,δ0)为视轴中心指向的赤经、赤纬,(A0、E0)为视轴中心对应的方位俯仰值,(Ai、Ei)为第i颗定标星对应的方位俯仰值;
将定标星的理想坐标(Xi,Yi)代入六常数模型求出望远镜的测量模型,通过最小二乘法,求出望远镜中心偏差等误差:
式中:(xi,yi)为匹配成功的定标星的像面坐标,
修正望远镜中心指向,实现精确定向A0z:
A0z=A0+ΔAi+a
A0为视轴中心对应的方位,a为六常数模型系数a
本发明的有益效果是:提供了一种可以全自动实现的机动望远镜自主定向的方法,即使在机动望远镜精确天文经纬度坐标未知的情况下也可实现,解决了船载、地面机动等机动布站望远镜的自主定向问题。
具体实施方式
发明内容部分的方案即是本发明的最佳实施方式。
本发明提供了一种可以全自动实现的望远镜自主定向方法,即使在机动望远镜天文经纬度未知的情况下也可实现。
在剔除非恒星目标环节,通过长短曝光图像结合的方式,剔除了图像中的卫星、撞星和其他假目标情况,提高了匹配效率和准确性;
在定标星选取环节,采用“回”形区域分区的方式确保了定标星之间的角距尽量大,以减少匹配冗余;
在恒星匹配环节,根据地面站的特点,将定标星的纬度排列作为匹配的输入条件之一,提高了匹配效率。通过角距与星等结合两轮匹配方案,解决了匹配冗余的难题。
在望远镜视轴中心指向解算环节,通过已知标准星的理想坐标和像面坐标建立望远镜测量模型,在望远镜站址精确天文坐标未知的情况下,依然可以解算出望远镜视轴中心指向。
Claims (1)
1.一种望远镜自主定向方法,其特征在于,包括以下步骤:
S1、获取天空图像:望远镜指向一个天区静止不动,连续获取长曝光天空背景图像和短曝光天空背景图像,分别统计两张图像上所有目标的像元数、像面坐标、灰度值,并记录曝光中心时刻;
S2、剔除非恒星目标:
统计长曝光图像中像元数大于阈值的目标的像面坐标,在短曝光图像中扣除以上目标;
对比长曝光图像和短曝光图像,若短曝光图像中有目标的位置,但长曝光图像中该位置没有目标,在短曝光图像中扣除该目标;
S3、将短曝光图像分割成“回”形区域,按照星像的平均灰度从大到小在四个“一”字型区域中的三个区域各选择一颗星作为定标星,分别记为N1、N2、N3,三个定标星俯仰角从高到低依次为ξ1、ξ2、ξ3,根据定标星的像面坐标计算出三颗星两两之间的角距,为σ12、σ23、σ13;
S4、检索恒星星库中对应赤纬从高到低分别为ξ1′、ξ2′、ξ3′的恒星作为标准星,构建出三个顶点分别位于纬度ξ1′、ξ2′、ξ3′上的所有三角形,计算出三角形的角距,形成配对用标准数据库,设用于配对的恒星编号分别为M1、M2、M3,三颗星的理论角距依次排列为γ12、γ23、γ13;所述恒星星库的构建方式为:将设定纬度范围内的恒星,按照纬度区域增加的方式进行分区,在各纬度区域,按星等增加的方式排列;
对定标星和标准星进行按以下条件进行匹配,确定出三颗定标星对应的标准星:
条件1:
min(σ12、σ23、σ13)≥min(γ12、γ23、γ13)+∈1
max(σ12、σ23、σ13)≤max(γ12、γ23、γ13)+∈1
条件2:
其中,∈1是设定误差门限1,∈2是设定误差门限2;
当符合要求的标准星超过一组时,根据星等进行再次匹配,将每组标准星按照星等从亮到暗排序,对应的赤纬分别为ξMmax′、ξMmid′、ξMmin′,对其两两做差,记录差值的正负性;对应的定标星按照星象平均灰度从高到低排序,对应的俯仰角分别为ξMmax、ξMmid、ξMmin,对其两两做差,记录差值的正负性,两组记录正负性的数组完全一致时,认为匹配成功,匹配成功后,通过查询定标星对应标准星在星库中的数据,得出定标星的赤经、赤纬;
S5、对于已知站点天文经纬度的测站,直接根据第i颗定标星的赤经、赤纬和地方时t求出定标星相对测站的理论方位角,再根据理论方位角和实测方位角之间的偏差ΔAi对方位编码器进行置数,实现定向;对于站点精确天文经纬度未知的测站,先根据站址的地理经纬度粗略估算出一个零位,再根据望远镜的误差模型进行修正,得出望远镜的精确零位,具体为:
根据测站粗略地理经纬度站址和第i颗定标星的赤经、赤纬(αi,δi)和地方时t粗略估算出测站中心点相对正北的偏差,为了提高精度,对多颗定标星的偏差值求平均值:
ΔAi=Ai-Ai'
式中:Ai'、Ai为第i颗定标星相对测站的理论方位角和实测方位角;
计算得出定标星的理想坐标(Xi,Yi):
式中:(αi,δi)为第i颗定标星的赤经、赤纬,(α0,δ0)为视轴中心指向的赤经、赤纬,(A0、E0)为视轴中心对应的方位俯仰值,(Ai、Ei)为第i颗定标星对应的方位俯仰值;
将定标星的理想坐标(Xi,Yi)代入六常数模型求出望远镜的测量模型,通过最小二乘法,求出望远镜中心偏差等误差:
式中:(xi,yi)为匹配成功的定标星的像面坐标,
修正望远镜中心指向,实现精确定向A0z:
A0z=A0+ΔAi+a
其中,A0为视轴中心对应的方位,a为六常数模型系数a。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210639586.0A CN114998423B (zh) | 2022-06-08 | 2022-06-08 | 一种望远镜自主定向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210639586.0A CN114998423B (zh) | 2022-06-08 | 2022-06-08 | 一种望远镜自主定向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114998423A CN114998423A (zh) | 2022-09-02 |
CN114998423B true CN114998423B (zh) | 2024-04-12 |
Family
ID=83033689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210639586.0A Active CN114998423B (zh) | 2022-06-08 | 2022-06-08 | 一种望远镜自主定向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114998423B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8767072B1 (en) * | 2010-03-26 | 2014-07-01 | Lockheed Martin Corporation | Geoposition determination by starlight refraction measurement |
CN111156989A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位方法 |
CN111486868A (zh) * | 2020-03-30 | 2020-08-04 | 中国人民解放军63636部队 | 基于地物特征的光电望远镜免方位标扩展标校法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110352420A (zh) * | 2018-09-13 | 2019-10-18 | 陈加志 | 一种基于图像识别的望远镜搜星方法、搜星装置及望远镜 |
-
2022
- 2022-06-08 CN CN202210639586.0A patent/CN114998423B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8767072B1 (en) * | 2010-03-26 | 2014-07-01 | Lockheed Martin Corporation | Geoposition determination by starlight refraction measurement |
CN111156989A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位方法 |
CN111486868A (zh) * | 2020-03-30 | 2020-08-04 | 中国人民解放军63636部队 | 基于地物特征的光电望远镜免方位标扩展标校法 |
Also Published As
Publication number | Publication date |
---|---|
CN114998423A (zh) | 2022-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106871927B (zh) | 一种无人机光电吊舱安装误差标校方法 | |
Clarke et al. | The development of camera calibration methods and models | |
CN104913780B (zh) | 集成gnss和ccd天顶筒的高精度垂线偏差快速测量方法 | |
CN109443359B (zh) | 一种地面全景影像的地理定位方法 | |
CN106500731A (zh) | 一种基于恒星模拟系统的船载经纬仪的标校方法 | |
CN111811538B (zh) | 一种基于常见星体的光电经纬仪定向标校方法 | |
CN111486868B (zh) | 基于地物特征的光电望远镜免方位标扩展标校法 | |
CN113538595B (zh) | 利用激光测高数据辅助提升遥感立体影像几何精度的方法 | |
CN111156989B (zh) | 基于指向自动测定的空间碎片实时天文定位方法 | |
CN114001756B (zh) | 一种小视场星敏感器外场地面寻星方法 | |
CN107589431B (zh) | 一种提高机载光电系统目标定位精度的校靶方法 | |
CN115511956A (zh) | 一种无人机成像定位方法 | |
CN110068313B (zh) | 一种基于投影变换的数字天顶仪定向方法 | |
CN114167900B (zh) | 一种基于无人机和差分gps的光电跟踪系统标校方法及装置 | |
CN115597586A (zh) | 一种基于大气偏振模式对称性的航向角提取方法 | |
CN114998423B (zh) | 一种望远镜自主定向方法 | |
WO2021135161A1 (zh) | 基于指向自动测定的空间碎片实时天文定位和测光方法 | |
CN110068312B (zh) | 一种基于球面三角形的数字天顶仪定位方法 | |
CN115218861B (zh) | 一种基于自动经纬仪的天文方位测量方法 | |
CN116309798A (zh) | 一种无人机成像定位方法 | |
CN116124094A (zh) | 基于无人机侦察图像与组合导航信息的多目标协同定位方法 | |
CN113109829B (zh) | 一种同步扫描交会测量敏感器的标定方法 | |
CN104458653A (zh) | 一种测量大天顶距处的大气折射值的方法和系统 | |
CN115326007A (zh) | 一种基于恒星观测的船载经纬仪测角时序检测方法 | |
CN111879298B (zh) | 一种用于月球无人值守望远镜的自动指星方法 |
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 |