劉忠奎，葛俊鋒，張羽鵬，張　興

示溫漆溫度自動(dòng)判讀技術(shù)研究

劉忠奎1，葛俊鋒2，張羽鵬1，張興1

（1.（中航工業(yè)）沈陽發(fā)動(dòng)機(jī)設(shè)計(jì)研究所，遼寧沈陽110015；2.華中科技大學(xué)自動(dòng)化學(xué)院，湖北武漢430074）

在總結(jié)國(guó)內(nèi)外示溫漆自動(dòng)判讀算法研究現(xiàn)狀的基礎(chǔ)上，詳細(xì)討論現(xiàn)有算法的優(yōu)缺點(diǎn)，進(jìn)一步研究示溫漆的自動(dòng)判讀技術(shù)。通過引用紅外特征的方法，并采用圖像融合將其與示溫漆顏色特征組合成新的圖像數(shù)據(jù)，然后提取出新的圖像數(shù)據(jù)的統(tǒng)計(jì)特征，采用機(jī)器學(xué)習(xí)算法自動(dòng)建立等溫線和像素點(diǎn)溫度識(shí)別模型，兩種模型相結(jié)合用于示溫漆溫度自動(dòng)判讀，達(dá)到提高示溫漆測(cè)溫準(zhǔn)確度的目的。

示溫漆；自動(dòng)判讀算法；顏色溫度特性曲線；顏色空間

doi：10.11857/j.issn.1674-5124.2015.09.005

0　引言

示溫漆是一種溫度敏感涂料，其涂敷在物體表面并隨物體溫度變化而改變顏色，實(shí)現(xiàn)物體溫度的測(cè)量。示溫漆作為一種接觸的、無損害的測(cè)溫手段，具有方便、成本低、不受測(cè)試位置限制、對(duì)試件無影響等優(yōu)點(diǎn)［1］，廣泛應(yīng)用在航空領(lǐng)域的溫度測(cè)量中。

示溫漆溫度判讀是示溫漆測(cè)溫的關(guān)鍵，目前主要采用人工判讀的方法［2］，但該方法判讀精度較低。因此，近年來基于圖像處理的示溫漆溫度自動(dòng)識(shí)別方法得到廣泛關(guān)注。

國(guó)外，最早提出的示溫漆自動(dòng)判讀方法是基于RGB空間顏色溫度特性曲線的溫度判別方法［3］，針對(duì)該方法的缺陷，又提出了消除特性曲線中歧義部分的馬爾科夫隨機(jī)場(chǎng)方法［4-5］。2008年Lempereur C等［6］提出了目前最好的示溫漆自動(dòng)判讀方法，即在Lab空間進(jìn)行圖像分割，找出與變色點(diǎn)對(duì)應(yīng)的區(qū)域，在HSI空間進(jìn)行像素顏色匹配，用最近鄰法確定每個(gè)區(qū)域中像素點(diǎn)的溫度。

目前，國(guó)內(nèi)示溫漆自動(dòng)判讀技術(shù)研究主要有在LUV空間建立顏色溫度曲線，利用歐氏距離最近判讀示溫漆像素點(diǎn)溫度的方法［7］，基于LUV空間顏色溫度曲線的點(diǎn)和區(qū)域溫度識(shí)別方法［8］，采用HSI空間顏色溫度曲線的K均值聚類、中位切割等示溫漆圖像分割及量化方法［9-10］，采用Lab空間的顏色溫度曲線判斷示溫漆單像素溫度的方法［11］，基于圖像處理及圖像融合算法的示溫漆溫度識(shí)別方法［12-13］以及基于等溫線溫度識(shí)別的示溫漆溫度自動(dòng)判讀方法［14］。

1　示溫漆自動(dòng)判讀的基本原理

示溫漆自動(dòng)判讀分兩步，第一步為學(xué)習(xí)，即由標(biāo)準(zhǔn)樣板構(gòu)建示溫漆顏色溫度模型；第二步為識(shí)別，即將試件上示溫漆的顏色與模型匹配，得到相應(yīng)的溫度。

