磁共振線圈合并系數(shù)計算方法、磁共振成像方法及其裝置的制造方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及磁共振成像領(lǐng)域，特別涉及一種磁共振線圈合并系數(shù)計算方法、磁共 振成像方法及其裝置。
【背景技術(shù)】
[0002] 在磁共振成像技術(shù)中，成像的速度是衡量成像方法的一個很重要標準。限制成像 速度的很重要因素是數(shù)據(jù)采集，以及k空間填充。一般的數(shù)據(jù)采集方式要采滿k空間數(shù)據(jù)， 然后才能進行重建得到圖像。磁共振并行采集重建技術(shù)，是利用線圈重組合并的方式，對欠 采樣的數(shù)據(jù)進行填補，利用填補完整的k空間數(shù)據(jù)進行重建。利用這樣的方式，可以根據(jù)需 求，只采集一部分k空間數(shù)據(jù)，不必采滿整個k空間。因此這樣的方法可以大大加快成像的 速度。
[0003] 比較常用的并行重建方法之一是GRAPPA。傳統(tǒng)的GRAPPA算法如圖1所示，黑色 實點代表為實際采集的k空間數(shù)據(jù)；白色空點為欠采樣需要填補的數(shù)據(jù)；灰色實點表示為 了計算線圈合并參數(shù)，而適量全采的校準數(shù)據(jù)。GRAPPA算法認為，圖中任意一個空心點可 以表示為周圍黑色實點的線性疊加，相當于對多個線圈的數(shù)據(jù)進行了合并。而線圈合并系 數(shù)Ii ij(第i個線圈，第j個位置，如圖1)可以通過黑色的實點擬合灰色點來確定。線圈合 并系數(shù)確定后，其他白色空心點即可根據(jù)求得的合并系數(shù)以及黑色實點計算填補得到，從 而重建得到完整的k空間數(shù)據(jù)。
[0004] 線圈合并系數(shù)又可稱為卷積核，在傳統(tǒng)方法中，卷積核的計算方向，只加入了相位 編碼方向，以及通道方向。為了優(yōu)化效果，近年來，很多方法引入了其他方向如頻率編碼方 向，以及k-t中時間方向等。這些維度的引入，使卷積核的信息增加，但是隨著過多的引入， 信息會出現(xiàn)冗余，而且由于噪聲的影響，使得計算的系數(shù)仍然受到影響；而卷積核過大，使 得要擬合的系數(shù)增多，也帶來計算的不穩(wěn)定性。而卷積核過小，則數(shù)據(jù)信息不足，計算過程 不夠準確。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明要解決的問題是提供一種磁共振并行成像線圈合并系數(shù)計算方法、磁共振 并行成像方法及其裝置，解決在卷積核選取過程中的由于卷積核選取過大造成信息量增 加，信息冗余，并且噪聲過多，擬合系數(shù)增多，計算不穩(wěn)定的影響，以及在卷積核選取過小 時，計算數(shù)據(jù)不足，結(jié)果不準確的問題。
[0006] 為了實現(xiàn)上述目的，本發(fā)明提供了一種磁共振并行成像線圈合并系數(shù)計算方法， 包括如下步驟：
[0007] a)采集k空間數(shù)據(jù)，所述k空間數(shù)據(jù)包括采集數(shù)據(jù)和校準數(shù)據(jù)；
[0008] b)對所述采集數(shù)據(jù)進行數(shù)據(jù)提煉得到提煉數(shù)據(jù)，所述采集數(shù)據(jù)在提煉方向上的維 數(shù)為n，所述提煉數(shù)據(jù)在所述提煉方向上的維數(shù)為t，所述t、η均為正整數(shù)且t〈n ;
[0009] c)由所述提煉數(shù)據(jù)和所述校準數(shù)據(jù)計算得到線圈合并系數(shù)。
[0010] 優(yōu)選的，所述步驟a)中k空間數(shù)據(jù)的數(shù)據(jù)方向包括以下任意一種或多種：頻率編 碼方向、相位編碼方向、通道方向、時間方向、3d掃描中的第二相位編碼方向或上述多個方 向合并形成的方向。
[0011] 優(yōu)選的，所述步驟b)中數(shù)據(jù)提煉方法為矩陣降維法。
[0012] 優(yōu)選的，所述步驟b)中矩陣降維法包括如下步驟：
[0013] 將所述米集數(shù)據(jù)的第一方向向量作為基礎(chǔ)向量，所述第一方向與提煉方向垂直；
[0014] 計算所述采集數(shù)據(jù)的協(xié)方差矩陣，并選取所述協(xié)方差矩陣歸一化正交特征矢量中 特征值較大的前t個值，組成提煉系數(shù)；
[0015] 由所述提煉系數(shù)計算得到提煉數(shù)據(jù)。
[0016] 優(yōu)選的，所述步驟b)中矩陣降維法包括如下步驟：
[0017] 將所述采集數(shù)據(jù)的第一方向向量作為基礎(chǔ)向量，所述第一方向與提煉方向垂直；
[0018] 按照相同數(shù)據(jù)方向?qū)⑺霾杉瘮?shù)據(jù)分類得到各類第一方向向量〇bk_Sxl、0bk_ Sx2、…、0bk_Sxi ;
[0019] 由所述各類第一方向向量計算各類樣本第一方向向量均值〇bk_Sml、0bk_Sm2、…、 0bk_Smi和總樣本第一方向向量均值0bk_Sm，并進而計算得到樣本類內(nèi)離散度矩陣0bk_ Ssl、0bk_Ss2、...、0bk_Ssi 以及樣本間離散度矩陣 0bk_Ssb ;
[0020] Obk^SsrkObk^Ssb的前t個最大特征值對應(yīng)的特征向量分別為wil、wi2、…、wit, 將wil、wi2至wit組合成一個矩陣作為所述提煉數(shù)據(jù)。
[0021] 進一步的，所述第一方向為所述采集數(shù)據(jù)的列向量方向或行向量方向。
[0022] 本發(fā)明還提供了一種磁共振并行成像方法，包括在所述的磁共振并行成像線圈合 并系數(shù)計算方法計算得到線圈合并系數(shù)后，根據(jù)所述線圈合并系數(shù)和所述采集數(shù)據(jù)，填補 得到完整的k空間數(shù)據(jù)，將所述完整的k空間數(shù)據(jù)變換到圖像域得到圖像。
[0023] 本發(fā)明還提供了一種磁共振并行成像裝置，包括：
[0024] 采集單元，適于采集k空間數(shù)據(jù)，所述k空間數(shù)據(jù)包括采集數(shù)據(jù)和校準數(shù)據(jù)，所述 k空間數(shù)據(jù)至少包含一個數(shù)據(jù)方向；
[0025] 計算單元，適于對所述采集數(shù)據(jù)進行數(shù)據(jù)提煉得到提煉數(shù)據(jù)，所述采集數(shù)據(jù)在提 煉方向上的維數(shù)為n，所述提煉數(shù)據(jù)在所述提煉方向上的維數(shù)為t，所述t、η均為正整數(shù)且 t〈n，并由所述提煉數(shù)據(jù)和所述校準數(shù)據(jù)計算得到線圈合并系數(shù)；
[0026] 填補單元，適于根據(jù)所述線圈合并系數(shù)和所述采集數(shù)據(jù)，填補得到完整的k空間 數(shù)據(jù)；
[0027] 成像單元，適于將所述完整的k空間數(shù)據(jù)變換到圖像域得到圖像。
[0028] 優(yōu)選的，所述計算單元適于：將所述采集數(shù)據(jù)的第一方向向量作為基礎(chǔ)向量，所述 第一方向與提煉方向垂直；計算所述采集數(shù)據(jù)的協(xié)方差矩陣，并選取所述協(xié)方差矩陣歸一 化正交特征矢量中特征值較大的前t個值，組成提煉系數(shù)；由所述提煉系數(shù)計算得到提煉 數(shù)據(jù)。
[0029] 優(yōu)選的，所述計算單元適于：將所述采集數(shù)據(jù)的第一方向向量作為基礎(chǔ)向量，所述 第一方向與提煉方向垂直；按照相同數(shù)據(jù)方向?qū)⑺霾杉瘮?shù)據(jù)分類得到各類第一方向向量 0bk_Sxl、0bk_Sx2、…、0bk_Sxi ;由所述各類第一方向向量計算各類樣本第一方向向量均 值0bk_Sml、0bk_Sm2、...、0bk_Smi和總樣本第一方向向量均值0bk_Sm，并進而計算得到樣 本類內(nèi)離散度矩陣〇bk_Ssl、Obk_Ss2、…、Obk_Ssi以及樣本間離散度矩陣Obk_Ssb ;Obk_ Ssr1=I=Obk^Ssb的前t個最大特征值對應(yīng)的特征向量分別為wiI、wi2、…、wit,將wiI、wi2 至wit組合成一個矩陣作為所述提煉數(shù)據(jù)。
[0030] 通過本發(fā)明所提供的技術(shù)方案，針對磁共振并行采集重建過程中選取的卷積核， 進行數(shù)據(jù)重提煉，使得提煉后的數(shù)據(jù)量減小，信號特征得到增強，利用提煉后的卷積核，進 行數(shù)據(jù)填補，得到新的卷積核比原有的得到了優(yōu)化。這樣一方面可以一定程度去除噪聲的 影響，使得信號特征強化，優(yōu)化系數(shù)計算的準確性，穩(wěn)定性；更進一步，可以使得卷積核的選 取更加容易。
【附圖說明】
[0031] 圖1為現(xiàn)有技術(shù)線圈合并系數(shù)計算示意圖；
[0032] 圖2為本發(fā)明磁共振并行成像線圈合并系數(shù)計算方法流程圖。
【具體實施方式】
[0033] 為了使本發(fā)明的上述目的、特征、優(yōu)點能夠更為顯而易懂，下面結(jié)合附圖和實施例 對本發(fā)明的【具體實施方式】作進一步描述。
[0034] 本發(fā)明提供了一種磁共振并行成像線圈合并系數(shù)計算方法，圖2為其流程圖，請 參見圖2,所述方法包括如下步驟：
[0035] S101，采集k空間數(shù)據(jù)，所述k空間數(shù)據(jù)包括采集數(shù)據(jù)和校準數(shù)據(jù)；
[0036] S102,對所述采集數(shù)據(jù)進行數(shù)據(jù)提煉得到提煉數(shù)據(jù)，所述采集數(shù)據(jù)在