專利名稱:磁共振成像裝置的制作方法
技術(shù)領域:
本發(fā)明的實施方式涉及磁共振成像(MRI =Magnetic Resonance Imaging)裝置。
背景技術(shù):
MRI是如下的攝像法利用拉莫爾頻率的高頻(RF :radio frequency)信號對置于靜磁場中的被檢體的原子核自旋進行磁激勵,并根據(jù)伴隨該激勵產(chǎn)生的核磁共振(NMR nuclear magnetic resonance) fW^SlilSl^o在該磁共振成像的領域中,公知有半傅里葉成像(HFI :half Fourier imaging) 法。HFI法是如下方法利用k空間中的數(shù)據(jù)的復對稱性,根據(jù)所收集的k空間數(shù)據(jù),對在k 空間上未收集數(shù)據(jù)的區(qū)域的數(shù)據(jù)進行填補從而圖像化。在HFI的現(xiàn)有的數(shù)據(jù)處理中,執(zhí)行根據(jù)所收集的k空間數(shù)據(jù)對在k空間上未被收集的數(shù)據(jù)進行填補的HFI處理。由此,生成整個空間的數(shù)據(jù)。HFI按照以下的步驟執(zhí)行。這里,設k空間的kz方向的高頻區(qū)域是k空間數(shù)據(jù)的非采樣區(qū)域。首先,對k空間的kz方向進行Homodyne濾波(Homodyne filter 零差濾波)fh和低通濾波(LPF:low pass filter) f 1這兩種濾波處理。Homodyne濾波fh是實施與對非采樣區(qū)域填補復共軛數(shù)據(jù)的處理等價的濾波處理的濾波。即,如式(1-1)和式(1- 所示,對 k空間數(shù)據(jù)K分別執(zhí)行Homodyne濾波處理fh和低通濾波處理Π,分別生成濾波處理后的 k 空間數(shù)據(jù) Kh (kx,ky, kz)、Kl (kx,ky, kz)。K (kx, ky, kz) fh — Kh (kx, ky, kz) (1-1)K (kx, ky, kz) fl — Kl (kx, ky, kz) (1-2)接著,進行三維(3Dthree dimensional)的快速傅里葉變換(FFT :fast Fourier Transform),使其成為r空間(實空間)的數(shù)據(jù)。S卩,如式(2_1)和式(2_2)所示,對濾波處理后的k空間數(shù)據(jù)Kh (kx, ky, kz)、Kl (kx, ky, kz)分別進行3D FFT,從而生成r空間數(shù)
據(jù)乂??!“,又,乙”論,又,乙)。FFT {Kh (kx, ky, kz)} — Vh (χ, y, ζ) (2-1)
FFT {ΚΙ (kx, ky, kz)} — Vl (χ, y, ζ) (2-2)接著,對Homodyne濾波處理fh后的r空間數(shù)據(jù)Vh進行相位校正處理和用于去除誤差量的REAL化處理。如式C3)所示,執(zhí)行相位校正處理和REAL化處理,生成相位校正后的r空間數(shù)據(jù)V(x,y,z)。V (χ, y, ζ) = REAL{Vh (χ, y, z)*V]*(x, y, z)/|Vl(x, y, ζ) |} (3)這里,νΓ表示Vl的復共軛(complex conjugation),REAL ()是取實部的函數(shù)。另外,如式(4)所示,使用加權(quán)系數(shù)α對將r空間數(shù)據(jù)V(x,y, ζ)進行3D快速逆傅里葉變換 (IFFT Jnverse Fast FourierTransform)得到的 k 空間數(shù)據(jù) IFFT{V(x,y,ζ)}和原始的 k 空間數(shù)據(jù)K(kx,ky, kz)進行加權(quán)運算。K(kx,ky,kz) = α IFFT {V (χ, y, ζ)} + (1- α ) K (kx, ky, kz) (4)
進而,加權(quán)運算的結(jié)果為用于再次計算r空間數(shù)據(jù)V(x,y,z)的k空間數(shù)據(jù)K(kx, ky,kz),為了提高k空間中的填補處理的精度,反復執(zhí)行1 4次左右式(1-1)、式(1-2)、 式0-1)、式0-2)、式(3)以及式(4)所示的HFI處理。即,反復進行濾波處理、3D FFT、相位校正處理以及加權(quán)運算處理。但是,在三維的HFI中,與二維的HFI相比,收集的數(shù)據(jù)量大,因此,重建所需的時間長。為了提高k空間中的填補處理的精度,需要反復執(zhí)行HFI處理,因此,重建時間進一步延長與HFI處理的反復次數(shù)相應的時間,數(shù)據(jù)處理時間增大。
發(fā)明內(nèi)容
本發(fā)明的實施方式是鑒于上述課題而完成的,其目的在于,縮短三維的HFI中的數(shù)據(jù)處理時間。本發(fā)明的實施方式提供一種磁共振成像裝置,該磁共振成像裝置具備數(shù)據(jù)收集單元,根據(jù)半傅里葉成像用的攝像條件,收集三維k空間中的磁共振信號;相位校正量計算單元,將比上述三維k空間中的磁共振信號少的磁共振信號作為用于求出第一相位校正量的k空間數(shù)據(jù),對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括基于用于求出上述第一相位校正量的k空間數(shù)據(jù)的相位校正、以及利用了 k空間的復對稱性對非采樣區(qū)域進行的數(shù)據(jù)填補,由此,求出上述第一相位校正量;以及圖像數(shù)據(jù)生成單元,對上述三維 k空間中的磁共振信號執(zhí)行如下處理該處理包括使用了基于上述第一相位校正量的第二相位校正量的相位校正、上述數(shù)據(jù)填補以及圖像重建處理,由此,生成磁共振圖像數(shù)據(jù)。
圖1是表示本發(fā)明的實施方式的磁共振成像裝置的構(gòu)成圖。圖2是圖1所示的計算機的功能框圖。圖3是表示在圖2所示的第一相位校正量計算部中執(zhí)行的第一 HFI處理的一個例子的圖。圖4是表示在圖2所示的HFI處理部中使用的Homodyne濾波特性的一個例子的圖。圖5是表示在圖2所示的HFI處理部中使用的LPF特性的一個例子的圖。圖6是表示在圖2所示的第一相位校正量計算部中執(zhí)行的第一 HFI處理的另一個例子的圖。圖7是示出在圖1所示的磁共振成像裝置中,表示通過HFI獲得被檢體的圖像數(shù)據(jù)時的流程的流程圖。
具體實施例方式本發(fā)明的實施方式的磁共振成像裝置具備數(shù)據(jù)收集單元、相位校正量計算單元以及圖像數(shù)據(jù)生成單元。數(shù)據(jù)收集單元根據(jù)半傅里葉成像用的攝像條件收集三維k空間中的磁共振信號。相位校正量計算單元將比上述三維k空間中的磁共振信號少的磁共振信號作為用于求出第一相位校正量的k空間數(shù)據(jù),對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括基于用于求出第一相位校正量的k空間數(shù)據(jù)的相位校正、以及利用k空間的復對稱性對非采樣區(qū)域的數(shù)據(jù)填補,由此,求出上述第一相位校正量。圖像數(shù)據(jù)生成單元對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括使用基于上述第一相位校正量的第二相位校正量的相位校正、上述數(shù)據(jù)填補以及圖像重建處理,由此生成磁共振圖像數(shù)據(jù)。參照
本發(fā)明的實施方式的磁共振成像裝置。圖1是表示本發(fā)明的實施方式的磁共振成像裝置的構(gòu)成圖。磁共振成像裝置20具備形成靜磁場的筒狀的靜磁場用磁鐵21 ;設于靜磁場用磁鐵21的內(nèi)部的勻磁線圈22、傾斜磁場線圈23以及RF線圈M。并且,磁共振成像裝置20具備控制系統(tǒng)25??刂葡到y(tǒng)25具備靜磁場電源沈、傾斜磁場電源27、勻磁線圈電源觀、發(fā)送器四、接收器30、順序控制器31以及計算機32??刂葡到y(tǒng)25的傾斜磁場電源27由X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y和Z軸傾斜磁場電源27z構(gòu)成。并且,計算機32具備輸入裝置33、顯示裝置34、運算裝置35以及存儲裝置36。靜磁場用磁鐵21與靜磁場電源沈連接,通過從靜磁場電源沈供給的電流而在攝像區(qū)域形成靜磁場。靜磁場用磁鐵21例如由超導線圈構(gòu)成。靜磁場用磁鐵21在勵磁時與靜磁場電源沈連接而被供給電流,但一旦被勵磁后,通常就成為非連接狀態(tài)。并且,有時也利用永久磁鐵構(gòu)成靜磁場用磁鐵21,而不設置靜磁場電源26。并且,在靜磁場用磁鐵21的內(nèi)側(cè)同軸地設有筒狀的勻磁線圈22。勻磁線圈22被構(gòu)成為與勻磁線圈電源28連接,從勻磁線圈電源28向勻磁線圈22供給電流,從而使靜磁場均勻化。傾斜磁場線圈23由X軸傾斜磁場線圈23x、Y軸傾斜磁場線圈23y和Z軸傾斜磁場線圈23z構(gòu)成,在靜磁場用磁鐵21的內(nèi)部形成為筒狀。在傾斜磁場線圈23的內(nèi)側(cè)設有臥床37作為攝像區(qū)域,被檢體P置于臥床37上。在RF線圈M具有內(nèi)置于雷達天線的RF 信號收發(fā)用的全身用線圈(WBC :whole body coil)、設于臥床37或被檢體P附近的RF信號的接收用的局部線圈等。并且,傾斜磁場線圈23與傾斜磁場電源27連接。傾斜磁場線圈23的X軸傾斜磁場線圈23x、Y軸傾斜磁場線圈23y和Z軸傾斜磁場線圈23z分別與傾斜磁場電源27的X 軸傾斜磁場電源27x、Y軸傾斜磁場電源27y和Z軸傾斜磁場電源27z連接。進而,從X軸傾斜磁場電源27x、Y軸傾斜磁場電源27y和Z軸傾斜磁場電源27z 分別向X軸傾斜磁場線圈23x、Y軸傾斜磁場線圈23y和Z軸傾斜磁場線圈23z供給電流。 由此,在攝像區(qū)域分別形成X軸方向的傾斜磁場Gx、Y軸方向的傾斜磁場Gy、Z軸方向的傾斜磁場(iz。RF線圈M與發(fā)送器四和接收器30中的至少一方連接。接收用的RF線圈M從發(fā)送器四接收RF信號并發(fā)送給被檢體P。接收用的RF線圈M接收伴隨被檢體P內(nèi)部的原子核自旋的RF信號的激勵產(chǎn)生的MR信號并將其賦予接收器30。另一方面,控制系統(tǒng)25的順序控制器31與傾斜磁場電源27、發(fā)送器四以及接收器30連接。順序控制器31存儲驅(qū)動傾斜磁場電源27、發(fā)送器四以及接收器30所需的控制信息,例如記述了要對傾斜磁場電源27施加的脈沖電流強度、施加時間、施加定時等動作控制信息的順序信息。并且,順序控制器31按照存儲的順序信息來驅(qū)動傾斜磁場電源27、發(fā)送器四以及接收器30,由此,產(chǎn)生X軸傾斜磁場Gx、Y軸傾斜磁場Gy、Z軸傾斜磁場 Gz以及RF信號。并且,順序控制器31構(gòu)成為,接收通過接收器30中的MR信號的檢波和A/ D(analog todigital)轉(zhuǎn)換得到的復數(shù)據(jù)即原始數(shù)據(jù)(raw data)并賦予給計算機32。發(fā)送器四根據(jù)從順序控制器31接收到的控制信息將RF信號賦予給RF線圈24。 接收器30對從RF線圈M接收到的MR信號進行檢波并執(zhí)行信號處理,并且,通過對其進行 A/D轉(zhuǎn)換,從而生成被數(shù)據(jù)化了的復數(shù)據(jù)即原始數(shù)據(jù),并將生成的原始數(shù)據(jù)賦予給順序控制器31。在RF線圈M由多個線圈要素構(gòu)成的情況下,既可以將從各線圈要素輸出的MR信號分別輸入到接收通道進行處理,也可以將從多個線圈要素中兩個以上的線圈要素輸出的MR 信號合成而輸入到接收通道進行處理。并且,利用運算裝置35執(zhí)行在計算機32的存儲裝置36中存儲的程序,由此,能夠在計算機32中實現(xiàn)各種功能。但是,不限于程序,也可以在磁共振成像裝置20中設置具有各種功能的電路。圖2是圖1所示的計算機32的功能框圖。計算機32中,通過運算裝置35執(zhí)行程序,由此,能夠作為攝像條件設定部40、順序控制器控制部41、k空間數(shù)據(jù)庫42、HFI處理部43以及圖像處理部44發(fā)揮功能。攝像條件設定部40根據(jù)來自輸入裝置33的指示信息,設定3D HFI用的攝像條件,將設定的攝像條件賦予給順序控制器控制部41。3D-HFI法是如下方法在3D的k空間中以在相位編碼(PE:phase encode)方向或分層編碼(SE :slice encode)方向非對稱的方式收集MR數(shù)據(jù),利用k空間中的數(shù)據(jù)的復對稱性,根據(jù)收集到的MR數(shù)據(jù)(采樣區(qū)域中的k空間數(shù)據(jù)),對未收集MR數(shù)據(jù)的非采樣區(qū)域的k空間數(shù)據(jù)進行填補,由此圖像化。因此,作為3D-HFI法用的攝像條件,設定與PE方向或SE方向上的正側(cè)或負側(cè)的高頻區(qū)域?qū)腗R數(shù)據(jù)未被收集的3D數(shù)據(jù)收集條件。另外,作為HFI用的脈沖序列,具有FASE (fast asymmetric spin echo對禾爾自旋回波、或fast advanced spin echo ‘決速高級自旋回波)序列等。順序控制器控制部41根據(jù)來自輸入裝置33的掃描開始指示信息,對順序控制器 31賦予包括脈沖序列的攝像條件,由此進行驅(qū)動控制。并且,順序控制器控制部41從順序控制器31接收原始數(shù)據(jù)而將其配置于在k空間數(shù)據(jù)庫42形成的k空間。因此,在根據(jù)3D-HFI用的攝像條件收集到三維k空間中的MR信號的情況下,不存在PE方向或SE方向的正側(cè)或負側(cè)的高頻區(qū)域的MR數(shù)據(jù)的k空間數(shù)據(jù)將被配置于在k空間數(shù)據(jù)庫42形成的k空間。HFI處理部43具有如下功能根據(jù)從k空間數(shù)據(jù)庫42獲得的k空間數(shù)據(jù),在半傅立葉QiaIf Fourier)法中執(zhí)行生成最終圖像數(shù)據(jù)所需的處理(HFI處理)。HFI處理是這樣的處理包括用于根據(jù)采樣區(qū)域中的k空間數(shù)據(jù)對k空間的非采樣區(qū)域中的k空間數(shù)據(jù)進行填補的濾波處理和相位校正處理。另外,相位校正處理是為了確保受相位分布的影響而實質(zhì)上不成立的復共軛的對稱性而實施的。用于算出相位校正量的相位分布能夠根據(jù)采樣區(qū)域中的全部或一部分k空間數(shù)據(jù)來推定。特別是從高精度地推定相位分布的觀點出發(fā),優(yōu)選根據(jù)k空間數(shù)據(jù)對稱采樣得到的低頻區(qū)域中的k空間數(shù)據(jù)推定相位分布。在該情況下,在相位校正之前先實施用于根據(jù)采樣區(qū)域中的k空間數(shù)據(jù)抽出低頻區(qū)域的k空間數(shù)據(jù)的LPF處理。但是,即便進行相位校正處理,也存在相位校正的錯誤、即復數(shù)信號的虛部也不是零的情況。因此,通過將由數(shù)據(jù)填補前后的k空間數(shù)據(jù)的加權(quán)運算得到的k空間數(shù)據(jù)作為初始數(shù)據(jù)反復進行數(shù)據(jù)的填補處理和相位校正處理的循環(huán)處理,能夠降低相位校正錯誤并提高數(shù)據(jù)填補的精度。除此之外,除了低頻區(qū)域的k空間數(shù)據(jù)外,有時通過使用對k空間數(shù)據(jù)進行非對稱采樣得到的高頻區(qū)域的k空間數(shù)據(jù)來推定相位分布,也能夠降低相位校正的錯誤。在該情況下,可以代替LPF,在相位校正之前先執(zhí)行用于根據(jù)采樣區(qū)域中的k空間數(shù)據(jù)抽出高頻區(qū)域和低頻區(qū)域雙方的k空間數(shù)據(jù)的濾波處理。另外,有時通過在非采樣區(qū)域的數(shù)據(jù)的填補前實施相位校正,也能夠降低相位校正錯誤。并且,還知道有循環(huán)處理中不包括數(shù)據(jù)填補用的濾波處理,在數(shù)據(jù)填補后反復進行包括相位校正和加權(quán)運算的循環(huán)處理的HFI處理。這樣,在HFI處理中存在如下各種變更為了推定相位校正用的相位分布而僅使用低頻區(qū)域的k空間數(shù)據(jù),或者還使用高頻區(qū)域的k空間數(shù)據(jù),或者將相位校正作為非采樣區(qū)域的數(shù)據(jù)的填補處理的前處理進行或作為后處理進行,或者是否進行相位校正錯誤的降低用的循環(huán)處理,或者在循環(huán)處理中包括數(shù)據(jù)填補用的濾波處理等各種變更。以下,作為一個例子,說明如下情況僅使用低頻區(qū)域的k空間數(shù)據(jù)推定相位分布,并且將相位校正作為非采樣區(qū)域的數(shù)據(jù)的填補處理的后處理,反復進行包括數(shù)據(jù)填補處理和相位校正的循環(huán)處理。當然,對于基于其他方法的HFI處理也可以同樣應用以后的處理。在該情況下,HFI處理通過反復進行預定次數(shù)的如下多個處理,即,用于進行k空間的非采樣區(qū)域中的數(shù)據(jù)填補的濾波處理、用于進行相位分布的推定用的低頻區(qū)域中的數(shù)據(jù)抽出的濾波處理、對實空間數(shù)據(jù)的相位校正處理、取相位校正后的實空間數(shù)據(jù)的實部的 REAL化處理、k空間數(shù)據(jù)的加權(quán)運算處理等,從而成為使實空間數(shù)據(jù)的虛部收斂于零的處理。因此,在HFI處理部43中具備如下功能除了上述處理外,還執(zhí)行構(gòu)成HFI處理的各種處理,例如將k空間數(shù)據(jù)轉(zhuǎn)換為實空間數(shù)據(jù)的FT處理或FFT處理、或?qū)⒆鳛閷嵖臻g數(shù)據(jù)的圖像數(shù)據(jù)轉(zhuǎn)換為k空間數(shù)據(jù)的IFT處理或IFFT處理等。特別是HFI處理部43構(gòu)成為不使用3D的采樣區(qū)域中收集到的MR數(shù)據(jù)(k空間數(shù)據(jù))的全部而是根據(jù)一部分MR數(shù)據(jù)進行相位分布的推定,由此,與以往相比能夠減少HFI 處理量及HFI處理時間。因此,HFI處理部43具有第一相位校正量計算部43A、第二相位校正量決定部43B以及診斷圖像生成部43C。第一相位校正量計算部43A具有如下功能通過抽出根據(jù)3D HFI用的攝像條件收集到的3D的k空間中的MR信號的一部分,從而將比收集到的MR信號少的MR信號設定為用于求出相位校正用的第一相位校正量的k空間數(shù)據(jù)的功能;以及通過第一 HFI處理求出第一相位校正量的功能,所述第一 HFI處理是這樣的處理對3D的k空間中的MR信號反復執(zhí)行包括利用k空間中的復對稱性向非采樣區(qū)域填補數(shù)據(jù)用的取復共軛的濾波處理、以及基于用于求出第一相位校正量的k空間數(shù)據(jù)的相位校正的處理。圖3是表示在圖2所示的第一相位校正量計算部43A中執(zhí)行的第一 HFI處理的一個例子的圖。
如圖3㈧所示,通過3D數(shù)據(jù)收集,來收集具有與實空間的PE方向、SE方向以及讀取(R0 readout)方向分別對應的k空間的kpe方向、kse方向以及kro方向的三軸方向的空間性的k空間數(shù)據(jù)。其中,在HFI中,收集在PE方向或SE方向非對稱的3D-k空間數(shù)據(jù)。圖3(A)示出將SE方向的負極側(cè)的高頻區(qū)域作為非采樣區(qū)域的例子。因此,只有SE方向的低頻區(qū)域和正極側(cè)的高頻區(qū)域成為采樣區(qū)域,在kpe ^ -Ka的區(qū)域填充k空間數(shù)據(jù)。PE方向、SE方向和RO方向分別根據(jù)診斷目的設定為χ軸、y軸、ζ軸以及傾斜方向的軸等任意的空間坐標軸。圖3㈧示出在實空間中將PE方向、SE方向和RO方向分別設定為y軸、ζ軸、χ軸的例子。因此,k空間中的kpe方向、kse方向以及kro方向分別為 ky方向、kz方向和kx方向。當對圖3 (A)所示的3D_k空間數(shù)據(jù)實施Homodyne濾波時,生成與在非采樣區(qū)域填補了 k空間數(shù)據(jù)得到的k空間數(shù)據(jù)等價的k空間數(shù)據(jù)。另外,Homodyne濾波是實施與對非采樣區(qū)域填補采樣區(qū)域的復共軛數(shù)據(jù)的處理等價的濾波處理的濾波。圖4是表示在圖2所示的HFI處理部43中使用的Homodyne濾波特性的一個例子的圖。在圖4㈧中,橫軸表示kz方向,縱軸表示Homodyne濾波的強度。并且,在圖4(B) 中,橫軸表示kz方向,縱軸表示kx方向。如圖3所示,在收集了在kz方向非對稱的3D_k空間數(shù)據(jù)的情況下,例如濾波處理可以使用具有圖4(A)所示的kz方向的強度分布、并且kx方向和ky方向的強度分布恒定的Homodyne濾波。即,濾波處理可以使用具有如下強度輪廓的Homodyne濾波在k空間數(shù)據(jù)未被收集的負極側(cè)的高頻區(qū)域-Kmax ^kz ^ -ka的范圍中強度為零、在收集了 k空間數(shù)據(jù)的低頻區(qū)域-Ka < kz < ka的范圍中強度為1. O、在收集了 k空間數(shù)據(jù)的正極側(cè)的高頻區(qū)域Ka彡kz彡Kmax的范圍中強度為2. O。圖4(A)所示的Homodyne濾波的強度分布示出了在恒定值之間直線變化的例子, 但濾波處理也可以使用以強度平滑變化的方式具有曲線形狀的強度分布的Homodyne濾波。當使用具有圖4㈧所示的強度分布的Homodyne濾波對在kz方向非對稱收集到的3D-k空間數(shù)據(jù)進行濾波處理時,得到圖4 (B)所示的Homodyne濾波處理后的k空間數(shù)據(jù)。進而,當對Homodyne濾波處理后的3D_k空間數(shù)據(jù)執(zhí)行3D-FFT時,圖3 (B)所示的多個2D-r空間數(shù)據(jù)作為分層圖像數(shù)據(jù)生成。另一方面,能夠從圖3 (A)所示的3D_k空間數(shù)據(jù)抽出一張2D_k空間數(shù)據(jù),并設定為用于求出相位校正用的第一相位校正量的k空間數(shù)據(jù)。從高精度地求出第一相位校正量的觀點出發(fā),選擇至少SNR(Signal to noise ratio)最好的編碼量為零時的2D_k空間數(shù)據(jù)來用于算出相位校正量是極為重要的。但是,為了算出第一相位校正量,需要推定非對稱方向的k空間數(shù)據(jù)的相位分布。因此,在非對稱方向上無法使編碼量為零。因此,當收集到在PE方向非對稱的MR信號時,優(yōu)選將SE方向上的編碼量為O的 2D-k空間數(shù)據(jù)作為用于求出第一相位校正量的k空間數(shù)據(jù)。另一方面,當收集到在SE方向非對稱的MR信號時,優(yōu)選將PE方向上的編碼量為O的2D-k空間數(shù)據(jù)作為用于求出第一相位校正量的k空間數(shù)據(jù)。圖3㈧示出收集了在SE方向非對稱的MR信號的例子。因此,在圖3㈧所示的例子中,PE量為0、即kpe = 0的2D-k空間數(shù)據(jù)被選擇為用于求出第一相位校正量的k空間數(shù)據(jù)。對用于求出第一相位校正量的2D_k空間數(shù)據(jù)實施LPF。由此,抽出2D_k空間數(shù)據(jù)的低頻區(qū)域中的k空間數(shù)據(jù),可以將抽出的低頻區(qū)域中的2D-k空間數(shù)據(jù)用于相位分布的推定。圖5是在圖2所示的HFI處理部43中使用的LPF特性的一個例子的圖。在圖5(A)中,橫軸表示kz方向,縱軸表示LPF的強度。并且,在圖5(B)中,橫軸表示kz方向,縱軸表示kx方向。如圖3所示,在收集了在kz方向非對稱的3D_k空間數(shù)據(jù)的情況下,濾波處理可以使用例如具有如圖5(A)所示的kz方向的強度分布、并且kx方向和ky方向的強度分布恒定的LPF。即,濾波處理可以使用具有如下強度輪廓的LPF 在收集了 k空間數(shù)據(jù)的低頻區(qū)域-Ka < kz < ka的范圍中強度為1. 0,在k空間數(shù)據(jù)未被收集的負極側(cè)的高頻區(qū)域-Kmax ^kz ^ -ka的范圍中、以及收集了 k空間數(shù)據(jù)的正極側(cè)的高頻區(qū)域Ka ^ kz ^ Kmax 的范圍中強度為零。圖5 (A)所示的LPF的強度分布示出了在恒定值之間直線變化的例子,但濾波處理也可以使用以強度平滑變化的方式具有曲線形狀的強度分布的LPF。使用具有圖5 (A)所示的強度分布的LPF進行用于求出第一相位校正量的2D_k空間數(shù)據(jù)的濾波處理時,得到圖5(B)所示的LPF處理后的k空間數(shù)據(jù)。S卩,抽出低頻區(qū)域中的2D-k空間數(shù)據(jù)用于推定相位分布。進而,當對LPF處理后的相位分布的推定用的2D_k空間數(shù)據(jù)執(zhí)行2D-FFT時,圖 3(C)所示的1張2D-r空間數(shù)據(jù)作為分層圖像數(shù)據(jù)生成。另外,通過利用相位分布的推定用的2D-r空間數(shù)據(jù)的復共軛數(shù)據(jù)除以絕對值,從而能夠算出低頻區(qū)域中的相位分布。接著,使用根據(jù)與PE = 0對應的2D_r空間數(shù)據(jù)算出的相位分布,進行Homodyne 濾波處理以及通過3D-FFT生成的多個2D-r空間數(shù)據(jù)的相位校正。其結(jié)果,能夠得到對非采用區(qū)域進行數(shù)據(jù)填補和相位校正后的多分層圖像數(shù)據(jù)。但是,當相位分布有誤差時,相位校正后的2D_r空間數(shù)據(jù)的虛部不會是零。因此, 將通過數(shù)據(jù)填補前后的k空間數(shù)據(jù)的加權(quán)運算得到的k空間數(shù)據(jù)作為初始數(shù)據(jù),反復執(zhí)行數(shù)據(jù)填補和相位校正。即,執(zhí)行使相位校正后的2D-r空間數(shù)據(jù)的虛部收斂于零的循環(huán)處理。具體而言,執(zhí)行用于通過去除相位校正后的多個2D_r空間數(shù)據(jù)的虛部而分別獲取實部從而去除誤差量的REAL化處理。接著,通過REAL化處理后的多個2D_r空間數(shù)據(jù)的2D IFFT生成數(shù)據(jù)填補后的整個頻率區(qū)域的3D k空間數(shù)據(jù)。接著,執(zhí)行數(shù)據(jù)填補前后的 3D-k空間數(shù)據(jù)的加權(quán)運算。作為該加權(quán)運算的結(jié)果得到的3D-k空間數(shù)據(jù)再次作為用于利用Homodyne濾波求出數(shù)據(jù)填補處理和第一相位校正量的與kpe = O對應的2D_k空間數(shù)據(jù)的抽出對象即初始數(shù)據(jù)使用。進而,反復進行上述循環(huán)處理預定次數(shù),該預定次數(shù)憑經(jīng)驗或通過模擬等預先決定,由此,能夠使相位校正后的r空間數(shù)據(jù)的虛部接近零。換言之,能夠?qū)⑹占降腗R數(shù)據(jù)的相位和循環(huán)處理后的數(shù)據(jù)相位的差作為適于相位校正的第一相位校正量來算出。在圖3所示的第一 HFI處理中,相位分布的推定用的數(shù)據(jù)不是使用3D_k空間數(shù)據(jù),而是使用單一的2D_k空間數(shù)據(jù)。因此,相位分布的推定所需的處理不是以往那種針對相位分布的推定用的3D-k空間數(shù)據(jù)的3D-FFT而是針對2D-k空間數(shù)據(jù)的2D-FFT。因此,能夠使算出第一相位校正量所需的數(shù)據(jù)處理量比以往減少。特別是在循環(huán)處理的次數(shù)多的情況下,由于反復進行FFT,所以數(shù)據(jù)處理量和數(shù)據(jù)處理的減少量變得顯著。另外,在使用多個線圈要素進行MR數(shù)據(jù)的收集的情況下,對與多個接收通道對應的多個3D-k空間數(shù)據(jù)分別執(zhí)行HFI處理。因此,能夠進一步減少相位分布的推定用的數(shù)據(jù)處理量和數(shù)據(jù)處理時間。圖3中示出了使用一張2D_k空間數(shù)據(jù)推定相位補正用的相位分布的例子,但也可以使用多張2D-k空間數(shù)據(jù)推定相位校正用的相位分布。在使用一張2D-k空間數(shù)據(jù)推定相位分布的情況下,推定一個作為相位校正對象的r空間數(shù)據(jù)的相位分布的代表值。與此相對,若使用多張2D_k空間數(shù)據(jù)推定相位校正用的相位分布,則可以推定作為相位校正對象的r空間數(shù)據(jù)的相位分布的多個代表值。因此,能夠進行更高精度的相位校正來生成最終的MR圖像數(shù)據(jù)。圖6是示出在圖2所示的第一相位校正量計算部43A中執(zhí)行的第一 HFI處理的另一個例子的圖。圖6所示的第一 HFI處理如(A)所示從3D_k空間數(shù)據(jù)抽出多張2D_k空間數(shù)據(jù), 并將其設定為用于求出第一相位校正量的k空間數(shù)據(jù),這點與圖3所示的第一 HFI處理不同。因此,對于相同的內(nèi)容省略說明。圖6示出將與kpe = _kpel,0,+kpel對應的三張2D_k空間數(shù)據(jù)設定為用于求出第一相位校正量的k空間數(shù)據(jù)的例子。在該情況下,抽出通過LPF得到的三張低頻區(qū)域中的2D-k空間數(shù)據(jù)。進而,通過對三張低頻區(qū)域的2D-k空間數(shù)據(jù)的3D-FFT,如(C)所示,生成三張2D-r空間數(shù)據(jù)用于推定相位分布。因此,能夠根據(jù)三張2D_r空間數(shù)據(jù),算出與相互不同的kpe方向的位置對應的三個相位分布的推定值。這相當于將(A)所示的3D的k空間在kpe方向分割為三部分,并針對每個分割區(qū)域推定相位分布。因此,作為相位校正的對象的2D-r空間數(shù)據(jù)使用根據(jù)與同一分割區(qū)域?qū)南辔环植纪贫ㄓ玫?D-r空間數(shù)據(jù)、即最接近的相位分布推定用的2D-r 空間數(shù)據(jù)推定的相位分布來進行相位校正。這樣,將在收集的MR數(shù)據(jù)對稱的編碼方向上編碼量相互不同的多個2D_k空間數(shù)據(jù)作為用于求出第一相位校正量的k空間數(shù)據(jù)的話,能夠更高精度地進行相位分布的推定和相位校正。即,通過將至少包括編碼量不是零的2D-k空間數(shù)據(jù)的多個2D-k空間數(shù)據(jù)作為用于求出第一相位校正量的k空間數(shù)據(jù),從而能夠求出與相互不同的編碼量對應的多個第一相位校正量。另一方面,通過將用于求出第一相位校正量的k空間數(shù)據(jù)作為比收集到的MR數(shù)據(jù)少的k空間數(shù)據(jù),從而能夠減少用于算出相位分布的3D-FFT的數(shù)據(jù)處理量和數(shù)據(jù)處理時間。因此,通過調(diào)整用于求出第一相位校正量的2D-k空間數(shù)據(jù)的數(shù)量和位置,從而能夠使處于折衷關(guān)系的相位校正的精度和數(shù)據(jù)處理量最佳化。第二相位校正量決定部4 具有這樣的功能根據(jù)第一相位校正量決定第二相位校正量,該第二相位校正量用于對根據(jù)3D-HFI用的攝像條件收集到的3D的k空間中的MR 信號執(zhí)行的三維相位校正。
例如第一相位校正量作為單一的值算出,并且,在收集的MR數(shù)據(jù)對稱的編碼方向上使第二相位校正量恒定的情況下,能夠?qū)⒌谝幌辔恍U恐苯記Q定為第二相位校正量。并且,在針對收集的MR數(shù)據(jù)對稱的編碼方向上不同的編碼量的每個算出多個第一相位校正量的情況下,也可以直接決定為與各編碼量對應的多個第二相位校正量。例如若是圖6所示的例子,則可以將與kpe方向的多個位置對應的多個第一相位校正量直接決定為與kpe方向的多個位置對應的第二相位校正量。另一方面,在設定與根據(jù)收集的MR數(shù)據(jù)對稱的編碼方向上的編碼量相應的多個第二相位校正量的情況下,可以基于在第一相位校正量計算部43A中算出的單一或多個第一相位校正量,求出與編碼量相應的多個第二相位校正量。例如在收集的MR數(shù)據(jù)對稱的編碼方向上,若第二相位校正量視為與編碼量成比例地呈線性地變化的話,則能夠通過使一個第一相位校正量呈線性變化、或多個第一相位校正量的插補或擬合(fitting)算出與任意的多個編碼量對應的多個第二相位校正量。例如若是圖3或圖6所示的例子,則可以算出與不同的多個kpe方向的位置對應的第二相位校正量。另外,對于用于計算第二相位校正量的系數(shù),可以通過經(jīng)驗或模擬來預先決定。診斷圖像生成部43C具有這樣的功能通過對根據(jù)3D-HFI用的攝像條件收集到的 3D的k空間中的MR信號執(zhí)行第二 HFI處理,從而生成診斷用的MR圖像數(shù)據(jù),該第二 HFI處理包括通過取三維的復共軛而對非采樣區(qū)域填補數(shù)據(jù)用的濾波處理、和使用了基于第一相位校正量的第二相位校正量的相位校正及圖像重建處理。即,在第二 HFI處理中,確定相位校正量。因此,不執(zhí)行相位分布的推定用的k空間數(shù)據(jù)的抽出和相位分布的推定。并且,未必一定要進行用于使相位校正錯誤減小的循環(huán)處理。并且,圖像處理部44具有這樣的功能對由HFI處理部43取得的HFI處理后的MR 圖像數(shù)據(jù)實施失真校正處理等必要的圖像處理,并使其顯示于顯示裝置擬。下面,對磁共振成像裝置20的動作進行說明。圖7是示出在圖1所示的磁共振成像裝置20中,表示通過HFI獲得被檢體P的圖像數(shù)據(jù)時的流程的流程圖。在步驟Sl中,攝像條件設定部40設定HFI用的攝像條件。S卩,通過FASE序列等的脈沖序列設定收集在SE方向或PE方向非對稱的3DMR數(shù)據(jù)的攝像條件。在步驟S2中,根據(jù)攝像條件設定部40設定的攝像條件進行MR數(shù)據(jù)的收集。艮口, 被檢體P設置于臥床37,在通過靜磁場電源沈勵磁的靜磁場用磁鐵21 (超導磁鐵)的攝像區(qū)域形成靜磁場。并且,從勻磁線圈電源觀向勻磁線圈22供給電流,從而使形成于攝像區(qū)域的靜磁場均勻化。進而,當從輸入裝置33對順序控制器控制部41賦予掃描開始指示時,順序控制器控制部41將從攝像條件設定部40取得的包括脈沖序列的攝像條件賦予順序控制器31。 順序控制器31根據(jù)從順序控制器控制部41賦予的脈沖序列驅(qū)動傾斜磁場電源27、發(fā)送器 29以及接收器30,從而在設置有被檢體P的攝像區(qū)域形成傾斜磁場,并且,從RF線圈M產(chǎn)生RF信號。由此通過被檢體P的內(nèi)部的磁共振產(chǎn)生的MR信號由RF線圈M接收而賦予給接收器30。接收器30從RF線圈M接收MR信號而執(zhí)行信號處理后,通過A/D轉(zhuǎn)換生成數(shù)字數(shù)據(jù)的MR信號即三維原始數(shù)據(jù)。
這里,根據(jù)HFI用的數(shù)據(jù)收集條件,以不收集三維k空間的一部分的MR數(shù)據(jù)的方式利用順序控制器控制部41進行控制。未被收集的k空間的數(shù)據(jù)通過后述的HFI處理填補。接收器30將生成的原始數(shù)據(jù)賦予給順序控制器31。順序控制器31將原始數(shù)據(jù)賦予給順序控制器控制部41,順序控制器控制部41將原始數(shù)據(jù)作為三維的k空間數(shù)據(jù)配置于在k空間數(shù)據(jù)庫42形成的k空間。在步驟S3中,HFI處理部43通過執(zhí)行第一 HFI處理來算出第一相位校正量。更具體而言,在步驟S31中,第一相位校正量計算部43A抽出在k空間數(shù)據(jù)庫42 中配置的三維的k空間數(shù)據(jù)中的用于求出第一相位校正量的單一或多個2D-k空間數(shù)據(jù)。從高精度地求出第一相位校正量的觀點出發(fā),選擇k空間數(shù)據(jù)對稱的方向的編碼量為零的2D k空間數(shù)據(jù)用于算出相位校正量較為重要。因此,這里,如圖3所示,以ky方向設定為PE方向、采樣在設定為SE方向的kz方向非對稱的k空間數(shù)據(jù)的情況為例進行說明。在該情況下,設定至少與kpe = ky = 0對應的PE方向的編碼量為零的2D_k空間數(shù)據(jù)用于算出第一相位校正量。接著,在步驟S32中,第一相位校正量計算部43A進行式(5_1)所示的Homodyne 濾波處理及式(5- 所示的LPF處理。K(kx,kz)fh — Kh(kx,kz) (5-1)KO (kx, kz) fl — Kl (kx, kz) (5-2)g卩,如式(5-1)所示,對構(gòu)成非對稱的3D_k空間數(shù)據(jù)的與Kx-Ky平面平行的多個 2D-k空間數(shù)據(jù)K(kx,kz)分別施加Homodyne濾波fh。由此,生成填補了非采樣區(qū)域中的數(shù)據(jù)的多個2D-k空間數(shù)據(jù)Kh (kx,kz)。另一方面,如式(5- 所示,對第一相位校正量計算用的單一或多個2D_k空間數(shù)據(jù)kO(kx,kz)施加LPF fl。由此,生成低頻區(qū)域中的相位分布推定用的2D-k空間數(shù)據(jù) kl (kx, kz) ο接著,在步驟S33中,如式(6-1)和式(6- 所示,第一相位校正量計算部43A對各濾波處理后的2D-k空間數(shù)據(jù)Kh (kx, kz),kl (kx, kz)實施FFT處理。由此,得到2D-R空間數(shù)據(jù) Vh (χ, ζ), Vl (χ, ζ)。FFT {Kh (kx, kz)} — Vh (χ, ζ) (6-1)FFT {Kl (kx, kz)} — Vl (x, z) (6-2)S卩,重建作為相位校正的對象的多個2D_r空間數(shù)據(jù)Vh (x,z)和用于推定相位校正用的相位分布的單一或多個2D-r空間數(shù)據(jù)Vl (χ, ζ)。接著,在步驟S34中,第一相位校正量計算部43Α使用根據(jù)式(6- 得到的2D_r 空間數(shù)據(jù)Vl (χ, ζ)算出的相位分布,進行各2D-r空間數(shù)據(jù)Vh (χ,ζ)的相位校正。接著,在步驟S35中,第一相位校正量計算部43Α判定從步驟S31到步驟S34的處理的反復次數(shù)被預先決定的次數(shù)的閾值N。閾值N憑經(jīng)驗設定為大約2-4。并且,處理的反復次數(shù)N也可以作為數(shù)據(jù)處理條件通過輸入裝置33的操作可變設定。在判定為循環(huán)處理的反復次數(shù)未達到閾值N的情況下,在步驟S36中,第一相位校正量計算部43Α進行取相位校正后的各2D-r空間數(shù)據(jù)的實部的REAL化處理。其結(jié)果,得到數(shù)據(jù)填補和相位校正后的只有實部信號的2D-r空間數(shù)據(jù)Vcor (χ,ζ)。式(7)表示該步驟S34和步驟S36的處理。Vcor (χ, ζ) = REAL {Vh (χ, ζ) *V] * (χ,ζ) / | Vl (χ, ζ) |} (7)其中,在式(7)中REALO示出取實部的處理,νΓ表示Vl的復共軛。另外,為了提高相位校正和k空間中的數(shù)據(jù)填補處理的精度,對通過上述處理生成的數(shù)據(jù)填補和相位校正后的2D-k空間數(shù)據(jù)、與數(shù)據(jù)填補和相位校正前的2D-k空間數(shù)據(jù)進行加權(quán)運算。進而,將通過加權(quán)運算得到的2D-k空間數(shù)據(jù)作為初始數(shù)據(jù),反復進行上述的2D-k空間數(shù)據(jù)的濾波處理、濾波處理后的k空間數(shù)據(jù)的FFT處理、相位校正及加權(quán)運算。因此,在步驟S37中,如式⑶所示,第一相位校正量計算部43A通過IFFT將數(shù)據(jù)填補和相位校正后的只有實部信號的2D-r空間數(shù)據(jù)Vcor (χ,ζ)轉(zhuǎn)換為2D_k空間數(shù)據(jù) Kcor (kx, kz)。Kcor (kx, kz) = IFFT {Vcor (x, z)} (8)接著,在步驟S38中,第一相位校正量計算部43A對數(shù)據(jù)填補和相位校正后的各 2D-k空間數(shù)據(jù)Kcor (kx,kz)、以及對應的數(shù)據(jù)填補和相位校正前的各2D_k空間數(shù)據(jù)K(kx, kz)進行加權(quán)運算。該加權(quán)運算的結(jié)果如式(9)所示,再次成為步驟S32中的濾波處理用的初始數(shù)據(jù)。K (kx, kz) = α Kcor (kx, kz) + (1- α ) K (kx, kz) (9)其中,式(9)中α是取O < α < 1的加權(quán)系數(shù)。加權(quán)系數(shù)α可以通過經(jīng)驗或模擬決定為適當?shù)闹?,典型的是,設定為與k空間數(shù)據(jù)的非對稱方向的位置kz對應的值。具體而言,設定加權(quán)系數(shù)α (kz),以使得在非采樣區(qū)域中數(shù)據(jù)填補和相位校正后的2D_k空間數(shù)據(jù)Kc0r(kX,kz)的比例大,在采樣區(qū)域中,數(shù)據(jù)填補和相位校正前的2D_k空間數(shù)據(jù)K (kx, kz)的比例大。并且,加權(quán)運算的加權(quán)系數(shù)α還可以作為數(shù)據(jù)處理條件通過輸入裝置33的操作
可變設定。當從步驟S31到步驟S38的循環(huán)處理反復進行反復次數(shù)的閾值N次時,在步驟S35 的判定中判定為是。由此,循環(huán)處理完成。當循環(huán)處理完成時,在步驟S39中,第一相位校正量計算部43Α根據(jù)循環(huán)處理的結(jié)果算出第一相位校正量。第一相位校正量如式(10)所示,可以作為r空間數(shù)據(jù)與包括相位校正的循環(huán)處理后的2D-r空間數(shù)據(jù)Vhfi (χ, ζ)的相位差Φ (ζ)的指數(shù)函數(shù)β 1 (ζ)求出, 所述r空間數(shù)據(jù)是變換作為第一 HFI處理的原始數(shù)據(jù)的相位校正前的2D-k空間數(shù)據(jù)K (kx, kz)得到的。β l(z) = exp {-i Φ (ζ) } = conjugate [FFT {K (Kx, Kz) } ] / | FFT {K (Kx, Kz)} I * I Vhfi (x,z) I /Vhfi* (x, z) (10)其中,式(10)中conjugate ()是取復共軛的函數(shù),Vhfi*表示Vhfi的復共軛。接著,在步驟S4中,第二相位校正量決定部4 根據(jù)第一相位校正量β I(Z)決定第二相位校正量β2。在PE方向即y軸方向使第二相位校正量β2為恒定的情況下,如式 (11)所示,可以使第一相位校正量β (ζ)為3D相位校正用的第二相位校正量β2(ζ)。β 2 (ζ) = β l(z) (11)另一方面,在y軸方向的位置、即與PE量對應地使第二相位校正量β 2變化的情況下,如式(1 所示,將與y軸方向的位置y對應的系數(shù)Y (y)乘以第一相位校正量β l(z) 得到的值作為第二相位校正量i32(y,z)。在使第二相位校正量i32(y,z)在y軸方向呈線性地變化的情況下,系數(shù)Y (y)成為與位置y成比例的值。β2(γ, ζ) = β 1(ζ) γ (y) (12)并且,在設定與多個位置ky = kyi(i為自然數(shù))對應的多個2D_k空間數(shù)據(jù)K0(kx, kyi,kz)用于計算第一相位校正量β 1的情況下,求出與第一相位校正量β 1的計算用的 2D-k空間數(shù)據(jù)KO (kx, kyi, kz)的數(shù)量i相應的第一相位校正量β l(yi, ζ) 0在該情況下, 可以直接使用多個第一相位校正量0 1(yi,ζ),或通過基于多個第一相位校正量0 1(yi, ζ)的外插或內(nèi)插等插補,決定與位置y對應的第二相位校正量i32(y,z)。接著,在步驟S5中,診斷圖像生成部43C對收集到的3D的k空間中的MR信號執(zhí)行包括基于Homodyne濾波fh的數(shù)據(jù)填補用的濾波處理和使用第二相位校正量β 2的相位校正的第二 HFI處理。具體而言,如式(1 所示,執(zhí)行對在步驟S2中收集的非對稱的3D_k空間數(shù)據(jù) K (kx, ky,kz)進行數(shù)據(jù)填補用的Homodyne濾波處理fh和作為針對Homodyne濾波處理后的 3D-k空間數(shù)據(jù)的圖像重建處理的3D-FFT。由此,作為3D-r空間數(shù)據(jù)V(x,y,z)重建3D-MR 圖像數(shù)據(jù)。V (x, y, ζ) = FFT {K (kx, ky, kz) fh} (13)接著,如式(14)所示,進行使用第二相位校正量β 2的3D_r空間數(shù)據(jù)V(x,y,ζ) 的相位校正及相位校正后的3D-r空間數(shù)據(jù)的REAL化處理。由此,生成最終的MR診斷圖像數(shù)據(jù) Vfinal(x, 1,ζ)。Vfinal (x, y, ζ) = gREAL {V (χ, y, ζ) β 2} (14)另外,g是用于調(diào)整圖像信號的振幅的任意系數(shù),例如設定為g = 2. O。該系數(shù)g 相當于增益,能夠通過經(jīng)驗或模擬等決定。然后,圖像處理部44對由HFI處理部43取得的MR診斷圖像數(shù)據(jù)Vfinal實施失真校正處理等必要的圖像處理,并使其顯示于顯示裝置34。由此,使用者能夠觀察顯示于顯示裝置34上的MR診斷圖像來進行診斷。如上所述,根據(jù)磁共振成像裝置20,與以往相比,在HFI處理中能夠減少推定相位校正用的相位分布所需的FFT的處理量。特別是若使相位分布的推定用的k空間數(shù)據(jù)作為一張2D-k空間數(shù)據(jù),則能夠通過2D-FFT算出相位校正量。其結(jié)果,能夠減少總體的數(shù)據(jù)處理量和數(shù)據(jù)處理時間?;蛘?,能夠使總體的數(shù)據(jù)處理量和數(shù)據(jù)處理時間為相同程度、同時增加循環(huán)處理的反復次數(shù),從而提高精度。另外,也可以不是對r空間數(shù)據(jù)而是對k空間數(shù)據(jù)進行相位校正。只要對k空間數(shù)據(jù)進行相位校正,則處理就變得簡單。另一方面,只要對r空間數(shù)據(jù)進行相位校正,就能夠提高精度。以上,對特定的實施方式進行了記載,但記載的實施方式只不過是一個例子,并不限定發(fā)明的范圍。這里記載的新方法和裝置可以以各種各樣的其他方式具體化。并且,在這里記載的方法和裝置的方式中,在不脫離發(fā)明主旨的范圍內(nèi),能夠進行各種省略、置換和變更。所附權(quán)利要求及其等同物包括于發(fā)明范圍和主旨中,包括各種方式和變形例。
權(quán)利要求
1.一種磁共振成像裝置,該磁共振成像裝置具備數(shù)據(jù)收集單元,根據(jù)半傅里葉成像用的攝像條件,收集三維k空間中的磁共振信號;相位校正量計算單元,將比上述三維k空間中的磁共振信號少的磁共振信號作為用于求出第一相位校正量的k空間數(shù)據(jù),對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括基于用于求出上述第一相位校正量的k空間數(shù)據(jù)的相位校正、以及利用了 k空間的復對稱性對非采樣區(qū)域進行的數(shù)據(jù)填補,由此,求出上述第一相位校正量;以及圖像數(shù)據(jù)生成單元,對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括使用了基于上述第一相位校正量的第二相位校正量的相位校正、上述數(shù)據(jù)填補以及圖像重建處理,由此,生成磁共振圖像數(shù)據(jù)。
2.根據(jù)權(quán)利要求1所述的磁共振成像裝置,其中,上述相位校正量計算單元構(gòu)成為,將包括編碼量不是零的二維k空間數(shù)據(jù)的多個二維 k空間數(shù)據(jù)作為用于求出上述第一相位校正量的k空間數(shù)據(jù),求出多個第一相位校正量。
3.根據(jù)權(quán)利要求1或2所述的磁共振成像裝置,其中,上述數(shù)據(jù)收集單元構(gòu)成為收集在相位編碼方向非對稱的磁共振信號,上述相位校正量計算單元構(gòu)成為,將分層編碼方向的編碼量為0的二維k空間數(shù)據(jù)作為用于求出上述第一相位校正量的k空間數(shù)據(jù)。
4.根據(jù)權(quán)利要求1或2所述的磁共振成像裝置,其中,上述數(shù)據(jù)收集單元構(gòu)成為收集在分層編碼方向非對稱的磁共振信號,上述相位校正量計算單元構(gòu)成為,將相位編碼方向的編碼量為0的二維k空間數(shù)據(jù)作為用于求出上述第一相位校正量的k空間數(shù)據(jù)。
5.根據(jù)權(quán)利要求1 4中的任一項所述的磁共振成像裝置,其中,上述圖像數(shù)據(jù)生成單元構(gòu)成為,將上述第一相位校正量作為上述第二相位校正量。
6.根據(jù)權(quán)利要求1 4中的任一項所述的磁共振成像裝置,其中,上述圖像數(shù)據(jù)生成單元構(gòu)成為,根據(jù)上述第一相位校正量求出與編碼量對應的多個第二相位校正量。
全文摘要
實施方式的磁共振成像裝置,具備數(shù)據(jù)收集單元、相位校正量計算單元以及圖像數(shù)據(jù)生成單元。數(shù)據(jù)收集單元根據(jù)半傅里葉成像用的攝像條件收集三維k空間中的磁共振信號。相位校正量計算單元將比上述三維k空間中的磁共振信號少的磁共振信號作為用于求出第一相位校正量的k空間數(shù)據(jù),對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括基于用于求出上述第一相位校正量的k空間數(shù)據(jù)的相位校正、以及利用k空間的復對稱性對非采樣區(qū)域進行的數(shù)據(jù)填補,由此,求出上述第一相位校正量。圖像數(shù)據(jù)生成單元對上述三維k空間中的磁共振信號執(zhí)行如下處理該處理包括使用了基于上述第一相位校正量的第二相位校正量的相位校正、上述數(shù)據(jù)填補以及圖像重建處理,由此生成磁共振圖像數(shù)據(jù)。
文檔編號A61B5/055GK102440778SQ20111030510
公開日2012年5月9日 申請日期2011年9月30日 優(yōu)先權(quán)日2010年9月30日
發(fā)明者高井博司 申請人:東芝醫(yī)療系統(tǒng)株式會社, 株式會社東芝