CN104778341B - 磁共振线圈合并系数计算方法、磁共振成像方法及其装置 - Google Patents

磁共振线圈合并系数计算方法、磁共振成像方法及其装置 Download PDF

Info

Publication number
CN104778341B
CN104778341B CN201410010303.1A CN201410010303A CN104778341B CN 104778341 B CN104778341 B CN 104778341B CN 201410010303 A CN201410010303 A CN 201410010303A CN 104778341 B CN104778341 B CN 104778341B
Authority
CN
China
Prior art keywords
data
obk
gathered
magnetic resonance
refinement
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
Application number
CN201410010303.1A
Other languages
English (en)
Other versions
CN104778341A (zh
Inventor
翟人宽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201410010303.1A priority Critical patent/CN104778341B/zh
Priority to CN201710368811.0A priority patent/CN107092805B/zh
Publication of CN104778341A publication Critical patent/CN104778341A/zh
Application granted granted Critical
Publication of CN104778341B publication Critical patent/CN104778341B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Epidemiology (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种磁共振并行成像线圈合并系数计算方法、磁共振并行成像方法及其装置，磁共振并行成像线圈合并系数计算方法包括如下步骤：a)采集k空间数据，所述k空间数据包括采集数据和校准数据；b)对所述采集数据进行数据提炼得到提炼数据；c)由所述提炼数据和所述校准数据计算得到线圈合并系数。通过本发明所提供的技术方案，针对磁共振并行采集重建过程中选取的卷积核，进行数据重提炼，使得提炼后的数据量减小，信号特征得到增强，利用提炼后的卷积核，进行数据填补，得到新的卷积核比原有的得到了优化。

Description

