CN103505210A - 一种实现水脂分离的磁共振成像方法和装置 - Google Patents

一种实现水脂分离的磁共振成像方法和装置 Download PDF

Info

Publication number
CN103505210A
CN103505210A CN201210218752.6A CN201210218752A CN103505210A CN 103505210 A CN103505210 A CN 103505210A CN 201210218752 A CN201210218752 A CN 201210218752A CN 103505210 A CN103505210 A CN 103505210A
Authority
CN
China
Prior art keywords
image
same phase
data
initial data
phase image
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
CN201210218752.6A
Other languages
English (en)
Other versions
CN103505210B (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.)
Siemens Shenzhen Magnetic Resonance Ltd
Original Assignee
Siemens Shenzhen Magnetic Resonance 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 Siemens Shenzhen Magnetic Resonance Ltd filed Critical Siemens Shenzhen Magnetic Resonance Ltd
Priority to CN201210218752.6A priority Critical patent/CN103505210B/zh
Priority to US13/927,359 priority patent/US9390476B2/en
Priority to KR1020130074429A priority patent/KR101644677B1/ko
Priority to JP2013135001A priority patent/JP6245861B2/ja
Publication of CN103505210A publication Critical patent/CN103505210A/zh
Application granted granted Critical
Publication of CN103505210B publication Critical patent/CN103505210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4818MR 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
    • G01R33/4824MR 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 using a non-Cartesian trajectory
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5615Echo 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/5617Echo 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 RF refocusing, e.g. RARE

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例提供了一种实现水脂分离的磁共振成像方法,包括:利用快速自旋回波刀锋伪影校正序列(TSE BLADE)采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;根据所述同相位图像的原始数据重建同相位图像,并根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;根据重建后的同相位图像和反相位图像计算水和脂肪的图像。由于本发明采用TSE BLADE序列采集k空间数据,从而继承了BLADE序列对刚体运动和脉动不敏感的优点,降低了对运动伪影的敏感程度,并且提高了图像的信噪比。

Description

