季云峰, 陸愛發(fā), 任 杰, 施之皓

(上海體育學(xué)院 中國乒乓球?qū)W院,上海 200438)

?

基于乒乓球機(jī)器人視覺系統(tǒng)的單色乒乓球旋轉(zhuǎn)三維速度測定

季云峰, 陸愛發(fā), 任 杰, 施之皓

(上海體育學(xué)院 中國乒乓球?qū)W院,上海 200438)

旋轉(zhuǎn)是乒乓球最核心的制勝因素,故獲得乒乓球的旋轉(zhuǎn)信息是乒乓球機(jī)器人視覺中必不可少的任務(wù)。提出基于單目視覺的測算乒乓球旋轉(zhuǎn)三維速度的方向及大小的算法,用于測量無標(biāo)記的乒乓球旋轉(zhuǎn),將其應(yīng)用于乒乓球機(jī)器人的視覺識(shí)別,獲得乒乓球的旋轉(zhuǎn)信息,對(duì)機(jī)器人后續(xù)運(yùn)動(dòng)規(guī)劃和控制決策起重要的作用。對(duì)所提算法進(jìn)行驗(yàn)證,得到與理論結(jié)果誤差較小的實(shí)驗(yàn)數(shù)據(jù),驗(yàn)證了方法的可行性。

乒乓球; 視覺系統(tǒng); 三維速度; 旋轉(zhuǎn); 機(jī)器人

Author’s address China Table Tennis College, Shanghai University of Sport, Shanghai 200438, China

乒乓球機(jī)器人是研究快速視覺感知和運(yùn)動(dòng)控制的良好平臺(tái)。自20世紀(jì)80年代Billingsley提出舉辦乒乓球機(jī)器人比賽的倡議以來,各國先后展開了乒乓球機(jī)器人方面的研究[1]。Andersson[2-3]、Acosta等[4],Matsushima等[5-7]先后設(shè)計(jì)了不同的機(jī)器人系統(tǒng),成功實(shí)現(xiàn)了人機(jī)對(duì)打。

實(shí)現(xiàn)乒乓球機(jī)器人的人機(jī)對(duì)打,主要需要機(jī)器視覺和運(yùn)動(dòng)控制的共同配合。先通過機(jī)器視覺找到乒乓球的空間位置,再通過對(duì)乒乓球進(jìn)行動(dòng)力學(xué)分析或經(jīng)驗(yàn)判斷預(yù)測乒乓球的運(yùn)動(dòng)軌跡,傳達(dá)給控制端,選擇一個(gè)合適的回球落點(diǎn),將乒乓球拍移動(dòng)到相應(yīng)位置,完成擊球。在乒乓球運(yùn)動(dòng)中,旋轉(zhuǎn)球是運(yùn)動(dòng)員使用最多也是最具特色的一項(xiàng)技術(shù),是運(yùn)動(dòng)員比賽制勝的重要法寶。一直以來很少有人根據(jù)旋轉(zhuǎn)運(yùn)動(dòng)過程對(duì)運(yùn)動(dòng)軌跡進(jìn)行分析,旋轉(zhuǎn)球基本未應(yīng)用到乒乓球機(jī)器人之中。目前無論是國外還是國內(nèi)乒乓球機(jī)器人,大多數(shù)的預(yù)測軌跡工作都是圍繞無旋球展開的,針對(duì)旋轉(zhuǎn)球的預(yù)測目前只集中在經(jīng)驗(yàn)學(xué)習(xí)方法中的模糊預(yù)測。Zhang等[8]根據(jù)受力情況采用近似的物理模型對(duì)乒乓球運(yùn)動(dòng)軌跡進(jìn)行預(yù)測,優(yōu)點(diǎn)是預(yù)測時(shí)間短,可以滿足實(shí)時(shí)性,但其中參數(shù)的不確定性導(dǎo)致不能精確建模,難以推廣到旋轉(zhuǎn)球的應(yīng)用中。

