細(xì)化算法在干涉圓環(huán)參數(shù)提取中的作用

摘" 要： 為解決光學(xué)干涉圖像亮暗不均和噪聲對(duì)干涉圓環(huán)半徑和圓心識(shí)別的影響問題，以Matlab仿真圖像為對(duì)象，研究細(xì)化算法在圓環(huán)參數(shù)提取中的作用。首先，在仿真干涉同心圓環(huán)圖像中加入噪聲和亮暗不均的因素；然后，在預(yù)處理基礎(chǔ)上依次應(yīng)用二值算法、細(xì)化算法得到二值圖和細(xì)化圖；最后，分別在預(yù)處理圖、二值圖和細(xì)化圖上進(jìn)行Hough變換提取半徑和圓心。對(duì)比研究發(fā)現(xiàn)，二值化基礎(chǔ)上的細(xì)化處理可以極大地減少亮暗不均的影響，抑制噪聲，不僅能夠消除粗大誤差，提高提取精度，而且使得識(shí)別結(jié)果變得更加穩(wěn)定。

關(guān)鍵詞： 圖像處理； 細(xì)化算法； 干涉圓環(huán)； Hough變換； 圓心識(shí)別； 圓環(huán)參數(shù)提取

中圖分類號(hào)： TN911.7?34； TP394.1； TH691.9" " " " "文獻(xiàn)標(biāo)識(shí)碼： A" " " " " " " 文章編號(hào)： 1004?373X（2024）17?0149?04

Role of thinning algorithm in interference ring parameter extraction

OUYANG Yefeng， TANG Jiali， ZHU Yuyang， FU Xuezhou， ZHANG Baowu

（College of Metrology and Measurement Engineering， China Jiliang University， Hangzhou 310018， China）

Abstract： In view of the uneven brightness of the optical interference images and the influence of noise on the identification of the radius and center of the interference rings， the role of thinning algorithm in the extraction of the ring parameters is studied by taking the Matlab simulation image as the object. Noise and uneven brightness are added to the simulated interference concentric ring image. On the basis of preprocessing， binary algorithm and thinning algorithm are applied in turn to obtain binary image and thinned image. Hough transform is performed on the preprocessed image， binary image and thinned image to extract the radius and center of the circle， respectively. The comparative study shows that the thinning processing based on binarization can reduce the influence of uneven brightness and suppress noise greatly， which can not only eliminate the gross error and improve the extraction accuracy， but also make the recognition results more steady.

Keywords： image processing; thinning algorithm; interference ring; Hough transform; centre identification; ring parameter extraction

0" 引" 言

干涉測量作為精密測量中常見的無損測量方式，在測角[1?2]、測距[3?4]、流體計(jì)量[5?6]等各領(lǐng)域被廣泛應(yīng)用。干涉測量法根據(jù)光源的選取、參考光路的設(shè)置以及干涉圖像的形成等不同的因素，分為各種不同的類別。但就其本質(zhì)而言，都是將被測量按一定物理關(guān)系轉(zhuǎn)化為干涉圖像信息，再通過提取這些信息反演出被測量。盡管在不同的應(yīng)用場景中測量形式各不相同，但在光學(xué)測量手段中，對(duì)干涉圖像的信息提取都是測量中極為重要的一環(huán)。

點(diǎn)光源非定域干涉為角度和厚度精密測量提供了一種新的可能。這種技術(shù)在實(shí)際應(yīng)用中的關(guān)鍵環(huán)節(jié)就是干涉同心圓環(huán)半徑和圓心的精確提取。得益于計(jì)算機(jī)技術(shù)的發(fā)展，圖像處理手段日益精進(jìn)，其處理速度和精度都有飛躍式的發(fā)展[7?10]，這為干涉圓環(huán)的參數(shù)提取算法實(shí)現(xiàn)提供了極大的便利。但是，在實(shí)際的干涉圖像采集中，亮暗分布不均是難以避免的問題，這一現(xiàn)象會(huì)影響圓環(huán)參數(shù)的提取精度，甚至可能會(huì)導(dǎo)致粗大誤差。針對(duì)這一問題，常用的解決方式為使用自適應(yīng)閾值二值算法對(duì)圖像進(jìn)行二值處理，然后，在二值圖上進(jìn)行Hough變換，實(shí)現(xiàn)圓心和半徑的識(shí)別。經(jīng)研究，在Hough變換之前先對(duì)二值圖進(jìn)行細(xì)化可以更好地減少亮暗分布不均造成的影響，大大提高圓環(huán)參數(shù)識(shí)別的精度。比較經(jīng)典的細(xì)化算法有Hilditch算法[11]、Zhang并行化算法[12]、Deutsch算法[13]、Yangtagai極值檢測算法[14]等。本文以點(diǎn)光源非定域干涉仿真圖像為對(duì)象，在加入噪聲和亮暗不均的因素后，對(duì)比研究了單獨(dú)預(yù)處理、預(yù)處理+二值化、預(yù)處理+二值化+細(xì)化處理后的Hough變換半徑和圓心提取結(jié)果，以此詳細(xì)考察了Zhang并行快速細(xì)化算法在干涉圓環(huán)參數(shù)提取中的作用。

