趙昆越,田漢民,郭 丹,王 錚,常衛(wèi)洪
(1.河北工業(yè)大學(xué)電子信息工程學(xué)院,天津 300401;2.天津市電子材料與器件重點(diǎn)實(shí)驗(yàn)室,天津 300401)
目前,在接觸角測(cè)量領(lǐng)域,應(yīng)用最廣泛,測(cè)量精度和效率最高的方法是圖像外形分析法[1],其原理為將液滴滴于固體樣品臺(tái),通過顯微鏡頭和相機(jī)拍攝液滴圖像,應(yīng)用數(shù)字圖像處理技術(shù)和相應(yīng)的測(cè)量算法計(jì)算接觸角。常用接觸角測(cè)量算法包含量高法[2]、量角法[3]、擬合法[4]、L-Y法等。其中量高法、量角法、擬合法基于一定的數(shù)學(xué)模型,即把液滴近似為球型的一部分,通過測(cè)量幾何參數(shù),例如寬度、高度、切線斜率等計(jì)算接觸角。
液滴圖像的拍攝質(zhì)量對(duì)于接觸角測(cè)量精度的影響很大。很多時(shí)候由于人為不正當(dāng)操作和硬件設(shè)備(樣品臺(tái),背景光源,攝像頭)的優(yōu)劣會(huì)導(dǎo)致測(cè)量誤差的產(chǎn)生,目前市場(chǎng)中不同廠家的接觸角測(cè)量?jī)x硬件的類別和規(guī)格不同,精度也不同。造成液滴圖像質(zhì)量差的因素很多,包含光源誤差、機(jī)械結(jié)構(gòu)定位誤差、環(huán)境誤差等。
精確定位固-液-汽三相交點(diǎn),尋找液滴邊界,對(duì)于接觸角測(cè)量的準(zhǔn)確性至關(guān)重要。液滴圖像呈現(xiàn)高亮度特點(diǎn),背景光源的作用是讓待測(cè)液滴和背景之間形成高明暗對(duì)比度,凸顯液滴輪廓,便于后續(xù)邊緣檢測(cè),光照強(qiáng)度分布對(duì)圖像的亮度和對(duì)比度起到關(guān)鍵作用,不均勻光照強(qiáng)度會(huì)使圖像傳感器曝光不足或過飽和[5],導(dǎo)致圖像局部過亮或過暗,丟失部分細(xì)節(jié),進(jìn)而產(chǎn)生測(cè)量誤差。
本文重點(diǎn)研究了不均勻背景光照對(duì)于接觸角測(cè)量誤差的影響,實(shí)驗(yàn)通過改變不均勻光照強(qiáng)度,分別采用量高法、量角法、五點(diǎn)擬合法測(cè)量接觸角,記錄并分析數(shù)據(jù),發(fā)現(xiàn)光照強(qiáng)度和接觸角測(cè)量誤差之間的關(guān)系,應(yīng)用Otsu二值化[6]、Canny邊緣檢測(cè)、三次樣條插值算法擬合邊緣,從圖像角度分析總結(jié)了圖像液滴邊緣隨光照改變的變化規(guī)律,并通過Mie光散射理論闡述其原理。
接觸角測(cè)量實(shí)驗(yàn)簡(jiǎn)易模型如圖1所示,1-CCD攝像頭(成像分辨率2 048×1 536),2-微量進(jìn)樣器,3-H2O2/H2O溶液液滴,4-載物臺(tái)(搭載材料為拋光后的硅片),5-背景光源(采用的是led單色平行冷光源)。
圖1 接觸角測(cè)量系統(tǒng)簡(jiǎn)易模型
首先在控制溫度不變的情況下,通過調(diào)節(jié)光源亮度改變液滴圖像光照強(qiáng)度分析液滴圖像灰度曲線,為盡量避免液滴蒸發(fā)帶來的影響,采用快速抓拍功能,每幀時(shí)間間隔1.5 s,如圖2所示,從1到5對(duì)液滴光照逐漸增強(qiáng)。液滴圖像1~5灰度分布如圖3所示。其中液滴與搭載液滴的硅片所占灰度級(jí)大致分布在0~75范圍內(nèi),背景光照灰度級(jí)分布在75~255范圍內(nèi),隨著光照強(qiáng)度增強(qiáng),見(1)處,液滴灰度峰值變小且向右逐漸偏移。(2)處峰的存在是由于光照較暗時(shí),在圖片1背景中產(chǎn)生的灰暗部分。
圖2 不同光照強(qiáng)度下的液滴圖像
圖3 不同亮度下液滴圖像的灰度直方圖
從灰度直方圖角度分析,假設(shè)一個(gè)像素點(diǎn)代表單位面積,統(tǒng)計(jì)1~5液滴圖像中灰度值分布在0~75之間的像素點(diǎn)數(shù)量,如表1所示,隨著光照增強(qiáng),灰度值分布在0~75的像素點(diǎn)數(shù)量逐漸減少,換言之,在圖像中液滴和硅片所占面積逐漸減小。
表1 在不同光照強(qiáng)度的圖像中0~75灰度值像素點(diǎn)數(shù)量
從直觀圖像角度來看,將1~5圖片首先應(yīng)用Otsu二值化方法進(jìn)行閾值分割[7],進(jìn)一步應(yīng)用canny邊緣檢測(cè)提取液滴輪廓,并將5張邊緣檢測(cè)后的圖像疊加,如圖4所示。隨著光照從1到5逐漸增強(qiáng),圖像中液滴邊緣由外向內(nèi)逐漸收縮,硅片水平面逐漸下降??梢钥闯鲭S著光照增強(qiáng),圖像中液滴逐漸發(fā)生形變,其所占面積比例逐漸減小。
圖4 不同光照強(qiáng)度的圖像中液滴邊緣疊加
圖像中液滴的形變必然帶來接觸角測(cè)量的誤差,然而其形變的原因是光在傳播過程中,光子與液滴邊界水分子相互作用發(fā)生散射,如圖5所示,光波被均勻非吸收介質(zhì)的球狀顆粒散射可由Mie散射理論形式給出[8],如式(1)反映了距離散射點(diǎn)r處的散射光強(qiáng),式(2)中振幅函數(shù)S1(θ)、S2(θ)表達(dá)式見式(3)和式(4),其中an和bn是與貝爾賽爾函數(shù)和漢克爾函數(shù)有關(guān)的函數(shù)。
圖5 光在液滴邊界散射示意圖
(1)
式中:θ為散射角;Isca為散射光強(qiáng);I0為入射光強(qiáng);λ為入射光波長(zhǎng);φ為偏振角。
I(θ,φ)=|S1(θ)|2sin2φ+|S2(θ)|2cos2φ
(2)
(3)
(4)
從式(1)中可以看出,當(dāng)波長(zhǎng)λ、散射角θ、偏振角φ以及散射距離r固定時(shí),散射光強(qiáng)Isca與入射光強(qiáng)I0成正比關(guān)系,同時(shí)文獻(xiàn)[9]中指出散射光強(qiáng)Isca與散射角θ的關(guān)系,見圖6,在顆粒直徑x和折射率m不變的情況下,隨著散射角θ增大,散射光強(qiáng)Isca逐漸減小,在0~10°范圍內(nèi)Isca下降率最高,在θ趨近于0°時(shí),Isca最大,幾乎等于入射光強(qiáng),Isca/I0趨近于1。
圖6 散射角與散射光強(qiáng)的關(guān)系
當(dāng)光源照射液滴在CCD上成像時(shí),圖像中某一點(diǎn)像素的灰度值h(X,Y)只與CCD感光面光強(qiáng)分布及變化有關(guān),如表達(dá)式(5),其中I(x,y)表示單位像素上光照平均值。
(5)