示溫漆自動(dòng)判讀的主要依據(jù)是示溫漆的標(biāo)準(zhǔn)樣板。目前，國(guó)內(nèi)外的標(biāo)準(zhǔn)樣板基本一致，通常有兩種形狀，一種是長(zhǎng)方形，可用于多變色示溫漆和單變色示溫漆的標(biāo)定，如圖1（a）所示；另一種是蝴蝶形狀，其阻抗呈梯形分布，從而在蝶形片上呈現(xiàn)示溫漆所有變化顏色，用于多變色示溫漆的標(biāo)定，如圖1（b）所示［6］。

圖1　示溫漆的標(biāo)準(zhǔn)樣板

可以看出：示溫漆每隔一定溫度就會(huì)有一個(gè)顏色變化明顯的變色點(diǎn)。在變色點(diǎn)兩邊，顏色變化具有階躍性，在兩個(gè)變色點(diǎn)之間，顏色變化具有漸變性，有些尚可區(qū)分，有些不易區(qū)分。

根據(jù)上述變色特點(diǎn)，示溫漆自動(dòng)判讀算法在構(gòu)建顏色溫度模型時(shí)通常有兩種方式：1）直接對(duì)示溫漆顏色和溫度關(guān)系進(jìn)行建模，通過擬合插值等方法獲得顏色溫度特性曲線；2）對(duì)示溫漆變色點(diǎn)兩側(cè)的顏色和溫度關(guān)系進(jìn)行建模，通過插值實(shí)現(xiàn)等溫線的自動(dòng)判讀。

2　現(xiàn)有方法的分析與比較

根據(jù)示溫漆顏色溫度模型和顏色匹配方法的不同，現(xiàn)有示溫漆自動(dòng)判讀方法可以歸納為3類，即基于顏色溫度特性曲線的方法，基于等溫線溫度識(shí)別的方法以及基于變色點(diǎn)區(qū)域分割的方法。

2.1基于顏色溫度特性曲線的方法

該方法關(guān)鍵是用樣條曲線插值法獲得顏色溫度特性曲線，然后根據(jù)離曲線距離最近的原則對(duì)像素點(diǎn)顏色對(duì)應(yīng)的溫度進(jìn)行判讀?？梢圆捎肦GB、HSI、Lab、YUV等多種顏色空間，實(shí)現(xiàn)步驟［7］：

1）根據(jù)示溫漆標(biāo)準(zhǔn)樣板建立色溫?cái)?shù)據(jù)庫，例如RGB值為（30，45，36）對(duì)應(yīng)溫度410℃。

2）建立RGB顏色溫度特征曲線，因RGB像素點(diǎn)個(gè)數(shù)為255×255×255，人工方式不可能測(cè)出所有RGB對(duì)應(yīng)溫度，所以需借助曲線擬合來識(shí)別未知RGB像素點(diǎn)的溫度。例如某型示溫漆的顏色溫度特征曲線如圖2所示。

圖2　示溫漆的顏色溫度特性曲線

3）進(jìn)行曲線擬合，用三次樣條插值法識(shí)別溫度。如圖2所示，為求未知點(diǎn)A的溫度T，根據(jù)式（1）對(duì)溫度特征曲線中所有點(diǎn)全局掃描，min（Δ）對(duì)應(yīng)點(diǎn)BI，點(diǎn)BI的溫度近似為A像素點(diǎn)的溫度。

此算法在不同的顏色空間都能實(shí)現(xiàn)，如RGB空間、LUV空間等，但不同的顏色空間對(duì)于該算法的性能有一定影響。表現(xiàn)為由于環(huán)境光改變引起圖像顏色偏差時(shí)，在LUV空間直接判讀仍能保持較高的識(shí)別準(zhǔn)確率，但在RGB空間的識(shí)別率明顯下降。

基于顏色溫度特性曲線進(jìn)行溫度識(shí)別的方法簡(jiǎn)單易行，但有以下缺陷：

