CN108050960A - 一种基于数字摄影测量技术的高精度转动测量方法 - Google Patents

一种基于数字摄影测量技术的高精度转动测量方法 Download PDF

Info

Publication number
CN108050960A
CN108050960A CN201711380705.0A CN201711380705A CN108050960A CN 108050960 A CN108050960 A CN 108050960A CN 201711380705 A CN201711380705 A CN 201711380705A CN 108050960 A CN108050960 A CN 108050960A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
mtd
mtr
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
CN201711380705.0A
Other languages
English (en)
Other versions
CN108050960B (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.)
Purple Mountain Observatory of CAS
Original Assignee
Purple Mountain Observatory of CAS
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 Purple Mountain Observatory of CAS filed Critical Purple Mountain Observatory of CAS
Priority to CN201711380705.0A priority Critical patent/CN108050960B/zh
Publication of CN108050960A publication Critical patent/CN108050960A/zh
Application granted granted Critical
Publication of CN108050960B publication Critical patent/CN108050960B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提供了一种基于数字摄影测量技术的高精度转动测量方法,并包括如下步骤:一、在转动机架的转子和定子上分别粘贴反光靶标,并将定子和转子上的反光靶标数量分别为Ns和Nr,其中,Ns和Nr均为正整数;二、在初始未发生转动的条件下,对转动机架上的反光靶标进行拍照,并在拍摄完毕后,根据拍摄的照片对各个反光靶标进行解算三维坐标;三、将转动机架的转子转动到Na组待测量角度上,分别进行拍摄并解算三维坐标,得到每一待测量角度下各个反光靶标的三维坐标;四、根据步骤二和步骤三得到的各个反光靶标的三维坐标

Description

