基于對偶圓錐曲線的橢圓中心精確定位方法

辛 睿， 武栓虎， 李愛娟

(煙臺大學(xué) 計算機(jī)與控制工程學(xué)院， 山東 煙臺 264005)

0 引 言

基于圓形模式的平面標(biāo)靶是近十多年來普遍采用的標(biāo)定工具，由于與棋盤模式平面標(biāo)靶[1]相比對圖像質(zhì)量要求較低[2]，在工業(yè)上具有較好的應(yīng)用價值，廣泛應(yīng)用于三維重建、動態(tài)跟蹤、機(jī)器人制導(dǎo)等方面。

在圓模式平面標(biāo)靶標(biāo)定過程中，由于存在鏡頭畸變和透視偏差，圓會投影成一個橢圓[2]。橢圓可以通過圖像強(qiáng)度或其梯度進(jìn)行估計，但是基于強(qiáng)度的方法[3]對包括光源、反射等非均勻照明非常敏感，而且算出的是準(zhǔn)確的橢圓中心，不是真實的圓心投影。與之相比基于梯度的方法受不均勻照明的影響較小，但往往需要精確提取亞像素級橢圓輪廓邊緣點(diǎn)，再將橢圓參數(shù)擬合到獲得的點(diǎn)集[4-6]。由于該方法依賴于一組事先提取的輪廓點(diǎn)，輪廓點(diǎn)的異常值會直接影響后續(xù)結(jié)果。而且在低對比度圖像中，有時會遇到很難沿橢圓提取每一個輪廓點(diǎn)的情況，會嚴(yán)重影響橢圓中心的提取。

針對以上問題，本文在研究圓模式識別方法的基礎(chǔ)上，提出一種直接利用橢圓輪廓附近的梯度矢量場精確獲取對偶圓錐曲線的橢圓中心的定位方法。該方法直接從包含特征的區(qū)域中計算出的圖像梯度估算橢圓中心特征的位置，避免了亞像素精度的橢圓輪廓邊緣點(diǎn)的提取，簡化了橢圓參數(shù)獲取過程并可以進(jìn)一步使用原始圖像數(shù)據(jù)來評估恢復(fù)的橢圓參數(shù)的不確定性。此外，通過模擬實驗和對比實驗評估該方法獲得橢圓中心的精度、對環(huán)境噪聲的抵抗性，證明該方法是可以用在相機(jī)標(biāo)定應(yīng)用中的高精度特征定位方法。

1 基于對偶圓錐曲線的特征點(diǎn)估計方法

1.1 對偶圓錐曲線基本性質(zhì)

橢圓輪廓通常由二次曲線進(jìn)行描述，在歐氏幾何中，二次曲線有3種主要類型：雙曲線、橢圓和拋物線。在經(jīng)典理論中，這3類曲線是不同方向的平面截二次錐面相得到的曲線，也稱為圓錐曲線[7]。在非齊次坐標(biāo)中，設(shè)p=[x,y]T為圓錐曲線上的點(diǎn)，則圓錐曲線方程為：

ax2+bxy+cy2+dx+ey+f=0,

(1)

其中，α=[a,b,c,d,e,f]T為圓錐曲線的系數(shù)，該方程可表示為矩陣形式，圓錐曲線系數(shù)矩陣為C：

pTCp=0,

(2)

(3)

設(shè)線l=[α,β,γ]T，由射影幾何中對偶原理，齊次點(diǎn)和線的基本結(jié)合方程是對稱的。則點(diǎn)圓錐曲線C有一個線圓錐曲線C*與之對應(yīng)，稱為對偶圓錐曲線、也稱作圓錐曲線的包絡(luò)[7](如圖1所示)，圓錐曲線C*的所有切線l均滿足：

lTC*l=0.

(4)

其中，對偶圓錐曲線系數(shù)矩陣C*為：

(5)

(a) 滿足pTCp=0的點(diǎn)p在一條點(diǎn)圓錐曲線上

(b) 滿足lTC*l=0的直線l是點(diǎn)圓錐曲線的切線

(b) The tangent setlof the point conic curve satisfieslTC*l=0

圖1 圓錐曲線與對偶圓錐曲線

Fig. 1 Conic curve vs dual conic curve

值得一提的是，設(shè)線l的極點(diǎn)為p0，則線l的極點(diǎn)可以由p0=C*l獲得。設(shè)無窮遠(yuǎn)線l的極點(diǎn)為p，那么對偶圓錐曲線的中心ec=(uc,vc)T，就是對偶圓錐曲線無窮遠(yuǎn)線的極點(diǎn)p?？梢杂蓀=C*l獲得[8]，具體公式如下：

(6)

其中，ξ為比例系數(shù)，l∞為無窮遠(yuǎn)線。由該公式可知，只要求得對偶圓錐曲線參數(shù)就可以直接獲得標(biāo)志點(diǎn)中心坐標(biāo)。對偶圓錐曲線參數(shù)的求取首先需要獲得一組切線方程，然后通過給定一組切線集合利用最小二乘法[9]進(jìn)行對偶圓錐曲線參數(shù)的估計。接下來對切線方程組的獲取進(jìn)行說明。

1.2 基于梯度圖像的切線估計方法

求取對偶圓錐曲線首先需要獲取切線方程組，切線方程組可以通過求取給定圖像的梯度獲得。實際上，特征邊緣的相鄰區(qū)域內(nèi)的像素處的梯度提供了穿過像素中心的法線方向。根據(jù)法線和梯度的性質(zhì)：法線方向垂直于切線方向，可以根據(jù)法線求得切線；梯度強(qiáng)的地方(特征邊緣)才存在切線，梯度為零的地方切線不存在。在法線不存在但切線存在時，設(shè)切線斜率為0，再根據(jù)具體像素點(diǎn)確定切線方程。