1" 點(diǎn)光源非定域干涉同心圓環(huán)

點(diǎn)光源的非定域相干疊加的原理如圖1a）所示，兩個(gè)相距為[d]的點(diǎn)光源[S1]和[S2]，它們各自發(fā)出的光波分別經(jīng)過[r1]和[r2]路徑后在空間[P]點(diǎn)相遇疊加。在[P0]點(diǎn)垂直于[x]軸的平面內(nèi)就可以觀察到一個(gè)同心圓環(huán)干涉圖形，如圖1b）所示。其中，干涉圓環(huán)的圓心即為圖1a）的[P0]點(diǎn)，任意一環(huán)的半徑記為[Ri]。為了對(duì)比，圖1c）給出了仿真圖片。另外，圖1d）和圖1e）分別給出了圖1b）和圖1c）中過圓心水平方向上的強(qiáng)度輪廓線。

對(duì)比圖1c）的仿真干涉圖像，圖1b）的實(shí)驗(yàn)圖片左右兩側(cè)的亮暗分布有著明顯的差異。這種差異會(huì)影響條紋分布的均勻性，如圖1d）所示，即暗區(qū)域內(nèi)的峰值明顯小于亮區(qū)域內(nèi)的峰值。另外，實(shí)驗(yàn)干涉圖像噪聲信號(hào)使圓環(huán)中心區(qū)域以及波峰（對(duì)應(yīng)亮環(huán)）和波谷（對(duì)應(yīng)暗環(huán)）分裂嚴(yán)重，或者存在飽和現(xiàn)象，極值點(diǎn)難以確定。這一干擾使干涉圓環(huán)參數(shù)（半徑和圓心）的識(shí)別精度受到影響。

2" 細(xì)化算法原理

細(xì)化算法是通過不斷連續(xù)剝離二值圖像的邊界元素，使得原圖像變?yōu)槿舾蓷l像素寬度為1的連通線的圖像處理方法。對(duì)比研究發(fā)現(xiàn)，在干涉圓環(huán)圖像細(xì)化處理方面，Zhang并行快速細(xì)化算法的綜合表現(xiàn)最好[15]。Zhang細(xì)化算法采用邏輯算術(shù)運(yùn)算進(jìn)行刪除判斷，使得迭代次數(shù)減小以提升執(zhí)行速度，其代價(jià)為細(xì)化后的圖像會(huì)殘留一些交叉和分支，不過，因?yàn)镠ough圓檢測采用梯度投票制，Zhang細(xì)化算法的殘留現(xiàn)象對(duì)結(jié)果的影響很小。

Zhang細(xì)化算法主要過程分為以下兩個(gè)步驟。

步驟1：建立如圖2所示的3×3滑動(dòng)窗口以遍歷全圖，每次滑動(dòng)讀入9個(gè)像素點(diǎn)，其中亮點(diǎn)值設(shè)為1，暗點(diǎn)為0。

步驟2：每滑動(dòng)一次窗口，判斷讀入的像素是否滿足以下兩組條件：

① [2≤N（P1）≤6]

② [S（P1）=1]

③ [P2×P4×P6=0]

④ [P4×P6×P8=0]

或：

① [2≤N（P1）≤6]

② [S（P1）=1]

③ [P2×P4×P8=0]

④ [P2×P6×P8=0]

其中：[N（P1）]表示以[P1]點(diǎn)為中心的窗口內(nèi)值為1的個(gè)數(shù)；[S（P1）]表示[P1]點(diǎn)的值。只要滿足兩組條件中的任意一組，則將[P1]點(diǎn)的像素進(jìn)行置0處理。

3" 結(jié)果分析

為了考察細(xì)化算法的作用效果，對(duì)圖1c）的仿真干涉同心圓環(huán)圖像進(jìn)行噪聲和亮暗不均效果修改（見圖3a）），依次進(jìn)行預(yù)處理（見圖3b））、二值化處理（見圖3c））和細(xì)化處理（見圖3d））。最后，遵循圖4的流程圖，對(duì)這些處理后的干涉圖像進(jìn)行Hough變換，完成圓環(huán)參數(shù)的識(shí)別，如圖5所示。其中，圖5a）是單獨(dú)預(yù)處理（方式一）后的識(shí)別結(jié)果，圖5b）是預(yù)處理+二值化（方式二）的識(shí)別結(jié)果；圖5c）是預(yù)處理+二值化+細(xì)化算法（方式三）后的識(shí)別結(jié)果。圓心處的點(diǎn)代表所識(shí)別圓環(huán)的圓心。

