步驟S103，根據(jù)第一子空間11和第二子空間12的對(duì)比度需求，確定除所述公共區(qū)間201之外的各個(gè)區(qū)間序列采集時(shí)每個(gè)相位編碼線所在的回波鏈位置。確定方法可通過如下迭代步驟進(jìn)行:
[0053](a)從除去所述公共區(qū)間之外所剩余的K空間Kl中，選取離原始K空間中心最近的區(qū)間Sn開始分配；
[0054](b)將回波鏈位置滿足min (En-E)的回波鏈數(shù)據(jù)填補(bǔ)到區(qū)間Sn，其中E為當(dāng)前對(duì)比度下有效TE所在的回波鏈位置，En為待確定的回波鏈位置；
[0055](c)從所述Kl空間中去除上所述Sn區(qū)間，作為當(dāng)前剩余的K空間，再根據(jù)前述步驟從剩余的回波鏈位置中合適的數(shù)據(jù)填補(bǔ)進(jìn)相應(yīng)的區(qū)間，直至所述原始K空間均已進(jìn)行填補(bǔ)。
[0056]在本實(shí)施例中，將除去所述公共區(qū)間201之后的原始K空間I記為K1，將所述Kl等分為4份(等分份數(shù)是以原回波鏈長(zhǎng)度5減去一個(gè)公共區(qū)間獲得)，從距離原始K空間I中心最近的區(qū)間202開始分配。對(duì)于第一子空間11，其有效TE所在回波鏈位置為第I個(gè)回波鏈211，在當(dāng)前剩余的4個(gè)回波鏈位置中，根據(jù)上述確定方法，11-11 < I 2-11 < I 4-11 < I 5-11，因此，所述區(qū)間202位置處應(yīng)填補(bǔ)回波鏈位置為第I個(gè)的回波鏈211 ;對(duì)于第二子空間12，其有效TE所在回波鏈位置為第5個(gè)回波鏈215，在當(dāng)前剩余的4個(gè)回波鏈位置中，根據(jù)上述確定方法，I 5-5 I < I 4-E I < I 2-5 I < 11-5 |所述區(qū)間202位置處應(yīng)填補(bǔ)回波鏈位置為第5個(gè)的回波鏈215。依次類推，對(duì)于第一子空間11，區(qū)間203位置處應(yīng)填補(bǔ)回波鏈位置為第2個(gè)的回波鏈212，此時(shí)碰巧區(qū)間204與區(qū)間205距離原始K空間I的中心距離一樣，因此，所述區(qū)間204與所述區(qū)間205填補(bǔ)的先后順序可隨意決定，在本實(shí)施例中，我們先就所述區(qū)間204進(jìn)行填補(bǔ)，而后對(duì)所述區(qū)間205進(jìn)行填補(bǔ)，因此，所述區(qū)間204位置處應(yīng)填補(bǔ)回波鏈位置為第4個(gè)的回波鏈214，所述區(qū)間205位置處應(yīng)填補(bǔ)回波鏈位置為第5個(gè)的回波鏈215 ;對(duì)于第二子空間12，區(qū)間203位置處應(yīng)填補(bǔ)回波鏈位置為第4個(gè)的回波鏈214，對(duì)于區(qū)間204和區(qū)間205，在第二子空間12中，我們?nèi)匀幌染退鰠^(qū)間204進(jìn)行填補(bǔ)，而后對(duì)所述區(qū)間205進(jìn)行填補(bǔ)，因此，所述區(qū)間204位置處應(yīng)填補(bǔ)回波鏈位置為第2個(gè)的回波鏈212，所述區(qū)間205位置處應(yīng)填補(bǔ)回波鏈位置為第I個(gè)的回波鏈211。
[0057]根據(jù)上述確定好的回波鏈位置需求，進(jìn)行數(shù)據(jù)的采集填補(bǔ)，如此可獲得各子空間序列數(shù)據(jù)。
[0058]然后再執(zhí)行步驟S104，以所述公共區(qū)間201的數(shù)據(jù)為校準(zhǔn)數(shù)據(jù)，利用磁共振并行采集K空間數(shù)據(jù)的方法，計(jì)算出并行采集2倍(此處的倍數(shù)為所需對(duì)比度的個(gè)數(shù))加速時(shí)所需的合并系數(shù)。由于并行采集重建方法非本專利所要求保護(hù)范圍，且已有不少并行采集重建相關(guān)現(xiàn)有技術(shù)，因此，此處不再對(duì)并行采集方法進(jìn)行贅述。
[0059]進(jìn)而再執(zhí)行步驟S105，將所述公共區(qū)間201的數(shù)據(jù)作為第一子空間11和第二子空間12的共享數(shù)據(jù)，分別填補(bǔ)在第一子空間11和第二子空間12相應(yīng)的相位編碼位置處，再根據(jù)所述合并系數(shù)對(duì)第一子空間11和第二子空間12中相對(duì)于原始K空間I中欠采集的數(shù)據(jù)進(jìn)行填補(bǔ)，以分別獲得滿采樣的第一子空間11和第二子空間12。
[0060]最后，執(zhí)行步驟S106，對(duì)所述滿采樣的第一子空間11和第二子空間12進(jìn)行圖像重建。此處的重建圖像重建方法為常規(guī)方法，故也不在贅述。
[0061]對(duì)應(yīng)于上述磁共振多對(duì)比度圖像重建方法，本實(shí)施例還提供了一種多對(duì)比度圖像重建系統(tǒng)。圖3是本實(shí)施例提供的多對(duì)比度圖像重建系統(tǒng)的結(jié)構(gòu)示意圖，如圖3所示，所述磁共振多對(duì)比度圖像重建系統(tǒng)包括:
[0062]子空間劃分單元301，用于根據(jù)所需對(duì)比度數(shù)量將原始K空間I劃分成多個(gè)子空間；
[0063]數(shù)據(jù)采集單元302，用于采集原始K空間I的數(shù)據(jù)信號(hào)；
[0064]數(shù)據(jù)分配單元303，用于確定各子空問所需分配的數(shù)據(jù)及分配方式，將所述數(shù)據(jù)采集單元302采集到的數(shù)據(jù)按前述確定好的分配方式分配至各子空間的對(duì)應(yīng)位置；
[0065]系數(shù)計(jì)算單元304，用于計(jì)算并行采集時(shí)所需的合并系數(shù)；根據(jù)所得合并系數(shù)，利用數(shù)據(jù)采集單元302對(duì)各子空間中數(shù)據(jù)缺失的位置進(jìn)行并行采集填補(bǔ)，獲得滿采樣子空間；
[0066]圖像重建單元305，用于將各個(gè)滿采樣子空間所得到的數(shù)據(jù)分別進(jìn)行重建處理，變換至圖像域，得到各對(duì)比度的重建圖像。
[0067]本實(shí)施例中，所述數(shù)據(jù)分配單元303還包括:
[0068]第一確定單元(圖中未示出)，用于確定公共區(qū)間所在位置；
[0069]第二確定單元(圖中未示出)，用于在除去公共區(qū)間外，根據(jù)各子空間的對(duì)比度需求，分別確定其余各相位編碼線所在的回波鏈位置，將各數(shù)據(jù)分配至相應(yīng)位置。
[0070]本發(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í)施例所作的任何簡(jiǎn)單修改、等同變化及修飾，均屬于本發(fā)明技術(shù)方案的保護(hù)范圍。
【主權(quán)項(xiàng)】
1.一種磁共振多對(duì)比度圖像重建方法，其特征在于，包括: 將原始K空間劃分成多個(gè)子空間，所述子空間的數(shù)量由所需對(duì)比度的數(shù)量決定； 確定原始K空間的公共區(qū)間及相應(yīng)填補(bǔ)的回波鏈位置； 確定除所述公共區(qū)間之外的各個(gè)區(qū)間序列采集時(shí)每個(gè)相位編碼線所在的回波鏈位置，根據(jù)所述回波鏈的位置需求進(jìn)行所采集數(shù)據(jù)的分配填補(bǔ)； 根據(jù)磁共振并行采集圖像重建的方法，計(jì)算并行采集時(shí)所需的合并系數(shù)； 將所述公共區(qū)間的數(shù)據(jù)作為各子空間的共享數(shù)據(jù)，分別填補(bǔ)在各子空間相應(yīng)的相位編碼位置，再根據(jù)所述合并系數(shù)對(duì)各子空間中相對(duì)于原始K空間中欠采集的數(shù)據(jù)進(jìn)行填補(bǔ)，以獲得多個(gè)不同的滿采樣子空間； 對(duì)所述各個(gè)滿采樣子空間進(jìn)行圖像重建。
2.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:所述公共區(qū)間所采集數(shù)據(jù)的回波鏈位置滿足min Σ (EP1-印)2，其中EPi為子空間有效TE所在的回波鏈位置，ep為公共區(qū)間所采集回波鏈位置。
3.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:按照對(duì)診斷對(duì)比度重要性的排列，以最重要對(duì)比度下采集所述回波鏈ep所在位置的區(qū)間作為所述公共區(qū)間。
4.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:與距離原始K空間中心最近的長(zhǎng)度等于原始K空間相位編碼長(zhǎng)度與所采集回波鏈長(zhǎng)度的比值的區(qū)間的相鄰區(qū)間為所述公共區(qū)間。
5.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:所述子空間在所述公共區(qū)間處只能填充相同的序列回波位置的數(shù)據(jù)。
6.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:除所述公共區(qū)間之外的每個(gè)子空間中各個(gè)區(qū)間序列采集時(shí)每個(gè)相位編碼線所在的回波鏈位置的確定通過如下迭代步驟進(jìn)行: 從除去所述公共區(qū)間之外所剩余的K空間Kl中，選取離原始K空間中心最近的區(qū)間Sn開始分配； 將回波鏈位置滿足min (En-E)的回波鏈數(shù)據(jù)填補(bǔ)到區(qū)間Sn，其中E為當(dāng)前某個(gè)子空間對(duì)比度下有效TE所在的回波鏈位置，En為待選的回波鏈位置； 從所述Kl空間中去除上所述Sn區(qū)間，作為當(dāng)前剩余的K空間，再根據(jù)前述步驟從剩余的回波鏈位置中合適的數(shù)據(jù)填補(bǔ)進(jìn)相應(yīng)的區(qū)間，直至所述原始K空間均已進(jìn)行填補(bǔ)。
7.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:所述子空間的數(shù)據(jù)采集，遵循第i個(gè)相位編碼線由第j=mod(i，n)+l個(gè)子空間進(jìn)行數(shù)據(jù)采集填補(bǔ)的方式，其中η代表所述子空間的個(gè)數(shù)，mod(i, η)代表i除以η的余數(shù)。
8.根據(jù)權(quán)利要求1所述的磁共振多對(duì)比度圖像重建方法，其特征在于:在計(jì)算所述合并系數(shù)時(shí)，以所述各子空間在所述公共區(qū)間的數(shù)據(jù)作為校準(zhǔn)數(shù)據(jù)。
9.一種磁共振多對(duì)比度圖像重建系統(tǒng)，其特征在于，包括: 子空間劃分單元，用于根據(jù)所需對(duì)比度數(shù)量將原始K空間劃分成多個(gè)子空間； 數(shù)據(jù)采集單元，用于采集原始K空間的數(shù)據(jù)信號(hào)； 數(shù)據(jù)分配單元，用于確定各子空間所需分配的數(shù)據(jù)及分配方式，將所述數(shù)據(jù)采集單元采集到的相應(yīng)數(shù)據(jù)分配至各子空間的對(duì)應(yīng)位置； 系數(shù)計(jì)算單元，用于計(jì)算并行采集時(shí)所需的合并系數(shù)，根據(jù)所得合并系數(shù)，利用所述數(shù)據(jù)采集單元對(duì)各子空間所數(shù)據(jù)缺失的位置進(jìn)行并行采集填補(bǔ)，獲得滿采樣子空間； 圖像重建單元，用于將各個(gè)滿采樣子空間所得到的數(shù)據(jù)分別進(jìn)行重建處理，變換至圖像域，得到各對(duì)比度的重建圖像。
10.根據(jù)權(quán)利要求9所述的磁共振多對(duì)比度圖像重建系統(tǒng)，其特征在于，所述數(shù)據(jù)分配單元還包括: 第一確定單元，用于確定所述公共區(qū)間所在位置； 第二確定單元，用于在除去所述公共區(qū)間外，根據(jù)各子空間的對(duì)比度需求，分別確定其余各相位編碼線所在的回波鏈位置，將各數(shù)據(jù)分配至相應(yīng)位置。
【專利摘要】一種磁共振多對(duì)比度圖像重建方法及系統(tǒng)，包括如下步驟：根據(jù)所需對(duì)比度數(shù)量將原始K空間劃分成多個(gè)子空間；確定原始K空間的公共區(qū)間及應(yīng)填補(bǔ)的回波鏈位置；確定除所述公共區(qū)間之外的各個(gè)區(qū)間應(yīng)采集的每個(gè)相位編碼線所在的回波鏈位置，根據(jù)此位置需求進(jìn)行數(shù)據(jù)的分配；計(jì)算并行采集時(shí)所需的合并系數(shù)；根據(jù)所述合并系數(shù)對(duì)各子空間中相對(duì)于原始K空間中欠采集的數(shù)據(jù)進(jìn)行填補(bǔ)，以獲得多個(gè)不同的滿采樣子空間；對(duì)所述各個(gè)滿采樣子空間進(jìn)行圖像重建。本發(fā)明通過組合回波鏈在K空間的填充方式，可實(shí)現(xiàn)在只采單個(gè)K空間的情況下，即可獲取多種對(duì)比度圖像，大大節(jié)省了多對(duì)比度圖像的獲取時(shí)間，提高診斷效率。
【IPC分類】G01R33-561
【公開號(hào)】CN104635189
【申請(qǐng)?zhí)枴緾N201310562638
【發(fā)明人】翟人寬
【申請(qǐng)人】上海聯(lián)影醫(yī)療科技有限公司
【公開日】2015年5月20日
【申請(qǐng)日】2013年11月13日