研究中，將篩選出包含橢圓特征的區(qū)域，從包含特征的區(qū)域中通過Scharr濾波器[10]的方式計算出圖像梯度。設(shè)點(diǎn)p=[x,y]T的梯度為G=[Gx,Gy]T，圖2(a)的符號表示一個3×3區(qū)域內(nèi)圖像點(diǎn)的灰度。通過使用圖2(a)中的3×3模板以及圖2(b)、圖2(c)中的Scharr算子[10]，以z5為中心的一個3×3鄰域?qū)x和Gy的近似如下式所示：

Gx=(3z7+10z8+3z9)-(3z1+10z2+3z3)，

(7)

Gy=(3z3+10z6+3z9)-(3z1+10z4+3z7),

(8)

通過Scharr濾波器可以快速得到橢圓輪廓附近的圖像梯度，也可以選擇其他濾波器進(jìn)行圖像梯度的獲取。當(dāng)p=[x,y]T的圖像梯度G=[Gx,Gy]T不為零時，可以定義通過像素中心的線的法線方向(法線方向垂直于切線方向)，得到切線方程l：

l=λ[Gx,Gy,-GTp]T.

(9)

其中，λ為任意比例系數(shù)。通過求取橢圓輪廓附近所有的圖像梯度，可以獲得該橢圓圖像的切線集合。

1.3 對偶圓錐曲線參數(shù)估計方法

通過求取的一組切線集合li,i為切線的下標(biāo)，i≥5。對偶圓錐曲線參數(shù)可以通過最小二乘法[9]進(jìn)行估計。設(shè)切線li=[xi,yi,ki]，用公式(10)作為近似擬合函數(shù)，在對偶圓錐曲線中加入約束A=1，則：

圖2 Scharr濾波器濾波系數(shù)

(10)

(11)

(12)

可由方程(12)解得:

ω=[B*,C*,D*,E*,F*]T.

由于D*,E*,F*與D,E,F只相差一個比例系數(shù)A。對公式(9)結(jié)果不產(chǎn)生影響，可由公式(6)求得對偶圓錐曲線中心。

1.4 算法步驟

本文提出的基于對偶圓錐曲線的橢圓中心精確定位方法的具體步驟如下：

(1)對圖像進(jìn)行預(yù)處理?；叶然?，均值濾波，再對圖像進(jìn)行二值化。

(2)通過1.2節(jié)方法提取二值圖像梯度圖像，篩選出包含橢圓特征的區(qū)域。

(3)利用梯度圖像上保留的輪廓法線信息通過1.2節(jié)方法獲得橢圓邊緣的切線集合。

(4)通過1.3節(jié)方法提取出每一幅圖片中的所有對偶圓錐曲線中心。

2 實驗與結(jié)果

實驗的目的是要驗證,無需亞像素精度的標(biāo)志點(diǎn)輪廓邊緣點(diǎn)的提取,通過本文提出的對偶圓錐曲線方法就可獲得較為精確的標(biāo)靶圓心對應(yīng)特征點(diǎn)。進(jìn)行了模擬噪聲實驗，驗證該方法受環(huán)境噪聲影響較小。

圖3 實驗圖像

在本實驗中，為了檢測本文方法的準(zhǔn)確性，在不添加任何噪聲的情況下，使用本文方法對200幅圖像進(jìn)行橢圓真實中心的識別，并與真值橢圓中心進(jìn)行比較算出橢圓中心的相對誤差。同時也對質(zhì)心法，點(diǎn)心法進(jìn)行了實驗對比，實驗結(jié)果見表1。由表1可知，本文方法取得了較為理想的效果。接下來，圖像被加性高斯噪聲污染，標(biāo)準(zhǔn)偏差從σn=2%，增加到圖像動態(tài)范圍的5%。對于每一個噪聲水平，在200個不同的單個橢圓圖像上評估橢圓位置的平均誤差。

表1 實驗結(jié)果對比

Tab. 1 Comparison of experimental resultspixel

實驗結(jié)果表明，在不添加任何噪聲的情況下，質(zhì)心法和點(diǎn)心法直接獲得的是橢圓中心，而對偶圓錐曲線方法更接近于真實的中心投影。但是由于質(zhì)心法在特征的所有像素上取平均值，因此具有很強(qiáng)的抗噪能力，而對偶圓錐曲線的方法隨著噪聲的增加，對梯度方向的估計變得不太精確，在σn=5%與點(diǎn)心法有了相似的精度，可以通過改進(jìn)梯度值的估計來提高本文方法的精度。最后，表1證實了這3種方法精度都在0.1pixel左右，具有較高的精度和抗噪性能。

3 結(jié)束語

特征點(diǎn)精確估計是相機(jī)標(biāo)定過程中非常重要的部分。圓模式平面標(biāo)靶是一常用的相機(jī)標(biāo)定工具，相對棋盤標(biāo)靶，具有制造容易、特征點(diǎn)恢復(fù)受光照和噪聲的影響較小的特點(diǎn)。但由于環(huán)境噪聲和像素化的影響，直接利用提取的橢圓輪廓往往會產(chǎn)生一些誤差，從而降低了標(biāo)定精度。本文提出了一種基于對偶圓錐曲線獲取標(biāo)志點(diǎn)中心的標(biāo)定方法，直接利用標(biāo)志點(diǎn)邊界附近的原始梯度信息，避免了不精確邊緣點(diǎn)集的提取，在對抗噪聲方面具有一定優(yōu)勢，并且直接應(yīng)用最小二乘法計算對偶圓錐曲線中心，可以簡化獲取橢圓真實中心的步驟的同時避免了部分噪聲點(diǎn)的干擾，具有較好的應(yīng)用價值。