一种实现水脂分离的磁共振成像方法和装置
技术领域
本发明涉及磁共振成像技术领域,特别涉及一种实现水脂分离的磁共振成像方法和装置。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)是利用磁共振现象进行成像的一种技术。磁共振现象的原理主要包括:包含单数质子的原子核,例如人体内广泛存在的氢原子核,其质子具有自旋运动,犹如一个小磁体,并且这些小磁体的自旋轴没有一定的规律,如果施加外在磁场,这些小磁体将按外在磁场的磁力线重新排列,具体为在平行于或反平行于外在磁场磁力线的两个方向排列,将上述平行于外在磁场磁力线的方向称为正纵向轴,将上述反平行于外在磁场磁力线的方向称为负纵向轴;原子核只具有纵向磁化分量,该纵向磁化分量既具有方向又具有幅度。用特定频率的射频(Radio Frequency,RF)脉冲激发处于外在磁场中的原子核,使这些原子核的自旋轴偏离正纵向轴或负纵向轴,产生共振,这就是磁共振现象。上述被激发的原子核的自旋轴偏离正纵向轴或负纵向轴之后,该原子核就具有了横向磁化分量。
停止发射射频脉冲后,被激发的原子核发射回波信号,将吸收的能量逐步以电磁波的形式释放出来,其相位和能级都恢复到激发前的状态,将原子核发射的回波信号经过空间编码等进一步处理即可重建图像。
由于人体内部脂肪中的氢原子核和水中的氢原子核所处的分子环境不相同,在使用相同的射频脉冲进行激发时,它们的共振频率也不相同。在不同的回波时间采集信号,脂肪组织和水表现出不同的相位以及信号强度。
狄克逊(Dixon)法是在磁共振成像中用以产生纯水质子图像的方法,其基本原理是分别采集水质子和脂肪质子的同相位(in Phase)和反相位(out phase)两种回波信号,这两种不同相位的回波信号通过运算,各自产生一幅纯水质子的图像和一幅纯脂肪质子的图像,从而在水质子图像上达到脂肪抑制的目的。Dixon法本身有一点Dixon法、两点Dixon法、三点Dixon法以及多点Dixon法等多种形式。
本领域具有多种与Dixon法相结合的k空间(k-space)数据采集方法,例如:笛卡尔(Cartesian)轨迹采集和径向(radial)或螺旋(spiral)轨迹采集等。通过研究可以发现,现有的笛卡尔轨迹采集方法虽然简单且节省时间,但是对刚体运动及脉动等运动非常敏感;而径向或螺旋轨迹采集方法会将运动伪影转化为重建后图像中的模糊,计算复杂并且耗时严重。因此,上述两类方法都不能消除刚体运动伪影。
此外,现有的三点Dixon法使用相位反卷叠技术计算水和脂肪图像,由于相位反卷叠固有的不稳定性,故而计算出的水和脂肪的图像可能调换,即理论上认为是水的图像,实际计算出来的却可能是脂肪的图像;而理论上认为是脂肪的图像,实际计算出来的却可能是水的图像;所以,在对要成像的物体进行多层扫描时,在有些扫描层面计算出的水和脂肪的图像可能发生调换,这就导致合成后的三维的水图像和脂肪图像出错。
发明内容
本发明实施例提出了一种实现水脂分离的磁共振成像方法,以降低成像过程中对运动伪影的敏感程度,提高信噪比。
本发明实施例提供了一种实现水脂分离的磁共振成像方法,该方法包括:
利用快速自旋回波刀锋伪影校正序列TSE BLADE采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;
根据所述同相位图像的原始数据重建同相位图像,并根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;以及
根据所重建的同相位图像和所重建的反相位图像计算水和脂肪的图像。
本发明实施例还提供了一种实现水脂分离的磁共振成像装置,该装置包括:
采集单元,用于利用快速自旋回波刀锋伪影校正序列TSE BLADE采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;
第一重建单元,用于根据所述同相位图像的原始数据重建同相位图像;
第二重建单元,用于根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;以及
计算单元,用于根据所重建的同相位图像和所重建的反相位图像分别计算水和脂肪的图像。
从上述技术方案可以看出,本发明实施例提供的实现水脂分离的磁共振成像方法,采用TSEBLADE序列采集k空间数据,从而继承了BLADE序列对刚体运动和脉动不敏感的优点,降低了对运动伪影的敏感程度,同时还提高了图像的信噪比。
附图说明
图1为依据本发明一实施的实现水脂分离的磁共振成像方法的流程示意图;
图2为BLADE轨迹的示意图;
图3A为依据本发明一实施例的实现水脂分离的磁共振成像方法中单极性读出梯度的序列图;
图3B为依据本发明一实施例的实现水脂分离的磁共振成像方法中双极性读出梯度的序列图;
图4为依据本发明一实施例的实现水脂分离的磁共振成像方法中对同相位图像的原始数据进行相位校正的流程示意图;
图5为依据本发明一实施例的实现水脂分离的磁共振成像方法中对反相位图像的原始数据进行相位校正的流程示意图;
图6为依据本发明一实施例的实现水脂分离的磁共振成像装置的组成结构示意图。
具体实施方式
下面将结合本发明中的附图,对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。应当理解,此处所描述的具体实施例仅仅用于解释本发明,并不用于限定本发明。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
通过研究可以发现,现有的笛卡尔轨迹采集方法虽然简单且节省时间,但是对刚体运动及脉动等运动非常敏感;而径向或螺旋轨迹采集方法会将运动伪影转化为重建后图像中的模糊,计算复杂并且耗时严重。因此,上述两类方法都不能消除刚体运动伪影。
此外,现有的三点Dixon法使用相位反卷叠技术计算水和脂肪图像,由于相位反卷叠固有的不稳定性,故而计算出的水和脂肪的图像可能调换,即理论上认为是水的图像,实际计算出来的却可能是脂肪的图像;而理论上认为是脂肪的图像,实际计算出来的却可能是水的图像;所以,在对要成像的物体进行多层扫描时,在有些扫描层面计算出的水和脂肪的图像可能发生调换,这就导致合成后的三维的水图像和脂肪图像出错。例如,对患者的头部进行扫描,从头顶到下颌一共扫描10个层面,分别为层面1、层面2……层面10,在利用现有的三点Dixon法计算该患者的头部水和脂肪的三维图像时,会先分别得到这10个层面的二维的水图像(从上到下分别为水图像1、水图像2……水图像10)和二维的脂肪图像(从上到下分别为脂肪图像1、脂肪图像2……脂肪图像10),然后再对这10个层面的二维的水图像和脂肪图像分别合成后得到三维的水图像和脂肪图像。如果在第8层面二维的水图像8和二维的脂肪图像8发生了调换,那么在这种情况下,在10个扫描层面的二维水图像中就混入了一张二维脂肪图像8,同样地,在10个扫描层面的二维脂肪图像中也混入了一张二维水图像8,这样,在利用10个扫描层面的二维的水图像和脂肪图像合成三维的水和脂肪图像时,就会导致合成后的水和脂肪的图像出错。
为了解决上述技术问题,本发明实施例提供了一种实现水脂分离的磁共振成像方法,即一种基于快速自旋回波刀锋伪影校正序列TSE BLADE的两点Dixon方法。在该方法中,先利用TSE BLADE序列采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;然后根据所述同相位图像的原始数据重建同相位图像,并根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;最后,再根据重建后的同相位图像和反相位图像计算水和脂肪的图像。
下面将结合附图详细描述本发明实施例所述的实现水脂分离的磁共振成像方法。参见图1,是本发明一实施例提供的实现水脂分离的磁共振成像方法的流程示意图,该方法具体包括以下步骤:
步骤101,磁共振成像设备利用快速自旋回波刀锋伪影校正序列(Turbo Spin Echo BLADE,TSEBLADE)采集一幅同相位图像的原始数据和一幅反相位图像的原始数据。
也即,在本步骤中,磁共振成像设备利用TSE BLADE序列采集一个同相位的回波作为同相位图像的原始数据以及采集一个反相位的回波作为反相位图像的原始数据。具体而言,在本例中,磁共振成像设备可以先采集同相位的回波再采集反相位回波,也可以先采集反相位的回波再采集同相位回波。
上述数据采集所使用的BLADE技术,也被称为螺旋浆(Periodically Rotated Overlapping ParallELLines with Enhanced Reconstruction,PROPELLER)技术,具体可参见James G.Pipe的论文“MotionCorrection With PROPELLER MRI:Application to head motion and free-breathing cardiac imaging”(Magnetic Resonance in Medicine,42:963-969,1999年11月)。由于BLADE序列对刚体运动和脉动不敏感,故而利用TSE BLADE序列来采集k空间数据,能够降低重建后的图像对运动伪影的敏感程度,同时还可以提高重建后的图像的信噪比。
采集每幅图像的原始数据的BLADE轨迹如图2所示。以N(N为正整数,图2中N取10)个数据带(strip)来采集K空间(K-space)数据,这些数据带沿圆周方向等角度或者非等角度旋转分布,每个数据带包括L(L为正整数,图2中L取9)行平行的数据线(line)。
在上述步骤101中,磁共振成像设备在采集回波信号时,既可采用单极性读出梯度(mono-polar readoutgradient)的方式,也可采用双极性读出梯度(bipolar readout gradient)的方式,来依次获取一幅同相位图像以及一幅反相位图像。
图3A和图3B示意性地给出了本发明实施例提供的方法在采集BLADE轨迹中各个数据带时的序列,其中图3A采用的是单极性读出梯度的方式,图3B采用的是双极性读出梯度的方式。在图3A和图3B中,RF表示射频脉冲,RO表示读出梯度(readout gradient),图中省略了选层梯度和相位编码梯度。
如图3A中的a部分所示,磁共振成像设备首先发射一个90度射频脉冲RF 0,然后再发射一个180度重聚相位射频脉冲RF 1。在距离90度射频脉冲RF 0二分之一个回波时间(time of echo,TE)至二分之三个回波时间之间,磁共振成像设备在读出梯度方向上施加读出梯度,读取一根数据线out 1。接着再发射一个180度重聚相位射频脉冲RF 2,得到第二个回波,并在读出梯度方向上施加读出梯度,再读取一根数据线out 2;重复上述操作,直至读取BLADE轨迹中所有的数据线,得到一幅反相位图像的原始数据。其中,数据线out 1、out 2、out 3……等构成一幅反相位图像的原始数据。
如图3A中的b部分所示,磁共振成像设备首先发射一个90度射频脉冲RF 0,然后再发射一个180度重聚相位射频脉冲RF 1。在距离90度射频脉冲RF 0一个回波时间时,磁共振成像设备在读出梯度方向上施加读出梯度,读取一根数据线in 1。接着再发射一个180度重聚相位射频脉冲RF 2,得到第二个回波,并在读出梯度方向上施加读出梯度,再读取一根数据线in 2;重复上述操作,直至读取BLADE轨迹中所有的数据线,得到一幅同相位图像的原始数据。其中,数据线in 1、in 2、in 3……等构成一幅同相位图像的原始数据。
如图3B所示,磁共振成像设备首先发射一个90度射频脉冲RF 0,然后再发射一个180度重聚相位射频脉冲RF 1。在距离90度射频脉冲RF 0一个回波时间时,磁共振成像设备在读出梯度方向上施加两个不同极性的读出梯度,分别读取两根数据线out 1和in 1。接着再发射一个180度重聚相位射频脉冲RF 2,得到第二个回波,并在读出梯度方向上施加两个不同极性的读出梯度,分别读取两根数据线out 2和in 2;重复上述操作,直至读取BLADE轨迹中所有的数据线,得到一幅反相位图像的原始数据和一幅同相位图像的原始数据。其中,数据线out 1、out 2、out 3……等构成一幅反相位图像的原始数据,数据线in 1、in 2、in 3……等构成另一幅同相位图像的原始数据。
需要说明的是,图3A和图3B只是示意性地给出了一种采集次序,本发明并不局限于此。例如,本发明可以先采集一个同相位的回波,然后采集一个反相位的回波,得到相应的原始数据;或者,先采集一个反相位的回波,然后再采集一个同相位的回波,得到相应的原始数据。
步骤102,磁共振成像设备根据所述同相位图像的原始数据重建同相位图像,并根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像。
在本步骤中,当重建同相位图像时,磁共振成像设备首先对同相位图像的各个数据带进行相位校正,然后,磁共振成像设备再对所述相位校正后的同相位图像的数据带依次进行旋转校正和平移校正,最后把所有校正后的数据带分配到笛卡尔坐标,并经过快速傅立叶变换得到一幅同相位图像。
其中,对同相位图像的各个数据带进行相位校正的流程如图4所示,具体包括以下步骤:
步骤401,利用窗函数(例如三角窗函数、金字塔窗函数)对同相位图像的原始数据的数据带进行窗操作得到第一处理结果;
步骤402,对所述第一处理结果进行二维快速傅立叶变换(Two-dimensional Fast FourierTransform,2D FFT)得到一个处理结果,不妨将该得到的处理结果称为所述同相位图像的窗数据;
步骤403,对所述同相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第二处理结果;
步骤404,从上述步骤403中得到的第二处理结果中去除从上述步骤402中得到的同相位图像的窗数据的相位,得到第三处理结果;
步骤405,对所述第三处理结果进行二维快速傅立叶逆变换(Two-dimensional inverse Fast FourierTransform,2D iFFT),从而得到相位校正后的同相位图像的数据带。
在本步骤中,当重建反相位图像时,磁共振成像设备也首先对反相位图像的各个数据带进行相位校正,然后,再对所述相位校正后的反相位图像的数据带进行旋转校正和平移校正,并经过快速傅立叶变换得到一幅反相位图像。其中,本申请的发明人对反相位图像的原始数据进行相位校正的方法提出了改进。如图5所示,在本发明实施例中对反相位图像的各个数据带进行相位校正的方法具体包括以下步骤:
步骤501,对反相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第一处理结果;
步骤502,利用窗函数对与上述反相位图像的原始数据的数据带对应的同相位图像的原始数据的数据带(即在k空间中与反相位图像的原始数据的数据带具有相同角度的同相位图像的原始数据的数据带)进行窗操作,得到第二处理结果;
步骤503,对上述步骤502中得到的第二处理结果进行二维快速傅立叶变换,得到所述同相位图像的原始数据的数据带的窗数据,称为相应的窗数据;
步骤504,从上述步骤501中得到的所述第一处理结果中去除从上述步骤503中得到的所述同相位图像的原始数据的数据带的窗数据的相位,得到第三处理结果;
步骤505,对上述步骤504中得到的第三处理结果进行二维快速傅立叶逆变换,从而得到相位校正后的反相位图像的数据带。
在上述过程中,利用同相位图像的数据带作为参考,对反相位图像的数据带进行相位校正,保留了反相位信息,从而在后续的处理过程中可以根据两点Dixon法进行水脂分离成像。
步骤103,磁共振成像设备根据上述重建后得到的一幅同相位图像和一幅反相位图像分别计算水图像和脂肪图像。
在本步骤中,在根据重建后的同相位图像和反相位图像计算水和脂肪的图像时,当所述重建后的同相位图像和反相位图像为二维图像时,利用二维区域增长算法(Two-dimensional region growingalgorithm),对重建后的同相位图像和反相位图像进行相位校正,根据相位校正后的同相位图像和反相位图像计算水和脂肪的图像。
同时,为了提高计算的可靠性,可以对要成像的物体进行多层次扫描,即同时得到反映多个扫描层面的多幅二维的同相位图像和反相位图像,然后根据所述多幅二维的同相位图像和反相位图像所对应的空间信息,把所述多幅二维的同相位图像以及反相位图像排列成一个三维的同相位图像以及一个三维的反相位图像。例如,对患者的头部进行扫描,从头顶到下颌一共扫描n个层面,分别为层面1、层面2……层面n,然后分别得到这n个层面的二维的同相位图像(从上到下分别为同相位图像1、同相位图像2……同相位图像n)和二维的反相位图像(从上到下分别为反相位图像1、反相位图像2……反相位图像n),按照这些图像的空间信息,即从头顶到下颌的排列顺序,将上述n幅二维的同相位图像排列成一个三维的同相位图像,将上述n幅二维的反相位图像排列成一个三维的反相位图像。然后,再对得到的所述三维同相位图像和三维反相位图像利用现有的各种方式来进行相位校正,例如三维区域增长算法(Three-dimensionalregion growing algorithm),这里不再一一赘述。也即,在本发明实施例中,当所述重建后的同相位图像和反相位图像为三维图像时,利用三维区域增长算法对重建后的同相位图像和反相位图像进行相位校正,根据相位校正后的同相位图像和反相位图像计算水和脂肪的图像。利用三维区域增长算法可以有效地解决传统三点Dixon法中由于采用相位反卷叠技术所带来的水和脂肪图像在不同扫描层的调换问题。
上述相位校正的目的是去除磁场不均匀所带来的相位误差。最后利用相位纠正后的三维的反相位图像以及三维的同相位图像,计算得到三维的水图像和脂肪图像。
图6为依据本发明一实施例的实现水脂分离的磁共振成像装置的组成结构示意图。该装置包括:
采集单元1,用于利用快速自旋回波刀锋伪影校正序列TSE BLADE采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;
第一重建单元2,用于根据所述同相位图像的原始数据重建同相位图像;
第二重建单元3,用于根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;以及
计算单元4,用于根据重建后的同相位图像和反相位图像分别计算水和脂肪的图像。
上述第一重建单元2包括:
第一相位校正子单元21,用于对所述同相位图像的原始数据的数据带进行相位校正;
第一旋转校正子单元22,用于对相位校正后的同相位图像的原始数据的数据带进行旋转校正;
第一平移校正子单元23,用于对旋转校正后的同相位图像的原始数据的数据带进行平移校正;
第一变换子单元24,用于对平移校正后的同相位图像的原始数据的数据带进行快速傅立叶变换。
其中,上述第一相位校正子单元21具体包括:
第一模块211,用于对所述同相位图像的原始数据的数据带进行窗操作得到第一处理结果;
第二模块212,用于对所述第一处理结果进行二维快速傅立叶变换,得到所述同相位图像的窗数据;
第三模块213,用于对所述同相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第二处理结果;
第四模块214,用于从所述第二处理结果中去除所述同相位图像的窗数据的相位得到第三处理结果;
第五模块215,用于对所述第三处理结果进行二维快速傅立叶逆变换得到所述校正结果。
上述第二重建单元3包括:
第二相位校正子单元31,用于对所述同相位图像的原始数据的数据带进行相位校正;
第二旋转校正子单元32,用于对相位校正后的同相位图像的原始数据的数据带进行旋转校正;
第二平移校正子单元33,用于对旋转校正后的同相位图像的原始数据的数据带进行平移校正;
第二变换子单元34,用于对平移校正后的同相位图像的原始数据的数据带进行快速傅立叶变换。
上述第二相位校正子单元31具体包括:
第六模块311,用于对反相位图像的原始数据的数据带进行二维快速傅立叶变换得到第一处理结果;
第七模块312,用于对与所述反相位图像的原始数据的数据带对应的同相位图像的原始数据的数据带进行窗操作得到第二处理结果;
第八模块313,用于对所述第二处理结果进行二维快速傅立叶变换,得到所述同相位图像的原始数据的相应数据带的窗数据;
第九模块314,用于从所述第一处理结果中去除所述同相位图像的原始数据的数据带的窗数据的相位得到第三处理结果;
第十模块315,用于对所述第三处理结果进行二维快速傅立叶逆变换得到所述校正结果。
从上述本发明实施例的技术方案可以看出,本发明实施例提供的实现水脂分离的磁共振成像方法,采用TSE BLADE序列采集k空间数据,从而继承了BLADE序列对刚体运动和脉动不敏感的优点,降低了对运动伪影的敏感程度,同时还提高了图像的信噪比。
更进一步,与现有的三点Dixon法相比,本发明利用二维区域增长算法或三维区域增长算法来重建图像,消除了在多层扫描中水和脂肪图像经常无规律的调换问题。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (12)

