CN107205686A - 去除在磁共振成像的重影伪影的方法及为此的磁共振装置 - Google Patents

去除在磁共振成像的重影伪影的方法及为此的磁共振装置 Download PDF

Info

Publication number
CN107205686A
CN107205686A CN201680002068.0A CN201680002068A CN107205686A CN 107205686 A CN107205686 A CN 107205686A CN 201680002068 A CN201680002068 A CN 201680002068A CN 107205686 A CN107205686 A CN 107205686A
Authority
CN
China
Prior art keywords
imaging
scanning line
odd
line imaging
magnetic resonance
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
Application number
CN201680002068.0A
Other languages
English (en)
Other versions
CN107205686B (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.)
Korea Advanced Institute of Science and Technology KAIST
Original Assignee
Korea Advanced Institute of Science and Technology KAIST
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 Korea Advanced Institute of Science and Technology KAIST filed Critical Korea Advanced Institute of Science and Technology KAIST
Publication of CN107205686A publication Critical patent/CN107205686A/zh
Application granted granted Critical
Publication of CN107205686B publication Critical patent/CN107205686B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • 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
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
    • 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
    • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)

Abstract

根据实施例，提供复原磁共振成像的方法。复原磁共振成像的方法，其步骤可包括：接收回波平面成像(echo planner imaging)数据；在回波平面成像数据的K‑空间数据，获得偶数次的扫描线成像和奇数次的扫描线成像；及利用偶数次的扫描线成像及奇数次的扫描线成像，复原被损失的部分。

Description

