X射線圖像灰度值與透照厚度的定量關(guān)系

郭文明,陳宇亮

(北京郵電大學(xué) 軟件學(xué)院, 北京 100876)

?

X射線圖像灰度值與透照厚度的定量關(guān)系

郭文明,陳宇亮

(北京郵電大學(xué) 軟件學(xué)院, 北京 100876)

摘要：X射線數(shù)字成像檢測技術(shù)中，射線圖像對于沿透照方向的厚度變化，表現(xiàn)為灰度值的變化，如何由灰度值的變化得到被檢物體實際厚度的變化成為業(yè)內(nèi)研究的重點?；贐eer定理，結(jié)合影響X射線機射線強度的因素和射線強度與數(shù)字圖像的轉(zhuǎn)換關(guān)系，推導(dǎo)出灰度值與透照厚度之間的定量關(guān)系；通過搭建試驗平臺，得到系列數(shù)據(jù)，對理論計算值和試驗數(shù)據(jù)測量值進行了比對驗證。結(jié)果表明：理論灰度值與實際測試相差誤差不超過灰度滿量程的±10%。

關(guān)鍵詞：X射線;數(shù)字成像;缺陷;灰度值;定量分析

射線檢測作為無損檢測的主要方法之一，與其他方法相比具有直觀、準確等優(yōu)點。射線檢測中垂直于射線透照方向的缺陷尺寸可精確測量[1]，但平行于射線透照方向的缺陷尺寸無法直接測量得到，而未焊透、根部內(nèi)凹等缺陷深度的測量又是實際應(yīng)用中的常見問題。

對于未焊透、根部內(nèi)凹和根部咬邊等缺陷的射線檢測，為了確定底片黑度與厚度(缺陷深度)的關(guān)系，通常采用定制的通用槽型對比試塊或小徑管環(huán)焊縫專用對比試塊模擬焊縫及未焊透深度，得到帶缺陷和對比試塊的底片。但是底片上的缺陷較小，而一般黑度計孔至少在2.5 mm以上，如何通過測量得到精確的缺陷深度成為難點。而在先進的射線數(shù)字成像檢測中，由于使用了計算機和數(shù)字處理技術(shù)，使得通過計算機軟件實現(xiàn)缺陷尺寸的準確測量成為可能。

目前針對X射線數(shù)字圖像灰度值(以下簡稱灰度值)的研究，主要著重于分析灰度值與管電壓和被測物體厚度等因素的定性關(guān)系，灰度值與各因素簡單描述為簡單的正比或反比關(guān)系[2-4]，而對于圖像灰度值與被測物體厚度之間的關(guān)系，并沒有給出明確的計算公式或?qū)?yīng)關(guān)系。筆者通過分析Beer定律以及灰度值的相關(guān)理論，推導(dǎo)出精確的灰度值計算公式以及灰度值與厚度的對應(yīng)關(guān)系，并通過試驗測試公式的準確性。

1圖像灰度值與透照厚度的定量分析

1.1　射線強度與透照厚度的關(guān)系

對于某一物質(zhì)，射線穿過物質(zhì)的衰減系數(shù)與物質(zhì)密度、原子序數(shù)和射線能量有關(guān)，一般射線檢測的衰減理論均是基于Beer定律[2]發(fā)展起來的。在現(xiàn)有的研究中，給出了光子與材料的各種相互作用截面參數(shù)、衰減系數(shù)、吸收系數(shù)等國際標準的數(shù)據(jù)表[5-7]，這些數(shù)據(jù)表同樣驗證了在理想的單色源情況下，穿過物體的衰減規(guī)律符合Beer定律，如下式：

(1)

式中：μ為射線穿過單位厚度物質(zhì)時的衰減系數(shù);T為物體的厚度;I為射線透射后的強度；I0為入射射線強度。射線強度指在單位時間內(nèi)、通過垂直于射線傳播方向上單位面積的X射線光子的能量[8-9]。

對于X射線的照射來說，X射線管的輻射強度與管電壓的平方成正比，與管電流成正比。而根據(jù)平方反比定律，強度I0與距離F的平方成反比，所以即可得出下式：

(2)

式中：K1為常數(shù)系數(shù)，可以定義為灰度值常數(shù)；A為X射線機的管電流；U為X射線機的管電壓值；F為X射線機到被測物體表面的焦距距離。