1）示溫漆在應(yīng)用中受到多種不定因素影響，很難有理想的成像圖像。該算法過于依賴RGB像素點(diǎn)，導(dǎo)致溫度識(shí)別可靠性差，準(zhǔn)確度低。

2）用一個(gè)RGB值表示一個(gè)標(biāo)準(zhǔn)樣板的溫度，存在一定不足。如果標(biāo)準(zhǔn)樣板的顏色差別比較小，該方法就容易出現(xiàn)誤識(shí)別。

另外，基于顏色溫度曲線的區(qū)域溫度識(shí)別方法，即先對(duì)示溫漆圖像進(jìn)行分割［8］，將顏色相近的像素分割在同一區(qū)域，這樣不同區(qū)域的顏色差別非常明顯，然后對(duì)各個(gè)區(qū)域的溫度進(jìn)行識(shí)別。但該方法識(shí)別的是分割區(qū)域的平均顏色對(duì)應(yīng)的溫度，溫度的分辨率依賴于示溫漆能夠明顯變化的顏色，通常由于示溫漆的變色點(diǎn)有限，該方法的溫度分辨率和溫度識(shí)別的準(zhǔn)確度都不高。

2.2基于等溫線溫度識(shí)別的方法

該方法首先識(shí)別各條等溫線的位置和數(shù)量，然后用K-Means聚類識(shí)別各等溫線附近區(qū)域顏色，最后與標(biāo)準(zhǔn)樣板對(duì)比，識(shí)別等溫線的溫度，具體實(shí)現(xiàn)步驟［14］：

1）提取等溫線。示溫漆圖像依次經(jīng)濾波、邊緣增強(qiáng)、改進(jìn)的Canny算子邊緣提取和雜點(diǎn)處理等步驟，準(zhǔn)確識(shí)別等溫線。

2）等溫線編號(hào)。通過全局掃描，把大于設(shè)定長(zhǎng)度的等溫線依次編號(hào)為line［i］，同時(shí)計(jì)算出各相鄰等溫線間最近距離distance［i］。等溫線長(zhǎng)度閥值選100左右為宜，這樣可減小等溫線數(shù)量，提高運(yùn)算速度。

3）跟蹤等溫線。逐一跟蹤等溫線line［i］，把其附近distance［i］范圍區(qū)域內(nèi)灰度素點(diǎn)值存入矩陣A、B。選擇灰度值而非RGB值的原因是RGB像素矩陣大計(jì)算負(fù)擔(dān)重，且因RGB與灰度值是線性變換的，具有一致性。

4）聚類法識(shí)別A、B聚類中心。用K-Means聚類法計(jì)算矩陣A、B各自的聚類中心a、b。a、b代表了等溫線附近區(qū)域的像素分布特征。

5）建立區(qū)域像素特征庫。獲取各示溫漆標(biāo)準(zhǔn)樣板組成該示溫漆標(biāo)準(zhǔn)樣板庫，用四叉樹法量化樣板，抽象出溫度數(shù)據(jù)庫。對(duì)于五變色示溫漆，其數(shù)據(jù)庫最小只有5組數(shù)據(jù)，相對(duì)基于RGB顏色溫度特征曲線的示溫漆溫度識(shí)別算法，數(shù)據(jù)庫數(shù)據(jù)量明顯減少。

6）溫度識(shí)別。根據(jù)式（2），利用插值原理列掃描示溫漆數(shù)據(jù)庫，計(jì)算點(diǎn)（a，b）與數(shù)據(jù)庫標(biāo)準(zhǔn)樣板矩陣的距離。若min（Δ）對(duì)應(yīng)點(diǎn)（aI，bI），則該點(diǎn)對(duì)應(yīng)的溫度即為該等溫線的溫度。

國(guó)內(nèi)大都采用顏色進(jìn)行示溫漆標(biāo)定和判讀，這種標(biāo)定和判讀誤差超過±50℃以上。而對(duì)等溫線進(jìn)行標(biāo)定和判讀，判讀準(zhǔn)確度可達(dá)±10℃。

