摘要：巡線機(jī)器人沿相線行走時(shí)必須探測(cè)識(shí)別各種障礙，并根據(jù)障礙類型規(guī)劃越障行為 。針對(duì)220 kV架空輸電線路的結(jié)構(gòu)特點(diǎn)，利用視覺傳感器，設(shè)計(jì)了基于結(jié)構(gòu)約束的障礙識(shí)別算法．算法利用圖像的邊緣信息，采用改進(jìn)的基于存在概率圖的圓/橢圓檢測(cè)方法和分層決策機(jī)制，以減少自然環(huán)境中光線變化和機(jī)器人運(yùn)動(dòng)對(duì)識(shí)別質(zhì)量的影響，滿足了巡線機(jī)器人的實(shí)時(shí)越障要求。實(shí)驗(yàn)室模擬線路和實(shí)際線路實(shí)驗(yàn)結(jié)果表明，算法能可靠地識(shí)別出復(fù)雜背景中的防震錘、懸垂線夾和耐張線夾等障礙物。 關(guān)鍵詞：障礙識(shí)別；圓檢測(cè)；巡線機(jī)器人；霍夫變換；輸電線路 Abstract：Power transmission line inspection robot must plan its behaviors to negotiate obstacles according to their typeswhen it is crawling along the power transmission line For this purpose．based on the structure of a 220 kV transmission line，a structure—constrained obstacle recognition algorithm is designed with machine vision sensors．The algorithm uses an improved existence-probability—map—based circle/ellipse detection method and a hierarchical decision mechanism to reduce the effects of illumination variation and robot motion on obstacle recognition quality，which satisfies the needs of real—time obstacle negotiation of inspection robot．The results of experiments with simulation and real transmission lines show that the algorithm can reliably recognize obstacles such as counterweight，strain clamp，and suspension clamp from cluttered back— ground． Keywords：obstacle recognition；circle detection；inspection robot；Hough transform；power transmission line 1引言（Introduction） 架空輸電線路巡線機(jī)器人能在線路帶電情況下沿相線或地線爬行，跨越防震錘、耐張線夾、懸垂線夾、桿塔等障礙，利用攜帶的檢測(cè)儀器代替人對(duì)桿塔、導(dǎo)線、避雷線、絕緣子、線路金具、線路通道等實(shí)施接近檢查。這種作業(yè)方式將大大減輕輸電線路維護(hù)人員的勞動(dòng)強(qiáng)度，提高巡檢效率與質(zhì)量，對(duì)保證輸電系統(tǒng)的安全可靠運(yùn)行具有重要意義。 巡線機(jī)器人多采用輪式行走機(jī)構(gòu)，架空輸電線路相線和地線上安裝的防震錘、絕緣子、懸垂線夾、耐張線夾等桿塔支撐附件便成為影響巡線機(jī)器人行走的障礙。因此，巡線機(jī)器人必須借助傳感器探測(cè)并識(shí)別這些障礙類型，然后根據(jù)障礙類型規(guī)劃越障行為，才能跨越桿塔，沿輸電線路大范圍自主行走。常用的傳感器有接近傳感器（紅外傳感器、超聲傳感器）和視覺傳感器（攝像機(jī)），Peters等[1]在一臺(tái)巡線機(jī)器人上配置了由多達(dá)34個(gè)接近傳感器組成的傳感陣列，采用基于粗糙集的神經(jīng)計(jì)算方法設(shè)計(jì)了障礙分類器，根據(jù)傳感陣列的測(cè)量數(shù)據(jù)進(jìn)行障礙識(shí)別分類，進(jìn)而規(guī)劃?rùn)C(jī)器人的越障行為。Peungsungwal等利用攝像機(jī)識(shí)別架空線路障礙，完成了實(shí)驗(yàn)線路兩端絕緣子的識(shí)別。 本文根據(jù)220kV架空輸電線路的結(jié)構(gòu)特點(diǎn)，利用視覺傳感器，采用改進(jìn)的幾何圖形檢測(cè)方法和分層決策機(jī)制，設(shè)計(jì)了基于結(jié)構(gòu)約束的障礙檢測(cè)和識(shí)別算法。 2 220kV巡線機(jī)器人及運(yùn)行環(huán)境分析（220kVpowertransmissionlineinspectionro—botanditsworkingenvironmentanaly-sis） 巡線機(jī)器人采用雙臂輪結(jié)構(gòu)，沿220kV相線滾功行走。圖1是220kV輸電線路的結(jié)構(gòu)示意圖，巡線機(jī)器人需要跨越的典型障礙有3種：（1）防震錘；（2）懸垂線夾和單串絕緣子；（3）耐張線夾、雙串絕緣子及引流跳線。防震錘用夾扣固定在懸垂線夾、耐張線夾的兩側(cè)，是巡線機(jī)器人接近桿塔時(shí)首先要跨越的障礙。上述3種典型障礙具有不同的結(jié)構(gòu)組成和空間布局，巡線機(jī)器人在越障時(shí)需采用相應(yīng)的越障策略。 攝像機(jī)安裝在巡線機(jī)器人的本體前端，光軸與機(jī)器人前進(jìn)方向的夾角約為30[sup]。[/sup]，如圖1所示。在此配置下，背景是天空，這樣可以降低背景的復(fù)雜度，而防震錘端面、懸垂線夾支撐絕緣子、耐張線夾支撐絕緣子在攝像機(jī)成像平面上的投影呈直線、圓、橢圓等典型幾何圖形。對(duì)攝像機(jī)攝取的環(huán)境圖像，經(jīng)圖像分析和理解，識(shí)別圖像中是否含有特定的障礙目標(biāo)。 3障礙識(shí)別算法原理及實(shí)現(xiàn)（Principleoftheobstaclerecognitionalgorithmanditsimplementation） 防震錘、懸垂線夾、耐張線夾和絕緣子串既無(wú)豐富的表面紋理，也無(wú)鮮明的顏色特征，且相互鉸接，因而難以分割成單獨(dú)的區(qū)域。針對(duì)障礙物的外觀、形狀和結(jié)構(gòu)特點(diǎn)，我們選擇相對(duì)簡(jiǎn)單的直線、圓（圓?。?、橢圓等幾何圖形基元作為障礙物存在的線索，結(jié)合輸電線路結(jié)構(gòu)對(duì)上述圖形基元位置的約束來(lái)識(shí)別不同類型的障礙物。由于3類障礙物之間存在較大差異，因此，我們把3類障礙物的識(shí)別問題簡(jiǎn)化為3個(gè)分別把防震錘、懸垂線夾、耐張線夾及絕緣子串從背景中分辨出來(lái)的兩類分類問題，從而簡(jiǎn)化分類器的設(shè)計(jì)。算法的兩個(gè)關(guān)鍵問題是：（1）從邊緣圖像中快速可靠地抽取直線、圓（圓?。┖蜋E圓等圖形基元；（2）在結(jié)構(gòu)約束下判斷檢測(cè)到的圖形基元源自于何種障礙。 3．1相線的識(shí)別與定位 防震錘、懸垂線夾、耐張線夾及絕緣子串以不同方式固定在相線上，因此相線的識(shí)別、定位和跟蹤對(duì)障礙的識(shí)別至關(guān)重要。首先，用高斯濾波器消除高頻噪音，并對(duì)圖像適度模糊化，以減少局部紋理邊緣，提高導(dǎo)線邊緣直線的檢出率；然后，采用Canny算子提取圖像的邊緣圖，圖像大小為320×240；最后，用漸進(jìn)概率霍夫變換PPHT（ProgressiveProbabilisticHoughTransfoITn）[3]抽取邊緣圖中長(zhǎng)度大于給定閾值T的所有直線段，每一直線L[sub]i[/sub]用以下參數(shù)表示：θ，傾角，直線與軸正方向夾角，如圖2所示；（x[sup]1[/sup]，y[sup]1[/sup]），上端點(diǎn)坐標(biāo)；（x[sub]2[/sub]，Y[sub]2[/sub]），下端點(diǎn)坐標(biāo)。 巡線機(jī)器人沿相線行走時(shí)，相線在圖像平面中的投影呈棒狀從圖像平面上邊框傾斜向下。機(jī)器人的抖動(dòng)對(duì)相線在圖像平面中軸方向的投影位置及傾角有一定影響，實(shí)驗(yàn)表明，相線投影位置及傾角的變化在確定的較小范圍內(nèi)。因此，我們采用基于規(guī)則的分類器對(duì)檢出的直線段分類，找出相線對(duì)應(yīng)的兩條邊緣直線L[sub]L[/sub]、L[sub]R[/sub]。設(shè)E=｛L[sub]i[/sub]，i=1，2，…，n｝表示由PPHT檢出的直線集合，C表示相線邊緣直線候選集合，分類器分類規(guī)則為： 如果任意直線段L[sub]i[/sub]∈E滿足條件：（1）傾角θ滿足70[sup]。[/sup]<θ<110[sup]。[/sup]；（2）上端點(diǎn)坐標(biāo)Y軸分量y[sup]1[/sup]x[sub]1i[/sub]，i≠j）；T[sub]1[/sub]、T[sub]2[/sub]根據(jù)相線投影寬度來(lái)確定，則L[sub]i[/sub]、L[sub]j[/sub]分別為相線的左、右邊緣直線L[sub]L[/sub]、L[sub]R[/sub]。 這樣就完成了相線的識(shí)別和定位，大致確定了障礙所在的區(qū)域，為后續(xù)的障礙識(shí)別奠定了基礎(chǔ)。

