CN114061382B - 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法 - Google Patents

一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法 Download PDF

Info

Publication number
CN114061382B
CN114061382B CN202111090114.6A CN202111090114A CN114061382B CN 114061382 B CN114061382 B CN 114061382B CN 202111090114 A CN202111090114 A CN 202111090114A CN 114061382 B CN114061382 B CN 114061382B
Authority
CN
China
Prior art keywords
observation station
vector
central axis
axis
mid
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
CN202111090114.6A
Other languages
English (en)
Other versions
CN114061382A (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.)
Pla 63875 Unit
Original Assignee
Pla 63875 Unit
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 Pla 63875 Unit filed Critical Pla 63875 Unit
Priority to CN202111090114.6A priority Critical patent/CN114061382B/zh
Publication of CN114061382A publication Critical patent/CN114061382A/zh
Application granted granted Critical
Publication of CN114061382B publication Critical patent/CN114061382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B35/00Testing or checking of ammunition
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C3/00Measuring distances in line of sight; Optical rangefinders

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供了一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,获取仿真需求段对应的目标偏航角和俯仰角;将经过姿态旋转后的中轴三维矢量投影至两个观测站,得到对应的二维矢量和像面倾角;根据观测站的焦距以及二维矢量得到两个观测站的物距,以及目标中轴对应的观测站中轴成像像长,进而得到图像倾角判读最大误差角量;计算所求中轴线矢量理论值和中轴线矢量因像面倾角误差引起的主导误差量;最终得到实际偏航角和俯仰角的量化精度。本发明更加方便、直观、准确,而且解决了关键参数未知的难题,为靶场姿态测量提供了比较可靠的事前布站测试、精度预估以及事后分析的理论支撑。

Description

