李 豪， 王德志， 魏永起， 姚 武

（1.寧夏大學(xué) 土木與水利工程學(xué)院，寧夏 銀川 750021； 2.同濟(jì)大學(xué) 先進(jìn)土木工程材料教育部重點(diǎn)實(shí)驗(yàn)室，上海 201804； 3.同濟(jì)大學(xué) 材料科學(xué)與工程學(xué)院，上海 201804）

發(fā)泡混凝土是一種通過(guò)人工方式引入微米至毫米尺度氣孔的多孔材料，因此具有優(yōu)良的保溫、隔熱性能.前期研究發(fā)現(xiàn)，發(fā)泡混凝土中孔結(jié)構(gòu)是影響其傳熱性能的重要因素［1-3］，但其微米至毫米尺度的孔隙無(wú)法使用傳統(tǒng)氮吸附法或壓汞法進(jìn)行精確測(cè)定［4］.因此，研究人員嘗試采用光學(xué)成像技術(shù)表征材料毫米級(jí)的氣孔［5］，但該方法缺少孔結(jié)構(gòu)圖像的分析方法.鑒于此，彭軍芝等［6］使用Matlab對(duì)發(fā)泡混凝土的孔結(jié)構(gòu)圖像進(jìn)行了量化分析，但其圖像中孔的邊緣與基體對(duì)比度較差，很難對(duì)孔結(jié)構(gòu)進(jìn)行精確識(shí)別.張雄等［5］使用基體著色并填充納米微粉的后處理方式，提高了孔壁和孔內(nèi)的對(duì)比襯度，但該方法需要較精細(xì)的后處理工藝，且檢測(cè)精確程度易受操作者制樣水平的影響［7］.因此建立一種既能簡(jiǎn)化制樣流程又可以進(jìn)一步提高對(duì)比襯度，同時(shí)還可以通過(guò)選擇合適圖像識(shí)別算法對(duì)泡沫混凝土孔結(jié)構(gòu)進(jìn)行檢測(cè)和識(shí)別的方法就顯得尤為必要.通過(guò)調(diào)研發(fā)現(xiàn)，在金屬機(jī)械零件表面劃痕檢測(cè)時(shí)引入低角度光源，可明顯提高待測(cè)表面與劃痕間的對(duì)比襯度［8］，但其缺少光源角度與成像之間關(guān)系的分析.這一問(wèn)題在田洪志［9］的研究中得以解決，其深入研究了單側(cè)光源不同角度對(duì)檢測(cè)金屬板表面缺陷的影響，確定了單側(cè)光源照明的最佳角度，但對(duì)于缺陷復(fù)雜的表面檢測(cè)仍不理想.馮超［10］開發(fā)了低適光性表面缺陷的檢測(cè)方法，提出了針對(duì)復(fù)雜表面缺陷的檢測(cè)解決方案.

本文借鑒零度光源在金屬表面缺陷檢測(cè)方面的應(yīng)用，嘗試將其拓展至發(fā)泡混凝土孔結(jié)構(gòu)檢測(cè)領(lǐng)域，以解決發(fā)泡混凝土孔特征檢測(cè)和識(shí)別過(guò)程中繁瑣的樣品后處理工藝以及由此可能帶來(lái)的圖像失真問(wèn)題，為發(fā)泡材料孔特征的深入研究提供快速且精準(zhǔn)的定量表征方法.

1 試驗(yàn)

1.1 設(shè)備

圖像采集設(shè)備為L(zhǎng)apsun的USB500工業(yè)相機(jī)，像素為500萬(wàn)，光學(xué)尺寸為1/2.5彩色；光學(xué)鏡頭選用Lapsun品牌180X-720X型號(hào)；圖像采集軟件為S-EYE軟件，獲得2 048 像素×1 536像素的彩色圖像；零度光源選用Lapsun生產(chǎn)的RL-146-00W，為滿足視 域選定內(nèi) 徑為104 mm.同 時(shí)，采 用IFV-C40DW同軸光源和IDR-F60/20DW的常光源進(jìn)行對(duì)比分析.

零度光源產(chǎn)生的是一種與試件截面夾角極小的光線，當(dāng)光線越過(guò)一側(cè)孔壁時(shí)，在另一側(cè)孔壁只有小部分的光線可以到達(dá)其孔壁上側(cè)區(qū)域，從而形成了一個(gè)與孔形貌相同的“黑體”，使采集的圖像中孔上壁和孔內(nèi)壁具有很強(qiáng)的區(qū)分度［11］，其結(jié)構(gòu)示意圖見圖1（圖中安裝孔公稱直徑和深度均為3 mm）.同軸光源是通過(guò)單向透光玻璃而產(chǎn)生對(duì)界面垂直方向成像的特征光源.