實(shí)際打球過程中,球拍與乒乓球、球桌與乒乓球之間都有摩擦力的作用,使得乒乓球所受合力偏離了球的中心而產(chǎn)生旋轉(zhuǎn)?，F(xiàn)階段也有文獻(xiàn)對(duì)乒乓球的旋轉(zhuǎn)進(jìn)行探討,Leung等[9]詳細(xì)分析了旋轉(zhuǎn)球與球拍和球臺(tái)碰撞過程,并創(chuàng)建了一種旋轉(zhuǎn)球與球拍和球臺(tái)的反彈模型。Akira等[10-12]在分析了旋轉(zhuǎn)球與球臺(tái)反彈模型的基礎(chǔ)上,提出相應(yīng)的旋轉(zhuǎn)球的飛行軌跡預(yù)測算法,但這些算法都未針對(duì)特定的旋轉(zhuǎn)類型和大小設(shè)計(jì)。

以上這些算法都是通過對(duì)乒乓球進(jìn)行受力分析預(yù)測乒乓旋轉(zhuǎn)球的軌跡,但該預(yù)測只能獲得旋轉(zhuǎn)的大概信息,到目前為止缺乏準(zhǔn)確計(jì)算無標(biāo)記乒乓球旋轉(zhuǎn)的算法。有學(xué)者在乒乓球上做標(biāo)記測算旋轉(zhuǎn)速度,Nakashima等[13-15]在球上做標(biāo)記通過高速視覺對(duì)球的旋轉(zhuǎn)進(jìn)行測算,但在乒乓球上做標(biāo)記會(huì)增加乒乓球的質(zhì)量,對(duì)乒乓球的圓度也會(huì)產(chǎn)生一定影響,在比賽中不可能使用添加標(biāo)記的乒乓球,故該方法并不能被廣泛應(yīng)用于測算乒乓球旋轉(zhuǎn)中。

到目前為止,并沒有學(xué)者提出一個(gè)精確的不做標(biāo)記的方法對(duì)球體旋轉(zhuǎn)的三維速度進(jìn)行測量。本文主要針對(duì)球體旋轉(zhuǎn)速度大小及方向問題,利用球體旋轉(zhuǎn)中的一些特征,提出一種基于單目攝像機(jī)的測量球體三維速度的方法,并用實(shí)驗(yàn)驗(yàn)證算法的有效性。

1 實(shí)驗(yàn)準(zhǔn)備

目前乒乓球機(jī)器人的視覺系統(tǒng)一般由多臺(tái)攝像機(jī)共同定位乒乓球的空間位置,事先對(duì)攝像機(jī)進(jìn)行標(biāo)定,根據(jù)乒乓球的特征對(duì)其進(jìn)行識(shí)別跟蹤。由于大部分乒乓球機(jī)器人的視覺系統(tǒng)都是針對(duì)無旋轉(zhuǎn)乒乓球進(jìn)行識(shí)別,故將乒乓球當(dāng)成一個(gè)白色球體進(jìn)行識(shí)別,獲得其中心的空間坐標(biāo)。

本文提出一種全新的測算乒乓球旋轉(zhuǎn)速度和方向的方法,主要利用乒乓球面的商標(biāo)相對(duì)于球的運(yùn)動(dòng)計(jì)算旋轉(zhuǎn)速度。同樣利用高速攝像機(jī)作為視覺系統(tǒng),高速攝像機(jī)的幀率保證能看到連續(xù)幾幀乒乓球商標(biāo)即可。在機(jī)器人視覺系統(tǒng)中,如果使用多目攝像機(jī)進(jìn)行乒乓球定位,會(huì)由于定標(biāo)、環(huán)境的影響造成誤差多倍放大,故本文采用的是基于單目攝像機(jī)測算旋轉(zhuǎn)速度的方法。由于乒乓球在空中是有旋轉(zhuǎn)的,故攝像機(jī)架在球桌四周皆可,只要幀率跟上球的旋轉(zhuǎn)速度即可,一般將攝像機(jī)放在裁判席附近。因?yàn)樵趯?shí)驗(yàn)中需要識(shí)別乒乓球的商標(biāo),所以要確保攝像機(jī)能夠拍到商標(biāo),故在實(shí)驗(yàn)中采用2臺(tái)攝像機(jī)放在乒乓球桌的兩邊,并且高度不同。在實(shí)驗(yàn)過程中只對(duì)可以識(shí)別商標(biāo)的那臺(tái)攝像機(jī)進(jìn)行算法處理,故覆蓋范圍較大,在大部分情況下可以識(shí)別商標(biāo),所以對(duì)于前后旋轉(zhuǎn)的乒乓球也可以進(jìn)行識(shí)別,本文的算法基本適用于任何旋轉(zhuǎn)速度的識(shí)別計(jì)算。

