本發(fā)明涉及磁共振成像領(lǐng)域,尤其涉及一種磁共振快速成像方法及系統(tǒng)。
背景技術(shù):
作為放射學主要的成像技術(shù)之一,磁共振成像被廣泛地應用在臨床醫(yī)療診斷。因為磁共振成像具有軟組織高對比度,沒有放射危害,可實現(xiàn)多對比機制成像,可選擇任意方位成像,可以進行溫度成像等優(yōu)勢。但是,相對于電子計算機斷層掃描(CT)、超聲掃描等成像技術(shù),磁共振掃描速度要慢很多,尤其是對適合用于做介入手術(shù)的開放式低場磁共振成像系統(tǒng)。
目前,常見的磁共振掃描加速手段主要具有如下缺陷:
1)K空間關(guān)鍵孔技術(shù)。該技術(shù)先采集一個具有全觀數(shù)的參考數(shù)據(jù)集,然后頻繁更新K空間中心低頻部分的數(shù)據(jù),從而加速掃描速度。因為只是低頻的中心數(shù)據(jù)得到更新,所以,犧牲了很大的空間分辨率。對于像介入針這樣的小尺度目標,尤其是針尖,則會存在較大的誤差。
2)壓縮感知技術(shù)。該技術(shù)利用磁共振圖像的稀疏性或者在變換域具有稀疏性,通過擬隨機模式欠采樣K空間數(shù)據(jù),然后非線性重建出圖像。該方法基于嚴謹?shù)臄?shù)學基礎,是完全可以在低采樣率下重建出高質(zhì)量的圖像,只要能夠滿足圖像稀疏性或者圖像變換域具有稀疏性,和K空間數(shù)據(jù)欠采樣具有隨機性兩個要求。但是,因為醫(yī)學磁共振圖像多數(shù)只是在變換域具有稀疏性(少數(shù)比如血管造影圖在圖像域就具有稀疏性),而且非線性重建也需要復雜耗時的數(shù)學計算,因此,難以應用在實際工作中。
技術(shù)實現(xiàn)要素:
本發(fā)明的目的是提供一種磁共振快速成像方法及系統(tǒng),不僅可以保持較高的空間分辨率,還可以加快掃描速度。
本發(fā)明的目的是通過以下技術(shù)方案實現(xiàn)的:
一種磁共振快速成像方法,包括:
對目標視野施加參考掃描,采集K空間原始數(shù)據(jù)中的全部觀數(shù)的數(shù)據(jù),作為K空間滑動完整數(shù)據(jù)集的起始點;
對目標視野施加更新掃描,并結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù),且每次采集的數(shù)據(jù)均包含K空間中高頻成分與低頻成分;
每次更新掃描之后,將K空間滑動完整數(shù)據(jù)集中與更新掃描采集的觀數(shù)數(shù)據(jù)相對應的相位編碼行進行替換更新;
對更新后的K空間滑動完整數(shù)據(jù)集進行二維傅里葉逆變換,并對變換結(jié)果的值作取絕對值和歸一化處理得到圖像矩陣并呈現(xiàn)相應圖像。
所述參數(shù)對包括:更新比與偏移量;
其中,所述更新比,指更新掃描采集的部分觀數(shù)占參考掃描采集的全部觀數(shù)的百分比;
所述偏移量,是用來表示更新掃描采集的部分觀數(shù)相對K空間中心位置的相對位置,定義為更新掃描采集的部分觀數(shù)中被K空間中心線分割的較小部分與較大部分的比值。
所述更新比為25%、33%、40%、45%或者50%;
所述偏移量為1:4、1:3、1:2或者1:1。
所述結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù)包括:
如果是奇數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間上半空間為主的部分觀數(shù)數(shù)據(jù);
如果是偶數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間下半空間為主的部分觀數(shù)數(shù)據(jù)。
所述進行二維傅里葉逆變換時,施加兩次90°相位偏移。
一種磁共振快速成像系統(tǒng),包括:
參考掃描模塊,用于對目標視野施加參考掃描,采集K空間原始數(shù)據(jù)中的全部觀數(shù)的數(shù)據(jù),作為K空間滑動完整數(shù)據(jù)集的起始點;
更新掃描模塊,用于對目標視野施加更新掃描,并結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù),且每次采集的數(shù)據(jù)均包含K空間中高頻成分與低頻成分;
替換模塊,用于每次更新掃描之后,將K空間滑動完整數(shù)據(jù)集中與更新掃描采集的觀數(shù)數(shù)據(jù)相對應的相位編碼行進行替換更新;
成像模塊,用于對更新后的K空間滑動完整數(shù)據(jù)集進行二維傅里葉逆變換,并對變換結(jié)果的值作取絕對值和歸一化處理得到圖像矩陣并呈現(xiàn)相應圖像。
所述參數(shù)對包括:更新比與偏移量;
其中,所述更新比,指更新掃描采集的部分觀數(shù)占參考掃描采集的全部觀數(shù)的百分比;
所述偏移量,是用來表示更新掃描采集的部分觀數(shù)相對K空間中心位置的相對位置,定義為更新掃描采集的部分觀數(shù)被K空間中心線分割的較小部分與較大部分的比值。
所述更新比為25%、33%、40%、45%或者50%;
所述偏移量為1:4、1:3、1:2或者1:1。
所述結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù)包括:
如果是奇數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間上半空間為主的部分觀數(shù)數(shù)據(jù);
如果是偶數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間下半空間為主的部分觀數(shù)數(shù)據(jù)。
所述進行二維傅里葉逆變換時,施加兩次90°相位偏移。
由上述本發(fā)明提供的技術(shù)方案可以看出,每次進行更新掃描時,所采集到的K空間部分觀數(shù)數(shù)據(jù)同時包含高頻和低頻成分,因而可以既保持較高的空間分辨率又可以提高時間分辨率,并加快掃描速度。
附圖說明
為了更清楚地說明本發(fā)明實施例的技術(shù)方案,下面將對實施例描述中所需要使用的附圖作簡單地介紹,顯而易見地,下面描述中的附圖僅僅是本發(fā)明的一些實施例,對于本領(lǐng)域的普通技術(shù)人員來講,在不付出創(chuàng)造性勞動的前提下,還可以根據(jù)這些附圖獲得其他附圖。
圖1為本發(fā)明實施例提供的一種磁共振快速成像方法的流程圖;
圖2為本發(fā)明實施例提供的參考掃描時采集K空間原始數(shù)據(jù)中的全部觀數(shù)數(shù)據(jù)的示意圖;
圖3為本發(fā)明實施例提供的奇數(shù)次更新掃描時采集以K空間上半空間為主的部分觀數(shù)數(shù)據(jù)的示意圖;
圖4為本發(fā)明實施例提供的偶數(shù)次更新掃描時采集以K空間下半空間為主的部分觀數(shù)數(shù)據(jù)的示意圖;
圖5為本發(fā)明實施例提供的一種磁共振快速成像系統(tǒng)的示意圖。
具體實施方式
下面結(jié)合本發(fā)明實施例中的附圖,對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述,顯然,所描述的實施例僅僅是本發(fā)明一部分實施例,而不是全部的實施例?;诒景l(fā)明的實施例,本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例,都屬于本發(fā)明的保護范圍。
圖1為本發(fā)明實施例提供的一種磁共振快速成像方法的流程圖。如圖1所示,其主要包括如下步驟:
步驟11、對目標視野施加參考掃描,采集K空間原始數(shù)據(jù)中的全部觀數(shù)的數(shù)據(jù),作為K空間滑動完整數(shù)據(jù)集的起始點。
本發(fā)明實施例中,參考掃描由脈沖掃描序列實現(xiàn),根據(jù)掃描目標的不同,選擇的脈沖序列會有所不同。而且針對不同的目標視野和序列,選擇的掃描參數(shù)也會有所不同。完整K空間相位編碼數(shù)目和頻率編碼數(shù)目(一次連續(xù)數(shù)據(jù)讀出)由目標視野大小以及分辨率大小具體確定。
如圖2所示,參考掃描時采集K空間原始數(shù)據(jù)中的全部觀數(shù)的數(shù)據(jù),參考掃描獲得的參考數(shù)據(jù)作為K空間滑動完整數(shù)據(jù)集的起始點,且此次掃描記為記作第零次掃描。
步驟12、對目標視野施加更新掃描,并結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù),且每次采集的數(shù)據(jù)均包含K空間中高頻成分與低頻成分。
本發(fā)明實施例中,參數(shù)對包括:更新比與偏移量;
其中,所述更新比,指更新掃描采集的部分觀數(shù)占參考掃描采集的全部觀數(shù)的百分比;
所述偏移量,是用來表示更新掃描采集的部分觀數(shù)相對K空間中心位置的相對位置,定義為更新掃描采集的部分觀數(shù)中被K空間中心線分割的較小部分與較大部分的比值。
所述更新比為25%、33%、40%、45%或者50%;
所述偏移量為1:4、1:3、1:2或者1:1。
本發(fā)明實施例中,根據(jù)更新掃描是奇數(shù)次還是偶數(shù)次來采集K空間中不同部分的數(shù)據(jù);具體來說:
如果是奇數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間上半空間為主的部分觀數(shù)數(shù)據(jù);示例性的,如圖3所示,此時,更新比為40%,偏移量為1:3。
如果是偶數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間下半空間為主的部分觀數(shù)數(shù)據(jù)。示例性的,如圖4所示,此時,更新比為40%,偏移量為1:3。
如圖3-圖4可以看到,更新的部分觀數(shù)是相鄰的,即連在一起的。而且奇數(shù)次更新掃描采集的部分數(shù)據(jù)與偶數(shù)次更新掃描采集的部分數(shù)據(jù)是關(guān)于K空間中心對稱的。部分觀數(shù)的具體數(shù)量以及具體相位編碼位置由參數(shù)對(更新比與偏移量)決定。一般而言,部分觀數(shù)既包含K空間中的高頻成分又包含低頻成分。
本領(lǐng)域技術(shù)人員可以理解,圖2-4中矩形框的中心點為基準點,距離該基準點較遠的數(shù)據(jù)則為高頻成分,較近則為低頻成分。
步驟13、每次更新掃描之后,將K空間滑動完整數(shù)據(jù)集中與更新掃描采集的觀數(shù)數(shù)據(jù)相對應的相位編碼行進行替換更新。
步驟14、對更新后的K空間滑動完整數(shù)據(jù)集進行二維傅里葉逆變換,并對變換結(jié)果的值作取絕對值和歸一化處理得到圖像矩陣并呈現(xiàn)相應圖像。
此步驟中,進行二維傅里葉逆變換時,施加兩次90°相位偏移。
本發(fā)明實施例的上述方案,每次進行更新掃描時,所采集到的K空間部分觀數(shù)數(shù)據(jù)同時包含高頻和低頻成分,因而可以既保持較高的空間分辨率又可以提高時間分辨率,并加快掃描速度。
通過以上的實施方式的描述,本領(lǐng)域的技術(shù)人員可以清楚地了解到上述實施例可以通過軟件實現(xiàn),也可以借助軟件加必要的通用硬件平臺的方式來實現(xiàn)?;谶@樣的理解,上述實施例的技術(shù)方案可以以軟件產(chǎn)品的形式體現(xiàn)出來,該軟件產(chǎn)品可以存儲在一個非易失性存儲介質(zhì)(可以是CD-ROM,U盤,移動硬盤等)中,包括若干指令用以使得一臺計算機設備(可以是個人計算機,服務器,或者網(wǎng)絡設備等)執(zhí)行本發(fā)明各個實施例所述的方法。
本發(fā)明的另一實施例還提供一種磁共振快速成像系統(tǒng),該系統(tǒng)可以用來實現(xiàn)前述實施例的成像方法。如圖5所示,該系統(tǒng)主要包括:
參考掃描模塊51,用于對目標視野施加參考掃描,采集K空間原始數(shù)據(jù)中的全部觀數(shù)的數(shù)據(jù),作為K空間滑動完整數(shù)據(jù)集的起始點;
更新掃描模塊52,用于對目標視野施加更新掃描,并結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù);
替換模塊53,用于每次更新掃描之后,將K空間滑動完整數(shù)據(jù)集中與更新掃描采集的觀數(shù)數(shù)據(jù)相對應的相位編碼行進行替換更新;
成像模塊54,用于對更新后的K空間滑動完整數(shù)據(jù)集進行二維傅里葉逆變換,并對變換結(jié)果的值作取絕對值和歸一化處理得到圖像矩陣并呈現(xiàn)相應圖像。
本發(fā)明實施例中,所述參數(shù)對包括:更新比與偏移量;
其中,所述更新比,指更新掃描采集的部分觀數(shù)占參考掃描采集的全部觀數(shù)的百分比;
所述偏移量,是用來表示更新掃描采集的部分觀數(shù)相對K空間中心位置的相對位置,定義為更新掃描采集觀數(shù)中被K空間中心線分割的較小部分與較大部分的比值。
本發(fā)明實施例中,所述更新比為25%、33%、40%、45%或者50%;
所述偏移量為1:4、1:3、1:2或者1:1。
本發(fā)明實施例中,所述結(jié)合預設的參數(shù)對,來采集K空間上半空間為主的部分觀數(shù)數(shù)據(jù)或者下半空間為主的部分觀數(shù)數(shù)據(jù)包括:
如果是奇數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間上半空間為主的部分觀數(shù)數(shù)據(jù);
如果是偶數(shù)次更新掃描,則結(jié)合預設的參數(shù)對,采集以K空間下半空間為主的部分觀數(shù)數(shù)據(jù)。
本發(fā)明實施例中,所述進行二維傅里葉逆變換時,施加兩次90°相位偏移。
需要說明的是,上述系統(tǒng)中包含的各個功能模塊所實現(xiàn)的功能的具體實現(xiàn)方式在前面的各個實施例中已經(jīng)有詳細描述,故在這里不再贅述。
所屬領(lǐng)域的技術(shù)人員可以清楚地了解到,為描述的方便和簡潔,僅以上述各功能模塊的劃分進行舉例說明,實際應用中,可以根據(jù)需要而將上述功能分配由不同的功能模塊完成,即將系統(tǒng)的內(nèi)部結(jié)構(gòu)劃分成不同的功能模塊,以完成以上描述的全部或者部分功能。
以上所述,僅為本發(fā)明較佳的具體實施方式,但本發(fā)明的保護范圍并不局限于此,任何熟悉本技術(shù)領(lǐng)域的技術(shù)人員在本發(fā)明披露的技術(shù)范圍內(nèi),可輕易想到的變化或替換,都應涵蓋在本發(fā)明的保護范圍之內(nèi)。因此,本發(fā)明的保護范圍應該以權(quán)利要求書的保護范圍為準。