平面回波成像序列图像的重建方法
【专利摘要】本发明提供了一种平面回波成像序列的图像重建方法，包括获取平面回波成像数据Si，并同时采集三条没有经过相位编码的参考回波信号R1、R2、R3；通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FSi，并用所述校正参数校正FSi，计算出校正后的平面回波成像数据；对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。本发明提供的平面回波成像序列的图像重建方法，在有偏离场存在的情况下，能够保持平面回波序列成像快速成像优点且有效地去除N/2伪影，同时还能够校正由于偏离场存在而导致的图像变形。
【专利说明】平面回波成像序列图像的重建方法
【【技术领域】】
[0001]本发明涉及磁共振成像领域，尤其是涉及一种平面回波成像序列图像的重建方法。
【【背景技术】】
[0002]目前在磁共振成像过程中，平面回波成像序列是一种快速的、只需要单个激发脉冲就可获得完整图像的磁共振成像序列。磁共振梯度系统周期性变换，产生一系列梯度回波；受激平面的图像通过对产生的回波序列使用傅里叶变换而得到。
[0003]在传统的平面回波成像中，序列的相位编码梯度和频率编码梯度分别如图1中表示的a和b，由于硬件系统的不完美，比如涡电流，数据采集延迟，磁场漂移等原因，在正和负的频率编码梯度下读出的数据存在差异，如回波中心不一致，k空间每条线的间隔不等等情况。由于回波序列中奇数和偶数回波时间的不匹配或者相位差导致在重建之后得到的正常图像在偏离实际的图像视野一半的位置出现低强度的附加图像，此即为N/2伪影，图像出现了整体漂移以及变形如图5所示。
[0004]为了尽可能地消除N/2伪影，现有技术中采用了一些校正方法，主要为两种:相位校正和单边读出梯度法。
[0005]相位校正法包括一维相位校正方法以及二维相位校正法。一维相位校正法，在读取数据之前，采集几个不加相位编码的回波信号，如图1中的回波信号C，根据第一个和第二个回波信号，算出它们相位差异，把这些相位差作为校正量来校正采集到的图像数据，即校正所有偶数项信号和奇数项信号之间的相位差。由于不是所有偶数项信号与奇数项信号之间的相位差都等于求得的校正量，所以该方法虽然简单，但是效果不好，并不能有效地消除N/2伪影，同时目前的一维相位校正法并不能同时校正平面回波图像的变形，这在低磁场的磁共振系统尤为明显。
[0006]二维相位校正法，同样是在读取图像数据前，采集一序列加相位编码的回波信号，利用这些采集到的数据纠正真正的图像数据。该方法消除N/2伪影的效果比较好，但是序列采集时间延长，平面回波成像序列失去了快速成像的优点。
[0007]单边读出梯度法只在一个极性的读出梯度时采集数据，该方法可以完全避免N/2伪影问题，但是不能纠正图像变形，浪费扫描时间。而且由于回波时间间隔比较长，图像在相位编码方向上的位移，以及因磁场不均匀导致的图像变形更为严重。
[0008]综上所述，现有技术中能够有效消除N/2伪影的校正方法采集时间较长或者不能够校正由于偏离场而导致的图像变形，采集时间较短的校正方法去除N/2伪影的效果比较差且同时也不能够校正图像的变形。因此，需要提供一种更为精准的平面回波成像序列的图像重建方法。
【
【发明内容】
】
[0009]为了解决现有技术中平面回波成像序列中消除N/2伪影效果差，不能校正图像变形，以及采集时间较长的问题，本发明提供了一种更为精准的平面回波成像序列的图像重
建方法。
[0010]一种平面回波成像序列的图像重建方法，包括如下步骤:
[0011]获取平面回波成像数据Si，并同时采集三条没有经过相位编码的参考回波信号RpR2、R3，所述三条参考回波信号分别为偶信号、奇信号以及偶信号；
[0012]通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；
[0013]将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FSi,并用所述校正参数校正FSi,计算出校正后的平面回波成像数据；
[0014]对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。
[0015]可选地，所述校正参数包括零次项相位偏差因子、一次项相位偏差因子
coir
Ψ\α °
[0016]可选地，对所述平面回波成像数据的偶数项信号或者奇数项信号进行校正，得到校正后的数据CFSi (η),
[0017]
【权利要求】
1.一种平面回波成像序列的图像重建方法，其特征在于，包括如下步骤: 获取平面回波成像数据Si,并同时采集三条没有经过相位编码的参考回波信号Rp R2,R3，所述三条参考回波信号分别为偶信号、奇信号以及偶信号； 通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数； 将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FSi,并用所述校正参数校正FSi,计算出校正后的平面回波成像数据； 对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。
2.如权利要求1所述平面回波成像序列图像重建方法，其特征在于，所述校正参数包括零次项相位偏差因子
3.如权利要求2所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号或者奇数项信号进行校正，得到校正后的数据CFSi(Ii)，
4.如权利要求2所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFSi(Ii)，
5.如权利要求1或2所述的平面回波成像序列图像重建方法，其特征在于，所述校正参数包括零次项相位偏移因子
6.如权利要求5所述的平面回波序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFSi(Ii)，
7.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFSi(Ii)，
8.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，所述一次项相位偏差因子
9.如权利要求8所述的平面回波成像序列图像重建方法，其特征在于，所述一次项相位偏差因子
10.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，所述零次项相位偏差因子C、所述零次项相位偏移因子与第一零次项相位、第二零次项相位
11.如权利要求10所述的平面回波成像序列图像重建方法，其特征在于，所述零次项相位偏差因子
12.如权利要求8所述的平面回波成像序列图像重建方法，其特征在于， 对所述参考回波信号Rp R2> R3傅里叶变换后复数共轭相乘，用来计算需要对所述平面回波成像数据进行校正的参数。
13.如权利要求12所述的平面回波成像序列图像重建方法，其特征在于， 将所述参考回波信号Rp R2> R3傅里叶变换得变换结果FRp FR2与FR3,
FR1 = FTid {RJ
FR2 = FTid {RJ
FR3 = FTid {R3} 然后将所述FR1与FR2、FR2与FR3复数共轭相乘得运算结果P1 (n)、P2 (η)，
P1 (η) = FR1 (η) X con j [FR2 (η)]
P2 (η) = FR3 (η) X con j [FR2 (η)] 其中η的取值范围为[l，Nro]； 分别取P1 (η)、P2 (η)的幅值最大值AMAXp AMAX2 ； 将所述P1 (η)及P2 (η)邻近两点的复数共轭相乘得运算结果DP1 (n)、DP2 (η)，
DP1 (η) = P1 (η+1) X con j [P1 (η)]
DP2 (η) = P2 (η+1) X con j [P2 (η)] 将所述DP1 (η)对应所述P1 (η)幅值大于τ X AMAX1的所有点相加得运算结果SDP1 (η)，
14.如权利要求13所述的平面回波成像序列图像重建方法，其特征在于， 对所述P?及匕⑷的相位进行校正，以消除<,和的影响，得到校正结果CP1 (η)、CP2(η)，
15.如权要求13或14中任一项所述的平面回波成像序列图像重建方法，其特征在于，所述预设参数τ = 0.5。
16.如权利要求13所述的平面回波成像序列图像重建方法，其特征在于， 若使用多通道采集平面回波数据，则将所述FR1与FR2、FR2与FR3复数共轭相乘然后进行多通道合并，包括:
【文档编号】G01R33/565GK104035059SQ201310072198
【公开日】2014年9月10日 申请日期:2013年3月6日 优先权日:2013年3月6日
【发明者】张卫国, 吴环, 蒋先旺 申请人:上海联影医疗科技有限公司