去除在磁共振成像的重影伪影的方法及为此的磁共振装置
技术领域
本发明涉及在磁共振成像(magnetic resonance imaging，MRI)装置适用回波平面成像(echo palnner imaging，EPI)方法，为了去除在获得的回波平面成像数据发生的重影伪影的方法。
背景技术
磁共振成像是在一定的磁场内，收集将排列的元素由无线电频率(RF)脉冲共振出的信号，获得人体内部成像的技术。作为以超高速获得磁共振成像的代表性方法，具有回波平面成像(echo planner imaging，EPI)技术。回波平面成像技术在40至100ms时间期间，可接收单一的自由感应衰减(FID；Free Induction Decay)信号中的磁共振(MagneticResonance)图像，可实现超高速成像技术。由回波平面成像的快速扫描技术，可减少因目标移动发生的缺陷，或磁共振成像(MRI；Magnetic Resonance Imaging)的问题。但是，回波平面成像技术生成大领域梯度，大概由1kHz的比率进行交换，所以，从硬件上不容易体现图像，并且得到目标的磁化率(magnetic susceptibility)的影响，容易发生图像的严重歪曲。
与由一定方向读取数据的其他方法不同，在回波平面成像技术，因读取k-空间偶数次的扫描线和奇数次的扫描线数据的方向不同，会发生各线间的位相差。因这些位相差在成像域生成重影伪影(ghost artifact)，且这妨碍回波平面成像的正确地分析，所以，用于去除其的方法被很多地研究。
代表性地具有，得到被复原的成像之前得到参考成像(reference scan)之后，在k-空间计算偶数次的扫描线和奇数次的扫描线间的位相差，补正对象成像位相差的方法。但是，这需要附加的成像，需要更长的时间。因此，使用这种方法时，为了弥补时间上的缺点，在附加成像上大多利用只得到3个扫描线进行补正的方法。但是，在这种情况也是受到磁共振成像装置内的线圈绕度的影响，具有不可完全地去除重影伪影的缺点。
此外，也有修改用于获得回波平面成像的序列，得到可抵消重影伪影(ghostartifact)的附加成像的方法，但是，这也需要得到本成像之外的附加成像，所以，具有需要附加时间的缺点。
发明内容
技术课题
本发明的目的是当复原磁共振成像时，不获得基准成像(reference scan)也能去除在磁共振成像装置的回波平面成像生成的重影伪影(ghost artifact)。
技术方案
根据实施例，提供一种复原磁共振成像的方法。复原磁共振成像的方法，其步骤可包括：接收回波平面成像(echo planner imaging)数据；在回波平面成像数据的K-空间数据，获得偶数次的扫描线成像和奇数次的扫描线成像；及利用偶数次的扫描线成像及奇数次的扫描线成像，复原被损失的部分。
在这种情况下，利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤，可包括以消灭滤波器存在于k-空间为基础，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，通过低链接矩阵的完成化算法，插入未被测定的傅里叶数据，来复原被损失的部分。
将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，复原被损失的部分的步骤，可包括在频率域利用所述消灭滤波器和k-空间的卷积演算的交换法则，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，变换为分块汉克尔矩阵的形态。
其中，复原磁共振成像的方法，还可包括结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，生成一个被复原的成像。
生成一个被复原的成像的步骤，可包括利用被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的平方和(sum of squares)，生成一个被复原的成像。
利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤，可包括将具有金字塔结构的各金字塔补片变更为分块汉克尔矩阵之后，利用多个最小链接矩阵完成化算法，复原偶数次的扫描线成像和奇数次的扫描线成像。
利用金字塔补片复原成像的步骤，通过矩阵完成化算法，复原偶数次的扫描线成像和奇数次的扫描线成像，且缩小金字塔补片大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，通过所述矩阵完成化算法，可复原k-空间的低频领域。
利用金字塔补片复原成像的步骤，在具有金字塔结构的各金字塔补片中，对于噪音包括在预先设定的基准以上的金字塔补片，适用预先设定的基准以下的停止基准，且在具有金字塔结构的各金字塔补片中，对于只切除低频率领域的金字塔补片，可适用预先设定的基准以上的停止基准。
其中，回波平面成像数据可包括单一线圈回波平面成像数据及多重线圈回波平面成像数据中至少一个。
利用偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤，可包括在所述偶数次的扫描线成像和所述奇数次的扫描线成像的k-空间，分别赋予加重值(weighting)构成分块汉克尔矩阵，复原被损失的部分。
根据其他实施例，可提供磁共振成像复原装置。磁共振成像复原装置可包括：处理单元，在回波平面成像数据的K-空间数据，获得偶数次的扫描线成像和奇数次的扫描线成像，且利用偶数次的扫描线成像及奇数次的扫描线成像，复原被损失(missing)的部分。
在这种情况下，处理单元以消灭滤波器存在于k-空间为基础，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，通过低链接矩阵的完成化算法，可插入未被测定的傅里叶数据，来复原被损失的部分。
另外，处理单元在频率域利用消灭滤波器和k-空间的卷积演算的交换法则，可将偶数次的扫描线成像和奇数次的扫描线成像，变换为分块汉克尔矩阵的形态。
处理单元结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，可生成一个被复原的成像。
处理单元利用被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的平方和(sum of squares)，可生成一个被复原的成像。
处理单元将具有金字塔结构的各金字塔补片变更为分块汉克尔矩阵之后，利用多个最小链接矩阵完成化算法，可复原偶数次的扫描线成像和奇数次的扫描线成像。
处理单元通过矩阵完成化算法复原偶数次的扫描线成像和奇数次的扫描线成像，且缩小所述金字塔补片大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，可通过矩阵完成化算法，复原k-空间的低频领域。
处理单元在具有金字塔结构的各金字塔补片中，对于噪音包括在预先设定的基准以上的金字塔补片，适用预先设定的基准以下的停止基准，且在具有金字塔结构的各金字塔补片中，对于只切除低频率领域的金字塔补片，可适用预先设定的基准以上的停止基准。
处理单元在偶数次的扫描线成像和奇数次的扫描线成像的k-空间，可分别赋予加重值(weighting)构成分块汉克尔矩阵，复原被损失的部分。
技术效果
根据实施例，在回波平面成像(EPI)不获得基准成像，也可去除磁共振成像的重影伪影。
根据实施例，不利用多重线圈，只利用一个线圈也可去除重影伪影。
根据实施例，重影伪影有效地被去除，在功能性磁共振成像(fMRI)分析，可得到更正确的结果。
附图说明
图1是示出根据实施例的磁共振成像复原装置的框图。
图2是示出根据实施例，复原磁共振成像方法的流程图。
图3是示出根据实施例，利用k-空间复原磁共振成像的方法。
图4是示出根据实施例，在信息不足的k-空间构成分块汉克尔矩阵的方法。
图5是示出根据实施例，分别利用偶数次的扫描线成像和奇数次的扫描线成像，构成分块汉克尔矩阵的结果。
图6是示出根据实施例，根据金字塔结构的各等级(level)的复原补片大小的变化。
图7是示出根据实施例，适用金字塔结构的复原过程的流程图。
图8A至8B是示出根据实施例，按照复原磁共振成像方法的结果。
具体实施方式
以下，参考附图对本发明的实施例进行详细地说明。
图1是示出根据实施例的磁共振成像复原装置的框图。
参考图1，磁共振成像复原装置100可包括通信单元110、处理单元120及存储单元130。
根据实施例，磁共振成像复原装置100可在回波平面成像(EPI)补正重影伪影(ghost artifact)。其中，重影伪影的去除利用消灭滤波器存在于K-空间，将由偶数次的扫描线成像和奇数次的扫描线成像分离的K-空间数据，构成为分块汉克尔矩阵，可变更成链接不足的问题。在这种情况下，磁共振成像复原装置100可通过多样的低链接矩阵完成化算法，解决链接不足的问题。
根据实施例，通信单元110可从磁共振成像摄影装置(未示出)，接收回波平面成像(echo planner imaging)数据。其中，通行单元也可从磁共振成像摄影装置，以无线接收回波平面成像数据，且也可由有线接收回拨平面成像数据。
根据一个侧面，回波平面成像数据可包括单一线圈回波平面成像数据及多重线圈回波平面成像数据中至少一个。
根据实施例，通信单元110可将被复原的磁共振成像输出在用户可确认的输出装置。
根据实施例，处理单元120在接收的回波平面成像数据，可获得在K-空间数据的偶数次的扫描线成像和奇数次的扫描线成像。
在这种情况下，处理单元120可利用偶数次的扫描线成像和奇数次的扫描线成像，复原被损失的部分。
处理单元120以消灭滤波器存在于k-空间为基础，将偶数次的扫描线成像和奇数次的扫描线成像，构成各个的分块汉克尔矩阵，变更为链接不足的问题，插值未被测定的傅里叶数据，可复原被损失的部分。其中，处理单元120的链接不足的问题，可通过多个低链接矩阵的完成化算法复原。
处理单元120在频率域利用消灭滤波器和k-空间的卷积演算的交换法则，可将所述偶数次的扫描线成像和所述奇数次的扫描线成像，分别变换成分块汉克尔矩阵的形态。
根据实施例，处理单元120结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，可生成一个被复原的成像。在这种情况下，处理单元120利用被复原的偶数次的扫描线成像和所述被复原的奇数次的扫描线成像的平方和(sum of squares)，可生成一个被复原的成像。
根据实施例，处理单元120将具有金字塔结构的各金字塔补片，变更为分块汉克尔矩阵之后，利用多个最小链接矩阵完成化算法，可复原偶数次的扫描线成像和奇数次的扫描线成像。
处理单元120通过矩阵完成化算法，复原偶数次的扫描线成像和奇数次的扫描线成像，缩小金字塔补片的大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，通过阵列完成化算法，可复原k-空间的低频率领域。
在这种情况下，处理单元120在具有金字塔结构的各金字塔补片中，对于噪音包括在预先设定的基准以上的金字塔补片，适用预先设定的基准以下的停止基准，且在具有金字塔结构的各金字塔补片中，对于只切除低频率领域的金字塔补片，可适用预先设定的基准以上的停止基准
处理单元120在偶数次的扫描线成像和奇数次的扫描线成像的k-空间，分别赋予加重值(weighting)构成分块汉克尔矩阵，可复原被损失的部分。其中，加重值可以是根据回波平面成像的稀疏变换(sparsifying transform)特征的加重值。因此，处理单元120在构成分块汉克尔矩阵之前，在k-空间数据赋予加重值，可将回波平面成像稀疏(sparse)地变换。
根据实施例，存储单元130可临时存储从通信单元110接收的回波平面成像，也可存储从处理单元120复原的成像。
图2是示出根据实施例，复原磁共振成像方法的流程图。
在步骤210中，磁共振成像复原装置可接收回波平面成像。
在步骤220中，磁共振成像复原装置在接收的回波平面成像的K-空间，可分离偶数次的扫描线成像和奇数次的扫描线成像。
在步骤230中，磁共振成像复原装置利用偶数次的扫描线成像和奇数次的扫描线成像，可复原被损失的部分。
在这种情况下，磁共振成像复原装置以消灭滤波器存在于k-空间为基础，将所述偶数次的扫描线成像和所述奇数次的扫描线成像构成分块汉克尔矩阵，可变更为链接不足的问题进行复原。其中，链接不足的问题可通过多个低链接矩阵完成化算法解决。根据实施例，磁共振成像复原装置将偶数次的扫描线成像和奇数次的扫描线成像构成为分块汉克尔矩阵，通过低链接矩阵的完成化算法，插值未被测定的傅里叶数据，可复原被损失的部分。
根据实施例，磁共振成像复原装置在频率域，利用消灭滤波器和k-空间的卷积演算的交换法则，可将偶数次的扫描线成像和奇数次的扫描线成像，变换成分块汉克尔矩阵的形态。
根据实施例，磁共振成像复原装置在偶数次的扫描线成像和奇数次的扫描线成像的k-空间，分别赋予加重值(weighting)构成分块汉克尔矩阵，可复原被损失的部分。其中，加重值可以是根据回波平面成像的稀疏变换(sparsifying transform)特征的加重值。因此，磁共振成像复原装置在构成分块汉克尔矩阵之前，可稀疏(sparse)回波平面成像。
磁共振成像复原装置将具有金字塔结构的各金字塔补片变更为分块汉克尔矩阵之后，利用多个最小链接矩阵完成化算法，可分别复原所述偶数次的扫描线成像和所述奇数次的扫描线成像。在这种情况下，磁共振成像复原装置通过矩阵完成化算法，分别复原偶数次的扫描线成像和奇数次的扫描线成像，且缩小所述金字塔补片的大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像值作为初始值，通过所述矩阵完成化算法，可复原所述k-空间的低频率领域。
根据实施例，磁共振复原装置结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，可生成一个成像。
根据实施例，磁共振成像复原装置利用被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的平方和(sum of squares)，可生成一个被复原的成像。
图3是示出根据实施例，利用k-空间复原磁共振成像的方法。
参考图3，根据实施例的磁共振成像复原装置可将回波平面成像311的k-空间312，分为偶数次的扫描线成像321，322和奇数次的扫描线成像331，332。
根据实施例，磁共振成像复原装置复原偶数次的扫描线成像321，322，可获得被复原的偶数次的扫描线成像341，342，且复原奇数次的扫描线成像331，332，可获得被复原的奇数次的扫描线成像351，352。
根据实施例，磁共振成像复原装置结合被复原的偶数次的扫描线成像341，342和被复原的奇数次的扫描线成像351，352，可获得重影伪影(ghost artifact)被去除的最终成像360。
利用数学式说明一系列过程如下。
回波平面成像可显示为如下的数学式1。
【数学式1】
其中，m是成像的强度(intensity)，Δf是由磁场不均匀(field inhomogeneity)生成的频率偏差(frequency offset)。x和y分别显示判读输出(read-out)、位相编码(phase encoding)的方向。TE是回波时间，ESP显示各回波之间的时间。N显示总回波个数，n显示各回波线的指数。
在数学式1的样式中，包括由回波线的位相差的重影伪影(ghost artifact)。即，随着n是偶数还是奇数，(﹣1)n的值变化，所以，偶数次的扫描线和奇数次的扫描线具有不同的极性(polarity)。将此作为虚拟的k-空间时，可显示如下的数学式2。
【数学式2】
其中，Sn，+(kx，ky)是只聚集偶数次的扫描线的k-空间，且Sn，-(kx，ky)是只聚集奇数次的扫描线的k-空间。因此，如此被区分的两个k-空间，因为各偶数次的扫描线和奇数次的扫描线全部空着，所以，可视为两倍的降采样。因此，复原磁共振成像的方法可被变换成填充k-空间空部分的问题。
去除数学式2的两个式，可显示如下的数学式3。
【数学式3】
其中，
在数学式3中，频率偏差(frequency offset)充分地小时，可接近如下的数学式4。
【数学式4】
因此，上述式如同在无线电频率偏差(RF offset)和图像相乘值的梯度(gradient)进行傅里叶变换(fourier transform)。经A(x，y)Δf(x，y)光滑(smooth)的性质，此值的密度稀疏(sparse)，所以，利用奇数次的线k空间和利用偶数次的线k空间的差异(Sn，Δ)稀疏(sparse)，且利用此构成的汉克尔矩阵具有低链接。因此，存在满足数学式5的消灭滤波器
【数学式5】
其中，是Sn，Δ的汉克尔矩阵，是2D消灭滤波器的向量形态。在这种情况下，Sn，Δ＝Sn，--Sn，+，所以，可获得如同数学式6的关系。
【数学式6】
其中，
因此，从只由偶数次的扫描线形成的k-空间和只由奇数次的扫描线形成的k-空间，分别构成分块汉克尔矩阵，连接两个矩阵可形成y矩阵，且为了填充空的数据部分，可利用y的低链接特性(low rankness)。
分别填充Sn，+和Sn，-的空着部分复原之后，结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，可生成没有重影伪影(ghost artifact)的最终成像。在这种情况下，被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的结合，可利用平方和(sum of square)方法。
根据实施例，利用多重线圈时，也构成如下数学式6的汉克尔矩阵，由相同的方式可复原。
【数学式7】
y(P)＝[y(1)…y(r)]
其中，y(i)意味着第i个线圈的汉克尔矩阵。
图4是示出根据实施例，在信息不足的k-空间构成分块汉克尔矩阵的方法。
根据实施例，在回波平面成像(EPI)的k空间410，可抽取各补片411、412、413、414、415。其中，在各补片可确认奇数次的扫描线被填满，且偶数次的扫描线空着。
在这种情况下，磁共振成像复原装置对各补片，可构成汉克尔矩阵。即，磁共振成像复原装置可构成对补片411的汉克尔矩阵421、对补片412的汉克尔矩阵422，424、对补片413的汉克尔矩阵423，425，427、对补片414的汉克尔矩阵426，428、对补片415的汉克尔矩阵429。通过此利用汉克尔矩阵的链接不足特性，可复原空部分的信息。
图5是示出根据实施例，分别利用偶数次的扫描线成像和奇数次的扫描线成像，构成分块汉克尔矩阵的结果。
参考图5，根据实施例回波平面成像数据511可由汉克尔矩阵521构成，且回波平面成像数据512可由汉克尔矩阵522构成。
在这种情况下，磁共振成像复原装置分别填充在各汉克尔矩阵521，522的空部分复原之后，通过平方和(sum of square)可生成没有重影伪影(ghost artifact)的最终成像。
图6是示出根据实施例，根据金字塔结构的各等级(level)的复原补片大小的变化。
参考图6，低链接完成化算法可通过金字塔结构化方式有效地被计算。低链接完成化算法可由独立性的补片被计算。使用如图6的金字塔结构时，可按阶段减小补片大小，且按阶段考虑信号对噪音比，也可不同地进行复原，这可减少反复性复原方式的低链接矩阵完成化算法的相互作用的次数。这些金字塔结构可被适用在kx-ky域。
根据实施例，可确认在第一步骤的补片610和第二步骤的补片620。在这种情况下，根据一个实施例，步骤越增加，可将补片x轴长度减为之前步骤的一般。其中，各补片的中线单元611一直包括低频部分，所以，随着步骤增加可进行没有噪音的更正确的复原。这具有减少低链接矩阵完成化算法的反复次数的效果，且各步骤的初始值使用之前步骤的复原数据，所以，可获得更正确的结果值。
根据其他实施例，补片的大小比起之前步骤，可将横着和竖着长度减少为更小。
图7是示出根据实施例，适用金字塔结构的复原过程的流程图。
在步骤710中，磁共振成像复原装置可在回波平面成像的k-空间，分离偶数次的扫描线成像和奇数次的扫描线成像。
在步骤720中，磁共振成像复原装置可决定金字塔补片。在这种情况下，磁共振成像复原装置通过矩阵完成化算法，复原偶数次的扫描线成像和奇数次的扫描线成像，可缩小金字塔补片的大小。
在步骤730中，磁共振成像复原装置可对偶数次的扫描线成像和奇数次的扫描线成像，构成k-空间的汉克尔矩阵。
在步骤740中，磁共振成像复原装置可通过低链接矩阵完成算法，再构成成像。在这种情况下，磁共振成像复原装置将以被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，可通过矩阵完成化算法，复原k-空间的低频领域。
在步骤750中，磁共振成像复原装置可判断金字塔等级是否达到所愿的等级。
在这种情况下，根据实施例，比起经过步骤750之前的补片，将补片大小减少1/2，直到达到所愿的金字塔等级为止，可反复步骤720至步骤740。
在适用金字塔结构的复原过程中，根据各等级按步骤考虑信号对噪音比，可不同地进行复原，且可减少反复性复原方式的低链接矩阵完成化算法的相互作用次数。
图8A至8B是示出根据实施例，按照复原磁共振成像方法的结果。
参考图8A，可确认根据已有的方式，示出的包括重影伪影的磁共振成像的侧面811、正面812、上面813，和根据实施例去除重影伪影的磁共振成像的侧面821、正面822、上面823。
在这种情况下，参考图8B，在根据已有方式的磁共振成像831，可确认重影伪影(箭头)被包括，且在根据实施例的磁共振成像，可确认在单一线圈图像832和多重线圈图形833都去除重影伪影的结果。
上述说明的装置可由硬件构成要素、软件构成要素、和/或硬件构成要素及软件构成要素的组合被体现。例如，说明的装置及构成要素，可利用类似处理器、控制器、算术逻辑单元ALU(arithmetic logic unit)、数字信号处理器(digital signal processor)、微型计算机、现场可编程阵列FPA(field programmable array)、可编程逻辑单元PLU(programmable logic unit)、微处理器、或实行指令(instruction)的其他任何装置、一个以上的范用计算机或特殊目的计算机被体现。处理装置可实行操作系统(OS)及该操作系统中所实行的一个以上的软件应用程序。此外，处理装置可应答软件的实行，来存取、存储、操作、处理、生成数据。为了便于理解，处理装置被说明是使用一个，但在相关技术领域中，具有通常知识的技术人员应理解，处理装置可包括多个处理元件(processing element)和/或多个类型的处理元件。例如，处理装置可包括多个处理器或一个处理器，以及一个控制器。此外，也可以是类似并行处理器(parallel processor)的其他处理配置(processingconfiguration)。
软件可包括计算机程序(computer program)，码(code)，命令(instruction)，或者这些中一个以上的组合，为了如愿的操作，可命令处理装置的构成，或单独的或者结合的(collectively)处理装置。软件及/或数据可以是，按处理装置分析或为了在处理装置提供命令或者数据、某些类型的机械、构成要素(component)、物理装置、虚拟装置(virtualequipment)、计算机存储媒体或装置，或者被传输的信号波(signal wave)永久的或者暂时的体现(embody)。软件分散在连接网络的计算机系统上，可以用分散的方法存储或实行。软件及数据可存储在一个以上的可分析计算机的记录媒体。
根据实施例的方法可通过多种计算机手段以可实行的程序指令形态被记录在计算机可读媒体中。计算机可读媒体可包括独立的或结合的程序指令、数据文件、数据结构等。媒体和程序指令可为了本发明被专门设计和创建，或为计算机软件技术人员熟知而应用。计算机可读媒体的例子包括：磁媒体(magnetic media)，如硬盘、软盘和磁带；光学媒体(optical media)，如CD-ROM、DVD；磁光媒体(magneto-optical media)，如光盘(flopticaldisk)；和专门配置为存储和实行程序指令的硬件装置，如只读存储器(ROM)、随机存取存储器(RAM)等。程序指令的例子，既包括由编译器产生的机器代码，也包括使用解释程序并可通过计算机被实行的高级语言代码。为实行实施例的运作，所述硬件装置可被配置为以一个以上的软件模来运作，反之亦然。
如上所示，本发明虽然由限定的实施例和附图进行了说明，但是，在本发明所属领域中具备通常知识的人均可以从此记载中进行各种修改和变形。例如，可通过与说明的方法不同的顺序来实行所说明的技术，和/或是通过与说明的方法不同的形态来结合或组合所说明的系统、结构、装置、电路等的构成要素，或是通过其他构成要素或同等事物来代替或置换也可获得适当结果。所以，其他实现、其他实施例及与专利请求范围均等的，也属于后述的专利申请范围的范围。

Claims (20)

1.一种复原磁共振成像的方法，其步骤包括：
接收回波平面成像数据；
在所述回波平面成像数据的K-空间数据，获得偶数次的扫描线成像和奇数次的扫描线成像；及
利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分。
2.根据权利要求1所述的复原磁共振成像的方法，其特征为，
利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤包括：
以消灭滤波器存在于k-空间为基础，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，通过低链接矩阵的完成化算法，插入未被测定的傅里叶数据，来复原被损失的部分。
3.根据权利要求2所述的复原磁共振成像的方法，其特征为，
将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，复原被损失的部分的步骤包括：
在频率域利用所述消灭滤波器和k-空间的卷积演算的交换法则，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，变换为分块汉克尔矩阵的形态。
4.根据权利要求1所述的复原磁共振成像的方法，还包括：
结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，生成一个被复原的成像。
5.根据权利要求4所述的复原磁共振成像的方法，其特征为，
所述生成一个被复原的成像的步骤包括：
利用所述被复原的偶数次的扫描线成像和所述被复原的奇数次的扫描线成像的平方和，生成一个被复原的成像。
6.根据权利要求1所述的复原磁共振成像的方法，其特征为，
利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤包括：
将具有金字塔结构的各金字塔补片变更为分块汉克尔矩阵之后，利用多个最少链接矩阵完成化算法，复原所述偶数次的扫描线成像和所述奇数次的扫描线成像。
7.根据权利要求6所述的复原磁共振成像的方法，其特征为，
所述利用金字塔补片复原成像的步骤，通过矩阵完成化算法，复原所述偶数次的扫描线成像和所述奇数次的扫描线成像，且缩小所述金字塔补片大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，通过所述矩阵完成化算法，复原所述k-空间的低频领域。
8.根据权利要求6所述的复原磁共振成像的方法，其特征为，
所述利用金字塔补片复原成像的步骤，在具有所述金字塔结构的各金字塔补片中，对于噪音包括在预先设定的基准以上的金字塔补片，适用预先设定的基准以下的停止基准，且在具有所述金字塔结构的各金字塔补片中，对于只切除低频率领域的金字塔补片，适用预先设定的基准以上的停止基准。
9.根据权利要求1所述的复原磁共振成像的方法，其特征为，
所述回波平面成像数据包括单一线圈回波平面成像数据及多重线圈回波平面成像数据中至少一个。
10.根据权利要求1所述的复原磁共振成像的方法，其特征为，
利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分的步骤包括：
在所述偶数次的扫描线成像和所述奇数次的扫描线成像的k-空间，分别赋予加重值构成分块汉克尔矩阵，复原被损失的部分。
11.一种计算机程序，用于执行权利要求1方法，与硬件结合并存储在媒体。
12.一种复原磁共振成像的装置，其包括：
处理单元，在回波平面成像数据的K-空间数据，获得偶数次的扫描线成像和奇数次的扫描线成像，且利用所述偶数次的扫描线成像及所述奇数次的扫描线成像，复原被损失的部分。
13.根据权利要求12所述的复原磁共振成像的装置，其特征为，
所述处理单元以消灭滤波器存在于k-空间为基础，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，构成为分块汉克尔矩阵，通过低链接矩阵的完成化算法，插入未被测定的傅里叶数据，来复原被损失的部分。
14.根据权利要求13所述的复原磁共振成像的装置，其特征为，
所述处理单元在频率域利用所述消灭滤波器和k-空间的卷积演算的交换法则，将所述偶数次的扫描线成像和所述奇数次的扫描线成像，变换为分块汉克尔矩阵的形态。
15.根据权利要求12所述的复原磁共振成像的装置，其特征为，
所述处理单元结合被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像，生成一个被复原的成像。
16.根据权利要求15所述的复原磁共振成像的装置，其特征为，
所述处理单元利用所述被复原的偶数次的扫描线成像和所述被复原的奇数次的扫描线成像的平方和，生成一个被复原的成像。
17.根据权利要求12所述的复原磁共振成像的装置，其特征为，
将具有金字塔结构的各金字塔补片变更为分块汉克尔矩阵之后，利用多个最少链接矩阵完成化算法，复原所述偶数次的扫描线成像和所述奇数次的扫描线成像。
18.根据权利要求17所述的复原磁共振成像的装置，其特征为，
所述处理单元通过矩阵完成化算法复原所述偶数次的扫描线成像和所述奇数次的扫描线成像，且缩小所述金字塔补片大小之后，将被复原的偶数次的扫描线成像和被复原的奇数次的扫描线成像的值作为初始值，通过所述矩阵完成化算法，复原所述k-空间的低频领域。
19.根据权利要求12所述的复原磁共振成像的装置，其特征为，
所述处理单元在具有所述金字塔结构的各金字塔补片中，对于噪音包括在预先设定的基准以上的金字塔补片，适用预先设定的基准以下的停止基准，且在具有所述金字塔结构的各金字塔补片中，对于只切除低频率领域的金字塔补片，适用预先设定的基准以上的停止基准。
20.根据权利要求12所述的复原磁共振成像的装置，其特征为，
所述处理单元在所述偶数次的扫描线成像和所述奇数次的扫描线成像的k-空间，分别赋予加重值构成分块汉克尔矩阵，复原被损失的部分。
CN201680002068.0A 2015-10-20 2016-02-05 去除在磁共振成像的重影伪影的方法及为此的磁共振装置 Active CN107205686B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR1020150146208A KR101697872B1 (ko) 2015-10-20 2015-10-20 자기 공명 영상에서의 고스트 아티팩트를 제거하는 방법 및 이를 위한 자기 공명 장치
KR10-2015-0146208 2015-10-20
PCT/KR2016/001279 WO2017069347A1 (ko) 2015-10-20 2016-02-05 자기 공명 영상에서의 고스트 아티팩트를 제거하는 방법 및 이를 위한 자기 공명 장치

