本發(fā)明涉及磁共振成像(MRI)技術(shù)領(lǐng)域，特別是一種擴(kuò)散加權(quán)的回波平面成像(DW-EPI)的圖像畸變校正方法及系統(tǒng)，以及磁共振成像設(shè)備。

背景技術(shù)：
擴(kuò)散成像是一種非常具有臨床價值的功能磁共振成像技術(shù)，它可以探測水分子在人體組織中擴(kuò)散的生物學(xué)特性，從而檢測許多常規(guī)醫(yī)學(xué)影像手段不能發(fā)現(xiàn)的病變。擴(kuò)散成像通常使用回波平面成像序列(EPI)來實現(xiàn)，EPI序列成像速度非?？?，但是容易受到各種偽影尤其是圖像畸變的影響。擴(kuò)散加權(quán)成像(DWI)是擴(kuò)散成像的一種比較重要的成像方法，其通常使用DW-EPI序列來實現(xiàn)，DW-EPI序列是通過在EPI序列中加入擴(kuò)散梯度來實現(xiàn)擴(kuò)散加權(quán)的。DW-EPI序列的擴(kuò)散梯度會在成像空間中引入渦流場，而渦流場這種額外磁場是圖像畸變產(chǎn)生的原因之一。由于在擴(kuò)散加權(quán)成像中通常需要采集若干具有不同擴(kuò)散方向和擴(kuò)散權(quán)重的單張圖像(即DW-EPI圖像)，要用到不同的擴(kuò)散梯度，引起不同的渦流場分布，因此這些單張圖像的畸變程度和類型并不完全相同。其中，圖像畸變的外觀既取決于梯度的幅度，即擴(kuò)散權(quán)重，又取決于梯度的方向。擴(kuò)散權(quán)重的大小大多通過所謂的“b值”來確定，b值為零時的圖像一般不存在畸變，b值不為零時的圖像存在有圖像畸變?？梢?，這些具有不同b值和不同擴(kuò)散方向的單張圖像是不重合的，如果將這些單張圖像不經(jīng)校正的相互組合和計算以便產(chǎn)生用于診斷所需的各向同性DWI以及表觀擴(kuò)散系數(shù)(ADC)圖等，就會造成DWI的模糊，以及ADC圖計算的錯誤等。抑制渦流場的手段有使用自屏蔽的梯度線圈和梯度脈沖的預(yù)加強(qiáng)。自屏蔽線圈包括一個主線圈和一個屏蔽線圈，屏蔽線圈在成像時產(chǎn)生的磁場的方向與主線圈相反，經(jīng)過特殊設(shè)計，使得這一對線圈在成像區(qū)域內(nèi)產(chǎn)生需要的梯度磁場，而在成像區(qū)域外的磁場盡可能接近零，這就從根本上減小了渦流場。而預(yù)加強(qiáng)技術(shù)是通過改變梯度波形的形狀，來抵消渦流場的影響。通過使用這些方法，可以有效減小渦流場的影響，但是并不能完全消除。而DW-EPI序列對渦流場非常敏感，為了提高成像的質(zhì)量，還需要使用特別方法來處理。

