運動圖像解析標定框架改進及算法

霍洪峰， 趙煥彬

(河北師范大學(xué) 體育學(xué)院, 河北省人體運動生物信息測評重點實驗室, 河北 石家莊 50024)

0 引 言

體育運動中，人體的運動多是復(fù)雜的三維空間運動，若要準確定量描述人體的運動特征，須進行三維立體攝像與解析。立體攝像測量方法是采用兩臺或多臺攝像機，從不同角度對同一研究對象進行同步拍攝，進而把所拍攝的平面影片或錄像帶進行數(shù)字化，獲得所需的人體三維運動的空間坐標，計算有關(guān)的運動學(xué)參數(shù)[1-3]。三維圖像解析是體育科研工作者獲取人體運動學(xué)特征常采用的研究方法[4-8]，人體運動的三維圖像解析技術(shù)關(guān)鍵是確定三維物體的實際坐標系到攝像機上圖像坐標系的映射關(guān)系。針對不同的成像模型學(xué)術(shù)界提出了一系列經(jīng)典的定標算法，主要有DLT線性算法、非線性算法、兩步算法和主動定標法等[9]，定標算法標定精度和可靠程度將直接影響到測量結(jié)果的精度[10]。許多學(xué)者對攝像機標定技術(shù)及算法進行了研究[11-12]，其中DLT算法被廣泛應(yīng)用于圖像解析，實際使用中具有如下的優(yōu)點：各相機光軸不需相交；相機位置可任意放置而不需測量；只需兩部相機即可獲取三維坐標；可使用更多的相機，使獲取的三維坐標精度更高等[13-14]。但在實際應(yīng)用中，由于算法的限制，選擇特征點時需要充滿整個空間，因此，一般選擇的標定物都比較復(fù)雜，國內(nèi)大多數(shù)體育科研單位和高校所使用標定框架為放射狀框。放射狀框架因其安裝與攜帶不便，標定球容易混淆及相互之間容易遮擋，方向桿固定等缺點，影響了運動學(xué)特征獲取的順利進行。

本文通過設(shè)計標定框架改進了放射狀框架存在的可操作性差等問題。基于DLT線性算法和張正友定標算法的優(yōu)點，提出了一種適用于體育運動圖像解析的定標算法，為編制三維解析系統(tǒng)時提高精度，實現(xiàn)人體運動特征的快速反饋提供可靠的理論與技術(shù)依據(jù)。

1 三維重構(gòu)算法

1.1 定標算法基礎(chǔ)

DLT算法是通過對標定框架的拍攝，然后求出攝影系數(shù)來實現(xiàn)的，即從攝影測量學(xué)的角度研究攝像機圖像和環(huán)境物體之間的關(guān)系，建立攝像機成像幾何的線性模型。由線性方程求解線性模型參數(shù)。

式中：(R,t)稱為外部參數(shù)，分別代表3×3的旋轉(zhuǎn)矩陣和3×1的平移向量;K稱為內(nèi)部參數(shù);(u0,v0)是主點(即光軸與圖像平面的交點)坐標;fu為圖像u軸尺度因子;fv為圖像v軸尺度因子。令P=K(R,t)，則P為一個3×4的矩陣，

這時有：

即

p11X+p12Y+p13Z+p14-p31Xu-p32Yu-

p33Zu-p34u=0

p21X+p22Y+p23Z+p24-p31Xv-p32Yv-

p33Zv-p34v=0

限定條件p34=1。為了提高精度，總是要使觀察值的個數(shù)多于求解聯(lián)立方程的個數(shù)，因此，選定空間的6個以上坐標點，將圖像坐標(ui,vi)和對應(yīng)的空間點坐標(X,Y,Z)代入上面的方程就能線性求解出攝像機參數(shù)(p11,p12,p13,p14,p21,p22,p23,p24,p31,p32,p33)[15]。

1.2 定標算法改進

從DLT標定方法出發(fā)，其中選擇的標定點限制在坐標平面上，利用位于X=0和Y=0兩平面上的標定球，將相應(yīng)的實際三維空間坐標和圖像二維坐標代入DLT線性算法中的線性方程求解出攝像機內(nèi)外參數(shù)，求解出攝像機拍攝的運動圖像上二維點的三維空間坐標。具體過程及數(shù)學(xué)推導(dǎo)如下：

(1) 建立空間坐標系XYOZ。

(2) 設(shè)置好攝像機各項參數(shù)，固定攝像機的位置，將標定板分別平行于XOZ平面和YOZ平面拍攝。

(3) 利用圖像分析程序得到圖像中特征點的坐標(xi,yi)，測量出此特征點的實際空間坐標(Xi,Yi,Zi)。

(4) 將(xi,yi)和(Xi,Zi)代入方程組：

解得p11,p13,p14,p21,p23,p24,p31,p33,p34。將(xi,yi)和(Yi,Zi)代入，得：