磁共振线圈合并系数计算方法、磁共振成像方法及其装置
技术领域
本发明涉及磁共振成像领域，特别涉及一种磁共振线圈合并系数计算方法、磁共振成像方法及其装置。
背景技术
在磁共振成像技术中，成像的速度是衡量成像方法的一个很重要标准。限制成像速度的很重要因素是数据采集，以及k空间填充。一般的数据采集方式要采满k空间数据，然后才能进行重建得到图像。磁共振并行采集重建技术，是利用线圈重组合并的方式，对欠采样的数据进行填补，利用填补完整的k空间数据进行重建。利用这样的方式，可以根据需求，只采集一部分k空间数据，不必采满整个k空间。因此这样的方法可以大大加快成像的速度。
比较常用的并行重建方法之一是GRAPPA。传统的GRAPPA算法如图1所示，黑色实点代表为实际采集的k空间数据；白色空点为欠采样需要填补的数据；灰色实点表示为了计算线圈合并参数，而适量全采的校准数据。GRAPPA算法认为，图中任意一个空心点可以表示为周围黑色实点的线性叠加，相当于对多个线圈的数据进行了合并。而线圈合并系数nij(第i个线圈，第j个位置，如图1)可以通过黑色的实点拟合灰色点来确定。线圈合并系数确定后，其他白色空心点即可根据求得的合并系数以及黑色实点计算填补得到，从而重建得到完整的k空间数据。
线圈合并系数又可称为卷积核，在传统方法中，卷积核的计算方向，只加入了相位编码方向，以及通道方向。为了优化效果，近年来，很多方法引入了其他方向如频率编码方向，以及k-t中时间方向等。这些维度的引入，使卷积核的信息增加，但是随着过多的引入，信息会出现冗余，而且由于噪声的影响，使得计算的系数仍然受到影响；而卷积核过大，使得要拟合的系数增多，也带来计算的不稳定性。而卷积核过小，则数据信息不足，计算过程不够准确。
发明内容
本发明要解决的问题是提供一种磁共振并行成像线圈合并系数计算方法、磁共振并行成像方法及其装置，解决在卷积核选取过程中的由于卷积核选取过大造成信息量增加，信息冗余，并且噪声过多，拟合系数增多，计算不稳定的影响，以及在卷积核选取过小时，计算数据不足，结果不准确的问题。
为了实现上述目的，本发明提供了一种磁共振并行成像线圈合并系数计算方法，包括如下步骤：
a)采集k空间数据，所述k空间数据包括采集数据和校准数据；
b)对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n；
c)由所述提炼数据和所述校准数据计算得到线圈合并系数。
优选的，所述步骤a)中k空间数据的数据方向包括以下任意一种或多种：频率编码方向、相位编码方向、通道方向、时间方向、3d扫描中的第二相位编码方向或上述多个方向合并形成的方向。
优选的，所述步骤b)中数据提炼方法为矩阵降维法。
优选的，所述步骤b)中矩阵降维法包括如下步骤：
将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；
计算所述采集数据的协方差矩阵，并选取所述协方差矩阵归一化正交特征矢量中特征值较大的前t个值，组成提炼系数；
由所述提炼系数计算得到提炼数据。
优选的，所述步骤b)中矩阵降维法包括如下步骤：
将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；
按照相同数据方向将所述采集数据分类得到各类第一方向向量Obk_Sx1、Obk_Sx2、…、Obk_Sxi；
由所述各类第一方向向量计算各类样本第一方向向量均值Obk_Sm1、Obk_Sm2、…、Obk_Smi和总样本第一方向向量均值Obk_Sm，并进而计算得到样本类内离散度矩阵Obk_Ss1、Obk_Ss2、…、Obk_Ssi以及样本间离散度矩阵Obk_Ssb；
Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量分别为wi1、wi2、…、wit，将wi1、wi2至wit组合成一个矩阵作为所述提炼数据。
进一步的，所述第一方向为所述采集数据的列向量方向或行向量方向。
本发明还提供了一种磁共振并行成像方法，包括在所述的磁共振并行成像线圈合并系数计算方法计算得到线圈合并系数后，根据所述线圈合并系数和所述采集数据，填补得到完整的k空间数据，将所述完整的k空间数据变换到图像域得到图像。
本发明还提供了一种磁共振并行成像装置，包括：
采集单元，适于采集k空间数据，所述k空间数据包括采集数据和校准数据，所述k空间数据至少包含一个数据方向；
计算单元，适于对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n，并由所述提炼数据和所述校准数据计算得到线圈合并系数；
填补单元，适于根据所述线圈合并系数和所述采集数据，填补得到完整的k空间数据；
成像单元，适于将所述完整的k空间数据变换到图像域得到图像。
优选的，所述计算单元适于：将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；计算所述采集数据的协方差矩阵，并选取所述协方差矩阵归一化正交特征矢量中特征值较大的前t个值，组成提炼系数；由所述提炼系数计算得到提炼数据。
优选的，所述计算单元适于：将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；按照相同数据方向将所述采集数据分类得到各类第一方向向量Obk_Sx1、Obk_Sx2、…、Obk_Sxi；由所述各类第一方向向量计算各类样本第一方向向量均值Obk_Sm1、Obk_Sm2、…、Obk_Smi和总样本第一方向向量均值Obk_Sm，并进而计算得到样本类内离散度矩阵Obk_Ss1、Obk_Ss2、…、Obk_Ssi以及样本间离散度矩阵Obk_Ssb；Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量分别为wi1、wi2、…、wit，将wi1、wi2至wit组合成一个矩阵作为所述提炼数据。
通过本发明所提供的技术方案，针对磁共振并行采集重建过程中选取的卷积核，进行数据重提炼，使得提炼后的数据量减小，信号特征得到增强，利用提炼后的卷积核，进行数据填补，得到新的卷积核比原有的得到了优化。这样一方面可以一定程度去除噪声的影响，使得信号特征强化，优化系数计算的准确性，稳定性；更进一步，可以使得卷积核的选取更加容易。
附图说明
图1为现有技术线圈合并系数计算示意图；
图2为本发明磁共振并行成像线圈合并系数计算方法流程图。
具体实施方式
为了使本发明的上述目的、特征、优点能够更为显而易懂，下面结合附图和实施例对本发明的具体实施方式作进一步描述。
本发明提供了一种磁共振并行成像线圈合并系数计算方法，图2为其流程图，请参见图2，所述方法包括如下步骤：
S101，采集k空间数据，所述k空间数据包括采集数据和校准数据；
S102，对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n；
S103，由所述提炼数据和所述校准数据计算得到线圈合并系数。
首先执行步骤S101，对k空间数据进行采集，在传统的GRAPPA技术中，其采集得到的数据分为采集数据和校准数据。其中，采集数据由并行采集方式采集得到，如图1所示，其采集数据在k空间中为隔行的数据，空白部分为未采集到的欠采集数据，而校准数据则是通过全采集得到的k空间中心区域部分的数据。在计算线圈合并系数时，上述采集数据和校准数据之间关系为：
Obk_S*Cft=Obk_Sacs [1]
其中Obk_S表示采集数据，Obk_Sacs表示校准数据，Cft表示GRAPPA算法中待求的线圈合并系数。
如背景资料部分所述，采集得到的k空间数据通常包含多个数据方向，传统数据方向为相位编码方向和通道方向，近年来为了扩展所述数据，还使用到了频率编码方向、k-t中的时间方向，以及在进行3d扫描时的第二相位编码方向。K空间数据的数据方向也可以是上述一种或多种方向合并形成的方向。对于本发明的实施方式，之后的描述将介绍数据方向的数目并不影响本发明的实施。
接下来执行S102，对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n。
这里，需要将采集数据Obk_S进行提炼，得到提炼数据Dbk_S，采集数据Obk_S和提炼数据Dbk_S之间的关系为：
Dbk_S=Obk_S*Rft [2]
其中，Rft为提炼系数。
确定提炼系数方法的思想为，将采集数据Obk_S中的第一方向向量作为基础向量，进行主成分分析提取，去向量之间的相关性，这里的第一方向向量可以为采集数据矩阵中每一列的列向量或者每一行的行向量。这里选择的方法为矩阵降维的数学方法，如PCA，KLT，LDA等方法中任意的一种方法。
KLT算法为特征提取比较常用的算法之一[参见文献：1.M.Turk and A.Pentland，“Eigenfaces for recognition，”J.Cogn.Neurosci.3，71-86(1991).2.R.Everson and L.“Sirovich，Karhunen-Loeve procedure for gappy data”Vol.12，No.8／August1995／J.Opt.Soc.Am.A]。
以KLT算法为例进行采集数据矩阵的降维，进而确定提炼系数Rft时，其实施过程为将采集数据Obk_S的列向量作为提取对象，该列向量相当于图1中对应的一列黑色点数据。首先求得采集数据Obk_S的协方差矩阵C_Obk_S，然后求得C_Obk_S的归一化正交特征矢量q(假定共n个)，选取特征值较大的前t(t<n)个q，组成Rft，此即提炼系数。上述过程中，是以列向量作为基础向量，这里同样可以将行向量作为基础向量。在以列向量作为基础向量时，其与列向量垂直的方向记为提取方向(行方向)，在提取方向上采集数据有n个维度，由于提炼系数Rft有t个维度，经过提取后得到的提取数据在提取方向(行方向)上也是t个维度，其中选取的t<n。
若以LDA算法为例，本发明技术的具体实施过程为：首先将采集数据Obk_S中的列向量进行分类，按照数据方向(如：频率编码方向和通道方向)相同的向量作为一类记为Obk_Sx1、Obk_Sx2、…、Obk_Sxi，所述采集矩阵共c个类，N个向量。第i类的所有列向量表示为Ri，其共Ni个列向量。其中，c、N、i、Ni均为正整数。
记各类样本均值向量为Obk_Sm1、Obk_Sm2、…、Obk_Smi，其计算过程为：
总样本均值Obk_Sm，其计算过程为：
样本类内离散度矩阵Obk_Ss1、Obk_Ss2、…、Obk_Ssi，其计算过程为：
样本类间离散度矩阵Obk_Ssb，其计算过程为：
上述各式中，x表示第i类列向量Ri中的一个列向量，c、N、i、Ni与之前规定相同，(…)T表示为矩阵的倒置。
记wi1，wi2…wit为矩阵Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量，将所有类求得的特征向量集合组合成一个矩阵即为提炼数据Dbk_S，利用式[2]可求得提炼系数Rft。
上述步骤中将列向量选取作为基础向量，与KLT算法相同的是这里也可以将行向量作为基础向量。在以列向量作为基础向量时，行方向即为提炼方向，在提炼方向上的数据维数为n，经过提炼之后的提炼数据的数据维数为t，t<n。
最后执行步骤S1OS，由所述提炼数据和所述校准数据计算得到线圈合并系数。
经过上一步骤的提炼，最终计算线圈合并系数基组包含提炼数据Dbk_S与校准数据Obk_Sacs这两组数据；
其提炼数据Dbk_S和校准数据Obk_Sacs相应关系为：
Dbk_S*Cft_new=Obk_Sacs [7]
这里即转换为了求解Cft_new的过程。
根据步骤S101采集得到的采集数据Obk_S和校准数据Obk_Sacs，以及相应公式[2]和[7]，可知：
Obk_S*Rft*Cft_new=Obk_Sacs [8]
为了计算简便，这里简化为：
Cft=Rft*Cft_new [9]
这样，则：
Obk_S*Cft=Obk_Sacs [10]
上述公式中Cft_new即为本发明技术方案得到的线圈合并系数，填补得到完整k空间的过程就是利用线圈合并系数Cft_new以及采集数据Obk_S进行。
通过以上对于本发明技术方案的描述可知，由于在步骤S102中，通过数据的提炼过程，提炼数据Dbk_S所对应的线圈合并系数Cft_new的维度相应得到了控制，所以在步骤S101选择采集数据的时候，可以简单的选择比较大的采集数据进行处理，不必根据数据的情况变换采集数据的大小。即使得选取采集数据更加简单，举例说明如下：普通的方法选取采集数据，假定为采集数据包含两个数据方向，为nx*ny(如3*4)大小；由于校准数据量有限，拟合过程中nx*ny不宜过大；但另一方面，我们希望nx*ny足够大，来涵盖尽量多的信息，这样需要平衡卷积核的大小。所以利用本方法，可以选取一个nx*ny(如30*4)比较大的一个采集数据，然后利用提炼的方法，将nx*ny中有用的信息提取出来(缩小nx*ny至3*4数据大小)，使其在提炼过程中，最终得到的提炼数据变小。
本发明在上述磁共振并行成像线圈合并系数计算方法的基础上，还提供了一种磁共振并行成像方法，包括由上述磁共振并行成像线圈合并系数计算方法计算得到的线圈合并系数，以及采集数据进行并行加速数据重建，填补得到完整的k空间数据，之后将k空间数据变换到图像域，得到磁共振图像。
本发明对应于磁共振成像方法还提供了一种磁共振并行成像装置，包括：
采集单元，适于采集k空间数据，所述k空间数据包括采集数据和校准数据；
计算单元，适于对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n，并由所述提炼数据和所述校准数据计算得到线圈合并系数；
填补单元，适于根据所述线圈合并系数和所述采集数据，填补得到完整的k空间数据；
成像单元，适于将所述完整的k空间数据变换到图像域得到图像。
其中，计算单元适于：将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；计算所述采集数据的协方差矩阵，并选取所述协方差矩阵归一化正交特征矢量中特征值较大的前t个值，组成提炼系数；由所述提炼系数计算得到提炼数据。
可选的，计算单元适于：将所述采集数据的第一方向向量作为基础向量，所述第一方向与提炼方向垂直；按照相同数据方向将所述采集数据分类得到各类第一方向向量Obk_Sx1、Obk_Sx2、…、Obk_Sxi；由所述各类第一方向向量计算各类样本第一方向向量均值Obk_Sm1、Obk_Sm2、…、Obk_Smi和总样本第一方向向量均值Obk_Sm，并进而计算得到样本类内离散度矩阵Obk_Ss1、Obk_Ss2、…、Obk_Ssi以及样本间离散度矩阵Obk_Ssb；Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量分别为wi1、wi2、…、wit，将wi1、wi2至wit组合成一个矩阵作为所述提炼数据。
上述磁共振成像方法和磁共振成像装置的具体实施过程可参考磁共振并行成像线圈合并系数计算方法的实施过程，这里不再一一赘述。
虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。