由于本文測算旋轉(zhuǎn)主要依靠乒乓球面的商標(biāo),故選取所拍視頻中的一張圖像進(jìn)行處理實(shí)驗(yàn),經(jīng)過二值化,腐蝕和膨脹,查找輪廓,以乒乓球中心點(diǎn)作為原點(diǎn)建立坐標(biāo)系,可以獲得商標(biāo)中心的坐標(biāo)(圖1)。在乒乓機(jī)器人視覺系統(tǒng)中,獲得幾幀圖片中乒乓球的商標(biāo)中心坐標(biāo),再利用這幾個(gè)中心坐標(biāo)間的關(guān)系測算乒乓球的旋轉(zhuǎn)速度和方向。

2 乒乓球旋轉(zhuǎn)速度測試方法

2.1 乒乓球旋轉(zhuǎn)模型 乒乓球在空中的旋轉(zhuǎn)必將繞某一個(gè)旋轉(zhuǎn)軸進(jìn)行旋轉(zhuǎn),故乒乓球上任意一點(diǎn)旋轉(zhuǎn)一周必將形成一個(gè)圓形。用單目攝像機(jī)從某一角度拍,這個(gè)角度不一定是正對(duì)旋轉(zhuǎn)軸,若將商標(biāo)中心作為標(biāo)記點(diǎn),形成的軌跡在單目視覺中將形成一個(gè)橢圓。本文需要基于單目攝像機(jī)對(duì)旋轉(zhuǎn)乒乓球進(jìn)行拍攝,利用旋轉(zhuǎn)乒乓球的旋轉(zhuǎn)特征計(jì)算乒乓球的旋轉(zhuǎn)速度的大小和方向;先利用乒乓球的旋轉(zhuǎn)建立一個(gè)旋轉(zhuǎn)球體模型,對(duì)旋轉(zhuǎn)乒乓球的特征進(jìn)行分析(圖2)。

圖1 獲得乒乓球商標(biāo)中心坐標(biāo)

Figure 1. The acquisition of the coordinate of table tennis brand’s center

圖2 乒乓球旋轉(zhuǎn)模型

假設(shè)乒乓球在空間坐標(biāo)系Oxyz中進(jìn)行旋轉(zhuǎn),圖2中的2個(gè)圓代表旋轉(zhuǎn)的乒乓球,其中實(shí)線圓代表乒乓球投影在xOy平面形成的圓,虛線圓代表球體上任意一個(gè)截面所形成的圓。2個(gè)方框代表2個(gè)圓形所在的截面,在空間中2個(gè)平面相交會(huì)產(chǎn)生一條交線,在圖中2個(gè)平面相交的交線即為中間那條黑色直線;而這2個(gè)圓形截面的交線,即為虛線圓的直徑,記為直徑m,屬于2個(gè)平面交線的一部分。

2.2 乒乓球三維速度分解 高速旋轉(zhuǎn)的乒乓球必將沿著某條通過球心的軸轉(zhuǎn)動(dòng),故可將求解高速旋轉(zhuǎn)的乒乓球的三維速度方向問題轉(zhuǎn)化為求解轉(zhuǎn)動(dòng)軸空間位置問題?，F(xiàn)假設(shè)此球體旋轉(zhuǎn)方向?yàn)檠刂摼€圓旋轉(zhuǎn),則此球體的轉(zhuǎn)動(dòng)軸為垂直于虛線圓所在截面的直徑m的中垂線所在直線,即為直線L。

如圖3所示,將高速旋轉(zhuǎn)的乒乓球投影到xOy平面上,可以得到球體在新的坐標(biāo)系x′O′y′的平面圖形。此時(shí),空間坐標(biāo)中的實(shí)線圓保持不變,而虛線圓將在平面中被投影成一個(gè)橢圓。此外,空間中的轉(zhuǎn)動(dòng)軸L將被投影到平面x′O′y′中,成為橢圓長軸m′的中垂線L′。故可將求解旋轉(zhuǎn)乒乓球轉(zhuǎn)動(dòng)軸L空間位置問題轉(zhuǎn)化為2個(gè)問題:①直線L′在平面x′O′y′的位置;②直線L′與平面x′O′y′的夾角。

圖3 旋轉(zhuǎn)乒乓球投影

2.3 求解旋轉(zhuǎn)軸方向

