陳樹越 路宏年
摘要:目的 討論數(shù)字射線成像技術(shù)在國內(nèi)外的發(fā)展. 方法 根據(jù)目前國內(nèi)外數(shù)字射線成像技術(shù)的發(fā)展,著重討論X射線轉(zhuǎn)換屏和 CCD 相機(jī)在成像系統(tǒng)中的地位,以及數(shù)字圖像處理在該領(lǐng)域的應(yīng)用. 結(jié)果X射線轉(zhuǎn)換屏在成像質(zhì)量上起著非常重要的作用;采用增加 CCD 內(nèi)的累積電荷量和多幀疊加可以提高信噪比;數(shù)字圖像處理技術(shù)在改善X射線圖像上是*的. 結(jié)論 數(shù)字射線成像技術(shù)在無損檢測和評估方面將有很好的發(fā)展前景.
關(guān)鍵詞:數(shù)字式射線照相;無損檢測;圖像增強(qiáng);圖像恢復(fù)
中圖分類號(hào):TG115.281
文獻(xiàn)標(biāo)識(shí)碼:A
Research and Development of Digital X-ray Radiography in Nondestructive Testing
CHEN Shu-yue,LU Hong-nian
(Dept. of Electronic Engineering and Information Science, North China Institute of Technology,Taiyuan 030051,China)
Abstract:Aim To discuss the development of digital radiography in internationality. Methods According to the development of digital radiography,to analyse chiefly the function of X-ray-to-light converters and CCD cameras in the imaging system,and the application of digital image processing in the field. Results X-ray-to-light converter plays an important role on image quality. To increase the accumulated charges in CCD and multiframe superposition can improve signal-noise ratio. Digital image processing is necessary for improvement image quaeity of X-ray radiograph. Conclusion Digital Radiography will develope widely in non-destructive testing and evaluation.
Key words:digital radiography;non-destructive testing;image enhancement;image restoration
0 引言
自從1895年倫琴發(fā)現(xiàn)X射線以來,X射線就被廣泛應(yīng)用于醫(yī)學(xué),工業(yè)技術(shù)和科學(xué)研究等領(lǐng)域,利用X射線的無損檢測(NDT)技術(shù)是其重要應(yīng)用之一. zui初的X射線無損檢測方法是膠片照相法,目前這種zui原始的成像方法仍然占有主導(dǎo)地位,其原因包括以下幾個(gè)方面:①高質(zhì)量膠片可存貯的信息量大,可達(dá)到很高的空間分辨率;②根據(jù)檢測要求可以選擇膠片參數(shù),如寬容度和顆粒度等;③人們具有長期使用膠片所積累的經(jīng)驗(yàn);④可以用數(shù)字化設(shè)備,如激光掃描儀和計(jì)算機(jī)接口的相機(jī)等,將膠片的黑度轉(zhuǎn)化為數(shù)字化圖像. 如ARCUSⅡ型AGFA掃描儀可掃描的膠片尺寸為255 mm×205 mm,動(dòng)態(tài)范圍zui大為12bits. 但是,膠片照相并不能滿足實(shí)時(shí)成像,實(shí)時(shí)檢測與評估的要求,這是其zui主要的缺點(diǎn). 而且,由于膠片照相需要高質(zhì)量的膠片,因此成本高. 在膠片的保存和管理上,又受保存年限的限制并且不象電子檔案那樣易于管理. 本世紀(jì)50年代X射線像增強(qiáng)管和電視監(jiān)視器的產(chǎn)生,出現(xiàn)了電視屏直接觀察法,它是實(shí)時(shí)成像的先導(dǎo).
至目前為止,利用X射線的無損檢測方法有許多,如熒光透視法,工業(yè)X射線電視法,閃光射線照相法,數(shù)字式射線照相法,射線實(shí)時(shí)成像,線掃描成像,自動(dòng)射線照相法(Autoradiography),立體成像,背散射檢測法(Backscatter Detection),計(jì)算機(jī)層析技術(shù)等等(亦用于其它射線檢測). 隨著計(jì)算機(jī)的發(fā)展,在圖像的獲取,分析和存貯上,包括圖像質(zhì)量的改善,圖像重建,自動(dòng)分析和識(shí)別以及數(shù)字圖像存貯技術(shù),都有很大進(jìn)展. 數(shù)字式X射線成像技術(shù)是在計(jì)算機(jī)技術(shù)的基礎(chǔ)上發(fā)展起來的,并且具有很大的發(fā)展?jié)摿Γㄟ^使用高靈敏度的CCD相機(jī)和X射線轉(zhuǎn)換屏,可實(shí)現(xiàn)高質(zhì)量的無損檢測,其動(dòng)態(tài)范圍可達(dá)到3000以上,空間分辨率可達(dá)25 lp/mm,厚度靈敏度可達(dá)0.1%或更高[1]. 而基于真空管的X射線實(shí)時(shí)成像系統(tǒng),其動(dòng)態(tài)范圍低于250,因此,反差靈敏度很低,難用于對小缺陷的檢測. 使用像增強(qiáng)管的視頻相機(jī)系統(tǒng),其空間分辨率只有大約(3~4) lp/mm,并且存在圖像畸形和陰影. 數(shù)字式X射線成像技術(shù)不使用膠片,而是利用可與計(jì)算機(jī)接口的CCD相機(jī)將X射線轉(zhuǎn)換屏上的可見光轉(zhuǎn)換成數(shù)字圖像,然后對圖像自動(dòng)分析和識(shí)別,大大提高了檢測效率,并將所建立的圖像檔案存貯于光盤上,這對X射線的無損檢測和評估具有十分重要的意義,并將成為今后的發(fā)展趨勢[2]. 現(xiàn)就作者曾從事的實(shí)際工作作一評述.
。
1 數(shù)字式X射線成像技術(shù)
數(shù)字式X射線成像系統(tǒng)可有兩種結(jié)構(gòu)形式[1],如圖1,圖2所示. 該系統(tǒng)的主要組成部分是X射線到可見光的轉(zhuǎn)換屏和CCD相機(jī),其性能直接影響成像質(zhì)量. 圖1的結(jié)構(gòu)是通過反射鏡用CCD相機(jī)獲取圖像,反射鏡防止了X射線對CCD相機(jī)的直接照射所造成的噪聲增大和縮短CCD相機(jī)壽命的影響,這種結(jié)構(gòu)形式可對視場大小進(jìn)行調(diào)節(jié),因此可以滿足對較大試件的檢測. 圖2的結(jié)構(gòu)用光纖直徑為(6~16) μm的光纖耦合,將可見光轉(zhuǎn)換屏與CCD相機(jī)進(jìn)行直接耦合,而不需要鏡頭,其等效數(shù)值孔徑優(yōu)于使用透鏡轉(zhuǎn)像,但有不可變視場的缺點(diǎn)
1——X 射線源;2——準(zhǔn)直;3——物體;4——不透明層;
5——可見光轉(zhuǎn)換屏;6——X 射線與光防護(hù)罩;
7——防護(hù)罩;8——反射鏡;9——CCD 相機(jī)
圖 1 使用反射鏡的數(shù)字式X射線成像系統(tǒng)[1]
1——X 射線源;2——準(zhǔn)直;3——物體;4——X 射線閘門;
5——X 射線防護(hù)罩;6——可見光轉(zhuǎn)換屏;7——X 射線防護(hù)
8——光纖耦合器;9——CCD 相機(jī);10——X 射線與光防護(hù);11——不透明層
圖 2 使用光纖耦合的數(shù)字式X射線成像系統(tǒng)[1]
性能優(yōu)良的X射線轉(zhuǎn)換成像屏在數(shù)字式X射線成像系統(tǒng)中起著非常重要的作用,可以說它是系統(tǒng)成像質(zhì)量的瓶頸. 目前已有許多種,包括不同材料和結(jié)構(gòu)形式的屏,可應(yīng)用于從低能到MeV以上的高能X射線檢測. 使用真空管的圖像增強(qiáng)器系統(tǒng)是較早的一種實(shí)時(shí)成像系統(tǒng),只有(3~4) lp/mm的空間分辨率,后來在使用CCD相機(jī)的實(shí)時(shí)成像系統(tǒng)中,使用較多的是多晶熒光屏來作為從X射線到可見光的轉(zhuǎn)換屏,這種屏由熒光物微晶粒與有機(jī)粘合劑混合并涂于塑料襯板而成,由于成像過程中粒子間光的散射,使其空間分辨率降低. 而玻璃屏 (如鋱玻璃屏\\,硫氧化釓屏等) 不存在這一問題. 新型閃爍晶體可以提高X射線的轉(zhuǎn)化效率,高密度的玻璃閃爍體能夠改善X射線的衰減損失率,提高空間分辨率和信噪比,它比通常所用的顆粒狀熒光屏更為*,它所產(chǎn)生的光比磷光體的光更強(qiáng). 在玻璃屏中,硫氧化釓屏對X射線的轉(zhuǎn)換率比鋱玻璃屏要高,使用它時(shí),空間分辨度可達(dá)(12~13) lp/mm,但在高能X射線下,由于二次散射其空間分辨率不超過3 lp/mm.
幾年前,美國加利福尼亞Lawrence Livermore國家實(shí)驗(yàn)室研究人員研制出密度為3.75 g/cm3的鋱激活硅酸釓玻璃屏[13]在150 kV射線能量下,空間分辨率可達(dá)(20~25) lp/mm,在9 MeV射線能量下達(dá)到2 lp/mm. 加利福尼亞Lockheed公司無損檢測實(shí)驗(yàn)室研制的新型玻璃纖維閃爍屏[12]可大大提高對X射線的吸收,并且可降低余輝和二次散射的影響. 尤其是在1 MeV以上的高能X射線下對降低二次散射提高空間分辨率能夠起到一定的效果. 這種屏比硫氧化釓屏的空間分辨率更高,可達(dá)到20 lp/mm,反差靈敏度為0.2%~0.5%,在1 MeV以上的X射線下空間分辨率為3 lp/mm,密度范圍在(2.85~3.83) g/cm3之間,以圖1和圖2兩種方式工作. 該公司通過美國*Wright實(shí)驗(yàn)室利用光纖閃爍屏(Fiber-optic Scintillator)研制的半自動(dòng)高分辨率實(shí)時(shí)成像檢測系統(tǒng)[14],在對飛機(jī)構(gòu)件中的腐蝕和裂紋等缺陷以及渦輪發(fā)動(dòng)機(jī)和復(fù)合材料的缺陷檢測中取得了很好的效果,對一件較大的飛機(jī)構(gòu)件可在(10~30) s內(nèi)完成檢測. 但是,該檢測系統(tǒng)于低能X射線的條件,對大型構(gòu)件的高能X射線高分辨率的實(shí)時(shí)成像系統(tǒng)目前尚未見報(bào)導(dǎo),而在高能下成像屏的空間分辨率在很大程度上影響整個(gè)檢測系統(tǒng)的成像質(zhì)量.
目前,CCD技術(shù)已取得很大發(fā)展,具有大規(guī)模像元的科學(xué)級(jí)數(shù)字CCD相機(jī),圖像規(guī)格包括1024×1024×12 bits,2048×2048×12 bits,3072×2048×12 bits和4096×4096×12 bits像元數(shù),有些成像面的面積可達(dá)2000 mm2或更高,可直接與閃爍屏配合使用,也可以通過光纖耦合或標(biāo)準(zhǔn)鏡頭連接使用. 使用這種相機(jī),可獲得*的分辨率,幾乎與膠片接近. CCD的動(dòng)態(tài)范圍是在一個(gè)CCD像元中滿勢阱電子容量與阱底噪聲電子之比,因此希望CCD像元中滿勢阱電子多,而阱底噪聲電子少. 若CCD滿勢阱容量為85000個(gè)電子,阱底噪聲電子為17,由此計(jì)算出的動(dòng)態(tài)范圍是5000∶1,與此相應(yīng)地由于材料密度或厚度引起的zui大反差靈敏度為0.02%,*可以滿足工業(yè)檢測的需要. 表1給出了三種型號(hào)的CCD相機(jī)的性能,我們目前使用的CCD相機(jī)是PXL6300型,這種相機(jī)可由軟件設(shè)置不同的讀出速率,冷卻溫度和增益,并且像元多,暗電流極小.
表1 三種型號(hào)的CCD相機(jī)的性能
型號(hào) | Thomson | Tektronix | PXL 6300 |
格式(像元數(shù)) | 1 024×1 024 | 2 048×2 048 | 3 072×2 048 |
像元尺寸/μm\\+2 | 19×19 | 24×24 | 9×9 |
滿勢阱電子數(shù) e | 290 000 | 235 000 | 85 000 |
讀出噪聲/e | 15 | 28 | 17 |
讀出速率/(kPixels*s\\+\\) | 200 | 500 | 500~2 000(可軟件選定) |
暗電流/(e/Pixel*s\\+\\) | 16.6 | 1.0 | 0.012 |
A/D 轉(zhuǎn)換 bits | 14 | 12 | 12 |
冷卻溫度/℃ | -45 | -40 | -5~-45(可軟件選定) |
在簡單的CCD系統(tǒng)中,CCD經(jīng)快速掃描后傳輸出來,如果傳輸給電視系統(tǒng),一般曝光時(shí)間為每幀33 ms,也可以經(jīng)A/D變換電路轉(zhuǎn)化成數(shù)字量存儲(chǔ),顯示和處理. 對于一般的CCD相機(jī),通常在給定的曝光時(shí)間內(nèi)所產(chǎn)生的電子數(shù)達(dá)不到滿勢阱電子,并且成像中的噪聲很高,如行,場的同步抖動(dòng)噪聲,熱噪聲,量化噪聲等. 而讀出放大器的噪聲是隨CCD的掃描速率而增加的. 如果要檢測被檢物厚度或密度的微小變化,那么就應(yīng)增加光量子的數(shù)量,在CCD內(nèi)增加有用電子的數(shù)量提高信噪比. 因此,在提高成像質(zhì)量上,可以采用增加CCD內(nèi)的累積電荷量,增加曝光來改善噪聲特性. 另外一種提高成像質(zhì)量的方法是采用多幅疊加,可以連續(xù)采集十幾幅到幾百幅的圖像進(jìn)行疊加(可由軟件或硬件實(shí)現(xiàn)),它能夠?qū)⑿旁氡忍岣?/span>N倍(N為疊加幅數(shù)),以改善圖像質(zhì)量. 但是,在疊加時(shí),只有當(dāng)每幀中的像元達(dá)到或接近滿勢阱時(shí)才有效,否則應(yīng)當(dāng)采用延長曝光時(shí)間來積累電荷. 另外,用模擬疊加比用數(shù)字疊加所獲圖像的精度更高,因?yàn)閿?shù)字疊加有量化誤差的影響. 這些方法在滿足實(shí)時(shí)性要求的情況下,對提高成像質(zhì)量是很有效的.
。
2 X射線檢測中的數(shù)字圖像處理技術(shù)
在60年代早期,人們就開始研究射線檢測中的圖像處理技術(shù),從60年代中期至70年代,首先應(yīng)用圖像處理算法來解決醫(yī)學(xué)成像檢測問題. 這些算法在其它射線檢測領(lǐng)域中也是非常實(shí)用的. 在70年代后期至80年代初期,隨著膠片圖像數(shù)字化設(shè)備的發(fā)展及微處理器的不斷成熟,射線檢測的數(shù)字圖像處理技術(shù)取得了很大的進(jìn)展. 到80年代中期,在獲取圖像的途徑上,全面開發(fā)了膠片數(shù)字化系統(tǒng)和數(shù)字射線成像技術(shù),并使其得到了實(shí)際應(yīng)用,數(shù)字化射線檢測已占到了總檢測量的26%[3]. 在此期間,主要的問題是如何獲得高分辨率的灰度圖像. 人們從電視系統(tǒng)中可獲得的數(shù)字化圖像陣列達(dá)到了1024×1024×8或×10bits. 隨著超大規(guī)模集成電路的發(fā)展,大大提高了計(jì)算機(jī)的運(yùn)算速度,存貯容量和顯示分辨率,而半導(dǎo)體 CCD 器件的發(fā)展又為數(shù)字式射線檢測技術(shù)注入了新的活力,使相機(jī)像元可達(dá)4096×4096×12bits或更高,并具有*的信噪比,在此期間圖像處理算法也層出不窮. 使得此項(xiàng)技術(shù)出現(xiàn)了飛躍發(fā)展.
圖像增強(qiáng)技術(shù)是使處理過的圖像看上去比未處理的圖像更加能夠顯示圖像的固有特征,它可以改善人或機(jī)器對一幅圖像的視覺效果. 目前還沒有一種準(zhǔn)則來定義所謂“圖像”效果,這是因?yàn)檫€沒有一種統(tǒng)一的圖像增強(qiáng)的理論來衡量所謂“圖像”,這必須根據(jù)某種準(zhǔn)則對實(shí)際對象采取不同方法. 在對射線檢測的圖像進(jìn)行處理時(shí),由于所檢測的缺陷類型,位置,尺寸,成像條件不同等諸因素的影響,造成圖像的視覺效果是千差萬別的,因此,圖像增強(qiáng)技術(shù)是針對對象的處理技術(shù). 一般來說,圖像增強(qiáng)技術(shù)可用來處理圖像對比度低,灰度分布不均或椒鹽噪聲,量子起伏噪聲等等.
在圖像增強(qiáng)技術(shù)中,常用點(diǎn)處理技術(shù),如灰度反轉(zhuǎn),灰度拉伸,按函數(shù)要求的直方圖修正等等. 另外,還通過參考某像素鄰域的灰度來校正該像素的灰度. 如果有噪聲的像素位置是已知的,那么就可以對其鄰域像素來加權(quán)平均來取代該噪聲的灰度. 而實(shí)際上,這種信息一般是未知的,因此,這種平均化方法可以用掩模窗去處理整幅圖像,這種運(yùn)算是用窗與圖像進(jìn)行卷積,因此,可在空間域或頻率域中處理,這種中值濾波的方法對噪聲的處理往往是有效的[4].
對于射線檢測中的圖像增強(qiáng),一般來說,下列幾種方法是很有效的:
(1)圖像的灰度拉伸和直方圖均衡可以擴(kuò)展圖像的動(dòng)態(tài)范圍,提高對比度;
(2)多幀采集后經(jīng)疊加和平均可以消除隨機(jī)噪聲,提高信噪比;
(3)高通濾波來突出圖像中缺陷的細(xì)節(jié);
(4)對比度和亮度的調(diào)節(jié),以及偽彩色處理可以改善視覺效果.
通過上述4種處理可以改善圖像質(zhì)量,即能夠提高圖像的對比度,清晰度和降低噪聲[5].
數(shù)字圖像恢復(fù)技術(shù)起始于60年代末,到達(dá)70年末期已經(jīng)形成較為成熟的一門圖像處理技術(shù). 美國科學(xué)家H C Andrews和B R Hunt對本領(lǐng)域發(fā)展作出了杰出的貢獻(xiàn)[6]. 圖像恢復(fù)的方法有許許多多,其中有兩種zui基本的方法,即逆濾波法和維納濾波法. 兩種方法相比,逆濾波法更為直接,簡單,但它的zui大缺點(diǎn)是容易出現(xiàn)病態(tài)解,若出現(xiàn)這種情況,則恢復(fù)出來的效果極差. 為克服這種現(xiàn)象,B R Hunt建議[7]事先研究一下光學(xué)傳遞函數(shù)是否存在零點(diǎn),若有零點(diǎn)則對其進(jìn)行校正,則可獲得較好的效果. 維納濾波恢復(fù)則不會(huì)出現(xiàn)病態(tài)解,但需要計(jì)算噪聲和信號(hào)的功率譜,因此計(jì)算量較大. B R Hunt指出[7],有75%的圖像恢復(fù)問題可以用這兩種zui簡單辦法來解決,而剩下的25%的恢復(fù)問題需用其它方法. 圖像恢復(fù)即為從降質(zhì)圖像中估計(jì)恢復(fù)圖像原貌的過程. 圖像恢復(fù)包括消除以下兩個(gè)方面的圖像降質(zhì)因素:
(1)系統(tǒng)降質(zhì)因素,如光學(xué)散焦,大氣湍流和運(yùn)動(dòng)模糊. 在射線檢測中的系統(tǒng)降質(zhì)因素,有由于射線源尺寸有限大小造成的圖像幾何不清晰度以及由于試件缺陷方位與射線方向有一定偏角造成的幾何畸變等;
(2)由于噪聲所引起的圖像統(tǒng)計(jì)規(guī)律上的降質(zhì)因素,如量子起伏噪聲,射線散射噪聲和量化噪聲等. 在射線檢測中,散射線對圖像的降質(zhì)是zui為突出的因素之一.
在圖像恢復(fù)過程中為消除這兩種形式的降質(zhì)因素所設(shè)計(jì)的濾波器通常是相沖突的,因?yàn)橄S機(jī)噪聲一般是低通過程,因此很難保證圖像的細(xì)節(jié)(通常含有較強(qiáng)的高頻成份)信息不受損失. 因而,高性能的濾波器是這兩方面的合理折衷. 圖像恢復(fù)的質(zhì)量在很大程度上取決于以下兩方面的因素:
。
(1)成像系統(tǒng)建模應(yīng)盡可能準(zhǔn)確,包括對系統(tǒng)的線性與非線性分析和系統(tǒng)空間頻率響應(yīng)移變性的分析,用理論和實(shí)驗(yàn)方法確定各個(gè)環(huán)節(jié)以及系統(tǒng)的點(diǎn)擴(kuò)展函數(shù),另外還要考慮到系統(tǒng)統(tǒng)計(jì)規(guī)律所引起的降質(zhì)因素,如研究隨機(jī)噪聲的統(tǒng)計(jì)特征,設(shè)計(jì)濾波器類型和參數(shù);
。
(2)根據(jù)觀察者或視覺系統(tǒng)確定圖像質(zhì)量準(zhǔn)則.
針對射線成像系統(tǒng)的圖像恢復(fù)方法,很重要的一個(gè)環(huán)節(jié)是確定系統(tǒng)點(diǎn)擴(kuò)展函數(shù)(PSF),它反映了成像系統(tǒng)的傳遞函數(shù),即系統(tǒng)的特性. 而決定點(diǎn)擴(kuò)展函數(shù)因素主要來源于兩個(gè)方面:一方面射線源的有限尺寸造成的幾何不清晰度和射線方向與缺陷存在的夾角造成的幾何畸變及圖像灰度分布的不均勻性;另一方面則是由于散射線對成像質(zhì)量的影響. 使得成像系統(tǒng)的對比度和清晰度降低. 對于這兩方面所造成的圖像降質(zhì)的恢復(fù)問題從 70 年代就開始了研究\\+\\{ 8,9\\]\\}. 對于高能X射線檢測,散射影響是zui突出的因素. 目前,在高能X射線下的圖像恢復(fù)研究已取得較好的效果[10,11],但針對不同對象和條件的圖像恢復(fù)問題有待進(jìn)一步研究和探索.
作者單位:華北工學(xué)院 電子信息工程系,山西 太原 030051
參考文獻(xiàn):
1 Bueno C,Barker M D. High-resolution digital radiography and three-dimentional computed tomography. X-Ray Detector Physics Applications Ⅱ,SPIE,Vol 2009:179~191.
2 Berger H. Trends in Radiologic NDT. Material Evaluation,Nov,1994,1248~1251.
3 DwyerⅢ S J. Radiographic image processing. Digital Image Processing,SPIE,Vol 528:134~136.
4 Huang T S. A fast two-dimension median filtering algorithm. IEEE Trans. on Acoustics,Speech,and Signal Processing,ASSP-27,1979,(1):13~18.
5 Burch S F. Digital enhancement of video images for NDT. NDT International,1987,20(1).
6 Anddrews H C,Hunt B R. Digital Image Restoration. Englewood Ciffs,NT:Prentice-Hall,1979.
7 Hunt B R. Image restoration. Digital Image Processing Techniques. Edited by Ekstrom. M P,Academic Press Inc,1984.
8 Hunt B R. Introduction to Restoration and Enhancement of Radiographic Images. Report No. LA 430,Los Alamos Scientific Lab,Los Alamos,New Mexico,1970.
9 Krusos G A. Restoration of radiologic images by optical spatial filtering. Optical Engineering,1974,13(3):208-218.
10 韓焱. 彈藥裝藥高能X射線數(shù)字圖像處理技術(shù):[博士學(xué)位論文]. 北京:北京理工大學(xué),1998.
11 劉小軍,高能X射線圖像的頻域預(yù)處理技術(shù):[碩士學(xué)位論文]. 太原:華北工學(xué)院,1998.
12 Bueno C,Buchanan R A. Luminescent glass design for high energy real-time radiography. SPIE,1990,1327:79~91.
13 Placious R C. High Density Glass Scintillator for Real-time X-Ray Inspection. Material Evaluation,November 1991.
14 Bueno C. Nodestructive evaluation of aircraft structures using high-resolution real-time radiography. SPIE,1995,2455,114~124.