Claims (5)

1.一种磁共振并行成像线圈合并系数计算方法，其特征在于，包括如下步骤：
a)采集k空间数据，所述k空间数据包括采集数据和校准数据；
b)对所述采集数据进行数据提炼得到提炼数据，所述采集数据在提炼方向上的维数为n，所述提炼数据在所述提炼方向上的维数为t，所述t、n均为正整数且t<n；
c)由所述提炼数据和所述校准数据计算得到线圈合并系数；
所述步骤b)中数据提炼方法为矩阵降维法；
所述步骤b)中矩阵降维法包括如下步骤：
将所述采集数据的第一方向向量作为基础向量，所述提炼方向与所述第一方向垂直；
计算所述采集数据的协方差矩阵，并选取所述协方差矩阵归一化正交特征矢量中特征值较大的前t个值，组成提炼系数；
由所述提炼系数结合所述采集数据计算得到提炼数据。
2.如权利要求1所述的磁共振并行成像线圈合并系数计算方法，其特征在于，所述步骤a)中k空间数据的数据方向包括以下任意一种或多种：频率编码方向、相位编码方向、通道方向、时间方向、3d扫描中的第二相位编码方向或上述多个方向合并形成的方向。
3.如权利要求1所述的磁共振并行成像线圈合并系数计算方法，其特征在于，所述步骤b)中矩阵降维法包括如下步骤：
将所述采集数据的第一方向向量作为基础向量，所述提炼方向与所述第一方向垂直；
按照相同数据方向将所述采集数据分类得到各类第一方向向量Obk_Sx1、Obk_Sx2、…、Obk_Sxi；
由所述各类第一方向向量计算各类样本第一方向向量均值Obk_Sm1、Obk_Sm2、…、Obk_Smi和总样本第一方向向量均值Obk_Sm，并进而计算得到样本类内离散度矩阵Obk_Ss1、Obk_Ss2、…、Obk_Ssi以及样本间离散度矩阵Obk_Ssb；
Obk_Ssi-1*Obk_Ssb的前t个最大特征值对应的特征向量分别为wi1、wi2、…、wit，将wi1、wi2至wit组合成一个矩阵作为所述提炼数据。
4.如权利要求1所述的磁共振并行成像线圈合并系数计算方法，其特征在于，所述第一方向为所述采集数据的列向量方向或行向量方向。
5.一种磁共振并行成像方法，其特征在于，由权利要求1-4任一项所述的磁共振并行成像线圈合并系数计算方法计算得到线圈合并系数后，根据所述线圈合并系数和所述采集数据，计算并填补得到完整的k空间数据，将所述完整的k空间数据变换到图像域得到图像。
CN201410010303.1A 2014-01-09 2014-01-09 磁共振线圈合并系数计算方法、磁共振成像方法及其装置 Active CN104778341B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201410010303.1A CN104778341B (zh) 2014-01-09 2014-01-09 磁共振线圈合并系数计算方法、磁共振成像方法及其装置
CN201710368811.0A CN107092805B (zh) 2014-01-09 2014-01-09 磁共振并行成像装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410010303.1A CN104778341B (zh) 2014-01-09 2014-01-09 磁共振线圈合并系数计算方法、磁共振成像方法及其装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN201710368811.0A Division CN107092805B (zh) 2014-01-09 2014-01-09 磁共振并行成像装置