一种基于数字摄影测量技术的高精度转动测量方法
技术领域
本发明属于图像信息技术处理领域,具体地涉及一种基于数字摄影测量技术的高精度转动测量方法。
背景技术
现代天文望远镜对指向精度的要求日益提高。在射电领域,大型射电天文望远镜(如FAST望远镜、德令哈13.7米口径毫米波望远镜等)的主波束宽度的典型值为10-50角秒,为了保证观测效率,通常要求望远镜指向精度达到×1角秒量级。
望远镜指向误差通常包含系统偏差部分和随机误差部分,其中系统偏差部分可以通过指向模型得以修正,修正后的残余误差则主要为随机指向误差,这部分误差主要由望远镜机架的重复定位精度决定。为了能够对望远镜机架的重复定位能力进行评估,需要对望远镜各轴系在给定转动指令条件下的实际转动角度进行精确的测量,测量精度往往需要达到1-2角秒甚至是亚角秒量级(×0.1角秒)。
在天文领域,用于高精度转动测量的设备和方法包括经纬仪、全站仪、激光跟踪仪、角摆仪等。
经纬仪和全站仪是直接的测角设备,其缺点是需要人为对准读数,引起的读数误差通常就达到了角秒量级,无法满足高精度转动测量的要求。此外,这类设备需要通过工装架设在机架上,随机架一起转动,这对小型机架的测量会造成不便。经纬仪正常工作时要求水平姿态,因此也会限制其在空间三维转动测量中的应用。
激光跟踪仪通常有较高的测距精度,但测角精度较差,极限精度通常在×1角秒量级。在使用激光跟踪仪测量转动角度时,通常将反射靶标置于待测的机架的转子上,通过测量靶标相对于激光跟踪仪的转动角度来推算机架的转动角度。由于激光跟踪仪并不位于机架的转动中心,在计算过程中有可能会将角度测量误差放大。
角摆仪的测量精度较高,可以达到1角秒甚至更高的角分辨率,但其缺点是需要配合激光干涉仪使用。此外,角摆仪只能测量平面内的转动角,无法实现空间三维任意方向轴上的转动角测量。
发明内容
本发明的目的在于针对现有技术的缺陷,提供一种基于数字摄影测量技术的高精度转动测量方法。
本发明的技术方案如下:一种基于数字摄影测量技术的高精度转动测量方法包括如下步骤:一、在转动机架的转子和定子上分别粘贴反光靶标,并将定子和转子上的反光靶标数量分别为Ns和Nr,其中,Ns和Nr均为正整数;二、在初始未发生转动的条件下,对转动机架上的反光靶标进行拍照,并在拍摄完毕后,根据拍摄的照片对各个反光靶标进行解算三维坐标,记作:其中,i=1,2,…,Nr;三、将转动机架的转子转动到Na组待测量角度上,分别进行拍摄并解算三维坐标,得到每一待测量角度下各个反光靶标的三维坐标,并记作:
其中,i=1,2,…,Nr,j=1,2,…,Na,Na为正整数,分别为反光
靶标的坐标数值;四、根据步骤二和步骤三得到的各个反光靶标的三维坐标拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度。
优选地,在步骤四中,拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度的步骤中,利用最小二乘法找到使得拟合函数达到极小值的转动角θj,拟合函数为:
其中,i=1,2,…,Nr,j=1,2,…,Na,|| ||为求模运算,R(θj)为坐标旋转矩阵,为每一个反光靶标的旋转中心位置矢量。
优选地,在步骤四中反光靶标的旋转中心位置矢量的拟合包括如下步骤:
拟合每一个反光靶标所在的旋转平面Aix+Biy+Ciz=1,
拟合每一个反光靶标的旋转轴方向矢量
拟合转动机架的旋转轴
拟合每一个反射靶标的旋转中心位置矢量
优选地,在步骤四中,R(θj)为坐标旋转矩阵,其定义为:
本发明提供的技术方案具有如下有益效果:
1、需要的测量设备简单,如数字摄影相机,对测量场景没有特殊要求,适用性广;
2、测量设备不需要架设在待测转动机构上,因此不需对待测转动机构做任何结构上的改造;
3、可以方便地实现空间三维任意方向轴上的转动角测量;
4、可以连续对多种转动姿态进行拍摄,拍摄完毕后通过数据处理可一次性获得多个转动姿态下的转动角;
5、反光靶标在待测转动机构上粘贴的位置可以随机分布,但为了获得更高的测量精度,需将靶标尽量均匀地分布在转动圆周方向上,并且尽可能远离转动轴;增加靶标的数量可以有效地提高角度测量精度;而且,在实际应用中,可以根据最终所需的测角精度来确定靶标数目。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
除非上下文另有特定清楚的描述,本发明中的元件和组件,数量既可以单个的形式存在,也可以多个的形式存在,本发明并不对此进行限定。本发明中的步骤虽然用标号进行了排列,但并不用于限定步骤的先后次序,除非明确说明了步骤的次序或者某步骤的执行需要其他步骤作为基础,否则步骤的相对次序是可以调整的。可以理解,本文中所使用的术语“和/或”涉及且涵盖相关联的所列项目中的一者或一者以上的任何和所有可能的组合。
本发明实施例提供的基于数字摄影测量技术的高精度转动测量方法可以用于任意转动机架,例如望远镜机架、斜轴式太赫兹天线机架等天文领域需要精确调整角度的机架设备。
需要说明的是,所述转动机架的包括定子和转子,所述定子为固定部分,所述转子为可转动部分,并且转子可以相对于定子进行转动,从而调节转动机架的角度调节。
具体地,所述基于数字摄影测量技术的高精度转动测量方法包括如下步骤:
一、在转动机架的转子和定子上分别粘贴反光靶标,并将定子和转子上的反光靶标数量分别为Ns和Nr,其中,Ns和Nr均为正整数。
在步骤一中,具体地,反光靶标的粘贴位置可以随机分布,但要求尽量均匀分布在转动圆周方向上,并且尽可能远离转动轴。
二、在初始未发生转动的条件下,对转动机架上的反光靶标进行拍照,并在拍摄完毕后,根据拍摄的照片对各个反光靶标进行解算三维坐标,记作:其中,i=1,2,…,Nr。
在步骤二中,具体地,将得到的各个反光靶标的三维坐标作为后续转动测量的参照值。
三、将转动机架的转子转动到Na组待测量角度上,分别进行拍摄并解算三维坐标,得到每一待测量角度下各个反光靶标的三维坐标,并记作:
其中,i=1,2,…,Nr,j=1,2,…,Na,Na为正整数,分别为反光靶标的坐标数值。
四、根据步骤二和步骤三得到的各个反光靶标的三维坐标拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度。
具体地,在步骤四的拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度的步骤中,利用最小二乘法找到使得拟合函数达到极小值的转动角θj,拟合函数为:
其中,i=1,2,…,Nr,j=1,2,…,Na,|| ||为求模运算,R(θj)为坐标旋转矩阵,为每一个反光靶标的旋转中心位置矢量。
而且,在步骤四中反光靶标的旋转中心位置矢量的拟合包括如下步骤:
拟合每一个反光靶标所在的旋转平面Aix+Biy+Ciz=1,
拟合每一个反光靶标的旋转轴方向矢量
拟合转动机架的旋转轴
拟合每一个反射靶标的旋转中心位置矢量
而且,R(θj)为坐标旋转矩阵,其定义为:
例如,利用本发明实施例提供的基于数字摄影测量技术的高精度转动测量方法对斜轴式太赫兹天线机架进行了重复定位精度测量。
斜轴式太赫兹天线机架高约1.8米,底部直径约700毫米。机架轴系由一个方位轴和一个斜轴构成,两轴之间成45度交角。在机架底部的三脚支架(定子)和顶部平台(转子)上分别粘贴了15个和36个反光靶标,其他的外壁部分均匀粘贴了用于图像拼接的编码标志。待测姿态分别为:方位轴0°,120°和240°,斜轴0°,90°和180°。
首先在计算机上对测量过程进行了仿真,仿真中假设靶标坐标测量的误差为随机误差,并假设其均方根为3μm。利用本专利提出的方法对测量流程进行模拟仿真,最终得到角度测量误差的均方根为0.4角秒。
接下来利用商用数字摄影测量系统对天线机架的重复定位精度进行了实测。通过测量得到方位轴和斜轴的重复定位精度分别为3.3角秒RMS和3.9角秒RMS。测量系统自身的坐标测量精度约为3μm RMS,角度重复测量精度优于0.5角秒RMS,与仿真结果相符。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (4)