Publications (2)

Publication Number Publication Date
CN107205686A true CN107205686A (zh) 2017-09-26
CN107205686B CN107205686B (zh) 2020-08-07

Family

ID=57992085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680002068.0A Active CN107205686B (zh) 2015-10-20 2016-02-05 去除在磁共振成像的重影伪影的方法及为此的磁共振装置

Country Status (5)

Country Link
US (1) US10627469B2 (zh)
EP (1) EP3366209B1 (zh)
KR (1) KR101697872B1 (zh)
CN (1) CN107205686B (zh)
WO (1) WO2017069347A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581253A (zh) * 2017-12-26 2019-04-05 上海联影医疗科技有限公司 用于磁共振成像的方法和系统
CN111580032A (zh) * 2019-02-15 2020-08-25 通用电气精准医疗有限责任公司 用于去除磁共振成像中峰值噪声的方法和系统

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102131505B1 (ko) * 2018-08-29 2020-07-07 한국과학기술원 뉴럴 네트워크를 이용한 에코평면 영상 고스트 아티팩트 제거 방법 및 그 장치
KR102066009B1 (ko) * 2018-09-14 2020-01-14 가천대학교 산학협력단 순환신경망을 사용한 의료영상기기의 영상복원 시스템
US20220236358A1 (en) * 2019-04-26 2022-07-28 Mayo Foundation For Medical Education And Research Model-Based Iterative Reconstruction for Magnetic Resonance Imaging with Echo Planar Readout
WO2021167124A1 (ko) * 2020-02-19 2021-08-26 서울대학교 산학협력단 Mri 데이터 및 추가 정보를 이용하는 변환 네트워크부를 포함하는 mri 데이터 변환장치 및 이를 이용한 mri 데이터 변환방법
KR102472546B1 (ko) * 2021-08-12 2022-11-30 주식회사 에어스메디컬 노이즈 투 노이즈 기법이 적용되는 자기 공명 영상 처리 장치 및 그 방법

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5910728A (en) * 1996-11-12 1999-06-08 Beth Israel Deaconess Medical Center Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays
US6697507B1 (en) * 1997-12-10 2004-02-24 Btg International Limited Ghost artifact reduction
CN1616985A (zh) * 2003-11-14 2005-05-18 西门子公司 基于常规ppa再现方法改进的磁共振断层造影成像设备
US20090085563A1 (en) * 2007-09-27 2009-04-02 Yoshitaka Bito Magnetic resonance imaging system
US20150091570A1 (en) * 2013-10-02 2015-04-02 Siemens Aktiengesellschaft Method and magnetic resonance apparatus for image reconstruction with multiple virtual coils

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008154887A (ja) * 2006-12-26 2008-07-10 Ge Medical Systems Global Technology Co Llc Ｍｒｉ装置
KR20120138478A (ko) * 2011-06-15 2012-12-26 삼성메디슨 주식회사 합성 구경을 이용한 초음파 영상 장치 및 방법
KR20140122991A (ko) * 2013-04-11 2014-10-21 삼성전자주식회사 의료 영상 장치 및 그 제어 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5910728A (en) * 1996-11-12 1999-06-08 Beth Israel Deaconess Medical Center Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays
US6697507B1 (en) * 1997-12-10 2004-02-24 Btg International Limited Ghost artifact reduction
CN1616985A (zh) * 2003-11-14 2005-05-18 西门子公司 基于常规ppa再现方法改进的磁共振断层造影成像设备
US20090085563A1 (en) * 2007-09-27 2009-04-02 Yoshitaka Bito Magnetic resonance imaging system
US20150091570A1 (en) * 2013-10-02 2015-04-02 Siemens Aktiengesellschaft Method and magnetic resonance apparatus for image reconstruction with multiple virtual coils

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581253A (zh) * 2017-12-26 2019-04-05 上海联影医疗科技有限公司 用于磁共振成像的方法和系统
CN109581253B (zh) * 2017-12-26 2021-05-18 上海联影医疗科技股份有限公司 用于磁共振成像的方法和系统
CN111580032A (zh) * 2019-02-15 2020-08-25 通用电气精准医疗有限责任公司 用于去除磁共振成像中峰值噪声的方法和系统

