胡銀風(fēng)，江 晟，宋 德

（長春理工大學(xué)物理學(xué)院，長春 130022）

0 引言

近年來，先進(jìn)制造技術(shù)的發(fā)展日新月異，高精度、快速、高質(zhì)量的工件檢測技術(shù)已成為發(fā)展先進(jìn)制造至關(guān)重要的一個環(huán)節(jié)，伴隨著電子、光學(xué)和計算機(jī)科學(xué)技術(shù)的不斷成熟，機(jī)器視覺在精密工件檢測技術(shù)中得到了飛速發(fā)展［1-2］。

在圖像處理過程中，邊緣檢測是最重要的環(huán)節(jié)之一，邊緣能最直接地反映物體輪廓信息，因此可以通過圖像邊緣信息計算工件的尺寸。邊緣提取常用的方法如運用一階差分的Krisch算子、Sobel 算子、Canny 算子等，以及二階微分算子Laplace［3-4］。為了能夠得到平滑且較為完整的圖像邊緣，羅朝陽等［5］提出利用數(shù)學(xué)方法中形態(tài)學(xué)處理方法提取圖像邊緣，但是該方法對噪聲很敏感且檢測過程中容易出現(xiàn)偽邊緣。包從望等［6］提出改進(jìn)Canny 邊緣檢測方法，運用3*3 有限差分方法并結(jié)合OTSU 的自適應(yīng)閾值提升算法的抗噪能力。曹義親等［7］提出將改進(jìn)后的LSD 直線檢測應(yīng)用到鋼軌表面邊界檢測中，提升了圖像的信噪比和檢測精度。蘇少輝等［8］利用雙目相機(jī)結(jié)合改進(jìn)Hough 變換檢測大尺度下矩形工件，重復(fù)定位誤差小于0.05mm。

針對三角形工件測量精度不高的問題，本文首先利用改進(jìn)Canny算子檢測圖像邊緣并結(jié)合改進(jìn)Hough 變換提取圖像直線像素，然后設(shè)置直線像素點到所屬直線的距離閾值排除非直線點信息，最終通過最小二乘法擬合剩余直線像素點提高檢測精度。

1 尺寸檢測算法流程

算法流程如圖1所示。

圖1 算法流程

1.1 圖像預(yù)處理

由于圖像采集受到環(huán)境光照和相機(jī)內(nèi)部硬件的影響，導(dǎo)致在數(shù)字化過程中和數(shù)據(jù)傳輸中會產(chǎn)生噪聲，而降低了圖像分辨率并對圖像邊緣像素檢測造成誤差，影響精密工件的檢測精度。因此，首先對采集到的圖像進(jìn)行預(yù)處理。目前常用的圖像預(yù)處理方法有高斯濾波、中值濾波、維納濾波、均值濾波、雙邊濾波［9］等。由于高斯濾波是對同一距離的像素進(jìn)行同一尺度濾波計算而沒有考慮空間像素的差異，導(dǎo)致邊緣出現(xiàn)模糊的情況，從而降低了工件測量尺寸。中值濾波則采用鄰域像素取中值的方式代替中心像素點去除圖像中的噪聲，該方法僅對椒鹽噪聲有較好的處理效果。雙邊濾波則較高斯濾波考慮了像素點的空間信息與值域信息，該方法對將要進(jìn)行濾波的像素區(qū)域先做分割與分類處理，再對不同所屬區(qū)域類別的像素計算不同權(quán)重，最后進(jìn)行鄰域加權(quán)求和處理，從而達(dá)到很好地保留邊緣并去除噪聲的效果。雙邊濾波結(jié)果H如公式（1）所示：

其中Gr，Gs，Wq如公式（2）～（4）所示：