圖1 零度光源結(jié)構(gòu)示意圖Fig.1 Structure diagram of zero-degree light source（size： mm）

1.2 試驗(yàn)方法

1.2.1 試件制備

用密度500 kg/cm3的發(fā)泡混凝土，制備100 mm×100 mm×100 mm的 試 件，并 進(jìn)行切割、打 磨、催粉.具體步驟為：（1）使用鋸條將泡沫混凝土切割；（2）將具有待測(cè)截面的試件使用6.5 μm的極細(xì)砂紙進(jìn)行打磨［5］，打磨時(shí)控制每2～3個(gè)相對(duì)運(yùn)動(dòng)周期清理1次砂紙上的粉末，每個(gè)試件截面的打磨過(guò)程持續(xù)20～30 min；（3）使用壓縮空氣型噴槍，將打磨完成后試件上的殘留粉末催出，該過(guò)程中保持待測(cè)截面與噴槍噴出的流體相對(duì)平行，盡量避免流體對(duì)試件的破壞.

1.2.2 圖像檢測(cè)

傳統(tǒng)圖像檢測(cè)需要對(duì)試件進(jìn)行著色、填粉、再次打磨等后處理，存在對(duì)其孔結(jié)構(gòu)的二次破壞問(wèn)題［11］，將未處理的試件分別使用無(wú)光源（自然光）、常光源、同軸光源以及零度光源對(duì)泡沫混凝土試件的同一截面區(qū)域進(jìn)行采樣，不同光源下泡沫混凝土的圖像見圖2.由圖2可見：與無(wú)光源（自然光）相比，光源的存在使工業(yè)相機(jī)能收集到更多光線，有光源的泡沫混凝土孔形貌特征更加顯著；同軸光源下的泡沫混凝土孔內(nèi)與孔壁的對(duì)比襯度較小，孔邊緣模糊；零度光源下泡沫混凝土孔內(nèi)與孔壁的對(duì)比襯度較大，孔特征更加顯著.

圖2 不同光源下泡沫混凝土的圖像Fig.2 Images of foamed concrete under different light sources

2 結(jié)果與討論

2.1 孔結(jié)構(gòu)視覺圖像的最小取樣數(shù)

發(fā)泡混凝土發(fā)泡時(shí)容易受到環(huán)境溫度及發(fā)泡劑分散程度等變量的影響，成型后孔的特征分布是不均衡的［6］，因此需要在孔特征統(tǒng)計(jì)之前，確定隨機(jī)獲取圖像的最小取樣數(shù)n.當(dāng)取樣數(shù)足夠大時(shí)（n≥30），樣品均值近似服從均值為μ、方差為σ2的正態(tài)分布，因此最小取樣數(shù)n實(shí)際上是在某置信程度及誤差程

度下對(duì)樣品均值的區(qū)間估計(jì)［5］，即：

式中：λ為分布參數(shù)，其值由置信程度確定，一般情況下混凝土的置信程度為95%，可查得當(dāng)置信程度為95%，n=30時(shí)，λ為1.96［12］；d為 估 計(jì) 誤 差，一 般 取d≤15%［5］.

隨機(jī)取2組各30張發(fā)泡混凝土的圖像，用Image-Pro Plus軟件進(jìn)行分析，計(jì)算平均孔隙率（體積分?jǐn)?shù)）的最小樣本數(shù)，結(jié)果見表1.由表1可見，2組發(fā)泡混凝土的n分別為0.145 7、0.159 9，取整得到n=1，即只需1個(gè)樣本圖像即可在保證率為95%的情況下，滿足待統(tǒng)計(jì)量估計(jì)誤差小于15%.

表1 最小樣本數(shù)Table 1 Minimum sample amount

2.2 不同制樣方法增強(qiáng)孔結(jié)構(gòu)對(duì)比襯度效果對(duì)比