解得p12,p13,p14,p22,p23,p24,p32,p33,p34。

(5) 保持攝像機參數(shù)不變，拍攝運動圖像，將圖像上關(guān)節(jié)點坐標(u，v)代入下式

其中，P=pij(i=1,2,3，j=1,2,3,4)。用兩個攝像機拍攝的圖像解方程解出關(guān)節(jié)點實際空間坐標(X,Y,Z)。

由于該算法采用平面點重建圖像的三維空間坐標，要求解的方程的未知數(shù)減少(由12個減少到8個)，計算較DLT算法簡單，且易編程實現(xiàn)。

2 雙平面標定框架設(shè)計

雙平面標定框架采取折疊和收縮設(shè)計，標定框架主體結(jié)構(gòu)采取輕質(zhì)鎂合金，框架兩垂直面的穿球桿采取上下交錯設(shè)計，兩垂直面上的標定球涂成兩種不同色彩且排列規(guī)則。

2.1 標定框架主要套件及功能

(1) 90°定位桿轉(zhuǎn)套(見圖1中1)。定位桿可來回旋轉(zhuǎn)，框架折疊時可將定位桿旋至與穿球桿平行。

(2) 90°定位桿接頭套(見圖1中2)?？蓙砘鼗瑒樱惯B桿分開和連接。定位桿控制一定長度，保證兩平面相互垂直。

(3) 凸輪夾與內(nèi)置定位銷(見圖1中3)。凸輪夾可快速鎖緊和放松，便于快速拆裝框架。主干拆裝處用定位銷定位，保證精度。

(4) 方向標定桿(見圖1中4)。方向標定桿可在水平方向360°旋轉(zhuǎn)，方便坐標變換。

(5) 合葉(見圖1中5)。便于框架兩平面折疊。

圖1 雙平面標定框架工程圖

2.2 雙平面標定框架優(yōu)點

目前，在國內(nèi)廣泛使用的是放射狀框架，與放射狀框架相比，雙平面框架具有攜帶方便，操作簡便，不遮擋，易識別等優(yōu)點。

(1) 攜帶方便。放射狀框架收好后需要兩個比較大的長方形箱子和一個用來裝基架的袋子，攜帶困難。雙平面框架兩個相互垂直的標定面采取可折疊的設(shè)計，收好后大小和一個普通的手提箱相差無幾，無論是平常攜帶還是進入比賽現(xiàn)場都非常方便。

(2) 操作簡便。放射狀框架需要將9根標桿按照對應(yīng)順序一一固定在標定基架上。雙平面框架由于安裝了90°定位桿轉(zhuǎn)套、90°定位桿接頭套、凸輪夾與內(nèi)置定位銷等套件，便于框架的折疊和快速拆裝，并能保

證精度。

(3) 不遮擋，易識別。當攝像機角度不合適時，放射性框架就會出現(xiàn)標定桿相互遮擋的問題，在固定好標定框架后還需要在紙上畫出標定桿的順序草圖用于標定球的識別。雙平面框架兩垂直面的穿球桿采取上下交錯設(shè)計，因此不存在標定球間的遮擋問題。兩垂直面上標定球采用不同顏色并且排列規(guī)則，便于標定球的識別。

(4) 便于三維空間坐標變換。方向標定桿可在水平方向上360°轉(zhuǎn)動，通過矢量夾角運算，可以方便地進行三維空間坐標變換。

3 雙平面框架精度分析

3.1 框架精度分析實驗設(shè)計

以雙平面框架為研究對象，應(yīng)用三維攝像測量原理，兩臺JVC9800攝像機(拍攝頻率50 Hz)A、B分別放置于如圖2所示位置，使A、B兩攝像機主光軸的連線到拍攝區(qū)域分別為15.6 m和16.4 m，機高1.15 m。分別在拍攝區(qū)域放置雙平面框架、標桿長740.0 mm，進行拍攝。實驗測試圖如圖2所示。

圖2 框架精度驗證測試圖

3.2 框架精度分析實驗結(jié)果

線性方程解出的攝像機參數(shù)如表1所示。

表1 線性方程解出的攝像機參數(shù)

經(jīng)雙平面框架標定及三維重構(gòu)后得到的標桿長度實際值、標桿長度計算值和相對誤差數(shù)據(jù)見表2。從實驗數(shù)據(jù)可看出,研制的雙平面標定框架在三維測量中相對誤差小于5%，能夠滿足人體運動參數(shù)測量的需要。

表2 框架精度分析試驗數(shù)據(jù)

4 結(jié) 語

研制了在實際測量中簡單易行的雙平面標定框架?；贒LT算法提出了一種適合雙平面標定框架的標定算法，計算較DLT算法簡單，且易編程實現(xiàn)。研制的雙平面標定框架在三維測量中相對誤差小于5%，能夠滿足人體運動學(xué)參數(shù)測量的需要。

[1] 蔡云飛 .人體運動三維仿真與分析系統(tǒng)[D].天津：河北工業(yè)大學(xué)，2004.

[2] Song L M，Wang M P，Lu L，etal. High precision camera calibration in vision measurement [J]. Optics & Laser Technology，2007，39(7)：1413-1420.

[3] 徐 艷，袁同春，柴業(yè)宏. 基于三維圖像解析的籃球投籃技術(shù)生物力學(xué)分析[J]. 合肥工業(yè)大學(xué)學(xué)報(自然科學(xué)版 ),2012,35(10):1428-1431.

XU Yan, YUAN Tong-chun, CHAI Ye-hong. Biodynamic analysis of basketball shooting skills based on 3-D view analysis[J]. Journal of Hefei University of Technology, 2012, 35(10): 1428-1431.

[4] 敬龍軍，王 鵬. 中外優(yōu)秀 20 km 競走運動員髖關(guān)節(jié)運動特征的三維運動學(xué)分析[J]. 北京體育大學(xué)學(xué)報, 2012,35(10): 122-126.

JING Long-jun，WANG Peng. 3-Dimensional kinematics analysis of hip joint motion characteristics of 20 km walking race athletes at home and abroad[J]. Journal of Beijing Sport University, 2012,35(10): 122-126.

[5] 霍洪峰，趙煥彬，葉 靜，等.中國運動生物力學(xué)研究發(fā)展現(xiàn)狀及發(fā)展趨勢[J].中國體育科技, 2009,45(3): 84-86.

HUO Hong-feng,ZHAO Huan-bin,YE Jing,etal. The current situation and developing trend of sports biomechanics in China[J]. China Sport Science and Technology, 2009,45(3): 84-86.

[6] 隋新梅.張峻旋轉(zhuǎn)推鉛球過渡階段主要關(guān)節(jié)角度變化特征分析[J].體育科研,2012,33(3):81-84.

SUI Xin-mei. Analysis of the variation characteristics of the main joint angles in the transition phase of Zhang Jun’s spinning shot put[J]. Sport Science Research, 2012,33(3):81-84.

[7] 文世林.我國兩名優(yōu)秀男子鉛球運動員旋轉(zhuǎn)推鉛球技術(shù)的三維運動學(xué)分析[D].北京：首都體育學(xué)院，2009.

[8] 周長濤. 男子跆拳道高水平運動員下劈技術(shù)動作三維運動學(xué)特征研究[J].山東體育學(xué)院學(xué)報，2012,28(2):79-84.

ZHOU Chang-tao. Three dimensional kinematic features of men’s taekwondo athletes’ axe kick[J]. Journal of Shandong Institute of Physical Education and Sports，2012, 28(2): 79-84.

[9] 毛劍飛，鄒細勇.改進的平面模板兩步法標定攝像機[J].中國圖象圖形學(xué)報，2004，9(7): 846-852.

MAO Jian-fei, ZOU Xi-yong. Improved two-stages camera calibration from a plane[J]. Journal of Image and Graphics，2004，9(7): 846-852.

[10] 李中偉，王從軍，史玉升. 3D 測量系統(tǒng)中的高精度攝像機標定算法[J]. 光電工程, 2008, 35(4): 58-63.

LI Zhong-wei，WANG Cong-jun，SHI Yu-sheng. High precision camera calibration algorithm for 3D measurement system[J]. Opto-Electronic Engineering, 2008, 35(4): 58-63.

[11] 陳 召，柏開祥，康 帆. 三維攝影框架的設(shè)計、 標定、 檢驗與應(yīng)用研究[J]. 武漢體育學(xué)院學(xué)報，2013，47(1)：65-69.

CHEN Zhao, BAI Kai-xiang, KANG Fan. Design, marking, checking and utilizing of 3D photography frame[J]. Journal of Wuhan Institute of Physical Education，2013，47(1)：65-69.

[12] 余 峰．基于圖像的三維造型算法研究[D].杭州：浙江大學(xué)，2004.

[13] John H, Challis, David G, Kerwin. Accuracy assessment and control point configuration when using the DLT for photogrammetry[J]. J Biomechanics, 1992, 35(9): 1053-1058.

[14] 楊年峰, 王季軍, 黃昌華，等. DLT算法中象平面誤差對三維重構(gòu)的影響[J]. 中國生物醫(yī)學(xué)工程學(xué)報，2000,19(2):121-129.

Yang Nian-feng, Wang Ji-jun, Huang Chang-hua,etal. The influences of image coordinate errors on the accuracy of 3-D coordinate reconstruction in DLT algorithm[J]. Chinese Journal of Biomedical Engineering，2000,19(2):121-129.

[15] 鄭秀瑗. 現(xiàn)代運動生物力學(xué)[M]. 北京：國防工業(yè)出版社, 2007:13-24.