CN110133558B - 一种磁共振动态成像方法、装置及可读介质 - Google Patents
一种磁共振动态成像方法、装置及可读介质 Download PDFInfo
- Publication number
- CN110133558B CN110133558B CN201810132873.6A CN201810132873A CN110133558B CN 110133558 B CN110133558 B CN 110133558B CN 201810132873 A CN201810132873 A CN 201810132873A CN 110133558 B CN110133558 B CN 110133558B
- Authority
- CN
- China
- Prior art keywords
- data
- space
- projection
- reconstruction
- cardiac cycle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 44
- 238000005070 sampling Methods 0.000 claims abstract description 100
- 230000000747 cardiac effect Effects 0.000 claims abstract description 87
- 238000000034 method Methods 0.000 claims abstract description 46
- 230000000241 respiratory effect Effects 0.000 claims abstract description 33
- 230000003190 augmentative effect Effects 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000013138 pruning Methods 0.000 claims 1
- 238000002595 magnetic resonance imaging Methods 0.000 abstract description 11
- 125000004122 cyclic group Chemical group 0.000 abstract 1
- 230000001133 acceleration Effects 0.000 description 5
- 230000002123 temporal effect Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 244000248162 Xanthoceras sorbifolium Species 0.000 description 1
- 238000011166 aliquoting Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- RGCLLPNLLBQHPF-HJWRWDBZSA-N phosphamidon Chemical compound CCN(CC)C(=O)C(\Cl)=C(/C)OP(=O)(OC)OC RGCLLPNLLBQHPF-HJWRWDBZSA-N 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及磁共振成像的技术领域,公开了一种磁共振动态成像方法、装置及可读介质,包括以下步骤:步骤1:接收到触发信号后,在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据,其中,所述黄金角径向采样轨迹采用投影数量满足斐波那契数的均匀径向采样轨迹来近似;步骤2:在径向采样的投影和心动周期这两个维度上采用a‑f BLAST方法进行图像重建。本发明采用特殊的均匀径向采样轨迹(投影数是斐波那契数)近似黄金角的径向采样轨迹,打破了采样模式是循环偏移的均匀径向采样的限制,能够应用于动态成像常用的黄金角径向采样轨迹。
Description
技术领域
本发明涉及磁共振成像的技术领域,特别涉及一种磁共振动态成像方法、装置及可读介质。
背景技术
磁共振成像是利用核磁共振原理,通过静磁场、射频场和梯度线圈等对人体组织成像。它具有丰富的组织对比度和较高的软组织分辨率,且对人体无害,成为了医学临床诊断的一种强有力的工具。但是,磁共振成像通常需要较长的扫描时间,尤其在心脏实时成像等动态领域不能满足临床应用的要求。成像速度慢是制约其快速发展和广泛应用的一大瓶颈,所以如何提高成像速度一直是磁共振学术界和工业界的关注热点。
目前的动态成像的加速技术根据原理可以分为三类:
第一类是仅仅利用空间冗余性,如Riederer SJ等人(Riederer SJ,Tasciyan T,Farzaneh F,Lee JN,Wright RC,Herfkens RJ.MR fluoroscopy:technicalfeasibility.Magn Reson Med 1988;8:1–15.)研究的MR fluoroscopy技术和van VaalsJJ等人(van Vaals JJ,Brummer ME,Dixon WT,et al.“Keyhole”method foraccelerating imaging of contrast agent uptake.J Magn Reson Imaging 1993;3:671–675.)研究的Keyhold技术。
第二类是先利用时间冗余性,再利用空间冗余性,如Kellman P等人(Kellman P,Epstein FH,McVeigh ER.Adaptive sensitivity encoding incorporating temporalfiltering(TSENSE)Magn Reson Med 2001;45:846–852.)研究的利用时间滤波器的自适应敏感度编码(TSENSE)技术,Breuer FA等人(Breuer FA,Kellman P,Griswold MA,JakobPM.Dynamic autocalibrated parallel imaging using temporal GRAPPA(TGRAPPA).Magn Reson Med 2005;53:981–985.)研究的动态广义自动校准部分并行采集(TGRAPPA)技术等。
第三类是同时利用时间和空间冗余性,如Madore B等人(Madore B,Glover GH,Pelc NJ.Unaliasing by Fourier-encoding the overlaps using the temporaldimension(UNFOLD),applied to cardiac imaging and fMRI.Magn Reson Med 1999;42:813–828.)研究的UNFOLD技术,Tsao J等人(Tsao J,Boesiger P,Pruessmann KP.k-tBLAST and k-t SENSE:dynamic MRI with high frame rate exploitingspatiotemporal correlations.Magn Reson Med 2003;50:1031–1042.)研究的k-t BLAST技术和Hansen M S等人(Hansen M S,Baltes C,Tsao J,et al.k‐t BLASTreconstruction from non‐Cartesian k‐t space sampling[J].Magnetic resonance inmedicine,2006,55(1):85-91.)研究的k-tSENSE技术,Jung H等人(Jung H,Sung K,NayakKS,Kim EY,Ye JC.k-t FOCUSS:a general compressed sensing framework for highresolution dynamic MRI.Magn Reson Med 2009;61:103–116.)研究的k-t FOCUSS技术等。
这些动态成像的加速技术根据采样轨迹又可以分为笛卡尔采样和非笛卡尔采样(如径向采样)。非笛卡尔采样可以有效利用梯度场,从而快速填充K空间。而且非笛卡尔采样对K空间中心过采样,重建图像对比度更高,对运动伪影更不敏感。但是,非笛卡尔采样采集到的数据,由于数据没有落在K空间规则分布的笛卡尔网格点上,其重建不能直接应用快速傅里叶变换。现有的重建方法多数是基于迭代的重建方法,如Hansen M S等人(Hansen MS,Baltes C,Tsao J,et al.k‐t BLAST reconstruction from non‐Cartesian k‐t spacesampling[J].Magnetic resonance in medicine,2006,55(1):85-91.)研究的非笛卡尔采样的k-t BLAST重建方法,重建速度较慢,不能满足实时成像的要求。
最近,Madison Kretzler等人(Kretzler,M.,Hamilton,J.,Griswold,M.&Seiberlich,N.a-f BLAST:A Non-Iterative Radial k-t BLAST Reconstruction inRadon Space.in ISMRM(2016).)在2016年ISMRM会议上提出了一种非迭代的径向采样动态成像加速方法(a-f BLAST)。A-f BLAST方法主要思想是将这些数据变换到a-f空间,然后在a-f空间使用BLAST算法进行重建。A-fBLAST是非迭代的重建方法,相比以往常用的迭代重建,如Hansen M S等人(Hansen M S,Baltes C,Tsao J,et al.k‐t BLAST reconstructionfrom non‐Cartesian k‐t space sampling[J].Magnetic resonance in medicine,2006,55(1):85-91.)研究的非笛卡尔采样的k-t BLAST重建方法,Pruessmann KP等人(Pruessmann KP,Weiger M,Bornert P,Boesiger P.Advances in sensitivity encodingwith arbitrary k-space trajectories.Magn Reson Med 2001;46:638–651.)研究的CGSENSE技术等,A-f BLAST非迭代的重建方法的重建速度大大加快。
A-f BLAST方法结合了径向采样快速填充K空间和非迭代快速重建的优点。但是,a-f BLAST方法要求的采样模式是循环偏移的均匀径向采样(uniform sampling in aninterleaved fashion),而动态成像的径向采样模式通常是黄金角(golden angle)采样。后者可以利用任意数量的投影进行一帧图像的重建,从而在动态成像中的任意时间点进行任意时间分辨率的重建。A-f BLAST的均匀径向采样模式限制了黄金角(golden angle)采样的使用,从而丧失了在任意时间点进行任意时间分辨率重建的良好性质。
发明内容
本发明要解决的技术问题是A-f BLAST的均匀径向采样模式限制了黄金角(golden angle)采样的使用,从而丧失了在任意时间点进行任意时间分辨率重建的良好性质。
为了解决上述技术问题,本发明中针对这一问题,通过研究a-f BLAST和黄金角(golden angle)采样的特点,将a-f BLAST方法拓展应用到黄金角(golden angle)采样的磁共振动态成像中。本发明在使用a-f BLAST方法加速磁共振成像的同时,保持黄金角(golden angle)采样在任意时间点进行任意时间分辨率重建的良好性质。
本发明的技术方案是这样实施的:提供一种磁共振动态成像方法,包括以下步骤:
步骤1:接收到触发信号后,在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据,其中,所述黄金角径向采样轨迹采用投影数量满足斐波那契数的均匀径向采样轨迹来近似;
步骤2:在径向采样的投影和心动周期这两个维度上采用a-f BLAST方法进行图像重建。
优选地,在所述步骤1之前还包括步骤:在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。
优选地,所述步骤1中,在不同的呼吸周期之间,分段采集nRS个呼吸周期,一个呼吸周期内至少包含4个心动周期,每个心动周期内利用nSG个投影用于一帧图像的重建,其中,TR×nSG×nPS<Tc,TR为序列的重复时间,nPS为每个心动周期成像的心脏时相,Tc为心动周期的平均时长。
优选地,所述步骤1中,在每个心动周期之内采用黄金角径向采样轨迹采集整个心动周期内的K空间数据,所述黄金角径向采样的投影角度在每个呼吸周期之内的4个心动周期之间有θf的角度偏移,即第k个心动周期黄金角径向采样的起始角度是
θf×(k-1),其中,k=0,1,2,3,以及
优选地,所述步骤1中,采集的K空间数据为其中,所述的大小为nR×(nSG·nRS)×nPS×4,所述kr表示每个投影的读出方向,kr∈{1,2,...,nR},所述kp表示用于重建一帧图像的投影,kp∈{1,2,...,nSG×nRS},所述tps表示帧,tPS∈{1,2,...,nPS},所述tc表示一个呼吸周期内的心动周期,tC∈{1,2,3,4}。
优选地,所述步骤2中,采用a-f BLAST方法进行图像重建包括以下步骤:
步骤2-1:寻找一个小于nSG×nPS的最大斐波那契数nFB,即nFB={1,1,2,3,5,8,13,21,34,...}∩{n|n<nSG×nPS},将所述数据第二维的最后(nSG×nPS-nFB)个投影丢弃,将数据修剪为其中,所述的数据大小为nR×nFB×nPS×4;
步骤2-4:将所述数据的投影(kp)和心动周期(tc)这两个维度的数据扩充为数据得到扩充后的数据其中,所述数据的大小为nFB×4,所述数据的大小为(4·nFB)×4,所述数据的大小为nR×(4·nFB)×nPS×4;
优选地,所述步骤2-5中,所述a-f BLAST数据重建包括如下步骤:
优选地,所述步骤S3包括:
S3-5:将作二维傅里叶变换,得到a-f BLAST重建后的二维数据对数据的不同x和tps,重复投影(kp)和心动周期(tc)这两个维度的数据进行a-f BLAST重建,得到重建后的数据所述数据的大小为nR×(4·nFB)×nPS×4。
本发明提供一种磁共振动态成像装置,包括:
数据采集模块,用于在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据;
a-f BLAST图像重建模块,用于在径向采样的投影和心动周期这两个维度上采用a-f BLAST方法进行图像重建。
优选地,所述数据采集模块还包括:触发模块,用于在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。
优选地,所述a-f BLAST图像重建模块包括:
数据扩展模块,用于对所述数据的投影(kp)和心动周期(tc)这两个维度的数据扩充为数据得到扩充后的数据其中,所述数据的大小为nFB×4,所述数据的大小为(4·nFB)×4,所述数据的大小为nR×(4·nFB)×nPS×4;
优选地,所述a-f BLAST数据重建模块包括:
本发明还提供一种计算机可读介质,该计算机可读介质具有存储在其中的程序,该程序是计算机可执行的以使计算机执行基于磁共振动态成像方法的各步骤。
实施本发明的有益效果主要有:
1、本发明通过研究a-f BLAST方法和黄金角(golden angle)采样的特点,将a-fBLAST方法拓展应用到黄金角(golden angle)采样的磁共振动态成像中。相比于现有技术a-f BLAST方法,采用特殊的均匀径向采样轨迹(投影数是斐波那契数)近似黄金角(goldenangle)的径向采样轨迹,打破了采样模式是循环偏移的均匀径向采样(uniform samplingin an interleaved fashion)的限制,能够应用于动态成像常用的黄金角(golden angle)径向采样轨迹。
2、本发明在使用a-f BLAST方法加速磁共振成像的同时,能保持黄金角(goldenangle)采样在任意时间点进行任意时间分辨率重建的良好性质。
3、本发明对于动态成像常用的黄金角(golden angle)径向采样轨迹的数据进行a-f BLAST重建,利用一个呼吸周期内的心动周期之间的时间冗余性进行加速成像,可以做到4倍加速的成像。
附图说明
为更好地理解本发明的技术方案,可参考下列的、用于对现有技术或实施例进行说明的附图。这些附图将对部分实施例或现有技术涉及的产品或方法进行简要的展示。这些附图的基本信息如下:
图1为一实施例中,一种磁共振动态成像方法流程图;
图2为一实施例中,一种磁共振动态成像方法数据采集方案示意图;
图3为一实施例中,一种采用a-f BLAST方法进行图像重建流程图;
图4为一测试例中,一个心动周期内的3帧NUFFT重建结果图;
图5为一测试例中,一个心动周期内的3帧a-f BLAST重建结果图。
具体实施方式
现在对本发明实施例中的技术方案或有益效果作进一步的展开描述,显然,所描述的实施例仅是本发明的部分实施方式,而并非全部。
本发明考虑到一个呼吸周期内至少包含4个心动周期,而且由于呼吸等人体运动导致这4个心动周期在同一时相的心脏形态也不一致,所以成像方案是将一个呼吸周期内的4个心动周期也重建出来。在每个心动周期之内采用黄金角(golden angle)的径向采样轨迹,径向采样的投影角度在每个呼吸周期之内的4个心动周期之间有偏移,而在不同的呼吸周期之间采用分段(segment)采集方式填充K空间。针对黄金角(golden angle)采样轨迹,采用特殊的均匀径向采样轨迹来近似,即均匀径向采样的投影(projection)数量满足斐波那契数。最后在径向采样的投影(projection)和心动周期(cardiac cycle)这两个维度上采用a-f BLAST方法进行重建,从而达到4倍加速的动态磁共振成像。
实施例一
如图1和图2所示,本实施例提供一种加速采集方案,采用黄金角(golden angle)径向采样轨迹、分段(segment)采集的采样方式采集数据。
如图1所示,本实施例提供的一种磁共振动态成像方法,包括以下步骤:
步骤1:接收到触发信号后,在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据,其中,所述黄金角径向采样轨迹采用投影数量满足斐波那契数的均匀径向采样轨迹来近似;步骤2:在径向采样的投影和心动周期这两个维度上采用a-f BLAST方法进行图像重建,利用一个呼吸周期内的心动周期之间的时间冗余性进行加速成像。
如图2所示,本实施例分段(segment)采集nRS个呼吸周期,每个呼吸周期之内采集4个心动周期,每个心动周期采用黄金角(golden angle)径向采样轨迹,相邻两个心动周期之间有θf的角度偏移。
具体的,所述步骤1中,在不同的呼吸周期之间,假设序列的重复时间是TR,分段采集nRS个呼吸周期,一个呼吸周期内至少包含4个心动周期,每个心动周期成像的心脏时相(phase)是nPS,即每个心动周期的成像有nPS帧。每个心动周期内利用nSG个投影用于一帧图像的重建,那么分段采集nRS个呼吸周期之后有nRS×nSG个投影用于一帧图像的重建,其中,TR×nSG×nPS<Tc,TR为序列的重复时间,nPS为每个心动周期成像的心脏时相,Tc为心动周期的平均时长。
具体的,所述步骤1中,首先检测到呼吸信号,然后在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。接收到触发信号后,即心电信号触发后,在每个心动周期之内采用黄金角径向采样轨迹采集整个心动周期内的K空间数据,所述黄金角径向采样的投影角度在每个呼吸周期之内的4个心动周期之间有θf的角度偏移,即第k个心动周期黄金角径向采样的起始角度是θf×(k-1),其中,k=0,1,2,3,
以及
根据本实施例中的采样模式,采集和储存数据。具体的,步骤1中采集的K空间数据为其中,所述的大小为nR×(nSG·nRS)×nPS×4,下标kr表示每个投影的读出方向(readout),kr∈{1,2,...,nR},下标kp表示用于重建一帧图像的投影(projection),kp∈{1,2,...,nSG×nRS},下标tps表示帧,tPS∈{1,2,...,nPS},下标tc表示一个呼吸周期内的心动周期,tC∈{1,2,3,4}。
具体的,如图3所示,所述步骤2中,采用a-f BLAST方法进行图像重建包括以下步骤:
步骤2-1:修建数据,首先寻找一个小于nSG×nPS的最大斐波那契数nFB,即nFB={1,1,2,3,5,8,13,21,34,...}∩{n|n<nSG<nPs,然后将采集的数据第二维的最后(nSG×nPS-nFB)个投影丢弃,也就是将数据修剪为其中,所述的大小为nR×(4·nFB)×nPS×4。
步骤2-2:重排列,将数据的第二维,即采用黄金角(golden angle)径向采样轨迹采集的投影重排列成均匀径向采样轨迹(nFB等分),得到数据其中,所述nFB是斐波那契数,所述的大小为nR×(4·nFB)×nPS×4。
步骤2-4:对所述数据的投影(kp)和心动周期(tc)这两个维度的数据扩充为所述数据的大小为nFB×4,所述数据的大小为(4·nFB)×4,其中,(4×(n-1)+k,k)位置的值与在(n,k)位置的值相同,其余位置的值为零,n=1,2,...,nFB,k=1,2,3,4。得到扩充后的数据所述数据的大小为nR×(4·nFB)×nPS×4。
步骤2-5:对数据的投影(kp)和心动周期(tc)这两个维度的数据进行a-f BLAST数据重建,得到重建后的数据所述数据的大小为nR×(4·nFB)×nPS×4。具体的,对数据的不同x和tps,重复投影(kp)和心动周期(tc)这两个维度的数据进行a-f BLAST重建,得到重建后的数据
其中,所述步骤2-5中的所述a-f BLAST数据重建包括如下步骤:
具体的,所述步骤S3包括:
S3-5:将作二维傅里叶变换,得到a-f BLAST重建后的二维数据对数据的不同x和tPS,重复投影(kp)和心动周期(tc)这两个维度的数据进行a-fBLAST重建,得到重建后的数据所述数据为nR×(4·nFB)×nPS×4。
本实施例中,将混叠的a-f空间进行a-f BLAST算法重建后,在心动周期tc做傅里叶逆变换,在读出方向(kr)和投影(kp)方向做NUFFT,得到a-f BLAST算法重建后的动态磁共振图像。
实施例二
本实施例提供一种磁共振动态成像装置,包括:
数据采集模块,用于在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据;a-f BLAST图像重建模块,用于在径向采样的投影和心动周期这两个维度上采用a-f BLAST方法进行图像重建。
具体的,所述a-f BLAST图像重建模块,包括:
数据扩展模块,用于对所述数据的投影(kp)和心动周期(tc)这两个维度的数据扩充为得到扩充后的数据其中,位置的值与在(n,k)位置的值相同,其余位置的值为零,n=12,...,nFB,k=1,2,3,4。其中,所述的大小为nFB×4,所述的大小为(4·nFB)×4,所述的大小为nR×(4·nFB)×nPS×4。
其中,所述a-f BLAST数据重建模块包括:混叠a-f空间模块,用于对作二维傅里叶逆变换,得到二维的混叠a-f空间低分辨率a-f空间模块,用于对数据在投影(kp)维度作插值和滤波操作,然后作傅里叶逆变换,得到低分辨率的a-f空间混叠a-f空间数据重建模块,用于将混叠的a-f空间进行a-f BLAST算法重建,得到数据其中,所述数据的大小为nR×(4·nFB)×nPS×4。
作为一优选实施例,所述数据采集模块还包括:触发模块,用于在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。
本实施例还提供一种计算机可读介质,该计算机可读介质具有存储在其中的程序,该程序是计算机可执行的以使计算机执行基于磁共振动态成像方法的各步骤。
测试例
本测试例分别使用NUFFT方法和a-f BLAST方法进行一个心动周期的动态磁共振图像重建,重建结果如图4和图5所示。
图4是直接使用NUFFT重建的结果,展示了一个心动周期内的3帧图像。每一帧图像使用了径向采样K空间的8个投影(kp)进行NUFFT重建,可以看到重建图像有明显的放射状伪影,这是使用比较少的投影进行重建时出现的典型伪影。
图5是使用a-f BLAST重建的结果,展示了一个心动周期内的3帧图像。每一帧图像也是使用径向采样K空间的8个投影(kp)进行重建,由于使用了a-f BLAST重建,加速因子是4,最后一个步骤对径向采样K空间做NUFFT时,相当于有8×4=32个投影(kp)做NUFFT,所以重建图像的质量大大提高,没有肉眼可见的放射状伪影。
最后需要指出的是,上文所列举的实施例,为本发明较为典型的、较佳实施例,仅用于详细说明、解释本发明的技术方案,以便于读者理解,并不用以限制本发明的保护范围或者应用。因此,在本发明的精神和原则之内所作的任何修改、等同替换、改进等而获得的技术方案,都应被涵盖在本发明的保护范围之内。
Claims (13)
1.一种磁共振动态成像方法,其特征在于,包括以下步骤:
步骤1:接收到触发信号后,在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据,其中,所述黄金角径向采样轨迹采用投影数量满足斐波那契数的均匀径向采样轨迹来近似;
步骤2:在径向采样的投影和心动周期这两个维度上采用a-f BLAST方法进行图像重建。
2.如权利要求1所述的磁共振动态成像方法,其特征在于,在所述步骤1之前还包括步骤:在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。
3.如权利要求1所述的磁共振动态成像方法,其特征在于:所述步骤1中,在不同的呼吸周期之间,分段采集nRS个呼吸周期,一个呼吸周期内至少包含4个心动周期,每个心动周期内利用nSG个投影用于一帧图像的重建,其中,TR×nSG×nPS<Tc,TR为序列的重复时间,nPS为每个心动周期成像的心脏时相,Tc为心动周期的平均时长。
6.如权利要求5所述的磁共振动态成像方法,其特征在于,所述步骤2中,采用a-fBLAST方法进行图像重建包括以下步骤:
步骤2-1:寻找一个小于nSG×nPS的最大斐波那契数nFB,即nFB={1,1,2,3,5,8,13,21,34,…}∩{n|n<nSG×nPS},将所述数据第二维的最后(nSG×nPS-nFB)个投影丢弃,将数据修剪为数据其中,所述数据的大小为nR×nFB×nPS×4;
步骤2-4:将所述数据的投影kp和心动周期tc这两个维度的数据填零扩充为得到填零扩充后的数据其中,所述数据的大小为nFB×4,所述数据的大小为(4·nFB)×4,所述数据的大小为nR×(4·nFB)×nPS×4;
8.如权利要求7所述的磁共振动态成像方法,其特征在于,所述步骤S3包括:
9.一种磁共振动态成像装置,其特征在于,包括:
数据采集模块,用于在不同的呼吸周期之间采用分段采集的采样方式填充K空间,在每个心动周期之内采用黄金角径向采样轨迹的采样方式采集K空间数据;
a-f BLAST图像重建模块,用于在径向采样的投影和心动周期这两个维度上采用a-fBLAST方法进行图像重建。
10.如权利要求9所述的磁共振动态成像装置,其特征在于,所述数据采集模块还包括:触发模块,用于在呼吸周期内检测心电信号的R波,接收到所述R波后产生触发信号。
11.如权利要求9所述的磁共振动态成像装置,其特征在于,所述a-f BLAST图像重建模块包括:
修建数据模块,用于将数据修剪为其中,所述为采集的K空间数据;kr表示每个投影的读出方向,kr∈{1,2,...,nR},kp表示用于重建一帧图像的投影,kp∈{1,2,...,nSG×nRS},tPS表示帧,tPS∈{1,2,...,nPS},tc表示一个呼吸周期内的心动周期,tC∈{1,2,3,4},所述数据的大小为nR×nFB×nPS×4;
数据扩展模块,用于对所述数据的投影kp和心动周期tc这两个维度的数据扩充为数据得到扩充后的数据其中,所述的大小为nFB×4,所述的大小为(4·nFB)×4,所述数据的大小为nR×(4·nFB)×nPS×4;
其中,在不同的呼吸周期之间,分段采集nRS个呼吸周期,一个呼吸周期内至少包含4个心动周期,每个心动周期内利用nSG个投影用于一帧图像的重建,其中,TR×nSG×nPS<Tc,TR为序列的重复时间,nPS为每个心动周期成像的心脏时相,Tc为心动周期的平均时长。
13.一种计算机可读介质,该计算机可读介质具有存储在其中的程序,该程序是计算机可执行的以使计算机执行权利要求1~8中任一项所述的基于磁共振动态成像方法的各步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810132873.6A CN110133558B (zh) | 2018-02-09 | 2018-02-09 | 一种磁共振动态成像方法、装置及可读介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810132873.6A CN110133558B (zh) | 2018-02-09 | 2018-02-09 | 一种磁共振动态成像方法、装置及可读介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133558A CN110133558A (zh) | 2019-08-16 |
CN110133558B true CN110133558B (zh) | 2021-04-30 |
Family
ID=67567589
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810132873.6A Active CN110133558B (zh) | 2018-02-09 | 2018-02-09 | 一种磁共振动态成像方法、装置及可读介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133558B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11741356B2 (en) * | 2019-02-08 | 2023-08-29 | Korea Advanced Institute Of Science & Technology | Data processing apparatus by learning of neural network, data processing method by learning of neural network, and recording medium recording the data processing method |
CN112881958B (zh) * | 2021-02-04 | 2022-02-25 | 上海交通大学 | 基于低秩与稀疏分解的磁共振介入成像方法、系统及介质 |
CN112967297B (zh) * | 2021-03-23 | 2024-05-03 | 中国科学院深圳先进技术研究院 | 磁共振氧十七代谢成像方法、装置、存储介质及终端设备 |
CN113866695B (zh) * | 2021-10-12 | 2022-08-26 | 上海交通大学 | 磁共振实时引导介入的图像采集与重建方法及系统 |
CN118078249A (zh) * | 2024-04-23 | 2024-05-28 | 中国医学科学院阜外医院 | 心脏磁共振成像水脂分离的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102008307A (zh) * | 2010-12-29 | 2011-04-13 | 中国科学院深圳先进技术研究院 | 磁共振弥散张量成像方法和系统 |
CN103845055A (zh) * | 2012-12-04 | 2014-06-11 | 上海联影医疗科技有限公司 | 心脏磁共振成像方法及系统 |
US20150309135A1 (en) * | 2014-04-25 | 2015-10-29 | New York University | System, method and computer-accessible medium for rapid real-time cardiac magnetic resonance imaging utilizing synchronized cardio-respiratory sparsity |
CN107110940A (zh) * | 2014-10-07 | 2017-08-29 | 西诺德牙科设备有限公司 | 用于mri径向或螺旋成像的方法 |
-
2018
- 2018-02-09 CN CN201810132873.6A patent/CN110133558B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102008307A (zh) * | 2010-12-29 | 2011-04-13 | 中国科学院深圳先进技术研究院 | 磁共振弥散张量成像方法和系统 |
CN103845055A (zh) * | 2012-12-04 | 2014-06-11 | 上海联影医疗科技有限公司 | 心脏磁共振成像方法及系统 |
US20150309135A1 (en) * | 2014-04-25 | 2015-10-29 | New York University | System, method and computer-accessible medium for rapid real-time cardiac magnetic resonance imaging utilizing synchronized cardio-respiratory sparsity |
CN107110940A (zh) * | 2014-10-07 | 2017-08-29 | 西诺德牙科设备有限公司 | 用于mri径向或螺旋成像的方法 |
Non-Patent Citations (1)
Title |
---|
径向K空间采样技术在胎儿心脏磁共振成像中的应用;张弘等;《中国医学物理学杂志》;20160630;第33卷(第6期);580-583 * |
Also Published As
Publication number | Publication date |
---|---|
CN110133558A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110133558B (zh) | 一种磁共振动态成像方法、装置及可读介质 | |
US10605878B2 (en) | System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial sampling and compressed sensing | |
CN107750338B (zh) | 使用星形叠层采集的mr成像方法和设备 | |
JP6513398B2 (ja) | 事前情報に制約される正則化を使用するmr画像再構成 | |
US6448771B1 (en) | Magnetic resonance method for forming a fast dynamic image | |
KR101939642B1 (ko) | 모션 분해 mri를 위한 시스템 및 방법 | |
EP3322997A1 (en) | Mr imaging with motion detection | |
JP2001161657A (ja) | 核磁気共鳴撮影装置 | |
JP2020523075A (ja) | 可変コントラストのスタック・オブ・スター収集を使用したmrイメージング | |
JP3699304B2 (ja) | 磁気共鳴撮像装置 | |
EP2726893B1 (en) | A method and system for rapid mri acquisition using tailored signal extraction modules | |
JP7446328B2 (ja) | ソフト動きゲーティングによる3dラジアル又はスパイラル収集を用いたmr画像 | |
JP3051374B2 (ja) | 磁気共鳴撮像装置 | |
US11016154B2 (en) | Magnetic resonance imaging apparatus and image reconstruction program | |
EP1444529B1 (en) | Magnetic resonance method for forming a fast dynamic image | |
US11269037B2 (en) | MR imaging using motion-dependent radial or spiral k-space sampling | |
US20230341493A1 (en) | Magnetic resonance imaging with prior knowledge and oversampling | |
EP1642149A1 (en) | Undersampled magnetic resonance imaging | |
CN112327233B (zh) | 多相位快速磁共振弹性成像采集与重建方法及系统 | |
Ji | Compressive sensing imaging with randomized lattice sampling: Applications to fast 3D MRI | |
WO2024081777A2 (en) | Methods for accelerated time-resolved multi-echo magnetic resonance imaging by augmenting temporal correlation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |