何文耀， 劉林威， 黃圳煜， 易子川， 周國(guó)富， 水玲玲*

(1. 華南師范大學(xué)華南先進(jìn)光電子研究院， 廣州 510006； 2. 電子科技大學(xué)中山學(xué)院電子信息學(xué)院，中山 5284021)

電泳電子紙(EPD)作為一種新型顯示技術(shù)[1]，具有視角廣[2]、可擦寫(xiě)[3]、類紙顯示[4]、超低功耗[5]、強(qiáng)光下可讀的特點(diǎn)[6]. 在分析EPD顯示效果時(shí)，采集的畫(huà)面容易受到光干擾，影響分析結(jié)果. 因此，有效克服光干擾對(duì)圖像采集過(guò)程至關(guān)重要，也是計(jì)算機(jī)視覺(jué)有待解決的難題[7].

針對(duì)采集圖像處理過(guò)程中的光干擾問(wèn)題，不同研究者提出了大量算法，例如：楊穎等[8]提出1種基于高斯頻域低通濾波的方法，解決字符圖像光照不均的問(wèn)題；WANG等[9]提出了一種提高EPD圖像穩(wěn)定性的算法;劉先明等[10]的研究有效改善了直流平衡問(wèn)題；HSU等[11]研究了橫向電場(chǎng)對(duì)粒子驅(qū)動(dòng)的影響，但在測(cè)試數(shù)據(jù)時(shí)存在光干擾. YI等[12]提出了減弱鬼影的算法. 但這些研究中均未避免光干擾對(duì)數(shù)據(jù)采集準(zhǔn)確性的影響. 關(guān)于EPD視頻播放及驅(qū)動(dòng)系統(tǒng)設(shè)計(jì)的大量研究[1,13-18]表明：在測(cè)量反射率及分析鬼影時(shí)未濾除光干擾，會(huì)導(dǎo)致測(cè)試數(shù)據(jù)不夠準(zhǔn)確.

為了解決這些問(wèn)題，本文提出了一種基于Gabor濾波器的EPD圖像濾波系統(tǒng). 對(duì)于單一灰階圖像，先通過(guò)Gabor濾波器降低圖像中的光干擾，再使用閉運(yùn)算進(jìn)行邊緣補(bǔ)償；對(duì)于彩色圖像，分別對(duì)紅(R)、黃(G)、藍(lán)(B)分量濾波，再重新組合成彩色圖像. 該系統(tǒng)在弱化光照干擾、提高圖像質(zhì)量方面具有良好的性能.

1 研究方法

本實(shí)驗(yàn)顯示屏采用由E-Ink公司生產(chǎn)的微膠囊EPD，以微控制器為主控單元，以攝像系統(tǒng)所采集的EPD畫(huà)面為研究對(duì)象，使用MATLAB R2016a對(duì)采集的EPD畫(huà)面進(jìn)行分析，以Gabor濾波器為主要濾波工具，濾除圖像光干擾. 當(dāng)EPD畫(huà)面為單一灰階圖像時(shí)，經(jīng)過(guò)Gabor濾波后，使用閉運(yùn)算進(jìn)行邊緣補(bǔ)償；當(dāng)EPD畫(huà)面為彩色圖像時(shí)，先提取R、G、B分量，再分別對(duì)這3個(gè)分量進(jìn)行Gabor濾波處理，從而濾除光干擾.

1.1 硬件設(shè)計(jì)原理

電泳電子紙顯示畫(huà)面通過(guò)攝像系統(tǒng)將采集的JPG或BMP數(shù)據(jù)傳輸?shù)轿⑻幚砥鬟M(jìn)行處理，數(shù)據(jù)保存為“.JPG”或者“.BMP”格式并存放于SD卡，再通過(guò)USB Slave接口與PC端連接，讀取存放于SD卡的圖片. 工作流程如圖1所示.

圖1 硬件系統(tǒng)工作流程

1.2 軟件設(shè)計(jì)原理

1.2.1 單一灰階圖像濾波原理 Gabor濾波器具有較強(qiáng)的空間位置和方向選擇性，并且能夠捕捉對(duì)應(yīng)于空間和頻率的局部結(jié)構(gòu)信息，通過(guò)調(diào)整Gabor濾波器的方向、基頻帶寬及中心頻率，能夠很好地兼顧信號(hào)在空域和頻域中的分辨能力；且對(duì)于光照、圖像的亮度變化及形態(tài)具有較強(qiáng)的魯棒性，通過(guò)提取圖像在各個(gè)尺度和方向上的紋理信息，對(duì)圖像進(jìn)行濾波處理，在一定程度上降低了光照變化和噪聲對(duì)圖像的影響，提供圖像對(duì)光照變化良好的適應(yīng)性.

將采集的帶光干擾的EPD圖像作為樣本，根據(jù)圖像RGB的加權(quán)和將圖像的RGB值轉(zhuǎn)化為灰度值，實(shí)現(xiàn)灰度圖像的轉(zhuǎn)化：

H=0.299R+0.587G+0.114B.

(1)

將圖像轉(zhuǎn)換為灰度圖像后，采用Gabor濾波器去除光噪聲，Gabor濾波函數(shù)如下：

(2)

其中Sx和Sy分別代表著X和Y方向的差異系數(shù)；Sx、Sy通過(guò)改變圖像的差異系數(shù)來(lái)控制濾波強(qiáng)度；U和V分別代表X和Y方向的中心頻率；U、V通過(guò)控制中心頻率來(lái)影響濾波寬度.

經(jīng)過(guò)Gabor濾波后，圖像邊緣的亮度驟減，產(chǎn)生邊緣失真. 系統(tǒng)設(shè)計(jì)了一種基于單一灰階圖像的邊緣處理方法：當(dāng)所拍攝的電子紙圖像為單一灰階時(shí)，經(jīng)過(guò)Gabor濾波后，先對(duì)圖像進(jìn)行膨脹處理，將與物體接觸的所有背景點(diǎn)合并到該物體中，使邊界向外部擴(kuò)張，填充圖像中的小孔與圖像邊緣處的微小凹陷部分；再對(duì)圖像進(jìn)行腐蝕處理. 腐蝕是消除邊界點(diǎn)，使邊界向內(nèi)部收縮的過(guò)程；先膨脹再腐蝕的過(guò)程稱為閉運(yùn)算. 具體定義為：

結(jié)構(gòu)元素B對(duì)圖像A的膨脹

(3)

集合A被結(jié)構(gòu)元素B腐蝕

AΘB={x:Bx?A}.

(4)

閉運(yùn)算

A·B=(A⊕B)ΘB.

(5)

1.2.2 彩色圖像濾波原理 當(dāng)采集到的EPD畫(huà)面為彩色圖像時(shí)，直接使用Gabor濾波會(huì)使圖像失真，系統(tǒng)設(shè)計(jì)了一種可以保留原圖像色彩的Gabor濾波方法，具體過(guò)程：

當(dāng)采集到的畫(huà)面為彩色圖像時(shí)，先提取彩色圖像的R、G、B分量：

R=I(:,:,1)；

(6)

G=I(:,:,2)；

(7)

B=I(:,:,3).

(8)

再對(duì)3組分量依次進(jìn)行Gabor濾波，將濾波后的3組分量R′、G′、B′重新組合成彩色圖像.

I(:,:,1)=R′；

(9)

I(:,:,2)=G′；

(10)

I(:,:,3)=B′.

(11)

2 結(jié)果與討論

2.1 系統(tǒng)硬件電路及軟件設(shè)計(jì)

本系統(tǒng)硬件：微處理器、PC端、USB Slave接口、復(fù)位電路、電源電路、聯(lián)合測(cè)試工作組(JTAG)、靜態(tài)隨機(jī)存儲(chǔ)器(SRAM)、按鍵、SD卡、LED燈、攝像頭和液晶顯示器(LCD)，系統(tǒng)硬件框架如圖2A所示. Gabor濾波器濾波流程：導(dǎo)入拍攝的圖像，判斷是否為彩色；若是則先提取R、G、B分量，設(shè)置濾波參數(shù)，再分別對(duì)這3個(gè)分量進(jìn)行Gabor濾波，最后重新合成彩色圖像；若圖片為單一灰階圖片，則先將圖像轉(zhuǎn)換為灰度圖，再設(shè)置濾波參數(shù)，通過(guò)Gabor濾波后進(jìn)行閉運(yùn)算；最后輸出圖像. 系統(tǒng)軟件流程如圖2B所示.

圖2 系統(tǒng)硬件框架與軟件流程

2.2 單一灰階圖像濾波

對(duì)于1張單一灰階圖像，在采集時(shí)受到由右下角射入的燈光照射影響，導(dǎo)致圖像亮度不均勻(圖3A). 通過(guò)對(duì)圖像進(jìn)行Gabor濾波后得到如圖3B所示的圖像.

圖3 Gabor濾波前后的圖像對(duì)比

圖像灰度的方差可以表示圖像的均勻性，通過(guò)MATLAB計(jì)算出濾波后的圖像灰度發(fā)現(xiàn)，其方差下降了39.9%，這說(shuō)明圖像的均勻性提升了39.9%，也說(shuō)明濾除了39.9%的光干擾.

