基于PCB圖像預(yù)處理的研究及其實(shí)現(xiàn)

繆 皓， 趙紅東， 蘇 周， 李 杰

（河北工業(yè)大學(xué) 信息工程學(xué)院，天津 300401）

0 引言

在信息產(chǎn)業(yè)中 PCB是一個(gè)不可缺少的重要支柱，PCB作為各種電子產(chǎn)品的基本零組件和集成各電子元器件的信息載體向著高性能化、高速化、輕薄短小化方向得到了快速發(fā)展，其技術(shù)和復(fù)雜程度已經(jīng)達(dá)到了一個(gè)非常高的水平。因此，隨著PCB領(lǐng)域的不斷擴(kuò)大，PCB的重要性也在進(jìn)一步提高[1]。系統(tǒng)在 PCB 待測灰度圖像的采樣、量化、以及傳輸過程中，由于光照的不均變化、電荷耦合元件（CCD，Charge-Coupled Device）攝像機(jī)自身的傳感器噪聲、模數(shù)采集（AD，Analog-to-Digital）過程的量化噪聲、感光過程產(chǎn)生的顆粒噪聲、以及人為因素的輕微抖動等，在傳輸、接收處理得到的圖像都不可避免受到內(nèi)部器件和外部環(huán)境的影響，從而使圖像質(zhì)量失真，信噪比下降。為了降低噪聲，可以使用平滑濾波器對待測圖像進(jìn)行濾波處理，但選擇不同大小的平滑濾波器對處理后圖像模糊程度不同，因此要提高圖像的質(zhì)量要求所使用的濾波器不僅能夠有效地去除噪聲并且能夠盡可能保留圖像的原貌[2]。

1 圖像增強(qiáng)

圖像增強(qiáng)是改進(jìn)圖像質(zhì)量的技術(shù)，相對于圖像識別的一種前期處理，根據(jù)圖像增強(qiáng)處理的空間不同，可分為空間域處理以及頻率域處理兩大類，前者包括對圖像的灰度級作用以及直方圖修正，都是直接對像素灰度值進(jìn)行處理；后者是對圖像的頻譜成分進(jìn)行分析操作，經(jīng)過傅里葉變換，對圖像頻譜的高頻低頻部分進(jìn)行處理，再經(jīng)逆傅里葉變換獲得所需圖像結(jié)果。

采集的 PCB 圖像由于外部曝光以及在信道傳輸過程中的其他干擾因素，使得圖像明暗度降低、噪聲加大，為有效的消除噪聲干擾并且使圖像的明暗對比增強(qiáng)，本文在PCB 待測圖像增強(qiáng)中選擇空間域的灰度變換處理和圖像平滑處理。

1.1 空間域的灰度變換

灰度變換作為圖像增強(qiáng)的一種重要手段，可以使圖像動態(tài)范圍加大、使圖像對比度擴(kuò)展，使圖像特征更加明顯來改善圖像顯示效果?；叶茸儞Q可分為線性變換和非線性變換。令原圖像 m(x,y)的灰度范圍為[a,b],線性變換后圖像 n(x,y)的灰度擴(kuò)展為[c,d],它們之間的關(guān)系為：

由于曝光過度或不足，CCD采集的PCB圖像灰度可能在一個(gè)很小的范圍內(nèi)變化，在電腦上可能看到的是沒有灰度層次以及模糊的圖像。采用線性變換可以對模糊圖像的各個(gè)像素灰度作線性拉伸，能夠使圖像視覺效果有效地得到改善。

為提高對圖像識別的后期處理以及特征提取，先對原始 PCB圖像進(jìn)行二值化的灰度處理，然后采用直方圖修正技術(shù)的圖像域方法使圖像均衡化[3]。

1.2 空間域圖像平滑

圖像平滑的目的是為了減少和消除圖像噪聲以改善圖像質(zhì)量，以便圖像分割、圖像識別等后續(xù)處理[4]。在空間域內(nèi)可以用鄰域平均來減少噪聲；在頻率域，由于發(fā)生高頻段的噪聲頻譜概率更大，可采用各種形式的低通濾波。在空間域圖像平滑主要有噪聲門限、鄰域平均、加權(quán)平均、中值濾波等方法。

趙曉霞提出一種偏微分的圖像增強(qiáng)方法，這種方法降噪采用全變分（TV，Total Variation）模型與直方圖修正技術(shù)來實(shí)現(xiàn)，這種方法能有效地保留邊緣也能夠使對比度增強(qiáng)[5]。

萬睿等人先對待測圖像進(jìn)行二值化處理，然后通過閾值分割處理，得到二值化閾值，在去除噪聲的同時(shí)使處理速度也得到了提高[6]。

1.2.1 噪聲門限

