CN101856228A - 磁共振成像装置 - Google Patents

磁共振成像装置 Download PDF

Info

Publication number
CN101856228A
CN101856228A CN201010154533A CN201010154533A CN101856228A CN 101856228 A CN101856228 A CN 101856228A CN 201010154533 A CN201010154533 A CN 201010154533A CN 201010154533 A CN201010154533 A CN 201010154533A CN 101856228 A CN101856228 A CN 101856228A
Authority
CN
China
Prior art keywords
mentioned
magnetic resonance
image
resonance signal
signal
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
CN201010154533A
Other languages
English (en)
Other versions
CN101856228B (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.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN101856228A publication Critical patent/CN101856228A/zh
Application granted granted Critical
Publication of CN101856228B publication Critical patent/CN101856228B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
    • 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • 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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308—Characterization of motion or flow; Dynamic imaging
    • G01R33/56316—Characterization of motion or flow; Dynamic imaging involving phase contrast techniques
    • 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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366—Perfusion 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/5607—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood 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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]

Landscapes

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

Abstract

提供一种磁共振成像装置，具备：收集单元，分别利用第1序列和第2序列收集在投放造影剂前后在成像区域内产生的磁共振信号，其中，第1序列是使射频激励后的磁化失相从而使关于在第1流速范围内流动的流体的第1磁共振信号分量与关于在比第1流速范围低的第2流速范围流动的流体的第2磁共振信号分量相比信号大幅度降低的第1序列，第2序列是设定为与造影剂的浓度对应的大小的序列；重建单元，根据利用第1以及第2序列分别收集到的磁共振信号，分别重建将流体在成像区域中的空间分布反映磁共振信号的大小而表示的第1图像以及第2图像；以及生成单元，根据第1图像以及第2图像生成表示投放造影剂前后的流体变化程度的第3图像。

Description

