CN1515912A - 磁共振图像刚性平移运动伪影的自动逆向修正方法 - Google Patents
磁共振图像刚性平移运动伪影的自动逆向修正方法 Download PDFInfo
- Publication number
- CN1515912A CN1515912A CNA03140278XA CN03140278A CN1515912A CN 1515912 A CN1515912 A CN 1515912A CN A03140278X A CNA03140278X A CN A03140278XA CN 03140278 A CN03140278 A CN 03140278A CN 1515912 A CN1515912 A CN 1515912A
- Authority
- CN
- China
- Prior art keywords
- space line
- line
- image
- actual motion
- matrix
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005481 NMR spectroscopy Methods 0.000 title abstract 2
- 230000033001 locomotion Effects 0.000 claims abstract description 120
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 239000013598 vector Substances 0.000 claims abstract description 45
- 238000012937 correction Methods 0.000 claims abstract description 36
- 238000013519 translation Methods 0.000 claims abstract description 23
- 239000000047 product Substances 0.000 claims description 17
- 239000012467 final product Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- LPLLVINFLBSFRP-UHFFFAOYSA-N 2-methylamino-1-phenylpropan-1-one Chemical compound CNC(C)C(=O)C1=CC=CC=C1 LPLLVINFLBSFRP-UHFFFAOYSA-N 0.000 description 2
- 241000132539 Cosmos Species 0.000 description 2
- 235000005956 Cosmos caudatus Nutrition 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 208000019155 Radiation injury Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000000386 athletic effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Abstract
本发明公开了一种磁共振图像刚性平移运动伪影的自动逆向修正方法，获取伪影图像的原始频域矢量矩阵数据，并确定K空间坐标；以Ky等于零为参照，从第一个K空间线开始，逐个修正每一个K空间线偏移的相位；对要修正的K空间线，确定是否需要进行逆向估计线间运动距离，若需要，采用逆向修正得出实际运动距离来修正该K空间线；若不需要，可直接得出其实际运动距离来修正该K空间线；最后将所有修正好的K空间线变换到图像域，得到修正好的图像。本发明使得因刚性平移而损害的数据经修正后能重建出不含运动伪影的高质量的图像。
Description
技术领域
本发明涉及一种磁共振(MR)图像运动伪影的修正方法，尤其是涉及一种磁共振(MR)图像刚性平移运动伪影的自动逆向修正方法。
背景技术
由于MR具有无放射线伤害、分辨率很高等优点，故在临床诊断中的应用越来越广泛。但是MR成像时间较长，病人的自主性运动(如呼吸、心跳等)和非自主性运动会在图像中产生运动伪影，严重损害图像质量，从而影响对疾病的诊断。对于运动伪影的消除，基本分为两大类，即前处理与后处理。前处理技术中，主要使用快速扫描、门控扫描、欠扫描、过扫描、半扫描或可校正的重建模式，这些模式基本上也是成功的。但可校正重建去伪影的模式需要对各类运动建立起数学模型，而此类特定数学模型及其他扫描重建模式的事先选择参数是非常困难，因为患者的运动方式是随机的，不可预测的。至于后处理去伪影法，目前对此问题的研究，大都是在刚性运动的范畴内。现有方法需要人工参与，没有实现自动修正。事实上，任何实际的复杂刚性运动都是由两种简单的运动：平移和旋转复合而成，因此对平移和旋转运动造成的运动伪影的研究和消除技术为进一步研究实际的复杂运动产生的运动伪影奠定了基础，具有重要的意义。
MR刚性平移运动伪影消除的有关理论和技术可参阅下面三篇描述论文。
1、德里.M，严.H，罗斯费尔德.D“一种二维平移运动伪影修正的改进方法”，IEEE Trans.Med.Ima.，vol.10，No.4，pp.548-553，Dec.1991
2、瑞查.宙洛斐.A等“一种MRI图像平面平移运动伪影修正的改进方法”，IEEE Trans.Med.Ima.，Vol.14，No.3，PP471-479，Sept.1995
3、瑞查.宙洛斐.A等“刚性运动引起的MR图像平面伪影的消除方法”，IEEE Tans.Med.Ima.，Vol.15，No.6，PP768-784，Dec.1996
国际上对刚性平移运动伪影的消除研究较少，上面三篇论文提出的对刚性运动伪影的消除方法需要在图像域手工标识出图像空间约束区域，不能完全自动修正。国内目前尚无这方面的研究。
发明内容
本发明的目的在于提供一种磁共振图像刚性平移运动伪影的自动逆向修正方法，使得因刚性平移而损害的数据经修正后能重建出不含运动伪影的高质量的图像。
本发明的目的可通过以下的技术措施来实现，依次包括以下步骤：
(1)获取含有刚性平移运动伪影的磁共振图像的原始频域矢量矩阵数据，即K空间数据；
(2)以频率编码方向为Kx方向，相位编码方向为Ky方向，把相位编码梯度Gy为零时定为Ky等于零确定K空间坐标；运用傅立叶变换的周期特性，将所获各原始频域矢量矩阵转换成在K空间坐标中以相位编码方向Ky等于零为对称轴的上下共扼对称的矩阵数据；并获得呈矢量共扼对称的矩阵中每个频域矢量的相位，即K空间坐标中每个K空间线的原相位；
(3)以K空间坐标中Ky等于零为参照K空间线，沿坐标轴向正方向的第一个K空间线开始，逐个修正每一个K空间线已偏移的相位；
(4)除设定为必须要逆向估计线间运动距离的第一个K空间线之外，对要修正的K空间线，首先确定该K空间线是否需要进行逆向估计线间运动距离dx和dy，所述的线间运动距离dx和dy为两个K空间线线间的运动距离，即所需修正的K空间线相对于前一已修正好的K空间线的运动距离；
(5)若需要逆向估计线间运动距离dx和dy的K空间线时，采用逆向修正得出所求的线间运动距离dx和dy，然后根据此线间运动距离dx和dy加上前一已被修正好的K空间线的实际运动距离x_和y_得出该被修正K空间线的实际运动距离x和y，据此实际运动距离x和y修正该K空间线；
(6)若不需要逆向估计线间运动距离dx和dy的K空间线时，可直接把前一已修正好的K空间线的实际运动距离x_和y_加上该被修正K空间线前两个已修正好K空间线的实际运动距离之差作为该被修正K空间线的实际运动距离x和y，然后据此实际运动距离x和y修正该K空间线；
(7)重复步骤4、5、6中的处理过程，修正全部沿坐标轴向正方向的K空间线，对Ky等于零为参照的向负方向的各个K空间线的修正，以Ky等于零为对称轴取正方向对称行的共扼矢量即可；
(8)将所有修正好的K空间线的频域矢量运用傅立叶变换的周期特性，变换成跟初始矩阵一致的非共扼对称矩阵，然后通过傅立叶反变换变换到图像域，并将矢量转化为向量，得到修正好的图像。
本发明中判断K空间线是否需要进行逆向估计线间运动距离dx和dy的具体步骤如下：
(1)先将该未修正该K空间线的矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，获得该图像的熵；
(2)再在前一已修正好的K空间线的实际运动距离x_和y_的基础上加上该被修正K空间线前两个已修正好K空间线的实际运动距离之差，得到该被修正K空间线相对于Ky等于零的实际运动距离x和y，据此修正好该K空间线，所述的实际运动距离为相对于Ky等于零的运动距离；将修正后所得矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，获得该图像的熵；
(3)将上述步骤2中获得的熵值与步骤1中获得的未修正该K空间线时的熵值比较，如果熵值增大，则此K空间线需要进行逆向估计线间运动距离dx和dy。
本发明步骤5中所述逆向修正得出所求的线间运动距离dx和dy具体步骤如下：
(1)在所需修正的K空间线的前一已被修正好的K空间线的实际运动距离x和y的基础上加上估计的线间运动距离dx′和dy′，得到估计的所需修正的K空间线的实际运动距离x′和y′，对所需修正的K空间线用每次估计的实际运动距离x′和y′来进行修正；
(2)将上述每次估计修正后所得矩阵再次运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，求出所得图像的熵；在所得所有的熵中找出最小值，该最小熵值所对应的估计的线间运动距离dx′和dy′就是所求的线间运动距离dx和dy。
为了使得估计的线间运动距离dx′和dy′估计得更加精确，本发明所述得出线间运动距离dx和dy的过程中，在修正需要逆向估计线间运动距离dx和dy的K空间线的同时，沿坐标轴向正方向对紧接该K空间线的余下未修正的K空间线中依次选取一部分也分别同时进行估计修正，估计修正时每个选取的K空间线的估计的实际运动距离分别为前一K空间线的估计的实际运动距离加上与估计的线间运动距离dx′和dy′，对估计修正每个选取的的K空间线用各自每次估计的实际运动距离x′和y′来进行估计修正。这样，在每次逆向迭代修正的过程中，就相当于有多个K空间线同时被估计修正了，以便使熵值之间的差值增大，免受积累误差和计算误差的影响，从而使得最小熵值所对应的估计的线间运动距离dx′和dy′更加接近所求的线间运动距离dx和dy的实际值。
本发明中所述根据所得实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向的实际运动距离x和所在K空间坐标Kx的乘积与Y方向的实际运动距离y和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小。
本发明中所述根据所得估计的实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向估计的实际运动距离x′和所在K空间坐标Kx的乘积与Y方向估计的实际运动距离y′和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小。
本发明中所述获得熵值的具体过程为：得出每个象素点灰度值Bj与Bmax的比值，其中Bmax为图像所有像素的灰度值平方和开根而得，然后用每个象素算得的比值乘以该比值的以e为底的对数值，将图像中每个象素所得的乘积相加，对所得的和取负值就得到该图像的熵。
本发明方法与现有技术相比具有以下的优点：
1、由于本发明方法能自动搜索出运动轨迹，因此不需知道运动的先验信息；即不需知道运动的形式和过程 2、由于本发明方法是对任意方向进行逆向的估计，因此不管是周期的还是任意的都可估计出来；3、由于本发明方法是对整个图像区域进行自动修正，不用人工给出感兴趣区(ROI)或人工勾画出近似模板图像。
附图说明
图1本发明的流程图；
图2模板原始图像；
图3模拟刚性平移运动后重建出的伪影图像；
图4运用本发明所提方法消除运动伪影后重建出的图像。
具体实施方式：
实施例一
含有刚性平移运动伪影的磁共振图像的原始频域矢量矩阵数据可从磁共振成像系统中直接获得，也可将原始模板图像根据MR半傅立叶成像策略模拟刚性平移运动来获得。如图2，大小为256*256的原始模板图像根据MR半傅立叶成像策略模拟刚性平移运动得到图3的伪影图像的频域矢量矩阵数据，模拟过程中x和y方向的运动范围在[-10，10](单位：像素)之间，对伪影图像的修正依次按如下步骤：
1、将所得频域矢量矩阵数据以频率编码方向为Kx方向，相位编码方向为Ky方向，把相位编码梯度Gy为零时定为Ky等于零确定K空间坐标；运用傅立叶变换的周期特性，将所获各原始频域矢量矩阵转换成在K空间坐标中以相位编码方向Ky等于零为对称轴的上下对称即的矩阵数据；并获得呈矢量共扼对称的矩阵中每个频域矢量的相位，即K空间坐标中每个K空间线的原相位；
2、以K空间坐标中Ky等于零为参照K空间线，沿坐标轴向正方向的第一个K空间线开始，逐个修正每一个K空间线已偏移的相位；
3、修正沿坐标轴向正方向的第一个K空间线，该第一个K空间线必须要逆向估计线间运动距离dx和dy，采用逆向修正得出所求的线间运动距离dx和dy，具体过程如下：在所需修正的第一个K空间线的前一已被修正好的K空间线的实际运动距离x和y的基础上加上估计的线间运动距离dx′和dy′，对应第一个所需修正的K空间线时前一已被修正好的K空间线的实际运动距离x_和y_为零，取dx′和dy′的估计值分别为[-3，-2，-1，0，1，2，3]，得到7*7共49个估计的所需修正的第一个K空间线的估计的实际运动距离x′和y′，对所需修正的第一个K空间线用每次估计的实际运动距离x′和y′进行估计修正，根据所得估计的实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向估计的实际运动距离x′和所在K空间坐标Kx的乘积与Y方向估计的实际运动距离y′和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小。
将每次估计修正后所得矩阵再次运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，求出49次估计修正后对应所得的49个估计修正图像的熵。在所得所有的熵中找出最小值，它所对应的估计的线间运动距离dx′和dy′就是所求的线间运动距离dx和dy。然后根据此线间运动距离dx和dy与前一已被修正好的K空间线的实际运动距离x和y(零)，得出该被修正第一个K空间线的实际运动距离x和y，据此修正该第一个K空间线。
根据所得实际运动距离x和y修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向的实际运动距离x和所在K空间坐标Kx的乘积与Y方向的实际运动距离y和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小256。
4、对要修正的第二个K空间线，首先判断它是否需要进行估计线间运动距离dx和dy，这里dx和dy为两个K空间线线间的运动距离，即所需修正的K空间线相对于前一已修正好的K空间线的运动距离。
判断步骤如下：先将该未修正的第二个K空间线的矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，计算该图像的熵；
再在前一已修正好的K空间线的实际运动距离x_和y_(此运动距离是相对于Ky等于零的运动距离)的基础上加上该被修正K空间线前两个已修正好K空间线的运动实际运动距离之差，得到该被修正K空间线相对于Ky等于零的实际运动距离x和y，据此修正好该第二个K空间线，将所得矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，计算该图像的熵；
上述获得熵值的具体过程为：得出每个象素点灰度值Bj与Bmax的比值，其中Bmax为图像所有像素的灰度值平方和开根而得，然后用每个象素算得的比值乘以该比值的以e为底的对数值，将图像中每个象素所得的乘积相加，对所得的和取负值就得到该图像的熵。
将上述得到的两个修正前后熵值比较，如果熵值增大，则此第二K空间线需要进行逆向估计线间运动距离dx和dy，其修正过程重复上述步骤3即可。
5、如果熵值减小，修正不需要估计线间运动距离dx和dy的K空间线，可直接把前一已修正好的K空间线的实际运动距离x_和y_(即上述已修正好的第一K空间线的实际运动距离)加上该被修正K空间线前两个已修正好K空间线的运动实际运动距离之差(即上述已修正好的第一K空间线的实际运动距离与参照空间线的零运动距离之差)作为该被修正K空间线的实际运动距离x和y，然后据此修正该第二K空间线。
根据所得实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向的实际运动距离x和所在K空间坐标Kx的乘积与Y方向的实际运动距离y和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小256。
6、重复步骤4、5中的处理过程，修正全部沿坐标轴向正方向的K空间线，对Ky等于零为参照的向负方向的各个K空间线的修正，以Ky等于零为对称轴取正方向对称行的共扼矢量即可；
7、将所有修正好的K空间线的频域矢量运用傅立叶变换的周期特性，变换成跟初始矩阵一致的非共扼对称矩阵，然后通过傅立叶反变换变换到图像域，并将矢量转化为向量，得到修正好的图像。
实施例二
步骤1、2与实施例一相同；
步骤3，为了使估计出的线间运动距离dx′和dy′估计得更加精确，根据运动一致性和连续性，在上述实施例一的步骤3中修正需要逆向估计线间运动距离dx和dy的第一K空间线的同时，对余下未修正的K空间线中紧接该第一个K空间线沿坐标轴向正方向依次选取即将被修正的七个K空间线也分别依次进行估计修正，选取余下未修正的K空间线的多少可任意，以满足图像的熵值变化不受积累误差和计算误差的影响为准。估计修正时每个选取的K空间线的估计的实际运动距离分别为前一K空间线的估计的实际运动距离加上估计的线间运动距离dx′和dy′，对估计修正每个选取的的K空间线用各自每次估计的实际运动距离x′和y′来进行修正。这样加上该被修正的第一个K空间线为8个，根据正负共扼对称，负方向也为8个，所以每次估计修正时同时对16个K空间线进行了估计修正。经上述每次估计修正后相应得到的49个矩阵在变换成图像后，所得的熵值之间的差值会变大，免受积累误差和计算误差的影响，从而使得最小熵值所对应的估计的线间运动距离dx′和dy′更加接近所求的线间运动距离dx和dy的实际值。
然后根据上述得出的更加接近实际值的线间运动距离dx和dy与前一已被修正好的K空间线的实际运动距离x和y(零)，得出该被修正第一个K空间线的实际运动距离x和y，据此修正该第一个K空间线。
余下步骤的处理过程与实施例一相同，得到修正好的图像(如图4)。
Claims (7)
1、一种磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：依次包括以下步骤：
(1)获取含有刚性平移运动伪影的磁共振图像的原始频域矢量矩阵数据，即K空间数据；
(2)以频率编码方向为Kx方向，相位编码方向为Ky方向，把相位编码梯度Gy为零时定为Ky等于零确定K空间坐标；运用傅立叶变换的周期特性，将所获各原始频域矢量矩阵转换成在K空间坐标中以相位编码方向Ky等于零为对称轴的上下共扼对称的矩阵数据；并获得呈矢量共扼对称的矩阵中每个频域矢量的相位，即K空间坐标中每个K空间线的原相位；
(3)以K空间坐标中Ky等于零为参照K空间线，沿坐标轴向正方向的第一个K空间线开始，逐个修正每一个K空间线已偏移的相位；
(4)除设定为必须要逆向估计线间运动距离的第一个K空间线之外，对要修正的K空间线，首先判断该K空间线是否需要进行逆向估计线间运动距离dx和dy，所述的线间运动距离dx和dy为两个K空间线线间的运动距离，即所需修正的K空间线相对于前一已修正好的K空间线的运动距离；
(5)若需要逆向估计线间运动距离dx和dy的K空间线时，采用逆向修正得出所求的线间运动距离dx和dy，然后根据此线间运动距离dx和dy加上前一已被修正好的K空间线的实际运动距离x_和y_得出该被修正K空间线的实际运动距离x和y，据此实际运动距离x和y修正该K空间线；
(6)若不需要逆向估计线间运动距离dx和dy的K空间线时，可直接把前一已修正好的K空间线的实际运动距离x_和y_加上该被修正K空间线前两个已修正好K空间线的实际运动距离之差作为该被修正K空间线的实际运动距离x和y，然后据此实际运动距离x和y修正该K空间线；
(7)重复步骤4、5、6中的处理过程，修正全部沿坐标轴向正方向的K空间线，对Ky等于零为参照的向负方向的各个K空间线的修正，以Ky等于零为对称轴取正方向对称行的共扼矢量即可；
(8)将所有修正好的K空间线的频域矢量运用傅立叶变换的周期特性，变换成跟初始矩阵一致的非共扼对称矩阵，然后通过傅立叶反变换变换到图像域，并将矢量转化为向量，得到修正好的图像。
2、根据权利要求1所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：在步骤4中判断K空间线是否需要进行逆向估计线间运动距离dx和dy的具体步骤如下：
(1)先将该未修正该K空间线的矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，获得该图像的熵；
(2)再在前一已修正好的K空间线的实际运动距离x_和y_的基础上加上该被修正K空间线前两个已修正好K空间线的实际运动距离之差，得到该被修正K空间线相对于Ky等于零的实际运动距离x和y，据此修正好该K空间线，所述的实际运动距离为相对于Ky等于零的运动距离；将修正后所得矩阵运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，获得该图像的熵；
(3)将上述步骤2中获得的熵值与步骤1中获得的未修正该K空间线时的熵值比较，如果熵值增大，则此K空间线需要进行逆向估计线间运动距离dx和dy。
3、根据权利要求1所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：步骤5中所述逆向修正得出所求的线间运动距离dx和dy具体步骤如下：
(1)在所需修正的K空间线的前一已被修正好的K空间线的实际运动距离x和y的基础上加上估计的线间运动距离dx′和dy′，得到估计的所需修正的K空间线的估计的实际运动距离x′和y′，对所需修正的K空间线用每次估计的实际运动距离x′和y′来进行修正；
(2)将上述每次修正后所得矩阵再次运用傅立叶变换的周期特性，变换成跟初始各矩阵一致的呈非共扼对称的矩阵，并将这些呈非共扼对称的矩阵通过傅立叶反变换变换到图像域，最后将矢量转换成向量而得图像，求出所得图像的熵；在所得所有的熵中找出最小值，该最小熵值所对应的估计的线间运动距离dx′和dy′就是所求的线间运动距离dx和dy。
4、根据权利要求3所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：所述得出线间运动距离dx和dy的过程中，在修正需要逆向估计线间运动距离dx和dy的K空间线的同时，沿坐标轴向正方向对紧接该K空间线的余下未修正的K空间线中依次选取一部分也分别依次进行估计修正，估计修正时每个选取的K空间线的估计的实际运动距离分别为前一K空间线的估计的实际运动距离加上估计的线间运动距离dx′和dy′，对估计修正每个选取的的K空间线用各自每次估计的实际运动距离x′和y′来进行估计修正。
5、根据权利要求1或2所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：所述根据所得实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向的实际运动距离x和所在K空间坐标Kx的乘积与Y方向的实际运动距离y和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小。
6、根据权利要求3或4所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：所述根据所得估计的实际运动距离修正K空间线的具体过程为：该K空间线数据的幅值保持不变，相位为在步骤2中所得的该K空间线原相位的基础上加上一个值，该值为X方向估计的实际运动距离x′和所在K空间坐标Kx的乘积与Y方向估计的实际运动距离y′和所在K空间坐标Ky的乘积的和，再乘上一个常量C；其中常量C为2π/N，N为图像大小。
7、根据权利要求2或3所述的磁共振图像刚性平移运动伪影的自动逆向修正方法，其特征在于：所述获得熵值的具体过程为：得出每个象素点灰度值Bj与Bmax的比值，其中Bmax为图像所有像素的灰度值平方和开根而得，然后用每个象素算得的比值乘以该比值的以e为底的对数值，将图像中每个象素所得的乘积相加，对所得的和取负值就得到该图像的熵。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB03140278XA CN1210671C (zh) | 2003-08-27 | 2003-08-27 | 磁共振图像刚性平移运动伪影的自动逆向修正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB03140278XA CN1210671C (zh) | 2003-08-27 | 2003-08-27 | 磁共振图像刚性平移运动伪影的自动逆向修正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1515912A true CN1515912A (zh) | 2004-07-28 |
| CN1210671C CN1210671C (zh) | 2005-07-13 |
Family
ID=34240271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB03140278XA Expired - Fee Related CN1210671C (zh) | 2003-08-27 | 2003-08-27 | 磁共振图像刚性平移运动伪影的自动逆向修正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1210671C (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1327809C (zh) * | 2005-03-28 | 2007-07-25 | 南方医科大学 | 在t1加权磁共振成像中用propeller采样方式消除运动伪影的方法 |
| CN101190128B (zh) * | 2006-11-30 | 2010-05-19 | Ge医疗系统环球技术有限公司 | 采集磁共振成像数据的方法和设备 |
| CN105701846A (zh) * | 2014-11-27 | 2016-06-22 | 株式会社东芝 | 医学图像处理装置和方法以及医学成像设备 |
| CN107607895A (zh) * | 2016-12-30 | 2018-01-19 | 上海联影医疗科技有限公司 | 磁共振成像方法和系统 |
| CN107621617A (zh) * | 2013-04-27 | 2018-01-23 | 上海联影医疗科技有限公司 | k空间运动伪影矫正装置 |
| CN108291950A (zh) * | 2015-12-03 | 2018-07-17 | 皇家飞利浦有限公司 | Sense-mri中的图像伪影的去除 |
| CN109765513A (zh) * | 2019-02-12 | 2019-05-17 | 上海联影医疗科技有限公司 | 用于产生磁共振场图的方法、装置及计算机设备 |
| US20220187406A1 (en) * | 2019-03-12 | 2022-06-16 | University Of Cincinnati | A system and method for motion correction of magnetic resonance image |
| US20220196771A1 (en) * | 2019-04-25 | 2022-06-23 | Insightec, Ltd. | Accelerated magnetic resonance thermometry |
-
2003
- 2003-08-27 CN CNB03140278XA patent/CN1210671C/zh not_active Expired - Fee Related
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1327809C (zh) * | 2005-03-28 | 2007-07-25 | 南方医科大学 | 在t1加权磁共振成像中用propeller采样方式消除运动伪影的方法 |
| CN101190128B (zh) * | 2006-11-30 | 2010-05-19 | Ge医疗系统环球技术有限公司 | 采集磁共振成像数据的方法和设备 |
| CN107621617B (zh) * | 2013-04-27 | 2019-12-20 | 上海联影医疗科技有限公司 | k空间运动伪影矫正装置 |
| CN107621617A (zh) * | 2013-04-27 | 2018-01-23 | 上海联影医疗科技有限公司 | k空间运动伪影矫正装置 |
| CN105701846A (zh) * | 2014-11-27 | 2016-06-22 | 株式会社东芝 | 医学图像处理装置和方法以及医学成像设备 |
| CN105701846B (zh) * | 2014-11-27 | 2019-01-15 | 东芝医疗系统株式会社 | 医学图像处理装置和方法以及医学成像设备 |
| CN108291950A (zh) * | 2015-12-03 | 2018-07-17 | 皇家飞利浦有限公司 | Sense-mri中的图像伪影的去除 |
| CN107607895A (zh) * | 2016-12-30 | 2018-01-19 | 上海联影医疗科技有限公司 | 磁共振成像方法和系统 |
| CN109765513A (zh) * | 2019-02-12 | 2019-05-17 | 上海联影医疗科技有限公司 | 用于产生磁共振场图的方法、装置及计算机设备 |
| CN109765513B (zh) * | 2019-02-12 | 2021-05-18 | 上海联影医疗科技股份有限公司 | 用于产生磁共振场图的方法、装置及计算机设备 |
| US20220187406A1 (en) * | 2019-03-12 | 2022-06-16 | University Of Cincinnati | A system and method for motion correction of magnetic resonance image |
| US11835612B2 (en) * | 2019-03-12 | 2023-12-05 | University Of Cincinnati | System and method for motion correction of magnetic resonance image |
| US20220196771A1 (en) * | 2019-04-25 | 2022-06-23 | Insightec, Ltd. | Accelerated magnetic resonance thermometry |
| US11921182B2 (en) * | 2019-04-25 | 2024-03-05 | Insightec Ltd. | Accelerated magnetic resonance thermometry |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1210671C (zh) | 2005-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110490832B (zh) | 一种基于正则化深度图像先验方法的磁共振图像重建方法 | |
| CN108460726B (zh) | 一种基于增强递归残差网络的磁共振图像超分辨重建方法 | |
| US7202663B2 (en) | Method for generating fast magnetic resonance images | |
| CN1210671C (zh) | 磁共振图像刚性平移运动伪影的自动逆向修正方法 | |
| CN110458756A (zh) | 基于深度学习的模糊视频超分辨率方法及系统 | |
| EP1980865A1 (en) | Rapid MRI dynamic imaging using MACH | |
| CN110827232B (zh) | 基于形态特征gan的跨模态mri合成方法 | |
| CN1686050A (zh) | 在t1加权磁共振成像中用propeller采样方式消除运动伪影的方法 | |
| CN113160380B (zh) | 三维磁共振影像超分辨重建方法、电子设备和存储介质 | |
| CN104181486A (zh) | 磁共振图像重建方法及装置 | |
| CN107038730A (zh) | 基于高斯尺度结构块分组的稀疏表示图像重建方法 | |
| US7778493B2 (en) | Pixelation reconstruction for image resolution and image data transmission | |
| CN112396672A (zh) | 一种基于深度学习的稀疏角度锥束ct图像重建方法 | |
| CN106952226A (zh) | 一种f‑msa超分辨率重建方法 | |
| CN1663526A (zh) | 磁共振图像重建时所产生Gibbs环状伪影的消除方法 | |
| Jiang et al. | Respiratory motion correction in abdominal MRI using a densely connected U-Net with GAN-guided training | |
| CN105678822A (zh) | 一种基于Split Bregman迭代的三正则磁共振图像重构方法 | |
| CA2468803A1 (en) | Method for compensating for effects of object motion in mri | |
| CN104181484A (zh) | 磁共振图像重建方法 | |
| CN107942271A (zh) | 基于迭代的speed快速磁共振成像方法 | |
| Iwamoto et al. | Super-resolution of MR volumetric images using sparse representation and self-similarity | |
| CN116758120A (zh) | 一种基于深度学习的3t mra到7t mra的预测方法 | |
| CN114596378A (zh) | 一种稀疏角ct伪影去除的方法 | |
| CN113192155A (zh) | 螺旋ct锥束扫描图像重建方法、扫描系统及存储介质 | |
| CN103860173B (zh) | 磁共振k空间重建方法与装置、磁共振成像方法与装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| ASS | Succession or assignment of patent right |
Owner name: GUANGDONG SOUTHERN MEDICAL UNIVERSITY ASSET MANAG Free format text: FORMER OWNER: NANFANG MEDICAL UNIV. Effective date: 20080725 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| C56 | Change in the name or address of the patentee |
Owner name: NANFANG MEDICAL UNIV. Free format text: FORMER NAME OR ADDRESS: NO.1 MMILITARY MEDICAL UNIVERSITY, CHINESE PLA |
|
| CP03 | Change of name, title or address |
Address after: No. 1023-1063 Sha Ha Road, Guangzhou Patentee after: SOUTHERN MEDICAL University Address before: Department of Biomedical Engineering, First Military Medical University, Guangzhou Tonghe Road, Guangdong, China Patentee before: No.1 Military Medical Univ., PLA |
|
| TR01 | Transfer of patent right |
Effective date of registration: 20080725 Address after: Floor two, attached building, science and technology building, Southern Medical University, Baiyun District, Guangzhou Patentee after: Guangdong Nanfang Medical University Asset Management Ltd. Address before: No. 1023-1063 Sha Ha Road, Guangzhou Patentee before: Southern Medical University |
|
| ASS | Succession or assignment of patent right |
Owner name: GUANGZHOU YI CHENG DIGITAL MEDICAL SYSTEM CO., LT Free format text: FORMER OWNER: GUANGDONG SOUTHERN MEDICAL UNIVERSITY ASSET MANAGEMENT CO., LTD. Effective date: 20080808 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20080808 Address after: Guangzhou economic and Technological Development Zone Science Town No. 11, No. 913, building B Patentee after: Guangzhou Yicheng Digital Medical System Co.,Ltd. Address before: Floor two, attached building, science and technology building, Southern Medical University, Baiyun District, Guangzhou Patentee before: Guangdong Nanfang Medical University Asset Management Ltd. |
|
| DD01 | Delivery of document by public notice |
Addressee: Chen Ni Document name: payment instructions |
|
| DD01 | Delivery of document by public notice | ||
| DD01 | Delivery of document by public notice |
Addressee: Chen Ni Document name: Notice of termination of patent right |
|
| DD01 | Delivery of document by public notice | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050713 Termination date: 20200827 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |