本發(fā)明涉及磁共振成像領(lǐng)域，特別涉及一種磁共振成像方法與裝置。

背景技術(shù)：
在磁共振成像技術(shù)中，成像的速度是衡量成像方法的一個(gè)很重要標(biāo)準(zhǔn)。限制成像速度的很重要因素是數(shù)據(jù)采集以及K空間填充。一般的數(shù)據(jù)采集方式要采滿K空間數(shù)據(jù)，然后才能進(jìn)行重建得到圖像。磁共振并行采集重建技術(shù)，是利用線圈重組合并的方式，對(duì)欠采樣的數(shù)據(jù)進(jìn)行填補(bǔ)，利用填補(bǔ)完整的K空間數(shù)據(jù)進(jìn)行重建。利用這樣的方式，可以根據(jù)需求，只采集一部分K空間數(shù)據(jù)，不必采滿整個(gè)K空間，由此可以大大加快成像的速度。比較常用的并行重建方法之一是廣義自校準(zhǔn)并行采集（GRAPPA，GeneRalizedAutocalibratingPartiallyParallelAcquisitions）。傳統(tǒng)的GRAPPA算法如圖1所示，黑色實(shí)點(diǎn)代表為實(shí)際采集的K空間數(shù)據(jù)；白色空點(diǎn)為欠采樣需要填補(bǔ)的數(shù)據(jù)；灰色實(shí)點(diǎn)表示為了計(jì)算線圈合并系數(shù)而適量全采樣的數(shù)據(jù)。GRAPPA算法認(rèn)為，圖1中任意一個(gè)白色空點(diǎn)可以表示為周圍黑色實(shí)點(diǎn)的線性疊加，相當(dāng)于對(duì)多個(gè)線圈的數(shù)據(jù)進(jìn)行了合并，而線圈合并系數(shù)nij（第i個(gè)線圈，第j個(gè)位置）可以通過黑色實(shí)點(diǎn)擬合灰色實(shí)點(diǎn)來確定。在線圈合并系數(shù)nij確定后，其他白色空點(diǎn)即可根據(jù)求得的線圈合并系數(shù)nij將線圈合并填補(bǔ)空白數(shù)據(jù)。需要說明的是，上述的黑色、白色及灰色點(diǎn)通常指的是某一方向（一般是讀方向，即頻率編碼方向）上的線數(shù)據(jù)集，如圖2所示的第i個(gè)線圈的K空間數(shù)據(jù)分布示意圖，其中，虛線代表圖1中的白色數(shù)據(jù)點(diǎn)，黑色粗實(shí)線代表圖1中的黑色數(shù)據(jù)點(diǎn)，灰色細(xì)實(shí)線代表圖1中的灰色數(shù)據(jù)點(diǎn)。按照傳統(tǒng)的方法，數(shù)據(jù)分為三種：應(yīng)采集數(shù)據(jù)、欠采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)。應(yīng)采集數(shù)據(jù)如圖2黑色粗實(shí)線所示，欠采集數(shù)據(jù)如圖2虛線所示，校準(zhǔn)數(shù)據(jù)如圖2灰色細(xì)實(shí)線所示，并且每一種數(shù)據(jù)都應(yīng)為整條線數(shù)據(jù)，校準(zhǔn)數(shù)據(jù)至少為一條線數(shù)據(jù)?，F(xiàn)有的磁共振成像方法，通常先通過上述并行采集重建方法對(duì)K空間域進(jìn)行重建，然后再通過某種變換操作（如傅里葉變換）將重建之后的K空間域變換至圖像域以獲得圖像。一般情況下，這種變換操作會(huì)消耗很多的時(shí)間，導(dǎo)致成像速度較慢，而且計(jì)算量較大。相關(guān)技術(shù)還可參考公開號(hào)為US2006184000A1的美國專利申請(qǐng)，該專利申請(qǐng)公開了一種磁共振成像快速廣義自校準(zhǔn)并行采集圖像重建算法。