3．2防震錘的識(shí)別 防震錘是巡線機(jī)器人接近桿塔時(shí)遇到的第一個(gè)障礙，它由兩個(gè)中空的圓筒連接而成，用夾扣固定在相線下部，如圖6所示。防震錘的識(shí)別基于兩個(gè)假設(shè)： 假設(shè)1 防震錘兩個(gè)可見端面各有一個(gè)孔孔洞，前大后小，投影到圖像中為局部暗區(qū)域，前端面孔洞面積較小，呈圓形；后端面孔洞面積較大，呈半月形。因此，可采用自適應(yīng)閾值分割方法，分割出上述兩個(gè)孔洞區(qū)域，作為防震錘存在的證據(jù)之一。 假設(shè)2 防震錘前端面邊緣呈現(xiàn)以孔洞中心為圓心的不規(guī)則上半圓弧，該半圓弧作為防震錘存在的證據(jù)之二。 在上述假設(shè)的基礎(chǔ)上，我們?cè)O(shè)計(jì)了防震錘識(shí)別算法，步驟如下： （1）低層圖像處理：根據(jù)已識(shí)別出的相線位置和寬度，確定含有防震錘的圖像窗口W[sub]r[/sub]。采用 Otsu[4]最優(yōu)閾值算法對(duì)圖像窗口W[sub]r[/sub]進(jìn)行二值化分割，并對(duì)得到的二值圖像進(jìn)行形態(tài)學(xué)開運(yùn)算濾波以消除圖像噪點(diǎn)，然后提取各區(qū)域質(zhì)心（x[sub]c[/sub]，y[sub]c[/sub]）、面積和輪廓等區(qū)域?qū)傩浴? （2）孔洞候選區(qū)域確定：由相線和防震錘的位置約束，確定代表防震錘存在的孔洞候選區(qū)域。 設(shè) θ[sub]L[/sub] 、θ[sub]R[/sub]（x[sub]1L[/sub]，y[sub]1L[/sub]）與 （x[sub]1R[/sub]，y[sub]1R[/sub]）分別表示相線左、右邊緣直線的傾角和上端點(diǎn)坐標(biāo)；Q是步驟（1）中分割出的任一區(qū)域，質(zhì)心坐標(biāo)為（x[sub]c[/sub]，y[sub]c[/sub]），則質(zhì)心到相線左、右邊緣直線的水平距離和質(zhì)心所在位置處相線寬度分別為： 如果滿足式（2），則區(qū)域 Q為端面孑L洞候選區(qū)域。 （3）端面孔洞區(qū)域確定：對(duì)候選區(qū)域施加端面幾何形狀約束。后端面暗孔洞面積較大，且分割出的區(qū)域邊界輪廓的下部應(yīng)呈現(xiàn)為半圓弧，因此對(duì)候選區(qū)域中最大面積區(qū)域的輪廓進(jìn)行半圓弧檢測(cè)，如果該區(qū)域下部輪廓為半圓弧，則進(jìn)一步支持該區(qū)域?yàn)楹蠖嗣婵锥磪^(qū)域的假設(shè)。 前端面為圓形，在邊緣圖中可提取與端面孔洞同心的上半圓弧斷續(xù)邊緣。因此，對(duì)于剩下候選區(qū)域，在邊緣圖中檢測(cè)是否存在以該候選區(qū)域質(zhì)心為圓心的上半圓弧。如存在，則進(jìn)一步確定該區(qū)域?yàn)榍岸嗣婵锥磪^(qū)域。 （4）如果關(guān)于前、后端面孔洞區(qū)域的兩條假設(shè)都成立，則確認(rèn)檢測(cè)到防震錘，同時(shí)估計(jì)出前、后端面圓參數(shù)。 如何檢測(cè)不規(guī)則邊緣在感知意義上是否為圓或半圓弧，是實(shí)現(xiàn)防震錘識(shí)別算法的關(guān)鍵，具體原理將在3.4節(jié)中詳細(xì)介紹。 3.3懸垂線夾識(shí)別 懸垂線夾把相線夾緊，掛在單絕緣子串下部，把相線架在空中。其本身體積小，特征難以提取。在圖1所示攝像機(jī)配置下，單絕緣子串的底部在圖像中呈橢圓形，主軸與軸基本平行，并且這種構(gòu)型在3種障礙中具有唯一性，這樣我們就可以把懸垂線夾識(shí)別問題轉(zhuǎn)化為橢圓檢測(cè)問題。文[5]給出了一種有效的橢圓檢測(cè)方法，計(jì)算量較大。在此基礎(chǔ)上，我們根據(jù)絕緣子成像特點(diǎn)對(duì)文[5]中的橢圓檢測(cè)算法進(jìn)行了改進(jìn)。 描述任意橢圓需5個(gè)未知參數(shù)，即中心坐標(biāo)（x[sub]0[/sub]， y[sub]0[/sub]）、長(zhǎng)半軸a、短半軸b和旋轉(zhuǎn)角α ，相應(yīng)的參數(shù)方程為： 如圖3所示，我們假定邊緣點(diǎn)對(duì)（x[sub]1[/sub]， y[sub]1[/sub]）、（x[sub]2[/sub]， y[sub]2[/sub] ）是待檢橢圓主軸的兩個(gè)頂點(diǎn)，則待檢橢圓的4個(gè)參數(shù)可按式（4）估計(jì) ： 設(shè)（ x，y）是待檢橢圓上的一個(gè)點(diǎn)，短半軸 b可用式（5）計(jì)算： 接下來(lái)，我們采用類似霍夫變換（HT）的投票法估計(jì)待檢橢圓的短半軸b。設(shè)A（i）是一個(gè)一維累加器數(shù)組，（min_b≤i≤max_b，min_b和max_b分別是預(yù)設(shè)b的最小、最大值）。W[sub]s[/sub]是邊緣圖像中包含待檢橢圓的最小矩形區(qū)域。對(duì)于中的每一邊緣點(diǎn)（x[sub]i[/sub]，y[sub]i[/sub]），用式（5）計(jì)算出對(duì)應(yīng)的短半軸b[sub]i[/sub]并取整，若滿足min_b≤b[sub]i[/sub]≤max_b，則累加器數(shù)組單元A（b[sub]i[/sub]）加1，即邊緣點(diǎn)（x[sub]i[/sub]，y[sub]i[/sub]）為待檢橢圓短半軸是bi投一票。投票結(jié)束后，統(tǒng)計(jì)累加器數(shù)組最大值，設(shè)為A（b[sub]m[/sub]），及對(duì)應(yīng)的短半軸b[sub]m[/sub]，那么圖像中存在由估計(jì)出的5個(gè)參數(shù)所描述的橢圓的概率為： 對(duì) P取閾值，如果 P>T[sub]e[/sub]，則圖像中存在由估計(jì)出的5個(gè)參數(shù)所描述的橢圓，否則不存在；T[sub]e[/sub]的大小視抽取的橢圓邊緣的質(zhì)量通過(guò)試驗(yàn)確定。