1.一种基于数字摄影测量技术的高精度转动测量方法,其特征在于:包括如下步骤:
一、在转动机架的转子和定子上分别粘贴反光靶标,并将定子和转子上的反光靶标数量分别为Ns和Nr,其中,Ns和Nr均为正整数;
二、在初始未发生转动的条件下,对转动机架上的反光靶标进行拍照,并在拍摄完毕后,根据拍摄的照片对各个反光靶标进行解算三维坐标,记作:其中,i=1,2,…,Nr;
三、将转动机架的转子转动到Na组待测量角度上,分别进行拍摄并解算三维坐标,得到每一待测量角度下各个反光靶标的三维坐标,并记作:
<mrow> <msubsup> <mover> <mi>P</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> <mi>j</mi> </msubsup> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
其中,i=1,2,…,Nr,j=1,2,…,Na,Na为正整数,分别为反光靶标的坐标数值;
四、根据步骤二和步骤三得到的各个反光靶标的三维坐标拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度。
2.根据权利要求1所述的一种基于数字摄影测量技术的高精度转动测量方法,其特征在于:在步骤四中,拟合计算在给定转动指令的条件下,转子相对于定子的相对转动角度的步骤中,利用最小二乘法找到使得拟合函数达到极小值的转动角θj,拟合函数为:
<mrow> <msub> <mi>f</mi> <mi>j</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>r</mi> </mrow> </munderover> <mo>|</mo> <mo>|</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msubsup> <mover> <mi>P</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> <mn>0</mn> </msubsup> <mo>-</mo> <msub> <mover> <mi>O</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>O</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msubsup> <mover> <mi>P</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> <mi>j</mi> </msubsup> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>,</mo> </mrow>
其中,i=1,2,…,Nr,j=1,2,…,Na,|| ||为求模运算,R(θj)为坐标旋转矩阵,为每一个反光靶标的旋转中心位置矢量。
3.根据权利要求2所述的一种基于数字摄影测量技术的高精度转动测量方法,其特征在于:在步骤四中反光靶标的旋转中心位置矢量的拟合包括如下步骤:
拟合每一个反光靶标所在的旋转平面Aix+Biy+Ciz=1,
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>C</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
拟合每一个反光靶标的旋转轴方向矢量
<mrow> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>g</mi> <mi>i</mi> </msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>A</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>B</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>C</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
<mrow> <msub> <mi>g</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <msubsup> <mi>A</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>B</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>C</mi> <mi>i</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mfrac> <mo>;</mo> </mrow>
拟合转动机架的旋转轴
<mrow> <mover> <mi>v</mi> <mo>^</mo> </mover> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>v</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>y</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>v</mi> <mi>z</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mi>r</mi> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>N</mi> <mi>r</mi> </mrow> </munderover> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>;</mo> </mrow>
拟合每一个反射靶标的旋转中心位置矢量
<mrow> <msub> <mover> <mi>O</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msub> <mover> <mi>S</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mover> <mi>v</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> <msub> <mover> <mi>S</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> </mrow> <msub> <mi>g</mi> <mi>i</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mover> <mi>v</mi> <mo>^</mo> </mover> <mo>,</mo> </mrow>
<mrow> <msub> <mover> <mi>S</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> </mrow> </mtd> <mtd> <mrow> <mi>N</mi> <mi>a</mi> <mo>+</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;CenterDot;</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>x</mi> <mi>i</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>y</mi> <mi>i</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <msubsup> <mi>z</mi> <mi>i</mi> <mi>j</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mi>a</mi> </mrow> </munderover> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>i</mi> <msup> <mi>j</mi> <mn>2</mn> </msup> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
4.根据权利要求3所述的一种基于数字摄影测量技术的高精度转动测量方法,其特征在于:在步骤四中,R(θj)为坐标旋转矩阵,其定义为:
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>+</mo> <msubsup> <mi>v</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>v</mi> <mi>x</mi> </msub> <msub> <mi>v</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>v</mi> <mi>z</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>v</mi> <mi>x</mi> </msub> <msub> <mi>v</mi> <mi>z</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>v</mi> <mi>y</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>v</mi> <mi>y</mi> </msub> <msub> <mi>v</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>v</mi> <mi>z</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>+</mo> <msubsup> <mi>v</mi> <mi>y</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>v</mi> <mi>y</mi> </msub> <msub> <mi>v</mi> <mi>z</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>v</mi> <mi>x</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>v</mi> <mi>z</mi> </msub> <msub> <mi>v</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>v</mi> <mi>y</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>v</mi> <mi>z</mi> </msub> <msub> <mi>v</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>v</mi> <mi>x</mi> </msub> <msub> <mi>sin&amp;theta;</mi> <mi>j</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>+</mo> <msubsup> <mi>v</mi> <mi>z</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>cos&amp;theta;</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
CN201711380705.0A 2017-12-20 2017-12-20 一种基于数字摄影测量技术的高精度转动测量方法 Active CN108050960B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711380705.0A CN108050960B (zh) 2017-12-20 2017-12-20 一种基于数字摄影测量技术的高精度转动测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711380705.0A CN108050960B (zh) 2017-12-20 2017-12-20 一种基于数字摄影测量技术的高精度转动测量方法