基于等溫線溫度識(shí)別方法的優(yōu)點(diǎn)為：

1）對(duì)等溫線進(jìn)行識(shí)別，而不是像素點(diǎn)，識(shí)別溫度的特征選擇不同，可以提高溫度識(shí)別的準(zhǔn)確度。

2）采用統(tǒng)計(jì)特征，而非單個(gè)像素的顏色特征對(duì)溫度進(jìn)行識(shí)別，可以減小光照對(duì)識(shí)別結(jié)果的影響。

該算法的不足為：

1）等溫線檢測(cè)采用Canny算子，雜點(diǎn)多，需要一定經(jīng)驗(yàn)值，以及后續(xù)處理。

2）特征采用單個(gè)平均灰度值過于簡(jiǎn)單，由于等溫線個(gè)數(shù)較少，不利于后續(xù)的插值計(jì)算。

3）對(duì)等溫線之間的區(qū)域沒有進(jìn)行進(jìn)一步的溫度識(shí)別，溫度分辨率不高。

2.3基于變色點(diǎn)區(qū)域分割的方法

該方法首先在Lab空間中利用種子點(diǎn)和區(qū)域增長(zhǎng)算法對(duì)示溫漆圖像進(jìn)行分割，通過經(jīng)驗(yàn)閾值，獲得與各段變色點(diǎn)間顏色相對(duì)應(yīng)的區(qū)域，同時(shí)排除與示溫漆顏色不能匹配的區(qū)域。然后在HSI空間中將各個(gè)區(qū)域中的像素顏色特征與標(biāo)準(zhǔn)樣板的顏色特征對(duì)比，對(duì)每個(gè)像素對(duì)應(yīng)的溫度進(jìn)行判讀，具體實(shí)現(xiàn)步驟［6］：

1）顏色空間變換。將RGB空間變換到Lab空間。

2）圖像分割。利用種子點(diǎn)和區(qū)域增長(zhǎng)算法對(duì)示溫漆圖像進(jìn)行分割，通過經(jīng)驗(yàn)閾值，獲得與各段變色點(diǎn)間顏色相近的區(qū)域。

3）變色點(diǎn)區(qū)域匹配。匹配分割區(qū)域顏色與各變色點(diǎn)的顏色信息，將各個(gè)區(qū)域與變色點(diǎn)對(duì)應(yīng)，同時(shí)排除與示溫漆顏色不能匹配的區(qū)域。匹配分割區(qū)域后排除了不能識(shí)別的區(qū)域，同時(shí)將各個(gè)區(qū)域像素溫度的識(shí)別限定在不同的溫度段，有利于提高后續(xù)溫度識(shí)別的準(zhǔn)確度。

4）顏色空間變換。將需要識(shí)別區(qū)域的像素從RGB空間變換到HSI空間，以便將亮度信息剔除，減少光照對(duì)溫度識(shí)別的影響。

5）像素點(diǎn)溫度識(shí)別。將各個(gè)區(qū)域中的像素點(diǎn)顏色特征與標(biāo)準(zhǔn)樣板的顏色特征進(jìn)行對(duì)比，選擇匹配度最高的樣板溫度作為對(duì)應(yīng)像素點(diǎn)的溫度。

基于變色點(diǎn)區(qū)域分割的方法優(yōu)點(diǎn)為：

1）利用不同的顏色空間特性進(jìn)行示溫漆圖像的分割（Lab）和溫度識(shí)別（HSI），提高算法的性能。

2）采用由粗到細(xì)的溫度識(shí)別方法。首先進(jìn)行區(qū)域分割，匹配區(qū)域溫度范圍并排除不能識(shí)別區(qū)域，其次在對(duì)應(yīng)溫度段內(nèi)進(jìn)一步識(shí)別區(qū)域內(nèi)的像素溫度，提高溫度識(shí)別的準(zhǔn)確度。

該算法的不足為：

