刁福林，陳 昊，李延豐，陳 曦，楊立余

(1. 南昌航空大學(xué) 信息工程學(xué)院，南昌 330063； 2. 無損檢測技術(shù)教育部重點實驗室(南昌航空大學(xué))，南昌 330063；3. 青島啤酒股份有限公司，山東 青島 266000)

0 引言

大多數(shù)金屬材料是由許多小晶粒組成，晶粒的大小影響著金屬的拉伸強度、韌性以及耐久度等常規(guī)力學(xué)性能[1]。晶粒圖像分析是檢測金屬材料性能的重要手段之一。因此，為了方便測量晶粒尺寸，準確而有效地提取出晶界非常重要。

現(xiàn)有的晶界提取技術(shù)中采用的Canny、Sobel和Kirsch等邊緣檢測算子雖然能夠有效定位到邊緣，但是容易受到諸多因素的干擾。文獻[2]使用Canny算子和Otsu算法相結(jié)合對現(xiàn)場采集的圖像進行實時處理；文獻[3]用中值濾波方法得到的閾值來代替原始的Sobel算子的閾值，同時利用非極大值抑制的方法對圖像作細化處理，從而提高對圖像紋理的識別率；文獻[4]改進傳統(tǒng)的Kirsch邊緣檢測方法，首先利用模糊數(shù)學(xué)方法消除圖像的背景隨機噪聲，得到了濾波圖像，最后用局部最大值圖像和利用最大熵算法自動選取的閾值得到邊緣圖像；文獻[5]提出用混度蟻群算法計算信息素矩陣得到閾值，進而獲得豐富的邊緣信息，能夠快速定位到圖像的邊緣并且較完整的提取出邊緣；文獻[6]將數(shù)學(xué)形態(tài)學(xué)中的腐蝕、膨脹等運算的組合，并依據(jù)不同的結(jié)構(gòu)元，提出了一種抗噪性能良好的邊緣檢測算法，有效地提取出圖像輪廓。

然而，這些傳統(tǒng)方法對高溫合金的晶界提取都難以得到理想的效果。由于目標金屬在制備和數(shù)據(jù)采集過程中，晶粒本身內(nèi)部會產(chǎn)生偽邊界，也會丟失部分晶界，進而在采集的金相圖像中往往存在劃痕、晶界模糊不清、晶粒內(nèi)部灰度分布不均勻等情況，這會干擾對晶粒的計算和分析。

用于晶界提取的經(jīng)典方法受噪聲影響嚴重，在復(fù)雜的圖像中難以取得理想的效果。為了獲取清晰完整的晶界，對存在的問題進行研究，提出先用拉普拉斯銳化進行處理，再做直方圖均衡化，去除圖像表面上的劃痕，用Otsu法進行分割，然后將處理后的圖像作最大連通域，最后做平滑等處理，從而得到較完整的晶界。由提取到的晶界圖計算出每一幅金相圖像晶粒的個數(shù)和對應(yīng)的面積等晶粒參數(shù)。

1 簡單晶界提取流程

金相圖像主要描繪晶粒和晶界，除了形狀外，每顆粒子都是相似的[5]。因為各個晶粒之間結(jié)構(gòu)的相互作用力不同，相鄰晶粒的邊界并不總是那么明顯，而且一部分的晶?？赡艽嬖谶B接它們的連接通道，使用傳統(tǒng)分割算法不能準確有效地對金相圖像進行分割。金相圖像中常見的干擾項有球狀析出相、平直晶界和劃痕等干擾項(圖1)。

這項研究的傳統(tǒng)算法流程如圖2所示。傳統(tǒng)的晶界提取方法是對圖像作預(yù)處理，去除大部分噪聲，然后利用金典的邊緣提取算法進行提取，最終獲得圖像的邊緣，這對應(yīng)一般的輪廓圖像會有較好的效果[7]。但對于類似于金相圖像的復(fù)雜圖像，金典邊緣提取算法提取的邊緣輪廓粗糙，抗噪聲能力弱且容易受到光照不均的影響，會引入虛假的邊緣分支[8]。因此，最后獲得的晶界圖像不可靠，不利于對晶粒的分析。

2 基于多種方法結(jié)合的晶界提取研究

2.1 圖像預(yù)處理