一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法
技术领域
本发明涉及兵器试验测试领域,特别是涉及靶场光学姿态测量领域。
背景技术
姿态参数在武器系统飞行稳定性、攻击精度等方面有重要作用。内测需要内置传感器,会提高测量成本,而且低高度区信号接收困难也是内测的劣势之一。光学姿态测量是重要的外测手段。目前光学姿态测量主要采用姿态测量系统来实施,一般基于中远距离(作用距离至少大于1公里)。姿态测量实施前需要设计测试方案,精度部分难以量化,主要原因在于传统精度分析采用偏导数分析法,涉及的具体成像测量参数未知,比如分站像面倾角、截距等,故精度评估仅靠测量距离所成大致像长进行大致评估,与实际测量场景相差较远,造成实际姿态测量方案设计缺乏必需的精度理论依据。
发明内容
为了克服现有技术的不足,本发明提供一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,利用两站图像信息和关键段弹道信息(姿态信息)获取目标姿态测量关键段的姿态精度,为靶场典型目标姿态交会测量提供精度评估方法。
本发明解决其技术问题所采用的技术方案包括以下步骤:
第一步,获取仿真需求段对应的目标偏航角ψ和目标俯仰角
Figure BDA0003267093230000011
第二步,将经过姿态旋转后的中轴三维矢量
Figure BDA0003267093230000012
投影至观测站1和观测站2,分别得到二维矢量(y1 z1)T和(y2 z2)T,并获取两个观测站对应的像面倾角α1、α2
第三步,根据观测站的焦距f1、f2,以及二维矢量(y1 z1)T、(y2 z2)T得到两个观测站的物距分别为u1和u2,以及长度为L的目标中轴对应的观测站中轴成像像长,进而得到图像倾角判读最大误差角量Δα;
第四步,计算所求中轴线矢量理论值l(l m n)T和中轴线矢量因像面倾角α误差引起的主导误差量Δl(l m n)T
第五步,计算实际偏航角ψ和俯仰角
Figure BDA0003267093230000013
的量化精度Δψ和
Figure BDA0003267093230000014
所述的第二步中,i取值为1或2,(Ai,Ei)为观测站i的光轴指向,
Figure BDA0003267093230000021
Figure BDA0003267093230000022
所述的第三步中,观测站Y向像长
Figure BDA0003267093230000023
观测站Z向像长
Figure BDA0003267093230000024
分站图像判读最大误差角量
Figure BDA0003267093230000025
所述的第四步中,
Figure BDA0003267093230000026
所述的第四步中,
Figure BDA0003267093230000027
所述的第五步中,
Figure BDA0003267093230000028
本发明的有益效果是:按照姿态模拟、透视投影、关键参数误差量获取、中轴线法重建、获取姿态角精度的过程,确立了中远距离下闭环反馈可执行的中轴姿态交会测量的精度评估方法。相较于传统偏导数方法,本发明不仅突出了关键参数的影响,更加方便、直观、准确,而且解决了关键参数未知的难题,为靶场姿态测量提供了比较可靠的事前布站测试、精度预估以及事后分析的理论支撑。
附图说明
图1是本发明的方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明以中轴线法为依据,根据靶场中远距离下姿态测量主要影响因素进行重点分析,直观准确地突出了光学姿态测量中方向判读误差对姿态测量精度的影响,并将目标中轴在空间向像面进行模拟投影,解决了传统方法关键参数未知的难题,整个仿真场景按照姿态模拟、透视投影、关键参数误差量获取、中轴线法重建、获取姿态角精度的过程进行,首次确立了可执行的闭环反馈的姿态交会仿真方法。本发明为靶场姿态测量提供了比较可靠的事前布站测试、精度预估以及事后分析提供理论支撑。本发明还能够拓展用于其他姿态测量相关领域。
本发明提供一种基于中远距离(适用距离不小于1公里)下中轴矢量交会姿态测量的精度模型,获取该精度模型包括下述步骤:
第一步:调整原中轴线交会法形式,将中轴线交会矢量叉积变换为矢量点积;
所有涉及坐标系以常规姿态处理算法中涉及的坐标系为基准。其中目标坐标系以目标中轴为XT轴,头部为正;基准坐标系以北向为X轴;摄像机坐标系以光轴指向为XC轴,平行于像面向上为YC轴;像面坐标系y、Z轴分别平行于像机坐标系YC轴、ZC轴。初始零姿态指目标坐标系平行于基准坐标系。
设分站1光轴指向为(A1,E1),点位坐标为(x1,y1,z1)、分站2光轴指向为(A2,E2),点位坐标为(x2,y2,z2);目标中轴在分站1的像面倾角为ɑ1弧度,截距为P1;分站1焦距为f1mm和分站2焦距为f2mm;目标中轴在分站2的像面倾角为ɑ2弧度,截距为P2。基准以常规中轴线法为准。则常规中轴线法定义中轴线矢量l(l m n)T为分站法线矢量l1(a1 b1 c1)T、l2(a2 b2c2)T的叉积。
Figure BDA0003267093230000031
其中:
Figure BDA0003267093230000032
变叉积形式为点积形式:
Figure BDA0003267093230000041
第二步:根据第一步获取所求中轴线矢量理论值l(l m n)T
中长远距离下姿态测量P/f量级约为0.001,忽略掉P/f项。则
Figure BDA0003267093230000042
第三步:获取中轴线矢量因分站像面倾角α误差引起的主导误差量Δl(ΔlΔmΔn)T
基于图像的中长远姿态测量中,主导误差为分站像面倾角ɑ提取误差,设分站2为主要判读误差来源,将l对α2求导,则
Figure BDA0003267093230000043
第四步:获取目标偏航角、俯仰角量化精度Δψ、
Figure BDA0003267093230000044
Figure BDA0003267093230000045
本发明还提供一种利用上述模型的中轴矢量交会姿态测量的精度预估仿真方法,包括下述步骤:
第一步:获取仿真需求段对应的目标偏航角ψ和俯仰角
Figure BDA0003267093230000046
第二步:将经过姿态旋转后的中轴三维矢量
Figure BDA0003267093230000047
投影至分站1和分站2成二维矢量(y1 z1)T、(y2 z2)T,并获取分站倾角α1、α2
Figure BDA0003267093230000048
Figure BDA0003267093230000051
第三步:根据分站焦距信息f1、f2,及第二步获取的二维矢量信息(y1 z1)T、(y2z2)T,弹道信息和分站点位信息可得分站物距u1米、u2米、可得目标中轴(长度L米)对应的分站中轴成像像长,计算获取图像倾角判读最大误差角量Δα;
分站Y向像长lY
Figure BDA0003267093230000052
分站Z向像长lZ
Figure BDA0003267093230000053
分站总像长
Figure BDA0003267093230000054
分站图像判读最大误差角量:
Figure BDA0003267093230000055
第四步:根据上述精度模型的第二步、第三步获取l(l m n)T、Δl(ΔlΔmΔn)T
第五步:根据上述精度模型的第四步获取实际偏航角ψ、俯仰角
Figure BDA0003267093230000056
的量化精度Δψ、
Figure BDA0003267093230000057
本发明的实施例包括以下步骤:
第一步:获取仿真需求段对应的目标理论偏航角ψ和俯仰角
Figure BDA0003267093230000058
具体理论姿态角信息如表1所示;
表1
Figure BDA0003267093230000059
第二步:根据分站关键点指向(AE)信息,如表2所示,将经过姿态旋转后的中轴三维矢量
Figure BDA00032670932300000510
投影至分站1和分站2成二维矢量(y1z1)T、(y2z2)T,并获取分站倾角α1、α2,结果如表3所示;
表2
Figure BDA0003267093230000061
表3
Figure BDA0003267093230000062
第三步:根据分站焦距信息1500mm,及第二步获取的二维矢量信息(y1 z1)T、(y2z2)T,弹道信息和分站点位信息可得分站物距u1米、u2米、可得目标中轴(长度2.45米)对应的分站中轴成像像长l2#,计算获取图像倾角判读最大误差角量,此处以分站2为例计算,结果如表4所示;
表4
Figure BDA0003267093230000063
第四步:根据发明内容1中第二、三步所述,计算获取关键点l(l m n)T,以分站2作为主要误差源获取Δl(ΔlΔmΔn)T,结果如表5所示;
表5
Figure BDA0003267093230000064
Figure BDA0003267093230000071
第五步:根据发明内容1中第四步所述,获取实际偏航角ψ、俯仰角
Figure BDA0003267093230000072
的量化精度Δψ、
Figure BDA0003267093230000073
结果如表6所示。
表6
Figure BDA0003267093230000074