1）識(shí)別像素點(diǎn)溫度時(shí)，采用的像素點(diǎn)鄰域特征的區(qū)分性不足。當(dāng)標(biāo)準(zhǔn)樣板顏色變化不明顯時(shí)，使用HSI空間中的色調(diào)（H）和飽和度（S）特征，不能有效識(shí)別像素點(diǎn)溫度。

2）在圖像分割中，需要用經(jīng)驗(yàn)閾值使得分割區(qū)域與變色點(diǎn)區(qū)域重合，這需要頻繁的人機(jī)交互，有時(shí)并不能保證獲得滿意的分割。

3　新方法研究

示溫漆自動(dòng)判讀技術(shù)的核心在于構(gòu)建顏色與溫度之間的關(guān)系模型，該模型需要對(duì)光照條件變化魯棒，難點(diǎn)在于如何解決示溫漆漸變色區(qū)域的溫度識(shí)別。為了進(jìn)一步提高示溫漆自動(dòng)判讀的準(zhǔn)確度，借鑒現(xiàn)有示溫漆自動(dòng)判讀方法中的優(yōu)點(diǎn)，新方法的研究重點(diǎn)如下。

1）從示溫漆彩色圖像處理，擴(kuò)展到示溫漆多光譜圖像處理。除了獲取顏色特征外，同時(shí)獲得紅外特征，紅外特征受可見光影響小，可以提高判讀的準(zhǔn)確性。

2）示溫漆等溫線兩側(cè)的顏色變化非常明顯，而等溫線之間顏色變化不明顯，因此對(duì)示溫漆圖像采用從粗到細(xì)的分析策略獲得溫度信息，先檢測(cè)等溫線，將示溫漆圖像劃分為顏色相近的不同區(qū)域，再對(duì)區(qū)域中的像素進(jìn)行溫度識(shí)別。

3）不能再用單個(gè)像素的顏色值作為示溫漆的特征表示方式，因其很容易受到光照條件干擾，應(yīng)該采用區(qū)域或者像素鄰域的統(tǒng)計(jì)特征，如直方圖、矩陣征等，通過歸一化減少光照影響，并獲得一個(gè)高維特征，具有更好的區(qū)分性。

4）現(xiàn)有方法通常采用擬合和插值來獲得顏色與溫度模型（即顏色溫度特性曲線），采用的特征易受干擾，區(qū)分性不強(qiáng)，該模型實(shí)際應(yīng)用效果不好。在采用高維統(tǒng)計(jì)特征的情況下，使用機(jī)器學(xué)習(xí)算法自動(dòng)構(gòu)建模型，將標(biāo)定樣板區(qū)域顏色直方圖統(tǒng)計(jì)特征作為輸入，通過相應(yīng)訓(xùn)練算法，獲得等溫線溫度識(shí)別模型，將標(biāo)準(zhǔn)樣板像素點(diǎn)與其領(lǐng)域顏色的組合特征作為輸入，通過相應(yīng)訓(xùn)練算法，獲得像素點(diǎn)溫度識(shí)別模型。

5）現(xiàn)有方法都采用歐氏距離來度量顏色與溫度的關(guān)系，其合理性存疑。在高維特征下，采用能夠區(qū)分漸變顏色特征的度量方法，找出特征與溫度之間更準(zhǔn)確的關(guān)系。

4　結(jié)束語

與其他測(cè)溫方法相比，示溫漆在航空測(cè)試領(lǐng)域具有不可替代的優(yōu)勢(shì)，隨著航空發(fā)動(dòng)機(jī)試驗(yàn)研究的不斷深人，對(duì)示溫漆測(cè)溫的精度要求越來越高，未來示溫漆測(cè)試技術(shù)的發(fā)展一方面可以從材料和工藝入手，研制變色靈敏、具有多變色點(diǎn)、變色間隔小、色差變化大的新型示溫漆，另一方面可以利用圖像采集器件，圖像處理，機(jī)器學(xué)習(xí)等領(lǐng)域的新方法、新技術(shù)研究新的示溫漆自動(dòng)判讀方法，實(shí)現(xiàn)更高的示溫漆判讀準(zhǔn)確度。