1.一种实现水脂分离的磁共振成像方法,该方法包括:
利用快速自旋回波刀锋伪影校正序列采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;
根据所述同相位图像的原始数据重建同相位图像,并根据所述同相位图像的原始数据以及所述反相位图像的原始数据重建反相位图像;以及
根据所重建的同相位图像和所重建的反相位图像计算水和脂肪的图像。
2.根据权利要求1所述的方法,其特征在于,所述根据所述同相位图像的原始数据重建同相位图像包括:
对所述同相位图像的原始数据的数据带依次进行相位校正、旋转校正、平移校正以及快速傅立叶变换得到所述同相位图像。
3.根据权利要求2所述的方法,其特征在于,所述对所述同相位图像的原始数据的数据带进行相位校正包括:
对所述同相位图像的原始数据的数据带进行窗操作得到第一处理结果;
对所述第一处理结果进行二维快速傅立叶变换,得到所述同相位图像的窗数据;
对所述同相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第二处理结果;
从所述第二处理结果中去除所述同相位图像的窗数据的相位得到第三处理结果;以及
对所述第三处理结果进行二维快速傅立叶逆变换得到校正结果。
4.根据权利要求1所述的方法,其特征在于,所述根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像包括:
对所述反相位图像的原始数据的数据带进行相位校正、旋转校正、平移校正以及快速傅立叶变换得到所述反相位图像。
5.根据权利要求4所述的方法,其特征在于,所述对所述反相位图像的原始数据的数据带进行相位校正包括:
对所述反相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第一处理结果;
对与所述反相位图像的原始数据的数据带对应的同相位图像的原始数据的数据带进行窗操作,得到第二处理结果;
对所述第二处理结果进行二维快速傅立叶变换,得到所述同相位图像的原始数据的数据带的窗数据;
从所述第一处理结果中去除所述同相位图像的原始数据的数据带的的窗数据的相位,得到第三处理结果;
对所述第三处理结果进行二维快速傅立叶逆变换,得到校正结果。
6.根据权利要求1所述的方法,其特征在于,所重建的同相位图像和所重建的反相位图像为二维图像;
所述根据所重建的同相位图像和所重建的反相位图像计算水和脂肪的图像包括:
利用二维区域增长算法,对所重建的同相位图像和所重建的反相位图像进行相位校正,根据相位校正后的同相位图像和反相位图像计算水和脂肪的图像。
7.根据权利要求1所述的方法,其特征在于,所重建的同相位图像和所重建的反相位图像为三维图像;
所述根据所重建的同相位图像和所重建的反相位图像计算水和脂肪的图像包括:
利用三维区域增长算法,对所重建的同相位图像和所重建的反相位图像进行相位校正,根据相位校正后的同相位图像和反相位图像计算水和脂肪的图像。
8.一种实现水脂分离的磁共振成像装置,该装置包括:
采集单元(1),用于利用快速自旋回波刀锋伪影校正序列采集一幅同相位图像的原始数据和一幅反相位图像的原始数据;
第一重建单元(2),用于根据所述同相位图像的原始数据重建同相位图像;
第二重建单元(3),用于根据所述同相位图像的原始数据以及反相位图像的原始数据重建反相位图像;以及
计算单元(4),用于根据所重建的同相位图像和所重建的反相位图像分别计算水和脂肪的图像。
9.根据权利要求8所述的装置,其特征在于,所述第一重建单元(2)包括:
第一相位校正子单元(21),用于对所述同相位图像的原始数据的数据带进行相位校正;
第一旋转校正子单元(22),用于对相位校正后的同相位图像的原始数据的数据带进行旋转校正;
第一平移校正子单元(23),用于对旋转校正后的同相位图像的原始数据的数据带进行平移校正;以及
第一变换子单元(24),用于对平移校正后的同相位图像的原始数据的数据带进行快速傅立叶变换。
10.根据权利要求9所述的装置,其特征在于,所述第一相位校正子单元(21)包括:
第一模块(211),用于对所述同相位图像的原始数据的数据带进行窗操作得到第一处理结果;
第二模块(212),用于对所述第一处理结果进行二维快速傅立叶变换,得到所述同相位图像的窗数据;
第三模块(213),用于对所述同相位图像的原始数据的数据带进行二维快速傅立叶变换,得到第二处理结果;
第四模块(214),用于从所述第二处理结果中去除所述同相位图像的窗数据的相位得到第三处理结果;以及
第五模块(215),用于对所述第三处理结果进行二维快速傅立叶逆变换得到校正结果。
11.根据权利要求8所述的装置,其特征在于,所述第二重建单元(3)包括:
第二相位校正子单元(31),用于对所述同相位图像的原始数据的数据带进行相位校正;
第二旋转校正子单元(32),用于对相位校正后的同相位图像的原始数据的数据带进行旋转校正;
第二平移校正子单元(33),用于对旋转校正后的同相位图像的原始数据的数据带进行平移校正;以及
第二变换子单元(34),用于对平移校正后的同相位图像的原始数据的数据带进行快速傅立叶变换。
12.根据权利要求11所述的装置,其特征在于,所述第二相位校正子单元(31)包括:
第六模块(311),用于对反相位图像的原始数据的数据带进行二维快速傅立叶变换得到第一处理结果;
第七模块(312),用于对与所述反相位图像的原始数据的数据带对应的同相位图像的原始数据的数据带进行窗操作得到第二处理结果;
第八模块(313),用于对所述第二处理结果进行二维快速傅立叶变换,得到所述同相位图像的原始数据的数据带的窗数据;
第九模块(314),用于从所述第一处理结果中去除所述同相位图像的原始数据的数据带的窗数据的相位得到第三处理结果;以及
第十模块(315),用于对所述第三处理结果进行二维快速傅立叶逆变换得到校正结果。
CN201210218752.6A 2012-06-28 2012-06-28 一种实现水脂分离的磁共振成像方法和装置 Active CN103505210B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201210218752.6A CN103505210B (zh) 2012-06-28 2012-06-28 一种实现水脂分离的磁共振成像方法和装置
US13/927,359 US9390476B2 (en) 2012-06-28 2013-06-26 Magnetic resonance imaging method and device achieving water/fat separation
KR1020130074429A KR101644677B1 (ko) 2012-06-28 2013-06-27 물/지방 분리를 달성하는 자기 공명 영상 방법 및 장치
JP2013135001A JP6245861B2 (ja) 2012-06-28 2013-06-27 水・脂肪の分離を実行する磁気共鳴画像化方法および装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210218752.6A CN103505210B (zh) 2012-06-28 2012-06-28 一种实现水脂分离的磁共振成像方法和装置