Publications (2)

Publication Number Publication Date
CN108050960A true CN108050960A (zh) 2018-05-18
CN108050960B CN108050960B (zh) 2020-06-02

Family

ID=62130227

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711380705.0A Active CN108050960B (zh) 2017-12-20 2017-12-20 一种基于数字摄影测量技术的高精度转动测量方法

Country Status (1)

Country Link
CN (1) CN108050960B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110132577A (zh) * 2019-05-29 2019-08-16 河北科技大学 齿轮箱传动性能检测系统及齿轮箱传动性能检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1597539A1 (de) * 2003-02-27 2005-11-23 Storz Endoskop Produktions GmbH Verfahren und optisches system zur vermessung der topographie eines messobjekts
US20070150194A1 (en) * 2003-03-31 2007-06-28 Gleb Chirikov Method for navigation with optical sensors, and a device utilizing the method
CN101226057A (zh) * 2008-02-01 2008-07-23 武汉朗视软件有限公司 一种数字近景摄影测量方法
CN101387501A (zh) * 2008-10-06 2009-03-18 天津大学 超大型工件圆形截面形状与方位测量装置及方法
CN101427155A (zh) * 2006-04-21 2009-05-06 法罗技术股份有限公司 具有能够旋转的反射镜的基于摄影机的六自由度标靶测量和标靶跟踪设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1597539A1 (de) * 2003-02-27 2005-11-23 Storz Endoskop Produktions GmbH Verfahren und optisches system zur vermessung der topographie eines messobjekts
US20070150194A1 (en) * 2003-03-31 2007-06-28 Gleb Chirikov Method for navigation with optical sensors, and a device utilizing the method
CN101427155A (zh) * 2006-04-21 2009-05-06 法罗技术股份有限公司 具有能够旋转的反射镜的基于摄影机的六自由度标靶测量和标靶跟踪设备
CN101226057A (zh) * 2008-02-01 2008-07-23 武汉朗视软件有限公司 一种数字近景摄影测量方法
CN101387501A (zh) * 2008-10-06 2009-03-18 天津大学 超大型工件圆形截面形状与方位测量装置及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
蔡名飞: "基于高频电压注入法的永磁同步电机转子初始位置检测", 《电力电子》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110132577A (zh) * 2019-05-29 2019-08-16 河北科技大学 齿轮箱传动性能检测系统及齿轮箱传动性能检测方法
CN110132577B (zh) * 2019-05-29 2020-11-06 河北科技大学 齿轮箱传动性能检测系统及齿轮箱传动性能检测方法