2.3.1 直線L′在平面x′O′y′的位置 先解決直線L′在平面x′O′y′位置問題,可轉(zhuǎn)化為求解橢圓長軸m′與x軸傾角及橢圓中心點(diǎn)坐標(biāo)問題。對(duì)于任意一個(gè)橢圓,都可以用其一般方程

Ax2+Bxy+Cy2+Dx+Ey+1=0

表示,總共有A、B、C、D、E5個(gè)系數(shù),即只要找到橢圓上面任何5個(gè)點(diǎn)的坐標(biāo)就可以解出橢圓的一般方程。在乒乓球旋轉(zhuǎn)測定中,可以利用乒乓球上面的商標(biāo)中心作為識(shí)別對(duì)象,記下5幀圖像中商標(biāo)中心位置即可。假設(shè)找到橢圓上任意5個(gè)點(diǎn)的坐標(biāo)(xn,yn),n=1,2,3,4,5,則可以得到以下五元一次方程組:

(1)

在方程組(1)中,x、y都是已知數(shù),未知數(shù)是系數(shù)A、B、C、D、E,它們都是一次的,故解出這個(gè)方程組的解是唯一的,利用5張圖像即可求出橢圓軌跡方程。

利用Matlab軟件編程求解方程(1)的解,獲得系數(shù)值,從而得到橢圓一般方程之后,開始計(jì)算橢圓長軸m′與x軸傾角及橢圓中心點(diǎn)坐標(biāo)。令

F(x,y)= Ax2+Bxy+Cy2+Dx+Ey+1

(2)

由于橢圓屬于中心型曲線,故中心點(diǎn)坐標(biāo)(x0,y0)滿足

(3)

代入F(x,y),有

(4)

即得橢圓中心坐標(biāo)為(x0,y0)。

(5)

令F(x,y)=0,將式(5)代入,可得

(6)

其中,

(7)

2(C-A)sinβcosβ+B(cos2β-sin2β)=0

(8)

利用二倍角公式,有

(9)

化簡得

(10)

即可得到橢圓長軸m′與x軸傾角β的值:

(11)

至此,計(jì)算出乒乓球的旋轉(zhuǎn)軸L在平面x′O′y′內(nèi)的投影L′與x、y軸夾角分別為:

2.3.2 直線L′與平面x′O′y′的夾角 接下來解決直線L′與平面x′O′y′夾角,在單目攝像機(jī)視角下,一個(gè)與攝像機(jī)視角平面傾角為α的平面中的一個(gè)圓會(huì)被識(shí)別成一個(gè)橢圓。其主要原因是將原本平面中圓投射到攝像機(jī)視角平面中,從而導(dǎo)致圖像變形,橢圓的短軸就是因?yàn)樵局睆酵渡涞叫缕矫嫦驴s短所形成的。

如圖4所示,假設(shè)攝像機(jī)視角平面為x1Oy1平面,則計(jì)算機(jī)視角平面與圓所在平面的夾角為α,投影形成的三角形的三條邊長分別為a、b、c,原本圓的直徑為c,投影到x1Oy1平面上變成橢圓短軸b,則有

(12)

此處的夾角α即為直線L′與平面x′O′y′夾角,三角形的邊長b和c即為橢圓的長軸和短軸,記橢圓長軸為a′,短軸為b′,可得

(13)

圖4 空間夾角轉(zhuǎn)換關(guān)系

故可將問題轉(zhuǎn)化為求解橢圓長軸a′和短軸b′問題。根據(jù)前文計(jì)算結(jié)果,可知橢圓中心坐標(biāo)為(x0,y0),假設(shè)一條直線過橢圓中心點(diǎn)與x軸夾角為φ,則這條直線的參數(shù)方程為

(14)

式中:t為參數(shù);φ為直線與x軸夾角。

令F(x)=0,將式(14)代入,可得

(Acos2φ+Bsinφcosφ+Csin2φ)t2+

[(2Ax0+By0+D)cosφ+

(Bx0+2Cy0+E)sinφ]t+H=0

(15)

式中,

(16)

將前文算出的中心坐標(biāo)值代入式(16),可得

(17)

(18)

易知

t1+t2=0

(19)

故可得

(20)

由韋達(dá)定理可得

(21)

故可得

(22)