由于X射線到達被檢物體表面的強度與管電壓、管電流、焦距有關(guān)，所以可將X射線強度看成是關(guān)于管電壓、管電流、焦距的函數(shù)：

(3)

1.2　線性衰減系數(shù)的確定

1.2.1衰減系數(shù)公式

射線管電壓的提高會造成射線有效能量的增加，而有效能量的增加會使射線波長變短，引起衰減系數(shù)下降[12-14]。根據(jù)光電效應(yīng)、原子物理學(xué)等理論，散射和吸收截面與衰減系數(shù)存在著換算關(guān)系，但在實際檢測中線性衰減系數(shù)隨著工件厚度變化而發(fā)生變化，而強濾波后的射線在X射線光譜中的分布應(yīng)為從最大強度附近到最短波長的范圍[15]。參考文獻[13] 給出了衰減系數(shù)和射線管電壓之間的關(guān)系：

(4)

式中:A,B,C,D,E,F為常數(shù)，隨透過材料的不同而不同，cm-1;V為射線的管電壓;K2為常數(shù)。

1.2.2鐵的衰減系數(shù)求解

求解A,B,C,D,E,F的六元一次方程,計算鐵的衰減系數(shù)。其中K2與被檢測材料的厚度有關(guān),其隨著材料厚度的增加而逐漸減小。在進行衰減系數(shù)測定時，對于連續(xù)譜的X射線來說會發(fā)生硬化現(xiàn)象，需要對衰減系數(shù)進行修正。對于給定厚度為x的衰減系數(shù)，μavex修正為理想單色源的衰減系數(shù)μave0時，需要加上修正項Δμ0;Δμ0的計算需結(jié)合Simpson公式，導(dǎo)出X射線修正模型的數(shù)值解法及修正方法，對修正后的衰減系數(shù)再作卷積反投影重構(gòu)，即可有效消除能譜硬化造成的影響[4-5]。由于Δμ0的計算公式冗長，文章不過多展開，具體方法參照文獻[5]。綜合式(4)和能量譜硬化修正模型，得出鐵在不同的射線管電壓下，衰減系數(shù)的變化情況，具體數(shù)值如表1所示。

表1　鐵在不同管電壓下的衰減系數(shù)

從表1可看出，入射能量隨著管電壓增大而增大，對應(yīng)的平均衰減系數(shù)減小。通過式(4)可知，被測材料的厚度增大，同樣會使對應(yīng)的平均衰減系數(shù)減小。這與文獻[5]的研究結(jié)論相一致，說明了表1的正確性。

1.3　灰度值與射線強度的關(guān)系

射線數(shù)字圖像是以灰度值的形式顯示的，即需要先把X射線強度轉(zhuǎn)化為數(shù)字信號(圖像的灰度值)。首先，將入射X射線光子轉(zhuǎn)換為電荷[10-11]，照射在像元上的X射線光子數(shù)量與像元轉(zhuǎn)換的電荷數(shù)近似存在線性關(guān)系，其具體公式如下：

(5)

式中：M為光通量；R為光電變換系數(shù)。

模數(shù)轉(zhuǎn)換器(A/D)將得到的電壓轉(zhuǎn)化為數(shù)字信號，形成數(shù)字圖像[15]。模擬數(shù)字轉(zhuǎn)換器的分辨率Q等于LSB電壓(輸出離散信號產(chǎn)生一個變化所需的最小輸入電壓的差值)：

(6)

式中：EFSR為總的電壓測量范圍；2D為電壓間隔數(shù)，D為ADC模塊的精度的位數(shù)，用比特作為單位。

在上述過程中，可直接用已經(jīng)量化的電壓值(即比特信息)作為圖像的灰度值，結(jié)合式(5)、(6)可推導(dǎo)出圖像灰度值的變化規(guī)律：

(7)

式中：f(R0,E0,D0)為關(guān)于R0,E0,D0參數(shù)的函數(shù)；R0為光電變換系數(shù)；E0為模數(shù)轉(zhuǎn)換器的轉(zhuǎn)換系數(shù)；D0為灰度值位數(shù)；G為圖像灰度值。

結(jié)合式(3)、(4)，可得到圖像灰度值的計算公式：

(8)