磁共振成像装置
(相关申请的交叉引用)
本申请基于2009年4月3日提交的在先的日本专利申请No.2009-091303并要求其优先权，其全部内容通过引用结合在本申请中。
技术领域
本发明涉及取得用于根据从被检体(subject)放射出的磁共振信号进行医用诊断的信息的技术。
背景技术
在dynamic susceptibility contrast-magnetic resonanceimaging(DSC-MRI：动态磁敏感性对比磁共振成像)或dynamiccontrast-enhanced MRI(DCE-MRI：动态对比增强磁共振成像)等使用血管内造影剂的灌流成像(perfusion imaging)中，由于血管成为诊断上的障碍，因此希望在解析映像上将其排除。能够对已经描出血管的图像通过阈值处理等的图像处理排除血管，但是，不只是需要多余的处理，阈值的最佳设定也是困难的。另外，在这种图像中，通常，为了重视信号噪声比(SNR)与时间分辨率，大多数情况下在低空间分辨率下进行成像。因此，根据部分容积(partial volume)效应，血管附近的实质部也被排除。另外，在去卷积(deconvolution)法等使用动脉输入功能(arterial input function)(AIF)的方法中，由于在解析中使用血管部分所以不能排除。自旋回波(SE)体系与梯度回波(GRE)体系相比血管信号抑制效果好，但是并不充分。
在使用GRE体系的平面回波成像(echo planar imaging)(EPI)法的功能磁共振成像(functional MRI)(fMRI)中，能够观察组织(包含毛细血管的组织)的活动变化，但是比较粗的血管内的流入(in-flow)效应或血氧水平依赖(blood oxygenation level dependent)(BOLD)效应，特别是在观察与血管毗邻的组织的活动变化上成为障碍。在不使用通常的EPI法的自旋卷绕(spin-warp)型的GRE中，特别是流动效应相对于原来的BOLD成为伪影。在图像处理上存在与上述相同的问题。另外，最近报告的diffusion-weighted(DW)-fMRI(参考“Le Bihanet al.PNAS 103，8263-8268，2006)，其主要目的是使用1000以上的大b值(b-factor)，试着根据组织内间质液或细胞内液等的D系数的差分离组织细胞内的D系数的变动。
像这样，以往，存在当描出像毛细血管中流动的血液或造影剂那样的以比较低的速度流动的流体时，像动脉或静脉中流动的血液或造影剂那样的以比较高的速度流动的流体也会被描出的问题。
发明内容
鉴于上述情况，希望能够将以比较低的速度流动的流体相对于以比其更高的速度流动的流体强调地描出。
根据本发明的第1实施方式提供一种磁共振成像装置，其特征在于，包括：收集单元，利用第1序列收集在投放造影剂前在成像区域内产生的磁共振信号，并且利用第2序列收集在投放上述造影剂后在上述成像区域内产生的磁共振信号，其中，第1序列是使射频激励后的磁化失相从而使关于在第1流速范围内流动的流体的第1磁共振信号分量与关于在比上述第1流速范围低的第2流速范围流动的上述流体的第2磁共振信号分量相比信号大幅度降低的序列，第2序列是设定为与上述造影剂的浓度对应的大小的序列；重建单元，根据利用上述第1以及第2序列分别收集到的上述磁共振信号，分别重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号大小而表示的第1图像以及第2图像；以及生成单元，根据上述第1图像以及上述第2图像生成表示了在投放上述造影剂后相对于投放前的上述流体的变化程度的第3图像。
根据本发明的第2实施方式提供一种磁共振成像装置，其特征在于，包括：收集单元，利用失相序列收集在成像区域内产生的磁共振信号，该失相序列是使射频激励后的磁化失相从而使在第1流速范围内流动的流体与在比上述第1流速范围低的第2流速范围流动的上述流体相比信号大幅度降低的序列；以及重建单元，根据通过上述收集单元所收集到的上述磁共振信号重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号的大小而表示的图像。
根据本发明的第3实施方式提供一种磁共振成像装置，其特征在于，包括：收集单元，使用使包含倾斜磁场失相脉冲的血管信号降低并较黑地描出血管的黑血法(black blood法)对投放阴性造影剂前的被检体进行成像，使用规定的成像方法对投放上述阴性造影剂后的被检体进行动态成像；重建单元，根据对上述造影剂投放前的被检体进行成像所取得的磁共振信号重建基准图像，并且根据对上述造影剂投放后的被检体进行动态成像所取得的磁共振信号重建多时相的造影图像；以及生成单元，根据上述基准图像以及上述造影图像生成诊断用图像。
根据本发明的第4实施方式提供一种磁共振成像装置，其特征在于，包括：收集单元，分别以第1条件和第2条件收集在成像区域内产生的磁共振信号，其中，第1条件是使关于在第1流速范围内流动的流体的第1磁共振信号分量的信号水平相对于关于在比上述第1流速范围低的第2流速范围流动的上述流体的第2磁共振信号分量的信号水平在规定方向上不同的条件，第2条件是使上述第1磁共振信号分量的信号水平相对于上述第2磁共振信号分量的信号水平在上述规定方向上不同，并且上述第1以及第2磁共振信号分量的信号水平都相对于上述第1条件在上述规定方向上不同，而且对于其变化量，上述第2磁共振信号分量相对于上述第1磁共振信号分量增大的条件；重建单元，根据分别以上述第1以及第2条件收集到的上述磁共振信号，分别重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号大小而表示的第1图像以及第2图像；以及生成单元，根据上述第1图像以及上述第2图像生成表示了在上述第1条件与上述第2条件下的信号水平的变化程度的第3图像。
在下面的描述中将说明本发明的其它目的和优点，部分内容可以从说明书的描述中明显看出或者可以通过实施本发明而理解。本发明的目的和优点可以通过下文中具体指出的实施方法和其组合实现和获得。
附图说明
结合在这里并构成说明书的一部分的附图描述本发明当前优选的实施方式，并且与上述的总的说明以及下面的优选实施方式的详细描述一同用来说明本发明的原理。
图1为表示一实施方式的磁共振成像装置(MRI装置)的概略结构的图。
图2为表示用于使横磁化失相的脉冲序列的一例的图。
图3为表示用于使纵磁化失相的脉冲序列的一例的图。
图4为表示使用阴性造影剂时的原信号强度-时间曲线的图。
图5为表示根据图4所示的原信号强度-时间曲线求出的浓度-时间曲线的图。
图6为表示CBF(脑血流)映像的一例的图。
图7为表示一般DSC-MRA时的原信号强度-时间曲线的图。
图8为表示根据图7所示的原信号强度-时间曲线求出的浓度-时间曲线的图。
图9为表示一般DSC-MRA时的CBF映像的一例的图。
图10为表示算出一般的动态成像中的诊断用图像的信号值的情况的图。
图11为表示算出一实施方式的动态成像中的诊断用图像的信号值的情况的第1例子的图。
图12为表示算出一实施方式的动态成像中的诊断用图像的信号值的情况的第2例子的图。
图13为表示算出一实施方式的动态成像中的诊断用图像的信号值的情况的第3例子的图。
图14为表示通过组合不同的b强度分别取得基准图像时的原信号强度-时间曲线的图。
图15为表示根据图14所示的原信号强度-时间曲线求出的浓度-时间曲线的图。
图16为表示背景相位校正的情况的图。
图17为表示余弦滤波的特性的图。
图18为表示进行余弦滤波处理时的摄像步骤的一例的流程图。
图19为表示进行余弦滤波处理时的摄像步骤的一例的流程图
图20为表示在背景相位校正以及余弦滤波处理后造影前后的血管信号的一例的图。
图21为表示在背景相位校正以及余弦滤波处理后造影前后的血管信号的一例的图。
图22为使用血管与背景的空间分布表示余弦滤波处理的效果的图。
具体实施方式
以下，参照附图说明本发明的一实施方式。
图1为表示本实施方式的磁共振成像装置(MRI装置)100的概略结构的图。
该MRI装置100具备：床部、静磁场产生部、倾斜磁场产生部、发送接收部以及控制·运算部。床部使载置的被检体200移动。静磁场产生部产生静磁场。倾斜磁场产生部产生用于向静磁场附加位置信息的倾斜磁场。发送接收部发送接收射频信号。控制运算部担负系统整体的控制以及图像重建。并且，MRI装置100中，作为上述各部的构成要素具有磁铁1、静磁场电源2、匀场线圈3、匀场线圈电源4、床板5、倾斜磁场线圈单元6、倾斜磁场电源7、RF线圈单元8、发送器9T、接收器9R、定序器(序列控制器)10、运算单元11、存储单元12、显示器13、输入器14、声音发生器15以及主计算机16。另外，MRI装置100与计测作为表示被检体200的心脏跳动的信号的ECG信号的心电计测部连接。
静磁场产生部包含磁铁1、静磁场电源2、匀场线圈3以及匀场线圈电源4。作为磁铁1，例如可以使用超导磁铁或常导磁铁。静磁场电源2向磁铁1供给电流。另外，作为磁铁1采用超导磁铁时，静磁场电源2可以省略。这样，静磁场产生部在送入被检体200的圆筒状的空间(诊断用空间)中产生静磁场B0。该静磁场B0的磁场方向与诊断用空间的轴方向(Z轴方向)大致一致。匀场线圈3接受来自主计算机16的控制下的匀场线圈电源4的电流供给，产生用于使静磁场均匀化的校正磁场。
床部将载置了被检体200的床板5送入诊断用空间，或从诊断用空间抽出。
倾斜磁场产生部包含倾斜磁场线圈单元6以及倾斜磁场电源7。倾斜磁场线圈单元6被配置在磁铁1的内侧。倾斜磁场线圈单元6具备用于产生互相正交的X轴方向、Y轴方向以及Z轴方向上的各倾斜磁场的3组线圈6x、6y、6z。倾斜磁场电源7在定序器10的控制下，向线圈6x、6y、6z分别供给用于产生倾斜磁场的脉冲电流。倾斜磁场产生部通过控制从倾斜磁场电源7向线圈6x、6y、6z供给的脉冲电流，合成作为物理轴的3轴(X轴、Y轴、Z轴)方向上的各倾斜磁场，任意设定互相正交的切片方向倾斜磁场Gs、相位编码方向倾斜磁场Ge以及读出方向(频率编码方向)倾斜磁场Gr所构成的逻辑轴方向上的各倾斜磁场。切片方向、相位编码方向以及读出方向上的各倾斜磁场Gs、Ge、Gr与静磁场B0重叠。
发送接收部包含RF线圈单元8、发送器9T以及接收器9R。RF线圈单元8被配置在诊断用空间中的被检体200的附近。发送器9T以及接收器9R与RF线圈单元8连接。发送器9T以及接收器9R在定序器10的控制下动作。发送器9T向RF线圈单元8供给用于产生核磁共振(NMR)的拉莫尔频率的射频(RF)电流脉冲。接收器9R获取由RF线圈单元8接收到的回波信号等MR信号(射频信号)，并对其实施前置放大、中频转换、相位检波、低频放大、或滤波等各种信号处理后，进行A/D转换，生成数字数据(原始数据)。
控制运算部包含定序器10、运算单元11、存储单元12、显示器13、输入器14、声音发生器15以及主计算机16。
定序器10具备CPU以及存储器。定序器10将从主计算机16发送来的脉冲序列信息存储到存储器中。定序器10的CPU按照存储器中存储的序列信息控制倾斜磁场电源7、发送器9T以及接收器9R的动作。定序器10的CPU一旦输入由接收器9R输出的原始数据，就将其转送到运算单元11。在此，序列信息是指使倾斜磁场电源7、发送器9T以及接收器9R按照一连串的脉冲序列动作所必需的全部信息，例如包含与向线圈6x、6y、6z施加的脉冲电流的强度、施加时间以及施加定时等相关的信息。
运算单元11通过定序器10输入由接收器9R输出的原始数据。运算单元11将输入的原始数据配置在内部存储器中设定的k空间(也称为傅立叶空间或频率空间)。运算单元11通过将配置在该k空间的数据作为对象进行二维或三维傅立叶转换而重建实空间的图像数据。另外，运算单元11根据需要也可以执行与图像相关的数据的合成处理或差分运算处理(也包含加权差分处理)。该合成处理包含对每个像素将像素值相加的处理或最大强度投影(MIP)处理、最小强度投影(minIP)处理等。另外，作为上述合成处理的其他例子，也可以在傅立叶空间上进行多个帧的轴的调整后，合成这些多个帧的原始数据，取得1帧的原始数据。另外，加法处理包含单纯加法处理、加法平均处理或加权加法处理等。
存储单元12存储重建后的图像数据或实施上述合成处理或差分处理后的图像数据。
显示器13在主计算机16的控制下显示应向用户提示的各种图像。作为显示器13，可以使用液晶显示器等显示装置。
输入器14输入操作者所希望的同步定时选择用的参数信息、扫描条件、脉冲序列、与图像合成或差分运算相关的信息等各种信息。输入器14将输入的信息发送到主计算机16。作为输入器14适当地具备鼠标或轨迹球等定位装置、模式切换开关等选择装置、或键盘等输入装置。
声音发生器15在由主计算机16发出指令时，将闭气(停止呼吸)开始以及闭气结束的信息作为声音发出。
主计算机16为了实现通过已有的MRI装置实现的各种动作而总括MRI装置100的各部的动作。
心电计测部包含ECG传感器17以及ECG单元18。ECG传感器17附着于被检体200的体表，将被检体200的ECG信号作为电信号(以下称为传感器信号)检测出。ECG单元18在对传感信号实施包含数字化处理的各种处理后，向主计算机16以及定序器10输出。作为该心电计测部，例如可以使用向量心电计。在执行与被检体200的心时相同步的扫描时，根据需要在定序器10中使用基于该心电计测部的传感器信号。
其次，针对以上构成的MRI装置100的动作进行说明。另外，MRI装置100能够进行通过已有的MRI装置实现的各种成像，但是省略对其说明。并且，在此，对通过DSC-MRI一边排除动脉或静脉等的血液流速比较高的血管(以下，简单地称为血管)一边描出包含毛细血管的组织(以下，简单地称为组织)时的动作进行说明。
(基本动作)
MRI装置100以根据运动探测梯度(motion probinggradient)(MPG)运动大的部分信号大幅度降低那样的适当的b值(b-factor)利用使RF激励后的横磁化或纵磁化失相(dephase)的脉冲序列进行摄影。
图2为表示用于使横磁化失相的脉冲序列的一例的图。图3为用于使纵磁化失相的脉冲序列的一例的图。
在图2所示的脉冲序列中，施加RF脉冲后，在使用于收集磁共振信号的收集部开始之前，施加MPG。此时，作为收集部中的脉冲序列，最好是GRE或EPI等。
在图3所示的脉冲序列中，在使收集部开始之前，施加MPG前置脉冲。MPG前置脉冲例如形成为90°脉冲-MPG-180°脉冲-MPG-90°脉冲那样的形状。并且，作为此时收集部中的脉冲序列，可以使用(快速场回波(fast field echo)(FFE)、稳态自由进动(steady state freeprecession)(SSFP)、或快速自旋回波(fast spin echo)(FSE)等。该图3所示的脉冲序列称为准备(preparation)方式。
本实施方式的特征在于通过使用了包含倾斜磁场失相脉冲的脉冲序列的黑血法(black blood法)对投放了阴性造影剂的被检体进行动态成像。黑血法是降低血管信号并较黑地描出血管的成像法。作为倾斜磁场失相脉冲也可以是一般使用的脉冲，但是在本实施方式中，作为倾斜磁场失相脉冲使用图2或图3所示的MPG脉冲。
由MPG脉冲所产生的失相量设为对血管以外的实质部以及其造影分量没有影响的程度的b值。b值与血流(造影剂)等流体的流速或其空间分布(体像素内的方向或方向变化的程度)有关。具体而言，b值例如在2-10[sec/nm2]的范围内，更优选是在2-4[sec/nm2]的范围内。
另外，在本实施方式中，对用于取得造影剂投放前后的两个图像的成像序列应用包含倾斜磁场失相脉冲的黑血法。然而，也可以采用对用于取得造影剂投放前的基准图像的成像序列应用包含倾斜磁场失相脉冲的黑血法、对用于取得投放造影剂后的图像的成像序列应用不包含倾斜磁场失相脉冲(仅施加倾斜磁场重相脉冲、或倾斜磁场重相脉冲和时相脉冲都没有施加)的黑血法的成像序列。
但是，描出包含毛细血管的组织时，一般使用造影剂。一边将b值设为上述那样的值一边使通过上述脉冲序列进行造影成像时的原信号强度-时间曲线(time-signal curve：TSC)成为如图4所示的曲线图。另外，tissue TSC是与组织相关的TSC，artery TSC是与动脉相关的TSC。但是，图4为表示使用阴性造影剂时的原信号强度-时间曲线的图。
关于tissue TSC的造影剂浓度C(t)以及关于artery TSC的造影剂浓度Ca(t)分别通过下式求出。但是，S(t)是tissue TSC的原信号强度。Sa(t)是artery TSC的原信号强度。S0、S(t)是造影剂投放前后的像。TE是回波时间。
C(t)＝ln[S0/S(t)]/TE
Ca(t)＝ln[S0/Sa(t)]/TE
因此，根据图4所示的原信号强度-时间曲线，求出图5所示的浓度-时间曲线。
据该图5所知，artery TSC的浓度峰值在tissue TSC的浓度峰值以下。其结果，脑血流(cerebral blood flow)(CBF)映像的动脉区域被抑制。图6为表示通过本实施方式所取得的CBF映像的一例的图。
与此相对，一般的DSC-MRA时，原信号强度-时间曲线、浓度-时间曲线以及CBF映像分别为如图7至图9所示的那样。即，arteryTSC的浓度峰值与tissue TSC的浓度峰值相比增大，CBF映像的动脉区域与组织区域相比被强调了。
像这样即使在不并用动态成像的成像中，根据本实施方式也能够描出与动脉相比强调了组织的图像。
(动态成像)
在动态成像中，取得根据在造影前收集到的磁共振信号重建的图像(以下，称为基准图像)与根据在造影开始后收集到的磁共振信号重建的图像(以下，称为动态图像)之比而生成医用诊断用图像。即，将关于某像素(pixel)的基准图像以及动态图像的信号值分别设为S0以及S(t)时，将作为S0/S(t)求出的值作为医用诊断用图像的信号值。
但是，动态成像中，一般而言，只要用于对基准图像以及动态图像分别成像的成像条件不同即可。即，将基准图像以及动态图像分别作为造影前后的图像是动态成像的一例。众所周知，动态成像中的成像条件根据例如应用的成像序列、失相量、或对被检体施加的负荷等多种要素而改变。
一般使用了造影剂的MR信号强度(振幅)通过下面的式(1)来表示。
S＝K·A(T1，ΔR1*，TR)·A(T2，TE)·A(ΔR2*，TE)·A(Dflow，ADC，b)  (1)
其中，式(1)中包含的各参数分别如下定义。
K是增益等硬性决定的系数，在一项研究内视为相同。
T1缓和项：A(T1，TR)＝1-exp[-TR(1/T1+ΔR1*)]
ΔR1*：基于造影剂的T1缓和速度
T2缓和项：A(T2，TE)＝exp[-TE/T2]
基于造影剂的T2缓和项：A(ΔR2*，TE)＝exp[-ΔR2*TE]
ΔR2*：基于造影剂的T2缓和速度
基于Flow与扩散的项：A(Dflow，ADC，b)＝exp[-b*(Dflow+ADC)]
Dflow：由血流产生的相位方差等价分量
ADC：扩散系数
(a)使用阴性造影效果时
在DSC-MRI中，视为通过造影剂投放，T1是一定的且ΔR2*随着时间变化而变化。
理想而言，如果求造影后的信号强度S相对于造影前的信号强度S0的比，则用式(2)来表示。
S/S0＝A(ΔR2*，TE)＝exp[-ΔR2*TE]      (2)
据该式(2)可知，包含含有b值的项，而其他项被取消。
然而，这是忽视了SNR的情况，在噪声存在下情况有所不同。
当在基准图像中SNR充分降低至噪声水平程度时，即使造影剂流入也几乎不会发生信号降低。通常，由于大多数解析通过绝对值图像进行的，因此噪声的平均值不会在正的某值以下。
在式(1)中将b＞0，将A(Dflow，ADC，b)减小接近0，预先使血管信号降低到噪声水平程度，则用式(2)求出的AIF成为像顶端被压扁了那样的特性。如果血管部分的信号强度S在造影前后完全达到噪声水平的话，则式(2)的解为1。即，造影了的部分也成为与没造影的部分相同。此时，即使信号不完全降低至噪声水平，作为结果，如果AIF的大小被降低到组织中的time-intensity curve(TIC，时间-强度曲线)的大小程度，则在解析时难以受到血管的影响。另外，在绝对值图像上如果预先实施空间性、时间性平滑，则在噪声下曲线在时间方向上的分散也减少。
即，作为b＝0的一般动态成像中，如图10所示，在基准图像上，在血管与组织中S0为相同程度，与此相对，在动态图像中，血管与组织相比S大幅度降低。因此，对于诊断用图像中的S，血管与组织相比大幅度变小，在造影剂浓度变换图像中血管被强调。
然而，在本实施方式中，如图11所示，对于基准图像中的S0，血管与组织相比接近噪声水平，因此对于随着造影剂流入的S的下降率，血管与组织相比变小。因此，对于诊断用图像中的S，血管与组织相比变小，从而降低血管的影响。
而且，如图12所示如果将b设定成可以使血管S0减小到与噪声水平相同程度的适当值，就可以将诊断用图像中的血管的S0设为噪声水平程度。其结果，在诊断图像中，能够与不含组织周围的毛细血管的周围组织一样地描出血管，进而降低血管的影响。
在b＞0时，不能取得正确的AIF。但是能够取得血管的影响少的灌流(perfusion)解析结果。另外，如果降低到某周围组织程度的水平，即使实施平滑，与b＝0的灌流解析结果相比，对血管周围组织的影响变小。组织的血液量CBV较少，为3％左右，因此100％血液所占的血管信号较大，为组织的100/3＝33倍左右。因此，特别是在像DCE-MRI或DSC-MRI那样使用血管造影剂的灌流测定方法中，容易产生血管信号的影响，降低血管信号的效果在提高精度上是很重要的。另外，成像条件也是，缩小距阵尺寸、增大切片厚度而使部分容积效应大的情况是重要的。
(b)使用阳性造影效果时
造影后的S相对于造影前的S0的比通过下面的式(3)来表示。
S/S0＝A(ΔR2*，TE)＝1-exp[-TR(1/T1+ΔR1*)]  (3)
即，对于血管部分的降低，在与噪声的关系上，与使用阴性造影剂的情况相同的讨论成立。
但是，在算出ΔR1*时需要测定或假定T1，但是这也可以根据在投放造影剂前以2阶段以上对TR成像所得的图像算出。
对于造影前后的信号强度的差，如下面的式(4)所示，b的项A(Dflow，ADC，b)的作用作为对整体的信号的一定权重而进行施加。因此，如果b＞0，则对于血管部分，乘以某个恒定的且比组织部分小的缩放比例的值，血管信号降低。
S-S0＝K·A(T2，TE)·A(ΔR2*，TE)·A(Dflow，ADC，b)·[{1-exp[-TR/T1]}-{1-exp[-TR(1/T1+ΔR1*)]}]                 (4)
对于基于用于流动失相的MPG的b值，当要完全BB化时变大，因此，只要AIF达到与组织TIC相同程度的大小就足够了。
在DSC-MRI或DCE-MRI中，与组织的T2*或T2值大致相等的TE(如果是脑部，TE＝50-80ms)是合适的，因此在TE间插入使血管信号降低得到的b＜10[sec/nm2]左右的小的MPG的时间是充裕的。
(c)对fMRI的应用
作为脑功能检查方法大多数情况下使用利用了BOLD的fMRI，使用以下的效果，即、通过随着由神经刺激引起的血流增加，脱氧血红蛋白(deoxy-hemoglobin)相对于氧合血红蛋白(oxy-hemoglobin)相对减少，信号强度比安静时提高。另一方面，通过在伴随血流速增加的薄块(slab)内纵磁化饱和的血液比安静时更多地流入也成为高信号。但是，流动效应在薄块内不均匀且依赖于流速等，未必是因为与血流增加相关，因此形成伪影，本来只想观察BOLD效果，而且是其的毛细血管水平的变化。
在本实施方式中，使用赋予了使比较粗的血管信号降低至噪声水平程度的b值(b＝2-10左右)的GRE体系序列。
由此，在能够进行三维(3D)收集的3D-GRE中也可以应用，对于时间分辨率也是，由于刺激周期为1分钟左右，所以每次30sec左右的收集时间也是充足的。在b＜10的程度下，由于存在毛细血管内的流动，因此不需要是高SNR。这样，如果根据利用这种本实施方式的序列收集到的磁共振信号实施fMRI，能够只对于BLOD效果，作为毛细血管水平的变化进行观察。非专利文献1中所公开的技术中，由于是不看血流而看细胞变化的DW-fMRI，所以b值大，毛细血管内的信号也大致变成0。
另外，在对fMRI应用本实施方式的序列时，GRE-EPI、spin warpGRE，gradient-and spin-echo(GRASE)、或asymmetric spinecho(ASE)等序列也能够应用。
(d)对每个时相组合不同的失相强度
共用于以上情况，在动态(dynamic)的全部时间相上，可以使用相同的MPG(b)强度的序列，如果也想使用AIF时也可以使用对每个时相区(block)追加了不同强度的MPG(b)的序列。例如，在每个初始时相中将基准图像设为b＝2-10，在以后的时相中设为b＝0，分别进行成像。即，对基准图像成像时的b值与对动态图像进行成像时的b值相比变大。并且，生成用于使用阈值等从基准图像排除血管的蒙片，也可以在以后的时相图像的区域选择中使用该蒙片。此时，如果使用重相(rephase)与失相(dephase)之差，则血管位置的鉴别进而变得容易。在以往的方法中，如果不是在造影后、且不检查全时相的图像的话就不清楚血管位置，但是如果是本实施方式，则可以在造影前或在活性化实验前鉴别血管位置，因此能够进行不存在血管伪影的实时的与基准图像相比的变化的动态图像显示或即时的解析映射(mapping)图像显示。另外，这些方法也可以在使用AIF的解析法中应用。进而，取得基准图像S0时，如果预先不仅追加b＞0还追加b＝0则可以与两者对应。
(e)重相的应用
每个时相中的跳动等变动带来血管信号的变动，会阻碍正确的计测。此时，上述失相应用例子也有利于血管的稳定化，但是相反，如果使用重相，则能够对每个时相测定出血流变动的影响少的稳定的AIF。相对于失相时的上述(a)至(d)的优点，血管伪影的排除是困难的，但是AIF的稳定取得成为不同的优点。
(f)基于图像处理的稳定的血管除去法
在上述(a)的说明中，叙述了通过以血管的信号值达到噪声程度的方式适当地设定b值，能够与周围组织一样描出血管的情况。然而，设定那样的适当的b值是不容易的。
因此，通过如图13所示那样对基准图像以及动态图像实施阈值处理，强制性地使血管的S0以及S(t)与阈值Sth一致。即，对基准图像、时间方向图像的全部图像以一定值赋予相当于原图像的血管部分的像素。具体而言，预先设定某阈值Sth，对基准图像以及动态图像的全部进行以下的式(5)以及式(6)的处理。
S0＜Sth 然后S0＝Sth       (5)
S＜Sth  然后S＝Sth        (6)
并且，如图13所示，求出阈值处理后的基准图像与动态图像之比而生成诊断用图像。
另外，阈值Sth可以根据相当于例如血管部分的像素的信号值来求出。即，例如，将基准图像以及动态图像的信号值S0、S(t)下的血管的信号值中的最大值作为阈值Sth。如果是DSC-MRI即使在时相方向不检索，只在为了造影而保证了在时相方向上为最大值的基准图像中预先抽出血管，也可以作为其蒙片内的血管部分的最大值求出阈值Sth。或者，作为基准图像以及动态图像的信号值S0、S(t)下的血管的信号值，也可以使用比一般取得的信号值更大地设定的固定值作为阈值Sth。
由此，即使血管的信号水平降低至噪声水平，在血管部分中S/S0＝1、即在表示造影剂浓度的C＝ln(S0/S)/TE中为相同图像值，也能够完全抑制血管。另外，阈值Sth能够与为了排除空气而使用的蒙片生成的阈值兼用。由于空气的信号强度比血管小，因此，通过上述处理必然排除空气。
另外，S0th、Sth以下的像素也可以不进行原信号强度映像的加工(置换为S0＝Sth、S＝Sth的处理)，而直接将浓度映像C的相符部分设为0。
(g)不同b强度(失相(b＞0)以及b＝0/重相)的组合
只在造影前利用b＝1-10[sec/nm2]左右的失相序列收集磁共振信号(最低一个点)，造影后利用b＝0或重相序列收集磁共振信号。并且，分别根据这样收集到的磁共振信号，分别重建2个基准图像S0d、S0。另外，通过b＝1-10[sec/nm2]左右的失相序列对多个点收集磁共振信号时，取其平均。进而，在non-AIF法的情况与AIF法的情况下，分别像以下那样进行处理。
(non-AIF法时)
针对全部像素分别判定下面的关系是否成立。
S0d/S(t)≤1
并且，对于该关系成立的像素设C(t)＝0，对于不成立的像素设C(t)＝ln[S0/S(t)]/TE。
由此，由于血管因失相效应S0d＝S0＜1，因此浓度C(t)为0，由于组织的S0d与S0不同，因此可以算出不为0的C(t)。
(AIF法时)
想要求出AIF时，不使用失相的基准图像，只要使用通常(标准、重相)的基准图像无条件地通过通常的下式对全部像素求出浓度即可。
C(t)＝ln[S0/S(t)]/TE
对于AIF、Ca(t)，在non-AIF法的情况以及AIF法的情况中是相同的，可使用峰值或面积、顶峰到达时间等从全部像素的C(t)中抽出。
图14以及图15为分别表示该情况中的原信号强度-时间曲线以及浓度-时间曲线的图。
这样，根据本实施方式，能够取得将包含血液比较低速流动的毛细血管的组织相对于血液高速流动的动脉等粗的血管强调地描出的图像，能够取得对观察毛细血管等有用的图像。
另外，通过与流动重相(flow-rephase)的组合，在使用了血管造影剂的灌流(perfusion)测定或诊断中能够降低成为阻碍的血管信号，通过与non-AIF法(例如，将与C(t)的最大倾斜的基准组织的比作为CBF指标，将与C(t)的曲线下方面积的基准组织的比作为CBV指标的方法)组合，即使不使用AIF也可以实现CBF，cerebral bloodvolume(CBV)，mean transit time(MTT)的半定量化。
另外，通过与流动重相的组合，改善所生成的AIF的稳定性，在DCE-MRI中能够改善基线(低浓度)部分中的SNR，在DSC-MRI中能够改善饱和效应引起的高浓度部分中的精度或稳定性。
另外，由于能够在注入造影剂以前的早期鉴别出血管部分，因此，也能够使鉴别出血管部分的灌流(perfusion)解析处理高速化。
进而，能够实现支配抑制了血管伪影的BOLD信号的fMRI。
本实施方式能够实施以下种种变形。
(1)即使是为了描出血液以外的淋巴液、脑脊髓液等流体的目的，也能够应用本发明。
可以代替失相序列，组合任意序列与利用了相位信息的后处理来应用。另外作为利用了相位信息的后处理，考虑与利用余弦滤波的相位信息对应的强调处理。
(2)也能够应用余弦滤波(COS-filter)处理。
如图16所示，背景相位校正后的血管信号Sv的实部(real-part)Re[Sv]与振幅|Sv|相比总是向负侧变小。因此，实部中的血管与背景之间的对比度(vessel-to-background contrast)Cvb(real)与振幅中的血管与背景之间的对比度Cvb(mag)相比，总是增大。其中，Cvb(mag)＝|Sv|-|Sb|。
图17为表示余弦滤波特性的图。
用于黑血法的余弦滤波是针对用于将BB图像中的血管向负方向强调的相位的滤波，将针对相位的余弦函数(cosine function)作为基准。
余弦滤波HB通过下式来表示。
HB＝2×(Mn-0.5)
在此，余弦滤波具有非对称型与对称型。如果Im[Scor]＜0或Φcor＜0的任何一个条件成立的话可以通过下式来确定非对称型时的M，如果上述条件都不成立则将其设为0。
M＝{cos(Φcor)+1}/2
对称型时的M可以无条件地通过下式来确定。
M＝{cos(Φcor)+1}/2
另外，n为强调因子。强调因子n为0以上的值，n越大则强调越强。
余弦滤波处理是使用了实部信息的信号强调处理。即，余弦滤波处理是通过基于实部信息通过上述特性使信号值变化而相对于背景信号强调血管信号的处理。BB图像中的余弦滤波处理是以使血管信号的信号值相对于背景信号向负方向增大的方式进行的。
相位校正后进行n＝1的余弦滤波处理与获取实部是等效的。进行n＞1的余弦滤波处理与向180°方向强调相位而获取实部是等效的。另外，在余弦滤波之前进行的使用了零差滤波(homodyne filter)的背景相位校正中，大多数情况下过小评价(underestimate)实际相位，进而，Cvb(real)向负侧变得最大，由于校正后的血管相位Φv是180°的情况，因此如果进行n＞1的余弦滤波处理就可以等价地使Φv接近180°。
如图18以及图19为表示进行余弦滤波处理时的摄像步骤的一例的流程图。
图18所示的摄像步骤表示在浓度变换前进行基于余弦滤波的血管遮蔽的例子的图。此时，依次为基准图像的收集(步骤Sa1)、动态图像的收集(步骤Sa2)、背景相位的校正(步骤Sa3)、余弦滤波处理(步骤Sa4)、阈值处理(步骤Sa5)、造影浓度变换(步骤Sa6)。
图19所示的摄像步骤表示在浓度变换后进行基于余弦滤波的血管遮蔽的例子。此时，依次为基准图像的收集(步骤Sb1)、动态图像的收集(步骤Sb2)、造影浓度变换(步骤Sb3)、背景相位的校正(步骤Sb4)、余弦滤波处理(步骤Sb5)、阈值处理(步骤Sb6)。
其次，针对使用阴性造影剂的DSC-MRI中的余弦滤波处理的效果进行说明。另外，以通过复数平面中的说明容易理解的方式通过实部信号与振幅之间的关系进行说明。
图20以及图21为表示在背景相位校正以及余弦滤波处理后的造影前后的血管信号的一例的图。另外，图20表示造影前后的血管信号彼此为相同极性的情况，图21表示相反极性的情况。
如图20、图21所示，当将造影前后的复数信号分别设为S0、S时，不论MPG是否存在，如果是同一条件，造影剂的磁化率效应引起的相移进一步增大，因此，即使在余弦滤波处理前、即在振幅中有|S0|＞|S|，也会有Re[S0]＜|S0|、Re[S]＜|S|。
图22为用血管与背景的空间分布图表示余弦滤波处理的效果的图。
在图22中，关于血管部分，使用虚线来表示基于振幅信息的信号值，使用实线来表示余弦滤波处理后的信号值。点划线表示余弦滤波处理后的基准图像与没有实施余弦滤波处理的造影图像的信号比。另外，图22的左侧是没有实施阈值处理的状态的分布图，右侧为实施了阈值处理后的分布图。此外，阈值处理中使用噪声水平程度的阈值。
从图22也可以看出，通过余弦滤波处理，基准图像中的血管的信号值变小。并且，与血管相关的信号的相位超过90度时，余弦滤波处理后的基准图像中的血管的信号值为负侧。像这样，通过应用余弦滤波处理可以增大血管与组织之间的信号值的差，因此容易基于阈值处理区别血管与组织。并且，在阈值处理中，作为阈值来，能够使用噪声水平或0。
即使将浓度变换后的ΔR2*也限定于由于ln为相同符号且＞0而能够计算出信号的情况，也有ln[Re[S0]/|S|]＜ln[|S0|/|S|]。即，与在造影前后都不应用余弦滤波处理的情况相比，只在造影前应用余弦滤波处理的情况的信号值变小，血管抑制效果增大。另外，ln＜0时，不能计算出信号值，但是此时，鉴于造影前后的相位差本身较大这一情况，是血管的概率高，因此可以视为血管而将浓度设为0等。如果与振幅一样MPG充分变大而可以降低到噪声水平程度，则可以省略使用了余弦滤波处理时的阈值处理。但是，为了发挥即使不提供大的MPG也能够降低的余弦滤波处理的优点，希望将相位强调到近180°而作为负信号进行噪声水平程度的阈值处理。
进而，即便在组织中，相移也会根据余弦滤波的强度或造影剂浓度变大或变小。因此，想要使余弦滤波处理对组织完全没有影响时，在通过阈值处理作为血管被抽出的范围以外最好原样地使用振幅图像。为此，将施加在振幅图像上的图像蒙片作为1，只将血管作为0即可，但是由于对于余弦滤波HB的增益，能够忽视相移的背景组织部能够大致为1，因此也可以容易地生成那样的血管抑制蒙片。余弦滤波本身的非线性大，针对0附近的相位，能够相对于稍微的变动而将增益维持在1附近。但是，如果需要的话，也可以将BB用余弦滤波本身做成使用相位的阈值阶段性改变后的形状。例如，通过下式来设定。
Hb＝1:|Φ|＜Φth1、＝0.5(1-Mn):Φth1＜|Φ|＜Φth2、＝0:Φth2＜|Φ|
其中，Φth1是划分背景组织与血管的阈值，Φth2是血管的上限以上且在静止组织中划分磁化率大的部分的阈值。进而也可以平滑连结相位方向断坡。
组合余弦滤波时，至少在基准图像收集中，既可以应用MPG，也可以不应用MPG。不应用MPG时，其优点是可以使用通常的脉冲序列。但是此时，由于在基准图像中也没有流入造影剂，因此由于不存在血管上的相移而不能期待进行余弦滤波处理的效果。因此，此时，通过对造影后的图像应用余弦滤波处理，生成抽出了血管的蒙片(背景组织为1，血管为0)，使用它来抑制血管部分。在余弦滤波处理中，即使在造影后，如果背景组织中的相移充分小的话，能够将余弦滤波设为1，即像振幅自身那样。因此，如果使用对动脉相限定了时间范围的造影剂浓度在时间方向上的最大值的映像、即峰高(peak height)的映像，由于只在动脉部分中不向组织移行的概率高，因此能够减小对组织的影响。
由于余弦滤波处理的应用阶段、即对振幅图像中的空间蒙片的应用阶段的目的在于，最终在浓度图像上抑制血管，因此只要注意实施方式中举出的负值或0的除法等的例外处理，也可以在原图像上，也可以在浓度变换后的图像上(如图19所示的步骤时)。此时，其优点是，不论造影剂的阴性/阳性，都能够使用共用的方法处理。
另外，余弦滤波处理不仅应用于阴性造影剂，也可同样应用于阳性造影剂的情况或f-MRI。
但是，与阳性情况不同的是，对原图像应用余弦滤波处理时，不对基准图像，而至少对造影后的数据实施余弦滤波处理。这与通过施加MPG抑制血管的情况相同。
对以上余弦滤波处理的效果概括如下：
·能够提高由于血管与背景的对比度的提高产生的抑制效果。
·能够省略由MPG产生的失相。
·希望只根据基准图像生成除血管以外的图像蒙片，即使在造影后也将该图像蒙片应用于原图像或浓度图像中。
在以上中，为了使BB图像中的血管与组织产生对比度或者为了强调其对比度而应用余弦滤波处理。然而，为了使WB图像中的血管与组织产生对比度，或者为了强调其对比度，也可以应用余弦滤波处理。此时所应用的余弦滤波与图17所示的特性相反，具有相位越接近0则增益越小的特性。
本领域技术人员容易想到其它优点和变更方式。因此，本发明就其更宽的方面而言不限于这里示出和说明的具体细节和代表性的实施方式。因此，在不背离由所附的权利要求书以及其等同物限定的总的发明构思的精神和范围的情况下，可以进行各种修改。

Claims (17)

1.一种磁共振成像装置，其特征在于，包括：
收集单元，利用第1序列收集在投放造影剂前在成像区域内产生的磁共振信号，并且利用第2序列收集在投放上述造影剂后在在上述成像区域内产生的磁共振信号，其中，上述第1序列是使射频激励后的磁化失相从而使关于在第1流速范围内流动的流体的第1磁共振信号分量与关于在比上述第1流速范围低的第2流速范围流动的上述流体的第2磁共振信号分量相比信号大幅度降低的序列，上述第2序列是设定为与上述造影剂的浓度对应的大小的序列；
重建单元，根据分别利用上述第1以及第2序列收集到的上述磁共振信号，分别重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号的大小而表示的第1图像以及第2图像；以及
生成单元，根据上述第1图像以及上述第2图像生成表示了在上述造影剂投放后相对于投放前的上述流体的变化程度的第3图像。
2.根据权利要求1所述的磁共振成像装置，其特征在于：
上述重建单元重建上述第1图像以及上述第2图像作为表示与在上述成像区域内产生的磁共振信号对应的信号值的空间分布的图像；
上述生成单元使用作为与上述造影剂投放前的上述第1磁共振信号分量的大小对应的信号值和与上述第2磁共振信号分量的大小对应的信号值之间的中间值所设定的阈值对上述第1图像以及上述第2图像中包含的信号值进行阈值处理，并根据进行该阈值处理后的上述第1图像以及上述第2图像生成上述第3图像。
3.根据权利要求1所述的磁共振成像装置，其特征在于：
上述收集单元除了利用上述第1序列之外，还利用使上述第1磁共振信号分量相对于上述第2磁共振信号分量的信号降低量比上述第1序列小的第3序列收集在投放上述造影剂前在成像区域内产生的磁共振信号；
上述重建单元根据利用上述第3序列收集到的上述磁共振信号重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号的大小而表示的第4图像；
上述生成单元根据上述第1图像、上述第2图像以及上述第4图像生成上述第3图像。
4.根据权利要求1所述的磁共振成像装置，其特征在于：
上述收集单元利用动作探测梯度脉冲使射频激励后的磁化失相，作为b值使用使在上述第1流速范围内流动的上述流体与在上述第2流速范围流动的上述流体相比信号大幅度降低的值。
5.根据权利要求1所述的磁共振成像装置，其特征在于：
上述收集单元使关于以上述第1流速范围内的流速流动的上述流体的磁共振信号分量的信号强度降低到关于上述流体以外的流体的磁共振信号分量的信号强度相同程度。
6.根据权利要求1所述的磁共振成像装置，其特征在于：
上述重建单元以将产生了更大的上述磁共振信号的上述流体强调地表示的方式重建上述第1图像或上述第2图像。
7.根据权利要求1所述的磁共振成像装置，其特征在于：
上述第1序列是通过包含上述倾斜磁场失相脉冲而使血管信号降低从而较黑地描出血管的黑血法的序列。
8.一种磁共振成像装置，其特征在于，包括：
收集单元，利用为了使在第1流速范围内流动的流体与在比上述第1流速范围低的第2流速范围流动的上述流体相比信号大幅度降低而使射频激励后的磁化失相的失相序列收集在成像区域内产生的磁共振信号；以及
重建单元，根据利用上述收集单元收集到的上述磁共振信号重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号的大小而表示的图像。
9.根据权利要求8所述的磁共振成像装置，其特征在于：
上述收集单元利用动作探测梯度脉冲使射频激励后的磁化失相，作为b值使用使在上述第1流速范围内流动的上述流体与在上述第2流速范围内流动的上述流体相比信号大幅度降低的值。
10.根据权利要求8所述的磁共振成像装置，其特征在于：
上述收集单元使关于以上述第1流速范围内的流速流动的上述流体的磁共振信号分量的信号强度降低到与关于上述流体以外的流体的磁共振信号分量的信号强度相同的程度。
11.根据权利要求8所述的磁共振成像装置，其特征在于：
上述重建单元以将产生了更大的上述磁共振信号的上述流体强调地表示的方式重建上述第1图像、上述第2图像或上述图像。
12.一种磁共振成像装置，其特征在于，包括：
收集单元，使用通过包含倾斜磁场失相脉冲而使血管信号降低从而较黑地描出血管的黑血法的脉冲序列对投放阴性造影剂前的被检体进行成像，使用规定的成像方法对投放上述阴性造影剂后的被检体进行动态成像；
重建单元，根据对上述造影剂投放前的被检体进行成像所取得的磁共振信号重建基准图像，并且根据对上述造影剂投放后的被检体进行动态成像所取得的磁共振信号重建多时相的造影图像；以及
生成单元，根据上述基准图像以及上述造影图像生成诊断用图像。
13.一种磁共振成像装置，其特征在于，包括：
收集单元，分别以第1条件和第2条件收集在成像区域内产生的磁共振信号，其中上述第1条件是使关于在第1流速范围内流动的流体的第1磁共振信号分量的信号水平相对于关于在比上述第1流速范围低的第2流速范围流动的上述流体的第2磁共振信号分量的信号水平在规定方向上不同的条件，上述第2条件是使上述第1磁共振信号分量的信号水平相对于上述第2磁共振信号分量的信号水平在规定方向上不同，并且上述第1以及第2磁共振信号分量的信号水平都相对于上述第1条件在上述规定方向上不同，而且对于其变化量，上述第2磁共振信号分量相对于上述第1磁共振信号分量增大的条件；
重建单元，根据以上述第1以及第2条件分别收集到的上述磁共振信号，分别重建将上述流体在上述成像区域中的空间分布反映上述磁共振信号的大小而表示的第1图像以及第2图像；以及
生成单元，根据上述第1图像以及上述第2图像生成表示了上述第1条件以及上述第2条件下的信号水平的变化程度的第3图像。
14.根据权利要求13所述的磁共振成像装置，其特征在于：
上述第1流速范围被设定为包含动脉以及静脉中的血液的流速，上述第2流速范围被设定为包含毛细血管中的血液的流速。
15.根据权利要求13所述的磁共振成像装置，其特征在于：
上述收集单元将第1以及第2条件设为使上述第2磁共振信号分量的信号水平相对于上述第1磁共振信号分量的信号水平降低的条件。
16.根据权利要求13所述的磁共振成像装置，其特征在于：
上述收集单元将第1以及第2条件设为使上述第2磁共振信号分量的信号水平相对于上述第1磁共振信号分量的信号水平增加的条件。
17.根据权利要求13所述的磁共振成像装置，其特征在于：
上述收集单元进行基于实部信息使上述第1磁共振信号分量的信号水平与上述第2磁共振信号分量的信号水平之差增大的余弦滤波处理。
CN2010101545337A 2009-04-03 2010-04-02 磁共振成像装置 Expired - Fee Related CN101856228B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2009-091303 2009-04-03
JP2009091303 2009-04-03

Publications (2)

Publication Number Publication Date
CN101856228A true CN101856228A (zh) 2010-10-13
CN101856228B CN101856228B (zh) 2013-04-24

Family

ID=42825666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101545337A Expired - Fee Related CN101856228B (zh) 2009-04-03 2010-04-02 磁共振成像装置

Country Status (3)

Country Link
US (1) US8344730B2 (zh)
JP (1) JP5624346B2 (zh)
CN (1) CN101856228B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102008292A (zh) * 2010-10-29 2011-04-13 中国人民解放军第二军医大学 术中脊髓监测系统
CN103208941A (zh) * 2012-01-13 2013-07-17 株式会社东芝 电源电路以及磁共振成像装置
CN103502832A (zh) * 2011-04-21 2014-01-08 皇家飞利浦有限公司 使用具有流动敏化梯度的受激回波脉冲序列的黑血mri
CN103733083A (zh) * 2011-08-17 2014-04-16 皇家飞利浦有限公司 减小磁共振成像期间预定体积中的射频发射场
CN104583795A (zh) * 2012-08-27 2015-04-29 皇家飞利浦有限公司 用于动态对比增强mri的适应性键孔压缩
CN104644172A (zh) * 2013-11-25 2015-05-27 深圳先进技术研究院 时间反转稳态进动快速成像序列中黑血的控制方法和系统
CN111025210A (zh) * 2019-12-20 2020-04-17 东软医疗系统股份有限公司 磁共振成像方法、装置、电子设备、存储介质

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5481117B2 (ja) * 2008-09-29 2014-04-23 株式会社東芝 磁気共鳴診断装置およびその制御方法
DE102010001552B4 (de) * 2010-02-03 2012-03-08 Siemens Aktiengesellschaft Verfahren zum Bestimmen eines Hintergrundphasenverlaufs, Magnetresonanzvorrichtung und Computerprogramm
WO2013138794A1 (en) * 2012-03-16 2013-09-19 Hologic, Inc. Mri display output reflecting contrast agent concentration as a function of time
WO2013144828A1 (en) * 2012-03-30 2013-10-03 Koninklijke Philips N.V. Simultaneous high spatial low temporal resolution magnetic resonance (mr) sequence for dynamic contrast enhanced (dce) magnetic resonance imaging (mri)
JP6266325B2 (ja) * 2013-12-06 2018-01-24 東芝メディカルシステムズ株式会社 画像解析装置、画像解析プログラム、及び、磁気共鳴イメージング装置
US10814019B2 (en) * 2014-06-30 2020-10-27 University Of Washington MRI signal suppression agents, compositions, and methods
DE102014223388B4 (de) * 2014-11-17 2016-07-28 Siemens Healthcare Gmbh Beschleunigung einer MR-Fingerprinting-Messung
KR101797141B1 (ko) * 2016-04-18 2017-11-13 성균관대학교산학협력단 자기 공명 영상 처리 장치 및 그 영상 처리 방법
US10297038B2 (en) * 2017-03-30 2019-05-21 Siemens Healthcare Gmbh Determination and presentation of flow transit curves
US10888246B2 (en) * 2017-07-18 2021-01-12 Synthetic Mr Ab Method and system for generating a contrast enhancement map
CN110391016B (zh) * 2019-07-10 2021-11-09 浙江大学 一种动态对比增强磁共振图像的分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6580937B2 (en) * 2000-12-30 2003-06-17 Ge Medical Systems Global Technology Co., Llc Method for optimal imaging of the peripheral vasculature emphasizing distal arterial visualization in a multi-station examination
CN101264018A (zh) * 2007-03-15 2008-09-17 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法
CN101330872A (zh) * 2005-12-15 2008-12-24 皇家飞利浦电子股份有限公司 一种在血管造影术中用于确定脉动程度的方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4663772A (en) * 1985-09-30 1987-05-05 Picker International, Inc. Bone mineral analysis phantom
JP3146033B2 (ja) * 1991-11-05 2001-03-12 株式会社東芝 磁気共鳴イメージング装置
JP4379559B2 (ja) * 1999-04-30 2009-12-09 株式会社日立メディコ 磁気共鳴イメージング装置
US6618605B1 (en) * 1999-09-08 2003-09-09 General Electric Company Method and apparatus for MR perfusion image acquisition using a notched RF saturation pulse
US6340887B1 (en) * 1999-09-21 2002-01-22 Picker International Inc. Multiple contrast FSE approach to black blood angiography with redundant and supplementary vascular information
JP4368988B2 (ja) * 1999-10-15 2009-11-18 Ｇｅヘルスケア・ジャパン株式会社 Ｍｒｉ装置
JP2001252262A (ja) * 2000-03-13 2001-09-18 Hitachi Medical Corp Ｍｒｉ差分画像処理方法及びｍｒｉ装置
JP4693227B2 (ja) * 2000-11-14 2011-06-01 株式会社東芝 磁気共鳴イメージング装置及び磁気共鳴イメージングのスキャン同期方法
US6804546B1 (en) * 2001-04-20 2004-10-12 Koninklijke Philips Electronics, N.V. Multiple contrast echo-planar imaging for contrast-enhanced imaging
US6559644B2 (en) * 2001-05-30 2003-05-06 Insightec - Txsonics Ltd. MRI-based temperature mapping with error compensation
US6980001B2 (en) * 2002-05-20 2005-12-27 The University Of Sheffield At Western Bank Methods & apparatus for magnetic resonance imaging
JP3434816B2 (ja) * 2002-11-06 2003-08-11 株式会社東芝 Ｍｒｉ装置
JP4537681B2 (ja) * 2003-09-24 2010-09-01 株式会社東芝 血流解析装置
JP2006158512A (ja) * 2004-12-03 2006-06-22 Hitachi Medical Corp 磁気共鳴イメージング装置
US9201129B2 (en) * 2006-09-13 2015-12-01 Kabushiki Kaisha Toshiba Magnetic-resonance image diagnostic apparatus and method of controlling the same
JP5288745B2 (ja) * 2006-09-13 2013-09-11 株式会社東芝 磁気共鳴イメージング装置
WO2008109783A2 (en) * 2007-03-06 2008-09-12 The Regents Of The University Of California Detecting spin perturbations using magnetic resonance imaging
JP5414335B2 (ja) * 2008-07-24 2014-02-12 株式会社東芝 磁気共鳴イメージング装置
AU2010232430B2 (en) * 2009-04-03 2014-08-28 The United States Of America, As Represented By The Secretary, Department Of Health And Human Services Magnetic microstructures for magnetic resonance imaging

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6580937B2 (en) * 2000-12-30 2003-06-17 Ge Medical Systems Global Technology Co., Llc Method for optimal imaging of the peripheral vasculature emphasizing distal arterial visualization in a multi-station examination
CN101330872A (zh) * 2005-12-15 2008-12-24 皇家飞利浦电子股份有限公司 一种在血管造影术中用于确定脉动程度的方法
CN101264018A (zh) * 2007-03-15 2008-09-17 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102008292B (zh) * 2010-10-29 2012-10-17 中国人民解放军第二军医大学 术中脊髓监测系统
CN102008292A (zh) * 2010-10-29 2011-04-13 中国人民解放军第二军医大学 术中脊髓监测系统
CN103502832B (zh) * 2011-04-21 2016-03-02 皇家飞利浦有限公司 使用具有流动敏化梯度的受激回波脉冲序列的黑血mri
CN103502832A (zh) * 2011-04-21 2014-01-08 皇家飞利浦有限公司 使用具有流动敏化梯度的受激回波脉冲序列的黑血mri
CN103733083B (zh) * 2011-08-17 2016-10-26 皇家飞利浦有限公司 采集成像区域内受检者的磁共振数据的系统及操作方法
CN103733083A (zh) * 2011-08-17 2014-04-16 皇家飞利浦有限公司 减小磁共振成像期间预定体积中的射频发射场
CN103208941B (zh) * 2012-01-13 2015-06-17 株式会社东芝 电源电路以及磁共振成像装置
CN103208941A (zh) * 2012-01-13 2013-07-17 株式会社东芝 电源电路以及磁共振成像装置
CN104583795A (zh) * 2012-08-27 2015-04-29 皇家飞利浦有限公司 用于动态对比增强mri的适应性键孔压缩
CN104583795B (zh) * 2012-08-27 2018-02-16 皇家飞利浦有限公司 具有适应性键孔压缩的动态对比增强mri系统和方法
CN104644172A (zh) * 2013-11-25 2015-05-27 深圳先进技术研究院 时间反转稳态进动快速成像序列中黑血的控制方法和系统
CN104644172B (zh) * 2013-11-25 2017-11-03 深圳先进技术研究院 时间反转稳态进动快速成像序列中黑血的控制方法和系统
CN111025210A (zh) * 2019-12-20 2020-04-17 东软医疗系统股份有限公司 磁共振成像方法、装置、电子设备、存储介质
CN111025210B (zh) * 2019-12-20 2022-02-01 东软医疗系统股份有限公司 磁共振成像方法、装置、电子设备、存储介质

Also Published As

Publication number Publication date
JP5624346B2 (ja) 2014-11-12
US8344730B2 (en) 2013-01-01
JP2010253256A (ja) 2010-11-11
CN101856228B (zh) 2013-04-24
US20100253342A1 (en) 2010-10-07

Similar Documents

Publication Publication Date Title
CN101856228B (zh) 磁共振成像装置
Kim et al. How we perform delayed enhancement imaging: HOW I DO…
US9002430B2 (en) System and method for combined time-resolved magnetic resonance angiography and perfusion imaging
CN101632584B (zh) 磁共振成像装置
US6381486B1 (en) Magnetic resonance angiography with vessel segmentation
US8958866B2 (en) System and method to analyze blood parameters using magnetic resonance imaging
EP1269212B1 (en) Magnetic resonance angiography with automated vessel segmentation
CN101354433B (zh) 记录和处理时间上相继的图像数据组的序列的方法及设备
US6556856B1 (en) Dual resolution acquisition of magnetic resonance angiography data with vessel segmentation
US10765349B2 (en) Magnetic resonance imaging device and method for calculating oxygen extraction fractions
US8332010B2 (en) Method for non-contrast enhanced magnetic resonance angiography
US9063207B2 (en) Magnetic resonance imaging apparatus and blood vessel image capturing method
US8099149B2 (en) MRI method for quantification of cerebral perfusion
CN101143094A (zh) 磁共振图像诊断装置及其控制方法
JP2002533144A (ja) コントラスト増幅型磁気共鳴血管造影法におけるコントラストの検出および誘導される再構築
JP6533571B2 (ja) 磁気共鳴イメージング装置及び画像作成方法
CN102266225A (zh) 磁共振成像装置
JP6708525B2 (ja) 磁気共鳴測定装置および画像処理方法
JPH05123312A (ja) 磁気共鳴イメージング装置
US20130088227A1 (en) Magnetic resonance-based method and system for determination of oxygen saturation in flowing blood
Rehwald et al. Cardiovascular MRI: its current and future use in clinical practice
Shao et al. Ultra-high field brain perfusion MRI
WO2001075469A1 (en) Magnetic resonance angiography with automated vessel segmentation
Cao Contrast enhanced MRI for the measurement of dynamic signal changes in the CSF and cerebral lymphatic vessels
WO2024081779A2 (en) Imaging slow flow dynamics of cerebrospinal fluid using magnetic resonance imaging

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160721

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130424

Termination date: 20190402