由式(22)可求得傾角為φ的以(x0,y0)為中點(diǎn)的弦的長度,因?yàn)闄E圓的長軸、短軸分別是其弦中最長和最短的兩根弦,當(dāng)Acos2φ+Bsinφcosφ+Csin2φ取得極值時(shí),對(duì)應(yīng)的弦即為橢圓的長軸和短軸。

令G=Acos2φ+Bsinφcosφ+Csin2φ,則

(23)

經(jīng)化簡,可得

G=(A+C)/2+

(24)

易知

(25)

故可得

(26)

由此可求得橢圓的長軸a′和短軸b′長度如下:

(27)

根據(jù)式(13)、(17)、(27)可得直線L′與平面x′O′y′夾角α的計(jì)算公式如下:

(28)

根據(jù)式(28),計(jì)算出乒乓球的旋轉(zhuǎn)軸L與空間坐標(biāo)系z(mì)軸夾角為:

(29)

至此,計(jì)算出乒乓球旋轉(zhuǎn)軸與空間3個(gè)軸的夾角,即算出乒乓球旋轉(zhuǎn)速度的三維方向。

2.4 求解旋轉(zhuǎn)速度大小 精確地計(jì)算速度大小的方法是先計(jì)算任意2幀圖像中同一點(diǎn)運(yùn)動(dòng)的弧線的距離,然后除以2幀圖像相隔時(shí)間,即得到球體轉(zhuǎn)速大小。其計(jì)算量會(huì)相對(duì)較大,在高速旋轉(zhuǎn)球體上,可以將兩點(diǎn)間的弧線距離近似為兩點(diǎn)間的直線距離,從而可以大幅度降低計(jì)算量,縮短計(jì)算時(shí)間。本文選取了橢圓上的5個(gè)點(diǎn),即同一個(gè)點(diǎn)連續(xù)5幀的運(yùn)動(dòng)軌跡。為了降低誤差,將選取這5個(gè)點(diǎn)中每相鄰2點(diǎn)之間距離的平均值。

根據(jù)前文可知,橢圓上任意5個(gè)點(diǎn)的坐標(biāo)為(xn,yn),n=1,2,3,4,5,記球體旋轉(zhuǎn)速度大小為v,每2幀之間間隔時(shí)間為t,物體實(shí)際大小與圖像中像素大小比值為k,則得出

(30)

式(30)可以計(jì)算出的v為選取坐標(biāo)點(diǎn)的線速度,可以根據(jù)線速度v計(jì)算球體旋轉(zhuǎn)的角速度。球體是繞旋轉(zhuǎn)軸進(jìn)行旋轉(zhuǎn)的,本文所選取的點(diǎn)的運(yùn)動(dòng)軌跡是垂直于旋轉(zhuǎn)軸的一個(gè)圓,但是這個(gè)圓不一定會(huì)經(jīng)過球心,所以不能簡單地利用球體的半徑計(jì)算角速度。

本文計(jì)算角速度利用所選取點(diǎn)運(yùn)動(dòng)軌跡圓的半徑,而在單目視覺中,該半徑即為橢圓軌跡的長軸a′。故利用式(28)和式(30)可算得球體旋轉(zhuǎn)的角速度計(jì)算公式為:

(31)

2.5 實(shí)驗(yàn)驗(yàn)證 為了驗(yàn)證設(shè)計(jì)算法的正確性,本文利用攝像機(jī)拍攝視頻進(jìn)行實(shí)驗(yàn)驗(yàn)證。將2臺(tái)攝像機(jī)擺放在乒乓球桌兩邊,用一個(gè)智能發(fā)球機(jī)進(jìn)行發(fā)球,發(fā)球的旋轉(zhuǎn)速度事先調(diào)試、測試好,控制在60 r/s,再利用2個(gè)攝像機(jī)同時(shí)進(jìn)行拍攝。攝像機(jī)的幀率采用100幀/s,最終用拍攝到乒乓球商標(biāo)的視頻進(jìn)行實(shí)驗(yàn)驗(yàn)證。

2臺(tái)攝像機(jī)放在球桌兩邊,且高度有一點(diǎn)差異,便于確保有一臺(tái)攝像機(jī)可以拍攝到乒乓球的商標(biāo)運(yùn)動(dòng)軌跡,再針對(duì)拍攝的視頻進(jìn)行算法代入,從而進(jìn)行實(shí)驗(yàn)驗(yàn)證。拍攝的視頻截圖如圖5所示,經(jīng)過二值化、去噪、腐蝕膨脹、查找輪廓燈一系列圖像處理操作,可找出乒乓球的商標(biāo)中心坐標(biāo)。