Also Published As

Publication number Publication date
KR101697872B1 (ko) 2017-01-18
EP3366209B1 (en) 2022-01-12
US20170299680A1 (en) 2017-10-19
EP3366209A4 (en) 2018-11-07
EP3366209A1 (en) 2018-08-29
US10627469B2 (en) 2020-04-21
CN107205686B (zh) 2020-08-07
WO2017069347A1 (ko) 2017-04-27

Similar Documents

Publication Publication Date Title
CN107205686A (zh) 去除在磁共振成像的重影伪影的方法及为此的磁共振装置
Ran et al. MD-Recon-Net: a parallel dual-domain convolutional neural network for compressed sensing MRI
Liu et al. MANTIS: model‐augmented neural network with incoherent k‐space sampling for efficient MR parameter mapping
WO2018223275A1 (zh) 基于深度卷积网的一维部分傅里叶并行磁共振成像方法
CN104661592B (zh) 磁共振成像装置以及图像处理装置
CN102043137B (zh) 全程放射grappa校准
CN108375746A (zh) 一种相位反卷折方法与设备
CN110133556B (zh) 一种磁共振图像处理方法、装置、设备及存储介质
CN106796716A (zh) 用于为低分辨率图像提供超分辨率的设备和方法
CN110211050A (zh) 用于稀疏图像重建的系统和方法
US20200341100A1 (en) Mri system and method using neural network for detection of patient motion
CN107576925B (zh) 磁共振多对比度图像重建方法和装置
CN105957117B (zh) 并行磁共振的图像重建方法、装置及并行磁共振成像系统
WO2020114329A1 (zh) 磁共振快速参数成像方法及装置
CN103027681B (zh) 用于重构并行获取的mri图像的系统
CN104749538A (zh) 一种并行磁共振成像相位处理方法
CN111353947A (zh) 磁共振并行成像方法及相关设备
Zhou et al. Parallel imaging and convolutional neural network combined fast MR image reconstruction: Applications in low‐latency accelerated real‐time imaging
CN103218795B (zh) 基于自适应双字典学习的部分k空间序列图像重构方法
CN108776318A (zh) 磁共振多对比度图像重建
CN105051564A (zh) 用于压缩感测MRI中的经改进的k-空间采样的方法
CN103584864A (zh) 一种磁共振成像方法和装置
WO2023045482A1 (zh) 基于长距离注意力模型重建的多层磁共振成像方法及装置
CN1327810C (zh) 磁共振图像重建时所产生Gibbs环状伪影的消除方法
CN105334479A (zh) 一种磁共振成像方法和装置

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