CN111678536A - 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法 - Google Patents

地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法 Download PDF

Info

Publication number
CN111678536A
CN111678536A CN202010422093.2A CN202010422093A CN111678536A CN 111678536 A CN111678536 A CN 111678536A CN 202010422093 A CN202010422093 A CN 202010422093A CN 111678536 A CN111678536 A CN 111678536A
Authority
CN
China
Prior art keywords
point
observation
angle
error
azimuth
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
CN202010422093.2A
Other languages
English (en)
Other versions
CN111678536B (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.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University 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 Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN202010422093.2A priority Critical patent/CN111678536B/zh
Publication of CN111678536A publication Critical patent/CN111678536A/zh
Application granted granted Critical
Publication of CN111678536B publication Critical patent/CN111678536B/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
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
    • 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

Abstract

本发明提供一种磁偏角及观瞄设备测角系统误差一并校准的标定实施方法,包括下列步骤:选择观测点及标定点并分别录取坐标;测量标定点的磁方位角和高低角;解算标定点的真方位角和高低角;消除观瞄设备系统误差及观测点磁偏角。此外,作为理论基础,还提供一种磁偏角及观瞄设备测角误差一并校准方法。同时,又给出一种基于定位精度确定标定点距离的基本方法。本发明提出基于标定点和观察者站立点地理坐标校准测量设备测角系统误差和本地磁偏角的模型,给出模型参数性能指标和标定实施步骤,能够解决观瞄设备真北方位角输出、系统误差校正等问题。本发明实施成本低、操作简单、效果好,具有很大的实用价值。

Description

地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定 方法
技术领域
本发明涉及空间目标定位、观测设备基准和测角系统误差校准领域,具体涉及地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法。
背景技术
近年来,随着低空空域的逐步开放和国内外低空有人、无人飞行器井喷式发展,低空目标探测需求激增。由于对低空目标探测的范围有限,需要大量配置对空观测设备进行探测。广泛分布的众多低空探测设备,需要在技术上确保时空一致性和误差的可控性,才能满足预警监视信息的有效、可用。
当前基于双天线的高精度GPS定位测向系统通过静态实验和动态实验分析系统定位、测角精度,虽然能达到经度的标准偏差在0.315m、纬度的标准偏差在0.196m,测向精度偏差在0.09°的精度级别。但观瞄设备受其尺寸、重量等便携因素影响,在实际野外观测校准活动中适用度不高,在实际观测活动中一般利用地磁指向作为寻北基准。但地磁受地理环境、地质矿藏、日月变化以及地壳运动等不可抗拒的因素影响,各地的地磁场强分布不均匀,且处于随时变化中。因此,利用地磁作为观测设备对目标实施测角时,需要将各点测量的磁方位角转化为真方位角,才能保证多用户方向基准的一致性。而某一地点、某一时刻的磁偏角,需要查询磁偏角及其年变率等数据库才能获取,给用户使用带来诸多不便。为此,迫切需要一种消除磁偏角的简易方法。
同时,任何观瞄设备都存在测高和测角的系统误差,特别是便携式简易观瞄设备,每次展开都会产生不同的系统误差,而其难以获知具体误差数据,在实际操作中系统误差的检测、消除的方法通常十分复杂繁琐,对检测设备的要求较高,成本和时效性代价昂贵。为此,迫切需要一种消除测角系统误差的建议方法。
发明内容
针对现有技术存在的问题,本发明提供一种磁偏角及观瞄设备测角误差一并校准方法,具体包括下列步骤:
磁偏角θ直接影响各个观测点的观瞄设备获取观测值的空间一致性;同时,观瞄设备的测角系统误差也需要消除;依据以下方法将二者一并消除;
第一,建立地面直角坐标系;以观测点A为坐标原点,其地理坐标为
Figure BSA0000209070230000021
Figure BSA0000209070230000022
λa分别表示A点的纬度和精度,以A点所在纬线的水平投影指东为Z轴,以所在点经线的水平投影指北为X轴,X轴为水平面向上的法线,即高程轴;AD为磁北方向线;
第二,在观测点A一定距离处选择标定点
Figure BSA0000209070230000023
利用观瞄设备或手持GPS/北斗终端、手机这样的其他定位手段,分别测量观测点A和标定点P的经纬坐标和高程,B为标定点P在观测点所在水平面XAZ上的投影;
第三,从B点做X轴的垂线相交于C点,依据测量的A和标P的经纬度坐标和高程,利用半正矢haversine算法求解出A、B、C三点两两之间的球面距离dab、dac和dbc;进而利用三角函数关系,计算出标定点P相对于观测点A的真方位角βz=∠CAB,即从经线北顺时针到观察线水平投影的夹角,和高低角εz=∠PAB;
第四,利用观瞄设备从观测点直接测量标定点P包含方位角系统误差的磁方位角βc=∠DAB,即从磁北顺时针到观察线水平投影的夹角,和包含高低角系统误差的高低角εC=∠PAB;
第五,观察标定点P所实际测量得的磁方位角βc减去计算所得的真方位角βz,就得到观瞄设备所在位置磁偏角和观瞄设备方位角系统误差之和θ;所测高低角εc减去计算高低角εZ得到高低角系统误差ψ;
以解算真方位角βz替代实际的真方位角,磁偏角为δ,则磁方位角β=βz-δ;测量设备方位角系统误差为γ,测量所测磁方位角为βc,则磁方位角为β=βc
即:βz-δ=βc-γ,
所以,实测磁方位角减去解算真方位角即为磁偏角和方位角系统误差之和:γ-δ=βcz
即:
θ=γ-δ=βcz (1)
以解算真高低角εz替代实际的真高低角,测量所得P点的高低角为εc,则
ψ=εcz (2)
第六,将包含磁偏角和方位角系统误差的方位角误差θ和高低角系统误差ψ输入观瞄设备相应置数区,用于在输出测量数据时消除观瞄设备的磁偏角和测角误差,使装置直接输出真方位角和高低角,且数值相对精确;
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc
还提供一种基于定位精度确定标定点距离的基本方法,具体包括下列步骤:
观测点A与标定点P之间的球面距离为dab,则由水平定位误差Δd引起的真方位角解算误差最大值为:
Δβ=arcsin(2Δd/dab) (3)
同理,当A与P两点在同一水平面,即标定点与观测点高程差为零(H=0)时,由垂直定位误差Δh引起的高低角解算误差最大值为:
Δε=arctan(Δh/dab) (4)
观测点A与标定点P之间的水平距离为dab、高程差为H,则AP之间距离为的
Figure BSA0000209070230000041
由定位高程误差引起的高低角解算误差的最大值近似为:
Figure BSA0000209070230000042
在本发明的一个实施例中,选择标定点距离较近时,尽量选择与观测点高差较大的目标点,有利于消除高低角解算误差。
此外,还提供一种磁偏角及观瞄设备测角系统误差一并校准的标定实施方法,利用所述方法消除观瞄设备系统误差及观测点磁偏角,具体有三种方式:一是观测设备具有数据处理和存储功能时,将包含磁偏角和方位角系统误差的方位角误差θ和观测设备高低角系统误差ψ解算出来装订存储在观瞄设备中,将每组测量的数据进行校准后发送;二是观测设备具有数据处理功能,不具有储存功能时,将校准数据θ和ψ发送至后台存储,由后台对各个观察哨实测高低角和磁方位角进行校准;三是观测设备既不具备解算功能,也不具备存储功能时,在后台安装所述方法相应的算法,在使用前,各观测点将包括经度、纬度、高程在内的观测点和标定点坐标、实测标定点高低角和方位角发送至后台,后台根据公式(15)、(16)解算出各观测点的校准数据θ和ψ并储存;具体如下:
(一)选择观测点及标定点并分别录取坐标
到达现地后,根据预案和现地实际情况选择标定点P和观测点A,首先,按照标志明显、便于观察和瞄准的要求选择标定点,测量、录入标定点经纬度坐标
Figure BSA0000209070230000043
和高程hp,其中
Figure BSA0000209070230000044
λp分别表示P点的纬度和经度,并做好标记,以便观察瞄准;其次,依据便于观察瞄准标定点、便于观察发现目标、便于隐蔽自身、便于通信联络、交通便利等因素确定观测点,测量、录入观测点经纬度坐标
Figure BSA0000209070230000051
和高程ha,其中
Figure BSA0000209070230000052
λa分别表示A点的纬度和经度;
(二)测量标定点的磁方位角和高低角
在确定标定点P和观测点A并完成定位数据录入之后,使用观测设备的方位和高低测量功能直接测出标定点的磁方位角和高低角的测量数据,并录入系统;其中,磁方位角βc包含方位角系统误差和磁偏角,高低角εc包含高低角系统误差;在测量标定点的方位角和高低角时,如果时间允许,采用经过多次测量取均值的方法,能够得到比较精准的磁方位角和高低角;
实施标定校准观测时,应保持观测设备稳定、标定点选取距离适中,数据测定准确,为后续精准测算方位角误差θ、高低角系统误差ψ提供数据支撑,其中方位角误差θ包含磁偏角系统误差和方位角系统误差;
(三)解算标定点的真方位角和高低角
根据第一步测定的观测点经纬度坐标
Figure BSA0000209070230000053
和高程ha,以及标定点经纬度坐标
Figure BSA0000209070230000054
和高程hp,依据公式(6-14),解算标定点的真方位角βz和高低角εZ;为介绍时方便计算,以下将椭圆形地球简化成球形地球,其半径近似为:
Figure BSA0000209070230000055
1.解算标定点P点真方位角
准确测得点A和点P的经纬度坐标后,要测算标定点P的真方位角βz=∠BAC,需要先解算出dab和dbc的距离,dab表示A点和B点之间的球面距离,dbc表示B点和C点之间的球面距离;直接使用求解球面两点间距离的半正矢haversine算法,解算出dbc如下:
Figure BSA0000209070230000061
Figure BSA0000209070230000062
式中,
Figure BSA0000209070230000063
点的地理坐标为
Figure BSA0000209070230000064
表示为B、C两点的纬度差,Δλbc=λbc表示B、C两点的经度差,abc和cbc为中间变量;
则,B、C两点间的球面距离为:
dbc=R·cbc (9)
其中,R为地球半径;
同理,可以直接解算得到A、B两点间的球面距离dab如下:
Figure BSA0000209070230000065
Figure BSA0000209070230000066
dab=R·cab (12)
式中,
Figure BSA0000209070230000067
表示为B、A两点的维度差,Δλab=λba表示B、A两点的经度差,aab和cab为中间变量,dab为A、B两点间的球面距离;
根据上式(9)和(12)解算得到的BC、AB的距离,根据三角函数之反正弦定理解算得到标定点P的真方位角为:
Figure BSA0000209070230000068
2.解算标定点P点高低角
由上公式(11)解算得到A、B两点间的距离为dab,根据三角函数之反正切函数定理,解算得标定点P的真实高低角εz如下:
Figure BSA0000209070230000069
式中,hpb为P点和B的高程差,在数值上等于P点和A点的高程差,能够通过GPS/北斗设备获得;
(四)消除观瞄设备系统误差及观测点磁偏角
在定位精度和标定点距离满足要求的情况下,认为依据定位数据解算的方位角βz和高低角εz是准确的,能够以此为基准消除测量设备的系统误差,且由于解算的方位角是真方位角,测量的方位角是磁方位角,消除方位角误差的同时也就完成了测量磁方位角转换为真方位角的工作;
根据公式(1),
θ=γ-δ=βcz (1)
其中,θ为磁偏角和方位角系统误差之和,γ为测量设备方位角系统误差,δ为磁偏角,βc为测量所得磁方位角,βz为解算的真方位角;
包含磁偏角和方位角系统误差的方位角误差θ为:
θ=βcz (15)
根据公式(2),高低角系统误差为ψ:
ψ=εcz (16)
其中,ψ为高低角系统误差,εc为测量所得P点的高低角,εz为解算真高低角;
根据观测点上述校准数据θ和ψ,对测量的每一组数据进行校准;
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc
本发明基于观瞄设备自身的GPS/北斗终端装置的定位功能,通过标定,解算出磁偏角与方位角系统误差之和,以及高低角系统误差,从而实现观测点磁偏角与观察设备测角系统误差一并校准和消除。
附图说明
图1示出定位点、标定点位置关系示意图;
图2示出真方位角、磁方位角、实测磁方位角关系示意图;
图3示出实测高低角与解算高低角偏差截面示意图;
图4示出解算真方位角误差受水平定位精度、标定点距离影响示意图;
图5示出解算高低角误差受高程定位精度、标定点距离、标定点高差影响示意图;
图6示出标定实施基本过程。
具体实施方法
一、磁偏角及观瞄设备测角误差一并校准基本原理
观瞄设备通常都具有定位和磁北测量功能。但是,地球上任意一点的磁北方向线与经线北方向线在在水平面上存在数字不等的夹角,称为磁偏角θ,这就直接影响了各个观测点获取观测值的空间一致性。同时,观瞄设备的测角系统误差也需要消除。可依据以下原理将二者一并消除。
第一,建立地面直角坐标系如图1所示。以观测点A为坐标原点,其地理坐标为
Figure BSA0000209070230000081
λa分别表示A点的纬度和精度,以A点所在纬线的水平投影指东为Z轴,以所在点经线的水平投影指北为X轴,X轴为水平面向上的法线,即高程轴;AD为磁北方向线。
第二,在观测点A一定距离处选择标定点
Figure BSA0000209070230000082
利用便携观瞄设备或手持GPS/北斗终端、手机等其他定位手段,分别测量观测点A和标定点P的经纬坐标和高程,B为标定点P在观测点所在水平面XAZ上的投影。
第三,从B点做X轴的垂线相交于C点,依据测量的A和标P的经纬度坐标和高程,利用半正矢(haversine)算法(wikipedia.Haversine formula[EB/OL].[2017-06-01].https://en.wikipedia.org/wiki/Haversine_formula.)求解出A、B、C三点两两之间的球面距离dab、dac和dbc;进而利用如公式(13)所示的三角函数关系,计算出标定点P相对于观测点A的真方位角βz=∠CAB(经线北顺时针到观察线水平投影的夹角)和高低角εz=∠PAB。具体方法见“(三)解算标定点的真方位角和高低角”。
第四,利用观瞄设备从观测点直接测量标定点P包含方位角系统误差的磁方位角βc=∠DAB(磁北顺时针到观察线水平投影的夹角)和包含高低角系统误差的高低角εC=∠PAB。
第五,观察标定点P所实际测量得的磁方位角βc减去计算所得的真方位角βz,就得到了观瞄设备所在位置磁偏角和观瞄设备方位角系统误差之和θ;所测高低角εc减去计算高低角εZ就得到了高低角系统误差ψ。
如图2所示,以解算真方位角βz替代实际的真方位角,磁偏角为δ,则磁方位角β=βz-δ;测量设备方位角系统误差为γ,测量所测磁方位角为βc,则磁方位角为β=βc
即:βz-δ=βc-γ,
所以,实测磁方位角减去解算真方位角即为磁偏角和方位角系统误差之和:γ-δ=βcz
即:
θ=γ-δ=βcz (1)
如图3所示,以解算真高低角εz替代实际的真高低角,测量所得P点的高低角为εc,则
ψ=εcz (2)
第六,将包含磁偏角和方位角系统误差的方位角误差θ和高低角系统误差ψ输入观瞄设备相应置数区,用于在输出测量数据时消除观瞄设备的磁偏角和测角误差,使装置直接输出真方位角和高低角,且数值相对精确。
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc
二、基于定位精度确定标定点距离的基本方法
从上述原理可以看出,这种标定方法校准后的剩余误差主要与算法、定位精度、操作误差和标定点与观测点距离有关;通常,在标定设备定位误差、算法误差和操作人员测量误差都应在用户能够接受的范围之内,且相对稳定。这样,标定的精度主要由标定点与观测点的距离决定,为此,提出了基于特定定位精度的标定点确定方法,为该方法的使用提供依据。
定位精度由定位设备的指标确定,在用户使用过程中通常是固定值,很显然,在忽略距离解算算法误差、人员测量误差的情况下,标定点选择的越远,标定精度越高。但是,标定点的选择过远会受到时间和各种保障条件制约,过近则无法满足精度要求。为此,标定点的合理选择十分重要。
观测点A与标定点P之间的球面距离为dab,如图4所示,则由水平定位误差Δd引起的真方位角解算误差最大值为:
Δβ=arcsin(2Δd/dab) (3)
同理,当A与P两点在同一水平面,即标定点与观测点高程差为零(H=0)时,由垂直定位误差Δh引起的高低角解算误差最大值为:
Δε=arctan(Δh/dab) (4)
例如,某地面观察哨观瞄设备的GPS/BD水平和垂直定位精度为Δd=Δh≤10m(陈杰生等,地面防空“观察哨”管理控制信息系统构建问题研究[J],空军工程大学学报,2019.03);则由水平定位误差引起的真方位角解算误差最大值,以及观测点与标定点高程差为零(H=0)时垂直定位误差引起的高低角解算误差最大值,随标定点距离变化情况见表1。
表1定位精度Δd=Δh=10m时标定点距离对解算测角误差的影响
Figure BSA0000209070230000111
观测点A与标定点P之间的水平距离为dab、高程差为H,如图5所示,则AP之间距离为的
Figure BSA0000209070230000112
与前述原理相似,则由定位高程误差引起的高低角解算误差的最大值为可近似为:
Figure BSA0000209070230000113
例如,前述地面观察哨观瞄设备的GPS/BD高程定位精度Δh也是≤10m,则由高程定位误差引起的高低角解算误差最大值随标定点距离dab、高程差H变化情况见表2。
表1定位精度Δh=10m时标定点距离、标定点高差对解算高低角误差的影响
Figure BSA0000209070230000114
虽然标定点选择的越远,校准精度越高,但对时间和其他保障条件的要求也越高。因此,校准精度的指标要求往往需要与观瞄设备的性能匹配,以校准误差不使测量误差显著增加即可。例如,需要地面观察哨观瞄设备概略指示目标,要求方位角输出测角精度≤5°,即可要求标定点距离大于300米;若要求高低角误差理想值为≤2°,即标要求标定点与观测点高差不小于30米。选择标定点距离较近时,尽量选择与观测点高差较大的目标点,有利于消除高低角解算误差。
三、磁偏角及观瞄设备测角系统误差一并校准的标定实施方法
利用该标定方法消除观瞄设备系统误差及观测点磁偏角,具体有三种方式:一是观测设备具有数据处理和存储功能时,可将包含磁偏角和方位角系统误差的方位角误差θ和观测设备高低角系统误差ψ解算出来装订存储在观瞄设备中,将每组测量的数据进行校准后发送;二是观测设备具有数据处理功能,不具有储存功能时,可将校准数据θ和ψ发送至后台存储,由后台对各个观察哨实测高低角和磁方位角进行校准;三是观测设备既不具备解算功能,也不具备存储功能时,在后台安装上述算法,在使用前,各观测点将观测点和标定点坐标(经度、纬度、高程),实测标定点高低角和方位角发送至后台,后台根据公式(15)、(16)解算出各观测点的校准数据θ和ψ并储存。本发明的技术方案和实施过程如图6所示。
(一)选择观测点及标定点并分别录取坐标
到达现地后,根据预案和现地实际情况选择标定点P和观测点A,如图1所示。首先,按照标志明显、便于观察和瞄准的要求选择标定点,测量、录入标定点经纬度坐标
Figure BSA0000209070230000121
和高程hp,其中
Figure BSA0000209070230000122
λp分别表示P点的纬度和经度,并做好标记,以便观察瞄准。其次,依据便于观察瞄准标定点、便于观察发现目标、便于隐蔽自身、便于通信联络、交通便利等因素确定观测点,测量、录入观测点经纬度坐标
Figure BSA0000209070230000123
和高程ha,其中
Figure BSA0000209070230000124
λa分别表示A点的纬度和经度。
(二)测量标定点的磁方位角和高低角
在确定标定点P和观测点A并完成定位数据录入之后,使用观测设备的方位和高低测量功能直接测出标定点的磁方位角和高低角的测量数据,并录入系统。其中,磁方位角βc包含方位角系统误差和磁偏角,高低角εc包含高低角系统误差。在测量标定点的方位角和高低角时,如果时间允许,可采用经过多次测量取均值的方法,能够得到比较精准的磁方位角和高低角。
实施标定校准观测时,应保持观测设备稳定、标定点选取距离适中,数据测定准确,为后续精准测算方位角误差θ(包含磁偏角和方位角系统误差)、高低角系统误差ψ提供数据支撑。
(三)解算标定点的真方位角和高低角
根据第一步测定的观测点经纬度坐标
Figure BSA0000209070230000131
和高程ha,以及标定点经纬度坐标
Figure BSA0000209070230000132
和高程hp,依据公式(6-14),解算标定点的真方位角βZ和高低角εZ。为介绍时方便计算,以下将地球长半径R1≈6378245m、短半径R2≈6356863m的地球椭圆近似简化成一个球,其半径近似为:
Figure BSA0000209070230000133
1.解算标定点P点真方位角
如图1所示,在准确测得点A和点P的经纬度坐标后,要测算标定点P的真方位角βz=∠BAC,则先解算出dab和dbc的距离,dab表示A点和B点之间的球面距离,dbc表示B点和C点之间的球面距离。可直接使用求解球面两点间距离的半正矢(haversine)算法,解算出dbc如下,其中haversine算法为本领域技术人员熟知:
Figure BSA0000209070230000134
Figure BSA0000209070230000135
式中,
Figure BSA0000209070230000136
点的地理坐标为
Figure BSA0000209070230000137
表示为B、C两点的纬度差,Δλbc=λbc表示B、C两点的经度差,abc和cbc为中间变量。
则,B、C两点间的球面距离为:
dbc=R·cbc (9)
其中,R为地球半径。
同理,可以直接解算得到A、B两点间的球面距离dab如下:
Figure BSA0000209070230000141
Figure BSA0000209070230000142
dab=R·cab (12)
式中,
Figure BSA0000209070230000143
表示为B、A两点的维度差,Δλab=λba表示B、A两点的经度差,aab和cab为中间变量,dab为A、B两点间的球面距离。
根据上式(9)和(12)解算得到的BC、AB的距离,根据三角函数之反正弦定理可以解算得到标定点P的真方位角为:
Figure BSA0000209070230000144
2.解算标定点P点高低角
由上公式(11)解算得到A、B两点间的距离为dab,根据三角函数之反正切函数定理,解算得标定点P的真实高低角εz如下:
Figure BSA0000209070230000145
式中,hpb为P点和B的高程差,在数值上等于P点和A点的高程差,可通过GPS/北斗设备获得。
(四)消除观瞄设备系统误差及观测点磁偏角
在定位精度和标定点距离满足要求的情况下,可认为依据定位数据解算的方位角βz和高低角εz是准确的,可以此为基准消除测量设备的系统误差,且由于解算的方位角是真方位角,测量的方位角是磁方位角,消除方位角误差的同时也就完成了测量磁方位角转换为真方位角的工作。
根据公式(1),包含磁偏角和方位角系统误差的方位角误差θ为:
θ=βcz (15)
根据公式(2),高低角系统误差为ψ:
ψ=εcz (16)
根据观测点上述校准数据θ和ψ,对测量的每一组数据进行校准。
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc
输出的校准值将大幅度提升设备观测数据的精度,方便观测数据信息入网和应用。
本发明优点如下:
第一,成本低、作用效果好。本发明不增加原有观瞄设备硬件,仅仅通过算法、软件程序及简单操作,就解决了困扰厂家和用户两大难题。一是解决了观瞄设备特别是望远镜等简易的观瞄设备系统误差难以获取和补偿的难题;二是各地磁偏角及其变化情况查询难,基于磁北标定的观瞄设备磁偏角校准的难题。
第二,数学模型成熟、简便,软件实现简单,使用灵活。本发明软件和算法即可全部在终端实现,直接输出校准系统误差后的真方位角,也可输出校准值在后台校准,还能将所有算法和处理软件统一配置在后台,终端只输出坐标、高程、标定角度等标定数据。
第三,使用操作简便易行,省时、省力。本发明可根据精度要求及时间限制、保障条件等灵活选择标点点位置,不需要增加人手、不需要多人配合,耗时很短。其中,标定点的选择和坐标、高程数据录取,可以在到达观测点前途径有显著标志处或人为留下标志,耗时通常在2分钟内,在观测点录取坐标和测量标定角度也不会超过2分钟。
上述优点使本发明在费效比、推广价值、使用前景等方面,具有一定的技术优势。

Claims (4)

1.一种磁偏角及观瞄设备测角误差一并校准方法,其特征在于,具体包括下列步骤:
磁偏角θ直接影响各个观测点的观瞄设备获取观测值的空间一致性;同时,观瞄设备的测角系统误差也需要消除;依据以下方法将二者一并消除;
第一,建立地面直角坐标系;以观测点A为坐标原点,其地理坐标为
Figure FSA0000209070220000011
Figure FSA0000209070220000012
λa分别表示A点的纬度和精度,以A点所在纬线的水平投影指东为Z轴,以所在点经线的水平投影指北为X轴,X轴为水平面向上的法线,即高程轴;AD为磁北方向线;
第二,在观测点A一定距离处选择标定点P
Figure FSA0000209070220000013
利用观瞄设备或手持GPS/北斗终端、手机这样的其他定位手段,分别测量观测点A和标定点P的经纬坐标和高程,B为标定点P在观测点所在水平面XAZ上的投影;
第三,从B点做X轴的垂线相交于C点,依据测量的A和标P的经纬度坐标和高程,利用半正矢haversine算法求解出A、B、C三点两两之间的球面距离dab、dac和dbc;进而利用三角函数关系,计算出标定点P相对于观测点A的真方位角βz=∠CAB,即从经线北顺时针到观察线水平投影的夹角,和高低角εz=∠PAB;
第四,利用观瞄设备从观测点直接测量标定点P包含方位角系统误差的磁方位角βc=∠DAB,即从磁北顺时针到观察线水平投影的夹角,和包含高低角系统误差的高低角εC=∠PAB;
第五,观察标定点P所实际测量得的磁方位角βc减去计算所得的真方位角βz,就得到观瞄设备所在位置磁偏角和观瞄设备方位角系统误差之和θ;所测高低角εc减去计算高低角εZ得到高低角系统误差ψ;
以解算真方位角βz替代实际的真方位角,磁偏角为δ,则磁方位角β=βz-δ;测量设备方位角系统误差为γ,测量所测磁方位角为βc,则磁方位角为β=βc
即:βz-δ=βc-γ,
所以,实测磁方位角减去解算真方位角即为磁偏角和方位角系统误差之和:γ-δ=βcz
即:
θ=γ-δ=βcz (1)
以解算真高低角εz替代实际的真高低角,测量所得P点的高低角为εc,则
ψ=εcz (2)
第六,将包含磁偏角和方位角系统误差的方位角误差θ和高低角系统误差ψ输入观瞄设备相应置数区,用于在输出测量数据时消除观瞄设备的磁偏角和测角误差,使装置直接输出真方位角和高低角,且数值相对精确;
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc-ψ。
2.基于定位精度确定标定点距离的基本方法,其特征在于,具体包括下列步骤:
观测点A与标定点P之间的球面距离为dab,则由水平定位误差Δd引起的真方位角解算误差最大值为:
Δβ=arcsin(2Δd/dab) (3)
同理,当A与P两点在同一水平面,即标定点与观测点高程差为零(H=0)时,由垂直定位误差Δh引起的高低角解算误差最大值为:
Δε=arctan(Δh/dab) (4)
观测点A与标定点P之间的水平距离为dab、高程差为H,则AP之间距离为的
Figure FSA0000209070220000031
由定位高程误差引起的高低角解算误差的最大值近似为:
Figure FSA0000209070220000032
3.如权利要求2所述的基于定位精度确定标定点距离的基本方法,其特征在于,选择标定点距离较近时,尽量选择与观测点高差较大的目标点,有利于消除高低角解算误差。
4.磁偏角及观瞄设备测角系统误差一并校准的标定实施方法,其特征在于,利用所述方法消除观瞄设备系统误差及观测点磁偏角,具体有三种方式:一是观测设备具有数据处理和存储功能时,将包含磁偏角和方位角系统误差的方位角误差θ和观测设备高低角系统误差ψ解算出来装订存储在观瞄设备中,将每组测量的数据进行校准后发送;二是观测设备具有数据处理功能,不具有储存功能时,将校准数据θ和ψ发送至后台存储,由后台对各个观察哨实测高低角和磁方位角进行校准;三是观测设备既不具备解算功能,也不具备存储功能时,在后台安装所述方法相应的算法,在使用前,各观测点将包括经度、纬度、高程在内的观测点和标定点坐标、实测标定点高低角和方位角发送至后台,后台根据公式(15)、(16)解算出各观测点的校准数据θ和ψ并储存;具体如下:
(一)选择观测点及标定点并分别录取坐标
到达现地后,根据预案和现地实际情况选择标定点P和观测点A,首先,按照标志明显、便于观察和瞄准的要求选择标定点,测量、录入标定点经纬度坐标P
Figure FSA0000209070220000041
和高程hp,其中
Figure FSA0000209070220000042
λp分别表示P点的纬度和经度,并做好标记,以便观察瞄准;其次,依据便于观察瞄准标定点、便于观察发现目标、便于隐蔽自身、便于通信联络、交通便利等因素确定观测点,测量、录入观测点经纬度坐标A
Figure FSA0000209070220000043
和高程ha,其中
Figure FSA0000209070220000044
λa分别表示A点的纬度和经度;
(二)测量标定点的磁方位角和高低角
在确定标定点P和观测点A并完成定位数据录入之后,使用观测设备的方位和高低测量功能直接测出标定点的磁方位角和高低角的测量数据,并录入系统;其中,磁方位角βc包含方位角系统误差和磁偏角,高低角εc包含高低角系统误差;在测量标定点的方位角和高低角时,如果时间允许,采用经过多次测量取均值的方法,能够得到比较精准的磁方位角和高低角;
实施标定校准观测时,应保持观测设备稳定、标定点选取距离适中,数据测定准确,为后续精准测算方位角误差θ、高低角系统误差ψ提供数据支撑,其中方位角误差θ包含磁偏角系统误差和方位角系统误差;
(三)解算标定点的真方位角和高低角
根据第一步测定的观测点经纬度坐标A
Figure FSA0000209070220000045
和高程ha,以及标定点经纬度坐标P
Figure FSA0000209070220000046
和高程hp,依据公式(6-14),解算标定点的真方位角βZ和高低角εZ;为介绍时方便计算,以下将椭圆形地球简化成球形地球,其半径近似为:
Figure FSA0000209070220000047
1.解算标定点P点真方位角
准确测得点A和点P的经纬度坐标后,要测算标定点P的真方位角βz=∠BAC,需要先解算出dab和dbc的距离,dab表示A点和B点之间的球面距离,dbc表示B点和C点之间的球面距离;直接使用求解球面两点间距离的半正矢haversine算法,解算出dbc如下:
Figure FSA0000209070220000051
Figure FSA0000209070220000052
式中,C
Figure FSA0000209070220000053
点的地理坐标为
Figure FSA0000209070220000054
表示为B、C两点的纬度差,Δλbc=λbc表示B、C两点的经度差,abc和cbc为中间变量;
则,B、C两点间的球面距离为:
dbc=R·cbc (9)
其中,R为地球半径;
同理,可以直接解算得到A、B两点间的球面距离dab如下:
Figure FSA0000209070220000055
Figure FSA0000209070220000056
dab=R·cab (12)
式中,
Figure FSA0000209070220000057
表示为B、A两点的维度差,Δλab=λba表示B、A两点的经度差,aab和cab为中间变量,dab为A、B两点间的球面距离;
根据上式(9)和(12)解算得到的BC、AB的距离,根据三角函数之反正弦定理解算得到标定点P的真方位角为:
Figure FSA0000209070220000058
2.解算标定点P点高低角
由上公式(11)解算得到A、B两点间的距离为dab,根据三角函数之反正切函数定理,解算得标定点P的真实高低角εz如下:
Figure FSA0000209070220000061
式中,hpb为P点和B的高程差,在数值上等于P点和A点的高程差,能够通过GPS/北斗设备获得;
(四)消除观瞄设备系统误差及观测点磁偏角
在定位精度和标定点距离满足要求的情况下,认为依据定位数据解算的方位角βz和高低角εz是准确的,能够以此为基准消除测量设备的系统误差,且由于解算的方位角是真方位角,测量的方位角是磁方位角,消除方位角误差的同时也就完成了测量磁方位角转换为真方位角的工作;
根据公式(1),
θ=γ-δ=βcz (1)
其中,θ为磁偏角和方位角系统误差之和,γ为测量设备方位角系统误差,δ为磁偏角,βc为测量所得磁方位角,βz为解算的真方位角;
包含磁偏角和方位角系统误差的方位角误差θ为:
θ=βcz (15)
根据公式(2),高低角系统误差为ψ:
ψ=εcz (16)
其中,ψ为高低角系统误差,εc为测量所得P点的高低角,εz为解算真高低角;
根据观测点上述校准数据θ和ψ,对测量的每一组数据进行校准;
观测目标方位角校正值为,
β=βc
观测目标高低角校正值为,
ε=εc-ψ。
CN202010422093.2A 2020-05-08 2020-05-08 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法 Active CN111678536B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010422093.2A CN111678536B (zh) 2020-05-08 2020-05-08 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010422093.2A CN111678536B (zh) 2020-05-08 2020-05-08 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法

Publications (2)

Publication Number Publication Date
CN111678536A true CN111678536A (zh) 2020-09-18
CN111678536B CN111678536B (zh) 2021-12-10

Family

ID=72451993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010422093.2A Active CN111678536B (zh) 2020-05-08 2020-05-08 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法

Country Status (1)

Country Link
CN (1) CN111678536B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113252073A (zh) * 2021-05-12 2021-08-13 无锡商业职业技术学院 一种应用于目标定位系统的现场标定方法及装置
CN113570649A (zh) * 2021-09-26 2021-10-29 南方电网数字电网研究院有限公司 基于三维模型的重力方向确定方法、装置、计算机设备

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6972902B1 (en) * 2004-09-28 2005-12-06 Pacific Telescope Corp. Telescope system having auto-tracking altitude-azimuthal mount and methods for calibrating same
US20050283988A1 (en) * 2004-06-11 2005-12-29 Yamaha Corporation Method and apparatus for measuring magnetic offset of geomagnetic sensor and portable electronic apparatus
US20080163504A1 (en) * 2007-01-05 2008-07-10 Smith John E Apparatus and methods for locating and identifying remote objects
CN103852058A (zh) * 2014-02-18 2014-06-11 中国人民解放军92859部队 一种基于陀螺经纬仪的超短基线磁偏角测量方法
CN107390277A (zh) * 2017-07-19 2017-11-24 中国地震局地球物理研究所 一种高精度自动化绝对地磁观测方法
CN206788370U (zh) * 2017-06-01 2017-12-22 中国人民解放军防空兵学院 一种自动化基线测量系统
CN108152838A (zh) * 2017-12-15 2018-06-12 中国航空工业集团公司洛阳电光设备研究所 一种基于观瞄测定目标位置的装置和方法
CN109541584A (zh) * 2018-12-29 2019-03-29 中国人民解放军空军工程大学 一种基于智能终端的低空飞行器侦察预警系统及方法
CN111006650A (zh) * 2019-11-22 2020-04-14 西安翔迅科技有限责任公司 地面观察哨侦察预警系统
CN111006651A (zh) * 2019-11-22 2020-04-14 西安翔迅科技有限责任公司 目标跟踪数据录取便携装置和观瞄跟踪设备

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050283988A1 (en) * 2004-06-11 2005-12-29 Yamaha Corporation Method and apparatus for measuring magnetic offset of geomagnetic sensor and portable electronic apparatus
US6972902B1 (en) * 2004-09-28 2005-12-06 Pacific Telescope Corp. Telescope system having auto-tracking altitude-azimuthal mount and methods for calibrating same
US20080163504A1 (en) * 2007-01-05 2008-07-10 Smith John E Apparatus and methods for locating and identifying remote objects
CN103852058A (zh) * 2014-02-18 2014-06-11 中国人民解放军92859部队 一种基于陀螺经纬仪的超短基线磁偏角测量方法
CN206788370U (zh) * 2017-06-01 2017-12-22 中国人民解放军防空兵学院 一种自动化基线测量系统
CN107390277A (zh) * 2017-07-19 2017-11-24 中国地震局地球物理研究所 一种高精度自动化绝对地磁观测方法
CN108152838A (zh) * 2017-12-15 2018-06-12 中国航空工业集团公司洛阳电光设备研究所 一种基于观瞄测定目标位置的装置和方法
CN109541584A (zh) * 2018-12-29 2019-03-29 中国人民解放军空军工程大学 一种基于智能终端的低空飞行器侦察预警系统及方法
CN111006650A (zh) * 2019-11-22 2020-04-14 西安翔迅科技有限责任公司 地面观察哨侦察预警系统
CN111006651A (zh) * 2019-11-22 2020-04-14 西安翔迅科技有限责任公司 目标跟踪数据录取便携装置和观瞄跟踪设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李飞战等: "不同磁差测定方法的研究", 《测绘技术装备》 *
田晋等: "地磁偏角测量GPS基准点误差分析", 《海洋测绘》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113252073A (zh) * 2021-05-12 2021-08-13 无锡商业职业技术学院 一种应用于目标定位系统的现场标定方法及装置
CN113570649A (zh) * 2021-09-26 2021-10-29 南方电网数字电网研究院有限公司 基于三维模型的重力方向确定方法、装置、计算机设备

Also Published As

Publication number Publication date
CN111678536B (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
US8411285B2 (en) Stationing an unleveled optical total station
CN105445774B (zh) 一种gnss与激光测距相结合的测量系统及测量方法
CN108051835B (zh) 一种基于双天线的倾斜测量装置及测量与放样方法
US5774826A (en) Optimization of survey coordinate transformations
US20120290199A1 (en) Apparatus, system and method for self orientation
US20120026322A1 (en) Method, tool, and device for determining the coordinates of points on a surface by means of an accelerometer and a camera
CN108278968A (zh) 一种车载扫描系统控制点检校方法
JP2005331499A (ja) 単又は複数不安定ゾーンを有するサイトにおける地上ベース測量方法および装置
CN108151698B (zh) 一种基于轴线相交法的天线旋转中心标校方法
CN111678536B (zh) 地面观察哨磁偏角及观瞄设备测角系统误差一并校准的标定方法
CN103697885B (zh) 自动补偿磁偏角的远程定位方法
CN111190204B (zh) 基于北斗双天线和激光测距仪的实时定位装置及定位方法
CN108168504B (zh) 一种基于拟合球心法的天线旋转中心标校方法
CN109883406A (zh) 基于更少点的全站仪概略位姿的计算方法及系统
CN104049269A (zh) 一种基于激光测距和mems/gps组合导航系统的目标导航测绘方法
CN112363191A (zh) 一种基于rtk的野外森林样地定位方法
CN216925591U (zh) 基于动态实时定位的便携式激光测量设备
US20190186913A1 (en) Method for verifying accuracy of rtk tilt compensation measurement
EP2040029A1 (en) A multi mode active surveying pole
CN105806304B (zh) 天线方向角测量方法及装置
Dumrongchai et al. Performance tests of geodetic receivers with tilt sensors in obstructed environments using the NRTK GNSS technique
CN207689674U (zh) 一种基于观瞄测定目标位置的装置
US20240069211A1 (en) Method for determining, using an optronic system, positions in a scene, and associated optronic system
KR101149348B1 (ko) 단독 gps 측량을 이용한 실시간 공간정보 위치 정확도 평가 시스템 및 평가방법
CN113447908B (zh) 一种固定式雷达探测目标的gps北斗坐标获取方法以及装置

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