董天驕， 于海濱， 安春雷， 史劍光， 江 曉

(1.杭州電子科技大學(xué) 電子信息學(xué)院，浙江 杭州 310018；2.中國極地研究中心，上海 200000)

0 引 言

冰芯光學(xué)檢測技術(shù)在保持冰芯固態(tài)下掃描冰芯，對冰芯剖面或冰芯切片的照片進(jìn)行分析，獲得冰芯的光學(xué)性質(zhì)，如灰度，晶體結(jié)構(gòu)等[1]。

目前，國外已經(jīng)研制出多套冰芯光學(xué)測定儀。日本在2003年就已經(jīng)研制出了一種冰芯光學(xué)測定儀，該冰芯光學(xué)測定儀應(yīng)用于Dome Fuji ice core project Ⅱ[2]。美國國家冰芯實驗室在2007年也研制出了冰芯光學(xué)測定儀[3]。該冰芯光學(xué)測定儀的掃描分辨率為0.05～1 mm。具有自動光學(xué)圖像采集和數(shù)據(jù)處理的功能。然而，這些設(shè)備只能在實驗室使用，無法抵抗極地的惡劣環(huán)境，難以在第一時間保存冰芯的原始信息[4]。

國內(nèi)光學(xué)影相技術(shù)研究在極地科研領(lǐng)域起步較晚[5]。國內(nèi)的光學(xué)影像儀應(yīng)用在冰芯檢測領(lǐng)域幾乎沒有成熟產(chǎn)品[6]，因此，本文設(shè)計并制作了冰芯光學(xué)特性二維重構(gòu)儀。儀器基于機(jī)電一體式設(shè)計，并引入全景圖像拼接技術(shù)，可實現(xiàn)自動光學(xué)掃描、冰芯全景圖重建等功能。同時，模塊化設(shè)計和智能變頻恒溫控制等關(guān)鍵技術(shù)的引入，使得本儀器可滿足冰川或極地惡劣環(huán)境的現(xiàn)場工作要求。最后，從控制精度、合成效果兩個方面進(jìn)行實驗分析，說明本儀器能夠穩(wěn)定可靠地應(yīng)用于極地科研工作。

1 系統(tǒng)設(shè)計

冰川或極地氣候惡劣，運(yùn)輸不便,因此，本儀器采用分模塊設(shè)計，便于運(yùn)輸和組裝。儀器包括工業(yè)相機(jī)、電機(jī)滑軌、發(fā)光臺、工作臺、電機(jī)控制器五部分。工業(yè)相機(jī)用于冰芯拍照，電機(jī)滑軌用于控制工業(yè)相機(jī)，發(fā)光臺用于冰芯打光。電機(jī)控制器用于控制電機(jī)。工作臺用于安裝所有模塊。工業(yè)相機(jī)和電機(jī)控制器裝在智能變頻恒溫箱中，用于保障在低溫環(huán)境可靠工作。系統(tǒng)結(jié)構(gòu)如圖1所示，計算機(jī)控制相機(jī)拍攝冰芯圖片，通過串口向電機(jī)驅(qū)動模塊發(fā)送前進(jìn)、后退和停止指令。計算機(jī)完成二維重構(gòu)的任務(wù)。電機(jī)驅(qū)動模塊控制步進(jìn)電機(jī)的啟停和方向。智能變頻恒溫箱保護(hù)系統(tǒng)的內(nèi)部電路。發(fā)光二極管(light-emitting diode,LED)燈負(fù)責(zé)給冰芯、標(biāo)尺打光。電源負(fù)責(zé)給整個系統(tǒng)供電。儀器結(jié)構(gòu)如圖2所示。

圖1 系統(tǒng)結(jié)構(gòu)

圖2 冰芯掃描機(jī)床結(jié)構(gòu)

2 測控系統(tǒng)

測控系統(tǒng)由自動光學(xué)掃描、圖像拼接兩部分組成。自動光學(xué)掃描是采用步進(jìn)電機(jī)和步進(jìn)電機(jī)控制器組成的一體化控制系統(tǒng)。自動光學(xué)掃描需要設(shè)定電機(jī)的運(yùn)動方向和運(yùn)動距離，從而滿足前后兩張圖像的重疊率；控制工業(yè)相機(jī)拍照，將照片保存到本地；掃描結(jié)束后自動返回初始位置。

執(zhí)行圖像拼接的操作，將拍攝的多張圖片合成一張完整的全景圖。系統(tǒng)流程如圖3所示。

圖3 系統(tǒng)流程

3 圖像拼接

實際工作中，冰芯樣品長度1 m左右，獲得一幅場景的全景圖像很難。如果擴(kuò)大圖像采集設(shè)備的視角，會導(dǎo)致圖像有很大的畸變，同時分辨率沒有提高。為了解決這個問題，本文采用基于特征點(diǎn)匹配的圖像拼接算法。具體拼接流程如圖4所示。

圖4 基于特征點(diǎn)的圖像拼接流程