技術(shù)實現(xiàn)要素：
有鑒于此，本發(fā)明一方面提出了一種DW-EPI的圖像畸變校正方法，另一方面提出了一種DW-EPI的圖像畸變校正系統(tǒng)，用以提高磁共振成像的圖像質(zhì)量。本發(fā)明還要提出一種磁共振成像設(shè)備。因此，本發(fā)明提供了一種擴(kuò)散加權(quán)的回波平面成像的圖像畸變校正方法，包括：利用組合序列采集得到復(fù)數(shù)個具有不同設(shè)定b值和不同設(shè)定擴(kuò)散方向的帶標(biāo)記圖像，其中所述組合序列包括擴(kuò)散加權(quán)的回波平面成像序列以及在其之前施加的標(biāo)記序列；利用擴(kuò)散加權(quán)的回波平面成像序列采集得到復(fù)數(shù)個具有與所述標(biāo)記圖像相一致的b值及擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像；將各組具有相同b值和相同擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像和帶標(biāo)記圖像相減，得到對應(yīng)各擴(kuò)散加權(quán)回波平面圖像的標(biāo)記圖像；其中，b值為零的標(biāo)記圖像為參考標(biāo)記圖像，b值不為零的標(biāo)記圖像為畸變標(biāo)記圖像；對b值不為零的每個擴(kuò)散加權(quán)回波平面圖像，利用與該擴(kuò)散加權(quán)回波平面圖像具有相同b值和相同擴(kuò)散方向的畸變標(biāo)記圖像和該擴(kuò)散加權(quán)回波平面圖像對應(yīng)的參考標(biāo)記圖像，計算得到對應(yīng)該擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，并利用所述伸縮系數(shù)和位移系數(shù)對所述擴(kuò)散加權(quán)回波平面圖像進(jìn)行校正?？蛇x地，計算得到對應(yīng)該擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)包括：從所述畸變標(biāo)記圖像中選取復(fù)數(shù)個畸變圖像數(shù)據(jù)列，對每個選取的畸變圖像數(shù)據(jù)列，在所述參考標(biāo)記圖像中選取位置對應(yīng)的參考圖像數(shù)據(jù)列；或者，從所述參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列，對每個選取的參考圖像數(shù)據(jù)列，在所述畸變標(biāo)記圖像中選取位置對應(yīng)的畸變圖像數(shù)據(jù)列；對當(dāng)前選取的參考圖像數(shù)據(jù)列及其對應(yīng)的畸變圖像數(shù)據(jù)列，從所述參考圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為參考標(biāo)記頻率；從所述畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為畸變標(biāo)記頻率；將所述畸變標(biāo)記頻率與所述參考標(biāo)記頻率的比值作為所述畸變圖像數(shù)據(jù)列相對于所述參考圖像數(shù)據(jù)列的伸縮系數(shù)；利用所述伸縮系數(shù)對所述畸變圖像數(shù)據(jù)列進(jìn)行校正，將校正后的畸變圖像數(shù)據(jù)列在設(shè)定方向上按照設(shè)定位移間隔依次進(jìn)行位移，并計算每次位移后的畸變圖像數(shù)據(jù)列與參考圖像數(shù)據(jù)列之間的相關(guān)值，將最大相關(guān)值對應(yīng)的位移量作為所述畸變圖像數(shù)據(jù)列相對于所述參考圖像數(shù)據(jù)列的位移系數(shù)；在計算得到復(fù)數(shù)個畸變圖像數(shù)據(jù)列所對應(yīng)的伸縮系數(shù)和位移系數(shù)后，分別對復(fù)數(shù)個伸縮系數(shù)和復(fù)數(shù)個位移系數(shù)進(jìn)行線性擬合和外推，得到對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，作為對應(yīng)所述擴(kuò)散加權(quán)回波平面圖像的各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)?？蛇x地，所述從畸變標(biāo)記圖像中選取復(fù)數(shù)個畸變圖像數(shù)據(jù)列包括：依次提取所述畸變標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一畸變圖像數(shù)據(jù)列。所述從參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列包括：依次提取所述參考標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一參考圖像數(shù)據(jù)列。本發(fā)明還提供了一種擴(kuò)散加權(quán)的回波平面成像的圖像畸變校正系統(tǒng)，包括：一個圖像采集單元，用于利用組合序列采集得到復(fù)數(shù)個具有不同設(shè)定b值和不同設(shè)定擴(kuò)散方向的帶標(biāo)記圖像，其中所述組合序列包括擴(kuò)散加權(quán)的回波平面成像序列以及在其之前施加的標(biāo)記序列；利用擴(kuò)散加權(quán)的回波平面成像序列采集得到復(fù)數(shù)個具有與所述標(biāo)記圖像相一致的b值及擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像；一個標(biāo)記圖像生成單元，用于將各組具有相同b值和相同擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像和帶標(biāo)記圖像相減，得到對應(yīng)各擴(kuò)散加權(quán)回波平面圖像的標(biāo)記圖像；其中，b值為零的標(biāo)記圖像為參考標(biāo)記圖像，b值不為零的標(biāo)記圖像為畸變標(biāo)記圖像；一個畸變系數(shù)計算單元，用于對b值不為零的各擴(kuò)散加權(quán)回波平面圖像中的當(dāng)前擴(kuò)散加權(quán)回波平面圖像，利用與所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像具有相同b值相同擴(kuò)散方向的畸變標(biāo)記圖像和所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像對應(yīng)的參考標(biāo)記圖像，計算得到對應(yīng)所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)；一個圖像校正單元，用于利用所述畸變系數(shù)計算單元得到的對應(yīng)當(dāng)前擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，對所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像進(jìn)行校正?？蛇x地，所述畸變系數(shù)計算單元包括：一個圖像數(shù)據(jù)列選取子單元，用于從所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像對應(yīng)的畸變標(biāo)記圖像中選取復(fù)數(shù)個畸變圖像數(shù)據(jù)列，對每個選取的畸變圖像數(shù)據(jù)列，在所述參考標(biāo)記圖像中選取位置對應(yīng)的參考圖像數(shù)據(jù)列；或者，從所述參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列，對每個選取的參考圖像數(shù)據(jù)列，在所述畸變標(biāo)記圖像中選取位置對應(yīng)的畸變圖像數(shù)據(jù)列；一個標(biāo)記頻率確定子單元，用于對所述圖像數(shù)據(jù)列選取子單元選取的當(dāng)前參考圖像數(shù)據(jù)列及其對應(yīng)的當(dāng)前畸變圖像數(shù)據(jù)列，從所述當(dāng)前參考圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為參考標(biāo)記頻率，從所述當(dāng)前畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為畸變標(biāo)記頻率；一個伸縮系數(shù)計算子單元，用于將所述標(biāo)記頻率確定子單元確定的所述畸變標(biāo)記頻率與所述參考標(biāo)記頻率的比值作為所述當(dāng)前畸變圖像數(shù)據(jù)列相對于所述當(dāng)前參考圖像數(shù)據(jù)列的伸縮系數(shù)；一個位移系數(shù)計算子單元，用于利用所述伸縮系數(shù)計算子單元計算的伸縮系數(shù)對所述當(dāng)前畸變圖像數(shù)據(jù)列進(jìn)行校正，將校正后的當(dāng)前畸變圖像數(shù)據(jù)列在設(shè)定方向上按照設(shè)定位移間隔依次進(jìn)行位移，并計算每次位移后的當(dāng)前畸變圖像數(shù)據(jù)列與當(dāng)前參考圖像數(shù)據(jù)列之間的相關(guān)值，將最大相關(guān)值對應(yīng)的位移量作為所述當(dāng)前畸變圖像數(shù)據(jù)列相對于所述當(dāng)前參考圖像數(shù)據(jù)列的位移系數(shù)；一個畸變系數(shù)擬合子單元，用于在所述伸縮系數(shù)計算子單元和所述位移系數(shù)計算子單元計算得到復(fù)數(shù)個畸變圖像數(shù)據(jù)列所對應(yīng)的伸縮系數(shù)和位移系數(shù)后，分別對復(fù)數(shù)個伸縮系數(shù)和復(fù)數(shù)個位移系數(shù)進(jìn)行線性擬合和外推，得到對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，作為對應(yīng)所述擴(kuò)散加權(quán)回波平面圖像的各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)。可選地，所述圖像數(shù)據(jù)列選取子單元依次提取所述畸變標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一畸變圖像數(shù)據(jù)列?；蛘撸鰣D像數(shù)據(jù)列選取子單元依次提取所述參考標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一參考圖像數(shù)據(jù)列。可選地，所述系統(tǒng)進(jìn)一步包括：一個序列組合單元，用于在擴(kuò)散加權(quán)的回波平面成像序列之前施加標(biāo)記序列，構(gòu)成組合序列；可選地，所述系統(tǒng)進(jìn)一步包括：一個處理器和一個計算機(jī)可讀存儲介質(zhì)；其中，所述序列組合單元、標(biāo)記圖像生成單元、畸變系數(shù)計算單元和圖像校正單元中的任一個或任意組合為存儲在所述計算機(jī)可讀存儲介質(zhì)中的程序代碼；所述處理器用于執(zhí)行所述計算機(jī)可讀存儲介質(zhì)中的程序代碼。本發(fā)明還提供了一種磁共振成像設(shè)備，包括如上所述的擴(kuò)散加權(quán)的回波平面成像的圖像畸變校正系統(tǒng)。從上述方案中可以看出，本發(fā)明實施例中利用由標(biāo)記序列和DW-EPI序列構(gòu)成的組合序列采集得到帶標(biāo)記圖像，并通過與正常采集圖像相減，得到純標(biāo)記圖像；之后利用與b值不為零的正常采集圖像對應(yīng)的參考標(biāo)記圖像和畸變標(biāo)記圖像，計算該b值不為零的正常采集圖像的畸變系數(shù)，之后利用該畸變系數(shù)對該b值不為零的正常采集圖像進(jìn)行校正，可以有效校正渦流引入的圖像畸變。此外，本發(fā)明實施例中以列為單位，從參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列，從畸變標(biāo)記圖像中選取復(fù)數(shù)個位置對應(yīng)的畸變圖像數(shù)據(jù)列，將渦流場引起的圖像畸的整體位移，剪切和伸縮變形轉(zhuǎn)換為只有位移和伸縮兩種形變，通過計算畸變圖像數(shù)據(jù)列相對于參考圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，確定整個畸變標(biāo)記圖像相對于參考標(biāo)記圖像的變形，簡化了整個計算過程。附圖說明下面將通過參照附圖詳細(xì)描述本發(fā)明的優(yōu)選實施例，使本領(lǐng)域的普通技術(shù)人員更清楚本發(fā)明的上述及其它特征和優(yōu)點，附圖中：圖1a和圖1b為本發(fā)明實施例中一個b值為零的DW-EPI圖像和一個b值不為零的DW-EPI圖像的示意圖。圖2為本發(fā)明實施例中DW-EPI成像的圖像畸變校正方法的示例性流程圖。圖3為本發(fā)明實施例中一個組合序列的示意圖。圖4a和圖4b分別為本發(fā)明實施例中一個DW-EPI圖像和一個帶標(biāo)記圖像的示意圖，圖4c為標(biāo)記圖像的示意圖。圖5為本發(fā)明實施例中根據(jù)參考標(biāo)記圖像和畸變標(biāo)記圖像計算對應(yīng)DW-EPI圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)的方法流程圖。圖6a為本發(fā)明實施例中一個參考標(biāo)記圖像及其對應(yīng)的一個畸變標(biāo)記圖像的示意圖。圖6b為圖6a中虛線所對應(yīng)的圖像數(shù)據(jù)列的一維曲線示意圖。圖6c為本發(fā)明實施例中圖6a所示參考圖像數(shù)據(jù)列的傅里葉變換譜曲線和畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線的示意圖。圖6d為對本發(fā)明實施例中圖6c所示譜曲線中小方框所示區(qū)域內(nèi)的波峰的擬合示意圖。圖7為本發(fā)明實施例中DW-EPI成像的圖像畸變校正系統(tǒng)的結(jié)構(gòu)示意圖。圖8為圖7中所示畸變系數(shù)計算單元的一個結(jié)構(gòu)示意圖。圖9a至圖9g為本發(fā)明示例一中對b＝1000的圖像進(jìn)行畸變校正的示意圖。其中，圖9a圖是b＝0的圖像；圖9b至圖9d是b＝1000的三個不同擴(kuò)散方向上的圖像；圖9e至圖9g分別是對圖9b至圖9d中各擴(kuò)散方向上的圖像校正后的圖像。圖10a和圖10b為本示例中在b值為1000時的各向同性擴(kuò)散加權(quán)圖像的校正示意圖。其中，圖10a為利用原始采集的圖9b至圖9d所示的單張DW-EPI圖像計算得到的各向同性DW-EPI圖像，圖10b為利用圖9e至圖9g所示校正過的單張DW-EPI圖像計算得到各向同性擴(kuò)散加權(quán)圖像。圖11a和圖11b為本示例中基于b值為1000時的校正前和校正后的DW-EPI圖像生成ADC的效果示意圖。圖12a至圖12g為本發(fā)明示例二中對b＝1000的圖像進(jìn)行畸變校正的示意圖。其中，圖12a圖是b＝0的圖像；圖12b至圖12d是b＝1000的三個不同擴(kuò)散方向上的圖像；圖12e至圖12g分別是對圖12b至圖12d中各擴(kuò)散方向上的圖像校正后的圖像。其中，附圖標(biāo)記如下：201-構(gòu)造組合序列202-采集帶標(biāo)記圖像和正常DW-EPI圖像203-獲取標(biāo)記圖像204-根據(jù)參考標(biāo)記圖像和畸變標(biāo)記圖像獲取對應(yīng)的DW-EPI圖像的畸變系數(shù)，并利用該畸變系數(shù)進(jìn)行圖像校正。501-選取畸變圖像數(shù)據(jù)列和參考圖像數(shù)據(jù)列502-根據(jù)所選取的圖像數(shù)據(jù)列計算對應(yīng)的畸變標(biāo)記頻率和參考標(biāo)記頻率503-計算伸縮系數(shù)504-計算位移系數(shù)505-對復(fù)數(shù)個伸縮系數(shù)和復(fù)數(shù)個位移系數(shù)進(jìn)行線性擬合和外推，得到所有圖像數(shù)據(jù)列的畸變系數(shù)701-序列組合單元702-圖像采集單元703-標(biāo)記圖像生成單元704-畸變系數(shù)計算單元705-圖像校正單元801-圖像數(shù)據(jù)列選取子單元802-標(biāo)記頻率確定子單元803-伸縮系數(shù)計算子單元804-位移系數(shù)計算子單元805-畸變系數(shù)擬合子單元具體實施方式下面首先對單張DW-EPI圖像畸變的類型進(jìn)行如下分析：在EPI序列中，額外的磁場ΔB在頻率編碼方向和相位編碼方向上引起的圖像畸變可以分別表達(dá)為式(1)和式(2)：其中，x為頻率編碼方向上的坐標(biāo)，y為相位編碼方向上的坐標(biāo)；Gx為頻率編碼梯度，為相位編碼梯度的平均值；τ為相位編碼梯度的持續(xù)時間；tesp為回波間距。實際應(yīng)用中，由于在頻率編碼方向上的變形通常很小，可以忽略，因此，一般只考慮在相位編碼方向上的變形。DW-EPI序列的擴(kuò)散梯度引入的渦流場即是額外磁場ΔB的一個主要來源，其是圖像畸變產(chǎn)生的原因。令表示渦流場在主磁場B0方向上的分量，則可以如式(3)所示展開為：其中，第一項b0(t)是與空間位置無關(guān)的，通常稱為零階項。第二項為線性項，或者稱為一階項，為矢量坐標(biāo)，分量為x，y，z，相應(yīng)地，矢量梯度的分量為gx、gy、gz，分別表示在x，y，z三個方向上渦流場的梯度。更高階的項通常不予考慮。其中，z為選層梯度方向上的坐標(biāo)。由于線性項中的z分量只會造成選層相位回聚的不完善，使得圖像整體信號強(qiáng)度下降，但是不會造成圖像畸變，因此本發(fā)明實施例中只考慮零階項和一階項中的x，y分量，式(3)可以重寫為式(4)：Be(x，y)＝b0+x·gx+y·gy(4)將式(4)代入式(2)的圖像畸變表達(dá)式中，可得到由渦流場引起的圖像畸變的表達(dá)式：從上式可以看出，零階渦流場會引起圖像整體的位移，沿著頻率編碼方向的一階項會使圖像產(chǎn)生剪切形變，沿著相位編碼方向的一階項使圖像產(chǎn)生拉伸形變。為此，本發(fā)明實施例中，可針對單張DW-EPI的變形類型對其進(jìn)行校正。為了便于實現(xiàn)，本發(fā)明實施例中，對每個單張DW-EPI圖像，以其中的每個單列圖像數(shù)據(jù)為單位，進(jìn)行縮放與位移的校正。令I(lǐng)0(x，y)表示b值為零的DW-EPI圖像，I1(x，y)表示b值不為零有畸變的DW-EPI圖像，分別如圖1a和圖1b所示。其中，x為頻率編碼方向，y為相位編碼方向。用I0(x)和I1(x)分別表示DW-EPI圖像中沿著相位編碼方向的一個圖像數(shù)據(jù)列，分別由圖1a和圖1b中的實線和虛線表示，根據(jù)上文的分析，渦流場引起的圖像畸變有整體位移、剪切和伸縮。如果逐列考慮的話，那么只有位移和伸縮兩種形變，只需把I1(x)進(jìn)行適量的伸縮和平移，就可以與I0(x)對齊：I0(x)＝S(x)·I1(x)+T(x)(6)其中，S(x)，T(x)分別表示第x圖像數(shù)據(jù)列的伸縮系數(shù)和平移系數(shù)。為了計算對應(yīng)每個圖像數(shù)據(jù)列的伸縮系數(shù)和平移系數(shù)，本發(fā)明實施例中考慮利用標(biāo)記圖像來計算該標(biāo)記圖像對應(yīng)的DW-EPI圖像的各圖像數(shù)據(jù)列的伸縮系數(shù)和平移系數(shù)。為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚，以下舉實施例對本發(fā)明進(jìn)一步詳細(xì)說明。圖2為本發(fā)明實施例中DW-EPI成像的圖像畸變校正方法的示例性流程圖。如圖2所示，該方法包括如下步驟：步驟201，在DW-EPI序列之前施加標(biāo)記序列，構(gòu)成組合序列。如圖3所示，圖3為本發(fā)明實施例中一個組合序列的示意圖，其中圖3的左側(cè)為一標(biāo)記序列，右側(cè)為DW-EPI序列。由于變形主要出現(xiàn)在相位編碼方向，因此本實施例中可采用只在相位編碼方向出現(xiàn)直線的標(biāo)記。步驟202，利用所述組合序列采集得到復(fù)數(shù)個具有不同設(shè)定b值和不同設(shè)定擴(kuò)散方向的帶標(biāo)記圖像；利用DW-EPI序列采集得到復(fù)數(shù)個具有與所述標(biāo)記圖像相一致的b值及擴(kuò)散方向的DW-EPI圖像。圖4a和圖4b分別示出了本發(fā)明實施例中一個DW-EPI圖像和一個帶標(biāo)記圖像，且二者具有相同的b值和相同的擴(kuò)散方向。步驟203，將各組具有相同b值和相同擴(kuò)散方向的DW-EPI圖像和帶標(biāo)記圖像相減，得到對應(yīng)各DW-EPI圖像的標(biāo)記圖像。其中，b值為零的標(biāo)記圖像可稱為參考標(biāo)記圖像，b值不為零的標(biāo)記圖像可稱為畸變標(biāo)記圖像。圖4c示出了本發(fā)明實施例中一個標(biāo)記圖像。該標(biāo)記圖像為圖4a所示DW-EPI圖像減去圖4b所示的帶標(biāo)記圖像后得到的標(biāo)記圖像。步驟204，對b值不為零的每個DW-EPI圖像，利用與該DW-EPI圖像具有相同b值相同擴(kuò)散方向的畸變標(biāo)記圖像和該DW-EPI圖像對應(yīng)的參考標(biāo)記圖像，計算得到對應(yīng)該DW-EPI圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，并利用所述伸縮系數(shù)和位移系數(shù)對所述DW-EPI圖像進(jìn)行校正。本步驟中，利用與該DW-EPI圖像具有相同b值相同擴(kuò)散方向的畸變標(biāo)記圖像和該DW-EPI圖像對應(yīng)的參考標(biāo)記圖像，計算得到對應(yīng)該DW-EPI圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)的具體實現(xiàn)方式可有多種。圖5示出了其中一種具體實現(xiàn)方法，如圖5所示，該方法包括如下步驟：步驟501從所述畸變標(biāo)記圖像中選取復(fù)數(shù)個畸變圖像數(shù)據(jù)列，對每個選取的畸變圖像數(shù)據(jù)列，在所述參考標(biāo)記圖像中選取位置對應(yīng)的參考圖像數(shù)據(jù)列。具體實現(xiàn)時，本實施例中無需對畸變標(biāo)記圖像的所有圖像數(shù)據(jù)列都計算畸變系數(shù)，即伸縮系數(shù)和位移系數(shù)，因為如果某一圖像數(shù)據(jù)列中只包含背景噪聲而沒有圖像，這種計算顯然是沒有意義的。另外，靠近邊緣的若干圖像數(shù)據(jù)列也是被排除在外的，因為這些圖像數(shù)據(jù)列中只有少量的圖像，如果用來計算位移系數(shù)和伸縮系數(shù)會引入很大的誤差。因此本步驟中可僅選取若干圖像數(shù)據(jù)列來進(jìn)行畸變系數(shù)的計算。圖6a中示出了本發(fā)明實施例中一個參考標(biāo)記圖像及其對應(yīng)的一個畸變標(biāo)記圖像，其中，上方的圖像為參考標(biāo)記圖像，下方的圖像為畸變標(biāo)記圖像。參考標(biāo)記圖像中虛線所對應(yīng)的位置為選取的一個參考圖像數(shù)據(jù)列，畸變標(biāo)記圖像中虛線所對應(yīng)的位置為選取的一個畸變圖像數(shù)據(jù)列。圖6b為圖6a中虛線所對應(yīng)的圖像數(shù)據(jù)列的一維曲線示意圖。具體實現(xiàn)時，本步驟中可依次提取所述畸變標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，如圖6c所示，圖6c中下方的曲線為圖6a中下方畸變標(biāo)記圖像中虛線所示圖像數(shù)據(jù)列的傅里葉變換譜曲線。若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰(如圖6c中下方曲線中小方框所示區(qū)域內(nèi)的波峰)的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一畸變圖像數(shù)據(jù)列?；蛘撸唧w實現(xiàn)時，本步驟中也可直接選取畸變標(biāo)記圖像中位于中間區(qū)域的若干圖像數(shù)據(jù)列作為畸變圖像數(shù)據(jù)列。當(dāng)然，本步驟中也可從所述參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列，并對每個選取的參考圖像數(shù)據(jù)列，在所述畸變標(biāo)記圖像中選取位置對應(yīng)的畸變圖像數(shù)據(jù)列。具體實現(xiàn)時，同樣可以依次提取所述參考標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一參考圖像數(shù)據(jù)列?；蛘?，直接選取參考標(biāo)記圖像中位于中間區(qū)域的若干圖像數(shù)據(jù)列作為參考圖像數(shù)據(jù)列。步驟502，對當(dāng)前選取的參考圖像數(shù)據(jù)列及其對應(yīng)的畸變圖像數(shù)據(jù)列，從所述參考圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為參考標(biāo)記頻率，從所述畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為畸變標(biāo)記頻率。仍以圖6a中所示的參考圖像數(shù)據(jù)列及其對應(yīng)的畸變圖像數(shù)據(jù)列為例，則參考圖像數(shù)據(jù)列的傅里葉變換譜曲線和畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線可如圖6c所示，圖6c中上方的曲線為圖6a中上方參考標(biāo)記圖像中虛線所示圖像數(shù)據(jù)列的傅里葉變換譜曲線，圖6c中下方的曲線為圖6a中下方畸變標(biāo)記圖像中虛線所示圖像數(shù)據(jù)列的傅里葉變換譜曲線。對應(yīng)的波峰擬合示意圖可如圖6d所示，圖6d中上方的波峰示意圖為對圖6c中上方曲線中小方框所示區(qū)域內(nèi)的波峰的擬合示意圖，相應(yīng)地，波峰的峰值位置為參考標(biāo)記頻率fref(x)；圖6d中下方的波峰示意圖為對圖6c中下方曲線中小方框所示區(qū)域內(nèi)的波峰的擬合示意圖，相應(yīng)地，波峰的峰值位置為畸變標(biāo)記頻率fdistorted(x)。步驟503，將所述畸變標(biāo)記頻率與所述參考標(biāo)記頻率的比值作為所述畸變圖像數(shù)據(jù)列相對于所述參考圖像數(shù)據(jù)列的伸縮系數(shù)。對于圖6a中所示的參考圖像數(shù)據(jù)列及其對應(yīng)的畸變圖像數(shù)據(jù)列，所述畸變圖像數(shù)據(jù)列相對于所述參考圖像數(shù)據(jù)列的伸縮系數(shù)為：S(x)＝fdistorted(x)/fref(x)。步驟504，利用所述伸縮系數(shù)對所述畸變圖像數(shù)據(jù)列進(jìn)行校正，將校正后的畸變圖像數(shù)據(jù)列在設(shè)定方向上按照設(shè)定位移間隔依次進(jìn)行位移，并計算每次位移后的畸變圖像數(shù)據(jù)列與參考圖像數(shù)據(jù)列之間的相關(guān)值，將最大相關(guān)值對應(yīng)的位移量作為所述畸變圖像數(shù)據(jù)列相對于所述參考圖像數(shù)據(jù)列的位移系數(shù)。對應(yīng)于步驟503中的S(x)可得到對應(yīng)的T(x)。步驟505，在計算得到復(fù)數(shù)個畸變圖像數(shù)據(jù)列所對應(yīng)的伸縮系數(shù)和位移系數(shù)后，分別對復(fù)數(shù)個伸縮系數(shù)和復(fù)數(shù)個位移系數(shù)進(jìn)行線性擬合和外推，得到對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，將所述對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)作為對應(yīng)所述擴(kuò)散加權(quán)回波平面圖像的各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)。考慮到具體實現(xiàn)時，由于圖像中必然會存在噪聲，因此畸變系數(shù)的計算必然也會引入一些誤差，所以計算得到的各伸縮系數(shù)和位移系數(shù)不是精確的，如果直接使用這些畸變系數(shù)對DW-EPI圖像相應(yīng)圖像數(shù)據(jù)列進(jìn)行校正的話，效果不會很好?？紤]到渦流場在整個成像空間的分布是線性的，所以各列的位移系數(shù)和伸縮系數(shù)也應(yīng)該滿足線性的關(guān)系，因此可以使用最小二乘法對S(x)和T(x)進(jìn)行擬合，然后外推到圖像中所有的圖像數(shù)據(jù)列，得到擬合和外推之后的結(jié)果Sfit(x)和Tfit(x)，以減小誤差。之后，可按照如下的方法計算出像素偏移圖：PSM(x，y)＝y(tǒng)·Sfit(x)+Tfti(x)(7)得到像素偏移圖后，就可以對所述存在畸變的DW-EPI圖像進(jìn)行校正，得到與b值為0的DW-EPI圖像配準(zhǔn)的圖像。以上對本發(fā)明實施例中的DW-EPI成像的圖像畸變校正方法進(jìn)行了詳細(xì)描述，下面再對本發(fā)明實施例中實施所述方法的一種DW-EPI成像的圖像畸變校正系統(tǒng)進(jìn)行詳細(xì)描述。圖7為本發(fā)明實施例中DW-EPI成像的圖像畸變校正系統(tǒng)的結(jié)構(gòu)示意圖。如圖7所示，該系統(tǒng)可包括：一個序列組合單元701、一個圖像采集單元702、一個標(biāo)記圖像生成單元703、一個畸變系數(shù)計算單元704和一個圖像校正單元705。其中，序列組合單元701用于在擴(kuò)散加權(quán)的回波平面成像序列之前施加標(biāo)記序列，構(gòu)成組合序列。圖像采集單元702用于利用所述組合序列采集得到復(fù)數(shù)個具有不同設(shè)定b值和不同設(shè)定擴(kuò)散方向的帶標(biāo)記圖像；利用擴(kuò)散加權(quán)的回波平面成像序列采集得到復(fù)數(shù)個具有與所述標(biāo)記圖像相一致的b值及擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像。標(biāo)記圖像生成單元703用于將各組具有相同b值和相同擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像和帶標(biāo)記圖像相減，得到對應(yīng)各擴(kuò)散加權(quán)回波平面圖像的標(biāo)記圖像；其中，b值為零的標(biāo)記圖像為參考標(biāo)記圖像，b值不為零的標(biāo)記圖像為畸變標(biāo)記圖像?；兿禂?shù)計算單元704用于對b值不為零的各擴(kuò)散加權(quán)回波平面圖像中的當(dāng)前擴(kuò)散加權(quán)回波平面圖像，利用與所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像具有相同b值相同擴(kuò)散方向的畸變標(biāo)記圖像和所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像對應(yīng)的參考標(biāo)記圖像，計算得到對應(yīng)所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)。圖像校正單元705用于利用所述畸變系數(shù)計算單元得到的對應(yīng)當(dāng)前擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，對所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像進(jìn)行校正。具體實現(xiàn)時，所述畸變系數(shù)計算單元704可有多種具體實現(xiàn)形式，圖8中示出了其中一種。如圖8所示，該畸變系數(shù)計算單元704可包括：一個圖像數(shù)據(jù)列選取子單元801、一個標(biāo)記頻率確定子單元802、一個伸縮系數(shù)計算子單元803、一個位移系數(shù)計算子單元804和一個畸變系數(shù)擬合子單元805。其中，圖像數(shù)據(jù)列選取子單元801用于從所述當(dāng)前擴(kuò)散加權(quán)回波平面圖像對應(yīng)的畸變標(biāo)記圖像中選取復(fù)數(shù)個畸變圖像數(shù)據(jù)列，對每個選取的畸變圖像數(shù)據(jù)列，在所述參考標(biāo)記圖像中選取位置對應(yīng)的參考圖像數(shù)據(jù)列；或者，從所述參考標(biāo)記圖像中選取復(fù)數(shù)個參考圖像數(shù)據(jù)列，對每個選取的參考圖像數(shù)據(jù)列，在所述畸變標(biāo)記圖像中選取位置對應(yīng)的畸變圖像數(shù)據(jù)列。標(biāo)記頻率確定子單元802用于對所述圖像數(shù)據(jù)列選取子單元選取的當(dāng)前參考圖像數(shù)據(jù)列及其對應(yīng)的當(dāng)前畸變圖像數(shù)據(jù)列，從所述當(dāng)前參考圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為參考標(biāo)記頻率，從所述當(dāng)前畸變圖像數(shù)據(jù)列的傅里葉變換譜曲線中確定最高峰的一側(cè)緊鄰的一波峰的位置，對所述波峰進(jìn)行位置擬合，將擬合后的波峰的峰值位置作為畸變標(biāo)記頻率。伸縮系數(shù)計算子單元803用于將所述標(biāo)記頻率確定子單元確定的所述畸變標(biāo)記頻率與所述參考標(biāo)記頻率的比值作為所述當(dāng)前畸變圖像數(shù)據(jù)列相對于所述當(dāng)前參考圖像數(shù)據(jù)列的伸縮系數(shù)。位移系數(shù)計算子單元804用于利用所述伸縮系數(shù)計算子單元計算的伸縮系數(shù)對所述當(dāng)前畸變圖像數(shù)據(jù)列進(jìn)行校正，將校正后的當(dāng)前畸變圖像數(shù)據(jù)列在設(shè)定方向上按照設(shè)定位移間隔依次進(jìn)行位移，并計算每次位移后的當(dāng)前畸變圖像數(shù)據(jù)列與當(dāng)前參考圖像數(shù)據(jù)列之間的相關(guān)值，將最大相關(guān)值對應(yīng)的位移量作為所述當(dāng)前畸變圖像數(shù)據(jù)列相對于所述當(dāng)前參考圖像數(shù)據(jù)列的位移系數(shù)?；兿禂?shù)擬合子單元805用于在所述伸縮系數(shù)計算子單元和所述位移系數(shù)計算子單元計算得到復(fù)數(shù)個畸變圖像數(shù)據(jù)列所對應(yīng)的伸縮系數(shù)和位移系數(shù)后，分別對復(fù)數(shù)個伸縮系數(shù)和復(fù)數(shù)個位移系數(shù)進(jìn)行線性擬合和外推，得到對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，將所述對應(yīng)整個畸變標(biāo)記圖像各畸變圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)作為對應(yīng)所述擴(kuò)散加權(quán)回波平面圖像的各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)。具體實現(xiàn)時，所述圖像數(shù)據(jù)列選取子單元801可以依次提取所述畸變標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一畸變圖像數(shù)據(jù)列；或者，也可直接選取畸變標(biāo)記圖像中位于中間區(qū)域的若干圖像數(shù)據(jù)列作為畸變圖像數(shù)據(jù)列。此外，所述圖像數(shù)據(jù)列選取子單元801也可依次提取所述參考標(biāo)記圖像中的各個圖像數(shù)據(jù)列，對當(dāng)前提取的圖像數(shù)據(jù)列，計算得到所述圖像數(shù)據(jù)列的傅里葉變換譜曲線，若所述譜曲線中最高峰的峰值與所述最高峰一側(cè)的一波峰的峰值之比大于設(shè)定的選取閾值，則將所述當(dāng)前提取的圖像數(shù)據(jù)列選取為一參考圖像數(shù)據(jù)列；或者，直接選取參考標(biāo)記圖像中位于中間區(qū)域的若干圖像數(shù)據(jù)列作為參考圖像數(shù)據(jù)列。具體實現(xiàn)時，圖像采集單元702可以為磁共振成像設(shè)備中用于實現(xiàn)圖像采集的所有元件的組合，例如可包括各梯度線圈及圖像采集元件等。序列組合單元701、標(biāo)記圖像生成單元703、畸變系數(shù)計算單元704和圖像校正單元705中的任一個或任意組合可以為硬件電路，也可以為程序代碼。當(dāng)序列組合單元701、標(biāo)記圖像生成單元703、畸變系數(shù)計算單元704和圖像校正單元705中的任一個或任意組合為程序代碼時，本發(fā)明實施例中的DW-EPI成像的畸變圖像校正系統(tǒng)可進(jìn)一步包括：處理器(圖中未示出)和計算機(jī)可讀存儲介質(zhì)(圖中未示出)。其中，所述處理器用于執(zhí)行所述計算機(jī)可讀存儲介質(zhì)中的程序代碼。所述序列組合單元701、標(biāo)記圖像生成單元703、畸變系數(shù)計算單元704和圖像校正單元705中的任一個或任意組合為存儲在所述計算機(jī)可讀存儲介質(zhì)中的程序代碼。本發(fā)明實施例中的磁共振成像設(shè)備，可包括上述任一具體實現(xiàn)形式的擴(kuò)散加權(quán)的回波平面成像的圖像畸變校正系統(tǒng)。本發(fā)明實施例中所描述的擴(kuò)散加權(quán)的回波平面成像(DW-EPI)的圖像畸變校正技術(shù)，適用于所有基于DW-EPI序列的成像技術(shù)，例如：擴(kuò)散加權(quán)成像(DWI)，表觀擴(kuò)散系數(shù)計算(ADCmapping)，擴(kuò)散張量成像(DTI)，纖維束跟蹤等。下面列舉兩個應(yīng)用本發(fā)明實施例中技術(shù)方案的具體示例。示例一：DW-EPI圖像的成像具體參數(shù)為：1.5T磁共振成像儀，TE/TR＝127/500msec；帶寬：1000Hz；視野(FoV)：200×200mm；成像矩陣：128×128；層厚：1；平均次數(shù)：4；擴(kuò)散成像使用的b值分別為：0，500，1000；擴(kuò)散模式：正交；擴(kuò)散方法：單極性。帶標(biāo)記圖像的成像具體參數(shù)與DW-EPI圖像的成像具體參數(shù)大部分相同，只是平均次數(shù)為1。圖9a至圖9g顯示了對b＝1000的圖像進(jìn)行畸變校正的示意圖。圖9a圖是b＝0的圖像；圖9b至圖9d是b＝1000的三個不同擴(kuò)散方向上的圖像，其中，圖9b為相位編碼方向上的圖像，圖9c為讀出方向上的圖像，圖9d為選層方向上的圖像。圖9e至圖9g分別是對圖9b至圖9d中各擴(kuò)散方向上的圖像校正后的圖像，其中，圖9e為對圖9b進(jìn)行校正后的圖像，圖9f為對圖9c進(jìn)行校正后的圖像，圖9g為對圖9d進(jìn)行校正后的圖像。通過比較校正前和校正后的圖像可以看出，校正后的圖像的邊緣與b＝0的圖像的輪廓線吻合很好，說明本發(fā)明實施例中的方案有效校正了渦流場造成的圖像畸變。圖10a和圖10b示出了本示例中在b值為1000時的各向同性擴(kuò)散加權(quán)圖像，其中，圖10a由原始采集的圖9b至圖9d所示的單張DW-EPI圖像計算得到，圖10b由圖9e至圖9g所示校正過的單張DW-EPI圖像計算得到。通過對各自右下角方框內(nèi)的邊緣區(qū)域的放大比較，可以看出，校正后的圖像可以有效降低各向同性擴(kuò)散加權(quán)圖像的模糊程度。圖11a和圖11b分別示出了本示例中利用原始各向同性擴(kuò)散加權(quán)圖像生成的ADC圖像和利用校正后的各向同性擴(kuò)散加權(quán)圖像生成的ADC圖像，其中，圖11a為利用原始各向同性擴(kuò)散加權(quán)圖像生成的ADC圖像，圖11b為利用校正后的各向同性擴(kuò)散加權(quán)圖像生成的ADC圖像。通過對其下方矩形框中的邊緣區(qū)域的放大比較，可以看出，由于有效校正了渦流引起的圖像畸變，利用校正后的各向同性擴(kuò)散加權(quán)圖像生成的ADC圖像更加精確。示例二DW-EPI圖像的成像具體參數(shù)為：3T磁共振成像儀，TE/TR＝108/5900msec；帶寬：1040Hz；視野(FoV)：235×235mm；成像矩陣：192×192；層厚：25；平均次數(shù)：4；擴(kuò)散成像使用的b值分別為：0，1000；擴(kuò)散模式：正交；擴(kuò)散方法：單極性。帶標(biāo)記圖像的成像具體參數(shù)與DW-EPI圖像的成像具體參數(shù)大部分相同，只是平均次數(shù)為1。圖12a至圖12g顯示了對b＝1000的圖像進(jìn)行畸變校正的結(jié)果。圖12a圖是b＝0的圖像；圖12b至圖12d是b＝1000的三個不同擴(kuò)散方向上的圖像，其中，圖12b為相位編碼方向上的圖像，圖12c為讀出方向上的圖像，圖12d為選層方向上的圖像。圖12e至圖12g分別是對圖12b至圖12d中各擴(kuò)散方向上的圖像校正后的圖像，其中，圖12e為對圖12b進(jìn)行校正后的圖像，圖12f為對圖12c進(jìn)行校正后的圖像，圖12g為對圖12d進(jìn)行校正后的圖像。通過比較校正前和校正后的圖像可以看出，校正后的圖像的邊緣與b＝0的圖像的輪廓線吻合很好，說明本發(fā)明實施例中的方案有效校正了渦流場造成的圖像畸變。通過比較校正前和校正后的圖像可以看出，校正后的圖像的邊緣與b＝0的圖像的輪廓線吻合很好，說明本發(fā)明實施例中的方案有效校正了渦流場造成的圖像畸變。本發(fā)明涉及磁共振成像中擴(kuò)散加權(quán)的回波平面成像，并且公開了一種圖像畸變校正方法，包括：在擴(kuò)散加權(quán)的回波平面成像序列之前施加標(biāo)記序列，構(gòu)成組合序列；利用所述組合序列采集得到復(fù)數(shù)個具有不同設(shè)定b值和不同設(shè)定擴(kuò)散方向的帶標(biāo)記圖像；利用擴(kuò)散加權(quán)的回波平面成像序列采集得到復(fù)數(shù)個具有與所述標(biāo)記圖像相一致的b值及擴(kuò)散方向的擴(kuò)散加權(quán)回波平面圖像；計算得到對應(yīng)該擴(kuò)散加權(quán)回波平面圖像各圖像數(shù)據(jù)列的伸縮系數(shù)和位移系數(shù)，并利用所述伸縮系數(shù)和位移系數(shù)對所述擴(kuò)散加權(quán)回波平面圖像進(jìn)行校正。本發(fā)明還公開了一種擴(kuò)散加權(quán)的回波平面成像的圖像畸變校正系統(tǒng)，以及包括該系統(tǒng)的磁共振成像設(shè)備。本發(fā)明能夠有效校正渦流引入的圖像畸變。以上所述僅為本發(fā)明的較佳實施例而已，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。