其中，Im為輸入圖像，In為輸出圖像，m代表目標(biāo)像素，n代表任意像素，σr表示值域標(biāo)準(zhǔn)差，σs表示空間域標(biāo)準(zhǔn)差，Gr，Gs表示值域信息和空間域信息，Wq表示窗口內(nèi)值域權(quán)重和。原圖和濾波效果分別如圖2 和如圖3所示。

圖2 原圖

圖3 濾波后圖像

1.2 邊緣檢測

邊緣檢測本質(zhì)上是加強(qiáng)輪廓像素，而邊緣是圖像灰度值變化最大的位置，因此通過計算梯度可以實現(xiàn)邊緣檢測。Canny 邊緣檢測［10-11］方法則能較好地判定邊緣位置。Canny邊緣檢測分為四步：

第一步：對輸入圖像進(jìn)行高斯去噪。由于高斯噪聲對邊緣產(chǎn)生模糊效果，因此在進(jìn)行邊緣檢測之前運用雙邊濾波改進(jìn)Canny邊緣檢測。

第二步：計算圖像梯度和幅值。運用改進(jìn)八方向Sobel 算子［12］進(jìn)行邊緣檢測。G0、G45、G90、G135、G180、G225、G270、G315分別表示八個方向角，因此計算梯度及梯度幅度值如公式（5）～（7）所示：

其中i，j表示圖像像素坐標(biāo)。

第三步：非極大值抑制。在得到圖像中每一點的梯度幅值和方向后，根據(jù)所計算信息去除所有非邊緣點，細(xì)化邊緣。

第四步：雙閾值確定邊緣位置。通過設(shè)置一個高閾值與低閾值，從而防止將偽邊緣信息判定為邊緣。

由于Canny檢測算法閾值由人工設(shè)定，本文采用OTSU 自適應(yīng)閾值分割改進(jìn)Canny 需要人工標(biāo)注閾值的問題。該算法事先假定閾值threshold能將像素分為兩類，并計算劃分為兩類像素的均值M1和M2，以及圖像全局均值MT，并計算被分為兩類像素的概率P1和P2，從而計算類間方差：

式（8）出現(xiàn)最大化的灰度級就是所計算的自適應(yīng)閾值。最終圖像邊緣檢測結(jié)果如圖4 所示。從圖4（a）中可見工件表面出現(xiàn)少許非邊緣信息，這將影響直線段信息的判斷情況，經(jīng)過自適應(yīng)閾值后，最終改進(jìn)結(jié)果如圖4（b）所示，較好地排除了非邊緣點信息，提高了后續(xù)直線檢測的準(zhǔn)確性。

圖4 Canny檢測對比圖

1.3 直線檢測

Hough變換直線檢測為現(xiàn)如今最有效的直線檢測方法，該方法在1962 年由Paul Hough 首次提出。Hough變換運用直角坐標(biāo)到極坐標(biāo)的映射關(guān)系找到在極坐標(biāo)中出現(xiàn)峰值的點，從而判斷出直線信息。

然而當(dāng)直線相對于圖像水平軸存在較大角度時，像素點會出現(xiàn)巨大的錯位現(xiàn)象，這將導(dǎo)致同一條直線會檢測出多條直線或是檢測不出直線的現(xiàn)象，影響直線長度尺寸標(biāo)定。因此，本文根據(jù)邊緣點位置在像素點八鄰域周圍進(jìn)行像素填充，進(jìn)而擴(kuò)大同一條直線像素點的檢測范圍。最終改進(jìn)圖像邊緣如圖6所示，與圖5所示的原直線邊緣相比，其邊緣像素得到擴(kuò)充，直線像素更加平緩。在圖6 邊緣圖像基礎(chǔ)上通過Hough 變換進(jìn)行直線檢測，未改進(jìn)邊緣圖像直線檢測如圖7（a）所示，該圖中存在一條直線多段檢測與直線漏檢問題，圖7（b）為改進(jìn)后邊緣圖像直線檢測結(jié)果，通過實驗證明該方法很好地解決了Hough直線檢測中漏檢或多檢的問題。