晶粒圖像在獲取過程中，由于光照等原因，導(dǎo)致圖像部分區(qū)域亮度不均勻，不利于晶界提取?？烧{(diào)節(jié)圖像的對比度和局部亮度實現(xiàn)局部灰度非線性變換，以增強晶界的對比度[6]。但晶界內(nèi)部的黑點邊界同時也被增強，需要進一步處理。

中值濾波有較強的降噪能力，可以起到平滑圖像的作用，還能盡可能的保留邊緣信息[9]。由于金相圖像中存在很多較小的黑點以及椒鹽噪聲，在平滑處理時會影響計算效果，中值濾波會直接忽略這點。選擇3×3像素進行平滑，并且與其他平滑方法相比，降低了一定的模糊程度。

2.2 拉普拉斯銳化

由于金相圖像在獲取、傳輸?shù)冗^程中受到傳遞函數(shù)的衰減影響，使得其中部分晶界有不同程度的模糊情況。為了增強晶界的邊緣和細節(jié)，用拉普拉斯銳化使圖像的灰度反差增大，加強晶界的邊緣區(qū)域，減弱非邊緣區(qū)域，并且保留了原圖的背景，相比于原圖像清晰度得到了明顯增強。其定義如式(1)所示，式中2f表示二階微分，f(x,y) 為像素位置。拉普拉斯可用過零點來確定突變的邊緣，或者是緩慢變換邊緣線的位置。它對孤立的點和線較為敏感，適用于邊緣檢測。

(1)

將計算得到的值替換原來點(x,y)出的灰度值，可確定晶界點的位置，然后由式(2)得到銳化圖像h(x)。

(2)

2.3 直方圖均衡化

經(jīng)過拉普拉斯銳化處理后，圖像的整體清晰度得到了提高，但是在一些位置依然不夠清晰，因此，使用直方圖均衡化來增加圖像局部對比度，其映射方法如式(3)所示。該方法是將圖像灰度直方圖中一個比較集中的灰度區(qū)間重新分配其他灰度區(qū)間的灰度值，使得整體灰度均勻化，實現(xiàn)對圖像非線性拉伸，有利于下一步的處理。

(3)

式中，n表示所有像素的和，ni為當前灰度級下像素數(shù)目，L表示可能存在的灰度級數(shù)。

2.4 晶界提取

Otsu法是按照圖像的灰度特點，將圖像分為背景部分和輪廓部分，分割效果由類間方差來決定：

w=m0m1(μ0-μ1)2

(4)

式中：m0為輪廓的像素點占整個圖像的比值，平均灰度記為μ0；m1表示背景區(qū)域像素點占整幅圖像的比值，平均灰度為μ1。

計算圖像的直方圖并歸一化，根據(jù)經(jīng)驗值假設(shè)分割閾值為h，通過遍歷整張圖像，最終求得閾值k使類間方差w取得最大值，記下此時的閾值，即為最佳閾值。此方法對低對比度和亮度有一定的魯棒性，計算簡單。

用該方法將均衡化后的圖像進行分割提取出晶界，然后對每個連通區(qū)域進行標記，并且求出對應(yīng)的大小，再求最大連通域的索引，最后獲取最大連通域的圖像：

f(i)=Fi-Fi-Δ/Fi-Δ

(5)

式中：Fi表示第i個連通域的面積；Δ表示微小變化量，一般取1；f(i)為在i處的Fi變化率，當f(i)取局部極小值時，那么Fi就為面積最大穩(wěn)定區(qū)域。

由上述方法處理過的晶界圖像邊緣處存在毛刺以及部分晶界線未連接的情況，直接用這些圖像來計算晶粒的個數(shù)和對應(yīng)面積會有較多的誤差，因此，有必要進行進一步平滑細化處理。

使用開運算方法對其進行處理?？梢杂脕硐∥矬w、與晶界分離的細小點，在不明顯改變晶粒面積的情況下能夠平滑晶界邊緣。再用閉運算處理未連接的鄰域晶界和填充微小的空洞，平滑金相圖像的晶界同時不影響計算晶粒的面積。最后再對平滑后的金相圖像作細化去枝處理，細化之前邊緣處有凸起的位置，在細化后得到的圖像會存在一些短線分支需要去除。根據(jù)短線比較短的特點，將閾值設(shè)置在4～6范圍內(nèi)，如果線條的尺度小于閾值，則認為是多余的線條，應(yīng)該去除；反之，如果線條尺度大于該閾值，則視為原圖像中本來存在的線條，不用去除。