［1］馬春武.示溫漆溫度自動(dòng)判讀與數(shù)字圖像處理研究［D］.南京：南京航空航天大學(xué)，2008.

［2］徐鳳花.示溫漆技術(shù)在航空發(fā)動(dòng)機(jī)高溫部件表面溫度測(cè)試上的應(yīng)用研究［D］.成都：電子科技大學(xué)，2009.

［3］Griffin A，Dittler J，Windeatt T，et al.Techniques for the interpretation of thermal paint coated samples［C］∥IEEE.Proceedings of Int Conf on Pattern Recognition 1996.Vienna：IEEE，1996：959-961.

［4］Griffin A，Windeatt T，Stoddart A J，et al.A Markov random fields approach for interpreting thermochromic paint［C］∥IEEE Proceedings of Int Conf on Pattern Recognition 1998.Brisbane：IEEE，1998：1741-1743.

［5］Smith M.Interpretation of thermal paints：Europe，EP0947813［P］.1999-03-19.

［6］Lempereur C，Andraland R，Prudhomme J Y.Surface temperature measurement on engine components by means of irreversible thermal coatings［J］.Measurement Science and Technology，2008（19）：105501.

［7］張志龍，曹承倜.示溫漆顏色溫度特性分析與溫度識(shí)別系統(tǒng)［J］.計(jì)算機(jī)自動(dòng)測(cè)量與控制，2001（3）：20-21.

［8］王美玲.基于彩色圖像處理的示溫漆溫度識(shí)別系統(tǒng)［D］.南京：南京航空航天大學(xué)，2009.

［9］林茂松，陳念年.示溫漆色彩量化及其圖像分割算法研究［J］.西南科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2006（2）：48-53.

［10］蔡茂蓉，龔巍.基于小波和模糊聚類的示溫漆彩色圖像分割算法［J］.微電子學(xué)與計(jì)算機(jī)，2010（5）：170-173.

［11］張興，薛秀生，陳斌.示溫漆在發(fā)動(dòng)機(jī)測(cè)試中的應(yīng)用與研究［J］.測(cè)控技術(shù)，2008（1）：21-23.

［12］于坤林.基于數(shù)字圖像處理的溫度場(chǎng)測(cè)量方法的研究［J］.長(zhǎng)沙航空職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2011（2）：31-33.

［13］于坤林，謝志宇.基于多傳感器圖像融合的溫度場(chǎng)測(cè)試系統(tǒng)［J］.現(xiàn)代電子技術(shù)，2013（20）：97-99.

［14］王榮華，杜平安，黃明鏡.基于等溫線溫度識(shí)別的示溫漆溫度自動(dòng)識(shí)別算法［J］.電子測(cè)量與儀器學(xué)報(bào)，2010（6）：542-547.

Study on the automatic tem perature recognition technology of thermal paint

LIU Zhongkui1，GE Junfeng2，ZHANG Yupeng1，ZHANG Xing1
（1.Shenyang Engine Design and Research Institute，Shenyang 110015，China；2.School of Automation，Huazhong University of Science and Technology，Wuhan 430074，China）

The advantages and disadvantages of existing algorithms for automatic temperature reading were discussed and the automatic temperature reorganization techniques of thermal paint were researched according to the summarization of current research at home and abroad on automatic temperature recognition algorithm of thermal paint.Infrared feature was introduced in combination with color feature to form new image data by image fusion.Afterwards，statistical features were extracted from these new data，and two temperature identification models for isotherm line and pixel were automatically created by machine-learning algorithm.The two models were used to automatically identify the thermal paint to improve its temperature-measuring accuracy.

thermal paint；automatic recognition algorithm；color-temperature curve；color space

A

1674-5124（2015）09-0020-04

2015-02-03；

2015-03-30

劉忠奎（1968-），男，遼寧沈陽市人，高級(jí)工程師，研究方向?yàn)楹娇瞻l(fā)動(dòng)機(jī)測(cè)試技術(shù)研究與應(yīng)用。