3.1 圖像預(yù)處理

拍攝的冰芯照片一般偏亮，因此，冰芯的原始圖片灰度主要集中在比較亮的區(qū)間。采用直方圖均衡化的方法將圖像進(jìn)行非線性拉伸，重新分配圖像像素值，使一定灰度范圍內(nèi)的像素數(shù)量大致相同，以增強(qiáng)局部的對比度而不影響整體的對比度，使圖像更加清晰。

3.2 圖像匹配區(qū)域定位

為了減少圖像拼接的運(yùn)算時間，需要確定兩幅圖像中配準(zhǔn)的區(qū)域。因為本儀器是水平掃描，所以在圖像配準(zhǔn)時，只需在相似區(qū)域內(nèi)進(jìn)行特征點(diǎn)的提取和配準(zhǔn)，從而縮短了圖像拼接的時間。

圖像匹配的時前后圖像的重疊率在30 %左右效果較好，考慮到步進(jìn)電機(jī)在運(yùn)動中存在丟包的現(xiàn)象，行進(jìn)距離不是非常精確。因此，本文選擇前后圖片40 %的區(qū)域作為圖像匹配區(qū)域，如圖5所示。

圖5 圖像匹配區(qū)域定位

3.3 圖像特征提取

確定圖像匹配區(qū)域后，就要在匹配區(qū)域內(nèi)進(jìn)行特征提取和匹配。因為本文中采集的圖像只存在平移變換，為了加快檢測速度，采用FAST特征提取。FAST角點(diǎn)的定義為:若某像素與其周圍鄰域內(nèi)足夠多的像素點(diǎn)相差較大，則該像素可能是角點(diǎn)

(1)

式中I(x)為圓周上任意一點(diǎn)的灰度，I(p)為圓心的灰度，εd為灰度差的閾值,N為滿足不等式的個數(shù)。如果N大于給定閾值，一般為周圍圓圈點(diǎn)的3/4，則認(rèn)為該像素是1個特征點(diǎn)。

3.4 圖像特征匹配

為了提高圖像匹配精度，采用RANSAC算法對變換矩陣進(jìn)行求解與精煉。

3.5 圖像融合

圖像匹配后，圖像間的變換關(guān)系唯一確定。但由于拍攝時光照、環(huán)境條件的差異，兩幅待拼接圖像的重疊區(qū)域可能會有較大的差異。如果直接對圖像進(jìn)行疊加拼接，則在拼接位置上會存在明顯的拼接縫以及重疊區(qū)域的模糊和失真的現(xiàn)象。

本文采用加權(quán)平均的融合方法進(jìn)行圖像平滑過渡。假設(shè)f1,f2是兩幅待拼接的圖像，f是融合后的圖像，則

(2)

式中d1,d2為權(quán)重值。在重疊區(qū)域中，d1由1漸變到0，d2由0漸變到1，由此實現(xiàn)了在重疊區(qū)域無縫拼接。

4 實驗與結(jié)果分析

實驗計算機(jī)采用Windows 7 64位系統(tǒng)、Intel i3處理器、8GB內(nèi)存；網(wǎng)絡(luò)工業(yè)相機(jī)分辨率為1 920像素×1 080像素。系統(tǒng)采用Microsoft Visual Studio 2013的開發(fā)環(huán)境，C#編寫PC端軟件，冰芯采用純凈水凍結(jié)而成的冰磚樣品。

1)步進(jìn)電機(jī)控制精度分析

控制方式為開環(huán)工作，選用6，12，16 cm，3個設(shè)定距離來分析運(yùn)行精度，每個重復(fù)3次，得到結(jié)果如表1所示。

表1 步進(jìn)電機(jī)運(yùn)行距離的誤差

從結(jié)果上可以看出，由于步進(jìn)電機(jī)存在丟幀的現(xiàn)象，步進(jìn)電機(jī)的實際運(yùn)行距離都會比設(shè)定距離小一點(diǎn)。但運(yùn)行誤差不超過10 %，不會對前后兩張圖片的重疊區(qū)域產(chǎn)生太大影響。因此，誤差可以接受。

2)冰芯圖像拼接分析

本文在4種亮度環(huán)境中對冰芯進(jìn)行圖像拼接，拼接后的冰芯全景如圖6所示,可以看出,拼接后的圖片完整的保存了整條冰芯，整幅圖細(xì)節(jié)清晰，過渡自然，達(dá)到設(shè)計要求。

圖6 冰芯全景

5 結(jié) 論

設(shè)計的冰芯光學(xué)特性二維重構(gòu)儀實現(xiàn)了光學(xué)掃描、重建冰芯全景等功能。本儀器全程自動控制，有效減輕了科研人員的負(fù)擔(dān)；重建后冰芯全景圖分辨率高達(dá)7 200像素×1 080像素，有效保存了冰芯的細(xì)節(jié)信息。儀器工作穩(wěn)定可靠，滿足科研工作的要求，為冰芯科研工作提供一種便捷高效的現(xiàn)場分析工具。