CN110068313B - 一种基于投影变换的数字天顶仪定向方法 - Google Patents

一种基于投影变换的数字天顶仪定向方法 Download PDF

Info

Publication number
CN110068313B
CN110068313B CN201910381422.0A CN201910381422A CN110068313B CN 110068313 B CN110068313 B CN 110068313B CN 201910381422 A CN201910381422 A CN 201910381422A CN 110068313 B CN110068313 B CN 110068313B
Authority
CN
China
Prior art keywords
star
zenith instrument
value
azimuth
orientation
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
Application number
CN201910381422.0A
Other languages
English (en)
Other versions
CN110068313A (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.)
Rocket Force University of Engineering of PLA
Original Assignee
Rocket Force University of Engineering 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 Rocket Force University of Engineering of PLA filed Critical Rocket Force University of Engineering of PLA
Priority to CN201910381422.0A priority Critical patent/CN110068313B/zh
Publication of CN110068313A publication Critical patent/CN110068313A/zh
Application granted granted Critical
Publication of CN110068313B publication Critical patent/CN110068313B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • 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

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)
  • Navigation (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及一种基于投影变换的数字天顶仪定向方法，采用具有星图与识别、定位解算软件***的通用数字天顶仪，特征在于包括四大步骤：测量天顶仪所在位置的天文坐标；构建恒星坐标的投影模型；方位角的计算；方位角A的修正。同现有技术相比，受光学***参数影响小，即使没有对光学***进行标校也能够准确解算方位，省去了非常困难的标校焦距、畸变等光学***参数的过程，也不会因为会受到外部因素的影响而随着时间发生相应的变化，且只需通过天文观测即可确定天顶仪方位并输出方位能够实现测站点的快速定向，环境适应性较强。实验验证，本发明能够高精度完成的定向精度在8＂以内，完全满足实际工作对于定向精度的要求。

Description

一种基于投影变换的数字天顶仪定向方法
技术领域
本发明属于大地天文测量技术领域，涉及一种天顶仪的定向方法，尤其是涉及一种基于投影变换的数字天顶仪定向方法。
背景技术
国内外研究数字天顶仪技术的主要有三家，德国Hannover大学、瑞士苏黎世理工大学、二炮装备研究院。三家所研究的数字天顶仪软硬件***基本相似，硬件结构主要由光学***、倾角测量与调平***、授时***、控制***与数据处理***构成。软件***主要包括星图与识别、天顶仪旋转轴解算和相关补偿等。在本发明以前的现有技术中，针对数字天顶仪的研究主要集中在通过硬件结构和精度改变来提高定位的精度上，其中基于星敏感器的多矢量方法，可以通过观测恒星实现定向，它的缺点在于，多矢量方法在定向时需要标校焦距、畸变等光学***参数，而这些***参数的标校往往是非常困难的，且参数会受到外部因素的影响而随着时间发生相应的变化，一般的定位精度在10＂以上，因此这种技术并不适合在数字天顶仪中使用。《数字天顶仪旋转轴解算方法》(发明专利申请号：201110406263.9)中采用了不同方位拍摄的星图，对数字天顶仪的旋转轴坐标提出了针对性的解算方法，但该方法并未就定位方法提供更具体的技术解决方案。至目前为止，数字天顶仪定向方法在国内外还未有专门的研究报道。
发明内容
针对上述现有技术状况，本发明的目的在于，提供一种基于投影变换的思路、在不对光学***进行标校的情况下准确解算方位、且通过天文观测确定天顶仪方位并输出方位的数字天顶仪定向方法。
现将本发明构思及技术解决方案叙述如下：
本发明的基本构思是，所提供的定向方法，应该是受光学***参数的影响很小，即使没有对光学***进行标校也能够准确解算方位，省去了非常困难的标校焦距、畸变等光学***参数的过程，同时也不会因为会受到外部因素的影响而随着时间发生相应的变化，且只需用本发明方法通过天文观测即可确定天顶仪方位并输出方位的的基于投影变换思路的数字天顶仪定向方法。
本发明一种基于投影变换的数字天顶仪定向方法，硬件***采用通用数字天顶仪，包括光学望远镜、CCD成像装置及旋转平台调平装置；软件***包括星图与识别、天顶仪定位解算和相关补偿，其特征在于：包括以下步骤：
步骤1：测量天顶仪所在位置的天文坐标
在进行天文定向前，先测量数字天顶仪所在位置的天文坐标：(α0,δ0)，在定位基础上对数字天顶仪进行定向，天顶仪的定位精度在0.3＂以内；
步骤2：构建恒星坐标的投影模型
步骤2.1：拍摄的星图做预处理，去除噪声；
步骤2.2：提取星图中恒星中心的像片坐标(xij,yij)，并结合恒星星表识别所拍摄到的恒星；
步骤2.3：根据拍摄时刻信息及星表中的相关参量，计算所识别恒星于拍摄历元在地球固联坐标系中的位置(αij,δij)；
步骤2.4：根据已经测量得到的天顶仪高精度天文坐标(α0,δ0)，计算恒星在切平面坐标系中的坐标：
Figure BDA0002053490120000031
步骤2.5：建立每幅图像的投影模型，利用最小二乘法解算参数
[ai1,bi1,ci1,ai2,bi2,ci2,ai3,bi3]
Figure BDA0002053490120000032
步骤3：方位角的计算
步骤3.1：利用计算得到的参数，分别计算每幅图像的方位角：
Figure BDA0002053490120000033
方位角αi的取值范围根据其正弦值与余弦值的正负情况确定：
Figure BDA0002053490120000034
式(1)到(4)式中，i表示图像编号，j表示恒星编号，i＝j＝1,2...；
步骤3.2：将每幅图像所计算得到的多个方位角去除异常值之后，将其均值作为最终的方位角α；
步骤3.3：通过公式(4)可以得到多个方位角的输出值，对输出值进行分析，通过对残差的处理筛选出含有粗大误差的星点，可得残差的均值为
Figure BDA0002053490120000035
残差的标准差
Figure BDA0002053490120000041
将方位角输出值按照大小顺序排列α(1)≤α(2)≤…≤α(n)；
步骤3.4：构建统计量
Figure BDA0002053490120000042
当g(1)≥g0(n,a)，认为测得值α(1)是异常数据，应剔除；当g(n)≥gn(n,a)，认为α(n)是异常数据，应剔除。对异常值踢出后求解均值，，将其均值作为最终的方位角α。
步骤4：.方位角A的修正
步骤4.1：经度修正公式：
Δα1＝Δαsinδ0 (6)
式(6)中，Δα＝α0-α1；δ0为定位后所得天文坐标的纬度值；
步骤4.2：纬度的修正
因纬度所产生的影响较小，大多数情况下可以忽略，但当定位所得天文坐标纬度值与所使用概略位置纬度值超过1°时，应用下式进行修正：
Δα2≈1/2ac(Δδ)2 (6)
步骤4.3：方位的输出
最终的输出方位角为：
α＝α+Δα1+Δα2 (7)
本发明通现有技术相比，能够实现测站点的快速定向，而且环境适应性较强。经实际实验验证，利用本发明能够高精度的完成定向，定向的精度在8＂以内，完全可以满足实际工作中对于定向精度的要求。
附图说明
图1：数字天顶仪的硬件***示意图
图2：数字天顶仪方位输出示意图
具体实施方式
现结合附图对本发明方法作进一步说明。
参见图1：本发明所用数字天顶仪主要由光学望远镜、CCD成像装置及旋转平台调平装置组成；本发明方法具体包括以下步骤：
步骤1：:测量天顶仪所在位置的天文坐标
在进行天文定向前，先测量天顶仪所在位置的天文坐标，得到较高精度的天文坐标：(α0,δ0)；
步骤2：.建立同一位置、不同图像的投影模型
步骤2.1：对拍摄的星图做预处理，去除噪声；
步骤2.2：提取星图中恒星中心的像片坐标(xij,yij)，并结合恒星星表识别所拍摄到的恒星；
步骤2.3：根据拍摄时刻、IAU2003模型及星表中的相关参量，计算所识别恒星于拍摄历元在地球固联坐标系中的位置(αij,δij)；
步骤2.4：根据已经测量得到的天顶仪高精度天文坐标(α0,δ0)，计算恒星在切平面坐标系中的坐标。
Figure BDA0002053490120000051
步骤2.5：建立每幅图像的投影模型，利用最小二乘法解算参数向量
[ai1,bi1,ci1,ai2,bi2,ci2,ai3,bi3]:
Figure BDA0002053490120000061
步骤3：方位角的计算
步骤3.1：利用计算得到的参数，分别计算每幅图像的方位角：
Figure BDA0002053490120000062
步骤3.2：方位角αi的取值范围根据其正弦值与余弦值的正负情况确定：
Figure BDA0002053490120000063
式(1)到(4)式中，i表示图像编号，j表示恒星编号，i＝j＝1,2...；
步骤3.3：通过公式(4)可以得到多个方位角的输出值，对输出值进行分析，通过对残差的处理筛选出含有粗大误差的星点，可得残差的均值为
Figure BDA0002053490120000064
残差的标准差
Figure BDA0002053490120000065
将方位角输出值按照大小顺序排列α(1)≤α(2)≤…≤α(n)；
步骤3.4：构建统计量
Figure BDA0002053490120000066
当g(1)≥g0(n,a)，认为测得值α(1)是异常数据，应剔除；当g(n)≥gn(n,a)，认为α(n)是异常数据，应剔除，对异常值踢出后求解均值，将其均值作为最终的方位角α；
步骤4：.方位角A的修正
步骤4.1：经度修正公式：
Δα1＝Δαsinδ0 (6)
式(6)中，Δα＝α0-α1；δ0为定位后所得天文坐标的纬度值。
步骤4.2：纬度的修正
因纬度所产生的影响较小，大多数情况下可以忽略。但是，当定位所得天文坐标纬度值与所使用概略位置纬度值相差比较大时，应用下式进行修正：
Δα2≈1/2ac(Δδ)2 (7)
步骤4.3：方位的输出
最终的输出方位角为：
α＝α+Δα1+Δα2 (8)
经实验验证，利用本发明能够实现快速的天文定向，定向的精度在8＂以内。由于该方法通过数字天顶仪在单一位置拍摄的数幅恒星星图进行定向，并采用了投影变换模型，因此该方法对外界的适应性较强。

Claims (1)

1.一种基于投影变换的数字天顶仪定向方法，硬件***采用通用数字天顶仪，包括光学望远镜、CCD成像装置及旋转平台调平装置；软件***包括星图与识别、天顶仪定位解算和相关补偿，其特征在于：包括以下步骤：
步骤1：测量天顶仪所在位置的天文坐标；在进行天文定向前，先测量数字天顶仪所在位置的高精度天文坐标(α0,δ0)，在定位基础上对数字天顶仪进行定向，天顶仪的定位精度在0.3＂以内；
步骤2：构建恒星坐标的投影模型；
步骤2.1：拍摄的星图做预处理，去除噪声；
步骤2.2：提取星图中恒星中心的像片坐标(xij,yij)，并结合恒星星表识别所拍摄到的恒星；
步骤2.3：根据拍摄时刻及星表中的相关参量，计算所识别恒星于拍摄历元在地球固联坐标系中的位置(αij,δij)；
步骤2.4：根据已经测量得到的天顶仪高精度天文坐标(α0,δ0)，计算恒星在切平面坐标系中的坐标：
Figure FDA0003071500690000011
步骤2.5：建立每幅图像的投影模型，利用最小二乘法解算参数[ai1,bi1,ci1,ai2,bi2,ci2,ai3,bi3]
Figure FDA0003071500690000021
式(1)、(2)式中，i表示图像编号，j表示恒星编号，i＝j＝1,2...；
步骤3：方位角α的计算；
步骤3.1：利用计算得到的参数，分别计算每幅图像的方位角：
Figure FDA0003071500690000022
方位角αi的取值范围根据其正弦值与余弦值的正负情况确定：
Figure FDA0003071500690000023
式(3)、(4)式中，i表示图像编号，j表示恒星编号，i＝j＝1,2...；
步骤3.2：将每幅图像所计算得到的多个方位角αi去除异常值之后，将其均值作为最终的方位角α′；
步骤3.3：通过公式(4)可以得到多个方位角αi的输出值，对输出值进行分析，通过对残差的处理筛选出含有粗大误差的星点，可得方位角的均值为：
Figure FDA0003071500690000024
残差的标准差
Figure FDA0003071500690000025
将方位角αi输出值按照大小顺序排列α1≤α2≤…≤αn；
步骤3.4：构建统计量
Figure FDA0003071500690000026
Figure FDA0003071500690000027
为方位角的均值；当g1≥g0(n,a)，认为测得值α1是异常数据，应剔除；当gn≥gn(n,a)，认为αn是异常数据，应剔除；对异常值剔除后求解均值，将其均值作为最终的方位角α′；
步骤4：方位角α的修正
步骤4.1：经度修正公式：
Δα1＝Δαsinδ0 (6)
式(6)中，Δα为经度偏差量；δ0为定位后所得天文坐标的纬度值；
步骤4.2：纬度的修正
因纬度所产生的影响较小，大多数情况下可以忽略，但当定位所得天文坐标纬度值与所使用概略位置纬度值大于1°时，应用下式进行修正：
Δα2≈1/2ac(Δδ)2 (7)
步骤5：方位的输出
最终的输出方位角α为：
α＝α′+Δα1+Δα2 (8)。
CN201910381422.0A 2019-05-08 2019-05-08 一种基于投影变换的数字天顶仪定向方法 Active CN110068313B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910381422.0A CN110068313B (zh) 2019-05-08 2019-05-08 一种基于投影变换的数字天顶仪定向方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910381422.0A CN110068313B (zh) 2019-05-08 2019-05-08 一种基于投影变换的数字天顶仪定向方法

Publications (2)

Publication Number Publication Date
CN110068313A CN110068313A (zh) 2019-07-30
CN110068313B true CN110068313B (zh) 2021-08-20

Family

ID=67370333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910381422.0A Active CN110068313B (zh) 2019-05-08 2019-05-08 一种基于投影变换的数字天顶仪定向方法

Country Status (1)

Country Link
CN (1) CN110068313B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112284333B (zh) * 2020-10-23 2022-06-17 中国科学院云南天文台 一种多功能天文经纬仪ccd照相机靶面偏转改正方法
CN114111762B (zh) * 2021-11-15 2023-06-06 北京航天计量测试技术研究所 一种基于双轴水平仪的单星定向方法
CN115406408B (zh) * 2022-10-31 2023-01-17 中国科学院长春光学精密机械与物理研究所 光电经纬仪垂直轴倾斜误差检测、修正方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521506A (zh) * 2011-12-09 2012-06-27 中国人民解放军第二炮兵装备研究院第五研究所 数字天顶仪旋转轴解算方法
CN103047985A (zh) * 2012-12-21 2013-04-17 中国科学院长春光学精密机械与物理研究所 一种空间目标的快速定位方法
CN108362268A (zh) * 2018-02-06 2018-08-03 中国人民解放军战略支援部队信息工程大学 一种基于视频测量的自动天文测量方法和测量***

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3778000B2 (ja) * 2001-04-11 2006-05-24 三菱電機株式会社 写真天頂望遠鏡
CN103017762B (zh) * 2012-12-21 2015-10-28 中国科学院长春光学精密机械与物理研究所 地基光电望远镜的空间目标快速捕获定位方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521506A (zh) * 2011-12-09 2012-06-27 中国人民解放军第二炮兵装备研究院第五研究所 数字天顶仪旋转轴解算方法
CN103047985A (zh) * 2012-12-21 2013-04-17 中国科学院长春光学精密机械与物理研究所 一种空间目标的快速定位方法
CN108362268A (zh) * 2018-02-06 2018-08-03 中国人民解放军战略支援部队信息工程大学 一种基于视频测量的自动天文测量方法和测量***

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
改进的数字天顶仪定位方法;张新帅等;《红外与激光工程》;20150430;第44卷(第4期);第1254-1259页 *

Also Published As

Publication number Publication date
CN110068313A (zh) 2019-07-30

Similar Documents

Publication Publication Date Title
CN108921901B (zh) 一种基于精密二轴转台和激光跟踪仪的大视场相机标定方法
Clarke et al. The development of camera calibration methods and models
CN110068313B (zh) 一种基于投影变换的数字天顶仪定向方法
CN109655079B (zh) 星敏感器测量坐标系到棱镜坐标系测量方法
CN112629431B (zh) 土木结构变形监测方法及相关设备
CN113538595B (zh) 利用激光测高数据辅助提升遥感立体影像几何精度的方法
CN100376883C (zh) 一种基于像素频率的星敏感器校准方法
CN106468544B (zh) 基于光电自准直仪的卫星高精度测角方法
CN109712201B (zh) 一种广域相机的定位能力标定装置及标定方法
CN116597013B (zh) 一种基于不同经纬度区域的卫星图像几何定标方法
CN110363758B (zh) 一种光学遥感卫星成像质量确定方法及***
CN111473802A (zh) 一种基于线阵推扫的光学传感器内方位元素定标方法
CN110501026B (zh) 基于阵列星点的相机内方位元素标定装置及方法
CN101655378A (zh) 一种数字摄影测量相机的内方位参数测量方法
CN103940590A (zh) 大口径光学镜头畸变的标定方法
CN108154535B (zh) 基于平行光管的摄像机标定方法
CN108447100B (zh) 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法
CN113218577A (zh) 一种星敏感器星点质心位置精度的外场测量方法
CN105387996B (zh) 多光轴地面恒星观测***光轴一致性检校方法
CN115222819A (zh) 一种基于机场大范围场景下多模信息参照的相机自标定和目标追踪方法
CN106289317A (zh) 一种单镜头数字航测相机的单机检校方法及装置
CN116123998A (zh) 多站点基于视频采集对空中炸点进行实时测量方法
CN114001756B (zh) 一种小视场星敏感器外场地面寻星方法
CN112802115B (zh) 一种多焦面拼接大视场离轴相机的几何标定方法及装置
CN101726316A (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