基于橢圓檢測(cè)的懸垂線夾識(shí)別步驟為：（1）提取邊緣圖像，盡量消除局部紋理形成的瑣碎邊緣點(diǎn)，以提高橢圓檢出率。 （2）沿相線由上至下搜索橢圓，并根據(jù)相線和懸垂線夾的位置約束及絕緣子投影橢圓可能的大小，確定搜索區(qū)域的寬度，如圖4所示。搜索時(shí)，找出在長(zhǎng)條矩形框中最左、最右的邊緣點(diǎn)，分別記為（x[sub]1[/sub]，y[sub]1[/sub]）、（x[sub]2[/sub]，y[sub]2[/sub]），然后按式（4）計(jì)算出假設(shè)橢圓的4個(gè)參數(shù)，并確定包圍該橢圓的最小矩形窗W[sup]s[/sup]，遍歷W[sup]s[/sup]內(nèi)所有邊緣點(diǎn)，按式（5）、（6）估計(jì)假設(shè)橢圓的短半軸b和存在概率P，按順序保存橢圓參數(shù)和相應(yīng)存在概率，｛（x[sub]0[/sub]，y[sub]0[/sub]），a，b，α，P｝。 （3）統(tǒng)計(jì)最大存在概率P[sub]m[/sub]，如果P[sub]m[/sub]>T[sub]e[/sub]，則表明從圖像中識(shí)別出懸垂線夾，位置由P對(duì)應(yīng)的橢圓參數(shù)給出。

