[发明专利]一种基于定焦数码相机的旋转全景摄影测量方法和装置

摘要

本发明提出了一种基于定焦数码相机的旋转全景摄影测量方法和装置，当难以在被测物体上或周围布设控制点，但可以在其它位置或被测物对面布设控制点的情况下，利用少量控制点，可解算被测点的物方坐标，是解决现实场景中相机视场角小、控制点稀少难题的技术方案，是在困难场景进行近景摄影测量的一种有效的手段；本发明设备简单、成本较低，能够解决实际工作中控制点稀少时的摄影测量问题。

主权项

一种基于定焦数码相机的旋转全景相机的摄影测量方法，其特征在于，包括以下步骤：步骤1：在全景三维控制场中选择M个摄站，在每一摄站用定焦数码相机拍摄含控制点的影像，共M张影像；步骤2：计算定焦数码相机的内方位元素(f,x0,y0)、畸变参数(K1,K2,P1,P2)；提取M张影像上控制点的像点坐标(x,y)，利用已知的控制点像点坐标(x,y)、控制点在全景三维控制场的物方空间坐标(X,Y,Z)和共线方程式，见公式(1)，采用最小二乘的平差方法，解求定焦数码相机的内方位元素(f,x0,y0)、外方位元素(Xs,Ys,Zs,a1,a2,a3,b1,b2,b3,c1,c2,c3)、畸变参数(K1,K2,P1,P2)，f,x0,y0和K1,K2,P1,P2的结果作为真值参与后续计算；公式(1)如下所示：x-x0-Δx=-fa1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)y-y0-Δy=-fa2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)---(1)]]>式中，(x0,y0)为像主点坐标，Δx,Δy为像点畸变差，Δx=(x-x0)(K1r2+K2r4)+P1(r2+2(x-x0)2)+2P2(x-x0)·(y-y0)Δy=(y-y0)(K1r2+K2r4)+P2(r2+2(y-y0)2)+2P1(x-x0)·(y-y0),]]>式中，K1,K2为径向畸变参数，P1,P2为偏心畸变参数，r为M张影像上控制点的像点到像主点的距离，r＝(x‑x0)2+(y‑y0)2，f为定焦数码相机的焦距，Xs,Ys,Zs为投影中心S在全景三维控制场物方空间坐标系(G‑X,Y,Z)中的坐标，a1,a2,a3,b1,b2,b3,c1,c2,c3为从像空间坐标系到全景三维控制场物方空间坐标系的旋转矩阵中的元素；步骤3：在全景三维控制场中选择N个摄站，在每一摄站上用旋转全景相机拍摄360°水平角的全景影像，共N张全景影像；步骤4：计算旋转全景相机的旋转参数(US,VS,WS)和RS；提取N张全景影像上控制点的像点坐标(x1,y1)；根据步骤2中的(f,x0,y0)、(K1,K2,P1,P2)、全景三维控制场中控制点的像点坐标(x1,y1)、全景三维控制场中控制点的物方空间坐标(X,Y,Z)和公式(2)，采用最小二乘的平差方法解求旋转中心在全景三维控制场物方空间坐标系中的坐标(XO,YO,ZO)、旋转坐标系在全景三维控制场物方空间坐标系中的初始角度α0和旋转全景相机的旋转参数；旋转全景相机的旋转参数包括定焦数码相机投影中心与旋转中心的偏移量(US,VS,WS)和相对姿态RS；公式(2)如下所示：x1-x0-Δx=-f(a1′cosα0-c1′sinα0)(X-XO)+b1′(Y-YO)+(a1′sinα0+c1′cosα0)(Z-ZO)-a1′US-b1′VS-c1′WS(a3′cosα0-c3′sinα0)(X-XO)+b3′(Y-YO)+(a3′sinα0+c3′cosα0)(Z-ZO)-a3′US-b3′VS-c3′WSy1-y0-Δy=-f(a2′cosα0-c2′sinα0)(X-XO)+b2′(Y-YO)+(a2′sinα0+c2′cosα0)(Z-ZO)-a2′US-b2′V-c2′WS(a3′cosα0-c3′sinα0)(X-XO)+b3′(Y-YO)+(a3′sinα0+c3′cosα0)(Z-ZO)-a3′US-b3′VS-c3′WS---(2)]]>式中，(x0,y0)为像主点坐标，Δx,Δy为像点畸变差，Δx=(x1-x0)(K1r12+K2r14)+P1(r12+2(x1-x0)2)+2P2(x1-x0)·(y1-y0)Δy=(y1-y0)(K1r12+K2r14)+P2(r12+2(y1-y0)2)+2P1(x1-x0)·(y1-y0),]]>式中，K1,K2为径向畸变参数，P1,P2为偏心畸变参数，r1为N张影像上控制点的像点到像主点的距离，r1＝(x1‑x0)2+(y1‑y0)2，f为定焦数码相机的焦距，X′,Y′,Z′为全景三维控制场内控制点的物方空间坐标，(US,VS,WS)为投影中心S在旋转坐标系(O‑U,V,W)中的坐标，a1′,a2′,a3′,b1′,b2′,b3′,c1′,c2′,c3′为从像空间坐标系到旋转坐标系的旋转矩阵中的元素，(XO,YO,ZO)为旋转坐标系与全景三维控制场物方空间坐标系的相对位置，α0为旋转坐标系在全景三维控制场物方空间坐标系中的初始角度；步骤5：在测量现场布设3个以上的控制点Q1,Q2,Q3……；步骤6：在测量现场选择T个摄站，在每一摄站上拍摄含Q1,Q2,Q3……和待测点的测量现场的全景影像，共T张全景影像；步骤7：计算测量现场中旋转全景相机的旋转中心在测量现场物方空间坐标系中的坐标(X′O,Y′O,Z′O)和旋转坐标系在测量现场物方空间坐标系中的初始角度α′0；提取测量现场的T张全景影像中控制点Q1,Q2,Q3……的像点坐标(x2,y2)；根据步骤2中的(f,x0,y0)、(K1,K2,P1,P2)、步骤4计算得到的全景相机的旋转参数US,VS,WS、RS、控制点Q1,Q2,Q3……的像点坐标(x2,y2)、控制点在测量现场物方空间坐标系中的坐标(XQ,YQ,ZQ)和公式(3)，采用最小二乘的平差方法解求旋转全景相机在测量现场的旋转中心坐标(X′O,Y′O,Z′O)、旋转坐标系在测量现场物方空间坐标系中的初始角度α′0；公式(3)如下所示：x2-x0-Δx=-f(a1′cosα0′-c1′sinα0′)(XP-XO′)+b1′(YP-YO′)+(a1′sinα0′+c1′cosα0′)(ZP-ZO′)-a1′US-b1′VS-c1′WS(a3′cosα0′-c3′sinα0′)(XP-XO′)+b3′(YP-YO′)+(a3′sinα0′+c3′cosα0′)(ZP-ZO′)-a3′US-b3′VS-c3′WSy2-y0-Δy=-f(a2′cosα0′-c2′sinα0′)(XP-XO′)+b2′(YP-YO′)+(a2′sinα0′+c2′cosα0′)(ZP-ZO′)-a2′US-b2′V-c2′WS(a3′cosα0′-c3′sinα0′)(XP-XO′)+b3′(YP-YO′)+(a3′sinα0′+c3′cosα0′)(ZP-ZO′)-a3′US-b3′VS-c3′WS---(3)]]>式中，(x2,y2)为像点坐标，(x0,y0)为像主点坐标，Δx,Δy为像点畸变差，K1,K2为径向畸变参数，P1,P2为偏心畸变参数，r2为T张全景影像中控制点Q1,Q2,Q3,……的像点到像主点的距离r2＝(x2‑x0)2+(y2‑y0)2，f为定焦数码相机的焦距，XP,YP,ZP为测量现场内控制点的物方空间坐标，(US,VS,WS)为投影中心S在旋转坐标系(O‑U,V,W)中的坐标，a′1,a′2,a′3,b′1,b′2,b′3,c′1,c′2,c′3为从像空间坐标系到旋转坐标系的旋转矩阵中的元素，(X′O,Y′O,Z′O)为旋转坐标系与测量现场物方空间坐标系的相对位置，α′0为旋转坐标系在测量现场物方空间坐标系的初始角度；步骤8：由α′0和旋转角度β＝(360/C)°，根据公式(4)计算每张像片的旋转角度α；公式(4)如下所示：α＝α′0+(i‑1)β (i＝1,2,…,C) (4)其中，C为转动定焦数码相机的次数；步骤9：计算待测点在测量现场物方空间坐标系中的坐标(XA,YA,ZA)；提取测量现场的T张全景影像上待测点的像点坐标(x3,y3)；利用已知的f,x0,y0、K1,K2,P1,P2、(US,VS,WS)、RS、(x3,y3)、(X′O,Y′O,Z′O)、α和公式(5)，采用最小二乘的平差方法解求待测点在测量现场物方空间坐标系中的坐标(XA,YA,ZA)；公式(5)如下所示：x3-x0-Δx=-f(a1′cosα0′-c1′sinα0′)(XA-XO′)+b1′(YA-YO′)+(a1′sinα0′+c1′cosα0′)(ZA-ZO′)-a1′US-b1′VS-c1′WS(a3′cosα0′-c3′sinα0′)(XA-XO′)+b3′(YP-YO′)+(a3′sinα0′+c3′cosα0′)(ZA-ZO′)-a3′US-b3′VS-c3′WSy3-y0-Δy=-f(a2′cosα0′-c2′sinα0′)(XA-XO′)+b2′(YA-YO′)+(a2′sinα0′+c2′cosα0′)(ZA-ZO′)-a2′US-b2′V-c2′WS(a3′cosα0′-c3′sinα0′)(XA-XO′)+b3′(YA-YO′)+(a3′sinα0′+c3′cosα0′)(ZA-ZO′)-a3′US-b3′VS-c3′WS---(5)]]>式中，(x3,y3)为像点坐标，(x0,y0)为像主点坐标，Δx,Δy为像点畸变差，K1,K2为径向畸变参数，P1,P2为偏心畸变参数，r3为T张全景影像上待测点的像点到像主点的距离，r3＝(x3‑x0)2+(y3‑y0)2，f为定焦数码相机的焦距，XA,YA,ZA为待测点在测量现场物方空间坐标系中的坐标，(US,VS,WS)为投影中心S在旋转坐标系(O‑U,V,W)中的坐标，a′1,a′2,a′3,b′1,b′2,b′3,c′1,c′2,c′3为从像空间坐标系到旋转坐标系的旋转矩阵中的元素，(X′O,Y′O,Z′O)为旋转坐标系与测量现场物方空间坐标系的相对位置，α′0为旋转坐标系在测量现场物方空间坐标系的初始角度；步骤10：结束。