利用MATLAB對(duì)濾波后的圖像作出三維模擬圖，以描述圖像像素間灰度的變化情況，其中x、y代表灰度圖像中某一像素點(diǎn)的位置，z代表當(dāng)前像素點(diǎn)的灰度. 經(jīng)過(guò)Gabor濾波后，圖像邊緣的灰度出現(xiàn)很大程度的衰減(圖4A). 通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，運(yùn)用閉運(yùn)算的形態(tài)學(xué)運(yùn)算處理后，圖像邊緣衰減現(xiàn)象明顯下降(圖4B).

為了驗(yàn)證本文方法的優(yōu)越性，在相同條件下將本文方法與其他濾波方法進(jìn)行對(duì)比分析. 通過(guò)灰度化和高斯低通濾波的方法處理后，圖像灰度仍不均勻且整體偏白，可見(jiàn)圖像表面的光干擾并沒(méi)有被濾除；圖像的右下角仍有明顯的入射光痕跡. 而通過(guò)本文方法處理的圖片灰度較均勻，且右下角入射光的痕跡已基本被濾除(圖5).

圖4 Gabor濾波和閉運(yùn)算得到的三維模擬圖

Figure 4 The 3-dimensional simulation curves obtained with Gabor filtering and closed operation

圖5 本文方法與其他方法的濾波效果對(duì)比

Figure 5 The filtering effect with the current method in compa-rison with that with other methods

2.3 彩色圖像濾波

當(dāng)原圖像為彩色圖像時(shí)，直接進(jìn)行Gabor濾波會(huì)使圖像失真. 通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，先將彩色圖像分解為R、G、B分量，對(duì)各分量依次進(jìn)行Gabor濾波，再按原順序重新組合成彩色圖片，使彩色圖像在濾除光干擾的前提下又不失真. 為了驗(yàn)證本文方法的優(yōu)越性，在相同條件下將本文方法與其他濾波方法相比. 通過(guò)高斯低通濾波方法處理后，圖片整體偏亮，且左上方及右下角亮度仍然較高，可見(jiàn)圖像表面的光干擾并沒(méi)有被濾除； 而通過(guò)本文方法處理后的圖片，灰度較均勻，且左上方及右下角的光干擾濾除效果明顯(圖6).

圖6 本文方法與其他方法的濾波效果對(duì)比

Figure 6 The filtering effect with the current method in comparison with that with other methods

圖像灰度的方差可表示圖像的均勻性(表1). 對(duì)于單一灰階圖像，濾波前圖像的灰度方差為80.2，通過(guò)Gabor濾波后灰度方差為48.2(下降了39.9%)，圖像的均勻性僅提升39.9%；通過(guò)灰度化后，灰度方差為128.9(上升了60.7%)，圖像的均勻性降低了60.7%. 通過(guò)高斯低通濾波后，灰度方差為79.9(下降了0.4%)，說(shuō)明圖像的均勻性僅提升0.4%. 對(duì)于彩色圖像，在濾波前圖像的灰度方差為3 000，通過(guò)Gabor濾波后，灰度方差為1 810(下降了39.7%)，說(shuō)明圖像的均勻性提升了39.7%； 通過(guò)高斯低通濾波后，灰度方差為2 886(下降0.4%)，說(shuō)明圖像的均勻性僅提升0.4%. 可見(jiàn)本文的方法可有效濾除EPD圖像采集過(guò)程中引入的光干擾.

表1 圖像灰度值的方差比較

3 結(jié)論

針對(duì)EPD顯示圖像采集易受光干擾的問(wèn)題，設(shè)計(jì)了一種基于Gabor濾波器的圖像濾波系統(tǒng). 對(duì)于受光干擾的單一灰階圖像，先通過(guò)Gabor濾波降低光照變化和噪聲對(duì)圖像的影響，再通過(guò)閉運(yùn)算進(jìn)行邊緣補(bǔ)償；對(duì)于受光干擾的彩色圖像，先提取圖像的R、G、B分量，分別對(duì)各分量進(jìn)行Gabor濾波，最后將3組分量重新組合成彩色圖像. 結(jié)果表明：該系統(tǒng)針對(duì)單一灰階圖像能有效濾除39.9%的光干擾, 針對(duì)彩色圖像能有效濾除39.7%的光干擾. 該系統(tǒng)的不足之處在于Gabor濾波器參數(shù)的設(shè)定仍需要手動(dòng)完成，后期可通過(guò)分析不同光干擾場(chǎng)景的數(shù)據(jù)來(lái)統(tǒng)計(jì)參數(shù)變化規(guī)律.