3．4 耐張線夾識(shí)別 支撐耐張線夾的兩串絕緣子底部在圖像中呈現(xiàn)被相線左、右分開的兩個(gè)圓。類似懸垂線夾的識(shí)別 ， 我們把耐張線夾的識(shí)別問題轉(zhuǎn)化為圓的檢測(cè)。任意圓可由圓心坐標(biāo)（x[sub]0[/sub]，y[sub]0[/sub]）和半徑 r三個(gè)參數(shù)描述，即： 霍夫變換 （HT）被廣泛應(yīng)用于幾何圖形檢測(cè) ， 文[6—8]提出了幾種基于霍夫變換的圓檢測(cè)方法。 針對(duì)絕緣子底部圓形邊緣不完整、不規(guī)則的情況，本文提出了一種基于存在概率圖（existence probability map）的圓檢測(cè)方法，能夠有效地將邊緣不清晰圓檢測(cè)出來(lái)。 設(shè)f[sub]e[/sub]（x，y）表示由Canny算子提取的邊緣圖像；P（x，y）為二維數(shù)組，保存圓存在概率圖，單元（x，y）的值表示邊緣圖像f[sub]e[/sub]（x，y）中存在以點(diǎn)（x，y）為圓心、半徑為Rx,y的圓的概率，R[sub]x,y[/sub]=R（x，y）；R（x，y）為二維數(shù)組，單元（x，y）保存邊緣圖像（x，y）中點(diǎn)（x，y）處可能存在的圓的半徑；Rmin表示待檢測(cè)圓的最小允許半徑，視具體問題而定；A（i）是一個(gè)一維累加器數(shù)組，算法原理如下： （1）選擇圓心（u，v）∈f[sub]e[/sub]（x，y），如圖5所示，確定位于該點(diǎn)的圓的可能最大半徑Rmax，計(jì)算包圍以（u，v）為圓心、Rmax為半徑的圓的最小矩形區(qū)域 （2）對(duì)于W[sub]r[/sub]中的每一邊緣點(diǎn)（x[sub]i[/sub]，y[sub]i[/sub]），計(jì)算該點(diǎn)到圓心（u，v）的歐氏距離D，如果R[sub]min[/sub] 將概率值P（u，v）和半徑。分別保存到圓存在概率圖P（x，y）和R（x，y）對(duì)應(yīng)點(diǎn)。 （4）對(duì)于邊緣圖像（x，y）中的所有點(diǎn)（x，y），重復(fù)步驟（1）～（3），最終生成圓存在概率圖及對(duì)應(yīng)的半徑圖。 （5）分析圓存在概率圖P（x，y）的峰值，確定檢出圓的數(shù)量及圓參數(shù)。圓存在概率圖中的每一個(gè)峰值點(diǎn)表示該處可能存在一個(gè)圓，一種峰值檢出方法是：設(shè)定圖中圓的數(shù)量為N，概率閾值為T[sub]p[/sub]，先找出P（x，y）的最大值： 如果P[sub]max[/sub]（x[sub]1[/sub]，y[sub]1[/sub]）>Tp，則認(rèn)為圖像f[sub]e[/sub]（x，y）中有以點(diǎn)（x[sub]1[/sub]，y[sub]1[/sub]）為圓心、半徑為R（x[sub]1[/sub]，y[sub]1[/sub]）的圓，然后，令圓存在概率圖P（x，y）中點(diǎn)（x[sub]1[/sub]，y[sub]1[/sub]）及其鄰域的概率值為0，鄰域的大小視邊緣質(zhì)量和圖中相鄰圓的距離而定。 重復(fù)上述峰值檢出過(guò)程，直至找夠Ⅳ個(gè)圓，或剩下峰值小于給定閾值后停止。

