基于多核CPU的微型CT并行錐束圖像重建算法應(yīng)用研究.pdf

1、20世紀80年代，普通醫(yī)學(xué)CT的空間分辨率就已達到毫米量級，但僅能粗糙地滿足宏觀病理學(xué)診斷的要求。如果醫(yī)學(xué)CT的空間分辨率提高到微米量級，即可顯示細胞和組織的微細結(jié)構(gòu)，從而真正滿足臨床病理學(xué)診斷金標準的要求。為此，學(xué)術(shù)界開始研發(fā)高分辨率的微型CT(Micro CT；也稱為顯微CT、微焦點CT)。Micro CT通常采用錐形X線束(Cone Beam)，采用錐形束不僅能夠獲得真正各向同性的容積圖像，提高空間分辨率及X射線利用率，而且在采集

2、相同3D圖像時速度遠遠快于扇形束CT。因此，錐束圖像重建算法成為Micro CT學(xué)術(shù)理論及工程技術(shù)界的研究熱點。 微型CT采用錐形X射線束(簡稱錐束)和二維檢測器取代傳統(tǒng)的扇形x射線束和一維檢測器，這樣不僅可以一次同時掃描多層組織，而且提高了獲取投影數(shù)據(jù)的速度和X射線的利用率，同時重建出的圖像的軸向分辨率也得到很大的提高。錐束圖像重建算法包括解析算法和迭代算法。其中，解析算法又分為近似算法和精確算法。相對于迭代算法，解析算法具有

3、更快的重建速度，并且占有較少的內(nèi)存資源，因此，成為實際CT系統(tǒng)中的主流算法。而在解析算法中，近似算法由于數(shù)學(xué)上簡單，實現(xiàn)起來容易，而且當(dāng)錐角比較小時，能夠取得較好的重建效果，所以在實際中有著廣泛的應(yīng)用。著名的Feldkamp算法(Kudo and Saito 1990，Yahand Leahy 1992，Wang et al 1993)就是典型的能夠處理部分錐束投影數(shù)據(jù)的近似處理算法。然而，當(dāng)檢測器陣列變寬、射束錐角變大時，其應(yīng)用性能總

4、體將明顯趨于惡化。但由于Feldkamp算法以及摹于Feldkamp算法的其它近似錐束圖像重建算法優(yōu)異的時間性能、良好的理論簡潔性和工程易實現(xiàn)性，它們?nèi)匀痪哂泻艽蟮墓こ虘?yīng)用價值，并構(gòu)成目前Micro CT工程的主流圖像重建算法，也是本課題的研究內(nèi)容。 由于CT圖像重建過程運算量大，尤其是反投影過程，在很多情況下，重建所需時間難以滿足實際應(yīng)用的需求，因此提高反投影重建速度一直是CT理論及工程技術(shù)界關(guān)注的熱點。國內(nèi)外的學(xué)者和機構(gòu)在這

5、方面做了大量的研究工作，提出很多快速反投影算法，但這些方法大都是基于正、余弦函數(shù)性質(zhì)的優(yōu)化方法，需要預(yù)先計算并存儲三角函數(shù)值，內(nèi)存開銷比較大。此外，采用并行處理技術(shù)可有效地提高圖像重建速度，縮短重建時間。本課題采用多核程序設(shè)計和并行處理技術(shù)，主要研究了并行錐束圖像重建算法，并且針對反投影問題，介紹了一種基于ADRT的快速反投影算法及其并行實現(xiàn)。雖然并行處理技術(shù)早在1990就被應(yīng)用于二維CT圖像重建，但它們都是在并行機或集群系統(tǒng)中實現(xiàn)的，

6、而在多核CPU上實現(xiàn)并行，能帶來更強的并行處理能力和更高的計算速度。 同時，仿真實驗是研究人員驗證其算法正確性和有效性、評價其算法性能的基本途徑，這在生物醫(yī)學(xué)圖像重建領(lǐng)域表現(xiàn)的尤為突出。本課題以3D Shepp-Logan模型為參考模型，在VC++6.0環(huán)境下開發(fā)微型CT錐束重建仿真實驗平臺，對各種圖像重建算法進行性能評測。該實驗平臺一方面能夠保障本課題算法的順利實現(xiàn)，更重要的，它是今后Micro CT工程技術(shù)人員必不可少的科研