喬佳亮，王明泉，韓得水

（中北大學(xué)a.動態(tài)測試技術(shù)重點實驗室;b.中北大學(xué)信息與通信工程學(xué)院，山西太原 030051）

基于多水平系統(tǒng)的錐束ART算法

喬佳亮a，b，王明泉a，b，韓得水a(chǎn)，b

（中北大學(xué)a.動態(tài)測試技術(shù)重點實驗室;b.中北大學(xué)信息與通信工程學(xué)院，山西太原 030051）

代數(shù)重建技術(shù)（ART）是X－射線CT圖像重建的一種重要迭代算法，投影的訪問順序?qū)ζ涫諗克俣扔兄匾绊憽?994年，Guan和Gordon提出了多水平訪問系統(tǒng)（MLS），該系統(tǒng)在投影數(shù)為2的整數(shù)次冪時效果最好。針對MLS僅適用于［0，π）上的均勻投影的問題，將MLS推廣到［0，2π）上，在特定的投影數(shù)情形下，分析投影數(shù)特點，提出了三等分方法，該方法對對稱投影訪問順序進(jìn)行了改進(jìn)。仿真實驗結(jié)果表明，推廣后的兩系統(tǒng)下的錐束ART算法在迭代中收斂速度都要快于隨機訪問系統(tǒng)（RAS）和順序訪問系統(tǒng)（SAS）下的錐束ART算法。

錐束ART算法;CT圖像重建;多水平訪問系統(tǒng);投影

ART算法是計算機斷層成像技術(shù)（CT）中用于圖像重建的一種重要方法，因其具有采用不完全投影數(shù)據(jù)也能實現(xiàn)較高成像質(zhì)量的優(yōu)點，在大型工業(yè)檢測和地質(zhì)勘探中應(yīng)用廣泛［1］。ART 主要缺點是計算量大，重建速度慢［2］。對此，很多學(xué)者從投影訪問順序、松弛因子的選擇和并行計算等方面改進(jìn)ART算法，提高迭代速度［3］。對于ART算法中的投影訪問順序系統(tǒng)，常用的有隨機訪問系統(tǒng)（RAS）［4］、多水平訪問系統(tǒng) （MLS）［5］等。

針對MLS僅適用于［0，π）上的均勻投影的問題，本文將MLS推廣到［π，2π）上，首先將其在［0，π）上的順序?qū)ΨQ地應(yīng)用在［π，2π）上，獲得［π，2π）上的投影訪問順序，鑒于MLS在投影數(shù)為2的整數(shù)次冪時效果最好，在360個投影角度數(shù)情形下，本文采用三等分法對對稱投影訪問順序進(jìn)行了改進(jìn)。仿真實驗結(jié)果表明，改進(jìn)后的投影訪問順序錐束ART算法在迭代中收斂速度要快于隨機訪問系統(tǒng)（RAS）和順序訪問系統(tǒng)（SAS）下的錐束ART算法。

1 多水平系統(tǒng)介紹

多水平系統(tǒng)（MLS）是一種較好的投影訪問順序系統(tǒng)，當(dāng)投影方向數(shù)為2的冪次方時，系統(tǒng)具有很好的效果［6］。設(shè)投影方向總數(shù)為2的整數(shù)次冪，首先被訪問的是第一層的數(shù)據(jù)，即投影在0（0°）和M/2（90°）方向最大程度正交。第三個被訪問的是第二層的數(shù)據(jù)，此時投影角為M/4（45°），然后是3M/4（135°）。以此類推，接下來被訪問的投影角分別為:M/8，5M/8，3M/8，7M/8等。表1給出了L=4時各層的訪問順序。

表1 MLS中M個投影方向的訪問順序

當(dāng)投影方向總數(shù)不是2的整數(shù)次冪，即M≠2L時，也可以使用MLS順序。在第一層時，可選任意兩個成90°（或接近的）夾角的投影角度，在第二層時選取上述兩個方向夾角的角平分線位置（或接近的位置），依次類推獲得相應(yīng)的投影角度。投影角度的訪問順序繼續(xù)依據(jù)表1，總層數(shù)為L=［lbM］+1（［］表示取整）［7］。

2 ［0，2π）的多水平系統(tǒng)的ART算法

2.1 對稱推廣

考慮投影射線的相關(guān)性，在［0，π）上MLS是一個效果很好的投影順序訪問系統(tǒng)。將［0，π）上的MLS推廣到［0，2π），可以將［0，π）上的 MLS 訪問順序中心對稱用在［π，2π）上，該順序系統(tǒng)在［0，π）和［π，2π）上分別達(dá)到了MLS的效果，本文稱此系統(tǒng)為MLS1。設(shè)投影總數(shù)為2的整數(shù)次冪。表2列出了L=5時MLS1各層的訪問順序。圖1給出了在［0，2π）上均勻分布的16個投影角度的MLS1訪問順序示意圖。當(dāng)投影總數(shù)不是2的整數(shù)次冪時，處理方法與［0，π）上MLS的處理方法相同。

表2 MLS1下M個投影的訪問順序

圖1 在［0，2π）上均勻分布的16個投影角度的MLS1訪問順序