在識(shí)別耐張線夾時(shí)，首先使用上述圓檢測(cè)方法檢出圖像中所有可能的圓，然后對(duì)檢出的圓實(shí)施結(jié)構(gòu)約束檢查，如果存在一對(duì)圓滿足下面兩個(gè)條件，則認(rèn)為已從圖像中識(shí)別出耐張線夾： （1）兩圓的半徑接近； （2）兩圓圓心約水平分布在相線左右兩側(cè)。 因算法要遍歷所有圖像點(diǎn)，計(jì)算量較大，考慮到耐張線夾的識(shí)別不需要精確估計(jì)圓參數(shù)，我們用最近鄰采樣法把原圖像從320×240縮小為80×60，來(lái)提高運(yùn)算速度，但這同時(shí)會(huì)造成圖像邊緣質(zhì)量下降影響圓檢測(cè)質(zhì)量。盡管如此，算法仍表現(xiàn)出較穩(wěn)定的性能。 4實(shí)驗(yàn)結(jié)果（Experimental results） 我們分別對(duì)實(shí)驗(yàn)室模擬線路和實(shí)際線路上的障礙進(jìn)行識(shí)別，圖6（1）是模擬線路防震錘、懸垂線夾耐張線夾的識(shí)別結(jié)果；圖6（2）是實(shí)際線路懸垂線夾耐張線夾的識(shí)別結(jié)果。從圖中可以看出，盡管因視角、雜亂背景和對(duì)比度較弱等因素造成提取的障礙物邊緣不規(guī)則完整，算法仍能檢出圖像中的圓/橢圓等圖形基元，再經(jīng)結(jié)構(gòu)約束下的檢驗(yàn)分類，最終將各種障礙物可靠地識(shí)別出來(lái)。 5結(jié)論（Conclusion） 從圖像中提取的邊緣點(diǎn)，在一定條件下包含了原圖像的全部信息，而興趣目標(biāo)也常常出現(xiàn)在圖像的邊緣處。因此，算法充分利用圖像邊緣信息，采用改進(jìn)的圓／橢圓檢測(cè)方法從中提取構(gòu)成障礙物的典型幾何圖形，并利用巡線機(jī)器人運(yùn)行環(huán)境的結(jié)構(gòu)約束，分層決策，縮小搜索空間，使算法簡(jiǎn)單有效?；诖嬖诟怕蕡D的圓／橢圓檢測(cè)方法，將邊緣點(diǎn)的“投票”結(jié)果變換為參數(shù)空間中圓／橢圓存在的相對(duì)可能性，提高了不清晰邊緣的圓／橢圓檢出率。本文設(shè)計(jì)的架空高壓輸電線路巡線機(jī)器人障礙視覺識(shí)別算法，基本滿足了巡線機(jī)器人越障要求。由于采用基于閾值的決策方法，算法雖然在較寬的范圍內(nèi)表現(xiàn)出穩(wěn)定一致的性能，但也存在對(duì)環(huán)境變化適應(yīng)能力不確定的問題。 參考文獻(xiàn) （References） [1] Peters J F，Ahn T C，Borkowskii M．Obstacle classification by a line—crawling robot：a rough neurocomputing approach[A]．Pro-ceedinss 0f the Third International Conference on Rough Sets and Current Tren（Is in Computing — Lecture Notes in Artificial Intelli• genee[C]．London，UK：Springer—Verlag，2002．594—601．Peungsungwal S，Pung～fi B，Chamnongthai K，et a1．Autonomous [2] robot for a power transmission line inspection[A]．Proceedings ofthe 2001 IEEE Iuternafional Symposium on Cimuits and Systems [C]．Piscataway，NJ，USA：IEEE，2001．121—124． [3] Galamhos C，Matas J，Kittler J．Progressive probabilistic Hough transform for line detection[A]．Proceedings of the 1999 IEEE Computer Society Conference on Computer Vision and Pattern Re— cognition[C]．Los Alamitos，CA，USA：IEEE Computer Society， 1999．554—560． [4]Otsu N，A threshold selection method from gray—level histogram[J]． IEEE Transactions on System， Man，an d Cybernetics， 1979，9 （1）：62—66． [5]Xie Y H，Ji Q，A new efficient ellipse detection method[A]，Pro— ceedings of the 16th International Conference on Pattern Recognition [C]．Los Alamitos，CA，USA：IEEE Computer Society，2002， 957—960． [6]Duda R 0，Hart P E．Use of the Hough transformation to detect lines and curves in pictures[J]．Communications of the ACM， 1972，15（1）：11一l5． [7]Kimme C，Ballard D H，Sklansky J．Finding circles by an array of accumulators[J]．Communicatim～s oftheACM，1975，18（2）：120 ～l22． [8] Leavem V F．Shape Detection in Computer Vision Using Hough Transform[M]．London：Springer-Verlag，1992．