壓縮感知磁共振成像基本原理及特點(diǎn)

壓縮感知磁共振成像基本原理及特點(diǎn)一、壓縮感知理論的核心基礎(chǔ)壓縮感知（CompressedSensing,CS）是一種突破傳統(tǒng)奈奎斯特采樣定理限制的信號采集與重構(gòu)理論，由Donoho、Candes等學(xué)者在2006年提出。其核心思想是：如果信號在某個(gè)變換域是稀疏的，那么就可以用遠(yuǎn)低于奈奎斯特采樣率的方式對信號進(jìn)行非自適應(yīng)采樣，并且能夠通過非線性重構(gòu)算法準(zhǔn)確地恢復(fù)出原始信號。這一理論為磁共振成像（MagneticResonanceImaging,MRI）的加速采集提供了重要的理論支撐。（一）信號稀疏性信號稀疏性是壓縮感知理論的前提條件。所謂稀疏性，是指信號可以通過某種變換（如傅里葉變換、小波變換、曲波變換等）表示為少數(shù)非零系數(shù)的線性組合。例如，一幅自然圖像在小波變換域中，大部分系數(shù)的絕對值都很小，只有少數(shù)大系數(shù)包含了圖像的主要信息，這樣的圖像就具有稀疏性。在磁共振成像中，人體組織的磁共振信號在k空間（傅里葉空間）通常具有一定的稀疏性，尤其是在經(jīng)過適當(dāng)?shù)淖儞Q后，這種稀疏性會更加明顯。例如，腦部磁共振圖像在小波變換域或壓縮感知變換域（如有限角度變換、局部余弦變換等）中可以被稀疏表示，這為壓縮感知在MRI中的應(yīng)用提供了可能。（二）非自適應(yīng)采樣傳統(tǒng)的采樣方式通常是基于奈奎斯特采樣定理，即采樣頻率必須大于信號最高頻率的兩倍，才能保證從采樣信號中無失真地恢復(fù)出原始信號。而壓縮感知采用的是非自適應(yīng)采樣，即采樣過程不需要根據(jù)信號的具體特征進(jìn)行調(diào)整，而是通過設(shè)計(jì)一種與稀疏變換基不相關(guān)的測量矩陣，對信號進(jìn)行隨機(jī)或偽隨機(jī)采樣。在磁共振成像中，傳統(tǒng)的k空間采樣是按照笛卡爾網(wǎng)格進(jìn)行的，需要采集大量的采樣點(diǎn)才能保證圖像的質(zhì)量。而壓縮感知MRI則采用非笛卡爾采樣方式，如隨機(jī)采樣、偽隨機(jī)采樣、螺旋采樣、徑向采樣等，只采集k空間中的部分采樣點(diǎn)，從而實(shí)現(xiàn)加速采集的目的。（三）非線性重構(gòu)算法由于壓縮感知采用了遠(yuǎn)低于奈奎斯特采樣率的采樣方式，因此不能直接使用傳統(tǒng)的線性重構(gòu)方法（如逆傅里葉變換）來恢復(fù)原始信號。壓縮感知的重構(gòu)過程是一個(gè)非線性優(yōu)化問題，需要通過求解一個(gè)凸優(yōu)化問題來實(shí)現(xiàn)信號的準(zhǔn)確恢復(fù)。常用的重構(gòu)算法包括基追蹤（BasisPursuit,BP）、正交匹配追蹤（OrthogonalMatchingPursuit,OMP）、壓縮匹配追蹤（CompressiveSamplingMatchingPursuit,CoSaMP）等。這些算法的基本思想是在稀疏變換域中尋找一個(gè)稀疏解，使得該解經(jīng)過測量矩陣變換后與測量值之間的誤差最小。在磁共振成像中，重構(gòu)算法的性能直接影響到圖像的質(zhì)量和重構(gòu)速度，因此選擇合適的重構(gòu)算法是壓縮感知MRI的關(guān)鍵環(huán)節(jié)之一。二、壓縮感知磁共振成像的基本原理壓縮感知磁共振成像（CS-MRI）是將壓縮感知理論應(yīng)用于磁共振成像領(lǐng)域的一種新技術(shù)，其基本原理是利用磁共振信號的稀疏性，通過非自適應(yīng)采樣方式采集k空間中的部分采樣點(diǎn)，然后利用非線性重構(gòu)算法從這些采樣點(diǎn)中恢復(fù)出完整的磁共振圖像。具體來說，CS-MRI的實(shí)現(xiàn)過程主要包括以下幾個(gè)步驟：（一）信號稀疏表示首先，需要對磁共振信號進(jìn)行稀疏表示。這通常通過選擇合適的稀疏變換基來實(shí)現(xiàn)，如小波變換、曲波變換、有限角度變換、局部余弦變換等。這些變換基能夠?qū)⒋殴舱裥盘栟D(zhuǎn)換到一個(gè)稀疏域，使得信號在該域中只有少數(shù)非零系數(shù)。例如，曲波變換具有良好的方向選擇性和各向異性，能夠有效地表示圖像中的邊緣和紋理信息，因此在磁共振圖像的稀疏表示中得到了廣泛的應(yīng)用。在實(shí)際應(yīng)用中，需要根據(jù)不同的成像部位和成像目的選擇合適的稀疏變換基，以獲得最佳的稀疏表示效果。（二）k空間采樣在傳統(tǒng)的MRI中，k空間采樣是按照笛卡爾網(wǎng)格進(jìn)行的，需要采集大量的采樣點(diǎn)才能保證圖像的質(zhì)量。而在CS-MRI中，采用非自適應(yīng)采樣方式采集k空間中的部分采樣點(diǎn)。常用的采樣方式包括隨機(jī)采樣、偽隨機(jī)采樣、螺旋采樣、徑向采樣等。其中，隨機(jī)采樣是最常用的一種采樣方式，它通過隨機(jī)選擇k空間中的采樣點(diǎn)來實(shí)現(xiàn)加速采集。偽隨機(jī)采樣則是基于某種偽隨機(jī)序列（如混沌序列、m序列等）來選擇采樣點(diǎn)，具有更好的可控性和重復(fù)性。螺旋采樣和徑向采樣則是通過改變采樣軌跡來實(shí)現(xiàn)加速采集，這些采樣方式能夠減少k空間采樣的冗余性，提高采樣效率。在選擇采樣方式時(shí)，需要考慮采樣點(diǎn)的分布、采樣密度、采樣軌跡等因素，以保證采樣信號能夠包含足夠的圖像信息，同時(shí)實(shí)現(xiàn)較高的加速倍數(shù)。（三）測量矩陣設(shè)計(jì)測量矩陣是壓縮感知理論中的一個(gè)重要組成部分，它用于將原始信號投影到一個(gè)低維空間中，得到測量值。在CS-MRI中，測量矩陣通常與k空間采樣方式相對應(yīng)。例如，當(dāng)采用隨機(jī)采樣方式時(shí)，測量矩陣是一個(gè)隨機(jī)矩陣，其中的元素只有0和1，1表示對應(yīng)的k空間采樣點(diǎn)被采集，0表示未被采集。測量矩陣需要滿足約束等距性（RestrictedIsometryProperty,RIP），即對于任意稀疏信號，測量矩陣在該信號的稀疏變換域中能夠保持信號的線性結(jié)構(gòu)，從而保證信號能夠被準(zhǔn)確重構(gòu)。在實(shí)際應(yīng)用中，由于測量矩陣的規(guī)模通常很大，直接驗(yàn)證其是否滿足RIP是非常困難的。因此，通常通過選擇合適的采樣方式和稀疏變換基，來間接保證測量矩陣滿足RIP條件。（四）圖像重構(gòu)圖像重構(gòu)是CS-MRI的核心環(huán)節(jié)，它通過求解一個(gè)凸優(yōu)化問題，從測量值中恢復(fù)出完整的磁共振圖像。常用的重構(gòu)算法可以分為兩類：一類是基于L1范數(shù)最小化的凸優(yōu)化算法，如基追蹤（BP）、梯度投影稀疏重構(gòu)（GradientProjectionforSparseReconstruction,GPSR）等；另一類是基于貪婪算法的迭代重構(gòu)算法，如正交匹配追蹤（OMP）、壓縮匹配追蹤（CoSaMP）、正則化正交匹配追蹤（RegularizedOrthogonalMatchingPursuit,ROMP）等。基于L1范數(shù)最小化的凸優(yōu)化算法的基本思想是將圖像重構(gòu)問題轉(zhuǎn)化為一個(gè)L1范數(shù)最小化問題，即在滿足測量值約束的條件下，尋找一個(gè)稀疏解，使得該解的L1范數(shù)最小。這類算法通常具有較高的重構(gòu)精度，但計(jì)算復(fù)雜度較高，重構(gòu)速度較慢。例如，基追蹤算法通過求解一個(gè)線性規(guī)劃問題來實(shí)現(xiàn)圖像重構(gòu)，它能夠保證在滿足一定條件下得到信號的稀疏解，但需要大量的計(jì)算時(shí)間和內(nèi)存資源。基于貪婪算法的迭代重構(gòu)算法則是通過迭代選擇與測量值最相關(guān)的原子（稀疏變換基中的元素），逐步逼近信號的稀疏表示。這類算法的計(jì)算復(fù)雜度較低，重構(gòu)速度較快，但重構(gòu)精度相對較低。例如，正交匹配追蹤算法通過每次選擇與殘差最相關(guān)的原子，并將其加入到支撐集中，然后通過最小二乘法求解支撐集中的系數(shù)，不斷迭代直到滿足停止條件。壓縮匹配追蹤算法則在正交匹配追蹤算法的基礎(chǔ)上進(jìn)行了改進(jìn)，通過每次選擇多個(gè)原子，提高了重構(gòu)速度和精度。在實(shí)際應(yīng)用中，為了提高重構(gòu)速度和精度，還可以將多種重構(gòu)算法相結(jié)合，或者采用一些加速技術(shù)，如并行計(jì)算、GPU加速等。此外，還可以通過引入先驗(yàn)知識（如圖像的邊緣信息、組織的解剖結(jié)構(gòu)等）來進(jìn)一步提高重構(gòu)質(zhì)量。例如，在腦部磁共振成像中，可以利用腦部的解剖結(jié)構(gòu)先驗(yàn)知識，對重構(gòu)過程進(jìn)行約束，從而得到更清晰、更準(zhǔn)確的腦部圖像。三、壓縮感知磁共振成像的特點(diǎn)（一）加速采集，縮短掃描時(shí)間傳統(tǒng)的磁共振成像通常需要較長的掃描時(shí)間，這不僅會導(dǎo)致患者的不適，還容易產(chǎn)生運(yùn)動偽影，影響圖像的質(zhì)量。而壓縮感知磁共振成像通過采用非自適應(yīng)采樣方式，只采集k空間中的部分采樣點(diǎn)，能夠顯著縮短掃描時(shí)間。例如，在腦部磁共振成像中，采用壓縮感知技術(shù)可以將掃描時(shí)間縮短到原來的1/2到1/10，甚至更短。這對于那些不能長時(shí)間保持靜止的患者（如兒童、老年患者、重癥患者等）來說具有重要的意義，能夠減少運(yùn)動偽影的產(chǎn)生，提高圖像的質(zhì)量。同時(shí)，縮短掃描時(shí)間還可以提高磁共振成像設(shè)備的使用效率，減少患者的等待時(shí)間。（二）減少運(yùn)動偽影運(yùn)動偽影是磁共振成像中常見的問題之一，它是由于患者在掃描過程中的運(yùn)動（如呼吸運(yùn)動、心臟運(yùn)動、頭部運(yùn)動等）導(dǎo)致的。傳統(tǒng)的磁共振成像由于掃描時(shí)間較長，患者更容易產(chǎn)生運(yùn)動，從而導(dǎo)致運(yùn)動偽影的產(chǎn)生。而壓縮感知磁共振成像由于掃描時(shí)間短，能夠減少患者的運(yùn)動時(shí)間，從而減少運(yùn)動偽影的產(chǎn)生。此外，壓縮感知的重構(gòu)算法具有一定的抗噪聲和抗偽影能力，能夠在一定程度上抑制運(yùn)動偽影的影響。例如，在腹部磁共振成像中，由于呼吸運(yùn)動的影響，傳統(tǒng)的磁共振成像容易產(chǎn)生明顯的運(yùn)動偽影，而采用壓縮感知技術(shù)可以在較短的時(shí)間內(nèi)完成掃描，減少呼吸運(yùn)動對圖像的影響，得到更清晰的腹部圖像。（三）提高圖像質(zhì)量壓縮感知磁共振成像在加速采集的同時(shí)，還能夠保證圖像的質(zhì)量。這是因?yàn)閴嚎s感知理論利用了信號的稀疏性，通過非線性重構(gòu)算法能夠從少量的采樣點(diǎn)中準(zhǔn)確地恢復(fù)出原始信號。與傳統(tǒng)的加速成像技術(shù)（如并行成像）相比，壓縮感知MRI具有更高的加速倍數(shù)和更好的圖像質(zhì)量。并行成像技術(shù)是利用多個(gè)接收線圈的空間敏感性來實(shí)現(xiàn)加速采集，但它對線圈的數(shù)量和性能要求較高，并且在加速倍數(shù)較高時(shí)容易產(chǎn)生明顯的偽影。而壓縮感知MRI則不受線圈數(shù)量的限制，能夠在較高的加速倍數(shù)下保持較好的圖像質(zhì)量。例如，在高分辨率磁共振成像中，傳統(tǒng)的成像方式需要采集大量的采樣點(diǎn)，掃描時(shí)間很長，而采用壓縮感知技術(shù)可以在較短的時(shí)間內(nèi)完成高分辨率圖像的采集，并且圖像的質(zhì)量能夠滿足臨床診斷的需求。（四）拓展磁共振成像的應(yīng)用領(lǐng)域壓縮感知磁共振成像的出現(xiàn)為磁共振成像的應(yīng)用領(lǐng)域帶來了新的拓展。例如，在動態(tài)磁共振成像中，傳統(tǒng)的成像方式由于掃描時(shí)間長，無法實(shí)時(shí)捕捉到組織的動態(tài)變化。而壓縮感知MRI可以在很短的時(shí)間內(nèi)完成多期相的圖像采集，實(shí)現(xiàn)動態(tài)磁共振成像，如心臟動態(tài)成像、腹部動態(tài)成像等。這對于研究組織的功能和代謝變化具有重要的意義。此外，壓縮感知MRI還可以用于磁共振血管成像（MagneticResonanceAngiography,MRA）、磁共振彌散成像（DiffusionWeightedImaging,DWI）、磁共振灌注成像（PerfusionWeightedImaging,PWI）等領(lǐng)域，為這些領(lǐng)域的研究和臨床應(yīng)用提供新的技術(shù)手段。例如，在磁共振血管成像中，采用壓縮感知技術(shù)可以在較短的時(shí)間內(nèi)完成血管圖像的采集，減少患者的對比劑用量，同時(shí)提高血管圖像的質(zhì)量和分辨率。（五）對硬件和算法的要求較高雖然壓縮感知磁共振成像具有很多優(yōu)點(diǎn)，但它也存在一些不足之處。首先，壓縮感知MRI對磁共振成像設(shè)備的硬件要求較高。由于需要進(jìn)行非自適應(yīng)采樣，設(shè)備需要具備快速切換梯度場和射頻脈沖的能力，以實(shí)現(xiàn)不同采樣軌跡的采集。此外，設(shè)備的信號采集系統(tǒng)需要具備較高的靈敏度和動態(tài)范圍，以保證采集到的信號質(zhì)量。其次，壓縮感知MRI的重構(gòu)算法計(jì)算復(fù)雜度較高，需要大量的計(jì)算資源和時(shí)間。雖然近年來隨著計(jì)算機(jī)技術(shù)的發(fā)展，GPU加速、并行計(jì)算等技術(shù)的應(yīng)用使得重構(gòu)速度有了很大的提高，但在一些實(shí)時(shí)性要求較高的應(yīng)用中，重構(gòu)速度仍然是一個(gè)需要解決的問題。此外，重構(gòu)算法的選擇和參數(shù)設(shè)置對圖像的質(zhì)量也有很大的影響，需要根據(jù)具體的應(yīng)用場景進(jìn)行優(yōu)化和調(diào)整。（六）存在一定的局限性壓縮感知磁共振成像的應(yīng)用也存在一定的局限性。首先，信號稀疏性是壓縮感知理論的前提條件，如果信號不具有稀疏性，或者稀疏性較差，那么壓縮感知MRI的重構(gòu)質(zhì)量就會受到影響。例如，一些復(fù)雜的組織或病變（如腫瘤組織、炎癥組織等）的磁共振信號稀疏性可能較差，采用壓縮感知技術(shù)進(jìn)行成像時(shí)，可能會導(dǎo)致圖像的模糊或偽影。其次，壓縮感知MRI的加速倍數(shù)也受到一定的限制。雖然理論上可以實(shí)現(xiàn)很高的加速倍數(shù)，但在實(shí)際應(yīng)用中，過高的加速倍數(shù)會導(dǎo)致圖像質(zhì)量的下降，因此需要在加速倍數(shù)和圖像質(zhì)量之間進(jìn)行權(quán)衡。此外，壓縮感知MRI的圖像質(zhì)量還受到采樣方式、測量矩陣設(shè)計(jì)、重構(gòu)算法等多種因素的影響，需要進(jìn)行綜合優(yōu)化和調(diào)整。四、壓縮感知磁共振成像的臨床應(yīng)用前景（一）在腦部成像中的應(yīng)用在腦部磁共振成像中，壓縮感知技術(shù)具有廣闊的應(yīng)用前景。例如，在腦部功能磁共振成像（fMRI）中，傳統(tǒng)的成像方式需要采集大量的時(shí)間點(diǎn)數(shù)據(jù)，掃描時(shí)間很長，而采用壓縮感知技術(shù)可以在較短的時(shí)間內(nèi)完成數(shù)據(jù)采集，提高時(shí)間分辨率，從而更好地捕捉腦部的功能活動變化。此外，在腦部彌散張量成像（DTI）中，壓縮感知技術(shù)可以減少彌散加權(quán)圖像的采集數(shù)量，縮短掃描時(shí)間，同時(shí)保證彌散張量參數(shù)的準(zhǔn)確性。這對于研究腦部的白質(zhì)纖維束結(jié)構(gòu)和功能具有重要的意義。在腦部腫瘤成像中，壓縮感知技術(shù)可以在較短的時(shí)間內(nèi)完成多參數(shù)磁共振成像（如T1加權(quán)成像、T2加權(quán)成像、彌散加權(quán)成像、灌注加權(quán)成像等），為腫瘤的診斷、分期和治療評估提供更全面、更準(zhǔn)確的信息。（二）在心血管成像中的應(yīng)用心血管磁共振成像（CMR）是一種無創(chuàng)性的心血管疾病診斷方法，但傳統(tǒng)的CMR掃描時(shí)間較長，容易受到心臟運(yùn)動和呼吸運(yùn)動的影響。壓縮感知技術(shù)在心血管成像中的應(yīng)用可以顯著縮短掃描時(shí)間，減少運(yùn)動偽影的產(chǎn)生。例如，在心臟電影成像中，采用壓縮感知技術(shù)可以在一個(gè)心動周期內(nèi)采集更多的圖像幀，提高時(shí)間分辨率，從而更準(zhǔn)確地評估心臟的功能和運(yùn)動情況。在心肌灌注成像中，壓縮感知技術(shù)可以減少對比劑的用量，同時(shí)提高灌注圖像的質(zhì)量和時(shí)間分辨率，更好地顯示心肌的灌注情況，為冠心病的診斷和治療評估提供更可靠的依據(jù)。（三）在腹部成像中的應(yīng)用腹部磁共振成像由于受到呼吸運(yùn)動的影響，傳統(tǒng)的成像方式需要患者配合屏氣，這對于一些不能屏氣的患者來說是非常困難的。而