Also Published As

Publication number Publication date
CN108050960B (zh) 2020-06-02

Similar Documents

Publication Publication Date Title
CN100504301C (zh) 一种天体敏感器测量基准转换方法及装置
CN103323855B (zh) 一种基线动态测量系统的精度获取方法
CN107543495A (zh) 航天器设备自动准直测量系统、准直方法与测量方法
CN106468544B (zh) 基于光电自准直仪的卫星高精度测角方法
CN104457688B (zh) 卫星上批量设备姿态角度矩阵的高精度自动化测量装置
CN108759798A (zh) 一种高精度航天器精度测量的实现方法
CN109459059B (zh) 一种星敏感器外场转换基准测定系统及方法
CN110873578B (zh) 一种基于转台传递的六面体棱镜和imu安装误差标定方法
CN103017686A (zh) 一种利用激光跟踪仪调整主镜光轴与水平轴垂直的方法
CN105823420B (zh) 一种角锥组合件回光能量中心坐标精密导出方法
CN106091966B (zh) 真空低温环境下的热变形测量方法
CN105444778B (zh) 一种基于成像几何反演的星敏感器在轨定姿误差获取方法
CN111707451B (zh) 干涉型成像光谱仪内方位元素及畸变标定方法
CN110766757B (zh) 一种具有二维指向镜的面阵相机几何成像模型标定方法
CN109087355A (zh) 基于迭代更新的单目相机位姿测量装置与方法
CN108106637A (zh) 一种分布式pos的精度检校方法和装置
CN108375383A (zh) 多相机辅助的机载分布式pos柔性基线测量方法和装置
CN105157729A (zh) 一种用于室内空间测量定位系统的调平与校准方法
Yang et al. Precision analysis of flatness measurement using laser tracker
CN108050960A (zh) 一种基于数字摄影测量技术的高精度转动测量方法
CN114167900A (zh) 一种基于无人机和差分gps的光电跟踪系统标校方法及装置
CN109470274A (zh) 一种车载光电经纬仪载车平台变形测量系统及方法
CN104330077B (zh) 一种基于两点对中模型的联合测量方法
CN104458653B (zh) 一种测量大天顶距处的大气折射值的方法和系统
CN113899324B (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