流程圖和/或方框圖中的流程和/或方框的結(jié)合。可提供運(yùn)些計算機(jī)程序 指令到通用計算機(jī)、專用計算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器W產(chǎn) 生一個機(jī)器，使得通過計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實 現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或多個方框中指定的功能的裝置。
[0186] 運(yùn)些計算機(jī)程序指令也可存儲在能引導(dǎo)計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備W特 定方式工作的計算機(jī)可讀存儲器中，使得存儲在該計算機(jī)可讀存儲器中的指令產(chǎn)生包括指 令裝置的制造品，該指令裝置實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一個方框或 多個方框中指定的功能。
[0187] 運(yùn)些計算機(jī)程序指令也可裝載到計算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計 算機(jī)或其他可編程設(shè)備上執(zhí)行一系列操作步驟W產(chǎn)生計算機(jī)實現(xiàn)的處理，從而在計算機(jī)或 其他可編程設(shè)備上執(zhí)行的指令提供用于實現(xiàn)在流程圖一個流程或多個流程和/或方框圖一 個方框或多個方框中指定的功能的步驟。
[0188] W上所述的具體實施例，對本發(fā)明的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳 細(xì)說明，所應(yīng)理解的是，W上所述僅為本發(fā)明的具體實施例而已，并不用于限定本發(fā)明的保 護(hù)范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本 發(fā)明的保護(hù)范圍之內(nèi)。
【主權(quán)項】
1. 一種磁共振溫度成像方法，其特征在于，所述方法包括： 基于高加速并行成像混疊控制CAIPRINHA方法，確定包含多個單激發(fā)脈沖的復(fù)合激發(fā) 脈沖； 在多個溫度成像時刻，分別使用所述復(fù)合激發(fā)脈沖同時激發(fā)一檢測對象的多層成像區(qū) 域，并分別采集所述多層成像區(qū)域的混疊信號數(shù)據(jù)； 根據(jù)所述混疊信號數(shù)據(jù)，通過并行成像算法計算得到在所述溫度成像時刻各層所述成 像區(qū)域的磁共振圖像數(shù)據(jù)； 根據(jù)所述溫度成像時刻的每層所述成像區(qū)域的磁共振圖像數(shù)據(jù)，通過溫度算法計算得 到所述溫度成像時刻的每層所述成像區(qū)域的溫度圖像數(shù)據(jù)； 其中，所述多個溫度成像時刻包含一溫度變化前的初始溫度時刻，所述溫度圖像數(shù)據(jù) 包含溫度變化后的所述溫度成像時刻相對于所述初始溫度時刻的溫度差。2. 如權(quán)利要求1所述的磁共振溫度成像方法，其特征在于，基于高加速并行成像混疊控 制CAIPRINHA方法，確定包含多個單激發(fā)脈沖的復(fù)合激發(fā)脈沖，包括： 利用CAIPRINHA方法確定各所述單激發(fā)脈沖的初始相位； 根據(jù)拉莫方程計算得到所述單激發(fā)脈沖的激發(fā)頻率； 根據(jù)所述多層成像區(qū)域的層數(shù)、所述單激發(fā)脈沖的初始相位、所述激發(fā)頻率及一設(shè)定 脈沖類型確定所述復(fù)合激發(fā)脈沖。3. 如權(quán)利要求1所述的磁共振溫度成像方法，其特征在于，采集所述多層成像區(qū)域的混 疊信號數(shù)據(jù)，包括： 根據(jù)所述復(fù)合激發(fā)脈沖，結(jié)合一選層梯度，獲取所述多層成像區(qū)域的回波信號數(shù)據(jù)； 根據(jù)所述回波信號數(shù)據(jù)，結(jié)合一相位編碼梯度和一讀出梯度，計算所述多層成像區(qū)域 的初始信號數(shù)據(jù)； 對所述初始信號數(shù)據(jù)進(jìn)行模數(shù)轉(zhuǎn)換處理，得到所述混疊信號數(shù)據(jù)。4. 如權(quán)利要求1所述的磁共振溫度成像方法，其特征在于，根據(jù)所述混疊信號數(shù)據(jù)，通 過并行成像算法計算得到各層所述成像區(qū)域的磁共振圖像數(shù)據(jù)，包括： 使用梯度回波GRE序列逐層采集所述多層成像區(qū)域的設(shè)定分辨率圖像； 利用所有所述設(shè)定分辨率圖像估算用于采集所述混疊信號數(shù)據(jù)的相控陣線圈的敏感 度空間分布數(shù)據(jù)； 利用所述敏感度空間分布數(shù)據(jù)和所述混疊信號數(shù)據(jù)，計算得到各層所述成像區(qū)域的磁 共振圖像數(shù)據(jù)。5. 如權(quán)利要求1所述的磁共振溫度成像方法，其特征在于，根據(jù)所述混疊信號數(shù)據(jù)，通 過并行成像算法計算得到各層所述成像區(qū)域的磁共振圖像數(shù)據(jù)，包括： 利用相控陣線圈中的多個線圈，以奈奎斯特脈沖頻率對所述成像區(qū)域的K空間中心數(shù) 據(jù)進(jìn)行全采樣，以獲取各所述線圈的K空間自校準(zhǔn)線ACS數(shù)據(jù)； 針對每個所述線圈，利用所述ACS數(shù)據(jù)和所述線圈采集到的所述混疊信號所對應(yīng)的K空 間數(shù)據(jù)，估算所述線圈的權(quán)重系數(shù)； 針對每個所述線圈，利用所述線圈的權(quán)重系數(shù)和所述線圈采集到的所述混疊信號所對 應(yīng)的K空間數(shù)據(jù)，擬合得到所述線圈未采集到的各層所述成像區(qū)域所對應(yīng)的K空間數(shù)據(jù)； 針對每個所述線圈，利用擬合得到的K空間數(shù)據(jù)、所述ACS數(shù)據(jù)以及采集到的所述混疊 信號數(shù)據(jù)所對應(yīng)的K空間數(shù)據(jù)，合并得到所述線圈所對應(yīng)的各層所述成像區(qū)域的完整K空間 數(shù)據(jù)； 對所述完整Κ空間數(shù)據(jù)進(jìn)行逆傅里葉變換，得所述線圈采集的各層所述成像區(qū)域的圖 像數(shù)據(jù)； 利用Β1權(quán)重線圈組合方法，將所有所述線圈采集的圖像數(shù)據(jù)合成為解混疊后的所述多 層成像區(qū)域的磁共振圖像數(shù)據(jù)，并獲取各層所述成像區(qū)域的磁共振圖像數(shù)據(jù)。6. 如權(quán)利要求1所述的磁共振溫度成像方法，其特征在于，根據(jù)所述溫度成像時刻的每 層所述成像區(qū)域的磁共振圖像數(shù)據(jù)，通過溫度算法計算得到所述溫度成像時刻的每層所述 成像區(qū)域的溫度圖像數(shù)據(jù)，包括： 根據(jù)所述初始溫度時刻的每層所述成像區(qū)域的磁共振圖像數(shù)據(jù)的相位和溫度變化后 的所述溫度成像時刻的每層所述成像區(qū)域的磁共振圖像數(shù)據(jù)的相位，計算得到相位差； 根據(jù)所述相位差，通過質(zhì)子共振頻率偏移PRFS方法，計算得到在所述溫度差，以得到溫 度變化后的所述溫度成像時刻的每層所述成像區(qū)域的溫度圖像數(shù)據(jù)。7. 如權(quán)利要求2所述的磁共振溫度成像方法，其特征在于，所述初始相位為：其中，η為所述成像區(qū)域的層數(shù)序號，Ν為所述多層成像區(qū)域的總層數(shù)，Μ為所述復(fù)合激 發(fā)脈沖的重復(fù)周期序號。8. 如權(quán)利要求4所述的磁共振溫度成像方法，其特征在于，所述混疊信號數(shù)據(jù)為：其中，x，y為像素位置，i為線圈序號，CKxj)為所述混疊信號;j為所述成像區(qū)域的層數(shù) 序號，N為所述多層成像區(qū)域的總層數(shù)，F(xiàn)OV為視野范圍為考慮CAIPRINHA方法 后第j層成像區(qū)域的圖像在上述混疊信號所對應(yīng)圖像內(nèi)的面內(nèi)位移，為線圈i在第j層所 述成像區(qū)域的敏感度空間分布數(shù)據(jù);P』為第j層所述成像區(qū)域受單激發(fā)脈沖激發(fā)后發(fā)出的 信號。9. 如權(quán)利要求6所述的磁共振溫度成像方法，其特征在于，所述溫度差與所述相位差之 間的關(guān)系為：其中，Δρ為所述相位差，γ為旋磁比，α為溫度系數(shù)，Bo為主磁場強(qiáng)度，△ T為所述溫度 差，TE為回波時間。10. 如權(quán)利要求5所述的磁共振溫度成像方法，其特征在于， 估算所述線圈的權(quán)重系數(shù)的關(guān)系式為：擬合得到所述線圈未采集到的各層成像區(qū)域所對應(yīng)的K空間數(shù)據(jù)的關(guān)系式為：其中，表示第j個線圈的ACS數(shù)據(jù)，ky表示相位編碼方向，△ ky為相位編碼的間隔，m 表示ACS數(shù)據(jù)《^^在冗空間的塊中的相對位置，L為相控陣線圈中的線圈總數(shù)，1為線圈的序 號，Nb為相位編碼方向塊的總個數(shù)，b為相位編碼塊的序號，11〇，13，1，111)為需要計算的線圈 的權(quán)重系數(shù)，R為加速因子，SKky-bRAky)為第1個線圈所鄰近ACS線的已采集到的K空間數(shù) 據(jù)，&為第j個線圈未采集的數(shù)據(jù)。11. 一種磁共振溫度成像系統(tǒng)，其特征在于，所述系統(tǒng)包括： 復(fù)合激發(fā)脈沖生成單元，用于基于高加速并行成像混疊控制CAIPRINHA方法，確定包含 多個單激發(fā)脈沖的復(fù)合激發(fā)脈沖； 混疊信號數(shù)據(jù)生成單元，用于在多個溫度成像時刻，分別使用所述復(fù)合激發(fā)脈沖同時 激發(fā)一檢測對象的多層成像區(qū)域，并分別采集所述多層成像區(qū)域的混疊信號數(shù)據(jù)； 磁共振圖像數(shù)據(jù)重建單元，用于根據(jù)所述混疊信號數(shù)據(jù)，通過并行成像算法計算得到 在所述溫度成像時刻各層所述成像區(qū)域的磁共振圖像數(shù)據(jù)； 溫度圖像數(shù)據(jù)生成單元，用于根據(jù)所述溫度成像時刻的每層所述成像區(qū)域的磁共振圖 像數(shù)據(jù)，通過溫度算法計算得到所述溫度成像時刻的每層所述成像區(qū)域的溫度圖像數(shù)據(jù)； 其中，所述多個溫度成像時刻包含一溫度變化前的初始溫度時刻，所述溫度圖像數(shù)據(jù) 包含溫度變化后的所述溫度成像時刻相對于所述初始溫度時刻的溫度差。12. 如權(quán)利要求11所述的磁共振溫度成像系統(tǒng)，其特征在于，所述復(fù)合激發(fā)脈沖生成單 元包括： 初始相位生成模塊，用于利用CAIPRINHA方法確定各所述單激發(fā)脈沖的初始相位； 激發(fā)頻率生成模塊，用于根據(jù)拉莫方程計算得到所述單激發(fā)脈沖的激發(fā)頻率； 復(fù)合激發(fā)脈沖生成模塊，用于根據(jù)所述多層成像區(qū)域的層數(shù)、所述單激發(fā)脈沖的初始 相位、所述激發(fā)頻率及一設(shè)定脈沖類型確定所述復(fù)合激發(fā)脈沖。13. 如權(quán)利要求11所述的磁共振溫度成像系統(tǒng)，其特征在于，所述混疊信號數(shù)據(jù)生成單 元包括： 回波信號數(shù)據(jù)生成模塊，用于根據(jù)所述復(fù)合激發(fā)脈沖，結(jié)合一選層梯度，獲取所述多層 成像區(qū)域的回波信號數(shù)據(jù)； 初始信號數(shù)據(jù)生成模塊，用于根據(jù)所述回波信號數(shù)據(jù)，結(jié)合一相位編碼梯度和一讀出 梯度，計算所述多層成像區(qū)域的初始信號數(shù)據(jù)； 混疊信號數(shù)據(jù)生成模塊，用于對所述初始信號數(shù)據(jù)進(jìn)行模數(shù)轉(zhuǎn)換處理，得到所述混疊 信號數(shù)據(jù)。14. 如權(quán)利要求11所述的磁共振溫度成像系統(tǒng)，其特征在于，所述磁共振圖像數(shù)據(jù)重建 單元包括： 設(shè)定分辨率圖像生成模塊，用于使用梯度回波GRE序列逐層采集所述多層成像區(qū)域的 設(shè)定分辨率圖像； 敏感度空間分布數(shù)據(jù)生成模塊，用于利用所有所述設(shè)定分辨率圖像估算用于采集所述 混疊信號數(shù)據(jù)的相控陣線圈的敏感度空間分布數(shù)據(jù)； 第一磁共振圖像數(shù)據(jù)生成模塊，用于利用所述敏感度空間分布數(shù)據(jù)和所述混疊信號數(shù) 據(jù)，計算得到各層所述成像區(qū)域的磁共振圖像數(shù)據(jù)。15. 如權(quán)利要求11所述的磁共振溫度成像系統(tǒng)，其特征在于，所述磁共振圖像數(shù)據(jù)重建 單元包括： ACS數(shù)據(jù)生成模塊，用于利用相控陣線圈中的多個線圈，以奈奎斯特脈沖頻率對所述成 像區(qū)域的K空間中心數(shù)據(jù)進(jìn)行全采樣，以獲取各所述線圈的K空間自校準(zhǔn)線ACS數(shù)據(jù)； 線圈權(quán)重系數(shù)生成模塊，用于針對每個所述線圈，利用所述ACS數(shù)據(jù)和所述線圈采集到 的所述混疊信號所對應(yīng)的K空間數(shù)據(jù)，估算所述線圈的權(quán)重系數(shù)； 未采集K空間數(shù)據(jù)生成模塊，用于針對每個所述線圈，利用所述線圈的權(quán)重系數(shù)和所述 線圈采集到的所述混疊信號所對應(yīng)的K空間數(shù)據(jù)，擬合得到所述線圈未采集到的各層所述 成像區(qū)域所對應(yīng)的K空間數(shù)據(jù)； 完整K空間數(shù)據(jù)生成模塊，用于針對每個所述線圈，利用擬合得到的K空間數(shù)據(jù)、所述 ACS數(shù)據(jù)以及采集到的所述混疊信號數(shù)據(jù)所對應(yīng)的K空間數(shù)據(jù)，合并得到所述線圈所對應(yīng)的 各層所述成像區(qū)域的完整K空間數(shù)據(jù)； 初始圖像數(shù)據(jù)生成模塊，用于對所述完整K空間數(shù)據(jù)進(jìn)行逆傅里葉變換，得所述線圈采 集的各層所述成像區(qū)域的圖像數(shù)據(jù)； 第二磁共振圖像數(shù)據(jù)生成模塊，用于利用B1權(quán)重線圈組合方法，將所有所述線圈采集 的圖像數(shù)據(jù)合成為解混疊后的所述多層成像區(qū)域的磁共振圖像數(shù)據(jù)，并獲取各層所述成像 區(qū)域的磁共振圖像數(shù)據(jù)。16. 如權(quán)利要求11所述的磁共振溫度成像系統(tǒng)，其特征在于，所述溫度圖像數(shù)據(jù)生成單 元包括： 相位差生成模塊，用于根據(jù)所述初始溫度時刻的每層所述成像區(qū)域的磁共振圖像數(shù)據(jù) 的相位和溫度變化后的所述溫度成像時刻的每層所述成像區(qū)域的磁共振圖像數(shù)據(jù)的相位， 計算得到相位差； 溫度圖像數(shù)據(jù)生成模塊，用于根據(jù)所述相位差，通過質(zhì)子共振頻率偏移PRFS方法，計算 得到在所述溫度差，以得到溫度變化后的所述溫度成像時刻的每層所述成像區(qū)域的溫度圖 像數(shù)據(jù)。
【專利摘要】本發(fā)明提供了一種磁共振溫度成像方法及系統(tǒng)，該方法包括：基于高加速并行成像混疊控制方法，確定包含多個單激發(fā)脈沖的復(fù)合激發(fā)脈沖；在多個溫度成像時刻，分別使用該復(fù)合激發(fā)脈沖同時激發(fā)多層成像區(qū)域，并分別采集該多層成像區(qū)域的混疊信號數(shù)據(jù)；根據(jù)該混疊信號數(shù)據(jù)，通過并行成像算法計算得到在該溫度成像時刻各層該成像區(qū)域的磁共振圖像數(shù)據(jù)；根據(jù)該溫度成像時刻的每層該成像區(qū)域的磁共振圖像數(shù)據(jù)，通過溫度算法計算得到所述溫度成像時刻的每層所述成像區(qū)域的溫度圖像數(shù)據(jù)；多個溫度成像時刻包含一溫度變化前的初始溫度時刻，該溫度圖像數(shù)據(jù)包含溫度變化后的該溫度成像時刻相對于該初始溫度時刻的溫度差。本發(fā)明能夠提高溫度成像效率。
【IPC分類】A61B5/055, A61B5/01
【公開號】CN105662357
【申請?zhí)枴緾N201511007976
【發(fā)明人】郭文莎, 劉新, 鄒超, 帖長軍
【申請人】中國科學(xué)院深圳先進(jìn)技術(shù)研究院
【公開日】2016年6月15日
【申請日】2015年12月29日