CN106023146B - 用于摄影测量中的场相关单边自标定光束平差方法 - Google Patents
用于摄影测量中的场相关单边自标定光束平差方法 Download PDFInfo
- Publication number
- CN106023146B CN106023146B CN201610298071.3A CN201610298071A CN106023146B CN 106023146 B CN106023146 B CN 106023146B CN 201610298071 A CN201610298071 A CN 201610298071A CN 106023146 B CN106023146 B CN 106023146B
- Authority
- CN
- China
- Prior art keywords
- field
- parameter
- calibration
- model
- self
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000003384 imaging method Methods 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 230000001419 dependent effect Effects 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000003044 adaptive effect Effects 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 11
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- YTAHJIFKAKIKAV-XNMGPUDCSA-N [(1R)-3-morpholin-4-yl-1-phenylpropyl] N-[(3S)-2-oxo-5-phenyl-1,3-dihydro-1,4-benzodiazepin-3-yl]carbamate Chemical compound O=C1[C@H](N=C(C2=C(N1)C=CC=C2)C1=CC=CC=C1)NC(O[C@H](CCN1CCOCC1)C1=CC=CC=C1)=O YTAHJIFKAKIKAV-XNMGPUDCSA-N 0.000 claims description 2
- 230000000903 blocking effect Effects 0.000 claims 1
- 230000000295 complement effect Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 22
- 238000002474 experimental method Methods 0.000 description 7
- 241000350481 Pterogyne nitens Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- 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)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种用于摄影测量中的场相关单边自标定光束平差方法,包括步骤:a)建立线性的针孔成像模型,并附加场相关的非线性畸变模型;b)确定场相关的单边自标定光束平差理论模型;c)建立场相关单边自标定成像模型的误差方程;d)求取误差方程中场相关成像模型对于外方位参数的偏导数;e)求取误差方程中场相关模型对于空间坐标的偏导数;f)求取误差方程中场相关模型对于单边内参数的偏导数;g)对误差方程雅可比矩阵中的各项进行自适应的比例调整;h)通过分块方式快速计算法方程中的各项;i)求解单边自标定中的各项参数;j)参数比例逆调整,消除比例变化;k)计算中误差,并进行单边自标定各项参数的不确定度估计。
Description
技术领域
本发明涉及一种场相关单边自标定光束平差方法,尤其涉及大尺寸多站位的静态摄影测量系的场相关单边自标定光束平差方法。
背景技术
摄影测量使用测量网络中多站位、多角度拍摄的被测物图片,利用光线束交会定位空间点坐标,大量冗余信息的使用降低了测量的***误差。相比于测量范围,摄影测量的普遍坐标测量相对误差在[1/100000,1/200000],长度测量相对误差在[1/50000,1/100000][1]。摄影测量除了具有非接触、成本低、高效率的优点外,在精度方面又不亚于其他测量设备。
摄影测量的高精度得益于两个关键因素:回光反射目标点的使用以及像面目标点的高精度定位;具有自标定能力的光束平差技术。自标定光束平差是指在测量过程中完成相机内参数标定的光束平差技术,自标定的意义在于:一方面降低了实验室标定的要求,另一方面使得相机内参数的标定结果与拍摄环境(温度、压力、光照、材质等)相适应,能够在复杂的测量环境中得到最精确的测量结果。
当人们最开始研究成像***时,已经注意到象差同放大率之间的关系,但是很长时间并未在摄影测量相机模型中应用。直至上世纪70年代,Duane Brown提出了简便的数学描述和标定方法,将相关模型引入摄影测量,Fryer和Fraser也认为物距对于畸变参数有不可忽略的影响。但是这种模型在自标定光束平差中的应用受到限制。
国内对于摄影测量的研究和开发没有深入探讨目标点空间距离对于其成像带来的影响,光束平差模型和算法也不考虑空间坐标、外方位参数与畸变参数之间的相关性。虽然平差结果给出较为满意的空间坐标误差估计,但事实上,相关研究指出,空间误差远大于平差结论。
发明内容
为了解决上述技术问题,本发明提供一种摄影测量相机的场相关单边自标定光束平差方法,包括步骤:a)建立线性的针孔成像模型,并附加场相关的非线性畸变模型,用于描述物空间点、外方位参数、内方位参数之间的数学关系;b)确定场相关单边自标定光束平差理论模型;c)建立场相关单边自标定光束平差理论模型的误差方程;d)求取误差方程中场相关成像模型对于外方位参数的偏导数;e)求取误差方程中场相关成像模型对于空间坐标的偏导数;f)求取误差方程中场相关成像模型对于单边内参数的偏导数;g)对误差方程雅可比矩阵中的各项进行自适应的比例调整,以降低法方程的病态;h)通过分块方式快速计算法方程中的各项;i)通过最小二乘技术迭代求解单边自标定光束平差理论模型中的各项参数;j)参数比例逆调整,消除步骤g)中的人为定义的比例变化;k)计算中误差,并进行单边自标定光束平差理论模型各项参数的不确定度估计。
优选地,所述步骤a)中建立线性成像模型的步骤为:
a1)确定物空间坐标系至相机空间坐标系之间的刚体旋转矩阵:
其中,Az,El,Ro为空间坐标系的三个旋转角;
a2)通过旋转矩阵和平移向量描述物空间至相机空间之间的刚体变换:
a3)在相机坐标系下,经小孔成像,将相机空间中的点投影至像面,得到对应的像点坐标:
优选地,所述步骤a)中建立非线性畸变模型的步骤为:
a4)计算位于某一坐标(x,y)处的像点畸变量:
其中,相机的主点偏移为xp和yp,对应此空间点的畸变参数为K1ss’,K2ss’,K3ss’,P1,P2,B1,B2。
a5)获取空间点在该站位相机像面上最终的像点坐标:
x=xl+xp-Δx
优选地,所述步骤b)中确定场相关单边自标定光束平差理论模型的步骤为:
b1)标定两个距离s1和s2上径向畸变参数;
b2)推导任意其他物距s′上的畸变参数:
其中,s是成像***对焦距离,以及是两个距离上径向畸变参数的标定结果,和分别是两个距离对应高斯成像公式的像距;s′为物空间点到相机xoy平面的距离。
b3)将两个距离s1和s2中任一距离的径向畸变参数作为已知量,对另一端的径向畸变参数以及其他参数进行平差。
优选地,所述步骤c)中建立场相关单边自标定成像模型的误差方程为:
其中,用(xij,yij)表示第i张相片对于空间第j点的成像坐标;Xij 0是所有与成像坐标(xij,yij)相关的参数,(vxij,vyij)是残余误差,Aij,Bij和Cij分别是步骤a5中获取的空间点在该站位相机像面上最终的像点坐标对于第i图片的外方位参数、第j空间点坐标以及相机成像参数的雅可比矩阵。
优选地,所述步骤d)至步骤f)中求取外方位参数、空间坐标以及内参数的偏导数的方法为:利用描述的场相关参数与畸变参数之间的直接关系,通过复合函数求导方法,利用中间量减少偏导数求解的复杂性和计算量。
优选地,所述步骤g中对误差方程雅可比矩阵中的各项进行自适应的比例调整的步骤为:
g1)在每次平差迭代过程开始,首先统计误差方程中对应单边内方位参数雅可比矩阵各列的数量级Mj及调整系数kj,计算方法如下:
g2)将对应单边内方位参数的雅可比矩阵的每一列计算结果乘以相应列的比例系数,得到调整后的雅可比矩阵为:
优选地,所述步骤h)通过分块方式快速计算法方程中的各项的步骤为:
h1)计算Aij,Bij和Cij,分别是步骤c中获取的成像模型对于第i图片的外方位参数、第j空间点坐标以及相机成像单边内参数的雅可比矩阵。
h2)将求解得到的Aij,Bij和Cij,按照如下公式对法方程中的各项进行有规则的分块描述:从而在对图像、目标点的索引过程中逐步得到法方程中的矩阵:
h3)获取最终的法方程的分块模型:
优选地,所述步骤i)通过最小二乘技术迭代求解单边自标定中的各项参数的步骤为:
i1)利用N22项的分块对角性,将法方程进一步改写为:
i2)利用上式第一、三等式:
求解外方位参数和单边内方位参数的增量:
i3)将i2)的计算结果带入i1)中求解空间坐标增量其中利用N22的分块对角性进行简化:
优选地,所述步骤j)参数比例逆调整,消除步骤g)中的人为定义的比例变化的方法为:确定未知量δ′与需要求解的δ之间的关系为:
优选地,所述步骤k)计算中误差,并进行单边自标定各项参数的不确定度估计的步骤为:
k1)计算单位权中误差:
其中,V是所有观测坐标的残余误差;n是观测坐标数目;u是未知变量数目。
k2)通过协方差矩阵的误差传播定律,计算外方位参数、内参数以及空间坐标的协方差矩阵:
其中,N为系数矩阵。
总结上述描述,本发明的用于摄影测量中的场相关单边自标定光束平差方法适用于摄影测量的大型稀疏矩阵的分块运算技术已经非常成熟,本发明的相关研究也展示出了场相关的畸变模型对于测量精度的有效提高。如何运用摄影测量传统光束平差中的分块运算技术实现场相关模型的扩展并取得更高精度的测量结果,在国内是一个空白。
本发明以准确的数学描述为基础,从算法上实现场相关的自标定光束平差,解决了以下技术问题:
1、分析场相关的相机成像模型,确定同畸变量相关的相机参数;
2、建立场相关的单边自标定误差方程,对误差方程、法方程中相关项进行准确描述;
3、分块光束平差算法、不确定度估计及自适应的数值比例调节。
附图说明
为了更清楚的说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明的用于摄影测量中的场相关单边自标定光束平差方法的流程图;
图2a为本发明的用于摄影测量中的场相关单边自标定光束平差方法的实验中测量场;
图2b为实验中测量场中的定向靶。
具体实施方式
通过参考示范性实施例,本发明的目的和功能以及用于实现这些目的和功能的方法将得以阐明。然而,本发明并不受限于以下所公开的示范性实施例,可以通过不同形式来对其加以实现。说明书的实质仅仅是帮助相关领域技术人员综合理解本发明的具体细节。
如图1所示,为本发明的用于摄影测量中的场相关单边自标定光束平差方法的流程图,步骤为:
101)建立线性的针孔成像模型,并附加场相关的非线性畸变模型,来描述物空间点、外方位参数、内方位参数之间的数学关系;
空间点通过坐标系转换、线性成像以及非线性畸变,投影成为一个像面点,其中,建立线性成像模型的步骤为:
a1)在世界坐标系下,令空间点坐标为[X Y Z]T,相机位置和指向通过[X0 Y0 Z0AzEl Ro]T描述,其中前三项表示相机位置坐标,后三项表示相机的方位、俯仰及旋转角,可以通过下式计算旋转矩阵:
其中,Az,El,Ro为空间坐标系的三个旋转角;
a2)通过旋转矩阵和平移向量描述物空间至相机空间之间的刚体变换:
a3)在相机坐标系下,经过成像***投影至像面上,得到线性投影坐标:
其中,建立非线性畸变模型的步骤为:
a4)计算位于某一坐标(x,y)处的像点畸变量:
其中,相机的主点偏移为xp和yp,对应此空间点的畸变参数为K1ss’,K2ss’,K3ss’,P1,P2,B1,B2。
a5)于是空间点在此站位相机像面上最终的像点坐标为:
x=xl+xp-Δx
102)确定场相关单边自标定光束平差理论模型;根据本发明的一个实施例,确定场相关单边自标定光束平差理论模型的步骤为:
b1)在回光反射共面直线阵列标定场内标定两个距离s1和s2上的径向畸变参数:
b2)推导任意其他物距s′上的畸变参数:
其中,s是成像***对焦距离,和分别是两个距离对应高斯成像公式的像距;s′为物空间点到相机xoy平面的距离。
b3)将s1和s2两端的径向畸变参数区别对待。本发明定义,固定一端的径向畸变参数而平差另一端参数的方法叫单边自标定光束平差。将s1端的径向畸变参数作为已知量,对s2端的径向畸变参数 以及其他参数进行平差。
103)建立场相关单边自标定光束平差理论模型的误差方程;
用(xij,yij)表示第i张相片对于空间第j点的成像坐标,可以得到如下面公式描述的误差方程。其中,Xij 0是所有与成像坐标(xij,yij)相关的参数,(vxij,vyij)是残余误差,Aij,Bij和Cij分别是成像模型(4)对于第i图片的外方位参数、第j空间点坐标以及相机成像参数的雅可比矩阵。
将公式(5)描述的畸变参数模型应用至自标定光束平差的关键在于将s1和s2两端的径向畸变参数区别对待和处理。本发明中定义,固定两端径向畸变参数的方法叫做非自标定光束平差;固定一端的径向畸变参数而平差另一端参数的方法叫单边自标定光束平差;两端的径向畸变参数都参与平差的方法叫做双边自标定光束平差。本发明中只介绍单边自标定光束平差。
由上可知,径向畸变参数同空间坐标、外方位参数紧密相关,下面的分析均将s1端的径向畸变参数作为已知量,对s2端的径向畸变参数 以及其他参数进行平差。
104)求取场相关成像模型的误差方程对于外方位参数的偏导数;
详细推导成像模型对于X0偏导数,其他项偏导数可以按照类似方法求解。各外方位参数及空间坐标通过s′与畸变参数产生联系,先求解畸变参数对于s′的偏导数:
可知:
通过公式(2)的描述,可以求解xl和yl对于X0的偏导数为:
于是,结合公式(4),观测值对于外方位参数中X0项的偏导数为:
同理,可以得到观测值对于其他外方位参数的偏导数:
线性项xl和yl对于各角度的偏导数描述比较复杂,可以参照相关文献。这里仅就Δx和Δy对于各角度的偏导数进行分析。
105)求取误差方程中场相关成像模型对于空间坐标的偏导数;
由于:
因此:
106)求取误差方程中场相关成像模型对于单边内参数的偏导数;
由上可知,参与平差的内参数是xp,yp,f,P1,P2,B1,B2以及S2单边径向畸变参数
107)对误差方程雅可比矩阵中的各项进行自适应的比例调整,以降低法方程的病态;
待求未知量在数量级上相差巨大,尤其是内参数,例如K1、P1、P2、B1和B2项一般具有10-5,K2项具有10-7,而K3具有10-11,这种数量级上的巨大差异带来了误差方程中矩阵A的病态,由于数值计算中普遍存在的机器精度、大数吃小数以及大数除小数等现象,会严重影响法方程的计算结果,往往会出现矩阵缺秩等问题。为了统一误差方程中各项的数量级,这里设计了自适应的比例调整技术。
在每次平差迭代过程开始,首先统计误差方程中各列的数量级Mj及调整系数kj,计算方法如下:
然后将每一列的计算结果乘以相应列的比例系数,得到调整后的误差方程矩阵:
108)通过分块方式快速计算法方程中的各项;
h1)对于第i图片的外方位参数、第j空间点坐标以及相机成像单边内参数分别求解偏导数:
h2)设一个具有m张图片的测量网络对于n个空间点进行摄影测量,那么误差方程组为:
对应的法方程为:
(ATA)δ=ATl (26)
鉴于A的稀疏性,ATA及ATl可以通过对于图片和目标点的索引进行有规则的分块描述:
h3)于是,对于(26)的求解演变成对于下式的求解:
109)通过最小二乘技术迭代求解单边自标定光束平差理论模型中的各项参数;
i1)充分利用N22项的分块对角性,将法方程进一步描述为:
i2)利用上式第一、三等式:
求解外方位参数和内参数的增量:
i3)将i2)的计算结果带入1)中求解空间坐标增量其中利用N22的分块对角性进行简化:
110)参数比例逆调整,消除步骤g)中的人为定义的比例变化;
未知量δ′与需要求解的δ之间的关系为:
111)计算中误差,并进行单边自标定光束平差理论模型各项参数的不确定度估计。
公式(32)、(33)和(34)计算的外方位参数、单边内参数以及空间坐标的不确定度可以通过其误差的协方差矩阵进行描述,误差协方差矩阵通过误差传播定律给出。单位权中误差为:
其中,V是所有观测坐标的残余误差;n是观测坐标数目;u是未知变量数目。
各未知量的协方差矩阵为:
其中,N是h3)描述的系数矩阵。
以下为实验验证本发明的用于摄影测量中的场相关单边自标定光束平差方法的技术效果。
实验为用一台AVT工业相机,具有1600万像素的全画幅传感器,搭配35mm尼康定焦镜头及商用闪光灯。数据通过千兆网络传输,并直接进入编写的单边场相关自标定光束平差软件算法进行处理,图2a为本发明的用于摄影测量中的场相关单边自标定光束平差方法的实验中测量场,如图所示,由编码点201a和普通点201b、定向靶203以及背景202组成测量场;其中如图2b所示,定向靶203的两端布置有编码点或/和普通点。
测量之前通过特殊标定场和特殊标定方法标定相机在4.489米处的各项畸变参数(主要是径向畸变参数),在单边场相关自标定光束平差的过程中计算自标定另一距离上的各项畸变参数,结果如下表所示。
表1场相关单边自标定光束平差相机内参数
距离1 | 距离2 | |
S(mm) | 4489.4 | 1909.5 |
xp(mm) | -0.1267 | -0.2057 |
yp(mm) | 0.0278 | 0.1973 |
f(mm) | 36.0997 | |
K1 | 6.6597E-05 | 6.5463E-05 |
K2 | -5.2602E-08 | -3.0186E-08 |
K3 | -1.1628E-11 | -4.6766E-11 |
P1 | 8.1816E-06 | 6.9022E-06 |
P2 | -1.2852E-05 | -3.4025E-06 |
B1 | -7.9364E-03 | 1.1967E-04 |
B2 | -1.7230E-04 | -2.8629E-05 |
可以看出,自标定结果体现出了不同距离上径向畸变参数的差别。空间点坐标不确定度估计的统计结果如表2所示。
表2场相关模型光束平差统计结果(单位:mm)
为了实现场相关及场无关自标定模型测量精度的比较,设计了相关实验。在测量场中不同位置放置等长度基准尺,用两种模型分别进行摄影测量,最后统计长度测量的不确定度,统计结果如下表所示。
表3不同光束平差模型长度测量不确定度的统计结果
场相关模型 | 场无关模型 | |
长度不确定度(mm) | 0.012 | 0.020 |
长度测量相对误差 | 1/79000 | 1/47500 |
本发明详细介绍了场相关的径向畸变模型在摄影测量中的应用方法。推导了场相关的单边自标定光束平差中的误差方程;解决了如何利用误差方程的稀疏性分块求解法方程及未知数改正数,在平差结束给出未知参数的不确定度估计;通过自适应的比例调节消除平差运算中的病态问题;最后通过实验验证了场相关单边自标定技术的可行性,并验证了这种模型对于摄影测量长度测量精度的提高。
以上只是本发明较佳的实例,并非来限制本发明实施范围,故凡依本发明申请专利范围所述的构造、特征及原理所做的等效变化或修饰,均应包括于本发明专利申请范围内。
Claims (7)
1.一种用于摄影测量中的场相关单边自标定光束平差方法,所述方法包括以下步骤:
a)建立线性的针孔成像模型,并附加场相关的非线性畸变模型,用于描述物空间点、外方位参数、内方位参数之间的数学关系;
b)确定场相关单边自标定光束平差理论模型;
c)建立场相关单边自标定光束平差理论模型的误差方程;
d)求取误差方程中场相关成像模型对于外方位参数的偏导数;
e)求取误差方程中场相关成像模型对于空间坐标的偏导数;
f)求取误差方程中场相关成像模型对于单边内参数的偏导数;
g)对误差方程雅可比矩阵中的各项进行自适应的比例调整,以降低法方程的病态;
h)通过分块方式快速计算法方程中的各项;
i)通过最小二乘方法迭代求解单边自标定光束平差理论模型中的各项参数;
j)参数比例逆调整,消除步骤g)中的人为定义的比例变化;
k)计算中误差,并进行单边自标定光束平差理论模型各项参数的不确定度估计;
所述步骤b)中确定场相关单边自标定光束平差理论模型的步骤为:
b1)标定两个距离s1和s2上径向畸变参数;
b2)推导任意其他物距s′上的畸变参数:
其中,s是成像***对焦距离,以及是两个距离上径向畸变参数的标定结果;
和分别是两个距离对应高斯成像公式的像距;
s′为物空间点到相机xoy平面的距离;
αs′为中间计算量;
b3)将两个距离s1和s2中任一距离的径向畸变参数作为已知量,对另一端的径向畸变参数以及其他参数进行平差;
其中,所述步骤c)中建立场相关单边自标定成像模型的误差方程为:
其中,用(xij,yij)表示第i张相片对于空间第j点的成像坐标;Xij 0是所有与成像坐标(xij,yij)相关的参数,(vxij,vyij)是残余误差,Aij,Bij和Cij分别是步骤a中建立非线性畸变模型时获取的空间点在站位相机像面上最终的像点坐标对于第i图片的外方位参数、第j空间点坐标以及相机成像参数的雅可比矩阵;
其中,所述步骤g)中对误差方程雅可比矩阵中的各项进行自适应的比例调整的步骤为:
g1)在每次平差迭代过程开始,首先统计误差方程中对应单边内方位参数雅可比矩阵各列的数量级Mj及调整系数kj,计算方法如下:
g2)将对应单边内方位参数的雅可比矩阵的每一列计算结果乘以相应列的比例系数,得到调整后的雅可比矩阵为:
2.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤a)中建立线性成像模型的步骤为:
a1)确定物空间坐标系至相机空间坐标系之间的刚体旋转矩阵:
其中,Az,El,Ro为空间坐标系的三个旋转角;
a2)通过旋转矩阵和平移向量描述物空间至相机空间之间的刚体变换:
其中,R为步骤a1)确定的刚体旋转矩阵;
a3)在相机坐标系下,经小孔成像,将相机空间中的点投影至像面,得到对应的像点坐标:
3.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤a)中建立非线性畸变模型的步骤为:
a4)计算位于某一坐标(x,y)处的像点畸变量:
其中,相机的主点偏移为xp和yp,对应此空间点的畸变参数为K1ss’,K2ss’,K3ss’,P1,P2,B1,B2;
a5)获取空间点在站位相机像面上最终的像点坐标:
4.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤h)通过分块方式快速计算法方程中的各项的步骤为:
h1)计算Aij,Bij和Cij,分别是步骤c中获取的成像模型对于第i图片的外方位参数、第j空间点坐标以及相机成像单边内参数的雅可比矩阵;
h2)将求解得到的Aij,Bij和Cij,按照如下公式对法方程中的各项进行有规则的分块描述:从而在对图像、目标点的索引过程中逐步得到法方程中的矩阵:
h3)获取最终的法方程的分块模型
5.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤i)通过最小二乘技术迭代求解单边自标定中的各项参数的步骤为:
i1)利用N22项的分块对角性,改写法方程;
表示为:
i2)求解外方位参数和单边内方位参数的增量;
i3)将i2)的计算结果带入i1)中求解空间坐标增量其中利用N22的分块对角性进行简化;
6.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤j)参数比例逆调整,消除步骤g)中的人为定义的比例变化的方法为:确定未知量δ′与需要求解的δ之间的关系为:
7.根据权利要求1所述的场相关单边自标定光束平差方法,其特征是:所述步骤k)计算中误差,并进行单边自标定各项参数的不确定度估计的步骤为:
k1)计算单位权中误差:
其中,V是所有观测坐标的残余误差;n是观测坐标数目;u是未知变量数目;
k2)通过协方差矩阵的误差传播定律,计算外方位参数、内参数以及空间坐标的协方差矩阵:
其中,N为系数矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610298071.3A CN106023146B (zh) | 2016-05-06 | 2016-05-06 | 用于摄影测量中的场相关单边自标定光束平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610298071.3A CN106023146B (zh) | 2016-05-06 | 2016-05-06 | 用于摄影测量中的场相关单边自标定光束平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106023146A CN106023146A (zh) | 2016-10-12 |
CN106023146B true CN106023146B (zh) | 2018-10-30 |
Family
ID=57081943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610298071.3A Active CN106023146B (zh) | 2016-05-06 | 2016-05-06 | 用于摄影测量中的场相关单边自标定光束平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106023146B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106846467B (zh) * | 2017-01-23 | 2020-05-05 | 阿依瓦(北京)技术有限公司 | 基于每个相机位置优化的实体场景建模方法和*** |
CN107192375B (zh) * | 2017-04-28 | 2019-05-24 | 北京航空航天大学 | 一种基于航拍姿态的无人机多帧图像自适应定位校正方法 |
CN107707821B (zh) * | 2017-09-30 | 2020-11-06 | 努比亚技术有限公司 | 畸变参数的建模方法及装置、校正方法、终端、存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102889882A (zh) * | 2012-09-03 | 2013-01-23 | 北京信息科技大学 | 一种基于光束平差的三维重建方法 |
CN103247053A (zh) * | 2013-05-16 | 2013-08-14 | 大连理工大学 | 基于双目显微立体视觉的零件精确定位方法 |
-
2016
- 2016-05-06 CN CN201610298071.3A patent/CN106023146B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102889882A (zh) * | 2012-09-03 | 2013-01-23 | 北京信息科技大学 | 一种基于光束平差的三维重建方法 |
CN103247053A (zh) * | 2013-05-16 | 2013-08-14 | 大连理工大学 | 基于双目显微立体视觉的零件精确定位方法 |
Non-Patent Citations (6)
Title |
---|
Accurate depth dependent lens distortion models: an application to;Luis Alvarez 等;《Journal of Mathematical Imaging & Vision》;20110131;第39卷(第1期);第75-85页 * |
AUTOMATIC CAMERA CALIBRATION IN CLOSE-RANGE PHOTOGRAMMETRY;Clive S. Fraser;《Photogrammetric Engineering & Remote Sensing》;20130430;第79卷(第4期);第381-388页 * |
利用二维DLT及光束法平差进行数字摄像机标定;张永军 等;《武汉大学学报 信息科学版》;20021230;第27卷(第6期);第566-569页 * |
基于分块雅可比矩阵的最小二乘状态估计算法;王珍意 等;《继电器》;20050801;第33卷(第15期);摘要 * |
基于广域信息的电力***谐波状态估计的研究;吕云清;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20120415;第2012年卷(第04期);摘要 * |
随对焦状态与物距变化的畸变模型及标定方法;董明利 等;《仪器仪表学报》;20131220;第34卷(第12期);第2653-2659页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106023146A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102663767B (zh) | 视觉测量***的相机参数标定优化方法 | |
US9857172B1 (en) | Method for implementing high-precision orientation and evaluating orientation precision of large-scale dynamic photogrammetry system | |
CN108871373B (zh) | 一种基于俯仰滚转台和非线性优化的星敏感器标定方法 | |
CN102800096B (zh) | 一种摄像机参数的鲁棒性估计算法 | |
CN110345921B (zh) | 立体视场视觉测量及垂轴像差和轴向像差校正方法及*** | |
CN108801218B (zh) | 大尺寸动态摄影测量***的高精度定向及定向精度评价方法 | |
CN106023146B (zh) | 用于摄影测量中的场相关单边自标定光束平差方法 | |
CN104697463B (zh) | 一种双目视觉传感器的消隐特征约束标定方法及装置 | |
CN110223355B (zh) | 一种基于双重极线约束的特征标志点匹配方法 | |
CN110009667A (zh) | 基于罗德里格斯变换的多视点云全局配准方法 | |
CN109325927B (zh) | 工业相机摄影测量图像亮度补偿方法 | |
CN101149836A (zh) | 一种三维重构的双摄像机标定方法 | |
CN111369607A (zh) | 一种基于图片解析的预制构件拼装匹配方法 | |
Caja et al. | Optical measuring equipments. Part I: Calibration model and uncertainty estimation | |
CN101907705A (zh) | 通用的多源遥感影像几何校正模型联合平差方法 | |
CN105678088B (zh) | 一种靶标测头的平差优化算法 | |
CN108871204B (zh) | 摄影测量中长度测量相对误差的绝对评价方法 | |
CN113393507B (zh) | 无人机点云与地面三维激光扫描仪点云配准方法 | |
CN109242911A (zh) | 一种基于分区域的双目相机基本矩阵计算方法 | |
CN107590517B (zh) | 一种基于形状信息的图像相似度度量方法及图像配准方法 | |
US7310440B1 (en) | Replacement sensor model for optimal image exploitation | |
Sun et al. | A new method of camera calibration based on the segmentation model | |
CN109508482A (zh) | 一种用于复杂曲面表面轮廓度误差不确定度的计算方法 | |
CN111899304B (zh) | 远心光路畸变中心定位方法 | |
Mustaffar et al. | Improving area-based matching by using surface gradients in the pixel co-ordinate transformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |