Sairam Geethanath, MSHyeon-Man Baek, PhDSandeep K. Ganji, BEYao Ding, MSElizabeth A. Maher, MDRobert D. Sims, MDChangho Choi, PhD,,2dMatthew A. Lewis, PhD,,Vikram D. Kodibagkar, PhD,,

壓縮感知可以加速醫(yī)用1H核磁代謝成像過程

Sairam Geethanath, MS1Hyeon-Man Baek, PhD2aSandeep K. Ganji, BE2bYao Ding, MS2bElizabeth A. Maher, MD2cRobert D. Sims, MD2aChangho Choi, PhD2a,2b,2dMatthew A. Lewis, PhD1,2a,2bVikram D. Kodibagkar, PhD1,2a,2b

1. 德州大學(xué)阿靈頓分校與德州大學(xué)西南醫(yī)學(xué)院聯(lián)合生物醫(yī)學(xué)工程研究生計劃項目，美國 德克薩斯；2. 德州大學(xué)西南醫(yī)學(xué)院a.放射學(xué)系; b. 放射學(xué)研究生計劃項目; c. 內(nèi)科與神經(jīng)學(xué)系; d. 高級成像研究中心，美國德克薩斯

目的回顧性評測使用壓縮感知（compressed sensing）技術(shù)進行不同程度的核磁成像過程（MR）加速，對所得結(jié)果保真度的影響。材料和方法本文中的實驗手段均由所在院校倫理審查委員會（Institutional Review Board, IRB）審查通過，采集MR成像數(shù)據(jù)前獲得了患者書面知情同意。本研究遵守《健康保險攜帶和責(zé)任法案》 ( Health Insurance Portability and Accountability Act，HIPAA）?；仡櫺詨嚎s感知技術(shù)被用于10名實驗對象的核磁共振成像過程中，并獲取了以下體素數(shù)據(jù)：600個來自6名健康對象的腦部；163個來自兩名腦瘤患者的腦部；36個來自于兩名前列腺癌患者。研究人員在加速常數(shù)分別為2，3，4，5，10的條件下進行核磁信號的重建，并通過均方根誤差（RMSE）法，代謝物圖譜（包括膽堿，肌酸，N乙酰天冬氨酸[NAA]，檸檬酸）以及統(tǒng)計學(xué)分析（包括基于單個體素的配對t檢驗，代謝圖譜上不同變量的單因素誤差分析，以及加速重建后與原始變量的比值）進行評估。結(jié)果使用不超過10的加速常數(shù)，重建核磁信號的保真度很高，RMSE很低（＜0.05）。在加速常數(shù)不超過5的情況下，重建信號的平均生化信號強度和生化熱點位置與原信號非常接近，二者之間沒有統(tǒng)計學(xué)意義上的顯著差異。在加速常數(shù)不超過5的重建信號中，膽堿/NAA以及（膽堿+肌酸）/檸檬酸的比值與原始信號并無顯著差異，在部分加速常數(shù)為10的重建信號中上述條件亦成立。結(jié)論本文闡明了一種最高能夠?qū)⒉蓸訒r間降低80%，同時保持信號信息量的醫(yī)用1H核磁生化成像過程。

核磁成像（MR）技術(shù)正在越來越廣泛地被應(yīng)用于臨床 診療過程中[1-3]。相關(guān)研究表明，惡性前列腺和腦部腫瘤均會表達高量的膽堿水平[4]。另有研究表明，腦部腫瘤會表達較低水平的N乙酰天冬氨酸和肌酸信號，因此可以使用膽堿 /NAA 比例作為腫瘤的生化指標(biāo)[4-7]。與前者類似，對于前列腺腫瘤的診斷也可以使用（膽堿 +肌酸）/檸檬酸比例作為生化指標(biāo)。研究表明，在惡性前列腺腫瘤中，檸檬酸水平會降低，而膽堿水平會升高[4]。

多體素（二維/三維）核磁共振成像相對于單體素核磁共振技術(shù)的優(yōu)勢已被廣泛認知[7-8]。使用該技術(shù)，醫(yī)護人員可以獲取更優(yōu)秀的腫瘤定位信息，并可以觀察到形變病變區(qū)域以外的異常表象。但是多體素核磁技術(shù)的主要技術(shù)缺陷——增長的采樣時間——是限制其臨床應(yīng)用的主要因素[7-9]。壓縮感知技術(shù)[10-11]提供了一種創(chuàng)新性的利用在適當(dāng)變換域（如小波變換域）下信號密度來對k-空間進行低密度取樣的技術(shù)，利用該技術(shù)可以降低重建過程對核磁信號數(shù)據(jù)點的需求數(shù)量。此技術(shù)已經(jīng)在腦部核磁共振成像、核磁共振血管造影[10]、徑向成像[11]、心臟電影成像[12]，以及其他動態(tài)對比高亮核磁共振成像手段[13]中得到廣泛應(yīng)用。壓縮感知技術(shù)應(yīng)用于核磁成像中，因為核磁信號數(shù)據(jù)在小波變換域的多個頻率和空間維度上均比較稀疏，變化也較小。過去的研究中，基于小波的分析手段，包括信號量化波形修正和去噪技術(shù)，均被用于核磁共振成像[14-17]。研究人員通過使用小波技術(shù)手段，利用核磁數(shù)據(jù)在此變換域中頻率和空間維度上的信號稀疏，已經(jīng)實現(xiàn)了快速核磁共振成像[18-19]。本研究中，研究人員回顧性評測了壓縮采樣技術(shù)重建信號的保真程度。

1 材料與方法

本研究回顧性地展示了6位健康志愿者，2位腦瘤（惡性神經(jīng)膠質(zhì)瘤）患者及 2 位前列腺癌患者的二維1H 核磁成像數(shù)據(jù)，各項協(xié)議均經(jīng)大學(xué)的倫理審查委員會批準(zhǔn)，并獲得了志愿者及患者的書面知情同意。本研究遵守 HIPAA 法案。所得磁共振參數(shù)詳見表1。

1.1 核磁共振成像數(shù)據(jù)及采樣

原始數(shù)據(jù)歸一化至 0～1，在 k 空間中沿相位維度對該數(shù)據(jù)進行低密度采樣，采樣模式由不同的加速因子（2，3，4，5，10）以及采樣密度決定。采樣模式為二維平面分布函數(shù)加權(quán)的隨機采樣，使得k空間中心附近采樣密度高于 邊 緣 部位[11]。二 維 平 面分布函數(shù) 的 中 心抽樣范圍取決于加速因子的倒數(shù)（加速因子越大，充分抽樣區(qū)域的面積越小）。

1.2 信號重建

將抽樣 k空間全部設(shè)定為 0，并經(jīng)過逆式傅立葉轉(zhuǎn)換來獲取對核磁共振數(shù)據(jù)的預(yù)估。核磁共振信號重建由對下述代價函數(shù)進行凸優(yōu)化完成（Lustig 等人[11]）。

該 信 號 重 建 過 程 通 過 Matlab（MathWorks，Natick，Mass）提供的非線性共軛梯度算法完成，其中m代表核磁共振數(shù)據(jù)，F(xiàn)u代表傅立葉轉(zhuǎn)換運算符，y代表已測 k空間數(shù)據(jù)，W代表微波轉(zhuǎn)換運算符，TV 代表總變差運算符，‖1和 ‖2為 L1 和 L2 各自的基準(zhǔn)運算符 ；λL1和 λTV為基于L1和總變差各自的正規(guī)化參數(shù)。本研究中，研究人員使用Daubechies 小波變換技術(shù)[20]對 ky維度上所有像素點的 kx-t數(shù)據(jù)矩陣進行編譯。λL1和 λTV的正規(guī)化值通過實驗分別定為 0.001 和 0.005。通過對所獲數(shù)據(jù)進行 8 次迭代操作，實驗人員得以獲取并比對代價函數(shù)的具體數(shù)值。

1.3 錯誤度量標(biāo)準(zhǔn)

信號重建誤差由下列公式計算所得的均方根誤差（RMSE）決定 ：

其中N為一幀磁共振成像數(shù)據(jù)組中數(shù)據(jù)點的總和，y為完整的 k 空間的重建數(shù)據(jù)，y'為低密度取樣 k 空間的重建數(shù)據(jù)。因為所用數(shù)據(jù)已經(jīng)經(jīng)過正態(tài)化至 0～1，所以上述方根誤差等于通用的正態(tài)方根誤差。

1.4 后期處理

磁共振成像數(shù)據(jù)組受基于 Java 的磁共振用戶界面中的下列最少處理步驟的影響 (jMRUI; http://www.mrui.uab.es/mrui/ mrui_Overview.shtml)[21]： 用 衍 射 控 象 法 移 除 現(xiàn) 有 的 平 截 偽影， 基線校正， 時域的 Hankel-Lanczos 奇異值分解過濾的殘余水峰和脂肪峰， 自動無序全相位校正， 使用量子估算算法（QUEST）[22]對相位核磁共振數(shù)據(jù)的實部進行分析，用以獲取生化物質(zhì)分布圖，在該圖中信號的強度代表 QUEST 算法結(jié)果中的譜圖信號分布。本研究中，對于腦部磁共振成像的案例（健康志愿者與癌癥患者），研究人員只分析了大腦內(nèi)部區(qū)域的數(shù)據(jù)。進一步的分析中，關(guān)于健康大腦、腦瘤和前列腺癌患者的體素總數(shù)分別為600，163，36。計算相關(guān)比例時，原始數(shù)據(jù)或壓縮傳感重建數(shù)據(jù)中的分母為0的像素點的數(shù)值

（表示缺乏適合量子估值的適當(dāng)應(yīng)用手段）設(shè)定為0。

表1 磁共振數(shù)據(jù)采集參數(shù)

1.5 統(tǒng)計學(xué)分析

每組信號重建獲取的生化組分分布圖以及相對應(yīng)的生化信號比值（膽堿 /NAA，健康大腦和腦瘤患者 ；膽堿 +肌酸 /檸檬酸，前列腺癌癥患者）均通過配對t檢驗進行統(tǒng)計學(xué)分析（Excel； Microsoft，Redmond，Wash），并被用來與原始數(shù)據(jù)（加速因子為1）進行比對。此外，研究人員還通過單因子重復(fù)分析測試與 Bonferroni多組分比對測試來進一步評測各種加速比例對核磁信號重建質(zhì)量的影響（GraphPad Prism ；GraphPad Software，La Jolla，Calif）。通過 Bonferroni檢驗獲取的t 值數(shù)據(jù)被轉(zhuǎn)化為P 值（Excel；Microsoft）。 如果P ＜ 0.05，則認為該組比較中存在明顯差異。

2 結(jié)果

圖1和圖2展示了典型健康大腦磁共振成像數(shù)據(jù)組的重建結(jié)果。圖1展示了對解剖圖中同一部位使用不同加速因子進行成像的核磁共振圖譜。壓縮感知的重建核磁信號與原始信號顯示出了很高的相似度。用于評價成像質(zhì)量的兩個體素的位置由綠框和紅框在解剖圖中標(biāo)出（圖 2）。覆蓋腦室的體素（綠框）數(shù)據(jù)與另一個體素（紅框）相比，NAA，肌酸和膽堿（圖 2右排）的濃度均有所降低，此結(jié)果與此前預(yù)期相符。圖2中在不同加速因子下，重建的成像結(jié)果與原始數(shù)據(jù)相似。壓縮感知重建數(shù)據(jù)可以保持原始數(shù)據(jù)的譜線形狀，而且在高加速因子下，譜線更加平滑。圖3展示了 NAA，肌酸和膽堿（譜峰覆蓋面積）的代謝分布圖譜和膽堿 /NAA 比值與加速因子之間的關(guān)系。原始和重建信號中的高亮區(qū)域出現(xiàn)在相似位置。而且在加速因子為10的時候，生化信號的強度比原始信號更強。

圖1 1倍，2倍，5倍和10倍加速度下的典型大腦磁共振光譜成像數(shù)據(jù)組的數(shù)據(jù)網(wǎng)格的重建。中心部位：核磁共振成像顯示值得注意的成像區(qū)域（黃框內(nèi)）。2個位置（紅框和綠框中的體素）的相關(guān)細節(jié)詳見圖2。

圖2 在兩個體素（圖1中紅色和綠色的框體范圍）中采集到的加速因子為1，2，5和10的核磁共振光譜成像數(shù)據(jù)。y軸上下限為-0.25～1。Cho=膽堿，Cr=肌酸。

圖3 加速因子為1，2，5和10的實驗中健康大腦的NAA，肌酸，膽堿和膽堿/NAA指標(biāo)的分布對照圖。圖1解剖圖中黃框表明的區(qū)域為圖譜上描繪的要注意的磁共振光譜成像區(qū)域。a.u.=arbitrary unit（任意單元）。

正常組織和惡性腫瘤組織的2個體素的原始數(shù)據(jù)和壓縮感知重建信號的腦瘤核磁共振成像數(shù)據(jù)組詳見圖4。加速因子在 10以內(nèi)的重建信號如實地保留了用來辨別正常腦組織與癌變組織的生化組分分布特點。圖4為前列腺癌核磁共振成像數(shù)據(jù)組中代表正常組織與惡性腫瘤組織的2個體素的重建信號。壓縮感知重建信號保持了正常前列腺與前列腺腫瘤病例的代謝物分布的特征。根據(jù)這些信號重建生成的生化組分分布圖譜與加速因子之間的相關(guān)關(guān)系詳見圖 5a 和圖 5b。對于腦瘤數(shù)據(jù)組，加速因子為 2 和 5 時，代謝產(chǎn)物分布圖和膽堿 /NAA 指標(biāo)圖譜與原始信號擁有相同的特征，而加速因子為 10的壓縮感知重建信號則相對于原始數(shù)據(jù)有較大的偏差。對于前列腺癌數(shù)據(jù)組，壓縮感知重建信號中的生化組分分布圖的信號強度有所降低，但（膽堿+肌酸）/檸檬酸指標(biāo)圖譜則保留了原始數(shù)據(jù)的所有特征，尤其是高亮區(qū)域的位置。

圖4 加速因子為1，2，5和10的實驗中下腦瘤（左）和前列腺癌（右）核磁共振光譜成像數(shù)據(jù)。在所有情況下，左列展示的是健康組織的數(shù)據(jù)，而右列展示的是腫瘤區(qū)域的數(shù)據(jù)。y軸上下限為-0.25～1。Cho=膽堿，Cit=檸檬酸，Cr=全部肌酸的甲基信號，Cr2=全部肌酸的亞甲基信號。

圖5 a 加速因子為1，2，5和10的實驗中腦瘤（NAA，肌酸，膽堿）和核磁共振光譜成像生化組分分布比較和比例圖譜。

圖5 b 前列腺癌（檸檬酸，肌酸，膽堿）比例圖譜描繪了前列腺癌數(shù)據(jù)中腦瘤數(shù)據(jù)和（膽堿+肌酸）／檸檬酸的膽堿/NAA指標(biāo)（CNI）。解剖圖中的黃框指出感興趣區(qū)域。a.u.=arbitrary unit（任意單元）。

加速因子 1，2，5，10 下的每項生化組分含量的平均值±標(biāo)準(zhǔn)差的計算結(jié)果詳見表2。在加速因子為2和 5時，壓縮感知重建信號的生化組分含量水平與原始信號比較一致，而加速因子為 10時，9種生化組分中有 6種較原始數(shù)據(jù)有著統(tǒng)計學(xué)上的顯著差異。不同生化組分之間的比例指標(biāo)分布圖也與原始信號相似，這說明加速因子5以內(nèi)的健康大腦/腦瘤數(shù)據(jù)和加速因子 10以內(nèi)的前列腺癌的數(shù)據(jù)的關(guān)鍵生物指標(biāo)值與原始信號相比沒有顯著差異。通過 Bonferroni檢驗分析單因子差異實驗數(shù)據(jù)得出的P值詳見表3。這些數(shù)值表明加速因子 1，2 和 5 之間聯(lián)系緊密，而加速因子 10 和其它加速因子之間的關(guān)系顯然沒有那么緊密，從而導(dǎo)致P值的下降。在這 5 種加速度因子下的 3 種數(shù)據(jù)類型的 RMSE 量化信號重建誤差詳見圖 6。與預(yù)期相同，RMSE 值隨著加速因子的增大而增大，但是加速因子不超過 10的情況下，RMSE 維持在 0.05 以下。

3 討論

本研究的目的是尋找一種可以在不丟失信息的情況下加速獲取1H 磁共振圖像的途徑，以便將1H 核磁共振成像廣泛應(yīng)用于臨床。通過使用壓縮感知信號重建，我們證明了1H 核磁共振成像的采集時間可能節(jié)省 80% 甚至更多，且根據(jù)臨床標(biāo)準(zhǔn)評估，幾乎不損失任何重要信息。

表2 健康大腦，腦瘤及前列腺癌代謝強度與體素數(shù)據(jù)比值作為加速函數(shù)

表3 Bonferroni對比試驗后單因子差異分析產(chǎn)生的P值

圖6 加速倍數(shù)為2，3，4，5和10的實驗中得出的核磁共振光譜成像數(shù)據(jù)RMSE值

重建光譜發(fā)出的噪音比原始光譜的小。這要歸功于重建過程中的小波靜噪功能和總變差因子的濾波效果。如上文所述，本研究中采用的信號重建方法能夠產(chǎn)生相似于原始信號的生化組分分布圖譜。核磁共振成像獲取的生化組分分布信息對診斷和預(yù)后過程很重要，而本研究中的信號重建過程能夠如實地保留這些信息。在原始信號中（加速因子 =1），即便是健康大腦的成像結(jié)果中，生化信號中均可觀察到過高的標(biāo)準(zhǔn)差，而這一現(xiàn)象是源于頭部線圈覆蓋了整個頭部，包括腦髓液部分（圖 2～3）。除了對這些圖譜的直接觀測，配對t檢驗和單因子差異分析也可以被用來比較特定部位的原始信號和不同加速因子下壓縮感知重建信號。當(dāng)我們對加速度和原型之間的保真度和微小差異而非顯著性差異進行測試的時候，將P值閾值設(shè)置于 0.05 時，比其更低的數(shù)據(jù)，如 0.01 更為嚴(yán)格，因為P 值小于 0.05 表明兩種信號之間存在顯著性差異，也就意味著信號重建算法的失敗。本文的結(jié)果同樣表明，對于各種生化指標(biāo)，例如膽堿／NAA和膽堿／檸檬酸，信號重建并不導(dǎo)致其發(fā)生變化。因此，信號重建也同樣維持了此類衍生參數(shù)的保真度。統(tǒng)計分析表明，加速因子為5的壓縮感知重建信號與原始信號相似，而且沒有統(tǒng)計學(xué)上的顯著差異。然而在加速因子為 10 時，重建信號的誤差及生物指標(biāo)值的差異有所增加。另一方面，加速因子 10 以內(nèi)的所有重建數(shù)據(jù)，其 RMSE 數(shù)值均在 0.05 以內(nèi)。由于數(shù)據(jù)在處理之前已經(jīng)正態(tài)化，這說明重建信號與原始信號之間的誤差不超過5%。將核磁共振成像時域數(shù)據(jù)轉(zhuǎn)換成生化組分分布是一個復(fù)雜的過程，涉及許多用戶自定義的步驟，這些步驟會潛在地影響最終生化組分分布的圖譜。研究人員可以通過比較原始數(shù)據(jù)和最少后處理的壓縮感知重建信號數(shù)據(jù)來客觀評估信號重建的保真度，但較高加速因子下的重建信號圖譜中的差異可能是自動化整體位相校正中的差異造成的。大腦部位的生化組分分布圖信號強度隨加速因子的提高而增強，前列腺部位則反之?；诒O(jiān)控參數(shù)，網(wǎng)格尺寸以及k空間信噪比等參數(shù)的不同，上述兩種趨勢都可能出現(xiàn)，因此上述結(jié)果并無參考意義。對于給定的監(jiān)控參數(shù)和加速因子并采用最優(yōu)化程序的情況下，通過 RMSE 值來評價的試驗結(jié)果誤差隨著加速因子的增加而增加。因此，對生化組分含量的估算數(shù)據(jù)將隨加速因子的增加而逐漸偏離原始數(shù)據(jù)。從實驗方法的實用性角度出發(fā)，文中所示方法測算所得的數(shù)據(jù)誤差與加速因子成單增或單減關(guān)系（不管指標(biāo)誤差在越來越多的不充分采樣中是增加還是減少），RMSE 值在最大為 10的加速因子下也被控制得很小。

其他幾種快速化學(xué)位移成像和光譜學(xué)成像方法也已經(jīng)發(fā) 展 成熟并付諸在 理 論[23]和 實 踐[24]上進行了 相 關(guān) 對 比。上 述 技術(shù)中，回波平面 光 譜成像[25-26]是一 種在臨床上[27]尚未被充分利用過的強大的光譜成像技術(shù)。本方法通過在一個頻域和一個空間維度上同時進行回聲平面采樣來加速成像過程，從而達到相當(dāng)于那個空間維度中像素點數(shù)的加速因子。例如，使用回波平面采樣技術(shù)采集 1 個 16×16 的二維矩陣，所用時間會是傳統(tǒng)化學(xué)位移成像手段的 1/16。然而，提高采樣速度或許會以信噪比為代價，而且所得的成像結(jié)果也會受 Nyquist偽影[28]的影響。本文中使用的壓縮感知技術(shù)中有選擇地忽略了部分特定相位的數(shù)據(jù)，該步驟在二維和三維化學(xué)位移成像序列或回波平面成像技術(shù)中可以比較容易地完成。在進一步的研究中，可以通過設(shè)計快速成像序列，并將其與壓縮感知技術(shù)進行聯(lián)用，以獲取更快的成像速度。最近的研究表明，在更多維度上使用壓縮感知技術(shù)可以更有效地利用信號分布的稀疏性，并因此提高重建信號的質(zhì)量[11]。對于核磁共振光譜成像數(shù)據(jù)，頻域維度上的數(shù)據(jù)稀疏程度比空間維度上的更高。為了獲取頻域和時域上的隨機采樣數(shù)據(jù)，研究人員需要設(shè)計一組包含快速變換梯度的復(fù)雜激活信號序列。對超極化13C 核磁共振壓縮感知數(shù)據(jù)的重建過程中已經(jīng)用到此種技術(shù)。然而，通過在相位編譯維度上進行壓縮感知而減少的采樣時間，要遠遠超出僅在臨時維度上進行壓縮感知能夠達到的程度。而且與此同時，除了回聲平面采樣技術(shù)，使用其他采樣手段的時候，信號收集加速比例和采樣壓縮比例之間并無直接的相關(guān)關(guān)系。在本研究中，我們致力于加快常規(guī)臨床核磁共振成像的信號收集過程，而該結(jié)論的應(yīng)用范圍并不僅限于上述討論的范疇。

在光譜成像領(lǐng)域，所有新型信號重建技術(shù)的一個重要考量是該技術(shù)對點擴散函數(shù)的影響，以及因此導(dǎo)致的相鄰空間體素信息污染的問題。Lustig 等人[11]詳細討論過點擴散函數(shù)和用于核磁共振成像壓縮感知的變換點擴散函數(shù)的性質(zhì)。上述研究中涉及的使用三維傅立葉變換的案例中，點分布函數(shù)與本文中的二維核磁譜圖非常類似。在本文中，該二維核磁譜圖在三個維度上進行取樣：兩個相位編譯的空間維度和一個采用完全模擬信號的臨時維度（無延遲/回聲）[11]。壓縮感知技術(shù)的優(yōu)勢在于，該技術(shù)通過采用隨機采樣的方法，避免了偽影的干擾（在回聲平面成像以及回聲平面頻域成像中常見），在整個成像區(qū)域中形成類似于噪音的信號。與之相反，偽影則是點分布函數(shù)的擴散導(dǎo)致的體素之間的信息干擾。

當(dāng)前研究的局限在于它是在此前獲得的臨床數(shù)據(jù)的基礎(chǔ)上回顧性完成的，而回顧性信號重建是決定我們研究方法可行性的重要的第一步，接下來的必要步驟是對健康的志愿者進行實時的信號壓縮感知，重建，并對比不同加速因子下的結(jié)果。其他可能的改進包括將該技術(shù)擴展至三維核磁共振光譜成像，以及改進信號重建程序以更好地利用空間和頻域的信號稀疏性，特別是三維核磁共振光譜成像過程。

總而言之，本研究表明，在加速因子不超過 5 的情況下，壓縮感知重建信號對原始信號的損失是微不足道的，同時還可節(jié)省最高可達 80% 的采樣時間。這將能夠促進更多核磁共振光譜成像手段應(yīng)用于臨床。本文中設(shè)計的加速技術(shù)也可以用于其他組織的成像，如乳腺。此外，由于減少了每輪采樣所用的時間，本技術(shù)也可以用于對更難以檢測的生化組分，如甘氨酸。

出版前本文公布情況：

10.1148/radiol.11111098 Content code: MR

[1] Daly PF,Cohen JS.Magnetic resonance spectroscopy of tumors and potential in vivo clinical applications:a review[J].Cancer Res,1989,49(4):770-779.

[2] Gillies RJ,Bhujwalla ZM,Evelhoch J,et al.Applications of magnetic resonance in model systems: tumor biology and physiology[J].Neoplasia,2000,2(1-2):139-151.

[3] Pirzkall A,Li XJ,Oh JM,et al.3D MRSI for resected highgrade gliomas before RT:tumor extent according to metabolic activity in relation to MRI[J].Int J Radiat Oncol Biol Phys, 2004,59(1):126-137.

[4] Kurhanewicz J,Vigneron DB,Nelson SJ.Three-dimensional magnetic resonance spectroscopic imaging of brain and prostate cancer[J].Neoplasia,2000,2(1-2):166-189.

[5] Nelson SJ,Vigneron DB,Dillon WP.Serialevaluation of patients with brain tumors using volume MRI and 3D1H MRSI[J]. NMR Biomed,1999,12(3):123-138.

[6] Nelson SJ.Multivoxel magnetic resonance spectroscopy of brain tumors[J].Mol Cancer Ther,2003,2(5):497-507.

[7] Howe FA,Opstad KS.1H MR spectroscopy of brain tumours and masses[J].NMR Biomed,2003,16(3):123-131.

[8] Ebel A,Soher BJ,Maudsley AA.Assessment of 3D proton MR echo-planar spectroscopic imaging using automated spectral analysis[J].Magn Reson Med,2001,46(6):1072-1078.

[9] Zhu H,Barker PB.MR spectroscopy and spectroscopic imaging of the brain[J].Methods Mol Biol,2011,711:203-226.

[10] Candes EJ,Romberg J,Tao T.Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information[J].IEEE Trans Inf Theory,2006,52(2):489-509.

[11] Lustig M,Donoho D,Pauly JM.Sparse MRI: the application of compressed sensing for rapid MR imaging[J].Magn Reson Med,2007,58(6):1182-1195.

[12] Jung H,Sung K,Nayak KS,et al.k-t FOCUSS: a general compressed sensing framework for high resolution dynamic MRI[J].Magn Reson Med,2009,61(1):103-116.

[13] Jim J,Tao L.Dynamic MRI with compressed sensing imaging using temporal correlations[A].In: 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro[C]. Paris: Institute of Electrical and Electronics Engineers,2008: 1613-1616.

[14] Serrai H,Senhadji L,de Certaines JD,et al.Time-domain quanification of amplitude, chemical shift,apparent relaxation time T2,and phase by wavelet-transform analysis: application to biomedical magnetic resonance spectroscopy[J].J Magn Reson, 1997,124(1):20-34.

[15] Serrai H,Nadal-Desbarats L,Poptani H,et al.Lactate editing and lipid suppression by continuous wavelet transform analysis: application to simulated and (1)H MRS brain tumor timedomain data[J].Magn Reson Med,2000,43(5):649-656.

[16] Cancino-De-Greiff HF,Ramos-Garcia R,Lorenzo-Ginori JV.Signal de-noising in magnetic resonance spectroscopyusing wavelet transforms[J].Concepts Magn Reson,2002,14(6): 388-401.

[17] Zarei S,Alireazie J,Babyn P,et al.MRS feature extraction: time frequency and wavelet analysis[A].In:The 2nd International Conference on Bioinformatics and Biomedical Engineering[C].S hanghai,China:ICBBE,2008;1863-1866.

[18] Fu Y,Serrai H.Fast magnetic resonance spectroscopic imaging (MRSI) using wavelet encoding and parallel imaging:in vitro results[J].J Magn Reson,2011,211(1):45-51.

[19] Hu S,Lustig M,Chen AP,et al.Compressed sensing for resolution enhancement of hyperpolarized 13C flyback 3D-MRSI[J]. J Magn Reson,2008,192(2):258-264.

[20] Antonini M,Barlaud M,Mathieu P,et al.Image coding using wavelet transform[J].IEEE Trans Image Process,1992, 1(2): 205-220.

[21] Naressi A,Couturier C,Castang I,et al.Java-based graphical user interface for MRUI,a software package for quantitation of in vivo/medicalmagnetic resonance spectroscopy signals[J].Comput Biol Med,2001,31(4):269-286.

[22] Ratiney H,Coenradie Y,Cavassila S,et al.Time-domain quantitation of1H short echo-time signals:background accommodation[J].MAGMA,2004,16(6):284-296.

[23] Pohmann R,von Kienlin M,Haase A.Theoretical evaluation and comparison of fast chemical shift imaging methods[J].J Magn Reson,1997,129(2):145-160.

[24] Zierhut ML,Ozturk-Isik E,Chen AP,et al.(1)H spectroscopic imaging of human brain at 3 Tesla:comparison of fast threedimensional magnetic resonance spectroscopic imaging techniques[J].J Magn Reson Imaging,2009,30(3):473-480.

[25] Mansfield P.Spatial mapping of the chemical shift in NMR[J]. Magn Reson Med,1984,1(3):370-386.

[26] Posse S,DeCarli C,Le Bihan D.Threedimensional echo-planar MR spectroscopic imaging at short echo times in the human brain[J].Radiology,1994,192(3):733-738.

[27] Mulkern RV,Panych LP.Echo planar spectroscopic imaging[J]. Concepts Magn Reson,2001,13(4):213-237.

[28] Du W,Du YP,Fan X,et al.Reduction of spectral ghost artifacts in high-resolution echo-planar spectroscopic imaging of water and fat resonances[J].Magn Reson Med,2003,49(6):1113-1120.

[29] Choi C,Ganji SK,Deberardinis RJ,et al.Measurement of glycine in the human brain in vivo by1H-MRS at 3 T: application in brain tumors[J].Magn Reson Med,2011,66(3):609-618.

10.3969/j.issn.1674-1633.2012.07.001

1674-1633(2012)07-0001-08

本文英文原版出自Radiology雜志2012年第262卷3月刊第985～994頁。翻譯及轉(zhuǎn)載均經(jīng)過北美放射學(xué)會許可。北美放射學(xué)會對在翻譯過程中出現(xiàn)的譯文不準(zhǔn)確現(xiàn)象概不負責(zé)。

2011-05-27

修回日期：2011-08-08

錄用日期：2011-09-12

定稿日期：2011-10-18

美國國立衛(wèi)生研究院基金項目支持（NCRR UL1RR024982， NCIR21CA132096-01A1）。

聯(lián)系人：V.D.K. (e-mail: vikram.kodibagkar@asu.edu)

文章內(nèi)容不代表美國聯(lián)邦政府與資助機構(gòu)的立場。

?RSNA 2012