CN109059866B - 基于图像的行星近景摄影测量相机安装参数的测量方法 - Google Patents

基于图像的行星近景摄影测量相机安装参数的测量方法 Download PDF

Info

Publication number
CN109059866B
CN109059866B CN201810879441.1A CN201810879441A CN109059866B CN 109059866 B CN109059866 B CN 109059866B CN 201810879441 A CN201810879441 A CN 201810879441A CN 109059866 B CN109059866 B CN 109059866B
Authority
CN
China
Prior art keywords
camera
coordinate system
pointing mechanism
installation parameters
reference mirror
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
CN201810879441.1A
Other languages
English (en)
Other versions
CN109059866A (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.)
National Astronomical Observatories of CAS
Original Assignee
National Astronomical Observatories 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 National Astronomical Observatories of CAS filed Critical National Astronomical Observatories of CAS
Priority to CN201810879441.1A priority Critical patent/CN109059866B/zh
Publication of CN109059866A publication Critical patent/CN109059866A/zh
Application granted granted Critical
Publication of CN109059866B publication Critical patent/CN109059866B/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
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/02—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

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

Abstract

本发明提供了一种基于图像的行星近景摄影测量相机安装参数的测量方法，包括以下步骤：步骤1，建立行星近景摄影测量相机安装参数测量涉及的坐标系；步骤2，确定安装参数类型；步骤3，搭建安装参数的测量环境；步骤4，获取解算安装参数所需的观测数据；步骤5，解算行星近景摄影测量相机的安装参数。本发明能够对行星近景摄影测量相机的安装参数进行高精度测量，可操作性强，理论严密，并且显著提高测量效率。

Description