通過對上述處理的每一個步驟的詳細描述，現(xiàn)將整個晶界提取的方法流程總結(jié)如下：

1)輸入金相圖像，對圖像進行調(diào)節(jié)亮度，實現(xiàn)局部灰度拉伸，增強晶界的對比度。

2)考慮到第一步在灰度拉伸的同時圖像中的噪聲和孤立黑點得到了增強，需要及時處理，利用3×3像素的中值濾波去除圖像中的小黑點以及椒鹽噪聲，在一定程度上平滑了圖像。

3)步驟2獲得的去噪后圖像中邊緣依然不夠清晰，根據(jù)拉普拉斯對輪廓信息有較好的增強作用的優(yōu)點，對圖像作拉普拉斯銳化處理，進一步增強晶界區(qū)域，弱化背景區(qū)域，使得邊緣和細節(jié)部分更清晰。

4)為了更好地分隔出完整的晶界，對步驟3得到的圖像再作直方圖均衡化，提升依然不夠清晰的區(qū)域，使整幅圖像的灰度均勻化。

5)通過遍歷圖像，計算出使類間方差取最大值的最佳閾值，使用該閾值對處理完的圖像進行分割，計算并標記各連通域，獲得最大連通域的索引，得到最大連通域圖像。

6)對晶界邊緣處有許多毛刺、部分區(qū)域有虛連接的情況，再對分割完的圖像進行平滑細化處理，最終獲得完整晶界圖像。

3 實驗與分析

為了獲得金相圖像，對GH4169、GH738和GH706材料進行研磨、拋光和腐蝕，然后將其置于金相顯微鏡下獲取實驗樣本來驗證方法的提取性能。實驗所用的圖片均是通過圖像顯微采集系統(tǒng)采集并整理獲得。

3.1 晶界提取對比試驗

以一幅具有代表性的GH738金相圖像為例(圖3a)，用不同方法進行晶界提取實驗分析。為了便于觀察，放大了每個處理圖像中的相同區(qū)域。圖3b是Sobel算子提取晶界的效果圖，從中可以看出圖像中有很多噪點，邊緣處有大量的毛刺。這是由于該算子只是檢測水平和豎直兩個方向的邊緣信息，計算得到的灰度值如果大于等于設(shè)定的閾值，就認為是邊緣點，這會導(dǎo)致將噪聲誤判為邊緣點。這種判別方式欠合理，適用于簡單紋理圖像處理，但對于復(fù)雜圖像的檢測，其效果并不理想。圖3c為傳統(tǒng)的分水嶺法的檢測結(jié)果圖，對比原圖可見丟失了許多重要輪廓信息。該方法是以數(shù)學(xué)形態(tài)學(xué)為基礎(chǔ)的一種分割方法，主要分2步實現(xiàn)：首先將圖像中每一個像素的灰度值按照從小到大排序，然后再從下往上實現(xiàn)淹沒過程并對每一局部受極小值影響的區(qū)域進行標注。分水嶺法通常處理的是梯度圖像，由于對噪聲很敏感，圖像中的噪聲點會影響其迭代標注，從而導(dǎo)致分割偏移，也容易出現(xiàn)過分割現(xiàn)象。圖3d是采用Roberts算子進行處理得到的晶界，提取的晶界不完整，該算子是根據(jù)局部差分算子來尋找邊界的一種簡單算子，局部差分算子是利用四鄰域?qū)蔷€方向相鄰的2個像素值之差，從實驗效果上可以看出，其對于確定邊緣位置較準，但是同樣容易受到噪聲的干擾。圖3e是哈爾小波處理后再閾值分割的效果圖，一些微弱邊界被去除，同時因為噪聲的干擾，圖中出現(xiàn)較多的黑點，直接影響到實驗結(jié)果。針對現(xiàn)有的方法容易受噪聲影響等問題，本研究提出的方法處理結(jié)果如圖3f所示，該方法很好地抑制了噪聲的干擾和偽邊緣的生成，較完整地提取出晶界。

3.2 不同算法的性能分析