圖5 圖像處理

經(jīng)過一系列圖像處理操作之后可以找到乒乓球的商標(biāo)相對(duì)于乒乓球中心點(diǎn)的坐標(biāo),根據(jù)模型驗(yàn)證,選取連續(xù)5幀圖像的商標(biāo)中心坐標(biāo)為:(67,-25),(31,-58),(2,-64),(-30,-58),(-61,-35),代入橢圓一般方程,用Matlab計(jì)算可得:

(32)

利用上面模型可以算出中心坐標(biāo)為(-1.107,-2.645),此外夾角如下:

(33)

(34)

由于選擇坐標(biāo)系中心點(diǎn)盡量靠近坐標(biāo)軸,故理想中的β值為0ο,測量中可能會(huì)由于中心不準(zhǔn)有一定誤差,在可接受范圍內(nèi)。

根據(jù)實(shí)際測量數(shù)據(jù),攝像機(jī)鏡頭中心離地面距離112.5cm,球體中心離地面87.5cm,鏡頭中心離球體中心距離為62.5cm,根據(jù)實(shí)際數(shù)據(jù)可測得攝像機(jī)拍攝平面與球體平面夾角α計(jì)算公式如下:

(35)

可計(jì)算出:

(36)

至此,可以驗(yàn)證算法中關(guān)于求解乒乓球旋轉(zhuǎn)三維速度方向問題基本正確。

再利用上面角速度模型可以計(jì)算出角速度為:

(37)

角速度與實(shí)際測量角速度之間的誤差為:

(38)

即

(39)

角速度計(jì)算誤差控制在1.83%,基本可以驗(yàn)證本算法對(duì)于乒乓球旋轉(zhuǎn)速度計(jì)算的正確性。至此,在誤差允許的范圍內(nèi),已用對(duì)照實(shí)驗(yàn)驗(yàn)證了本算法的正確性。

以上是針對(duì)一般情況下的實(shí)驗(yàn)驗(yàn)證,特殊情況分為2種情況:① 當(dāng)M軸與攝像機(jī)的主光軸是平行的,這時(shí)形成的運(yùn)動(dòng)軌跡將不是一個(gè)橢圓,此時(shí)算出的軌跡方程將會(huì)是一個(gè)圓的軌跡方程;再利用公式可以算出此時(shí)L線與攝像機(jī)光軸的角度為90°,即垂直。② 當(dāng)商標(biāo)處于旋轉(zhuǎn)軸上時(shí),商標(biāo)相對(duì)于球體中心的坐標(biāo)是不變的,此時(shí)獲得的5個(gè)x、y坐標(biāo)是一樣的,代入方程組中有無窮多解,這種情況下無法獲得軌跡方程。這是最特殊的情況,此時(shí)L線與攝像機(jī)光軸線是平行的,即夾角為0°,這應(yīng)特別加以說明。

3 結(jié)論及展望

本文利用乒乓球商標(biāo)的運(yùn)動(dòng)軌跡參數(shù)方程,根據(jù)一些空間轉(zhuǎn)換關(guān)系獲得乒乓球旋轉(zhuǎn)的三維速度及大小,并用實(shí)驗(yàn)驗(yàn)證了算法的有效性。在后續(xù)研究中,主要從以下2個(gè)方面進(jìn)行深入:① 算法中軌跡點(diǎn)的數(shù)量要減少,盡量嘗試用更少的點(diǎn)進(jìn)行旋轉(zhuǎn)速度計(jì)算,從而降低對(duì)攝像機(jī)拍攝視頻的要求;② 在圖像處理算法上深入學(xué)習(xí)研究,爭取提高圖像識(shí)別精度,減少處理時(shí)間,滿足機(jī)器人視覺系統(tǒng)的實(shí)時(shí)性要求。

[1] Zhang Z,Xu D,Yu J.Research and latest development of ping-pong robot player[C]∥Proceedings of the 7thWorld Congress on Intelligent Control and Automation.Chongqing:IEEE,2008:4881-4886

[2] Andersson R L.A robot ping-pong player:Experiment in real time intelligent control[M].Cambridge,CA:MIT Press,1988:11-19