基于图像的行星近景摄影测量相机安装参数的测量方法
技术领域
本发明涉及一种数据处理方法，特别是一种基于图像的行星近景摄影测量相机安装参数的测量方法。
背景技术
近景摄影测量是摄影测量的一个分支，是利用对300m以内近距离目标摄影所获取的图像来确定其形态、几何位置和大小的技术。在月球探测、火星探测等行星探测任务中，近景摄影测量相机是重要的科学载荷，包括我国嫦娥三号月球探测任务、美国火星探险漫游者(MER)、火星科学实验室(MSL)等多个月球与火星探测任务均携带了全景相机、导航相机等近景摄影测量相机，我国即将实施的嫦娥五号和首次火星探测任务也将携带近景摄影测量相机。通过获取着陆区、巡视区周围场景的立体图像，基于摄影测量原理，可以恢复拍摄区域的地形，据此可开展行星表面地形地貌研究、地质构造解析、着陆器或巡视器定位、巡视探测路径规划等科学研究和工程应用。
行星近景摄影测量相机一般安装在行星探测器上可转动的指向机构上(如巡视器桅杆)，利用摄影测量原理进行地形恢复的关键是获取每幅图像的外方位元素，具体涉及相机在指向机构上的安装参数(包括旋转矩阵和偏移矩阵)、相机指向机构的转动角度、指向机构在探测器上的安装参数、以及探测器在行星表面上的实时位置与姿态等参数。其中，相机指向机构的转动角度、探测器的位置与姿态将由行星上相关传感器测量并以遥测参数的形式返回地面，而相机与指向机构的安装参数是相对固定的常数，必须在任务开始前对安装参数进行精确测量，确保图像外方位元素的解算精度及行星表面地形的恢复精度。
安装参数测量的主要目的是建立相机坐标系和行星探测器坐标系之间的转换关系，通常需要借助光学立方镜(基准镜)实现相关坐标系过渡。传统的安装参数测量方法需要利用高精度经纬仪等仪器设备对基准镜进行自准直，建立传递坐标系，实现各坐标之间的转换。然而，由于坐标系的概念比较抽象，基准镜自准直过程中随机误差较大，极大的影响了参数测量效率和精度。
发明内容
(一)要解决的技术问题
鉴于上述技术问题，本发明提供了一种基于图像的行星近景摄影测量相机安装参数的测量方法。本发明能够对行星近景摄影测量相机的安装参数进行高精度测量，可操作性强，理论严密，并且显著提高测量效率。
(二)技术方案
根据本发明的一个方面，提供了一种基于图像的行星近景摄影测量相机安装参数的测量方法，包括以下步骤：
步骤1，建立行星近景摄影测量相机安装参数测量涉及的坐标系；
步骤2，确定安装参数类型；
步骤3，搭建安装参数的测量环境；
步骤4，获取解算安装参数所需的观测数据；
步骤5，解算行星近景摄影测量相机的安装参数。
(三)有益效果
从上述技术方案可以看出，本发明至少具有以下有益效果：能够对行星近景摄影测量相机的安装参数进行高精度测量，可操作性强，理论严密，并且显著提高测量效率。
附图说明
图1为本发明实施例行星近景摄影测量相机安装参数的测量方法的流程图。
图2为本发明实施例嫦娥五号着陆器本体坐标系示意图。
图3为本发明实施例近景摄影测量相机摄测坐标系示意图。
图4A和图4B为本发明实施例嫦娥五号全景相机安装参数测量观测场景示意图。
【主要部件】
1-第一标记点；
2-第二标记点；
3-第三标记点。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
在本发明的示例性实施例中，提供了一种基于图像的行星近景摄影测量相机安装参数的测量方法。图1为本发明实施例行星近景摄影测量相机安装参数的测量方法的流程图。如图1所示，本发明的测量方法包括以下步骤：步骤1，建立行星近景摄影测量相机安装参数测量涉及的坐标系；步骤2，确定安装参数类型；步骤3，搭建安装参数的测量环境；步骤4，获取解算安装参数所需的观测数据；步骤5，解算行星近景摄影测量相机的安装参数。
在步骤1中，涉及的坐标系包括行星探测器本体坐标系、行星探测器基准镜坐标系、相机指向机构基准镜坐标系、相机指向机构偏航零位坐标系、相机指向机构俯仰零位坐标系、左、右相机摄测坐标系等与相机安装参数测量相关的坐标系。
进一步地，相机指向机构偏航零位坐标系与相机指向机构俯仰零位坐标系应保证其中一个坐标轴指向重合，另外两个坐标轴指向相互平行。
在步骤2中，安装参数类型包括两坐标系之间的偏移矩阵和旋转矩阵，其中，偏移矩阵为两坐标系原点的偏移量，旋转矩阵由一坐标系各轴指向相对于另一坐标系各轴指向的方向余弦构成。特殊的，如存在两坐标系原点重合的情况，则所述的安装参数仅包含旋转矩阵。
在步骤3中，安装参数的测量环境通过以下方式搭建：在观测场地中安置行星近景摄影测量相机、全站仪等仪器设备；在行星近景摄影测量相机前方布设若干标志牌；在相机指向机构的关节部位布设标志点。
在步骤4中，解算安装参数所需的观测数据通过以下方式获取：按照预先设定的步进，控制相机指向机构进行偏航、俯仰方向的转动，利用全站仪等设备测量转动前、后指向机构关节标志点的三维坐标；利用行星近景摄影测量相机拍摄相机指向机构转动到规定方向后的地面标志牌；利用全站仪等设备分别从两个相互垂直的方向对行星探测器基准镜、相机指向机构基准镜进行自准直；完成每个方向自准直观测后，利用全站仪等设备测量地面标志牌的坐标。
在步骤5中，安装参数的解算方法包括：利用行星近景摄影测量相机拍摄的立体像对，解算左、右相机拍摄的图像的外方位元素，构建相机摄测坐标系；利用相机指向机构的关节标志点的三维坐标，拟合相机指向机构偏航轴、俯仰轴的空间指向，解算两轴公垂线，构建相机指向机构偏航零位坐标系和相机指向机构俯仰零位坐标系；根据基准镜自准直方向的观测结果，构建行星探测器基准镜坐标系和相机指向机构基准镜坐标系，并解算坐标系之间的偏移矩阵和旋转矩阵。
进一步地，解算左、右相机拍摄的图像的外方位元素时，采用摄影测量空间后方交会的方法。
拟合相机指向机构偏航轴、俯仰轴的空间指向时，采用基于标志点三维坐标序列应位于同一个圆上的最小二乘拟合方法。
行星探测器基准镜坐标系和相机指向机构基准镜坐标系之间的转换基于七参数转换原理，所需控制点坐标为地面标志牌的观测坐标。
实施例：
本发明适用于对月球探测器近景摄影测量相机的安装参数测量，为了描述的方便，本发明以嫦娥五号任务全景相机为例对安装参数的测量方法进行具体阐述。
首先介绍一下嫦娥五号任务及全景相机的组成和工作原理。
我国探月工程三期嫦娥五号任务(CE-5)将实施月球采样并返回地球，其科学目标包括采样区现场就位探测和月球样品的分析研究。其中，采样区现场就位探测包括采样区月表形貌和地质构造调查等内容，为月壤取样对象选取提供依据，并建立现场探测数据与实验室分析数据的联系。CE-5着陆器上搭载了全景相机，由安装在相机指向机构上的两台相机组成。月面工作过程中，相机指向机构带动全景相机进行俯仰和偏航运动，左、右全景相机获取的图像构成一个立体像对，利用摄影测量原理，可恢复采样区地形，可实现采样过程图像信息支持、采样区地形地貌研究、地质构造解析等科学研究。
本发明中，安装参数是指建立全景相机坐标系和嫦娥五号着陆器本体坐标系之间的转换关系，进而求解每幅全景相机图像的外方位元素，恢复采样区地形所需的相关参数，该参数应在任务开始前进行精确测量。
如图1所示，本发明用于行星近景摄影测量相机安装参数的测量方法包含以下步骤：
步骤1，建立行星近景摄影测量相机安装参数测量涉及的坐标系。
对于嫦娥五号全景相机安装参数的测量，所涉及的坐标系主要包括：CE-5着陆器本体坐标系(如图2所示)、着陆器基准镜坐标系、相机指向机构基准镜坐标系、相机指向机构偏航零位坐标系、相机指向机构俯仰零位坐标系、左、右全景相机摄测坐标系(如图3所示)等六种坐标系。其定义分别如表1～表6所示。
表1 CE-5着陆器本体坐标系定义
Figure BDA0001754122760000051
表2着陆器基准镜坐标系定义
Figure BDA0001754122760000052
表3相机指向机构基准镜坐标系定义
Figure BDA0001754122760000053
表4相机指向机构偏航零位坐标系定义
Figure BDA0001754122760000054
Figure BDA0001754122760000061
表5相机指向机构俯仰零位坐标系定义
Figure BDA0001754122760000062
表6全景相机摄测坐标系定义
Figure BDA0001754122760000063
优选地，为了保证安装参数测量方法的可操作性强，相机指向机构偏航零位坐标系、相机指向机构俯仰零位坐标系设计为其中一个轴系重合(Y轴)，另外两个轴系相互平行的形式，参见表4、表5。
步骤2，确定安装参数类型。
针对步骤1建立的坐标系，安装参数类型包括：
右全景相机摄测坐标系相对于左全景相机摄测坐标系的偏移矩阵和旋转矩阵；
左全景相机摄测坐标系相对于相机指向机构俯仰零位坐标系的偏移矩阵和旋转矩阵；
相机指向机构俯仰零位坐标系相对于相机指向机构偏航零位坐标系的偏移矩阵和旋转矩阵；
相机指向机构偏航零位坐标系相对于相机指向机构基准镜坐标系的偏移矩阵和旋转矩阵；
相机指向机构基准镜坐标系相对于着陆器基准镜坐标系的偏移矩阵和旋转矩阵；
着陆器基准镜坐标系相对于着陆器本体坐标系的偏移矩阵和旋转矩阵。
步骤3，搭建安装参数的测量环境。
如图4A和图4B所示，安装参数的测量环境通过以下方式搭建：
在试验场地中安置CE-5全景相机和全站仪；
进行全站仪设站并建立观测坐标系；
在全景相机前方布设25个标志牌(如图4A所示)；
全景相机及相机指向机构上粘贴3个十字丝标记点，作为标志点，其中第一标记点1为相机指向机构俯仰方向的关节标志点，第二标记点2和第三标记点3均为相机指向机构偏航方向的关节标志点(如图4B所示)。
步骤4，获取解算安装参数所需的观测数据。
解算安装参数所需的观测数据通过以下方式获取：
步骤4.1，在相机指向机构偏航角为0°的情况下控制俯仰机构按步进5°由俯仰角0°旋转至俯仰角-90°，利用全站仪观测第一标记点1的坐标位置，共获取19组观测值，见表7所示。
表7第一标记点1的坐标位置观测值列表
标记点观测号 X(m) Y(m) Z(m)
1001 501.7879 500.0947 501.3551
1002 501.7881 500.0941 501.3473
1003 501.7879 500.0925 501.3397
1004 501.7887 500.0904 501.3321
1005 501.7876 500.0876 501.3249
1006 501.7879 500.0843 501.3178
1007 501.7882 500.0803 501.3112
1008 501.7885 500.0758 501.3049
1009 501.7889 500.0707 501.2988
1010 501.7887 500.065 501.2933
1011 501.7898 500.0592 501.2885
1012 501.7903 500.0525 501.284
1013 501.7897 500.0462 501.2805
1014 501.791 500.0395 501.2775
1015 501.791 500.0321 501.2749
1016 501.7919 500.0247 501.2728
1017 501.7924 500.0171 501.2715
1018 501.7932 500.0092 501.271
1019 501.7935 500.0015 501.271
相机指向机构每次转动到位后，全景相机对地面的标志牌进行拍图，共获取19对全景相机图像。
步骤4.2，在指向机构俯仰角为0°的情况下控制偏航机构按步进5°由偏航角80°旋转至偏航角-20°，利用全站仪观测第二标记点2的坐标位置，共获取21组观测值，见表8所示。
表8第二标记点2的坐标位置观测值列表
标记点观测号 X(m) Y(m) Z(m)
2001 501.9727 500.1017 501.3219
2002 501.9574 500.1032 501.3219
2003 501.9414 500.103 501.322
2004 501.9256 500.1017 501.3221
2005 501.9105 500.099 501.3224
2006 501.8955 500.0949 501.3226
2007 501.8807 500.0895 501.3229
2008 501.8662 500.0829 501.3232
2009 501.8526 500.075 501.3235
2010 501.8399 500.066 501.3239
2011 501.8285 500.056 501.3242
2012 501.8183 500.0459 501.3247
2013 501.8087 500.0342 501.3251
2014 501.7992 500.0215 501.3255
2015 501.7921 500.008 501.3259
2016 501.7863 499.9938 501.3264
2017 501.7793 499.9791 501.327
2018 501.7754 499.9639 501.3275
2019 501.7758 499.9485 501.3278
2020 501.7788 499.9326 501.328
2021 501.7733 499.9172 501.3288
保持相机指向机构俯仰角为0°，控制偏航机构按步进5°由偏航角90°旋转至偏航角-5°，利用全站仪观测第三标记点3的坐标位置，共获取20组观测值，见表9所示。
表9第三标记点3的坐标位置观测值列表
Figure BDA0001754122760000081
Figure BDA0001754122760000091
相机指向机构每次转动到位后，全景相机对地面的标志牌进行拍图，共获取41对全景相机图像。
步骤4.3，利用全站仪观测地面的标志牌的坐标位置，共获取20组观测值，见表10所示。
表10地面标志牌的坐标位置观测值列表
标志牌序号 X(m) Y(m) Z(m)
1 501.1591 509.4953 499.46
2 501.1586 503.7175 499.5813
3 502.2309 510.9434 499.4384
4 502.1854 503.8142 499.5835
5 501.6295 503.7872 499.5849
6 501.79 501.4719 500.0502
7 501.0005 504.9357 499.5627
8 501.4505 506.2992 499.5295
9 500.7354 506.2266 499.5289
10 502.144 509.4884 499.4783
11 501.7192 502.0903 499.9009
12 501.1041 510.8997 499.4304
12 500.1808 509.3579 499.461
13 502.0612 502.0853 499.9028
14 501.5232 504.9939 499.5633
15 502.1891 505.0395 499.5693
16 500.4362 507.7501 499.5018
17 501.2939 507.8236 499.505
18 502.2343 506.265 499.5369
19 502.169 502.8808 499.5921
19 502.2276 507.8282 499.5099
20 499.9868 510.8462 499.4171
21 501.2457 502.7871 499.591
22 501.736 502.8525 499.594
23 501.3853 502.0671 499.9
23 501.1591 509.4953 499.46
24 501.1586 503.7175 499.5813
24 502.2309 510.9434 499.4384
25 502.1854 503.8142 499.5835
步骤4.4，保持观测环境不变，利用全站仪对相机指向机构基准镜的一面进行自准直，记录自准直后全站仪的方位角和俯仰角读数，并观测基准镜平面内十字丝中心坐标。
换站，利用全站仪对相机指向机构基准镜的另一垂直面进行自准直，记录自准直后全站仪的方位角和俯仰角读数，并观测基准镜平面内十字丝中心坐标。
换站前、后分别利用全站仪观测至少5个标志牌的坐标位置，用于建立全站仪换站前后观测坐标系的转换关系。
步骤4.5，按照步骤4.4的方法对着陆器基准镜进行自准直并观测标志牌的坐标位置。
步骤5，解算行星近景摄影测量相机的安装参数。
解算安装参数的方法包括：
步骤5.1，构建左、右全景相机摄测坐标系。
利用步骤4.1、4.2所拍摄的图像，以及步骤4.3所测量的标志牌的坐标位置，根据摄影测量空间后方交会原理，解算相机指向机构不同俯仰角、偏航角状态下左、右相机图像的外方位元素，包括投影中心位置(摄测坐标系原点)、相机在观测坐标系中的姿态，进而构建左、右相机摄测坐标系。
步骤5.2，构建相机指向机构俯仰零位坐标系和相机指向机构偏航零位坐标系。
首先，根据步骤4.1获取的俯仰方向标志点第一标记点1的观测坐标序列，解算俯仰轴的空间指向(记为VFY_Z)，具体方法为基于最小二乘法拟合计算第一标记点1的观测坐标位置序列所构成的球面方程和平面方程，该球面与平面相交构成的圆曲线即为第一标记点1绕俯仰轴转动的轨迹，则过该圆圆心的圆平面法线即为指向机构俯仰轴的空间指向。
记第一标记点1的观测坐标位置序列所构成的球面方程表达式为：
(x-a)2+(y-b)2+(z-c)2＝R2
展开可得：
x2+y2+z2+a2+b2+c2-2ax-2by-2cz＝R2
令A＝2a,B＝2b,C＝2c,D＝a2+b2+c2-R2，则有：
Ax+By+Cz-D＝x2+y2+z2
对于每一个第一标记点1的观测坐标位置，都能列出一个如上方程，因此，对于N个观测坐标位置所列的方程组可表达为如下矩阵形式：
Figure BDA0001754122760000111
由此，可基于最小二乘原理解算A,B,C,D，进而解算球面球心位置和半径(a,b,c,R)。
同时，记第一标记点1的观测坐标位置序列所构成的平面方程表达式为：
A′x+B′y+C′z-1＝0
同理，对于N个观测坐标位置所列的方程组可表达为如下矩阵形式：
Figure BDA0001754122760000112
由此，可基于最小二乘原理解算(A′,B′,C′)，即解算平面的法向量。
至此，所求圆平面的圆心即为球心到该平面的垂足，问题转化为求平面外一点到平面的垂线问题。圆平面的圆心坐标(a‘,b’,c’)及其半径R‘为：
Figure BDA0001754122760000113
进一步，根据圆心坐标解算结果和第一标记点1任意两个位置的观测坐标，构建矢量V1和V2，则俯仰轴的空间指向为两矢量的叉积，即：
VFY_Z＝V1×V2
其次，按照同样的方法，根据步骤4.2获取的相机指向机构偏航方位标志点第二标记点2和第三标记点3的观测坐标序列，解算偏航轴的空间指向(记为VPH_X)。
最后，解算相机指向机构俯仰轴和偏航轴的公垂线方向，作为相机指向机构俯仰零位坐标系和相机指向机构偏航零位坐标系的公共轴系(记为VPH_Y)，公垂线分别与俯仰轴、偏航轴的交点作为指向机构俯仰零位坐标系、偏航零位坐标系的原点，按右手系构建指向机构俯仰零位坐标系和偏航零位坐标系。
VFY_Y＝VPH_Y＝VFY_Z×VPH_X
VPH_Z＝VPH_X×VPH_Y
VFY_X＝VFY_Y×VFY_Z
步骤5.3，构建相机指向机构基准镜坐标系和着陆器基准镜坐标系。
根据步骤4.4获取的对相机指向机构基准镜两个相邻、相互垂直面自准直后的方位、俯仰角度，以及十字丝中心的观测坐标，解算相机指向机构基准镜相邻两面的法线。
根据步骤4.4观测的地面标志牌坐标，解算换站前、后全站仪观测坐标系的转换七参数，实现两条法线的坐标系统一，进而解算两法线的交点，作为相机指向机构基准镜坐标系的原点，按右手系构建相机指向机构基准镜坐标系。
根据步骤4.5的观测结果，按照同样的方法构建着陆器基准镜坐标系。
步骤5.4，根据探测器设计值构建着陆器本体坐标系。
步骤5.5，解算各类安装参数。
其中，右全景相机摄测坐标系相对于左全景相机摄测坐标系的偏移矩阵和旋转矩阵，应根据步骤5.1中构建的相机指向机构不同俯仰、方位角状态下左、右全景相机摄测坐标系进行解算，解算结果取各状态的加权平均值。
其中，左全景相机摄测坐标系相对于相机指向机构俯仰零位坐标系的偏移矩阵和旋转矩阵，应根据步骤5.1中构建的指向机构俯仰角0°、方位角0°状态下左全景相机摄测坐标系，与步骤5.2中构建的相机指向机构俯仰机构零位坐标系进行解算。
其中，相机指向机构俯仰机构零位坐标系相对于相机指向机构偏航机构零位坐标系的偏移矩阵和旋转矩阵，应根据步骤5.2中构建的相机指向机构俯仰机构、偏航机构零位坐标系进行解算。
其中，相机指向机构偏航机构零位坐标系相对于相机指向机构基准镜坐标系的偏移矩阵和旋转矩阵，应根据步骤5.2中构建的相机指向机构偏航机构零位坐标系，与步骤5.3中构建的相机指向机构基准镜坐标系进行解算。
其中，相机指向机构基准镜坐标系相对于着陆器基准镜坐标系的偏移矩阵和旋转矩阵，应根据步骤5.3中构建的相机指向机构基准镜坐标系与着陆器基准镜坐标系进行解算。
其中，着陆器基准镜坐标系相对于着陆器本体坐标系的偏移矩阵和旋转矩阵，应根据步骤5.3中构建的着陆器基准镜坐标系，与步骤5.4中构建的着陆器本体坐标系进行解算。
至此，可以得到嫦娥五号全景相机的安装参数为：
右相机到左相机偏移矩阵＝(262.81mm,0mm,0mm)
Figure BDA0001754122760000131
左相机到相机指向机构俯仰零位坐标系偏移矩阵＝(-9.01mm,-137.33mm,28.40mm)
Figure BDA0001754122760000132
相机指向机构俯仰零位坐标系到偏航零位坐标系偏移矩阵＝(0.00mm,0.00mm,84.21mm)
Figure BDA0001754122760000133
相机指向机构偏航零位坐标系到相机指向机构基准镜偏移矩阵＝(99.83mm,0.05mm,184.74mm)
Figure BDA0001754122760000134
相机指向机构基准镜到着陆器基准镜偏移矩阵＝(0mm,0mm,0mm)
Figure BDA0001754122760000141
着陆器基准镜到着陆器本体偏移矩阵＝(0mm,0mm,0mm)
Figure BDA0001754122760000142
至此，已经结合附图对本实施例进行了详细描述。依据以上描述，本领域技术人员应当对本发明基于图像的行星近景摄影测量相机安装参数的测量方法有了清楚的认识。本发明提出了一种几何光学测量与立体影像相结合的行星近景摄影测量相机安装参数的测量方法，结合行星近景摄影测量相机的特点解算高精度的相机安装参数，解算精度满足后续图像外方位元素解算和地形恢复的需求。
需要说明的是，在附图或说明书正文中，未绘示或描述的实现方式，均为所属技术领域中普通技术人员所知的形式，并未进行详细说明。此外，上述对各元件和方法的定义并不仅限于实施例中提到的各种具体结构、形状或方式，本领域普通技术人员可对其进行简单地更改或替换。
还需要说明的是，除非特别描述或必须依序发生的步骤，上述步骤的顺序并无限制于以上所列，且可根据所需设计而变化或重新安排。并且上述实施例可基于设计及可靠度的考虑，彼此混合搭配使用或与其他实施例混合搭配使用，即不同实施例中的技术特征可以自由组合形成更多的实施例。
在此提供的算法和显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述，构造这类系统所要求的结构是显而易见的。此外，本发明也不针对任何特定编程语言。应当明白，可以利用各种编程语言实现在此描述的本发明的内容，并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于图像的行星近景摄影测量相机安装参数的测量方法，其特征在于，包括以下步骤：
步骤1，建立行星近景摄影测量相机安装参数测量涉及的坐标系，包括：行星探测器本体坐标系、行星探测器基准镜坐标系、相机指向机构基准镜坐标系、相机指向机构偏航零位坐标系、相机指向机构俯仰零位坐标系、左、右相机摄测坐标系；
步骤2，确定安装参数类型；
步骤3，搭建安装参数的测量环境，包括在相机指向机构的关节部位布设标志点；
步骤4，获取解算安装参数所需的观测数据，包括按照预先设定的步进，控制相机指向机构进行偏航、俯仰方向的转动，利用全站仪测量转动前、后相机指向机构关节标志点的三维坐标；
步骤5，解算行星近景摄影测量相机的安装参数。
2.根据权利要求1所述的测量方法，其特征在于，相机指向机构偏航零位坐标系与相机指向机构俯仰零位坐标系的其中一个坐标轴指向重合，另外两个坐标轴指向相互平行。
3.根据权利要求1所述的测量方法，其特征在于，所述安装参数类型包括两坐标系之间的偏移矩阵和旋转矩阵，其中，偏移矩阵为两坐标系原点的偏移量，旋转矩阵由一坐标系各轴指向相对于另一坐标系各轴指向的方向余弦构成。
4.根据权利要求1所述的测量方法，其特征在于，步骤3包括以下子步骤：
在观测场地中安置行星近景摄影测量相机、全站仪；
在行星近景摄影测量相机前方的地面上布设若干标志牌。
5.根据权利要求1所述的测量方法，其特征在于，步骤4包括以下子步骤：
利用行星近景摄影测量相机拍摄相机指向机构转动到规定方向后的标志牌；
利用全站仪分别从两个相互垂直的方向对行星探测器基准镜、相机指向机构基准镜进行自准直；
完成每个方向自准直观测后，利用全站仪测量标志牌的坐标。
6.根据权利要求1所述的测量方法，其特征在于，步骤5包括以下子步骤：
利用行星近景摄影测量相机拍摄的立体像对，解算左、右相机拍摄的图像的外方位元素，构建左、右相机摄测坐标系；
利用相机指向机构的关节标志点的三维坐标，拟合相机指向机构偏航轴、俯仰轴的空间指向，解算两轴公垂线，构建相机指向机构偏航零位坐标系和相机指向机构俯仰零位坐标系；
根据基准镜自准直方向的观测结果，构建行星探测器基准镜坐标系和相机指向机构基准镜坐标系；
解算坐标系之间的偏移矩阵和旋转矩阵。
7.根据权利要求6所述的测量方法，其特征在于，解算左、右相机拍摄的图像的外方位元素时，采用摄影测量空间后方交会的方法。
8.根据权利要求6所述的测量方法，其特征在于，拟合相机指向机构偏航轴、俯仰轴的空间指向时，采用基于标志点三维坐标序列应位于同一个圆上的最小二乘拟合方法。
9.根据权利要求6所述的测量方法，其特征在于，行星探测器基准镜坐标系和相机指向机构基准镜坐标系之间的转换基于七参数转换原理，所需控制点坐标为标志牌的观测坐标。
CN201810879441.1A 2018-08-03 2018-08-03 基于图像的行星近景摄影测量相机安装参数的测量方法 Active CN109059866B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810879441.1A CN109059866B (zh) 2018-08-03 2018-08-03 基于图像的行星近景摄影测量相机安装参数的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810879441.1A CN109059866B (zh) 2018-08-03 2018-08-03 基于图像的行星近景摄影测量相机安装参数的测量方法

Publications (2)

Publication Number Publication Date
CN109059866A CN109059866A (zh) 2018-12-21
CN109059866B true CN109059866B (zh) 2020-12-18

Family

ID=64833139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810879441.1A Active CN109059866B (zh) 2018-08-03 2018-08-03 基于图像的行星近景摄影测量相机安装参数的测量方法

Country Status (1)

Country Link
CN (1) CN109059866B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110095659B (zh) * 2019-04-30 2021-06-22 北京卫星环境工程研究所 深空探测巡视器通讯天线指向精度动态测试方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007021781A2 (en) * 2005-08-10 2007-02-22 Space Adventures, Ltd. Method, apparatus, and system for private lunar exploration
US8159357B1 (en) * 2009-03-30 2012-04-17 Philip Onni Jarvinen Means to prospect for water ice on heavenly bodies
CN104111061A (zh) * 2014-06-30 2014-10-22 中国电力工程顾问集团中南电力设计院 一种基础资料缺乏地区测量控制点成果获取方法
CN105374067A (zh) * 2015-10-10 2016-03-02 长安大学 一种基于pal相机的三维重建方法及其重建系统
CN105509733A (zh) * 2015-11-30 2016-04-20 上海宇航系统工程研究所 非合作空间圆目标的相对位姿测量方法
CN105678783A (zh) * 2016-01-25 2016-06-15 西安科技大学 折反射全景相机与激光雷达数据融合标定方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007021781A2 (en) * 2005-08-10 2007-02-22 Space Adventures, Ltd. Method, apparatus, and system for private lunar exploration
US8159357B1 (en) * 2009-03-30 2012-04-17 Philip Onni Jarvinen Means to prospect for water ice on heavenly bodies
CN104111061A (zh) * 2014-06-30 2014-10-22 中国电力工程顾问集团中南电力设计院 一种基础资料缺乏地区测量控制点成果获取方法
CN105374067A (zh) * 2015-10-10 2016-03-02 长安大学 一种基于pal相机的三维重建方法及其重建系统
CN105509733A (zh) * 2015-11-30 2016-04-20 上海宇航系统工程研究所 非合作空间圆目标的相对位姿测量方法
CN105678783A (zh) * 2016-01-25 2016-06-15 西安科技大学 折反射全景相机与激光雷达数据融合标定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
嫦娥3号巡视探测器桅杆的转动角度误差检校;张烁;《测绘科学》;20140930;第39卷(第9期);第1-6页 *
嫦娥三号月面巡视探测器导航定位技术研究;马友青;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20150615;第15-25页 *
马友青.嫦娥三号月面巡视探测器导航定位技术研究.《中国博士学位论文全文数据库 工程科技Ⅱ辑》.2015, *