Claims (6)

1.一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,包括以下步骤:
第一步,提供一种基于中远距离下中轴矢量交会姿态测量的精度模型,获取该精度模型包括下述步骤:
步骤1.1:调整原中轴线交会法形式,将中轴线交会矢量叉积变换为矢量点积;
所有涉及坐标系以常规姿态处理算法中涉及的坐标系为基准,其中目标坐标系以目标中轴为XT轴,头部为正;基准坐标系以北向为X轴;摄像机坐标系以光轴指向为XC轴,平行于像面向上为YC轴;像面坐标系y、Z轴分别平行于像机坐标系YC轴、ZC轴,初始零姿态指目标坐标系平行于基准坐标系,
设观测站1光轴指向为(A1,E1),点位坐标为(x1,y1,z1)、观测站2光轴指向为(A2,E2),点位坐标为(x2,y2,z2);目标中轴在观测站1的像面倾角为ɑ1弧度,截距为P1;观测站1焦距为f1mm和观测站2焦距为f2mm;目标中轴在观测站2的像面倾角为ɑ2弧度,截距为P2,基准以常规中轴线法为准,则常规中轴线法定义中轴线矢量l(l m n)T为观测站法线矢量l1(a1 b1c1)T、l2(a2 b2 c2)T的叉积,
其中:
变叉积形式为点积形式:
步骤1.2:根据第一步获取所求中轴线矢量理论值l(lm n)T
中长远距离下姿态测量P/f量级约为0.001,忽略掉P/f项,则
步骤1.3:获取中轴线矢量因观测站像面倾角α误差引起的主导误差量Δl(ΔlΔmΔn)T
基于图像的中长远姿态测量中,主导误差为观测站像面倾角ɑ提取误差,设观测站2为主要判读误差来源,将l对α2求导,则
步骤1.4:获取目标偏航角、俯仰角量化精度Δψ、
基于精度模型,获取仿真需求段对应的目标偏航角ψ和目标俯仰角
第二步,将经过姿态旋转后的中轴三维矢量投影至观测站1和观测站2,分别得到二维矢量(y1 z1)T和(y2 z2)T,并获取两个观测站对应的像面倾角α1、α2
第三步,根据观测站的焦距f1、f2,以及二维矢量(y1 z1)T、(y2 z2)T得到两个观测站的物距分别为u1和u2,以及长度为L的目标中轴对应的观测站中轴成像像长,进而得到图像倾角判读最大误差角量Δα;
第四步,计算所求中轴线矢量理论值l(lm n)T和中轴线矢量因像面倾角α误差引起的主导误差量Δl(lm n)T
第五步,计算实际偏航角ψ和俯仰角的量化精度Δψ和
2.根据权利要求1所述的基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,所述的第二步中,(A1,E1)为观测站1的光轴指向,(A2,E2)为观测站2的光轴指向,
3.根据权利要求1所述的基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,所述的第三步中,观测站Y向像长观测站Z向像长观测站图像判读最大误差角量
4.根据权利要求1所述的基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,所述的第四步中,
5.根据权利要求1所述的基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,所述的第四步中,
6.根据权利要求1所述的基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法,其特征在于,所述的第五步中,
CN202111090114.6A 2021-09-17 2021-09-17 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法 Active CN114061382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111090114.6A CN114061382B (zh) 2021-09-17 2021-09-17 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111090114.6A CN114061382B (zh) 2021-09-17 2021-09-17 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法