[3] Andersson R L.Aggressive trajectory generator for a robot ping-pong player[J].IEEE Control System Magazine,1989,9(2):15-21

[4] Acosta L,Rodrigo J,Mendez J A,et al.Ping-pong player prototype,a PC-based,low-cost,ping pong robot[J].IEEE Robotics and Automation Magazine,2003,10(4):44-52

[5] Matsushima M,Hashimoto T,Trkeuchi M,et al.A learning approach to robotic table tennis robotics[J].IEEE Transactions on Robotics and Automation,2005,21(4):767-771

[6] Matsushima M,Hashimoto T,Miyazaki F.Learning to the robot table tennis task-ball control and rally with a human[C]∥Proceedings of IEEE Conference on System,Man and Cybernetics.Washington D C,USA:IEEE,2003:2962-2969

[7] Miyazaki F,Matsushima M,Takeuchi M.Learning todynamically manipulate:a table tennis robot controls a ball and rallies with a human being[M].Advance in Robot Control.Berlin,Germany:Springer Berlin Heidelberg,2006:317-341

[8] Zhang Zhengtao,Xu De,Yu Junzhi.Visual measurement and prediction of ball trajectory for table tennis robot[J].IEEE Transactions on Instrumentation and Measurement,2010,59(12):3195-3205

[9] Leung H,Wang H Q,Cao C X.An improved locally weighted regression for a converter re-vanadium prediction modeling[C]∥The 6thWorld Congress on Intelligent Control and Automation.Dalian:[s.n.],2006:1515-1519

[10] Huang Yanlong,Xu De,Tan Min,et al.Trajectory prediction of spinning ball for ping-pong player robot[C]∥2011 IEEE/RSJ Int Conf on Intelligent Robots and Systems.San Francisco:[s.n.],2011:3434-3439

[11] Chen Xiaopeng,Tian Ye,Huang Qiang,et al.Dynamic model based ball trajectory prediction for a robot ping-pong player[C]∥The 2010 IEEE Int Conf on Robotics and Biomimetics.Tianjin:[s.n.],2010:603-608

[12] Akira N,Yosuke K,Yuki O,et al.Modeling of rebound phenomenon between ball and racket rubber[C]∥ICROSSICE Int Joint Conf.Fukuoka:[s.n.],2009:2295-2300

[13] Nakashima A,Yuki Ogawa.Robotic table tennis based on physical models of aerodynamics and rebounds[C]∥Conference:Robotics and Biomimetics(ROBIO),IEEE International Conference,2011:2348-2354

[14] Nakashima A,Tsuda Y,Liu C,et al.A real-time measuring method of translational/rotational velocity of a table tennis ball[C]∥Proceeding of 5thIFAC Symposium on Mechatronic Systems,2010:732-738

[15] Liu C,Hayakawa Y,Nakashima A.A registration algorithm for on-line measuring the rotational velocity of a table tennis ball[C]∥Proceeding of IEEE/RSJ International Conference on Intelligent Robots and Systems,2011:428-438

3D Velocity Measurement of Single-color Table Tennis Rotation Based on Visual System of Table Tennis Robot

JI Yunfeng, LU Aifa, REN Jie, SHI Zhihao

Rotation is the key winning factor to table tennis, thus acquiring the rotating information of the table tennis ball is an essential task of the table tennis robot vision. This paper presents the algorithm based on the monocular vision that can calculate 3D velocity of high-speed rotating sphere, which can be used to calculate the rotating velocity of table tennis ball with no marks. The algorithm can be applied to the vision recognition of table tennis robot to acquire the rotating information of table tennis ball, which plays an important role in motion planning and control decision of robot at the subsequent stages. The experimental results of the proposed algorithm are verified with little error, which proves the feasibility of the proposed method.

table tennis; visual system; 3D velocity; rotation; robot

2016-11-16；

2017-01-03

上海市科委科研計(jì)劃項(xiàng)目(15490503100);國家體育總局科研課題(2014B072)

季云峰(1990-),男,安徽合肥人,上海體育學(xué)院博士研究生;Tel.:15021307809,E-mail:ji_yunfeng2015@163.com 通信作者簡介:陸愛發(fā)(1958-),男,上海人,上海體育學(xué)院副教授;Tel.:13816860417,E-mail:luaifa2009@163.com

G846

A

1000-5498(2017)03-0083-06

DOI 10.16099/j.sus.2017.03.013