Publications (2)

Publication Number Publication Date
CN103505210A true CN103505210A (zh) 2014-01-15
CN103505210B CN103505210B (zh) 2015-09-16

Family

ID=49778240

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210218752.6A Active CN103505210B (zh) 2012-06-28 2012-06-28 一种实现水脂分离的磁共振成像方法和装置

Country Status (4)

Country Link
US (1) US9390476B2 (zh)
JP (1) JP6245861B2 (zh)
KR (1) KR101644677B1 (zh)
CN (1) CN103505210B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105809662A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN105809661A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN107153169A (zh) * 2017-07-04 2017-09-12 大连锐谱科技有限责任公司 一种稳态进动梯度多回波水脂分离成像方法
CN107997763A (zh) * 2017-12-26 2018-05-08 中国科学院生物物理研究所 一种水脂分离磁共振成像方法
CN110068783A (zh) * 2019-04-29 2019-07-30 华东师范大学 一种消除脂肪伪影的化学交换饱和转移成像方法
WO2020006959A1 (zh) * 2018-07-02 2020-01-09 清华大学 基于平面回波成像的磁共振水脂分离和定量方法及装置
CN111381204A (zh) * 2018-12-28 2020-07-07 西门子(深圳)磁共振有限公司 基于二维快速自旋回波的磁共振成像方法和装置
CN111856366A (zh) * 2019-04-24 2020-10-30 西门子(深圳)磁共振有限公司 k空间数据采集装置及方法、磁共振成像装置及方法
CN111948590A (zh) * 2020-07-13 2020-11-17 上海东软医疗科技有限公司 磁共振成像方法及装置、电子设备、存储介质
CN112285623A (zh) * 2019-07-25 2021-01-29 西门子(深圳)磁共振有限公司 磁共振成像水脂图像分离方法、装置及磁共振成像系统

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6467344B2 (ja) * 2012-09-04 2019-02-13 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ディクソン水脂肪分離を用いるプロペラ
JP2017506976A (ja) * 2014-03-05 2017-03-16 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 動き補正、水−脂肪分離、及び磁場不均一性情報の推定を用いるmriプロペラ
JP7208796B2 (ja) * 2016-06-02 2023-01-19 コーニンクレッカ フィリップス エヌ ヴェ ディクソン法による水/脂肪分離を用いたmr撮像
WO2018114554A1 (en) * 2016-12-20 2018-06-28 Koninklijke Philips N.V. Dixon-type water/fat separation mr imaging
WO2018187338A1 (en) * 2017-04-03 2018-10-11 Board Of Regents, The University Of Texas System System and method for robust mr imaging with prepared contrast using cartesian acquisition with spiral reordering (caspr)
KR102062766B1 (ko) 2017-12-11 2020-01-06 삼성전자주식회사 물 지방 분리 영상을 획득하는 방법 및 그 자기 공명 영상 장치
DE102018208569A1 (de) 2018-05-30 2019-12-05 Siemens Healthcare Gmbh Aufnahme zweier Magnetresonanz-Bilder
US10884087B2 (en) * 2019-03-01 2021-01-05 Canon Medical Systems Corporation Apparatuses and methods to improve chemical species separation in magnetic resonance imaging

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030062900A1 (en) * 2001-07-04 2003-04-03 Berthold Kiefer Magnetic resonance tomography apparatus and method for separating fat and water images by correction of phase values dependent on a noise phase
CN1827038A (zh) * 2005-02-28 2006-09-06 西门子(中国)有限公司 磁共振多通道成像水脂分离重建算法
US20110254547A1 (en) * 2010-04-15 2011-10-20 Reeder Scott B Method for separating magnetic resonance imaging signals using spectral distinction of species
CN102232830A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种磁共振成像水脂分离方法
CN102232831A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法
WO2014037870A1 (en) * 2012-09-04 2014-03-13 Koninklijke Philips N.V. Propeller with dixon water fat separation

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002052005A (ja) * 2000-08-10 2002-02-19 Hitachi Medical Corp 磁気共鳴イメージング方法
WO2004097387A2 (en) * 2003-04-25 2004-11-11 Case Western Reserve University Dixon techniques in spiral trajectories with off-resonance correction
US7227359B2 (en) * 2003-11-26 2007-06-05 Boards Of Regents, The University Of Texas System Method and apparatus for phase-sensitive magnetic resonance imaging
US8314614B2 (en) * 2007-05-03 2012-11-20 Koninklijke Philips Electronics N.V. PROPELLER MRI with phase correction
WO2009012115A1 (en) * 2007-07-13 2009-01-22 Board Of Regents, The University Of Texas System Methods for efficient and improved phase-sensitive mri
DE102011007850B4 (de) * 2011-04-21 2015-05-13 Siemens Aktiengesellschaft Verfahren zur getriggerten Akquisition von Messdaten mittels eines Magnetresonanzgeräts unter Unterdrückung von unerwünschten Signalen, Magnetresonanzgerät, Computerprogramm und elektronisch lesbarer Datenträger

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030062900A1 (en) * 2001-07-04 2003-04-03 Berthold Kiefer Magnetic resonance tomography apparatus and method for separating fat and water images by correction of phase values dependent on a noise phase
CN1827038A (zh) * 2005-02-28 2006-09-06 西门子(中国)有限公司 磁共振多通道成像水脂分离重建算法
US20110254547A1 (en) * 2010-04-15 2011-10-20 Reeder Scott B Method for separating magnetic resonance imaging signals using spectral distinction of species
CN102232830A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种磁共振成像水脂分离方法
CN102232831A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法
WO2014037870A1 (en) * 2012-09-04 2014-03-13 Koninklijke Philips N.V. Propeller with dixon water fat separation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JINGFEI MA: "Breath-hold water and fat imaging using a dual-echo two-point Dixon technique with an efficient and robust phase-correction algorithm", 《MAGNETIC RESONANCE IN MEDICINE》 *
燕树林等: "《全国医用设备使用人员(CT、MR、DSA)上岗考试指南》", 30 June 2009, 《军事医学科学出版社》 *
翁得河: "Dixon水脂分离法介绍", 《世界医疗器械》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105809661B (zh) * 2014-12-30 2018-08-07 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN105809661A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN105809662A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN105809662B (zh) * 2014-12-30 2018-09-04 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN107153169B (zh) * 2017-07-04 2020-03-27 南京拓谱医疗科技有限公司 一种稳态进动梯度多回波水脂分离成像方法
CN107153169A (zh) * 2017-07-04 2017-09-12 大连锐谱科技有限责任公司 一种稳态进动梯度多回波水脂分离成像方法
CN107997763A (zh) * 2017-12-26 2018-05-08 中国科学院生物物理研究所 一种水脂分离磁共振成像方法
WO2020006959A1 (zh) * 2018-07-02 2020-01-09 清华大学 基于平面回波成像的磁共振水脂分离和定量方法及装置
CN111381204A (zh) * 2018-12-28 2020-07-07 西门子(深圳)磁共振有限公司 基于二维快速自旋回波的磁共振成像方法和装置
CN111381204B (zh) * 2018-12-28 2022-07-12 西门子(深圳)磁共振有限公司 基于二维快速自旋回波的磁共振成像方法和装置
CN111856366A (zh) * 2019-04-24 2020-10-30 西门子(深圳)磁共振有限公司 k空间数据采集装置及方法、磁共振成像装置及方法
CN111856366B (zh) * 2019-04-24 2023-05-26 西门子(深圳)磁共振有限公司 k空间数据采集装置及方法、磁共振成像装置及方法
CN110068783A (zh) * 2019-04-29 2019-07-30 华东师范大学 一种消除脂肪伪影的化学交换饱和转移成像方法
CN110068783B (zh) * 2019-04-29 2021-02-26 华东师范大学 一种消除脂肪伪影的化学交换饱和转移成像方法
CN112285623A (zh) * 2019-07-25 2021-01-29 西门子(深圳)磁共振有限公司 磁共振成像水脂图像分离方法、装置及磁共振成像系统
CN111948590A (zh) * 2020-07-13 2020-11-17 上海东软医疗科技有限公司 磁共振成像方法及装置、电子设备、存储介质
CN111948590B (zh) * 2020-07-13 2023-06-16 上海东软医疗科技有限公司 磁共振成像方法及装置、电子设备、存储介质

