磁共振成像裝置以及磁共振成像方法【專利摘要】實(shí)施方式的磁共振成像裝置具備數(shù)據(jù)收集單元和圖像生成單元。數(shù)據(jù)收集單元從被檢體收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的磁共振數(shù)據(jù)。圖像生成單元,根據(jù)向上述磁共振數(shù)據(jù)的非采樣區(qū)域進(jìn)行的零填充后的第一k空間數(shù)據(jù),生成實(shí)空間的振幅圖像數(shù)據(jù),執(zhí)行將上述振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并在進(jìn)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理、或者利用將上述濾波用的窗函數(shù)進(jìn)行變換而得到的實(shí)空間的函數(shù)進(jìn)行的上述振幅圖像數(shù)據(jù)的卷積處理,由此生成磁共振圖像數(shù)據(jù)?!緦@f明】磁共振成像裝置以及磁共振成像方法【
技術(shù)領(lǐng)域:
】[0001]本發(fā)明的實(shí)施方式涉及磁共振成像(MR1:MagneticResonanceImaging)裝置以及磁共振成像方法。【
背景技術(shù):
】[0002]MRI是一種攝像方法,將置于靜磁場中的被檢體的原子核自旋用拉莫爾頻率的高頻(RF:radiofrequency,射頻)信號(hào)進(jìn)行磁激勵(lì),根據(jù)伴隨該激勵(lì)而產(chǎn)生的核磁共振(NMR:nuclearmagneticresonance)信號(hào)對(duì)圖像進(jìn)行重構(gòu)。[0003]作為MRI的成像法之一,公知有AFI(AsymmetricFourierImaging,非對(duì)稱傅里葉成像)法。AFI法是這樣一種方法,以在k空間(k-space)中在波數(shù)方向上成為非對(duì)稱的方式將數(shù)據(jù)采樣,利用根據(jù)采樣后的自身數(shù)據(jù)而推定的相位分布進(jìn)行相位修正后,對(duì)圖像數(shù)據(jù)進(jìn)行重構(gòu)。根據(jù)AFI法,能夠生成與根據(jù)在k空間中對(duì)稱采樣的數(shù)據(jù)而生成的圖像數(shù)據(jù)同等的圖像數(shù)據(jù)。[0004]AFI法中,提出了Margosian法、FIR(finiteimpulseresponse)法、MoFIR(ModifiedFIR)法、POCS(projectionontoconvexsets)法或hybrid法等各種各樣的方法。此外,公知有非AFI法的、以在k空間中成為非對(duì)稱的方式將數(shù)據(jù)采樣、在無數(shù)據(jù)的部分進(jìn)行O填充(Ο-filling)后進(jìn)行傅里葉變換(FT:Fouriertransform)的非對(duì)稱數(shù)據(jù)采樣法。已知在Ο-filling中,k空間數(shù)據(jù)的非對(duì)稱度越大,在圖像中產(chǎn)生越顯著的模糊。[0005]在Margosian法中,在對(duì)非對(duì)稱地米樣的k空間數(shù)據(jù)乘以作為窗(window)函數(shù)的homodynefilter后,通過FT生成與非對(duì)稱k空間數(shù)據(jù)對(duì)應(yīng)的實(shí)空間(r_space)數(shù)據(jù)。另一方面,在非對(duì)稱地采樣的k空間數(shù)據(jù)中,根據(jù)對(duì)稱地采樣的k空間的中心附近的低頻區(qū)域中的k空間數(shù)據(jù)來推定相位分布。并且,利用所推定的相位分布,進(jìn)行與非對(duì)稱k空間數(shù)據(jù)對(duì)應(yīng)的實(shí)空間數(shù)據(jù)的相位修正。[0006]POCS法是Margosian法的改良法,是在Margosian法之后進(jìn)行POCS循環(huán)處理的手法。POCS循環(huán)處理是通過反復(fù)進(jìn)行REAL化處理、合成處理以及相位修正處理而使虛部的變化收斂于閾值以下的處理,REAL化處理是僅將相位修正后的實(shí)空間數(shù)據(jù)的實(shí)部留下而將虛部設(shè)為零的處理,合成處理是將使REAL化了的實(shí)空間數(shù)據(jù)的相位復(fù)原后通過逆傅里葉變換(IFT:inverseFouriertransform)得到的k空間數(shù)據(jù)的非采樣部和原數(shù)據(jù)的采樣部合成的處理,相位修正處理是針對(duì)通過合成處理后的k空間數(shù)據(jù)的FT而得到的實(shí)空間數(shù)據(jù)進(jìn)行的處理。POCS法基于這樣的原理,即:若相位修正是完全的,則實(shí)空間數(shù)據(jù)的虛部成為零。根據(jù)POCS法,通過反復(fù)進(jìn)行幾次POCS循環(huán)處理,能夠降低在Margosian法中因homodynefilter處理而產(chǎn)生的相位修正的錯(cuò)誤。[0007]另一方面,在FIR法中,在對(duì)非對(duì)稱地采樣而得的k空間數(shù)據(jù)乘以homodynefilter之前進(jìn)行相位修正。即,在FIR法中,在通過非對(duì)稱k空間數(shù)據(jù)的FT而生成的實(shí)空間數(shù)據(jù)的相位修正后,通過IFT將相位修正后的實(shí)空間數(shù)據(jù)變換為k空間數(shù)據(jù)。并且,對(duì)相位修正后的k空間數(shù)據(jù)乘以homodynefilter。該FIR法中,與Margosian法相比,數(shù)據(jù)處理時(shí)間延長2次的FT的量,但由于在homodynefilter處理前進(jìn)行相位修正,所以能夠降低由homodynefilter處理引起的相位修正的錯(cuò)誤。[0008]MoFIR法是FIR法的改良法,根據(jù)不僅包含k空間數(shù)據(jù)被對(duì)稱采樣了的低頻區(qū)域、還包含非對(duì)稱采樣了的部分在內(nèi)的全部的k空間數(shù)據(jù),推定在相位修正用中使用的相位分布。即,F(xiàn)IR法中,根據(jù)k空間的中心附近的僅低頻區(qū)域的k空間數(shù)據(jù),推定相位修正用的低頻區(qū)域的相位分布,相對(duì)于此,MoFIR法中,根據(jù)非對(duì)稱地采樣而得的全部k空間數(shù)據(jù),推定相位修正用的相位分布。因此,MoFIR法中,與FIR法相比較,雖然與原來的相位分布不同,但是能夠推定更高頻區(qū)域的相位分布。因此,根據(jù)MoFIR法,能夠降低Margosian法或FIR法中的由homodynefilter處理引起的相位修正的錯(cuò)誤。[0009]另一方面,在作為最單純的針對(duì)非對(duì)稱采樣數(shù)據(jù)的重構(gòu)法的Ο-filling中,在圖像中產(chǎn)生模糊。但是,在對(duì)稱地采樣的數(shù)據(jù)的70%以上被采樣那樣的非對(duì)稱度比較小的非對(duì)稱采樣的情況下,圖像的模糊成為許容范圍。此外,在進(jìn)行Ο-filling的情況下,不僅不需要特別的處理,而且也不會(huì)出現(xiàn)AFI中的因相位的過修正而引起的偽像(artifact)。因此,0-fi11ing在非對(duì)稱度比較小的情況下依然被大量應(yīng)用。[0010]hybrid法是將o-filling及AFI法組合而得到的手法。即,hybrid法是將通過Ο-filling得到的Ο-fill圖像和通過AFI得到的AFI圖像進(jìn)行加權(quán)加算合成的手法。更具體而言,對(duì)權(quán)重進(jìn)行調(diào)整,使得在Ο-fill圖像與AFI圖像之間相位差或振幅差大的部分成為Ο-fill圖像,小的部分成為AFI圖像。[0011]現(xiàn)有技術(shù)文獻(xiàn)[0012]專利文獻(xiàn)[0013]專利文獻(xiàn)1:日本特開2010-217981號(hào)公報(bào)【
發(fā)明內(nèi)容】[0014]發(fā)明要解決的技術(shù)問題[0015]根據(jù)上述那樣的AFI法,若將k空間數(shù)據(jù)的成為非對(duì)稱的方向作為k空間數(shù)據(jù)的讀出方向,則能夠縮短回波時(shí)間(TE:echotime)。另一方面,若將k空間數(shù)據(jù)的成為非對(duì)稱的方向作為編碼方向,則導(dǎo)致攝像時(shí)間的縮短。但是,在AFI法中,希望兼顧進(jìn)一步的攝像時(shí)間的縮短和畫質(zhì)的提高。[0016]因此,本發(fā)明的目的在于,提供一種磁共振成像裝置以及磁共振成像方法,能夠根據(jù)在k空間中非對(duì)稱地采樣而得的MR數(shù)據(jù),更高速地生成更高精度的圖像數(shù)據(jù)。[0017]用于解決技術(shù)問題的手段[0018]本發(fā)明的實(shí)施方式的磁共振成像裝置具備數(shù)據(jù)收集單元和圖像生成單元。數(shù)據(jù)收集單元從被檢體收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的磁共振數(shù)據(jù)。圖像生成單元,根據(jù)向上述磁共振數(shù)據(jù)的非采樣區(qū)域進(jìn)行的零填充后的第一k空間數(shù)據(jù),生成實(shí)空間的振幅圖像數(shù)據(jù),執(zhí)行將上述振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并在進(jìn)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理、或者利用將上述濾波用的窗函數(shù)進(jìn)行變換而得到的實(shí)空間的函數(shù)進(jìn)行的上述振幅圖像數(shù)據(jù)的卷積處理,由此生成磁共振圖像數(shù)據(jù)。[0019]此外,本發(fā)明的實(shí)施方式的磁共振成像方法具有:從被檢體收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的磁共振數(shù)據(jù)的步驟;以及根據(jù)向上述磁共振數(shù)據(jù)的非采樣區(qū)域進(jìn)行的零填充后的第一k空間數(shù)據(jù)生成實(shí)空間的振幅圖像數(shù)據(jù),執(zhí)行將上述振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并在執(zhí)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理、或者利用將上述濾波用的窗函數(shù)進(jìn)行變換而得到的實(shí)空間的函數(shù)進(jìn)行的上述振幅圖像數(shù)據(jù)的卷積處理,由此生成磁共振圖像數(shù)據(jù)的步驟?!緦@綀D】【附圖說明】[0020]圖1是本發(fā)明的實(shí)施方式的磁共振成像裝置的構(gòu)成圖。[0021]圖2是圖1所示的計(jì)算機(jī)的功能框圖。[0022]圖3是說明圖2所示的AFI處理部中執(zhí)行的包含k空間中的濾波的AFI處理的方法的圖。[0023]圖4是表示式(2)所示的窗函數(shù)Hlw(k)的常數(shù)K1及K2的定義的圖。[0024]圖5是表示式(2)所示的低頻區(qū)域的提取用的窗函數(shù)H1ot(k)的一例的圖。[0025]圖6是表示式(I)所示的Ο-filling用的窗函數(shù)Hwtole(k)的一例的圖。[0026]圖7是表示對(duì)稱的窗函數(shù)Hh_的例子的圖。[0027]圖8是表示非對(duì)稱的窗函數(shù)Hh_的例子的圖。[0028]圖9是說明圖2所示的AFI處理部中執(zhí)行的包含實(shí)空間中的卷積處理的AFI處理的方法的圖。[0029]圖10是說明接續(xù)圖3所示的數(shù)據(jù)處理而執(zhí)行使相位誤差降低的循環(huán)處理的情況下的AFI處理的方法的圖。[0030]圖11是將通過圖3所示的AFI處理而生成的一維(ID:onedimensional)仿真數(shù)據(jù)與通過以往的AFI處理而生成的ID仿真數(shù)據(jù)進(jìn)行比較的圖?!揪唧w實(shí)施方式】[0031]參照附圖對(duì)本發(fā)明的實(shí)施方式的磁共振成像裝置以及磁共振成像方法進(jìn)行說明。[0032]圖1是本發(fā)明的實(shí)施方式的磁共振成像裝置的構(gòu)成圖。[0033]磁共振成像裝置20具備:形成靜磁場的筒狀的靜磁場用磁體21、在該靜磁場用磁體21的內(nèi)部設(shè)置的勻場線圈22、梯度磁場線圈23以及RF線圈24。[0034]此外,磁共振成像裝置20具備控制系統(tǒng)25??刂葡到y(tǒng)25具備靜磁場電源26、梯度磁場電源27、勻場線圈電源28、發(fā)送器29、接收器30、序列控制器31以及計(jì)算機(jī)32。控制系統(tǒng)25的梯度磁場電源27包括X軸梯度磁場電源27x、Y軸梯度磁場電源27y以及Z軸梯度磁場電源27z。此外,[0035]計(jì)算機(jī)32具備輸入裝置33、顯示裝置34、運(yùn)算裝置35以及存儲(chǔ)裝置36。[0036]靜磁場用磁體21與靜磁場電源26連接,具有利用從靜磁場電源26供給的電流在攝像區(qū)域形成靜磁場的功能。另外,靜磁場用磁體21多由超導(dǎo)線圈構(gòu)成,在勵(lì)磁時(shí)與靜磁場電源26連接而被供給電流,而一旦被激勵(lì)后通常被設(shè)為非連接狀態(tài)。此外,也存在由永久磁體構(gòu)成靜磁場用磁體21、而不設(shè)置靜磁場電源26的情況。[0037]此外,在靜磁場用磁體21的內(nèi)側(cè),在同軸上設(shè)有筒狀的勻場線圈22。勻場線圈22與勻場線圈電源28連接,從勻場線圈電源28向勻場線圈22供給電流而將靜磁場均勻化。[0038]梯度磁場線圈23包括X軸梯度磁場線圈23x、Y軸梯度磁場線圈23y以及Z軸梯度磁場線圈23z,在靜磁場用磁體21的內(nèi)部形成為筒狀。在梯度磁場線圈23的內(nèi)側(cè)設(shè)置診視床37作為攝像區(qū)域,在診視床37上設(shè)置被檢體P。RF線圈24有在架臺(tái)中內(nèi)置的RF信號(hào)的收發(fā)用的全身用線圈(WBC:wholebodycoil)、設(shè)置在診視床37或被檢體P附近的RF信號(hào)的接收用的局部線圈等。[0039]此外,梯度磁場線圈23與梯度磁場電源27連接。梯度磁場線圈23的X軸梯度磁場線圈23x、Y軸梯度磁場線圈23y以及Z軸梯度磁場線圈23z分別與梯度磁場電源27的X軸梯度磁場電源27x、Y軸梯度磁場電源27y以及Z軸梯度磁場電源27z連接。[0040]并且,構(gòu)成為,利用從X軸梯度磁場電源27x、Y軸梯度磁場電源27y以及Z軸梯度磁場電源27z分別向X軸梯度磁場線圈23x、Y軸梯度磁場線圈23y以及Z軸梯度磁場線圈23z供給的電流,能夠在攝像區(qū)域分別形成X軸方向的梯度磁場Gx、Y軸方向的梯度磁場Gy、Z軸方向的梯度磁場Gz。[0041]RF線圈24與發(fā)送器29及接收器30的至少一方連接。發(fā)送用的RF線圈24具有從發(fā)送器29接受RF信號(hào)并向被檢體P發(fā)送的功能,接收用的RF線圈24具有接收伴隨著被檢體P內(nèi)部的原子核自旋的RF信號(hào)的激勵(lì)而發(fā)生的NMR信號(hào)并提供給接收器30的功能。[0042]另一方面,控制系統(tǒng)25的序列控制器31與梯度磁場電源27、發(fā)送器29以及接收器30連接。序列控制器31具有存儲(chǔ)序列信息的功能、和按照所存儲(chǔ)的規(guī)定的序列來使梯度磁場電源27、發(fā)送器29以及接收器30驅(qū)動(dòng)從而產(chǎn)生X軸梯度磁場Gx、Y軸梯度磁場Gy、Z軸梯度磁場Gz以及RF信號(hào)的功能,序列信息記載了使梯度磁場電源27、發(fā)送器29以及接收器30驅(qū)動(dòng)所需要的控制信息、例如應(yīng)向梯度磁場電源27施加的脈沖電流的強(qiáng)度或施加時(shí)間、施加定時(shí)等動(dòng)作控制信息。[0043]此外,序列控制器31接受原始數(shù)據(jù)(rawdata)并提供給計(jì)算機(jī)32,原始數(shù)據(jù)是通過接收器30中的NMR信號(hào)的檢波以及A/D(analogtodigital:模數(shù))變換而得到的復(fù)數(shù)(complex)數(shù)據(jù)。[0044]因此,發(fā)送器29具備根據(jù)從序列控制器31接受的控制信息、將RF信號(hào)提供給RF線圈24的功能,另一方面,接收器30具備對(duì)從RF線圈24接受的NMR信號(hào)進(jìn)行檢波且執(zhí)行所需的信號(hào)處理并且進(jìn)行A/D變換、從而生成作為被數(shù)字化后的復(fù)數(shù)數(shù)據(jù)的原始數(shù)據(jù)的功能、以及將生成的原始數(shù)據(jù)提供給序列控制器31的功能。[0045]進(jìn)而,磁共振成像裝置20具備取得被檢體P的ECG(electrocardiogram,心電圖)信號(hào)的ECG單元38。構(gòu)成為,由ECG單元38取得的ECG信號(hào)經(jīng)由序列控制器31向計(jì)算機(jī)32輸出。另外,還能取得將搏動(dòng)作為脈波信息來表示的脈波同步(PPG:peripheralpulsegating)信號(hào),來替代將搏動(dòng)作為心跳信息來表示的ECG信號(hào)。[0046]此外,通過使運(yùn)算裝置35執(zhí)行在計(jì)算機(jī)32的存儲(chǔ)裝置36中保存的程序,使計(jì)算機(jī)32具備各種功能。但是,也可以替代程序的至少一部分而對(duì)磁共振成像裝置20設(shè)置具有各種功能的特定的電路。[0047]圖2是圖1所示的計(jì)算機(jī)32的功能框圖。[0048]計(jì)算機(jī)32的運(yùn)算裝置35通過執(zhí)行在存儲(chǔ)裝置36中保存的程序而作為攝像條件設(shè)定部40、數(shù)據(jù)處理部41以及AFI處理?xiàng)l件設(shè)定部42發(fā)揮功能。此外,存儲(chǔ)裝置36作為k空間數(shù)據(jù)庫43以及圖像數(shù)據(jù)庫44發(fā)揮功能。數(shù)據(jù)處理部41具有AFI處理部41A。[0049]攝像條件設(shè)定部40具有如下功能:根據(jù)來自輸入裝置33的指示信息,設(shè)定包含脈沖序列的攝像條件,將設(shè)定的攝像條件向序列控制器31輸出,由此驅(qū)動(dòng)控制序列控制器31。[0050]特別是,攝像條件設(shè)定部40能夠設(shè)定從被檢體P收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的MR數(shù)據(jù)時(shí)的AFI的攝像條件。AFI法是一種圖像重構(gòu)方法,根據(jù)在二維(2D:twodimensional)或三維(3D:threedimensional)的k空間的至少I軸方向的波數(shù)方向上非對(duì)稱地采樣而得的MR數(shù)據(jù),近似地生成與根據(jù)對(duì)稱地采樣而得的MR數(shù)據(jù)的圖像數(shù)據(jù)相近的圖像數(shù)據(jù)。AFI能夠應(yīng)用于包含磁共振血管攝影法(MRA:magnetiCresonanceangiography)、擴(kuò)散強(qiáng)調(diào)成像(DW1:diffusionweightedimaging)的各種各樣的成像。作為用于在k空間中進(jìn)行AFI的攝像條件,如果是2D的采樣,則存在k空間數(shù)據(jù)的讀出(readout)方向及相位編碼(phaseencode)方向中的一方向k的波數(shù)方向上的采樣區(qū)域(-Kc≤k≤Kmax).采樣區(qū)域的邊界K。例如能夠按照來自AFI處理?xiàng)l件設(shè)定部42的設(shè)定信息進(jìn)行可變設(shè)定。非采樣區(qū)域可以是k方向的正側(cè)以及負(fù)側(cè)的任一個(gè)。此外,作為3D的區(qū)域也能夠設(shè)定采樣區(qū)域和非采樣區(qū)域。這里以作為2D區(qū)域設(shè)定采樣區(qū)域和非采樣區(qū)域、并且將負(fù)側(cè)作為非采樣區(qū)域的情況為例進(jìn)行說明。[0051]數(shù)據(jù)處理部41具有從序列控制器31接收通過AFI用的條件而收集的原始數(shù)據(jù)并作為k空間數(shù)據(jù)配置在形成于k空間數(shù)據(jù)庫43的k空間中的功能、以及通過包含針對(duì)k空間數(shù)據(jù)的FT的圖像重構(gòu)處理來生成圖像數(shù)據(jù)的功能。診斷用的圖像數(shù)據(jù)能夠顯示在顯示裝置34上。此外,還能夠根據(jù)需要將診斷用的圖像數(shù)據(jù)、中間性的圖像數(shù)據(jù)等實(shí)空間數(shù)據(jù)保存在圖像數(shù)據(jù)庫44中。[0052]因此,在k空間數(shù)據(jù)庫43中保存k空間數(shù)據(jù),另一方面,在圖像數(shù)據(jù)庫44中保存實(shí)空間數(shù)據(jù)。[0053]數(shù)據(jù)處理部41的AFI處理部41A具有執(zhí)行AFI處理的功能,AFI處理是這樣一種處理,即:根據(jù)在k空間中在波數(shù)方向上非對(duì)稱地采樣而得的k空間數(shù)據(jù),不伴隨相位修正處理地生成與從對(duì)稱地采樣而得的k空間數(shù)據(jù)生成的圖像數(shù)據(jù)大致同等的圖像數(shù)據(jù)。即,不同于以往的伴隨相位修正處理的AFI處理,AFI處理部41A執(zhí)行不伴隨相位修正處理的AFI處理。[0054]AFI處理能夠按如下順序來執(zhí)行,即:通過向從非對(duì)稱采樣區(qū)域收集的MR數(shù)據(jù)的非采樣區(qū)域進(jìn)行的o-filling生成第一k空間數(shù)據(jù),通過第一k空間數(shù)據(jù)的FT生成振幅圖像數(shù)據(jù),通過振幅圖像數(shù)據(jù)的IFT變換為第二k空間數(shù)據(jù),通過第二k空間數(shù)據(jù)的濾波以及濾波后的第二k空間數(shù)據(jù)的FT生成圖像數(shù)據(jù)。因此,AFI處理部41A根據(jù)向MR數(shù)據(jù)的非采樣區(qū)域進(jìn)行o-filling后的第一k空間數(shù)據(jù)生成實(shí)空間的振幅圖像數(shù)據(jù),將振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并通過進(jìn)行濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理生成MR圖像數(shù)據(jù)。[0055]圖3是說明圖2所示的AFI處理部41A中執(zhí)行的包含k空間中的濾波的AFI處理的方法的圖。[0056]若按照將2D的k空間的負(fù)側(cè)作為非采樣區(qū)域的AFI法的攝像條件收集MR數(shù)據(jù),則收集在k空間中在波數(shù)方向上非對(duì)稱的MR數(shù)據(jù)SOTig(k)。如圖3所示那樣在成為非對(duì)稱的方向是k空間的kx方向的情況下,MR數(shù)據(jù)SOTig(k)成為將滿足-K。蘭k蘭Kmax的頻率k作為參數(shù)的數(shù)據(jù)群。其中,K。是取正值的非采樣區(qū)域的截止頻率,Kmax是采樣區(qū)域的最大頻率。因此,(KKcXKmax的關(guān)系成立。[0057]接著,執(zhí)行向非對(duì)稱的MR數(shù)據(jù)SOTig(k)的非采樣區(qū)域的0-filling。0-filling能夠通過對(duì)非對(duì)稱的MR數(shù)據(jù)SOTig(k)乘以式(I)定義的窗函數(shù)Hwh-(k)來執(zhí)行。[0058]Hwhole(k)=Hlow(k):-Kmax=k=O[0059]=l:0<k^Kmax[0060](I)[0061]式(I)中的函數(shù)Hlw(k)是從非對(duì)稱的MR數(shù)據(jù)SOTig(k)提取低頻區(qū)域中的對(duì)稱部分的窗函數(shù)。在將提取低頻區(qū)域的窗函數(shù)用Gaussian定義的情況下,能夠如式(2)那樣表不。[0062]Hlow(k)=1:1k=Kc-K1[0063]=exp[(-1n2){(k_(Kc-K1))/K2}2]=Kc-K1C|k|=Kmax[0064](2)[0065]式(2)中的K1及K2是用于使成為提取對(duì)象的低頻區(qū)域平滑地變化的常數(shù)。[0066]圖4是表示式(2)所示的窗函數(shù)Hlw(k)中的常數(shù)K1及K2的定義的圖。[0067]圖4中,橫軸表示k空間的頻率k,縱軸表示窗函數(shù)Hlw(k)的增益。如圖4所示,能夠?qū)⑹勾昂瘮?shù)H1ot(k)的增益從最大值I平滑地變化到O的頻率的寬度定義為常數(shù)K1,另一方面,能夠?qū)⑹勾昂瘮?shù)H1ot(k)的增益從最大值I平滑地變化到1/2的頻率的寬度定義為常數(shù)K2。[0068]圖5是表示式(2)所示的低頻區(qū)域的提取用的窗函數(shù)H1ot(k)的一例的圖,圖6是表示式(I)所示的Ο-filling用的窗函數(shù)Hwtole(k)的一例的圖。[0069]圖5及圖6的各橫軸表示k空間的頻率k,各縱軸表示窗函數(shù)Hlw(k),Hwhole(k)的增益。利用圖5及圖6中例示那樣的低頻區(qū)域的提取用的窗函數(shù)HlOT(k)以及0-filling用的窗函數(shù)Hwtole(k),能夠執(zhí)行非對(duì)稱的MR數(shù)據(jù)SOTig(k)的0-filling。[0070]若非對(duì)稱的MR數(shù)據(jù)SOTig(k)的0-filling完成,則生成第一k空間數(shù)據(jù)SOTig(k)XHwhole(k)。接著對(duì)0-filling后的第一k空間數(shù)據(jù)SOTig(k)XHwhole(k)執(zhí)行FT。由此,作為實(shí)空間數(shù)據(jù)而生成振幅圖像數(shù)據(jù)SQ_filled(r)。接著,通過對(duì)0-filling后的振幅圖像數(shù)據(jù)SQ_filled(r)執(zhí)行IFT,生成第二k空間數(shù)據(jù)S’(k)。[0071]0-filling后的振幅圖像數(shù)據(jù)SQ_fille(1(r)的相位成分成為O。即,振幅圖像數(shù)據(jù)S0_filled(r)是僅實(shí)部的數(shù)據(jù),虛部成為O。另一方面,在作為以往的AFI法的一個(gè)方法的MoFIR法中,利用全部的相位的數(shù)據(jù)執(zhí)行實(shí)空間數(shù)據(jù)的相位修正。因此,相位修正后的實(shí)空間數(shù)據(jù)的相位成分也成為O。[0072]因此,可以認(rèn)為,與Ο-filling后的振幅圖像數(shù)據(jù)SQ_filledCr)對(duì)應(yīng)的第二k空間數(shù)據(jù)S’(k)的分布與MoFIR法的相位修正后的k空間數(shù)據(jù)的分布同等。因此,能夠?qū)⑴cΟ-filling后的振幅圖像數(shù)據(jù)SQ_filled(r)對(duì)應(yīng)的第二k空間數(shù)據(jù)S’(k)作為MoFIR法的相位修正后的k空間數(shù)據(jù)來處理。[0073]但是,MoFIR法中,在相位修正后,為了對(duì)高頻成分的劣化進(jìn)行修正,乘以Homodyne-high-passfilter。Homodyne-high-passfilter是增益非對(duì)稱的AFI處理用的窗函數(shù)。具體而言,Homodynehigh-passfilter是非采樣區(qū)域的增益為O、采樣區(qū)域中的對(duì)稱部分的增益為I倍、符號(hào)與非采樣區(qū)域相反的采樣部分的增益為2倍的濾波器。通常,Homodyne-high-passfilter的增益以平滑變化的方式被決定。[0074]因此,通過對(duì)與Ο-filling后的振幅圖像數(shù)據(jù)SQ_filled(r)對(duì)應(yīng)的第二k空間數(shù)據(jù)S’(k)實(shí)施與Homodynehigh-passfilter等價(jià)的濾波處理,能夠?qū)崿F(xiàn)高頻成分的劣化即模糊的降低。由此,能夠生成具有與通過MoFIR法生成的MR圖像數(shù)據(jù)同等的畫質(zhì)的MR圖像數(shù)據(jù)。[0075]第二k空間數(shù)據(jù)S’(k)的濾波處理能夠通過對(duì)第二k空間數(shù)據(jù)S’(k)乘以窗函數(shù)Hh_來執(zhí)行。作為窗函數(shù)Hh_,可以考慮非對(duì)象的函數(shù)和對(duì)稱的函數(shù)。[0076]通過振幅圖像數(shù)據(jù)SQ_filled(r)的IFT而得到的第二k空間數(shù)據(jù)S’(k)的特性為復(fù)共軛。即,第二k空間數(shù)據(jù)S’(k)的實(shí)部是偶函數(shù),虛部是奇函數(shù)。因此,作為窗函數(shù)Hh_,采用非對(duì)象的函數(shù)或是采用對(duì)稱的函數(shù)也能得到相同的結(jié)果。[0077]第二k空間數(shù)據(jù)S’(k)的非采樣區(qū)域中也存在具有進(jìn)行了全采樣的情況下的1/2的強(qiáng)度的數(shù)據(jù)。另一方面,符號(hào)與非采樣區(qū)域相反的采樣部分中的數(shù)據(jù)的強(qiáng)度也成為進(jìn)行了全采樣的情況下的1/2的強(qiáng)度。[0078]因此,在第二k空間數(shù)據(jù)S’(k)的濾波處理中使用的模糊的修正用的窗函數(shù)Hh_適合于作為將第二k空間數(shù)據(jù)S’(k)的信號(hào)強(qiáng)度進(jìn)行修正以使其成為進(jìn)行了全采樣的情況下的強(qiáng)度的函數(shù)。由此,能夠抑制與由相位修正引起的錯(cuò)誤相當(dāng)?shù)腻e(cuò)誤導(dǎo)致的偽像,并且能夠降低由Ο-filling引起的模糊。[0079]圖7是表示對(duì)稱的窗函數(shù)Hh_的例子的圖。[0080]圖7中,橫軸表示k空間的頻率k,縱軸表示窗函數(shù)Hh_(k)的增益。如圖7所示,能夠利用窗函數(shù)Hh_(k)執(zhí)行對(duì)第二k空間數(shù)據(jù)S’(k)的濾波,窗函數(shù)Hh_(k)被設(shè)定為,與非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分和非采樣區(qū)域構(gòu)成的高頻側(cè)的區(qū)域的增益更大。即,通過對(duì)第二k空間數(shù)據(jù)S’(k)乘以相對(duì)于k空間的中心對(duì)稱的窗函數(shù)Hh_(k),能夠修正第二k空間數(shù)據(jù)S’(k)的高頻成分的劣化。[0081]利用相對(duì)于k空間的中心而線對(duì)稱的窗函數(shù)Hh_(k)的情況下,由于上述的理由,將高頻側(cè)的區(qū)域中的增益設(shè)定為低頻側(cè)的區(qū)域中的增益的2倍以下是適當(dāng)?shù)?。圖7所示的例子中,與高頻側(cè)的區(qū)域?qū)?yīng)的窗函數(shù)Hh_(k)的增益的最大值Hniax被設(shè)定為理論上的理想值2,與低頻側(cè)的區(qū)域?qū)?yīng)的增益被設(shè)定為I。[0082]將線對(duì)稱的窗函數(shù)Hh_(k)作為Hhigh.sym(k)時(shí),線對(duì)稱的窗函數(shù)Hhigh.sym(k)能夠通過式(3)求出。[0083]Hhigh.sym(k)=Hlow(k):|k|=Kc[0084]=Hmax-(Hmax-1)XHlow(k):Kc<|k|^Kfflax[0085](3)[0086]其中,1〈H—^2。此外,式(3)中的函數(shù)H1ot(k)是式(2)中定義的低頻區(qū)域的提取用的窗函數(shù)Hlw(k)。因此,線對(duì)稱的窗函數(shù)Hhigh.sym(k)利用Gauss函數(shù)而成為k=±K。前后的增益的過渡部分平滑變化的函數(shù)。通過將這樣平滑地變化的函數(shù)作為對(duì)稱或非對(duì)稱的窗函數(shù)Hh_(k),能夠?qū)⒌诙空間數(shù)據(jù)S’(k)修正為連續(xù)的k空間數(shù)據(jù)。[0087]圖8是表示非對(duì)稱的窗函數(shù)Hh_的例子的圖。[0088]圖8中,橫軸表示k空間的頻率k,縱軸表示窗函數(shù)Hh_(k)的增益。如圖8所示,能夠利用窗函數(shù)Hh_(k)執(zhí)行對(duì)第二k空間數(shù)據(jù)S’(k)的濾波,窗函數(shù)Hh_(k)被設(shè)定為,與非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分構(gòu)成的高頻側(cè)的區(qū)域的增益更大。即,通過對(duì)第二k空間數(shù)據(jù)S’(k)乘以相對(duì)于k空間的中心非對(duì)稱的窗函數(shù)Hh_(k),能夠修正第二k空間數(shù)據(jù)S’(k)的高頻成分的劣化。[0089]在使用關(guān)于k空間的中心非對(duì)稱的窗函數(shù)Hh_(k)的情況下,由于上述的理由,將高頻側(cè)的采樣部分的增益設(shè)定為低頻側(cè)的采樣部分的增益的2倍以上4倍以下、并且將非采樣區(qū)域的增益設(shè)定為O是適當(dāng)?shù)?。圖8所示的例子中,高頻側(cè)的采樣部分的增益的最大值Hmax設(shè)定為理論上的理想值4,低頻側(cè)的采樣部分的增益設(shè)定為1,非采樣區(qū)域的增益設(shè)定為O。[0090]若將非對(duì)稱的窗函數(shù)Hh_(k)作為Hhighasyni(k),則非對(duì)稱的窗函數(shù)Hhighasyni(k)能夠通過式(4)求出。[0091]Hhigh.asym(k)=Hlow(k):-Kmax=k=O[0092]=Hmax-(Hmax-1)XHlow(k):0<k^Kfflax[0093](4)[0094]其中,2^Hfflax^4。此外,在式(4)中也利用式(2)中定義的低頻區(qū)域的提取用的窗函數(shù)Hlw(k)。因此,非對(duì)稱的窗函數(shù)Hhigh.asym(k)也通過Gauss函數(shù)而成為平滑變化的函數(shù)。另外,圖8以及式(4)所示的非對(duì)稱的窗函數(shù)Hhigh.asym(k)是在k空間的負(fù)側(cè)設(shè)置非采樣區(qū)域的情況下的圖以及式。因此,在k空間的正側(cè)設(shè)置非采樣區(qū)域的情況下,極性反轉(zhuǎn)。[0095]若利用對(duì)稱或非對(duì)稱的窗函數(shù)Hh_(k)的第二k空間數(shù)據(jù)S’(k)的高頻成分的修正完成,則能夠通過FT生成MR圖像數(shù)據(jù)SeOT(r)作為實(shí)空間數(shù)據(jù)。這里生成的MR圖像數(shù)據(jù)SeOT(r)成為具有與以往的MoFIR法中進(jìn)行了相位修正以及高頻成分的修正的MR圖像數(shù)據(jù)同等的畫質(zhì)的圖像數(shù)據(jù)。MR圖像數(shù)據(jù)Scot(r)能夠如式(5-1)或式(5-2)那樣表示。[0096]Scor(r)=FT{S(k)XHhighsym(k)}(5-1)[0097]Scor(r)=FT{S(k)XHhigh.asym(k)}(5-2)[0098]另外,還能通過對(duì)振幅圖像數(shù)據(jù)StlHd(r)執(zhí)行與將振幅圖像數(shù)據(jù)Sew(r)變換為第二k空間數(shù)據(jù)S’(k)并在進(jìn)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理等價(jià)的實(shí)空間中的卷積處理,由此實(shí)現(xiàn)不伴隨相位修正的AFI處理。即,還能夠不在k空間而是在實(shí)空間中執(zhí)行AFI處理。[0099]圖9是說明圖2所示的AFI處理部41A中執(zhí)行的包含實(shí)空間中的卷積處理的AFI處理的方法的圖。另外,圖9中,對(duì)于與圖3所示的處理同等的處理省略說明。[0100]在實(shí)空間中執(zhí)行卷積處理的情況下,如圖9所示,利用將k空間中的濾波用的對(duì)稱或非對(duì)稱的窗函數(shù)Hh_(k)進(jìn)行FT而得到的實(shí)空間的函數(shù),執(zhí)行振幅圖像數(shù)據(jù)Cr)的卷積處理即可。即,能夠利用與k空間用的窗函數(shù)Hh_(k)等價(jià)的實(shí)空間中的函數(shù)執(zhí)行卷積處理。[0101]另外,實(shí)空間中的濾波函數(shù)FT(Hh_)成為解析式。因此,能夠通過事先將k空間用的窗函數(shù)Hh_(k)離散化并進(jìn)行FT,將實(shí)空間中的濾波函數(shù)FT(Hh_)表(table)化。由此,能夠不需要實(shí)空間中的濾波函數(shù)FT(Hh_)的計(jì)算。即,能夠通過對(duì)表進(jìn)行參照求出濾波函數(shù)FT(Hhomo)0[0102]利用了實(shí)空間中的濾波函數(shù)FT(Hh_)的卷積處理的MR圖像數(shù)據(jù)SeOT(r)的計(jì)算能夠通過式(6-1)或式(6-2)執(zhí)行。[0103]Scor(r)=Convolve(Sj(r),Hhighsym(r)}=Convolve[S,(r),FT{Hhighsym(k)}](6-1)[0104]Scor(r)=Convolve{S,(r),Hhigh.asym(r)}=Convolve[S,(r),FT(Hhighasym(k)}](6-2)[0105]其中,式(6-1)及式(6-2)中Convolve表示{}內(nèi)的函數(shù)間的卷積處理的算子(operator)ο[0106]進(jìn)而,還能夠?qū)⑹雇ㄟ^圖3或圖9所示的處理而生成的MR圖像數(shù)據(jù)Scot(r)中殘留的相位誤差降低的循環(huán)處理作為附加的后處理來執(zhí)行。[0107]圖10是說明后續(xù)于圖3所示的數(shù)據(jù)處理而執(zhí)行使相位誤差降低的循環(huán)處理的情況下的AFI處理的方法的圖。[0108]循環(huán)處理能夠由如下處理構(gòu)成,即:對(duì)MR圖像數(shù)據(jù)ScotCr)取實(shí)空間數(shù)據(jù)的實(shí)部的REAL化處理;使REAL化處理后的實(shí)部的相位向與相位的修正方向相反的方向移位的處理;將使相位向與相位的修正方向相反的方向移位的處理后的實(shí)空間數(shù)據(jù)被進(jìn)行IFT而得到的k空間數(shù)據(jù)中的采樣區(qū)域(-K。^k^Kmax)所對(duì)應(yīng)的數(shù)據(jù)置換為采樣區(qū)域的MR數(shù)據(jù)Sorig(k)的處理;以及將置換后的k空間數(shù)據(jù)FT為實(shí)空間數(shù)據(jù)并使相位向相位的修正方向移位從而取得更新后的MR圖像數(shù)據(jù)的處理。該循環(huán)處理能夠執(zhí)行一次或多次。在執(zhí)行多次循環(huán)處理的情況下,成為使MR圖像數(shù)據(jù)ScotCr)的虛部接近于O的收斂計(jì)算。[0109]另外,如圖10中例示的那樣,利用采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域中的相位分布ΦlOTt,執(zhí)行使實(shí)部的相位向與相位的修正方向相反的方向移位的逆相位修正處理,另一方面,利用與米樣區(qū)域?qū)?yīng)的相位分布執(zhí)行使相位向相位的修正方向移位的相位修正處理,已經(jīng)確認(rèn)到這樣能夠改善畫質(zhì)。因此,逆相位修正處理成為對(duì)實(shí)空間數(shù)據(jù)乘以exp(J^low)的處理,相位修正處理成為對(duì)實(shí)空間數(shù)據(jù)乘以expCj^whole)的處理。[0110]低頻區(qū)域中的相位分布Φ-能夠根據(jù)將低頻側(cè)的區(qū)域中的MR數(shù)據(jù)SOTig(k)進(jìn)行變換而得到的實(shí)空間數(shù)據(jù)S1ot(r)來求出。具體而言,對(duì)MR數(shù)據(jù)SOTig(k)乘以窗函數(shù)H1ot(k),來提取低頻區(qū)域中的對(duì)稱部分,通過對(duì)提取出的對(duì)稱部分進(jìn)行FT而能夠生成與對(duì)稱部分對(duì)應(yīng)的實(shí)空間數(shù)據(jù)S1ot(r)。這樣,通過式(7),能夠求出低頻區(qū)域中的相位分布Φ1Μ。[0111]exp{-jΦlow(χ)}=Slow*(r)/Isiow(r)|(7)[0112]其中,式(7)中*表示復(fù)共軛。[0113]同樣地,與采樣區(qū)域?qū)?yīng)的相位分布也能夠根據(jù)將采樣區(qū)域中的MR數(shù)據(jù)SOTig(k)進(jìn)行變換而得到的實(shí)空間數(shù)據(jù)來求出。具體而言,對(duì)MR數(shù)據(jù)SOTig(k)乘以0-filling用的窗函數(shù)Hwh-(k),提取采樣區(qū)域中的數(shù)據(jù),通過對(duì)提取出的數(shù)據(jù)進(jìn)行FT而能夠生成振幅圖像數(shù)據(jù)Sew(r)。這樣,通過式(8),能夠求出采樣區(qū)域中的相位分布ΦΛ&。[0114]exp{-jOwhole(X)}=S0_filled*(r)/|S0_filled(r)(8)[0115]其中,式(8)中*表示復(fù)共軛。[0116]將對(duì)逆相位修正處理后的實(shí)空間數(shù)據(jù)進(jìn)行IFT而得到的k空間數(shù)據(jù)中的與采樣區(qū)域?qū)?yīng)的數(shù)據(jù)置換為采樣區(qū)域中的MR數(shù)據(jù)SOTig(k)的處理,能夠通過利用將采樣區(qū)域中的數(shù)據(jù)切取的數(shù)據(jù)切取函數(shù)Hmwge(k)的運(yùn)算來進(jìn)行。具體而言,能夠通過將對(duì)MR數(shù)據(jù)SOTig(k)乘以數(shù)據(jù)切取函數(shù)Hmwge(k)而切取出的k空間數(shù)據(jù)、和對(duì)將逆相位修正處理后的實(shí)空間數(shù)據(jù)進(jìn)行IFT而得到的k空間數(shù)據(jù)乘以函數(shù){l-HMge(k)}從而切取出的k空間數(shù)據(jù)進(jìn)行加法運(yùn)算的處理來執(zhí)行置換處理。[0117]接著,說明磁共振成像裝置20的動(dòng)作以及作用。[0118]首先,預(yù)先在診視床37上設(shè)置被檢體P,在通過靜磁場電源26而勵(lì)磁后的靜磁場用磁體21的攝像區(qū)域中形成靜磁場。此外,根據(jù)通過勻場用的預(yù)掃描(pre-scan)收集的數(shù)據(jù),從勻場線圈電源28向勻場線圈22供給規(guī)定的電流,在攝像區(qū)域中形成的靜磁場被均勻化。[0119]并且,序列控制器31按照在攝像條件設(shè)定部40中設(shè)定的AFI用的脈沖序列來驅(qū)動(dòng)梯度磁場電源27、發(fā)送器29以及接收器30,由此,在設(shè)置了被檢體P的攝像區(qū)域中形成梯度磁場,并且從RF線圈24產(chǎn)生RF信號(hào)。[0120]因此,被檢體P的內(nèi)部的由磁共振產(chǎn)生的MR信號(hào)通過RF線圈24而被接收并被提供給接收器30。接收器30從RF線圈24接受MR信號(hào),生成作為數(shù)字?jǐn)?shù)據(jù)的MR信號(hào)的原始數(shù)據(jù)。接收器30將生成的原始數(shù)據(jù)提供給序列控制器31。序列控制器31將原始數(shù)據(jù)提供給數(shù)據(jù)處理部41,數(shù)據(jù)處理部41將原始數(shù)據(jù)配置在形成于k空間數(shù)據(jù)庫43的k空間中。[0121]另外,在k空間中配置的原始數(shù)據(jù)通過AFI用的攝像條件而被收集,所以成為非對(duì)稱的采樣區(qū)域-K。寫k寫Kmax的MR數(shù)據(jù)SOTig(k)。[0122]接著,數(shù)據(jù)處理部41利用包含針對(duì)從k空間取得的k空間數(shù)據(jù)進(jìn)行的FT在內(nèi)的圖像重構(gòu)處理,生成MR圖像數(shù)據(jù)。具體而言,AFI處理部41A執(zhí)行圖3或圖9所示那樣的不伴隨相位修正處理的AFI處理。或者,在AFI處理部41A中執(zhí)行圖10所示的將包含相位修正處理的循環(huán)處理作為后處理的AFI處理。由此,能夠以短的數(shù)據(jù)處理時(shí)間,生成與從對(duì)稱地采樣而得的k空間數(shù)據(jù)生成的MR圖像數(shù)據(jù)大致同等的MR圖像數(shù)據(jù)Scot(r)。[0123]并且,在數(shù)據(jù)處理部41中,在執(zhí)行了針對(duì)MR圖像數(shù)據(jù)Scot(r)的必要的圖像處理之后,能夠顯示在顯示裝置34上或保存在圖像數(shù)據(jù)庫44中。[0124]即,以上那樣的磁共振成像裝置20,將對(duì)按照AFI用的攝像條件而收集的非對(duì)稱的MR數(shù)據(jù)執(zhí)行Ο-filling而得到的振幅圖像數(shù)據(jù)SQ_filled(r),作為以往的MoFIR法的相位修正后的數(shù)據(jù)來處理。[0125]因此,根據(jù)磁共振成像裝置20,在AFI處理中不需要以往所必須的相位修正處理。此外,由于不需要相位修正處理,所以即使不收集相位修正用的原始數(shù)據(jù)或復(fù)數(shù)圖像數(shù)據(jù),也能生成具有與MoFIR法的畫質(zhì)同等畫質(zhì)的MR圖像數(shù)據(jù)。此外,能夠與以往相位修正所需要的計(jì)算時(shí)間相應(yīng)地縮短圖像重構(gòu)處理時(shí)間。即,磁共振成像裝置20中,能夠維持與MoFIR法同等的精度并高速地重構(gòu)圖像數(shù)據(jù)。[0126]圖11是將通過圖3所示的AFI處理生成的ID仿真數(shù)據(jù)與通過以往的AFI處理生成的ID仿真數(shù)據(jù)進(jìn)行比較的圖。[0127]圖11中,各橫軸表示I軸方向的位置X,各縱軸表示位置X的圖像的相對(duì)信號(hào)強(qiáng)度S。此外,圖11(A)是對(duì)非對(duì)稱的k空間數(shù)據(jù)僅執(zhí)行o-filling而得到的仿真數(shù)據(jù),圖11(B)是通過以往的MoFIR從非對(duì)稱的k空間數(shù)據(jù)重構(gòu)的仿真數(shù)據(jù),圖11(C)是利用對(duì)稱的窗函數(shù)Hhigh.sym(k)以圖3所示的順序從非對(duì)稱的k空間數(shù)據(jù)重構(gòu)的仿真數(shù)據(jù),圖11(D)是利用非對(duì)稱的窗函數(shù)Hhigh.asym(k)以圖3所示的順序從非對(duì)稱的k空間數(shù)據(jù)重構(gòu)的仿真數(shù)據(jù)。[0128]根據(jù)圖11,能夠確認(rèn)到,用本發(fā)明的AFI處理方法重構(gòu)的(C)以及(D)所示的仿真數(shù)據(jù)的模糊相比于僅執(zhí)行0-filling而得到的(A)所示的仿真數(shù)據(jù)得到了改善。此外,參照?qǐng)D11的(B)、(C)以及(D),能夠確認(rèn)到,無論是采用對(duì)稱的窗函數(shù)Hhigksyni(k)還是采用非對(duì)稱的窗函數(shù)Hhigh.asyni(k),根據(jù)本發(fā)明的AFI處理方法,相比于通過MoFIR重構(gòu)的仿真數(shù)據(jù),都特別能夠得到不劣化的數(shù)據(jù)。[0129]進(jìn)而,利用ID仿真確認(rèn)到,通過執(zhí)行圖10所示那樣的使相位誤差降低的循環(huán)處理,與不執(zhí)行循環(huán)處理的情況相比,均方根誤差(RMSE:rootmeansquareerror)變小。[0130]以上對(duì)特定的實(shí)施方式進(jìn)行了記載,但所記載的實(shí)施方式只不過是一例,并不限定發(fā)明的范圍。這里記載的新的方法以及裝置能夠以各種其他形態(tài)具體化。此外,這里記載的方法以及裝置的形態(tài)中,在不脫離發(fā)明的主旨的范圍內(nèi),能夠進(jìn)行各種省略、置換以及變更。隨附的權(quán)利要求的范圍以及其等同物包含在發(fā)明的范圍及主旨中,包含這樣的各種形態(tài)以及變形例。[0131]例如,上述的實(shí)施方式中,主要對(duì)于僅在具有kx軸及ky軸的2D的k空間(kx,ky)的kx方向上進(jìn)行非對(duì)稱的k空間數(shù)據(jù)的采樣的例子進(jìn)行了描述,但是在kx方向及ky方向中的I個(gè)方向或2個(gè)方向上進(jìn)行非對(duì)稱的k空間數(shù)據(jù)的采樣的情況下,也能用相同的方法進(jìn)行不伴隨相位修正處理的AFI處理。進(jìn)而,在具有kx軸、ky軸以及kz軸的3D的k空間(kx,ky,kz)中,在I方向、2方向或3方向上進(jìn)行非對(duì)稱的k空間數(shù)據(jù)的采樣的情況下,也能以相同的方法進(jìn)行不伴隨相位修正處理的AFI處理。[0132]在對(duì)2個(gè)方向以上的方向進(jìn)行非對(duì)稱的采樣的情況下,對(duì)于非對(duì)稱的多個(gè)采樣方向,分別將ID的窗函數(shù)按不同方向通過直積來生成,使生成的多個(gè)ID的窗函數(shù)相互相乘,由此能夠生成2D以上的窗函數(shù)。更具體而言,如式(9-1)及式(9-2)所示那樣分別對(duì)于kx軸方向、ky軸方向以及kz軸方向生成ID的0-filling用的窗函數(shù)Hwhau(k)以及用于修正高頻成分的劣化的窗函數(shù)Hh_(k),并使其相互相乘即可。[0133]Hwhole(kx,ky,kz)=Hwhole(kx)XHwhole(kx)XHwhole(kx)(9-1)[0134]Hh_(kx,ky,kz)=Hh_(kx)XHhomo(ky)XHhomo(kz)(9-2)[0135]此外,對(duì)應(yīng)于非對(duì)稱的采樣方向以1D、2D或3D進(jìn)行在AFI處理中最初執(zhí)行的FT以及IFT即可。另一方面,通過對(duì)全部的方向執(zhí)行最后的FT而能夠得到AFI處理后的圖像數(shù)據(jù)?!緳?quán)利要求】1.一種磁共振成像裝置,具備:數(shù)據(jù)收集單元,從被檢體收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的磁共振數(shù)據(jù);以及圖像生成單元,根據(jù)向上述磁共振數(shù)據(jù)的非采樣區(qū)域進(jìn)行的零填充后的第一k空間數(shù)據(jù),生成實(shí)空間的振幅圖像數(shù)據(jù),執(zhí)行將上述振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并在進(jìn)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理、或者利用將上述濾波用的窗函數(shù)進(jìn)行變換而得到的實(shí)空間的函數(shù)進(jìn)行的上述振幅圖像數(shù)據(jù)的卷積處理,由此生成磁共振圖像數(shù)據(jù)。2.如權(quán)利要求1記載的磁共振成像裝置,將用于執(zhí)行針對(duì)上述第二k空間數(shù)據(jù)的上述濾波的窗函數(shù)作為如下窗函數(shù),該窗函數(shù)設(shè)定為:與上述非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,上述非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分和非采樣區(qū)域構(gòu)成的高頻側(cè)的區(qū)域的增益更大。3.如權(quán)利要求2記載的磁共振成像裝置,將用于執(zhí)行針對(duì)第二k空間數(shù)據(jù)的上述濾波的窗函數(shù)作為如下窗函數(shù),該窗函數(shù)設(shè)定為:上述高頻側(cè)的區(qū)域的增益為上述低頻側(cè)的區(qū)域的增益的2倍以下。4.如權(quán)利要求1記載的磁共振成像裝置,將用于執(zhí)行針對(duì)上述第二k空間數(shù)據(jù)的上述濾波的窗函數(shù)作為如下窗函數(shù),該窗函數(shù)設(shè)定為:與上述非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,上述非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分構(gòu)成的高頻側(cè)的區(qū)域的增益更大。5.如權(quán)利要求4記載的磁共振成像裝置,將用于執(zhí)行針對(duì)上述第二k空間數(shù)據(jù)的上述濾波的窗函數(shù)作為如下窗函數(shù),該窗函數(shù)設(shè)定為:上述高頻側(cè)的區(qū)域的增益為上述低頻側(cè)的區(qū)域的增益的2倍以上4倍以下,并且上述非采樣區(qū)域的增益被設(shè)定為零。6.如權(quán)利要求2記載的磁共振成像裝置,作為上述窗函數(shù),使用平滑變化的函數(shù)。7.如權(quán)利要求1記載的磁共振成像裝置,將用于執(zhí)行上述卷積處理的上述實(shí)空間的函數(shù)作為將如下窗函數(shù)進(jìn)行變換而得到的函數(shù),該窗函數(shù)設(shè)定為:與上述k空間中的非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,上述非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分和非采樣區(qū)域構(gòu)成的高頻側(cè)的區(qū)域的增益更大。8.如權(quán)利要求1記載的磁共振成像裝置,將用于執(zhí)行上述卷積處理的上述實(shí)空間的函數(shù)作為將如下窗函數(shù)進(jìn)行變換而得到的函數(shù),該窗函數(shù)設(shè)定為:與上述k空間中的非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域相比較,上述非對(duì)稱的采樣區(qū)域中的非對(duì)稱的采樣部分構(gòu)成的高頻側(cè)的區(qū)域的增益更大。9.如權(quán)利要求1記載的磁共振成像裝置,上述圖像生成單元執(zhí)行一次或多次如下處理:對(duì)上述磁共振圖像數(shù)據(jù)取實(shí)空間數(shù)據(jù)的實(shí)部的處理;使上述實(shí)部的相位向與相位的修正方向相反的方向移位的處理;將使上述相位向與相位的修正方向相反的方向移位的處理后的實(shí)空間數(shù)據(jù)被變換而得到的k空間數(shù)據(jù)中的上述采樣區(qū)域所對(duì)應(yīng)的數(shù)據(jù),置換為上述采樣區(qū)域中的上述磁共振數(shù)據(jù)的處理;以及將上述置換后的k空間數(shù)據(jù)變換為實(shí)空間數(shù)據(jù)并使相位向相位的修正方向移位,由此取得更新后的上述磁共振圖像數(shù)據(jù)的處理。10.如權(quán)利要求9記載的磁共振成像裝置,上述圖像生成單元構(gòu)成為,利用根據(jù)將上述非對(duì)稱的采樣區(qū)域中的對(duì)稱的采樣部分構(gòu)成的低頻側(cè)的區(qū)域中的上述磁共振數(shù)據(jù)進(jìn)行變換而得到的實(shí)空間數(shù)據(jù)而求出的上述低頻側(cè)的區(qū)域中的相位分布,執(zhí)行使上述相位向與相位的修正方向相反的方向移位的處理,另一方面,利用根據(jù)將上述非對(duì)稱的采樣區(qū)域中的上述磁共振數(shù)據(jù)進(jìn)行變換而得到的實(shí)空間數(shù)據(jù)而求出的上述采樣區(qū)域所對(duì)應(yīng)的相位分布,執(zhí)行使上述相位向相位的修正方向移位的處理。11.一種磁共振成像方法,具有:從被檢體收集在k空間中在波數(shù)方向上非對(duì)稱的采樣區(qū)域所對(duì)應(yīng)的磁共振數(shù)據(jù)的步驟;以及根據(jù)向上述磁共振數(shù)據(jù)的非采樣區(qū)域進(jìn)行的零填充后的第一k空間數(shù)據(jù)生成實(shí)空間的振幅圖像數(shù)據(jù),執(zhí)行將上述振幅圖像數(shù)據(jù)變換為第二k空間數(shù)據(jù)并在進(jìn)行了濾波后變換為實(shí)空間數(shù)據(jù)的數(shù)據(jù)處理、或者利用將上述濾波用的窗函數(shù)進(jìn)行變換而得到的實(shí)空間的函數(shù)進(jìn)行的上述振幅圖像數(shù)據(jù)的卷積處理,由此生成磁共振圖像數(shù)據(jù)的步驟。【文檔編號(hào)】A61B5/055GK103747728SQ201380002760【公開日】2014年4月23日申請(qǐng)日期:2013年7月29日優(yōu)先權(quán)日:2012年8月4日【發(fā)明者】木村德典申請(qǐng)人:株式會(huì)社東芝,東芝醫(yī)療系統(tǒng)株式會(huì)社