對某一檢測物體來說，由于函數(shù)f(R0,E0,D0)的值與檢測系統(tǒng)有關(guān)，在固定檢測系統(tǒng)后，函數(shù)f(R0,E0,D0)應(yīng)該是一個常數(shù)，簡化為常數(shù)系數(shù)K。函數(shù)I0(A,U,F)雖然是關(guān)于射線強度的函數(shù)，但在X射線機發(fā)射器固定的情況下，可將函數(shù)I0(A,U,F)的具體形式換算成式(2)，得到試驗驗證公式(9)。

(9)

從式(9)可看出:數(shù)字圖像灰度值的大小受到材料衰減系數(shù)、管電壓、管電流、曝光時間和焦距的影響，同時與檢測系統(tǒng)的硬件特性有很大關(guān)系。

2計算結(jié)果與試驗驗證

2.1　公式系數(shù)K求解

為了求解公式系數(shù)K，采用CMOET 225 kV X射線機和以色列vidisco非晶硅數(shù)字探測器(像素127 μm)搭建射線檢測系統(tǒng)平臺，選取不同厚度的鋼板(7,10,12 mm)作為試驗用試塊(試塊尺寸為210 mm×150 mm)，通過改變管電壓，管電流和焦距來獲得不同數(shù)字圖像。試驗主要影響因素包括管電壓的不穩(wěn)定和硬化現(xiàn)象。表2給出了鋼板厚度為7 mm時，采用不同的管電壓(92,95,97 kV)、管電流(5,4,3 mA)、焦距(700,800,900 mm)時得到的圖像灰度值。由于衰減系數(shù)μ隨著管電壓而改變，故可通過查閱表1來確定衰減系數(shù)μ。將不同厚度下的管電壓、管電流、焦距和圖像灰度值的數(shù)值代入式(9)，可求出K值。在X射線系統(tǒng)固定的情況下，公式對于任何情況下的X射線都滿足，那么對于上述的數(shù)據(jù)來說，K值的波動不會太大。綜合另外兩組厚度為10，12 mm的數(shù)據(jù)，求出K值的平均值為428.030 1。從計算結(jié)果來看，K值誤差波動范圍在-3.4%～3.0%之間，大體趨近一個常數(shù)。

表2　不同透照條件下公式系數(shù)K的計算結(jié)果

2.2　試驗驗證

試驗平臺同上，采用管電流5 mA、焦距800 mm,改變管電壓(95,100,105 kV)，對6～12 mm厚的鋼制試塊進行灰度值的理論計算與試驗測試結(jié)果的比較。試驗主要影響因素同上。結(jié)果見表3,變化趨勢如圖1所示。

表3　不同管電壓下灰度值的試驗與計算結(jié)果比較

圖1　不同管電壓下灰度值和透照厚度的關(guān)系曲線

從表3可看出，在運用式(9)對圖像的灰度值進行計算和試驗驗證的過程中，在一定的透照條件下，求解得出的理論值與試驗值的誤差在±10%間。而從圖1可看出，在管電壓固定的情況下，連續(xù)計算一組數(shù)據(jù)的灰度值與試驗值的曲線大致吻合。產(chǎn)生這種現(xiàn)象的主要原因有：

① 數(shù)字圖像灰度值與數(shù)字射線檢測系統(tǒng)直接相關(guān)，不同的數(shù)字探測器其響應(yīng)不同，即便是同一數(shù)字探測器隨著時間推移也會發(fā)生相應(yīng)退化，而造成數(shù)字圖像灰度值變化。② X射線機是連續(xù)的射線源，在通過修正模型對射線的硬化現(xiàn)象進行校正的過程中，會出現(xiàn)計算偏差。③ 隨著鋼板厚度的增加，材料的衰減系數(shù)也會發(fā)生變化，衰減系數(shù)很大程度影響了灰度值的計算。④ 式(9)中K值是通過大量試驗數(shù)據(jù)求出的平均值，對于個別數(shù)據(jù)的驗證仍會出現(xiàn)較大的誤差，但是對于大量數(shù)據(jù)的驗證，從圖1的曲線中可以看出，試驗值和計算值比較吻合。

3結(jié)論

基于射線衰減理論，分析不同透照條件下引起射線圖像灰度值變化的影響因素，圖像灰度值與X射線機的能量、材料的厚度和衰減系數(shù)的定量關(guān)系，推導(dǎo)出不同鋼板厚度和圖像灰度值之間的定量關(guān)系，并搭建試驗平臺予以驗證。通過對試驗得到的多組數(shù)據(jù)的分析計算表明：對于一定厚度的工件，其灰度的理論值與試驗值存在一定的偏差，但在實際的射線檢測中，可基于公式初步估算厚度與灰度的關(guān)系，為后續(xù)開展基于灰度值精確測量缺陷深度奠定基礎(chǔ)。