Publications (2)

Publication Number Publication Date
CN104778341A CN104778341A (zh) 2015-07-15
CN104778341B true CN104778341B (zh) 2017-08-22

Family

ID=53619802

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201710368811.0A Active CN107092805B (zh) 2014-01-09 2014-01-09 磁共振并行成像装置
CN201410010303.1A Active CN104778341B (zh) 2014-01-09 2014-01-09 磁共振线圈合并系数计算方法、磁共振成像方法及其装置

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201710368811.0A Active CN107092805B (zh) 2014-01-09 2014-01-09 磁共振并行成像装置

Country Status (1)

Country Link
CN (2) CN107092805B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107582057B (zh) * 2017-09-30 2020-10-27 上海联影医疗科技有限公司 磁共振成像方法及装置
CN110146835B (zh) * 2019-05-22 2021-09-07 山东颐邦齐鲁医生集团管理有限公司 一种基于并行成像的自导航磁共振图像重建方法及装置
CN111582453B (zh) * 2020-05-09 2023-10-27 北京百度网讯科技有限公司 生成神经网络模型的方法和装置
CN111951344B (zh) * 2020-08-09 2022-08-02 昆明理工大学 一种基于级联并行卷积网络的磁共振图像重建方法
CN112557980B (zh) * 2020-11-02 2022-05-03 上海东软医疗科技有限公司 磁共振图像矫正方法、装置、介质和电子设备
WO2022218438A1 (en) * 2021-04-16 2022-10-20 Shanghai United Imaging Healthcare Co., Ltd. Calibration methods and systems for imaging field

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5887074A (en) * 1996-12-13 1999-03-23 Siemens Corporate Research, Inc. Local principal component based method for detecting activation signals in functional MR images
EP2063392A1 (en) * 2007-11-24 2009-05-27 Barco NV Image processing of medical images
CN102366323A (zh) * 2011-09-30 2012-03-07 中国科学院自动化研究所 一种基于pca和gca的磁共振脑成像因果连接强度的检测方法
CN102749600A (zh) * 2012-05-30 2012-10-24 苏州安科医疗系统有限公司 一种磁共振多通道图像的合成方法
CN103185878A (zh) * 2011-12-27 2013-07-03 上海联影医疗科技有限公司 磁共振图像并行采集以及图像重建方法
CN103430038A (zh) * 2011-03-17 2013-12-04 皇家飞利浦有限公司 没有图像劣化的较快的逐个信道重建的mri方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7202663B2 (en) * 2003-11-12 2007-04-10 Iovivo Corporation Method for generating fast magnetic resonance images
US7511495B2 (en) * 2005-04-25 2009-03-31 University Of Utah Systems and methods for image reconstruction of sensitivity encoded MRI data
US20080278165A1 (en) * 2007-04-18 2008-11-13 Yu Li Method and apparatus for reconstruction of an image in image space using basis functions (RIB) for partially parallel imaging
DE102009014461B4 (de) * 2009-03-23 2011-06-22 Siemens Aktiengesellschaft, 80333 Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquistionstechnik

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5887074A (en) * 1996-12-13 1999-03-23 Siemens Corporate Research, Inc. Local principal component based method for detecting activation signals in functional MR images
EP2063392A1 (en) * 2007-11-24 2009-05-27 Barco NV Image processing of medical images
CN103430038A (zh) * 2011-03-17 2013-12-04 皇家飞利浦有限公司 没有图像劣化的较快的逐个信道重建的mri方法
CN102366323A (zh) * 2011-09-30 2012-03-07 中国科学院自动化研究所 一种基于pca和gca的磁共振脑成像因果连接强度的检测方法
CN103185878A (zh) * 2011-12-27 2013-07-03 上海联影医疗科技有限公司 磁共振图像并行采集以及图像重建方法
CN102749600A (zh) * 2012-05-30 2012-10-24 苏州安科医疗系统有限公司 一种磁共振多通道图像的合成方法