Also Published As

Publication number Publication date
CN109059866A (zh) 2018-12-21

Similar Documents

Publication Publication Date Title
CN110296691B (zh) 融合imu标定的双目立体视觉测量方法与系统
CN111156998B (zh) 一种基于rgb-d相机与imu信息融合的移动机器人定位方法
Johnson et al. Precise image-based motion estimation for autonomous small body exploration
CN110095659B (zh) 深空探测巡视器通讯天线指向精度动态测试方法
CN107223275A (zh) 多路传感数据融合的方法和系统
CN109146958B (zh) 一种基于二维图像的交通标志空间位置测量方法
CN112184786B (zh) 一种基于合成视觉的目标定位方法
CN106525007A (zh) 分布交互式测绘通用机器人
CN106595608A (zh) 分布交互通用测绘仪
Lo et al. The direct georeferencing application and performance analysis of uav helicopter in gcp-free area
CN112461204A (zh) 卫星对动态飞行目标多视角成像联合计算航行高度的方法
CN106595600B (zh) 立体测绘卫星姿态测定系统低频误差补偿方法及系统
Burkard et al. User-aided global registration method using geospatial 3D data for large-scale mobile outdoor augmented reality
CN113607188B (zh) 基于经纬仪叉丝成像的多视场星敏感器的标定系统及方法
CN109506662B (zh) 一种小天体着陆初始对准方法、其相对导航基准确定方法及装置
CN109059866B (zh) 基于图像的行星近景摄影测量相机安装参数的测量方法
JPH11514434A (ja) イメージデータを用いてカメラの位置及び配向を決定する方法と装置
CN112785686A (zh) 一种基于大数据的林区地图构建方法及可读存储介质
RU2723199C1 (ru) Способ и система определения ориентации космического аппарата в пространстве с автономной коррекцией эффекта аберрации света
CN103033181B (zh) 确定第三目标的光学目标定位器
Wu Photogrammetry: 3-D from imagery
Johnston Ground object geo-location using UAV video camera
Wang et al. A mobile stereo vision system with variable baseline distance for three-dimensional coordinate measurement in large FOV
CN112712559B (zh) 基于NED坐标系向量旋转的SfM点云校正方法
CN113759365B (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