ASTM-C457-06《Standard test method for microscopical determination of parameters of the air-void system in hardened concrete》根 據(jù) 體 視 學(xué) 原理，使用線法和網(wǎng)格法對(duì)拍攝圖像的平均孔隙率進(jìn)行統(tǒng)計(jì)，而后研究人員對(duì)制樣和拍攝手法進(jìn)行改進(jìn)［13-14］，設(shè)法增加孔內(nèi)壁與孔壁之間的對(duì)比度，使圖像中孔的特征更加顯著.由文獻(xiàn)［6］可知，在試樣制備的基礎(chǔ)上對(duì)其進(jìn)一步處理，宏觀處理發(fā)泡混凝土孔結(jié)構(gòu)的傳統(tǒng)方法有2種，分別為將孔壁著色和孔壁著色-孔內(nèi)填充［5］.3種處理方法的制樣工序流程對(duì)比見圖3.由圖3可見：零度光源法共需要47～65 min；與孔壁著色、孔壁著色-孔內(nèi)填充相比，零度光源法分別節(jié)約了10～15、20～30 min.傳統(tǒng)方法若使用掃描電鏡（SEM）識(shí)別孔隙，在制樣工序后還需要進(jìn)行鍍金（30 s）、裝樣（2 min）、抽真空（3 min）以及對(duì)焦掃描（10 min）.傳統(tǒng)方法在著色階段要求顏色涂抹均勻以及著色材料不能與孔內(nèi)壁接觸，且填充后剝離白色粉末時(shí)，要求刮片不能破壞孔壁，因此其后處理對(duì)操作要求較高；而零度光源輔助拍攝的方法，省略了后處理部分，整個(gè)流程簡(jiǎn)便快速.

圖3 3種處理方法的制樣工序流程對(duì)比Fig.3 Comparison of sample preparation process of three processing methods

3種處理方法的采樣結(jié)果見圖4.由圖4可見：孔壁著色法的位置A（見圖4（g））有墨水沁入孔壁，將孔內(nèi)部分區(qū)域著色，導(dǎo)致孔的特征出現(xiàn)偏差；而零度光源法呈現(xiàn)出的孔內(nèi)部和孔壁區(qū)分完整，輪廓清晰；孔壁著色-孔內(nèi)填充法在位置B、C處（見圖4（i））均出現(xiàn)了孔相連通的狀況，原因可能是發(fā)泡混凝土中孔隙率偏大，相鄰孔的孔壁較薄，在使用孔壁著色-孔內(nèi)填充法刮去粉末的時(shí)侯，使孔壁遭到了破壞.

圖4 3種處理方法的采樣結(jié)果Fig.4 Sampling results of three processing methods

2.3 孔結(jié)構(gòu)圖像優(yōu)化

在進(jìn)一步提取圖像信息前，一般需要對(duì)原始圖像進(jìn)行降噪等優(yōu)化處理，而判斷圖像優(yōu)劣的標(biāo)準(zhǔn)是圖像對(duì)目標(biāo)特征的表達(dá)程度［15］.發(fā)泡混凝土的目標(biāo)特征是孔的邊緣特征以及截面孔的總體特征.首先使用灰度化方法將原始圖像變?yōu)閱瓮ǖ佬畔⒕仃?，其原理是將圖像相同位置的不同通道加權(quán)平均，得到［0-255］的單通道信息［16］.采用人眼對(duì)紅綠藍(lán)顏色的敏感程度進(jìn)行加權(quán)平均［16］，得到灰度圖像g(x，y)：

式中：Red（x，y）、Green（x，y）、Blue（x，y）分別為紅色、綠色、藍(lán)色通道像素矩陣；x、y分別為矩陣行、列的位置.

3種處理方法得到的發(fā)泡混凝土灰度化、降噪處理效果圖見圖5.由圖5可見，灰度化處理后，在孔上壁中有一些點(diǎn)狀和塊狀噪點(diǎn).將大鄰域的中值濾波［17］得到的圖像u（x，y）作為噪聲，取原圖的灰度圖像與大鄰域中值濾波圖像的差值，并將差值小于0的像素替代為0，得到降噪后的圖像γ（x，y）：

該方法本質(zhì)上是自動(dòng)依據(jù)圖像的整體特征選取全局圖像閾值.由圖5還可見：孔壁著色法灰度化后在A、B位置有明顯的由于著色導(dǎo)致的反光，因此在灰度圖中位置A、B處的像素特征與孔內(nèi)像素特征幾乎相同；孔壁著色-孔內(nèi)填充法在A位置存在無(wú)法去除的填充粉末，導(dǎo)致其像素特征與孔內(nèi)相同；零度光源法經(jīng)過(guò)灰度處理后，孔壁上出現(xiàn)噪聲點(diǎn)、塊，經(jīng)過(guò)降噪處理后噪聲明顯減小.