參考文獻：

[1]SIM L M, TAN T C, MON A A, et al. X-ray absorption-based technique to measure the thickness of multi-layered structures[J]. Ndt & E International, 2009, 42(4):291-296.

[2]PENG Guang-han, YANG Xue-heng, HAN Zhong, et al. Hardening correction model of energy spectrum for continuous spectrum X-ray ICT [J].Spectroscopy and Spectral Analysis,2005,25(11):1880-1883.

[3]XU Guang-duo, WANG Li-qiang, TONG Jian-min. Research on calibration model of X-ray thickness gauge[J].Atomic Energy Science and Technology,2014,5:925-929.

[4]CHEN Ping, HAN Yan, PAN Jin-xiao. Research on X-ray multispectrum imaging based on variable energy [J].Spectroscopy and Spectral Analysis,2013,32(5):1383-1387.

[5]PENG G H, CAI X H, YANG X H. Numerical analysis of hardening correction model of energy spectrum for X-ray TICT[J]. Nondestructive Testing, 2007.

[6]CULLEN D E. A simple model model of photon transport[J].Nucl InstrumMetthodsPhys Res,1995,B 101(1):499-510.

[7]KOLKOORI S, WROBEL N, ZSCHERPEL U, et al. A new X-ray backscatter imaging technique for non-destructive testing of aerospace materials[J]. Ndt & E International, 2014,70:41-52.

[8]SELTZER S M. Calculation of photon mass enery-transfer and mass enery-absorbtioncoefficient[J].Radiat Res,1993,136(2):147-170.

[9]REZA H P R, REZAIE A H, SADEGHI S H H, et al. A density-based fuzzy clustering technique for non-destructive detection of defects in materials[J]. Ndt& E International,2007,40(4):337-346.

[10]CASALTA S, DAQUINO G G, METTEN L, et al. Digital image analysis of X-ray and neutron radiography for the inspection and the monitoring of nuclear materials[J]. Ndt & E International,2003,36(5):349-355.

[11]趙肖東,蔡建剛,陳明,等. 航空零部件X射線數(shù)字成像檢測系統(tǒng)[J]. 無損檢測, 2013, 35(4):57-59.

[12]王瑋,李公平,潘小東,等. 工業(yè)CT錐束X射線能譜及強度的分布模擬[J].無損檢測,2014,36(3):9-13.

[13]馮長春,周煒. X射線衰減系數(shù)與射線管電壓的關(guān)系初探[J].無損探傷,1997,19(1):11-13.

[14]高闊,于虹,郭鐵橋. 射線數(shù)字成像(DR)技術(shù)在電力工業(yè)檢測中的應(yīng)用[J]. 無損檢測, 2013, 35(11):76-78.

[15]連廣，王庭魁.材中X射線衰減規(guī)律的驗證及吸收系數(shù)測定[J].林業(yè)科技，1990(6):28-30.

Quantitative Relationship between Gray Value of X-ray Image and Radiographic Thickness

GUO Wen-ming, CHEN Yu-liang

(School of Software Engineering,Beijing University of Posts and Telecommunications, Beijing 100876,China)

Abstract：In the technology of X-ray digital imaging detection, the X-ray imaging shows changes in gray value for thickness variations along the transillumination direction. The way how to get the object's actual thickness changes by grey value change becomes the focus of the research. In this paper, based on Beer theorem and combining with the factors influencing the intensity of X-ray machine and ray intensity and digital image transformation relations, the quantitative relationship between the grey value and thickness is deduced. By setting up the experimental platform to get series of experimental data, the theoretical calculation value and experimentally measured values are compared to verify. The results show that the difference between theoretical gray value and actual grey value is not more than gray full range of ±10%.

Key words:X-ray; Digital imaging; Detection; Gray value; Quantitative analysis

中圖分類號：TG115.28

文獻標志碼：A

文章編號：1000-6656(2016)02-0014-04

DOI:10.11973/wsjc201602004

作者簡介：郭文明(1967-),男,副教授,主要從事圖像處理的教學(xué)與科研工作。

基金項目：國家重大科學(xué)儀器設(shè)備開發(fā)專項資助項目(2013YQ240803)

收稿日期：2015-06-23