本发明属于图像信息技术处理领域，具体地涉及一种基于数字摄影测量技术的高精度转动测量方法。
背景技术：
现代天文望远镜对指向精度的要求日益提高。在射电领域，大型射电天文望远镜(如fast望远镜、德令哈13.7米口径毫米波望远镜等)的主波束宽度的典型值为10-50角秒，为了保证观测效率，通常要求望远镜指向精度达到×1角秒量级。
望远镜指向误差通常包含系统偏差部分和随机误差部分，其中系统偏差部分可以通过指向模型得以修正，修正后的残余误差则主要为随机指向误差，这部分误差主要由望远镜机架的重复定位精度决定。为了能够对望远镜机架的重复定位能力进行评估，需要对望远镜各轴系在给定转动指令条件下的实际转动角度进行精确的测量，测量精度往往需要达到1-2角秒甚至是亚角秒量级(×0.1角秒)。
在天文领域，用于高精度转动测量的设备和方法包括经纬仪、全站仪、激光跟踪仪、角摆仪等。
经纬仪和全站仪是直接的测角设备，其缺点是需要人为对准读数，引起的读数误差通常就达到了角秒量级，无法满足高精度转动测量的要求。此外，这类设备需要通过工装架设在机架上，随机架一起转动，这对小型机架的测量会造成不便。经纬仪正常工作时要求水平姿态，因此也会限制其在空间三维转动测量中的应用。
激光跟踪仪通常有较高的测距精度，但测角精度较差，极限精度通常在×1角秒量级。在使用激光跟踪仪测量转动角度时，通常将反射靶标置于待测的机架的转子上，通过测量靶标相对于激光跟踪仪的转动角度来推算机架的转动角度。由于激光跟踪仪并不位于机架的转动中心，在计算过程中有可能会将角度测量误差放大。
角摆仪的测量精度较高，可以达到1角秒甚至更高的角分辨率，但其缺点是需要配合激光干涉仪使用。此外，角摆仪只能测量平面内的转动角，无法实现空间三维任意方向轴上的转动角测量。
技术实现要素：
本发明的目的在于针对现有技术的缺陷，提供一种基于数字摄影测量技术的高精度转动测量方法。
本发明的技术方案如下：一种基于数字摄影测量技术的高精度转动测量方法包括如下步骤：一、在转动机架的转子和定子上分别粘贴反光靶标，并将定子和转子上的反光靶标数量分别为ns和nr，其中，ns和nr均为正整数；二、在初始未发生转动的条件下，对转动机架上的反光靶标进行拍照，并在拍摄完毕后，根据拍摄的照片对各个反光靶标进行解算三维坐标，记作：
靶标的坐标数值；四、根据步骤二和步骤三得到的各个反光靶标的三维坐标
优选地，在步骤四中，拟合计算在给定转动指令的条件下，转子相对于定子的相对转动角度的步骤中，利用最小二乘法找到使得拟合函数达到极小值的转动角θ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均为正整数。
在步骤一中，具体地，反光靶标的粘贴位置可以随机分布，但要求尽量均匀分布在转动圆周方向上，并且尽可能远离转动轴。
二、在初始未发生转动的条件下，对转动机架上的反光靶标进行拍照，并在拍摄完毕后，根据拍摄的照片对各个反光靶标进行解算三维坐标，记作：
在步骤二中，具体地，将得到的各个反光靶标的三维坐标作为后续转动测量的参照值。
三、将转动机架的转子转动到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μmrms，角度重复测量精度优于0.5角秒rms，与仿真结果相符。
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。