2.2 改進(jìn)算法

圖2 在［0，2π）上均勻分布的12個投影角度的MLS2訪問順序

3 實驗結(jié)果

將MLS1，MLS2，RAS，SAS應(yīng)用到三維物體重建。在實驗中，選用三維的Shepp－Logan模型，大小為128×128×128，設(shè)置參數(shù)為物距650，相距350，λ=0.2，投影射線數(shù)目為360×300×300，其中360為投影角度數(shù)，每隔1°取一次投影，300×300為探測器大小。

圖3給出了4個系統(tǒng)下錐束ART算法一次迭代重建所得物體分別在x=0，y=0，z=0處的截面圖，其中第1列為x=0截面圖，第2列為y=0的截面圖，第3列為z=0的截面圖。與［0，π）上的重建一樣，RAS的迭代中每次所用的隨機數(shù)組都是不同的。

可以明顯看出，推廣后的MLS1和MLS2的效果優(yōu)于SAS，為了更清楚地顯示出MLS1，MLS2與RAS之間的差異，表3、表4分別給出了4種系統(tǒng)下重建的歸一化均方距離d1、平均均方誤差d3隨迭代次數(shù)k的變化。

圖3 SAS，RAS，MLS1，MLS2 的 ART 算法一次迭代所得物體的截面圖

表3 4種系統(tǒng)下重建的歸一化均方距離d1

表4 4種系統(tǒng)下重建的平均均方誤差d3

表3、表4表明，在前5次的迭代中，MLS2優(yōu)于MLS1，MLS1優(yōu)于RAS，3個系統(tǒng)都優(yōu)于SAS。

4 結(jié)論

本文介紹了基于MLS的ART算法，并驗證了其在三維錐束ART算法中相對于RAS和SAS的優(yōu)越性。同時針對MLS僅適用于［0，π）上的均勻投影的問題，本文將MLS推廣到0，2[）π上，在特定的投影數(shù)情形下，分析投影數(shù)特點，提出了三等分方法，該方法對對稱投影訪問順序進(jìn)行了改進(jìn)。從實驗結(jié)果可以看出，無論是基于何種情況，文中的方法都優(yōu)于RAS和SAS，而且該方法在采用不完全投影數(shù)據(jù)的情況下也能成像出質(zhì)量較高的圖像，進(jìn)一步驗證了該方法確實是一種有效的圖像重建算法。

:

［1］干宗良，粱秀聚.自適應(yīng)稀疏約束圖像超分辨力重建方法［J］.電視技術(shù)，2012，36（14）:19－23.

［2］周斌，王連堂，王俊杰，等.代數(shù)重建算法中一種快速投影系數(shù)計算方法［J］.計算機工程與應(yīng)用，2008，44（25）:46－47.

［3］張順利，張定華，王凱，等.一種基于ART算法的快速圖像重建技術(shù)［J］.核電子學(xué)與探測技術(shù)，2007，27（3）:479－483.

［4］王宏鈞，路宏年，傅健.代數(shù)重建技術(shù)中投影序列選擇次序的研究［J］.光學(xué)技術(shù)，2006，32（3）:389－391.

［5］VAN D M C.Iterative methods in image reconstruction［D］.Netherlands:Utrecht University，1992.

［6］張順利，張定華，趙歆波，等.工業(yè)CT圖像重建的ART算法研究［J］.無損檢測，2007，29（8）:453－456.

［7］孫開一.計算機模擬斷層掃描技術(shù)［J］.科技資訊，2008（4）:103－104.

Multilevel Scheme for Cone－beam ART Algorithm

QIAO Jialianga，b，WANG Mingquana，b，HAN Deshuia，b

（a.Dynamic Testing Technology Key Laboratories，North University of China;b.School of Information and Communication Engineering，North University of China，Taiyuan 030051，China）

ART is one important iterative algorithm of CT image reconstruction using X－ray，while the sequential access of the projection is very important to the convergence rate.In 1994，Guan and Gordon raised multilevel scheme（MLS）which is best performance when the projection number is power of two.According to MLS only applies for the circumstance when projection is uniform in［0，π），then extend it to［0，2π）.Analyzing the features of the projection number in the certain projection number，the trisection method is proposed in this paper.This method is to promote the first method with trisection.The results of experiments show that convergence rate of cone－beam ART algorithms with the above two systems perform better than that with the RAS and SAS.

cone－beam ART algorithm;CT image reconstruction;MLS;projection

TN911.73

A

【本文獻(xiàn)信息】喬佳亮，王明泉，韓得水.基于多水平系統(tǒng)的錐束ART算法［J］.電視技術(shù)，2013，37（13）.

國家自然科學(xué)基金項目（6171177）;山西省工業(yè)攻關(guān)項目基金項目（20110321073）;中北大學(xué)2011青年基金項目（201261413752398）

喬佳亮（1988— ），碩士生，研究方向為信號與信息處理;

王明泉（1970— ），博士生導(dǎo)師，研究方向為信息獲取與處理技術(shù)、圖像處理與識別;

韓得水（1989— ），碩士生，研究方向為信號與信息處理。

責(zé)任編輯:時 雯

2012－10－23