本研究提出的方法雖然能夠提取更完整的晶界，但直接觀察還過于主觀?；诰Я€數(shù)和匹配率，分析各種算法用于晶界提取的性能(表1)。該金相圖像在實驗中不能直接用于比較，因為它不是一個標準的圖像。實驗中使用一款軟件用于提取晶界，雖然存在一些誤差，但在人工修飾的基礎(chǔ)上，對晶界的提取的誤差近乎為0，因此，使用軟件加手動標注晶粒的圖像作為數(shù)據(jù)的客觀標準具有一定的參考意義。表1中晶粒數(shù)量代表晶界提取效果的判斷，通過比較可以看出，Sobel法和哈爾小波處理方法不能完全提取出所有晶界，分水嶺法會出現(xiàn)過度提取現(xiàn)象。由不同算法提取的晶界與手動提取的晶界之間的錯誤提取率描述的是提取到錯誤晶界的比率，由表1可以看出：Sobel算子和Roberts算子因為有較多的偽邊界，導(dǎo)致錯誤率較高；分水嶺法的錯誤率也較高，這是由于過度分割，提取出了很多邊界，既包含正確的邊界，也包含錯誤的邊界。在計算錯誤率時，只考慮錯誤提取晶界的個數(shù)與手工提取的數(shù)量的比值。本研究提出的方法錯誤率遠低于其他算法，具有優(yōu)于其他方法的性能。

表1 不同方法用于晶界提取的性能比較Table 1 Performance comparison of different methods for grain boundary extraction

3.3 不同材料晶界提取的應(yīng)用實例及分析

為了進一步驗證所提出方法的可靠性，選取了幾種不同型號，且在不同放大倍數(shù)和腐蝕的情況下獲得的GH4169、GH738、GH706金相圖像為例進行測試，用于評判本研究所提出的晶界提取算法的準確性和穩(wěn)定性。圖4分別是GH4169、GH738、GH706材料晶界提取的效果圖，從圖中觀察比較可知，雖然材料、拍照環(huán)境和圖像尺寸不同，但本研究提出的方法提取的晶界仍然比較完整，可以達到較高的精度和較好的穩(wěn)定性，因此，對于不同類型的金相圖像的晶界提取具有一定的實用性，同時也為其他金屬材料的晶界提取提供了良好的思路。由于大多數(shù)金屬材料也同樣是由若干小晶粒組成，理論上使用本研究所提出的方法也具有同樣的效果，但具體問題還需要具體分析。

使用本研究方法提取出完整的邊界之后，根據(jù)圖中晶粒的輪廓，利用橢圓擬合法，以橢圓方程為模板盡可能地逼近晶界，使得對應(yīng)的橢圓方程盡量滿足邊緣上的數(shù)據(jù)，并求出橢圓方程的長短軸參數(shù)。由式(6)求出圓度和費雷特直徑。

P=∑r/N·Rd=A/l

(6)

式中：P為圓度，r是每個角的曲率半徑，N為晶粒角曲率半徑的個數(shù)，R表示晶粒內(nèi)最大內(nèi)接圓的半徑，A表示平均面積，l為晶粒平均截距。晶粒中棱角越多，則圓度越低，反之則越好。對同一種型號的金屬材料不同放大倍數(shù)下的金相圖像，計算出的橢圓長短軸、圓度以及長短軸之比，獲得的數(shù)據(jù)如表2所示。表2中列出了GH4169、GH738和GH706的參數(shù)，通過完整晶界的提取可以進一步計算出相關(guān)的參數(shù)，從而對材料進行分析。以GH4169為例：平均長軸計算結(jié)果為132.23 μm，與實際參數(shù)的誤差為±7.21 μm；平均短軸測試結(jié)果為97.11 μm，與實際值的誤差為±4.12 μm；圓度、長短軸之比、平均面積和平均晶粒直徑等參數(shù)誤差值都控制在較小的范圍內(nèi)。從這些數(shù)據(jù)可以看出，通過本研究中的方法提取出的晶粒計算出的參數(shù)與實際值非常接近，對衡量晶粒質(zhì)量有一定的實用價值。

表2 不同材料的特征參數(shù)Table 2 Characteristic parameters of different materials

4 結(jié)論

1)本研究提出了一種簡單高效、可靠的晶界提取方法。先將圖像進行預(yù)處理，然后使用Otsu法進行分割，最后再對提取出來的晶界作邊界平滑和細化去枝。該方法有效抑制了偽邊界的產(chǎn)生，同時得到較好的圖像效果。

2)對GH4149、GH706、GH738晶界提取結(jié)果進行對比試驗，驗證了該方法能有效提取晶界，并提高了晶界的完整性，且適用于各種金相圖像的檢測分析。