圖5 原直線邊緣圖像

圖6 改進(jìn)后直線邊緣圖像

圖7 直線檢測圖像

1.4 尺寸檢測

由于直線段邊緣出現(xiàn)毛刺或缺陷等問題，導(dǎo)致最終檢測出的直線像素點無法滿足真實直線段數(shù)據(jù)的要求，因此在得到直線像素點后通過最小二乘法擬合所有屬于同一條直線像素點信息，提高檢測精度。

由于在邊緣檢測中增加了八鄰域連通區(qū)間像素擴(kuò)充，這將導(dǎo)致存在一些非直線點像素位置被計入該條直線，因此需要計算直線所屬像素點到其直線距離，并設(shè)置閾值，當(dāng)距離數(shù)值大于閾值則判定為非直線點像素，且不計入最小二乘法擬合，從而提高擬合精度。

直線點擬合方程由公式（9）所示：

其中，（xi,yi）表示直線像素點位置，a，b表示待擬合直線參數(shù)。最終通過計算待求參數(shù)偏導(dǎo)數(shù)得出最佳參數(shù)估計值：

針對三角形直線區(qū)域線段長度標(biāo)定，需計算每條直線與其他直線的交點，從而判斷三角形區(qū)域三個頂點位置，進(jìn)而計算各個點的歐式距離，如公式（12）所示，其中i= 1,2,3;j= 1,2,3分別表示不同直線段編號。最終計算結(jié)果Di表示每段直線標(biāo)定長度。最終檢測結(jié)果如圖8所示。圖8中A、B、C分別表示三角形檢測交點。

圖8 直線交點檢測圖像

2 結(jié)果與結(jié)論

測試相機(jī)選用?？低?00萬像素工業(yè)面陣工業(yè)相機(jī)型號：MV-CS060-10GM/C，傳感器型號：Sony IMX178，像元尺寸：2.4μm*2.4μm，分辨率：3072*2048。

拍攝鏡頭選用DTCM110-36，0.2 倍雙側(cè)遠(yuǎn)心鏡頭。

實驗測試設(shè)備選用處理器：Intel（R）Core（TM）i5-10210U CPU @ 1.60GHz，內(nèi)存：16.00 GB（2133 MHz），顯卡：NVIDIA GeForce MX250。

本文將圖8 中AC、AB、BC段作為尺寸標(biāo)定直線，BC段加工尺寸為5 mm，AB與AC尺寸為5.6 mm。在測量中通過對該圖像進(jìn)行旋轉(zhuǎn)、重復(fù)拍攝等方法進(jìn)行了10 次測量，得到最終標(biāo)定尺寸誤差。直線長度平均誤差達(dá)到0.0021 mm，-0.0013 mm和0.0027mm，最大誤差小于0.007mm，10 次檢測時間共計21.45 s，平均檢測時間為2.145 s。有部分測量誤差數(shù)據(jù)較大或較小，該類數(shù)據(jù)可能是拍攝過程中光照和拍攝原因?qū)D像質(zhì)量產(chǎn)生了影響。

表1 直線和圓形工件測量結(jié)果單位：mm

3 結(jié)語

為更精確地檢測三角工件長度尺寸，提出基于改進(jìn)Canny 邊緣檢測和改進(jìn)Hough 直線檢測算法，依據(jù)仿真實驗表明：

（1）結(jié)合OTSU 的自適應(yīng)閾值分割方法和八方向Sobel 算子計算梯度改進(jìn)Canny 邊緣檢測，提升了邊緣檢測精度；

（2）運用Hough 變換檢測三角直線區(qū)域，并采用八連通區(qū)域增加直線邊緣像素，提高了直線判斷準(zhǔn)確性，最終通過歐氏距離判斷屬于每一條直線的邊緣像素點，利用最小二乘法擬合該類邊緣點提高了工件尺寸檢測精度。