圖5 3種處理方法得到的發(fā)泡混凝土灰度化、降噪處理效果圖Fig.5 Effect figures of gray and noise reduction of foamed concrete by three processing methods

2.4 不同制樣方法孔結(jié)構(gòu)特征統(tǒng)計(jì)對(duì)比

不同的制樣方法得到的孔結(jié)構(gòu)圖像有明顯的差異，本節(jié)分別從截面的孔邊緣特征、孔徑分布、孔隙率等維度進(jìn)行更深入的探究.分別采用Laplacian、Kirsch及Canny算 子［18］對(duì) 孔邊緣進(jìn)行 對(duì)比分析，結(jié)果見圖6.由圖6可見：Canny算子對(duì)孔邊緣的提取較好；與Canny算子相比，Laplacian算子和Kirsch算子處理后圖像的孔邊緣雜亂不清晰，存在更多不連續(xù)的“偽邊”且孔輪廓模糊；孔壁著色法用Canny算子處理（見圖6（g））后，由于著色部分反光，導(dǎo)致了灰度圖像出現(xiàn)“偽孔”及“偽邊”；孔壁著色-孔內(nèi)填充法用Canny算子處理后（見圖6（h）），孔壁較薄部分的孔邊緣輪廓線相連接，整體成為了1個(gè)孔隙；零度光源法用Canny算子處理后圖像的孔邊緣清晰.

用Image Pro Plus軟件計(jì)算發(fā)泡混凝土的孔徑分布和面積率，結(jié)果見圖7.由圖7可見：孔壁著色法共有1 536個(gè)孔，平均孔隙率為52.73%，半徑為0～90 μm的 孔占比81.58%，由于 孔壁著色時(shí)墨沁入了孔內(nèi)，從而造成原始圖像的邊緣不顯著，將單個(gè)的孔統(tǒng)計(jì)為多個(gè)破碎狀的小孔（見圖6（g）、圖7（g））；孔壁著色-孔內(nèi)填充法共有967個(gè)孔，平 均 孔 隙 率 為62.52%，半 徑 為10～90 μm的孔占比76.42%，與孔壁著色法相比破碎狀小孔數(shù)目明顯下降，但其最大的孔半徑是876 μm，這是因?yàn)槠浯罂讈?lái)源于多個(gè)孔的聯(lián)通，即相鄰的薄壁孔被統(tǒng)計(jì)為1個(gè)孔（見圖6（h）、圖7（h））；零度光源法共有677個(gè)孔，平均孔隙率為43.86%，其沒有出現(xiàn)大量破碎狀小孔以及由于薄壁而呈1個(gè)聯(lián)通孔的情況（見圖7（i））.綜上，零度光源法采集的圖像孔隙邊緣清晰且薄壁保持完整，可以觀察半徑為10～900 μm孔的形貌特征以及統(tǒng)計(jì)整體特征.

圖6 3種處理方法得到的發(fā)泡混凝土孔邊緣識(shí)別對(duì)比Fig.6 Comparison of pore edge recognition of foamed concrete by three processing methods

圖7 3種處理方法得到的發(fā)泡混凝土孔徑分布和面積率Fig.7 Pore size distribution and area percentage of foamed concrete by three processing methods

3 結(jié)論

（1）針對(duì)發(fā)泡混凝土孔測(cè)試中的制樣問(wèn)題，選用零度光源法，以非接觸式形式增加了孔壁與孔內(nèi)對(duì)比度，與傳統(tǒng)方法相比，零度光源法對(duì)孔壁的破壞小，孔特征保存完整.

（2）通過(guò)統(tǒng)計(jì)學(xué)方法確定了在固定尺度下發(fā)泡混凝土孔結(jié)構(gòu)分析所需的最小樣本數(shù)，從而增強(qiáng)了零度光源法在統(tǒng)計(jì)學(xué)意義上的科學(xué)性.

（3）優(yōu)化的均值濾波算法可以將零度光源法獲取圖像中的噪聲去除，并且不影響孔隙圖像的細(xì)部特征，從而得到相對(duì)清晰且完整的孔特征圖像.

（4）針對(duì)零度光源輔助拍攝的發(fā)泡混凝土孔結(jié)構(gòu)圖像，使用Canny算法，可以精確地提取發(fā)泡混凝土孔邊緣圖像特征.