Publications (2)

Publication Number Publication Date
CN114061382A CN114061382A (zh) 2022-02-18
CN114061382B true CN114061382B (zh) 2023-05-12

Family

ID=80233767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111090114.6A Active CN114061382B (zh) 2021-09-17 2021-09-17 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法

Country Status (1)

Country Link
CN (1) CN114061382B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8033221B2 (en) * 2007-08-13 2011-10-11 Raytheon Company System and method for sensing proximity
CN105910607B (zh) * 2016-04-07 2018-11-13 国家测绘地理信息局卫星测绘应用中心 基于地面控制的卫星长周期姿态误差修正方法
CN107314718B (zh) * 2017-05-31 2018-11-13 中北大学 基于磁测滚转角速率信息的高速旋转弹姿态估计方法
CN108225276B (zh) * 2018-01-11 2020-06-23 中国人民解放军战略支援部队航天工程大学 一种单星成像目标运动特性反演方法及系统
CN108444449B (zh) * 2018-02-02 2019-03-08 中国科学院西安光学精密机械研究所 一种对具有平行线特征的目标空间姿态测量方法
CN112598617B (zh) * 2020-11-13 2023-09-01 北京遥测技术研究所 一种基于虚拟平台的外弹道光学测量精度分析方法

Also Published As

Publication number Publication date
CN114061382A (zh) 2022-02-18

Similar Documents

Publication Publication Date Title
CN110296691A (zh) 融合imu标定的双目立体视觉测量方法与系统
CN109297436B (zh) 双目线激光立体测量基准标定方法
CN102472609B (zh) 位置和姿势校准方法及设备
EP1870856B1 (en) Information-processing method and apparatus for calculating information regarding measurement target on the basis of captured images
CN105423954B (zh) 一种基于视觉测量的柔壁喷管的测量方法
CN103559711A (zh) 基于三维视觉系统图像特征和三维信息的运动估计方法
CN111678521B (zh) 一种移动机器人定位精度的评估方法和系统
CN110849331B (zh) 基于三维点云数据库模型的单目视觉测量与地面试验方法
CN112729221A (zh) 一种飞机舵面偏转角度测量方法
CN110595479B (zh) 一种基于icp算法的slam轨迹评估方法
CN112815843B (zh) 一种3d打印过程中工件表面打印偏差的在线监测方法
CN113689473B (zh) 基于变分原理的河道表面流速计算方法及装置
CN109636927A (zh) 一种飞行器姿态测量算法训练与鉴定的系统及方法
CN111915685B (zh) 一种变焦摄像机标定方法
CN111307146A (zh) 一种基于双目相机和imu的虚拟现实头戴显示设备定位系统
CN112710303A (zh) 由运动平台运动引起目标在视场中姿态角θ变化的确定方法
CN109035343A (zh) 一种基于监控相机的楼层相对位移测量方法
CN108257184B (zh) 一种基于正方形点阵合作目标的相机姿态测量方法
CN113393413B (zh) 基于单目与双目视觉协同的水域面积测量方法和系统
CN114061382B (zh) 一种基于中远距离下中轴矢量交会姿态测量的精度预估仿真方法
CN107945166B (zh) 基于双目视觉的待测物体三维振动轨迹的测量方法
CN112712566A (zh) 基于结构参数在线校正的双目立体视觉传感器测量方法
CN107421476A (zh) 一种空间孔位测量基准误差补偿方法
CN107515012B (zh) 基于单轴旋转机构的动态视觉测量系统校准装置及方法
CN113240749B (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