如圖5所示，按方式一和方式二處理后識(shí)別的圓心存在離散現(xiàn)象，識(shí)別到的圓環(huán)與實(shí)際環(huán)也存在明顯的偏差，而方式三多次識(shí)別的圓心基本重合，并且識(shí)別的半徑與條紋中心基本吻合。為進(jìn)一步分析，對(duì)不同處理方法所得的數(shù)據(jù)進(jìn)行列表對(duì)照分析，見表1。

表1中，真實(shí)值由仿真數(shù)據(jù)和數(shù)學(xué)公式直接計(jì)算得到。對(duì)比表中數(shù)據(jù)可以看出，亮暗不均和噪聲的影響對(duì)干涉圓環(huán)的識(shí)別有很大的影響，甚至可能會(huì)出現(xiàn)粗大誤差，如方式一和方式二第1環(huán)的圓心識(shí)別偏差較大，而方式一第2、3環(huán)與方式二第2、3環(huán)的半徑識(shí)別偏差較大。而方式三可以有效提升圓心半徑提取結(jié)果的精度和穩(wěn)定性，其識(shí)別的圓心和半徑的誤差均在2 pix以下。

4" 結(jié)" 語

本文基于實(shí)驗(yàn)現(xiàn)象，在仿真干涉同心圓環(huán)中主動(dòng)添加亮暗不均和噪聲，在三種圖像處理基礎(chǔ)上分別利用Hough變換進(jìn)行圓環(huán)半徑和圓心的識(shí)別，對(duì)比研究了不同處理方式下的識(shí)別結(jié)果，以此驗(yàn)證細(xì)化算法在減少亮暗不均的影響、抑制噪聲方面的突出作用。在二值化處理后添加細(xì)化處理過程，識(shí)別的圓心和半徑的誤差均在2 pix以下，不僅識(shí)別精度更高，而且結(jié)果更穩(wěn)定。研究結(jié)果說明，細(xì)化算法在干涉圓環(huán)參數(shù)提取過程中有非常重要的作用。

注：本文通訊作者為張寶武。

參考文獻(xiàn)

[1] 王漢斌，黃垚，薛梓，等.激光小角度基準(zhǔn)測量不確定度評(píng)定[J].計(jì)量學(xué)報(bào)，2019，40（6）：970?974.

[2] 張寶武，崔建軍，歐陽燁鋒，等.雙臂非共程點(diǎn)光源非定域干涉小角度測量方法[J].計(jì)量學(xué)報(bào)，2023，44（8）：1202?1207.

[3] 包為政，張福民，曲興華.基于等光頻細(xì)分重采樣的調(diào)頻干涉測距方法[J].激光技術(shù)，2020，44（1）：1?6.

[4] DZYUBENKO M I， KOLENOV I V， PELIPENKO V P， et al. Laser rangefinder of relatively eye?safe wavelength range [J]. Radiofizika i elektronika， 2021， 26（1）： 41?48.

[5] 周國全，段玉婕，朱增旺.基于等傾干涉與圖像處理算法的流體折射率測量傳感系統(tǒng)[J].物理與工程，2021，31（6）：17?24.

[6] 李虎，郭子龍，楊文婷，等.空芯光纖多模干涉型光纖液位傳感技術(shù)研究[J].激光技術(shù)，2022，46（1）：120?124.

[7] 梅啟升，王敏，周群.多噪聲干涉條紋圖像的檢測方法[J].激光與光電子學(xué)進(jìn)展，2019，56（12）：103?108.

[8] 高銘鍔，袁鶴，徐敏，等.激光干涉條紋特征提取及定中測量技術(shù)研究[J].光子學(xué)報(bào)，2022，51（2）：311?318.

[9] 劉逸飛，蘇亞，姚曉天，等.OCT無創(chuàng)血糖檢測圖像處理最優(yōu)化方法研究[J].激光技術(shù)，2023，47（2）：178?184.

[10] 蒲亮，葉玉堂，吳云峰，等.基于經(jīng)驗(yàn)?zāi)B(tài)分解的激光干涉信號(hào)處理[J].激光技術(shù)，2011，35（3）：299?301.

[11] 陳家鳳，陳海清，任海霞.基于小波變換的最小二乘相位解纏算法[J].光學(xué)技術(shù)，2007（4）：613?616.

[12] 崔鳳奎，王曉強(qiáng)，張豐收，等.二值圖像細(xì)化算法的比較與改進(jìn)[J].洛陽工學(xué)院學(xué)報(bào)，1997（4）：50?54.

[13] DEUTSCH E S. Thinning algorithms on rectangular， hexagonal， and triangular arrays [J]. Communications of the ACM， 1972， 15（9）： 827?837.

[14] 戴福隆，王朝陽.條紋圖象的數(shù)字化自動(dòng)分析處理技術(shù)之一：條紋中心法[J].光子學(xué)報(bào)，1999（8）：700?706.

[15] 李偉.光干涉圖像處理中的關(guān)鍵算法研究[D].西安：西安科技大學(xué)，2013.