專利名稱：磁共振成像方法
技術領域：
本發明涉及一種用于形成對象的圖像的磁共振(MR)方法，其中在k空間中采集一組非線性軌跡，然而，所述軌跡組的密度遠低于對應于對象大小的密度。借助于一個或多個接收天線對沿著這些軌跡的信號進行采樣，由這些信號并且基于接收天線組的空間靈敏度分布得出磁共振圖像。本發明尤其涉及一種磁共振成像方法，其中借助于接收天線系統采集磁共振信號并基于磁共振信號重建磁共振圖像。
國際申請WO 01/73463公開了這種磁共振成像方法。
在這種已知的磁共振成像方法中，通過沿著k空間中的軌跡掃描采集磁共振信號。已知的磁共振成像方法在選擇穿過k空間將遵循的采集軌跡方面提供了高自由度。尤其是，可以使用導致k空間中不規則采樣圖案(pattern)的采集軌跡，特別是螺旋形軌跡。
本發明還涉及一種MR設備及執行這種方法的計算機程序產品。
通常，在諸如SENSE(Pruessmann)或SMASH(Sodickson)的平行成像中，分別通過k空間或圖像空間的笛卡爾(Cartesian)網格化執行圖像重建。
在US-A-2003/0122545中描述了一種磁共振成像方法，其中二次采樣度選擇成，使得用于接收由單個RF激發引起的磁共振信號(回波)序列的后繼采集時間少于MR信號的衰減時間。優選地，執行k空間的分段掃描，分段數目和在每個分段中的所掃描的線的數目可調節，并且掃描預定總數目的線。使用少量分段使得接收完整磁共振圖像的磁共振信號的采集時間少于相關動態過程的處理時間。盡管所述方法被描述成在k空間中使用直線掃描軌跡，也可以使用其它軌跡，例如諸如圓弧或螺線的曲線。然而，在這種情況中，需要對磁共振信號進行更復雜的頻率和相位編碼。該文獻未描述對k空間中連續非笛卡爾軌跡的具體解決方案。
在Magn.Reson.Med 10(2002)中M.Bydder等人的論文中提到，通過將問題看作大的線性方程系統，用于重建來自多個線圈的不充分采樣(under-sampled)的k空間數據的部分平行成像技術可以與任意采集方案(例如笛卡爾，螺旋等)一起使用。然而，對于實際應用而言，直接求解這個系統的計算成本很高。目前為止，對于在k空間中連續非笛卡爾軌跡(尤其例如螺旋采樣)上粗糙地不充分采樣的數據的快速重建過程缺乏實際解決方案。
本發明的目的在于進一步減少在磁共振成像方法中從所采集磁共振信號進行重建所涉及的計算工作量。
該目的由本發明的磁共振成像方法實現，其中-借助于接收天線系統通過多個信號通道采集磁共振信號-該接收天線系統具有一靈敏度分布-通過不充分采樣采集磁共振信號-對k空間中的各個取向扇形區域，在規則采樣網格上從不充分采樣采集的磁共振信號重采樣規則重采樣的磁共振信號-所述重采樣涉及不充分采樣的采集磁共振信號與網格化核函數(gridding kernel)的卷積-網格化核函數取決于-所關心扇形區域的取向，和-接收天線系統的靈敏度分布，以及-根據重采樣磁共振信號重建磁共振圖像。
本發明基于下面的理解。為了實現快速重建，使用諸如快速傅立葉變換(FFT)的技術。作為輸入，這些技術需要數據在規則采樣網格上被采樣。而且，k空間中各種各樣的采集軌跡，尤其是螺旋形軌跡和包括螺旋形片段的軌跡，由k空間的各扇形區域中的軌跡的(幾乎)平行的軌跡片段精確地或至少相當地近似。通常，扇形區域可以是具有穿過k空間原點的主軸的k空間區域。這種扇形區域在角邊界之間，即在k矢量的各個最小和最大模數之間延伸到k空間的邊界，并且由從k空間原點徑向延伸的徑向邊界限制。扇形區域可以是從k空間原點或穿過k空間原點延伸到k空間邊界的完整扇形。在二維中，扇形區域是穿過k空間原點點對稱地延伸的扇形或平面扇形或扇形片段，在三維中，扇形區域是k空間中的錐形或錐形的部分。根據本發明，對單個扇形區域單獨執行重建，該重建涉及重新網格化以重采樣采集的磁共振信號，從而在規則網格上重采樣磁共振信號。重新網格化包括與網格化核函數卷積。網格化核函數取決于所關心扇形區域的取向，以解釋圖像空間中合適的方向，由于磁共振圖像的像素值和k空間中重新采樣磁共振信號之間的傅立葉關系將在所述圖像空間中發生混疊。另外，網格化核函數包括接收天線系統的靈敏度分布以計入由采集的磁共振信號的不充分采樣導致的混疊。基于扇形區域的取向和基于靈敏度分布導出網格化核函數不需很多的計算量。由于相對小數量的像素或體素導致混疊，計算量顯著減少。因此，僅相對低維度的矩陣需要矩陣求逆。在諸如笛卡爾方形網格的規則網格上的實際重新采樣僅包括與網格化核函數的卷積，這僅花費很少的計算量。然后通過僅占用很短計算時間的FFT技術執行磁共振圖像的最后重建。
看來，對網格化核函數的k空間中的扇形區域的取向的依賴關系涉及十分平滑的變化。因此，特定取向的網格化核函數對k空間中相當寬的扇形也準確有效。
本發明的另一目的是提供一種磁共振成像方法，該方法能夠實現對粗糙地不充分采樣的k空間中非笛卡爾采樣(尤其是沿螺旋軌跡的)進行快速重建。本發明的另一目的是提供一種執行該方法的系統和計算機程序產品。
這個目的通過根據特別在權利要求1、2和3中所述的本發明的磁共振成像方法實現。
這些目的通過權利要求1所述的方法、權利要求6所述的MR設備和權利要求21所述的計算機程序產品來實現。
本發明的主要優點在于得到了代表傅立葉域中線圈靈敏度信息的SPACERIP和非笛卡爾SENSE的公式。由于所需傅立葉項的數量少，線性系統非常稀疏，從而允許有效求解方程式。從而，使得螺旋掃描在高度不充分采樣下是可行的，從而實現非常快速的采集和重建。
本發明的主要方面在于k空間中的非笛卡爾軌跡可以通過虛構平行切線的坐標系局部描述，以便執行類似SENSE或SMASH的二次采樣，所述切線局部形成笛卡爾網格。如果整個k空間被在360°的角度上均勻分割的射線再次分割，則得到連續的局部笛卡爾網格系統。那么k空間的這些部分被局部重建并作為整體轉換為圖像。
本發明的這個和其它優點在從屬權利要求和后面的描述中披露，在后面的描述中參考附圖描述本發明的示例性實施例。附圖中示出

圖1是k空間中不充分采樣的螺旋軌跡，圖2是與圖1相同的螺旋軌跡，其中有用于具有角θ的半徑周圍的區域的假想平行掃描線，圖3是對應圖2所示的假想笛卡爾采樣圖案的圖像中的折疊點，圖4是用于執行根據本發明的方法的設備，以及圖5是圖4所示設備的電路圖。
本發明的原理在圖1中，已經示出了螺旋掃描軌跡，其是在單次激發(singleshot)EPI中的單螺旋臂2。點3表示笛卡爾網格，其具有正確成像包含待成像對象的視野(FOV)所需的密度。實際密度也可以稍高(所謂的“過網格化”)，對應于稍大于對象的區域。螺旋臂2已經根據SENSE方法用約為2的不充分采樣因數被粗糙地不充分采樣。這能夠立即從圖1看出，因為臂的螺旋部分之間的距離是大約兩個點3的距離。從尼奎斯特準則的觀點來看，這是不充分采樣。然而，如果那個軌跡的信號已經由至少兩個具有不同空間靈敏度圖案(pattern)的接收天線或線圈采樣，圖像仍然能夠被重建。圖像的重建需要解出具有N個未知數的大致N個方程的方程組，其中N是大約采樣點的數目乘以線圈的數目的量值，或大約所得到圖像中象素的數目(即，N為約104到105)。這意味著應當解出N個方程式，這不能通過直接的矩陣反演(matrix inversion)得到。因此，一些作者提議使用迭代解，這需要大約10次迭代，每次涉及昂貴的計算網格化操作。
根據本發明算法的概要在圖2中，在螺旋臂2中示出了具有角θ的徑向線5，其橫穿螺旋臂2。在徑向線5和螺旋臂2之間的交點繪出了切線6a、6b、6c和6d，這些切線顯示了螺旋臂2的相鄰部分大致平行且等距。在此示出了阿基米德螺旋(Archimedic spiral)。但是也可以使用其它螺旋函數。這種情況可以從平行成像中的笛卡爾方法得知如果這些等距線6a到6b覆蓋了整個k空間，那么重建將更輕松。在圖像空間中，不連續數目個對象點將被“簡單地”相疊，如在圖3中為簡便以兩個點所示的那樣。
原則上，該問題能夠通過“常規”SENSE重建而被解決。這可以寫成接收天線信號mk(X，Y)“加權”函數ak(X，Y)的和。這也可以在傅立葉域中寫成p(X,Y)=Σcoils kF-1{αk(kx,ky)⊗μk(kx,ky)},]]>或 (1)p(X,Y)=F-1{Σcoils kαk(kx,ky)⊗μk(kx,ky)}---(2)]]>其中μk(kx，ky)是沿假想等距線6a到6b的測量數據，α是ak(X，Y)的傅立葉變換。
應當注意，方程式(1)或(2)確切描述了在沒有不充分采樣(SENSE，SMASH)的情況下為正常螺旋成像進行的相同操作。此時，ak(X，Y)的意思是“網格化核函數”，其實質上是框的傅立葉變換(但是與平滑邊緣一起逐漸減小以防止ak(X，Y)具有大的支集(support))。
在目前情況中，ak(X，Y)的形狀不是錐形框，而是一“重建函數”，其主要取決于所有接收天線或線圈的線圈靈敏度圖案，取決于SENSE方法的折疊距離，最終部分取決于對象形狀(由于正則化)。然而，由于線圈靈敏度函數被期望在空間中是平滑的函數，函數ak(X，Y)也被期望在空間中是平滑的。因此，網格化函數αk(X，Y)被期望具有相對小的支集。假定12*12到16*16笛卡爾點的支集將是足夠的(其中對于正常成像的網格化，4*4的支集通常足夠)。
所獲得的網格化函數αk(X，Y)可以完美地應用于來自一組平行等距線的重建數據，所述等距線相對所需網格具有角度。然而，在這種情況中，沿著螺旋臂采樣數據，而不是沿著線。這意味著所獲得的網格化核函數僅對嚴格定位在具有角θ的半徑上的點有效。嚴格地說，網格化核函數ak(X，Y)應當對無窮大情況進行計算。然而，線圈靈敏度圖案通常是空間的平滑函數。這意味著如果折疊方向略微改變，加權函數ak(X，Y)(并且因此，網格化函數αk(X，Y))將不會明顯改變。“折疊方向”由折疊點的線之間的角限定，或等同地，由假想平行線6a到6d的取向限定。因此，所獲得的網格化函數能夠應用在具有角θ的半徑周圍的預定區域中。這允許對有限數量的(例如10或20個)半徑計算αk(X，Y)。
得出的算法假定在整個有關區域已知線圈靈敏度，而且對對象的存在有一些認識(如同在笛卡爾SENSE中一樣)。假定在二維k空間中的螺旋軌跡。那么相關的問題僅僅在于螺旋臂之間的距離。根據本發明的重建將通過以下步驟進行1.選擇笛卡爾網格。距離(或密度)應當符合所研究對象的大小。
2.在k空間上限定多個等距半徑(例如10或20個)。
3.對每個半徑，確定與螺旋臂相切的線6a到6b的方向，由此切線之間的距離應當與角度無關。
4.對所獲得的折疊點集合，計算SENSE重建矩陣A＝(Sh·Ψ-1·S+R-1)-1·Sh·Ψ-1，其中S是接收天線或線圈靈敏度矩陣，Ψ是線圈通道之間的噪聲協方差，R是正則化矩陣，Sh表示S的厄米共軛。將矩陣A的片段組合成每個接收天線或線圈k的函數ak(X，Y)。(應當注意的是該步驟是常規SENSE重建的一部分)5.在圖像區域外，將函數ak(X，Y)設為零(補零)；然后執行這個函數的傅立葉變換，成為αk(X，Y)，或者使用在足夠精細的網格上內插αk(X，Y)的值的其它合適的方法，以用于隨后的卷積和再采樣。這個函數預期僅對于kx和ky的小值是相關，從而其外面部分可以被丟棄。
6.采集數據。
7.執行采樣密度補償。
8.對于沿著螺旋軌跡的每個點，確定最靠近那個采樣點的半徑。
9.使用為每個線圈的最近半徑計算的網格化函數αk(X，Y)，以及每個線圈的采集數據樣本，執行網格化操作。這個操作是常規螺旋掃描重建過程的一部分，僅αk的范圍可能更大。
10.對線圈元件上的數據求和。
11.對笛卡爾網格點進行傅立葉變換。
應當注意，對于動態掃描(或其中對相同幾何位置采集多個數據集的任何類型的掃描)，步驟1-5必須只執行一次。
算法的細化R1.可以選擇不對應于正常FOV(或“所研究對象的大小”)而是對應于稍大FOV的網格，從而網格稍密。已知在用戶選擇的FOV和更大或“過網格化”FOV之間的空白中不包含對象。朝向過網格化FOV的邊緣，正則化值R逐漸逼近零。在這種方式中，避免了函數ak(X，Y)的不連續，導致αk(X，Y)的更小支集。
R2.函數ak(X，Y)可以通過首先將它們與共用整形函數(shapingfunction)相乘被預處理，以便減少αk(X，Y)的支集。這可以是錐形窗口函數，或與例如靈敏度的平方和相乘，以避免在所有線圈不靈敏的空間中的點上存在極大的值。
R3.可以采用兩個最接近的半徑并且它們之間的網格化核函數可以被內插。在更有效的方式中，兩個半徑都被網格化并且其結果被內插。
R4.k空間的最中心區域可以通過另一方法(例如，直接求逆)重建，并且網格化結果與中心區域的可選重建混和。
R5.函數ak(X，Y)被分成規定數量的子函數，對于這些子函數，相應的 的支集特別小，其總和等于或近似于完整視野中的原始函數ak(X，Y)。這使得可以應付ak(X，Y)的急劇變化，這種急劇變化尤其發生在減小的視野的邊緣。在這種情況下，需要對每個子函數集進行單獨的網格化。與充分擴大局部卷積核的尺寸以表現急劇變化相比，這在計算上是有吸引力的。一個自然選擇是分別采用每個減小的視野，并用一些類型的錐形窗口函數處理具有或沒有周期重復的函數。
R6.由于類似原因，分配給相鄰半徑的樣本集合可以被單獨網格化和變換。這允許為完整FOV提供誤差“空間”，非常類似于傳統的過采樣在完整FOV邊緣提供的誤差“空間”。
R7.所提出的算法可以用于產生一圖像，然后使用該圖像初始化已知的迭代重建算法之一(例如，共軛梯度方法)。這可以顯著減少這些方法實現滿意的圖像質量所需的迭代次數。同樣，對于合理的參數設定(即，αk(kx，ky)的有限支集和半徑的有限數目)，這個策略可以消除可能伴隨所提出算法的任何偽影。
方法的擴展原則上，每個k空間采樣點的局部鄰域和二次采樣的局部程度可以被單獨考慮。在這種情況下，對具有類似局部性質的點的集合執行根據本發明的方法的步驟1-5，所述點集可以在k空間中任意分布。這將尤其允許將所提出的算法也應用于可變密度螺旋和傳統徑向采集。
圖4所示的設備是MR設備，該設備包括用于產生穩定、均勻磁場的四個線圈51的系統，該磁場的場強的大小大約是從零點幾特斯拉(Tesla)到幾個特斯拉。相對于z軸同心布置的線圈51可以設置在球表面52上。待檢查的患者60處于定位在這些線圈內的工作臺54上。為了產生在z方向延伸并在該方向線性變化的磁場(該磁場在下文中也被稱為梯度場)，在球表面52上提供作為多個接收天線的4個線圈53。還具有產生梯度場的4個線圈57，該梯度場也在x方向(垂直)延伸。由4個線圈55產生在z方向延伸并在y方向(垂直于圖面)具有梯度的梯度磁場，線圈55可以和線圈57相同，但是被布置成在空間中相對其偏轉90°。在此僅示出這4個線圈中的兩個。
由于用于產生梯度磁場的3個線圈系統53、55、57中的每個相對球面對稱布置，球中心處的場強由線圈51的穩定、均勻磁場唯一確定。還提供了產生基本均勻的RF磁場的RF線圈61，該RF磁場垂直于穩定、均勻磁場的方向(即，垂直于z方向)延伸。在每個RF脈沖期間，該RF線圈從RF發生器接收RF調制電流。RF線圈61也可以用于接收在檢查區域中產生的自旋共振信號。
如圖5所示，在MR設備中接收的MR信號由單元70放大并在基帶中換位(transpose)。由此獲得的模擬信號通過模數轉換器71轉換成數字值序列。模數轉換器71受控制單元69控制使得它僅在讀出相期間產生數字數據字。模數轉換器71之后是傅立葉變換單元72，該傅立葉變換單元對由MR信號數字化獲得的采樣值序列進行一維傅立葉變換，執行快速使得在下個MR信號接收到之前傅立葉變換終止。
由傅立葉變換如此產生的原始數據被寫入存儲器73，其存儲容量足夠存儲若干原始數據集。合成單元74從這些原始數據集以所述方式產生合成圖像；這個合成圖像被儲存在存儲器75中，該存儲器的儲存容量足夠儲存大量連續合成圖像80。這些數據集在不同時刻被計算，這些時刻之間的間隔優選比采集一數據集所需的測量時間小。執行連續圖像合成的重建單元76根據由此采集的數據集產生MR圖像，所述MR圖像被儲存。這些MR圖像代表預定時刻的檢查區域。從數據由此獲得的MR圖像系列適當地再現檢查區域中的動態過程。
單元70-76受控制單元69控制。如向下的指示箭頭所指示的，控制單元也對梯度線圈系統53、55和57中的電流時間和由RF線圈61產生的RF脈沖的中心頻率、帶寬和包絡強加變化。存儲器73和75以及重建單元76中的MR圖像存儲器(未示出)可以由具有足夠容量的單個存儲器實現。傅立葉變換單元72、合成單元74和重建單元76可以由適合運行根據上述方法的計算機程序的數據處理器實現。
權利要求
1.一種磁共振成像方法，其中-借助于接收天線系統通過多個信號通道采集磁共振信號-該接收天線系統具有靈敏度分布-在不充分采樣的情況下采集磁共振信號-對k空間中各個取向的扇形區域，在規則采樣網格上從不充分采樣采集的磁共振信號重采樣規則重采樣的磁共振信號-所述重采樣涉及不充分采樣采集的磁共振信號與網格化核函數的卷積-該網格化核函數取決于-所關心扇形區域的取向，和-接收天線系統的靈敏度分布，以及-根據重采樣的磁共振信號重建磁共振圖像。
2.根據權利要求1的磁共振成像方法，其中磁共振信號通過沿著非線性、尤其是螺旋形軌跡掃描k空間被采集。
3.根據權利要求1的磁共振成像方法，為了形成對象的圖像，其中-從二次采樣的磁共振信號并基于多個接收天線的空間靈敏度分布得到磁共振圖像，-施加RF脈沖和梯度序列，該序列對應于k空間中的包括至少一個基本非線性軌跡的一組軌跡，其中所述軌跡組的采樣密度顯著小于與對象尺寸對應的正常采樣密度，-沿著所述軌跡組的每個信號在至少兩個不同接收天線位置被采樣，得到多個接收天線信號，-通過與網格化核函數卷積，將所述信號的數據從所述軌跡組轉換成笛卡爾網格，以重建圖像，并且由此-網格化核函數是針對笛卡爾網格為每個天線加權的圖案的傅立葉變換，并且-網格化核函數圖案在k空間中的一個區域和另一區域之間不同。
4.根據權利要求1的方法，其中加權圖案由于以下原因被獲得-對k空間的每個單獨區域，指定一組平行等距線，這些線局部匹配所述軌跡組，-確定圖像空間中重疊點的圖案，其對應于k空間中的該平行等距線組，-在圖像空間中，計算每個天線的加權圖案，該圖案近似唯一地對應于k空間的單獨區域中的所述平行等距線的圖案。
5.根據權利要求2的方法，其中該軌跡組的至少部分對應于阿基米德螺旋，k空間中的區域由它們在k空間中的方位角限定。
6.根據權利要求4或5的方法，其中根據二次采樣數據和接收天線的空間靈敏度分布的笛卡爾方程組的反演，計算天線的加權圖案。
7.根據權利要求6的方法，其中所述笛卡爾方程組的反演被公式化為A＝(Sh·Ψ-1，S+R-1)-1·Sh·Ψ-1，其中A是重建矩陣，S是接收天線靈敏度矩陣(sij)，其中sij是天線i在點的重疊組的第j個點上的空間靈敏度分布，Ψ是天線之間的噪聲協方差，R是正則化矩陣，Sh表示S的厄米共軛。
8.根據權利要求7的方法，其中選擇網格化核函數以當正常FOV覆蓋所研究對象的大小時對應較大的FOV，并且該較大FOV的空白和正常FOV之間的正則化矩陣R的值設為零。
9.根據權利要求8的方法，其中從重建矩陣A得到的用于每個天線的網格化核函數圖案與共用整形函數相乘，該共用整形函數包括錐形窗口函數或每個天線的靈敏度平方和。
10.根據權利要求5-9的方法，其中穿過螺旋軌跡的兩個最接近半徑之間的網格化核函數被內插。
11.根據權利要求10的方法，其中兩個半徑都被網格化，并且它們的結果被內插。
12.根據權利要求1-11之一的方法，其中k空間的最中心區域通過直接反演以完全采樣密度被重建，網格化重建方法的結果與以完全采樣密度重建的結果混和。
13.根據權利要求7-12之一的方法，其中從重建矩陣A得到的用于每個天線的網格化核函數圖案被分成規定數量的子函數，對于所述子函數，k空間中對應函數的支集趨向零，以便丟棄網格化核函數圖案中的急劇變化，而每個子函數被單獨網格化。
14.根據權利要求10的方法，其中分配給相鄰半徑的樣本集被單獨網格化和變換。
15.由權利要求1-14之一的方法產生的圖像的用途，其目的是初始化圖像重建的傳統迭代算法。
16.一種用于從多個信號獲得MR圖像的磁共振成像設備，包括-主磁體，-用于在患者的預定區域內激發自旋的發射天線，-用于在主磁場的有限均勻區域內采樣信號的多個接收天線，-承載患者的工作臺，-用于連續移動工作臺穿過主磁體孔的裝置，-用于從二次采樣磁共振信號并基于每個所述接收天線位置的空間靈敏度分布得到磁共振圖像的裝置，-用于施加RF脈沖和梯度序列的裝置，該序列對應于k空間中的包括至少一個基本非線性軌跡的一組軌跡，其中所述軌跡組的密度顯著小于與對象尺寸對應的密度，-用于在至少兩個不同接收天線位置沿著所述軌跡組采樣每個信號以得到多個接收天線信號的裝置，-用于通過與網格化核函數卷積將所述信號的數據從所述軌跡組轉換成笛卡爾網格而重建圖像的裝置，并且由此-網格化核函數對于每個天線是特定的，-網格化核函數圖案在k空間中的一個區域和另一區域之間不同，且-網格化核函數是針對笛卡爾網格為每個天線加權的圖案的傅立葉變換。
17.根據權利要求16的設備，還包括-用于獲得加權圖案的裝置，其包括-用于為k空間的每個單獨區域指定一組平行等距線的裝置，這些線局部匹配所述軌跡組，-用于確定圖像空間中重疊點的圖案的裝置，該圖案對應于k空間中的該平行等距線組，以及-用于在圖像空間中計算每個天線的加權圖案的裝置，該圖案近似唯一地對應于k空間的單獨區域中的所述平行等距線的圖案。
18.根據權利要求17的設備，還包括通過它們在k空間中的方位角限定k空間中的區域的裝置，而所述軌跡組的至少部分對應于等距螺線。
19.根據權利要求17或18的設備，還包括用于計算天線的加權圖案的裝置，其根據二次采樣數據和接收天線空間靈敏度分布的笛卡爾方程組的反演進行此計算。
20.根據權利要求19的設備，所述用于計算所述笛卡爾方程組的反演的裝置基于公式A＝(Sh·Ψ-1·S+R-1)-1·Sh·Ψ-1，其中A是重建矩陣，S是接收天線靈敏度矩陣(sij)，其中sij是天線i在點的重疊組的第j個點上的空間靈敏度分布，Ψ是天線之間的噪聲協方差，R是正則化矩陣，Sh表示S的厄米共軛。
21.一種儲存在計算機可用介質上的通過磁共振方法形成圖像的計算機程序產品，包括用于使計算機控制以下執行的計算機可讀程序裝置-通過主磁體產生主磁場，-通過發射天線在患者預定區域中激發自旋，-在多個接收天線位置在主磁場的有限均勻區域內采樣多個信號，-連續移動支承患者的工作臺穿過主磁體的孔，-從二次采樣磁共振信號并基于每個所述接收天線位置的空間靈敏度分布得到磁共振圖像，-施加RF脈沖和梯度序列，該序列對應于k空間中的包括至少一個基本非線性軌跡的一組軌跡，其中所述軌跡組的密度顯著小于與對象尺寸對應的密度，-在至少兩個不同接收天線位置沿著所述軌跡組采樣每個信號，得到多個接收天線信號，-通過與網格化核函數卷積將所述信號的數據從所述軌跡組轉換成笛卡爾網格而重建圖像，并且由此-網格化核函數對于每個天線是特定的，-網格化核函數圖案在k空間中的一個區域和另一區域之間不同，并且-網格化核函數是針對笛卡爾網格為每個天線加權的圖案的傅立葉變換。
全文摘要
本發明描述了一種新穎的磁共振成像方法和設備，其中從二次采樣磁共振信號和基于各接收天線的空間靈敏度分布得到圖像。施加RF脈沖和梯度序列，該序列對應k空間中的包括至少一個基本非線性軌跡的一組軌跡，其中所述軌跡組的密度顯著小于與對象尺寸對應的密度。在至少兩個不同接收天線位置沿著所述軌跡組采樣每個信號。通過用網格核函數卷積將所述信號的數據轉換成笛卡爾網格而重建圖像，由此網格化核函數對每個天線是特定的，在k空間中的一個區域和另一區域之間不同，并且是針對笛卡爾網格為每個天線加權的圖案的傅立葉變換。
文檔編號G01R33/561GK1918480SQ200580004408
公開日2007年2月21日 申請日期2005年2月3日 優先權日2004年2月10日
發明者M·富德爾, H·埃格斯 申請人:皇家飛利浦電子股份有限公司