CN111457896B - 一种单像空间后方交会非迭代方法 - Google Patents
一种单像空间后方交会非迭代方法 Download PDFInfo
- Publication number
- CN111457896B CN111457896B CN202010309567.2A CN202010309567A CN111457896B CN 111457896 B CN111457896 B CN 111457896B CN 202010309567 A CN202010309567 A CN 202010309567A CN 111457896 B CN111457896 B CN 111457896B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- axis
- point
- image space
- control points
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000003384 imaging method Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims description 29
- 239000004576 sand Substances 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种单像空间后方交会非迭代方法,与传统后方交会方法相比本发明方法完全不依赖初值,可以适用于任意摄影成像条件;与目前的非迭代方法相比,本发明方法只需要求解一个二元二次方程组,恢复外方位参数时不需要额外求解方程组,在求解过程中可以避免矩阵奇异导致的不稳定因素,更加稳健高效。
Description
技术领域
本发明涉及一种单像空间后方交会非迭代方法,是对现有的现有单像空间后方交会中非迭代方法的改进,属于摄影测量技术领域。
背景技术
在摄影测量中单像空间后方交会起着非常重要的作用。针对单像空间后方交会的研究,可以分为迭代方法和非迭代方法。采用迭代方法进行后方交会,特别是倾斜角度较大时,常存在法方程病态、收敛错误和不收敛问题,这是由于没有较好的初始值造成的,无法确保最终能够收敛到正确解。非迭代方法能够求解准确初始值,但往往不能进行平差。采用非迭代方法求解外方位参数初值,再结合迭代方法对所有控制点统一进行平差,不仅能解决收敛错误问题的同时,而且能保证参数求解的精度。目前,三个控制点情况下,所用非迭代方法都需要将方程先转化为二元二次方程组,求解相关参数后再通过求解多个线性方程组恢复外方位元素,而线性方程组的求解常会遇到法方程病态问题,增加了求解过程的不稳定性。
发明内容
本发明的目的是提供一种新的求解外方位元素的非迭代方法,以解决上述现有单像空间后方交会中非迭代方法的缺陷。
本发明所述的一种单像空间后方交会的非迭代方法,采用以下技术解决方案:
S1:从所用控制点中选择三个控制点,使其在大地坐标系中构成三角形面积最大;
S2:使用步骤S1中的三个控制点A、B、C与摄影中心S构成的角锥体如图1所示,设置待求解参数k1和k2,使得向量且平面A'B'C'与平面ABC平行,根据摄影成像的几何关系构建二元二次方程组,并求解k1和k2的实数解:由摄影成像必然满足共线条件可知至少存在1组实数解,而二元二次方程组至多有4组实数解,所以k1和k2的实数解共有1至4组;
S3:根据步骤S2求解的k1和k2计算相似角锥体的相关向量,并根据几何关系恢复像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的旋转矩阵ROS和线元素POS。其中,大地坐标系O-XoYoZo为控制点采集时使用的地面坐标系,像空间坐标系S-XSYSZS以摄影中心S为原点,以焦平面水平向左为X轴方向,以焦平面水平向上为Y轴方向,以垂直焦平面指向摄影中心S的方向为Z轴方向,旋转矩阵TOS和线元素POS的含义为:若空间中任意一点在像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo中的坐标分别为PS和PO,则上角标T表示对矩阵转置;
S4:从所有控制点中选择三个控制点,使其构成的三角形与步骤S1中三角形的夹角最大,并按步骤S1-步骤S3解算相应的旋转矩阵和线元素同样共有1至4组。选择POS和的距离最小的一组记为和将和转化为对应的三个角元素和最终输出摄影成像的三个线元素和三个角元素其中,和表示大地坐标系O-XoYoZo先沿平移,再依次绕Y轴、X轴、Z轴逆时针旋转的角度为φ、ω、κ时得到像空间坐标系S-XSYSZS,和的关系如下:
其中步骤S2中所述的二元二次方程为:
其中步骤S3中根据k1和k2计算旋转矩阵ROS和线元素POS的步骤为:
1)求解辅助坐标系A-XtYtZt和像空间坐标系S-XSYSZS之间的旋转矩阵RSA和线元素PSA,若空间中任意一点在像空间坐标系S-XSYSZS和辅助坐标系A-XtYtZt中的坐标分别为PS和PA,则如图1所示,辅助坐标系A-XtYtZt以点A为坐标原点,以点A到点B方向为X轴方向,垂直平面ABC指向摄影中心S的方向为Z轴方向,在平面ABC内垂直AB方向为Y轴方向来构成右手坐标系:
3)计算像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的旋转矩阵ROS和线元素POS:ROS=ROARSA T、POS=POA-ROSPSA。
本发明的积极效果在于:
提供了一种新的单像空间后方交会的非迭代方法,与传统后方交会方法相比本发明方法完全不依赖初值,可以适用于任意摄影成像条件;与目前的非迭代方法相比,本发明方法只需要求解一个二元二次方程组,恢复外方位参数时不需要额外求解方程组,在求解过程中可以避免矩阵奇异导致的不稳定因素,更加稳健高效。
附图说明:
图1为本发明角锥体还原法示意图;
图2为本发明实施例1苏黎世市政厅(西侧视图)。
具体实施方式
通过以下实施例进一步举例描述本发明,并不以任何方式限制本发明,在不背离本发明的技术解决方案的前提下,对本发明所作的本领域普通技术人员容易实现的任何改动或改变都将落入本发明的权利要求范围之内。
实施例1
下面以计算一景实际影像的外方位为例来说明本发明的具体实施过程
如图2所示,为ISPRS发布的一组近景摄影测量数据中的一景影像,使用相机型号为Olympus C1400L,控制点分布如图2所示红色圆点标记的位置,具体数据如表1所示。
相机的内方位元素:像素大小为(W,H)=(1280像素,1024像素)、幅面大小为(w,h)=(8.245mm,6.600mm)、焦距f=8.595mm、像主点为(x0,y0)=(3.932mm,3.192mm);
表1.图2中控制点数据
本发明的具体实施过程如下:
S1:从所用控制点中选择三个控制点A(101)、B(103)、C(122),使其在大地坐标系中构成三角形面积最大;
S2:使用步骤S1中的三个控制点A(101)、B(103)、C(122)与摄影中心S构成的角锥体如图1所示,已知的控制点数据为:
求得k1和k2的实数解,如表2所示:
表2.k1和k2的实数解
编号 | k1 | k2 |
1 | 0.905 | 1.083 |
2 | 0.958 | 0.987 |
3 | 0.986 | 0.636 |
4 | 1.143 | 1.372 |
S3:对于步骤S2求解的任意一组k1和k2的实数解,计算相似角锥体的相关向量,并根据几何关系恢复像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的变换参数(旋转矩阵TOS和线元素POS),具体步骤如下:
1)求解辅助坐标系A-XtYtZt和像空间坐标系S-XSYSZS之间的旋转矩阵RSA和线元素PSA,若空间中任意一点在像空间坐标系S-XSYSZS和辅助坐标系A-XtYtZt中的坐标分别为PS和PA,则如图1所示,辅助坐标系A-XtYtZt以点A为坐标原点,以点A到点B方向为X轴方向,垂直平面ABC指向摄影中心S的方向为Z轴方向,在平面ABC内垂直AB方向为Y轴方向来构成右手坐标系。
3)计算像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的旋转矩阵ROS和线元素POS:ROS=ROARSA T、POS=POA-ROSPSA。
对于所有k1和k2的实数解,计算的TOS和POS如表3所示:
表3.使用控制点A、B、C求解的TOS和POS
S4:从所有控制点中选择三个控制点A1(102)、B1(121)、C1(123),使其构成的三角形与步骤S1中三角形的夹角最大,并按步骤S1-步骤S3解算相应的旋转矩阵和线元素同样共有1至4组,如表4所示。
结论:采用欧拉角法、四元数法使用表1的所有控制点进行整体的光束平差,是否采用本发明方法计算初始值情况下,收敛情况、迭代次数和整体平差后的外方位元素的对比如表5所示:
表5不同方法的后方交会结果
表5中可以看出,在后方交会中对所有控制点进行光束平差时,本发明方法可以提供良好的初始值,不仅能保证算法的最终收敛性,而且可以加快收敛速度。
Claims (1)
1.一种单像空间后方交会的非迭代方法,其特征在于包括以下步骤:
S1:从所用控制点中选择三个控制点,使其在大地坐标系中构成三角形面积最大;
S2:使用步骤S1中的三个控制点A、B、C与摄影中心S构成的角锥体,设置待求解参数k1和k2,使得向量且平面A'B'C'与平面ABC平行,根据摄影成像的几何关系构建二元二次方程组,并求解k1和k2的实数解:由摄影成像必然满足共线条件可知至少存在1组实数解,而二元二次方程组至多有4组实数解,所以k1和k2的实数解共有1至4组;
S3:根据步骤S2求解的k1和k2计算相似角锥体的相关向量,并根据几何关系恢复像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的旋转矩阵ROS和线元素POS;其中,大地坐标系O-XoYoZo为控制点采集时使用的地面坐标系,像空间坐标系S-XSYSZS以摄影中心S为原点,以焦平面水平向左为X轴方向,以焦平面水平向上为Y轴方向,以垂直焦平面指向摄影中心S的方向为Z轴方向,旋转矩阵TOS和线元素POS的含义为:若空间中任意一点在像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo中的坐标分别为PS和PO,则上角标T表示对矩阵转置;
S4:从所有控制点中选择三个控制点,使其构成的三角形与步骤S1中三角形的夹角最大,并按步骤S1-步骤S3解算相应的旋转矩阵和线元素同样共有1至4组;选择POS和的距离最小的一组记为和将和转化为对应的三个角元素和最终输出摄影成像的三个线元素和三个角元素其中,和表示大地坐标系O-XoYoZo先沿平移,再依次绕Y轴、X轴、Z轴逆时针旋转的角度为φ、ω、κ时得到像空间坐标系S-XSYSZS,和的关系如下:
其中步骤S2中所述的二元二次方程为:
其中步骤S3中根据k1和k2计算旋转矩阵ROS和线元素POS的步骤为:
1)求解辅助坐标系A-XtYtZt和像空间坐标系S-XSYSZS之间的旋转矩阵RSA和线元素PSA,若空间中任意一点在像空间坐标系S-XSYSZS和辅助坐标系A-XtYtZt中的坐标分别为PS和PA,则辅助坐标系A-XtYtZt以点A为坐标原点,以点A到点B方向为X轴方向,垂直平面ABC指向摄影中心S的方向为Z轴方向,在平面ABC内垂直AB方向为Y轴方向来构成右手坐标系:
3)计算像空间坐标系S-XSYSZS和大地坐标系O-XoYoZo之间的旋转矩阵ROS和线元素POS:ROS=ROARSA T、POS=POA-ROSPSA。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309567.2A CN111457896B (zh) | 2020-04-20 | 2020-04-20 | 一种单像空间后方交会非迭代方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010309567.2A CN111457896B (zh) | 2020-04-20 | 2020-04-20 | 一种单像空间后方交会非迭代方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111457896A CN111457896A (zh) | 2020-07-28 |
CN111457896B true CN111457896B (zh) | 2021-09-14 |
Family
ID=71681303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010309567.2A Expired - Fee Related CN111457896B (zh) | 2020-04-20 | 2020-04-20 | 一种单像空间后方交会非迭代方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111457896B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883325B (zh) * | 2021-01-15 | 2023-09-22 | 中国测绘科学研究院 | 一种单像空间后方交会数值计算方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281034A (zh) * | 2008-05-16 | 2008-10-08 | 南京师范大学 | 基于空间直角关系的建筑物单影像三维测量方法 |
CN102889882A (zh) * | 2012-09-03 | 2013-01-23 | 北京信息科技大学 | 一种基于光束平差的三维重建方法 |
CN104613942A (zh) * | 2015-01-30 | 2015-05-13 | 北京林业大学 | 一种像对相对外方位元素的解析方法 |
CN106248014A (zh) * | 2016-08-23 | 2016-12-21 | 中国人民解放军信息工程大学 | 一种基于单相片的三维坐标测量方法及装置 |
CN106918327A (zh) * | 2017-03-25 | 2017-07-04 | 北京林业大学 | 一种无人机航拍光束法平差摄影测量方法 |
US9857172B1 (en) * | 2017-09-25 | 2018-01-02 | Beijing Information Science And Technology University | Method for implementing high-precision orientation and evaluating orientation precision of large-scale dynamic photogrammetry system |
CN107806861A (zh) * | 2017-09-14 | 2018-03-16 | 中国人民解放军92859部队 | 一种基于本质矩阵分解的倾斜影像相对定向方法 |
CN108344398A (zh) * | 2017-01-22 | 2018-07-31 | 北京林业大学 | 一种应用3 个控制点获取地面定点多片3d 摄影外方位元素的算法 |
CN108680144A (zh) * | 2018-05-17 | 2018-10-19 | 北京林业大学 | 一种单片摄影测量标定地面点位的方法 |
CN109115218A (zh) * | 2018-07-06 | 2019-01-01 | 航天星图科技(北京)有限公司 | 一种无地面控制点的遥感影像定位方法 |
CN110160503A (zh) * | 2019-04-04 | 2019-08-23 | 武汉大学 | 一种顾及高程的无人机景观匹配定位方法 |
CN110285827A (zh) * | 2019-04-28 | 2019-09-27 | 武汉大学 | 一种距离约束的摄影测量高精度目标定位方法 |
CN110793508A (zh) * | 2019-11-11 | 2020-02-14 | 中国人民解放军战略支援部队信息工程大学 | 三维控制网中控制点的观测数据处理方法与装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8239220B2 (en) * | 2006-06-08 | 2012-08-07 | Injury Sciences Llc | Method and apparatus for obtaining photogrammetric data to estimate impact severity |
US20170091945A1 (en) * | 2015-09-25 | 2017-03-30 | Ecole Polytechnique Fed De Lausanne (Epfl) | Point and sensor estimation from images |
-
2020
- 2020-04-20 CN CN202010309567.2A patent/CN111457896B/zh not_active Expired - Fee Related
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281034A (zh) * | 2008-05-16 | 2008-10-08 | 南京师范大学 | 基于空间直角关系的建筑物单影像三维测量方法 |
CN102889882A (zh) * | 2012-09-03 | 2013-01-23 | 北京信息科技大学 | 一种基于光束平差的三维重建方法 |
CN104613942A (zh) * | 2015-01-30 | 2015-05-13 | 北京林业大学 | 一种像对相对外方位元素的解析方法 |
CN106248014A (zh) * | 2016-08-23 | 2016-12-21 | 中国人民解放军信息工程大学 | 一种基于单相片的三维坐标测量方法及装置 |
CN108344398A (zh) * | 2017-01-22 | 2018-07-31 | 北京林业大学 | 一种应用3 个控制点获取地面定点多片3d 摄影外方位元素的算法 |
CN106918327A (zh) * | 2017-03-25 | 2017-07-04 | 北京林业大学 | 一种无人机航拍光束法平差摄影测量方法 |
CN107806861A (zh) * | 2017-09-14 | 2018-03-16 | 中国人民解放军92859部队 | 一种基于本质矩阵分解的倾斜影像相对定向方法 |
US9857172B1 (en) * | 2017-09-25 | 2018-01-02 | Beijing Information Science And Technology University | Method for implementing high-precision orientation and evaluating orientation precision of large-scale dynamic photogrammetry system |
CN108680144A (zh) * | 2018-05-17 | 2018-10-19 | 北京林业大学 | 一种单片摄影测量标定地面点位的方法 |
CN109115218A (zh) * | 2018-07-06 | 2019-01-01 | 航天星图科技(北京)有限公司 | 一种无地面控制点的遥感影像定位方法 |
CN110160503A (zh) * | 2019-04-04 | 2019-08-23 | 武汉大学 | 一种顾及高程的无人机景观匹配定位方法 |
CN110285827A (zh) * | 2019-04-28 | 2019-09-27 | 武汉大学 | 一种距离约束的摄影测量高精度目标定位方法 |
CN110793508A (zh) * | 2019-11-11 | 2020-02-14 | 中国人民解放军战略支援部队信息工程大学 | 三维控制网中控制点的观测数据处理方法与装置 |
Non-Patent Citations (4)
Title |
---|
Iterative Algorithm of Space Resection Using Rodrigues Matrix;Huaien Zeng;《2010 2nd Conference on Environmental Science and Information Application Technology》;20100718;第191-194页 * |
USING PARALLEL PROJECTION MATHEMATICS TO ORIENT AN OBJECT RELATIVE TO A SINGLE IMAGE;Kyle,S;《The Photogrammetric Record》;20040331;第19卷(第105期);第38-49页 * |
单应性几何下的后方交会直接解法;黄旭 等;《遥感学报》;20160525;第431-440页 * |
单张大角度影像后方交会初值确定方法;徐振亮 等;《辽宁工程技术大学学报(自然科学版)》;20140731;第33卷(第7期);第951-954页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111457896A (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112434709B (zh) | 基于无人机实时稠密三维点云和dsm的航测方法及系统 | |
CN109272570B (zh) | 一种基于立体视觉数学模型的空间点三维坐标求解方法 | |
CN109064516B (zh) | 一种基于绝对二次曲线像的相机自标定方法 | |
CN109272574B (zh) | 基于投影变换的线阵旋转扫描相机成像模型构建方法和标定方法 | |
CN112270698B (zh) | 基于最邻近曲面的非刚性几何配准方法 | |
CN112258588A (zh) | 一种双目相机的标定方法、系统及存储介质 | |
CN108776989A (zh) | 基于稀疏slam框架的低纹理平面场景重建方法 | |
CN111457896B (zh) | 一种单像空间后方交会非迭代方法 | |
CN112150561A (zh) | 一种多相机标定方法 | |
CN113920275B (zh) | 三角网格构建方法、装置、电子设备及可读存储介质 | |
CN114463521B (zh) | 一种面向空地影像数据融合的建筑目标点云快速生成方法 | |
CN112862683A (zh) | 一种基于弹性配准和网格优化的邻接图像拼接方法 | |
CN110660099B (zh) | 基于神经网络的遥感影像处理的有理函数模型拟合方法 | |
CN113160335A (zh) | 一种基于双目视觉的模型点云及三维表面重建方法 | |
CN114202632A (zh) | 网格线性结构恢复方法、装置、电子设备及存储介质 | |
CN114170402B (zh) | 隧洞结构面提取方法、装置 | |
CN113902853A (zh) | 一种人脸三维重建方法、装置、电子设备及存储介质 | |
CN117156282A (zh) | 一种基于无人机航摄相片投影的全景图像生成方法及系统 | |
CN117372244A (zh) | 基于二维阵列表示的大场景特征立体匹配方法 | |
CN110503604B (zh) | 一种基于高精度pos的航空面阵影像实时正射拼接方法 | |
CN103295211B (zh) | 婴儿图像合成方法及装置 | |
CN114485574B (zh) | 基于卡尔曼滤波模型的三线阵影像pos辅助对地定位方法 | |
CN115797460A (zh) | 一种水下双目标定方法 | |
CN114998496A (zh) | 一种基于场景航拍图及稀疏点云的正射影像快速生成方法 | |
CN114663789A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210914 |