技術(shù)實(shí)現(xiàn)要素：
本發(fā)明要解決的問題是現(xiàn)有磁共振成像方法中成像速度較慢且計(jì)算量較大。為解決上述問題，本發(fā)明技術(shù)方案提供一種磁共振成像方法，包括：通過第一變換操作將欠采集的K空間域變換至中間域，所述中間域位于K空間域與圖像域之間；分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算；在所述中間域中，基于第二應(yīng)采集數(shù)據(jù)以及所述線圈合并系數(shù)進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)；所述第二應(yīng)采集數(shù)據(jù)由對(duì)在K空間域中采集的第一應(yīng)采集數(shù)據(jù)進(jìn)行所述第一變換操作后得到；通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域變換至圖像域以獲得圖像?？蛇x的，所述分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算是在所述中間域進(jìn)行，所述校準(zhǔn)數(shù)據(jù)為對(duì)在K空間域中采集的第一校準(zhǔn)數(shù)據(jù)進(jìn)行所述第一變換操作后得到的第二校準(zhǔn)數(shù)據(jù)?？蛇x的，所述分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算是在K空間域進(jìn)行，所述校準(zhǔn)數(shù)據(jù)為在K空間域中采集的第一校準(zhǔn)數(shù)據(jù)?？蛇x的，所述數(shù)據(jù)線性加權(quán)加和操作中引入或未引入頻率編碼方向的數(shù)據(jù)?？蛇x的，所述第一變換操作為頻率編碼方向的一維傅里葉變換，所述第二變換操作為相位編碼方向的一維傅里葉變換。為解決上述問題，本發(fā)明技術(shù)方案還提供一種磁共振成像裝置，包括：第一變換單元，適于通過第一變換操作將欠采集的K空間域變換至中間域，所述中間域位于K空間域與圖像域之間；計(jì)算單元，適于分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算；填補(bǔ)單元，適于在所述中間域中，基于第二應(yīng)采集數(shù)據(jù)以及所述線圈合并系數(shù)進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)；所述第二應(yīng)采集數(shù)據(jù)由對(duì)在K空間域中采集的第一應(yīng)采集數(shù)據(jù)進(jìn)行所述第一變換操作后得到；第二變換單元，適于通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域變換到圖像域以獲得圖像。與現(xiàn)有技術(shù)相比，本發(fā)明技術(shù)方案至少具有以下優(yōu)點(diǎn)：先通過第一變換操作將欠采集的K空間域變換至中間域，在所述中間域中完成填補(bǔ)欠采集數(shù)據(jù)的操作，再通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)的中間域變換至圖像域以實(shí)現(xiàn)圖像的并行采集重建，由于第一變換操作中無需對(duì)填補(bǔ)的欠采集數(shù)據(jù)進(jìn)行變換，因此減少了第一變換操作在計(jì)算上的開銷，總體上加快了計(jì)算速度，從而提高了磁共振成像的速度。附圖說明圖1是現(xiàn)有技術(shù)的GRAPPA算法的示意圖；圖2是第i個(gè)線圈的K空間數(shù)據(jù)分布示意圖；圖3是K空間域轉(zhuǎn)換至圖像域的示意圖；圖4是本發(fā)明實(shí)施方式提供的磁共振成像方法的流程示意圖；圖5是頻率編碼方向影響下K空間域的數(shù)據(jù)分布示意圖；圖6是本發(fā)明實(shí)施方式提供的磁共振成像裝置的結(jié)構(gòu)示意圖。具體實(shí)施方式為使本發(fā)明的上述目的、特征和優(yōu)點(diǎn)能夠更為明顯易懂，下面結(jié)合附圖對(duì)本發(fā)明的具體實(shí)施方式做詳細(xì)的說明。在以下描述中闡述了具體細(xì)節(jié)以便于充分理解本發(fā)明。但是本發(fā)明能夠以多種不同于在此描述的其它方式來實(shí)施，本領(lǐng)域技術(shù)人員可以在不違背本發(fā)明內(nèi)涵的情況下做類似推廣。因此本發(fā)明不受下面公開的具體實(shí)施方式的限制。如背景技術(shù)所述，現(xiàn)有的磁共振成像方法，通過某種變換操作將并行采集重建后的K空間域變換至圖像域的過程中，由于這種變換操作會(huì)消耗很多的時(shí)間，導(dǎo)致成像速度較慢，而且計(jì)算量較大。發(fā)明人經(jīng)研究后分析如下：現(xiàn)有的并行采集重建K空間的過程中，計(jì)算線圈合并系數(shù)和填補(bǔ)欠采集數(shù)據(jù)的操作全部是在K空間域中進(jìn)行的，最終要得到圖像，還需要對(duì)填補(bǔ)好的K空間域進(jìn)行某種變換操作（如傅里葉變換）使其轉(zhuǎn)換到圖像域，該變換操作的對(duì)象包含采集的數(shù)據(jù)（即應(yīng)采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)）以及計(jì)算填補(bǔ)的數(shù)據(jù)（即填補(bǔ)的欠采集數(shù)據(jù)）。由于該變換操作過程中需要變換的數(shù)據(jù)量較大，會(huì)消耗很多的時(shí)間，所以如何盡量減少空間變換操作，是提高圖像產(chǎn)生速度的一個(gè)很重要的手段。從K空間域轉(zhuǎn)換到圖像域的過程，一般可以分解為如圖3所示的兩個(gè)子過程：（a）通過變換操作1將K空間域轉(zhuǎn)換至某個(gè)中間域；（b）通過變換操作2將該中間域轉(zhuǎn)換至圖像域。在實(shí)際實(shí)施時(shí)，一般都會(huì)將變換操作1和變換操作2合為一個(gè)變換操作過程，若該變換操作過程以二維傅里葉變換為例，則圖3所示的中間域，可以認(rèn)為是對(duì)K空間域進(jìn)行了一維傅里葉變換（對(duì)應(yīng)變換操作1）后的數(shù)據(jù)空間，將得到的數(shù)據(jù)空間再次進(jìn)行一維傅里葉變換（對(duì)應(yīng)變換操作2）便轉(zhuǎn)換至圖像域。因此，發(fā)明人考慮，如果在K空間域不進(jìn)行對(duì)欠采集數(shù)據(jù)進(jìn)行填補(bǔ)的操作，而是通過變換操作1將欠采集的K空間域先變換至某個(gè)非K空間域的數(shù)據(jù)空間，該數(shù)據(jù)空間即上述的中間域，然后在所述中間域中實(shí)現(xiàn)對(duì)欠采集數(shù)據(jù)進(jìn)行填補(bǔ)的操作，再通過變換操作2將完成填補(bǔ)操作之后的中間域變換至圖像域，這樣便可以從總體上減少變換操作過程中的計(jì)算量，提高磁共振成像的速度，理由如下：由于一般的并行重建方法（如GRAPPA）是在K空間域?qū)η凡杉瘮?shù)據(jù)進(jìn)行填補(bǔ)，填補(bǔ)的欠采集數(shù)據(jù)、應(yīng)采集數(shù)據(jù)以及校準(zhǔn)數(shù)據(jù)都需要進(jìn)行某種變換操作（如二維傅里葉變換），才能轉(zhuǎn)換為圖像域的數(shù)據(jù)；而將填補(bǔ)欠采集數(shù)據(jù)的過程放在非K空間域（即本發(fā)明實(shí)施方式中所述的中間域）進(jìn)行，則填補(bǔ)的欠采集數(shù)據(jù)是所述中間域中的數(shù)據(jù)，該數(shù)據(jù)不需要進(jìn)行從K空間域到中間域的變換操作1，這樣將節(jié)省浪費(fèi)在變換操作1上的計(jì)算開銷，由此能加快計(jì)算速度，進(jìn)而提高磁共振成像的速度。基于上述分析，本發(fā)明實(shí)施方式提供一種磁共振成像方法，在非K空間域的某一個(gè)數(shù)據(jù)空間內(nèi)，進(jìn)行數(shù)據(jù)線性加權(quán)加和操作以完成對(duì)欠采集數(shù)據(jù)的填補(bǔ)，從而實(shí)現(xiàn)磁共振圖像的并行采集重建。該方法能夠在一方面加速成像采集速度的同時(shí)，另一方面還可以提高重建圖像的速度。圖4是本發(fā)明實(shí)施方式提供的磁共振成像方法的流程示意圖。如圖4所示，所述磁共振成像方法包括：步驟S101，通過第一變換操作將欠采集的K空間域變換至中間域，所述中間域位于K空間域與圖像域之間；步驟S102，分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算；步驟S103，在所述中間域中，基于第二應(yīng)采集數(shù)據(jù)以及所述線圈合并系數(shù)進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)；所述第二應(yīng)采集數(shù)據(jù)由對(duì)在K空間域中采集的第一應(yīng)采集數(shù)據(jù)進(jìn)行所述第一變換操作后得到；步驟S104，通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域變換至圖像域以獲得圖像。下面以具體實(shí)施例對(duì)上述磁共振成像方法作詳細(xì)說明。在具體實(shí)施時(shí)，需要先采集K空間域的數(shù)據(jù)，所采集的數(shù)據(jù)包括應(yīng)采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)，為了區(qū)別在中間域中對(duì)應(yīng)的應(yīng)采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)，因此將在K空間域中采集的應(yīng)采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)分別稱為第一應(yīng)采集數(shù)據(jù)和第一校準(zhǔn)數(shù)據(jù)。在K空間域采集到所述第一應(yīng)采集數(shù)據(jù)和第一校準(zhǔn)數(shù)據(jù)之后，執(zhí)行步驟S101，通過第一變換操作將欠采集的K空間域變換至中間域，所述中間域位于K空間域與圖像域之間。本實(shí)施例中，步驟S101中所述的“欠采集的K空間域”是指采集到所述第一應(yīng)采集數(shù)據(jù)和第一校準(zhǔn)數(shù)據(jù)之后形成的K空間域，該K空間域既區(qū)別于“未采集的K空間域”，也區(qū)別于“完成填補(bǔ)的K空間域”或者“并行采集重建后的K空間域”。因此，所述欠采集的K空間域中的數(shù)據(jù)包括在K空間域中采集的第一應(yīng)采集數(shù)據(jù)和第一校準(zhǔn)數(shù)據(jù)。本實(shí)施例中，所述第一變換操作具體為頻率編碼方向的一維傅里葉變換，所述中間域具體為對(duì)所述欠采集的K空間域進(jìn)行一維傅里葉變換之后得到的數(shù)據(jù)空間。欠采集的K空間域變換至中間域的過程，具體是對(duì)所述欠采集的K空間域中的數(shù)據(jù)進(jìn)行所述第一變換操作的過程。本發(fā)明實(shí)施方式中，將對(duì)所述第一應(yīng)采集數(shù)據(jù)進(jìn)行第一變換操作之后得到的數(shù)據(jù)稱為第二應(yīng)采集數(shù)據(jù)，將對(duì)所述第一校準(zhǔn)數(shù)據(jù)進(jìn)行第一變換操作之后得到的數(shù)據(jù)稱為第二校準(zhǔn)數(shù)據(jù)。因此，所述中間域中的數(shù)據(jù)包括所述第二應(yīng)采集數(shù)據(jù)和第二校準(zhǔn)數(shù)據(jù)。執(zhí)行步驟S102，分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算。本實(shí)施例中，步驟S102中所述的分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算是在所述中間域進(jìn)行的，此時(shí)，步驟S102在步驟S101之后執(zhí)行，步驟S102中所述的校準(zhǔn)數(shù)據(jù)具體為所述第二校準(zhǔn)數(shù)據(jù)。在其他實(shí)施例中，所述分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算也可以在K空間域進(jìn)行，此時(shí)，步驟S102在步驟S101之前執(zhí)行，步驟S102中所述的校準(zhǔn)數(shù)據(jù)具體為所述第一校準(zhǔn)數(shù)據(jù)。若所述分離出校準(zhǔn)數(shù)據(jù)以及計(jì)算線圈合并系數(shù)的過程在K空間域進(jìn)行，則可以采用現(xiàn)有的實(shí)現(xiàn)方式進(jìn)行，此處不再詳細(xì)描述；若所述分離出校準(zhǔn)數(shù)據(jù)以及計(jì)算線圈合并系數(shù)的過程在所述中間域進(jìn)行，其具體實(shí)現(xiàn)仍然可以參照在K空間域的實(shí)現(xiàn)方式進(jìn)行，因?yàn)榘l(fā)明人發(fā)現(xiàn)，雖然經(jīng)過所述第一變換操作，應(yīng)采集數(shù)據(jù)和校準(zhǔn)數(shù)據(jù)的數(shù)據(jù)內(nèi)容發(fā)生了變化，但第二應(yīng)采集數(shù)據(jù)和第二校準(zhǔn)數(shù)據(jù)之間的對(duì)應(yīng)關(guān)系依然保持不變，據(jù)此擬合出的線圈合并系數(shù)也與在K空間域中所計(jì)算出的線圈合并系數(shù)是相同的，這樣通過步驟S102計(jì)算出線圈合并系數(shù)后，在后續(xù)步驟（步驟S103）便能夠基于該線圈合并系數(shù)以及第二應(yīng)采集數(shù)據(jù)在中間域進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)，即在中間域?qū)η凡杉瘮?shù)據(jù)進(jìn)行填補(bǔ)是可以實(shí)施的。本實(shí)施例中，無論所述數(shù)據(jù)線性加權(quán)加和操作中是否引入頻率編碼方向的數(shù)據(jù)（即表明所述數(shù)據(jù)線性加權(quán)加和操作是否受到頻率編碼方向的影響），對(duì)于在中間域計(jì)算出的線圈合并系數(shù)并無影響，均與在K空間域中計(jì)算出的線圈合并系數(shù)相同，下面將對(duì)此進(jìn)行詳細(xì)說明。1）不引入頻率編碼方向?qū)?shù)據(jù)線性加權(quán)加和操作的影響，有：設(shè)計(jì)算的基矩陣（對(duì)應(yīng)于圖1中的黑色數(shù)據(jù)點(diǎn)）為Kb，校準(zhǔn)線矩陣（對(duì)應(yīng)于圖1中的灰色數(shù)據(jù)點(diǎn)）為L，線圈合并系數(shù)為N，則有關(guān)系：Kb·N=L（1）則經(jīng)過頻率編碼方向的一維傅里葉變換后，進(jìn)入另一線性空間（中間域）其基矩陣變換為Fb，即：Fb=F·Kb（2）校準(zhǔn)線矩陣變換為Fl，即：Fl=F·L（3）其中，F(xiàn)為傅里葉變換矩陣。在式子（1）兩邊左乘F，可得：Fb·N=Fl（4）所以在中間域里線圈合并系數(shù)N不變，如果式子（4）對(duì)應(yīng)的方程為超定的，則可以求得：N=(FbH·Fb)-1FbH·Fl（5）利用式子（5）得到的線圈合并系數(shù)N，可以在中間域?qū)η凡杉瘮?shù)據(jù)進(jìn)行填補(bǔ)，即填補(bǔ)如圖1所示的白色空缺點(diǎn)。2）引入頻率編碼方向?qū)?shù)據(jù)線性加權(quán)加和操作的影響，有：圖5是頻率編碼方向影響下K空間域的數(shù)據(jù)分布示意圖。引入頻率編碼方向?qū)?shù)據(jù)線性加權(quán)加和操作的影響，以圖5為例說明。引入頻率編碼方向的影響，在如圖5所示的K空間域中，將會(huì)額外多出兩個(gè)基矩陣，為Kb1和Kb2，分別為Kb向上向下平移一個(gè)ΔRo，ΔRo表示頻率編碼方向的步長。若記：則有：Kb1=T1·Kb（6）Kb2=T2·Kb（7）則在K空間域存在關(guān)系：Kb1·N1+Kb2·N2+Kb·N3=L（8）可以將式子（8）轉(zhuǎn)為：因?yàn)椋涸賹⑹阶樱?）、（7）代入式子（9），可以得到：(T1T2E)·Kb·N=L（10）經(jīng)過頻率編碼方向的一維傅里葉變換之后進(jìn)入中間域，（2）、（3）式仍然成立，對(duì)于（10）式兩邊左乘F，有：(F·T1F·T2F)·Kb·N=F·L（11）由于F為滿秩矩陣，所以存在F-1，則記：M1=F·T1·F-1（12）M2=F·T2·F-1（13）可知：F·T1=M1·F（14）F·T2=M2·F（15）將式子（14）、（15）代入式子（11），可得：(M1·FM2·FF)·Kb·N=F·L（16）可以將式子（16）轉(zhuǎn)為：(M1M2E)·F·Kb·N=F·L（17）將式子（2）、（3）代入式子（17），可得：(M1M2E)·Fb·N=Fl（18）式子（18）中，M1，M2，E，F(xiàn)b，F(xiàn)l都為已知量，由此可知，線圈合并系數(shù)N相對(duì)于K空間域不變。記：Fβ=(M1M2E)·Fb（19）則有：N=(FβH·Fβ)-1·FβH·Fl（20）利用式子（20）得到的線圈合并系數(shù)N，可以在中間域?qū)η凡杉瘮?shù)據(jù)進(jìn)行填補(bǔ)，即填補(bǔ)如圖5所示的白色空缺點(diǎn)。由于傅里葉變換以及T1、T2矩陣的一些特性，M1和M2都為對(duì)角陣，若記：則有：在實(shí)際操作中，M1·Fb的操作，可以簡化為Fb第一行至最后一行分別乘一個(gè)常數(shù)（對(duì)應(yīng)于M1對(duì)角線上的數(shù)值，如第j行則乘以Ej）。將K空間域轉(zhuǎn)換至中間域以及計(jì)算出線圈合并系數(shù)之后，執(zhí)行步驟S103，在所述中間域中，基于第二應(yīng)采集數(shù)據(jù)以及所述線圈合并系數(shù)進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)；所述第二應(yīng)采集數(shù)據(jù)由對(duì)在K空間域中采集的第一應(yīng)采集數(shù)據(jù)進(jìn)行所述第一變換操作后得到。如前所述，無論是否受到頻率編碼方向的影響，均不影響步驟S103中所述的數(shù)據(jù)線性加權(quán)加和操作的結(jié)果，所述數(shù)據(jù)線性加權(quán)加和操作的具體實(shí)施可以參照公式（4）及公式（18）以及現(xiàn)有的在K空間域中GRAPPA算法填補(bǔ)操作的實(shí)現(xiàn)方式，此處不再詳細(xì)描述。需要說明的是，現(xiàn)有技術(shù)中，對(duì)于欠采集數(shù)據(jù)的填補(bǔ)在K空間域進(jìn)行，線圈合并系數(shù)N可以通過下式得到：N=(KbH·Kb)-1KbH·Fl(21)由此可知，與在中間域計(jì)算出的線圈合并系數(shù)是相同的。此外，通過步驟S103對(duì)欠采集數(shù)據(jù)進(jìn)行填補(bǔ)后數(shù)據(jù)，與現(xiàn)有技術(shù)在K空間域?qū)η凡杉瘮?shù)據(jù)進(jìn)行填補(bǔ)后數(shù)據(jù)有所區(qū)別，前者為中間域中的數(shù)據(jù)，而后者為K空間域中的數(shù)據(jù)，對(duì)前者進(jìn)行所述第一變換操作可以得到后者。完成數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)對(duì)欠采集數(shù)據(jù)的填補(bǔ)之后，執(zhí)行步驟S104，通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域變換至圖像域以獲得圖像。本實(shí)施例中，步驟S104中所述的第二變換操作具體為相位編碼方向的一維傅里葉變換。完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域中的數(shù)據(jù)，再經(jīng)過一維傅里葉變換，即可以變換至圖像域，實(shí)現(xiàn)磁共振成像。對(duì)應(yīng)于上述磁共振成像方法，本發(fā)明實(shí)施方式還提供一種磁共振成像裝置。圖6是本發(fā)明實(shí)施方式提供的磁共振成像裝置的結(jié)構(gòu)示意圖。如圖6所示，所述磁共振成像裝置包括：第一變換單元101，適于通過第一變換操作將欠采集的K空間域變換至中間域，所述中間域位于K空間域與圖像域之間；計(jì)算單元102，適于分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算；填補(bǔ)單元103，適于在所述中間域中，基于第二應(yīng)采集數(shù)據(jù)以及所述線圈合并系數(shù)進(jìn)行數(shù)據(jù)線性加權(quán)加和操作，以實(shí)現(xiàn)欠采集數(shù)據(jù)的填補(bǔ)；所述第二應(yīng)采集數(shù)據(jù)由對(duì)在K空間域中采集的第一應(yīng)采集數(shù)據(jù)進(jìn)行所述第一變換操作后得到；第二變換單元104，適于通過第二變換操作將完成欠采集數(shù)據(jù)填補(bǔ)之后的中間域變換到圖像域以獲得圖像。具體實(shí)施時(shí)，所述計(jì)算單元102分離出校準(zhǔn)數(shù)據(jù)，并進(jìn)行線圈合并系數(shù)的計(jì)算可以在所述中間域進(jìn)行，所述校準(zhǔn)數(shù)據(jù)為對(duì)在K空間域中采集的第一校準(zhǔn)數(shù)據(jù)進(jìn)行所述第一變換操作后得到的第二校準(zhǔn)數(shù)據(jù)。此外，所述計(jì)算單元102分離出校準(zhǔn)數(shù)據(jù)，以進(jìn)行線圈合并系數(shù)的計(jì)算也可以在K空間域進(jìn)行，所述校準(zhǔn)數(shù)據(jù)為在K空間域中采集的第一校準(zhǔn)數(shù)據(jù)。所述數(shù)據(jù)線性加權(quán)加和操作可以在未受到頻率編碼方向的影響下進(jìn)行，也可以在受到頻率編碼方向的影響下進(jìn)行，在上述兩種情況下，數(shù)據(jù)線性加權(quán)加和操作所依據(jù)的線圈合并系數(shù)是相同的。在具體實(shí)施例中，所述第一變換單元101進(jìn)行的第一變換操作為頻率編碼方向的一維傅里葉變換，所述第二變換單元104進(jìn)行的第二變換操作為相位編碼方向的一維傅里葉變換。本發(fā)明實(shí)施方式提供的磁共振成像裝置的具體實(shí)施可參考上述磁共振成像方法的實(shí)施，在此不再贅述。本領(lǐng)域技術(shù)人員可以理解，實(shí)現(xiàn)上述實(shí)施例中的磁共振成像裝置的全部或部分是可以通過程序來指令相關(guān)的硬件來完成，所述的程序可以存儲(chǔ)于計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中，所述存儲(chǔ)介質(zhì)可以是ROM、RAM、磁碟、光盤等。本發(fā)明雖然已以較佳實(shí)施例公開如上，但其并不是用來限定本發(fā)明，任何本領(lǐng)域技術(shù)人員在不脫離本發(fā)明的精神和范圍內(nèi)，都可以利用上述揭示的方法和技術(shù)內(nèi)容對(duì)本發(fā)明技術(shù)方案做出可能的變動(dòng)和修改，因此，凡是未脫離本發(fā)明技術(shù)方案的內(nèi)容，依據(jù)本發(fā)明的技術(shù)實(shí)質(zhì)對(duì)以上實(shí)施例所作的任何簡單修改、等同變化及修飾，均屬于本發(fā)明技術(shù)方案的保護(hù)范圍。