基于貝耶斯方法的自發(fā)熒光斷層成像重建方法研究.pdf

1、醫(yī)學(xué)影像技術(shù)的不斷進(jìn)步，極大地推動(dòng)了現(xiàn)代醫(yī)學(xué)和生物學(xué)的發(fā)展。作為一門新興的前沿交叉學(xué)科，分子影像可以在活體狀態(tài)下對(duì)生物過(guò)程進(jìn)行細(xì)胞和分子水平的定性和定量研究，為研究生物體生長(zhǎng)發(fā)育、疾病的早期檢測(cè)、藥物定量評(píng)估及動(dòng)力學(xué)變化等提供有效的方法和工具，具有廣泛的應(yīng)用前景。自發(fā)熒光斷層成像(BLT)由于靈敏度高、無(wú)創(chuàng)傷、成本低、操作簡(jiǎn)單以及背景噪聲低等特性而成為分子影像領(lǐng)域的一個(gè)重要分支。
　　 BLT的目的是利用在生物體體表獲得的熒光圖

2、像重建熒光光源在生物體體內(nèi)的分布以反映生物體細(xì)胞分子水平變化。但熒光光子在到達(dá)生物體體表前，由于生物組織具有強(qiáng)散射、低吸收的特點(diǎn)，導(dǎo)致熒光光子在生物組織中不再沿直線傳輸，而是經(jīng)歷了大量的散射過(guò)程，使得BLT本身具有病態(tài)性，即真實(shí)光源存在但不能重建熒光光源或重建光源不唯一、不準(zhǔn)確。此外，缺少外部的激發(fā)光源和生物體的非勻質(zhì)特性進(jìn)一步加劇了熒光光源重建的難度。同時(shí)，由于在生物體體表獲得的數(shù)據(jù)有限，而要求的解數(shù)量非常大，導(dǎo)致重建光源不唯一且對(duì)噪

3、聲比較敏感。因此，如何降低BLT的病態(tài)性，唯一、準(zhǔn)確地重建熒光光源成為當(dāng)前研究的熱點(diǎn)。目前，國(guó)際上的相關(guān)研究也剛剛起步，但尚未取得突破性進(jìn)展，仍有大量的理論和技術(shù)問(wèn)題需要解決。由于BLT重建算法研究新進(jìn)展決定著B(niǎo)LT技術(shù)的應(yīng)用潛力，所以深入開(kāi)展重建算法的研究是非常有必要的。本文的主要研究工作就是圍繞重建算法展開(kāi)的，完成的主要工作包括：
　　 ①提出了基于貝耶斯方法的BLT重建算法以重建熒光光源分布。BLT光源重建是典型的病態(tài)問(wèn)題

4、，一般情況下，重建的光源不唯一，只有融入盡可能多的先驗(yàn)信息才能確保重建光源的唯一。但是如何融入更多的先驗(yàn)信息需要進(jìn)一步的研究和探索?？紤]到生物體特別是小鼠的形狀是一個(gè)不規(guī)則的，基于體素的方法對(duì)生物體進(jìn)行網(wǎng)格剖分實(shí)現(xiàn)起來(lái)比較復(fù)雜，而且在生物組織邊界處會(huì)產(chǎn)生較大的誤差。此外，假定生物體的生物組織是勻質(zhì)的，會(huì)導(dǎo)致重建光源不準(zhǔn)確。針對(duì)上述問(wèn)題，本文基于貝耶斯方法，在有限元方法的基礎(chǔ)上，提出了三維BLT重建算法。該算法首先用馬爾可夫隨機(jī)場(chǎng)先驗(yàn)?zāi)Ｐ?/p>

5、估計(jì)未知光源密度分布情況，同時(shí)融入先驗(yàn)可行光源區(qū)域信息、生物組織的結(jié)構(gòu)信息和光學(xué)特性參數(shù)信息，有效地降低了BLT的病態(tài)性。其次，建立了關(guān)于未知光源密度與噪聲模型的最大后驗(yàn)概率函數(shù)。通過(guò)對(duì)該函數(shù)優(yōu)化求解，就可重建出光源的分布。有限元方法的使用使得該算法適用于任意復(fù)雜形狀的生物體，增加了該算法的實(shí)用性。最后，利用micro-CT獲得的非勻質(zhì)小鼠和基于蒙特卡洛方法產(chǎn)生的前向仿真數(shù)據(jù)對(duì)算法有效性進(jìn)行了驗(yàn)證。真實(shí)仿體實(shí)驗(yàn)進(jìn)一步證明了該算法的可行性

6、。
　　 ②提出了一種自適應(yīng)有限元重建算法以提高重建精度和效率。該算法是在上一章算法的基礎(chǔ)上發(fā)展起來(lái)的。由于利用有限元方法進(jìn)行光源重建時(shí)存在著重建質(zhì)量和離散網(wǎng)格單元大小之間的矛盾。一般來(lái)說(shuō)，當(dāng)離散網(wǎng)格單元越細(xì)，重建結(jié)果越好，但BLT的病態(tài)性越嚴(yán)重，而且重建時(shí)間越長(zhǎng)；離散網(wǎng)格單元越粗，重建誤差越大，BLT病態(tài)性越小，重建時(shí)間也越短。自適應(yīng)有限元方法可以有效地解決這幾者之間的矛盾。在自適應(yīng)有限元方法的基礎(chǔ)上，結(jié)合貝耶斯方法可以精確地

7、重建出熒光光源，并且可以提高算法的計(jì)算效率。利用非勻質(zhì)仿體對(duì)算法進(jìn)行了有效性的驗(yàn)證，光源位置和強(qiáng)度得到了很好的重建。
　　 ③提出了一種優(yōu)化可行光源區(qū)域策略以降低BLT的病態(tài)性。在最初的研究中，利用先驗(yàn)可行光源區(qū)域降低BLT的病態(tài)性，并取得了不錯(cuò)的結(jié)果。但是，當(dāng)在生物體中存在多光源或光源位于生物體比較深時(shí)，推斷的先驗(yàn)可行光源區(qū)域可能會(huì)存在誤差。針對(duì)這個(gè)問(wèn)題，本文基于多光譜信息發(fā)展了優(yōu)化可行光源區(qū)域策略，該策略可以快速地選擇可行光

8、源區(qū)域，有效地降低BLT的病態(tài)性。最后，利用非勻質(zhì)的仿體實(shí)驗(yàn)對(duì)算法進(jìn)行了驗(yàn)證，重建的光源分布與真實(shí)光源非常接近。
　　 ④提出了融合多種先驗(yàn)信息的稀疏貝耶斯重建算法以提高重建質(zhì)量。為降低BLT的病態(tài)性，提高重建質(zhì)量，需要盡可能多的融入先驗(yàn)信息。本文提出了融入多光譜信息、優(yōu)化可行光源區(qū)域、生物組織的結(jié)構(gòu)信息、光學(xué)特性參數(shù)信息以及光源稀疏特性的重建算法，多種先驗(yàn)信息的使用有效地改善了重建結(jié)果?？紤]到計(jì)算效率，利用基于點(diǎn)的有限元方法降

9、低數(shù)據(jù)量。此外，優(yōu)化可行光源區(qū)域的使用進(jìn)一步地提高了算法的重建速度。通過(guò)micro-CT獲得的非勻質(zhì)小鼠對(duì)算法進(jìn)行了驗(yàn)證，結(jié)果表明了算法的優(yōu)勢(shì)和潛力。進(jìn)一步，利用小鼠形狀的仿體采集真實(shí)實(shí)驗(yàn)數(shù)據(jù)對(duì)算法進(jìn)行了可行性分析。
　　 ⑤本文利用主成分分析方法和基于有限內(nèi)存矩陣發(fā)展了兩種快速的重建算法。多光譜方法的引入雖然可以提高重建結(jié)果的質(zhì)量，但是也帶來(lái)了數(shù)據(jù)災(zāi)難，極大地影響了重建速度。本文在優(yōu)化可行光源區(qū)域的基礎(chǔ)上，通過(guò)利用主成分分析方