Also Published As

Publication number Publication date
CN104778341A (zh) 2015-07-15
CN107092805B (zh) 2020-08-04
CN107092805A (zh) 2017-08-25

Similar Documents

Publication Publication Date Title
CN104778341B (zh) 磁共振线圈合并系数计算方法、磁共振成像方法及其装置
CN110287932B (zh) 基于深度学习图像语义分割的道路阻断信息提取方法
CN110414377B (zh) 一种基于尺度注意力网络的遥感图像场景分类方法
Myllymäki et al. GET: Global envelopes in R
CN106796716B (zh) 用于为低分辨率图像提供超分辨率的设备和方法
KR101599219B1 (ko) 3차원 의료 영상에서 랜드마크들의 자동 등록을 위한 장치 및 방법
CN108647741A (zh) 一种基于迁移学习的图像分类方法和系统
KR20200032651A (ko) 3차원 영상 재구성 장치 및 그 방법
Zhang et al. LU-NET: An improved U-Net for ventricular segmentation
CN107622229A (zh) 一种基于融合特征的视频车辆重识别方法与系统
CN102651124B (zh) 基于冗余字典稀疏表示和评价指标的图像融合方法
CN106127263A (zh) 基于三维特征提取的人脑磁共振图像分类识别方法及系统
CN107679574A (zh) 超声图像处理方法及系统
CN106650572A (zh) 一种指纹图像的质量评估方法
CN104020430A (zh) 磁共振成像运动伪影的校正方法及系统
CN113283419B (zh) 基于注意力的卷积神经网络指针式仪表图像读数识别方法
CN104361318A (zh) 一种基于弥散张量成像技术的疾病诊断辅助系统及方法
US20210150284A1 (en) Classification model building apparatus and classification model building method thereof
CN109859189A (zh) 一种基于深度学习的年龄估计方法
CN102422324A (zh) 年龄估计设备、方法和程序
CN115439654B (zh) 动态约束下的弱监督农田地块精细化分割方法及系统
KR102604216B1 (ko) 전역 어텐션 기반 컨볼루션 네트워크를 이용한 이미지 분할 방법 및 장치
CN110853048A (zh) 基于粗、精训练的mri图像分割方法、装置和存储介质
CN104200441A (zh) 基于高阶奇异值分解的磁共振图像去噪方法
CN115131580A (zh) 基于注意力机制的空间目标小样本识别方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

CP01 Change in the name or title of a patent holder
CP02 Change in the address of a patent holder

Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee before: Shanghai Lianying Medical Technology Co.,Ltd.

CP02 Change in the address of a patent holder