噪聲門限方法是一種常見的噪聲消除辦法，它對噪聲的執(zhí)行效率高，并且消除簡單易行。它對圖像平滑時(shí)，第一是門限閾值的設(shè)定，門限閾值的設(shè)定直接影響濾波效果以及圖像細(xì)節(jié)，然后根據(jù)圖像的特性，對每個(gè)像素進(jìn)行順序的檢測，根據(jù)公式與其鄰域的所有像素值進(jìn)行對比，判斷是否該像素為噪聲。若不是噪聲，輸出像素原值，若為噪聲，則輸出為鄰域內(nèi)其它像素灰度的平均值。該方法中門限閾值 T 的選擇至關(guān)重要，T選擇太大或太小，都會或多或少導(dǎo)致噪聲平滑不夠或者圖像變得模糊[7]。

1.2.2 中值濾波

傳統(tǒng)中值濾波算法運(yùn)算量主要集中在窗口數(shù)據(jù)的排序上，為了減少排列次數(shù)更快的得到中值，提出了一種快速中值濾波算法，將陣列分解為一維陣列進(jìn)行運(yùn)算，先取陣列水平各行的中間值，再取水平各行中間值的中值，作為濾波最后的輸出。

該方法可表示為：

其中Mkn表示濾波窗口的第k行第n列的元素，該快速算法是將二維陣列簡化為一維陣列進(jìn)行比較中值運(yùn)算，此二維陣列中含有n×n個(gè)像素，簡化的一維陣列為n+1個(gè)僅含n個(gè)像素。這樣，相比于傳統(tǒng)中值濾波的 n2(n2-1)/2次運(yùn)算，快速算法只需進(jìn)行n( n2-1)/2次運(yùn)算，速度得到很大提高，復(fù)雜度減少為原來的1/n。但是和普通的中值濾波一樣，濾波窗口太小,圖像細(xì)節(jié)得到更好的保護(hù), 但會限制濾除噪聲的能力；濾波窗口太大,噪聲可以得到很好的濾除, 但圖像細(xì)節(jié)輪廓會有很大的損失，使圖像模糊不清[8]。

1.2.3 快速加權(quán)中值濾波

為了解決既要降噪, 又要保護(hù)圖像細(xì)節(jié)這一矛盾本文提出了一種快速加權(quán)中值濾波算法,不僅濾波速度上提高了，噪聲得到很好的濾除，還可以保護(hù)更多的圖像細(xì)節(jié)[9]。通過加權(quán),使在窗口中的中心象素比重增加,對于輸入:(I1, I2,… ,IN)，輸出加權(quán)中值濾波:

r代表加權(quán)，MED{}表示取中值運(yùn)算函數(shù),W=(W1, W2,… ,WN)表示加權(quán)窗口的權(quán), 并且規(guī)定 W= ( 1, 1,…, 1) 時(shí)為標(biāo)準(zhǔn)窗口。規(guī)定：

快速加權(quán)中值濾波中值運(yùn)算 MED的運(yùn)算過程是: 先對加權(quán)后的輸出中的Nc個(gè)數(shù)按照順序升冪排序,排序后輸出的Nc個(gè)數(shù)的第T個(gè)數(shù)就是最后的中值輸出。

2 MATLAB圖像增強(qiáng)處理

MATLAB支持包括索引圖像、灰度圖像、二值圖像、RGB圖像和多幀圖像陣列的5種圖像類型[10]；

采集的PCB圖像可以對其先增強(qiáng)對比度，再進(jìn)行去噪，若運(yùn)用該方法，則圖像原有的噪聲在增強(qiáng)對比的過程中也被大大增加了，使后續(xù)圖像去噪處理得不到很好的效果。

因此本文采用對圖像先去噪后增強(qiáng)，經(jīng)采集得到的 PCB 圖像上，可能存在噪聲需要去噪，光源強(qiáng)度不夠，圖像可能整體偏暗，首先用RGB2GRAY函數(shù)將采集到的圖片轉(zhuǎn)換為灰度圖像（256色），對含有椒鹽噪聲的PCB圖像分別用本文提出的方法以及傳統(tǒng)中值濾波的方法進(jìn)行比較，去噪后再利用灰度變換對PCB圖像進(jìn)行增強(qiáng)處理,圖1就是用這種方法處理后的結(jié)果。

圖1 各種濾波算法實(shí)現(xiàn)結(jié)果比較

從圖 1（a）和圖 1（c）可以看出對于有椒鹽噪聲的PCB圖像，傳統(tǒng)的3×3中值濾波處理孤立點(diǎn)的噪聲雖然具有很好的效果，但是PCB圖像上的走線，器件明顯模糊，圖像的細(xì)節(jié)輪廓，質(zhì)量不能達(dá)到很好的效果。圖1（e）是對帶有椒鹽噪聲的PCB圖像進(jìn)行快速加權(quán)中值濾波的仿真圖，去噪具有明顯效果，而且濾波速度也得到了提高，圖像的細(xì)節(jié)，走線都能夠很好的保持。

圖1（b）、圖1（d）、圖1（f）是PCB圖像經(jīng)去噪處理前后的灰度直方圖，圖中橫坐標(biāo)為灰度等級, 縱坐標(biāo)為各個(gè)灰度等級對應(yīng)的像素?cái)?shù)。其中像素值 255是椒鹽噪聲，比較兩幅圖可以明顯看出,兩種濾波方法都能將噪聲很好的去除，但是中值濾波把許多灰度等級都濾除掉了，例如150-200的灰度等級，使邊界變得模糊，而且波峰波谷都遭到了破壞，加權(quán)中值濾波去能克服這些缺點(diǎn)，灰度等級，灰度動態(tài)范圍都能夠很好的保留，并且對于波峰、波谷的位置都給予了很好的保留，因此圖像許多細(xì)節(jié)都看的比較清晰，為了更好地便于分析，圖2給出了圖像經(jīng)濾波變換前后的頻譜。

圖2 濾波處理前后的頻譜圖

3幅圖是經(jīng)傅里葉變換得到的頻譜，為了便于分析，根據(jù)傅里葉變換的周期性，把零頻移到中間。幅值譜的明暗程度反映了各頻率成分的幅值大小，圖像的能量主要集中在低頻段（中央部分）。

盡管高頻段包含少量的能量，但包含了圖像的重要信息，圖像的邊緣信息屬于高頻信息，同樣，噪聲的灰度級變化迅速，也是高頻信息。從圖中可以看出，兩種濾波方法都在一定程度上抑制圖像的高頻成分，能夠有效濾除高頻段的噪聲,但是低頻分量,也就是灰度變化平緩的部分-圖像輪廓信息遭到了破壞，從直方圖也可以看出變化緩慢的150～200灰度級遭到了破壞，而加權(quán)中值能夠很好的保護(hù)輪廓信息，為了使細(xì)節(jié)更突出圖3是增強(qiáng)后圖像。

圖3 增強(qiáng)后的效果

從灰度處理后的直方圖可以看出，直方圖占據(jù)了整個(gè)圖像灰度值允許的范圍，增加了圖像灰度動態(tài)范圍，也增加了圖像的對比度，在圖像視覺上有了較大的反差，使細(xì)節(jié)更加突出。

3 結(jié)語

本文主要對采集的PCB圖像整體偏暗，對比度差，噪聲大進(jìn)行圖像的預(yù)處理，主要包括空間域的濾波技術(shù)以及圖像增強(qiáng)的灰度變換，由于傳統(tǒng)的中值濾波受濾波窗口尺寸的影響較大，處理后圖像的細(xì)節(jié)變得模糊不清，通過一種改進(jìn)的加權(quán)中值濾波算法，根據(jù)灰度直方圖以及頻譜圖分析了處理后圖像質(zhì)量的好壞，結(jié)果表明其濾波速度及質(zhì)量大大超過傳統(tǒng)的中值濾波，顯著地提高了PCB圖像走線，元器件等邊緣，整體圖像輪廓的保護(hù)能力，最后通過灰度變換使圖像的細(xì)節(jié)變得更加清晰，提高PCB圖像的識別率的后期處理。

[1]胡秀亮,高志,葉健敏,等.基于圖像處理技術(shù)的 PCB 缺陷自動檢測系統(tǒng)[J].圖形圖像與多媒體,2010,9(01):24-27.

[2]喬鬧生,趙華君.保留灰度級的 PCB 圖像增強(qiáng)研究[J].計(jì)算機(jī)工程與應(yīng)用,2009,45(09):181-183.

[3]包建華.印刷電路板圖件識別研究中的圖像預(yù)處理實(shí)現(xiàn)[J].電子工程師,2006,8(08):35-37.

[4]李玉峰,葛雯.小波域含噪圖像的壓縮方法研究[J].通信技術(shù),2010,43(01):129-130.

[5]CASELLES V,LISANI J L,MOREL J I,et al.Shapepreserving Local Contrast Enhancemen[C].USA: [s.n.],1997:220-230.

[6]趙曉霞,王明泉,李高亮.一種基于偏微分方程的 PCB圖像增強(qiáng)方法[J].電視技術(shù),2012,36(03):33-35.

[7]林越偉.PCB 質(zhì)量檢測的算法研究[D].重慶:重慶大學(xué)光電工程學(xué)院,2007.

[8]董漢磊,徐力平,高穎敏,等.基于 Bior小波變換和中值濾波的圖像去噪[J].通信技術(shù),2009,42(11):171-172.

[9]鄧秀勤,熊勇.用于圖像處理的加權(quán)中值濾波算法[J].計(jì)算機(jī)技術(shù)與發(fā)展,2009,3(03):46-48.

[10]楊帆.數(shù)字圖像處理與分析[M].北京:北京航空航天大學(xué)出版社,2007:9-15.

[11]趙凱,賈玉韜,王丹.基于 IXP425的防火墻硬件平臺設(shè)計(jì)[J].信息安全與通信保密,2007(01):111-113.

[12]李侃.粗糙集文本過濾系統(tǒng)研究及硬件實(shí)現(xiàn)[J].信息安全與通信保密,2006(03):44-45,48.