Also Published As

Publication number Publication date
JP2014008409A (ja) 2014-01-20
US9390476B2 (en) 2016-07-12
JP6245861B2 (ja) 2017-12-13
US20140003694A1 (en) 2014-01-02
KR101644677B1 (ko) 2016-08-01
CN103505210B (zh) 2015-09-16
KR20140001779A (ko) 2014-01-07

Similar Documents

Publication Publication Date Title
CN103505210A (zh) 一种实现水脂分离的磁共振成像方法和装置
CN103477238B (zh) 采用来自先验采集的约束的压缩感测mr图像重建
CN107750338B (zh) 使用星形叠层采集的mr成像方法和设备
CN102435966B (zh) 三维磁共振成像方法及系统
JP4951264B2 (ja) シェルk空間サンプリング軌跡を使用したアンダーサンプリング3DMRI
CN103257333B (zh) 一种磁共振成像中的水脂分离成像方法及装置
CN102232831B (zh) 一种实现水脂分离的磁共振成像方法
JP6467344B2 (ja) ディクソン水脂肪分離を用いるプロペラ
US10852376B2 (en) Magnetic resonance imaging method and device
CN105074491B (zh) 具有使用压缩感知的图像重建的动态mri
CN1918480A (zh) 磁共振成像方法
CN1683939A (zh) 在磁共振断层造影中融合加权成像的运动校正的多点方法
JP2010221026A (ja) k空間上の空間周波数領域フィルタリングを利用したMRIシステムにおけるB1マッピング
CN103380384B (zh) 使用交错的水参考扫描的具有自动相位和b0校正的磁共振波谱成像方法和装置
CN102232830A (zh) 一种磁共振成像水脂分离方法
CN101470179A (zh) 磁共振成像中失真校准的方法和装置
CN108776318A (zh) 磁共振多对比度图像重建
CN112014782B (zh) 磁共振成像方法、装置及计算机存储介质
CN105473067A (zh) 用于降低mr成像中声音噪声水平的系统和方法
CN105738846A (zh) K空间数据采集方法及其磁共振成像方法
CN103584864A (zh) 一种磁共振成像方法和装置
CN104395773A (zh) 贯穿平面的导航器
WO2016055462A1 (en) Propeller mr imaging with artefact suppression
CN102096055B (zh) 一种用于磁共振成像非均匀采样数据的快速精确重建方法
CN105902270A (zh) 一种超快速单次激发四点Dixon水脂分离成像方法

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