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

磁共振成像装置 Download PDF

Info

Publication number
CN104582567B
CN104582567B CN201380043226.3A CN201380043226A CN104582567B CN 104582567 B CN104582567 B CN 104582567B CN 201380043226 A CN201380043226 A CN 201380043226A CN 104582567 B CN104582567 B CN 104582567B
Authority
CN
China
Prior art keywords
image
control unit
data
marker
perfusion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201380043226.3A
Other languages
English (en)
Other versions
CN104582567A (zh
Inventor
周相之
宫崎美津惠
星野勉
Original Assignee
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 Medical Systems Corp filed Critical Toshiba Medical Systems Corp
Publication of CN104582567A publication Critical patent/CN104582567A/zh
Application granted granted Critical
Publication of CN104582567B publication Critical patent/CN104582567B/zh
Active 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/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0044—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
    • 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/0263—Measuring blood flow using NMR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • 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/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

Abstract

实施方式所涉及的磁共振成像装置具备序列控制部和显示控制部。上述序列控制部关于长度不同的多个上述规定时间，执行在对标记区域施加了标记脉冲之后经过规定时间后收集心脏的切片的数据的数据收集。上述显示控制部将关于多个上述规定时间收集到的多个上述数据的分析结果显示于显示部，以使得能够识别上述规定时间与上述切片的关系。

Description

磁共振成像装置
技术领域
本发明的实施方式涉及磁共振成像装置。
背景技术
以往，在基于磁共振成像的心脏的摄像法之一中，存在被称为心肌灌流(Perfusion、灌流)的方法。心肌灌流是通过高速地对T1强调图像进行连续摄像，从而作为时间变化观察基于钆(Gd(gadolinium))类的造影剂的心肌的染色，对向心肌的血液供给进行评估的方法。
现有技术文献
非专利文献
非专利文献1：Kim,et al.,“Relationship of MRI DelayedContrastEnhancement to Irreversible Injury,Infarct Age,and ContractileFunction,”Circulation Journal of the American Heart Association,pages1991-2002(November9,1999)
发明内容
本发明要解决的问题在于，提供一种能够恰当地显示由非造影得到的心肌灌流图像的分析结果的磁共振成像装置。
实施方式所涉及的磁共振成像装置具备序列控制部和显示控制部。上述序列控制部关于长度不同的多个上述规定时间执行从对标记区域施加标记脉冲开始经过了规定时间之后收集心脏的切片的数据的数据收集。上述显示控制部将关于多个上述规定时间收集到的多个上述数据的分析结果以使得能够识别上述规定时间与上述切片的关系的方式显示于显示部。根据上述的磁共振成像装置能够恰当地显示由非造影得到的心肌灌流图像的分析结果。
附图说明
图1是构成为提供非造影动态MRI心肌灌流的分析以及视觉化的、示例的实施方式所涉及的MRI系统的高水准概略框图。
图2A是表示本实施方式所涉及的MRI系统所使用的、标记添加、标记删除MRI数据收集序列的第1个示例的概略描绘图。
图2B是表示本实施方式所涉及的MRI系统所使用的、标记添加、标记删除MRI数据收集序列的第2个示例的概略描绘图。
图3是表示本实施方式中的标记脉冲的施加区域的图。
图4是表示在本实施方式中，为了生成标记添加以及标记删除的图像的强度图像而使用复数计算的、标记添加以及标记删除的图像的像素单位的减法的图。
图5A是表示本实施方式中的不同的BBTI间的位置对准前的左心室心肌的强度图像的图。
图5B是表示本实施方式中的不同的BBTI间的位置对准后的左心室心肌的强度图像的图。
图6A是表示本实施方式中的、对标记添加或者标记删除的图像施加的分割的图。
图6B是表示本实施方式中的、与差分图像合成的对象部位的图。
图6C是表示本实施方式中的、与除去被分割的容积的外侧的非心肌信号的清晰的图像合成的对象部位的图。
图7A是本实施方式中的、将分割前的LV切片图像在1个显示面板上一起显示为BBTI的时间周期的函数的示例的灌流视觉化的描绘图。
图7B是本实施方式中的、将分割后的LV切片图像在1个显示面板上一起显示为BBTI的时间周期的函数的示例的灌流视觉化的描绘图。
图7C是本实施方式中的、将分割后的LV切片图像在一个显示面板一起显示为BBTI的时间周期的函数的示例的灌流视觉化的描绘图。
图8A是本实施方式中的、左心室的冠状动脉区域的分割的描绘图。在本图中，使用米国心脏协会(AHA：American Heart Association)的6分割模型，得到左心室被分割的示例的冠状动脉区域的颜色区分显示。
图8B是本实施方式中的、左心室的冠状动脉区域的分割的描绘图。在本图中，使用AHA的6分割模型，得到左心室被分割的示例的冠状动脉区域的颜色区分显示。
图8C是本实施方式中的、左心室的冠状动脉区域的分割的描绘图。在本图中，使用AHA的6分割模型，得到左心室被分割的示例的冠状动脉区域的颜色区分显示。
图8D是本实施方式中的、左心室的冠状动脉区域的分割的描绘图。在本图中，使用AHA的6分割模型，得到左心室被分割的示例的冠状动脉区域的颜色区分显示。
图9A是本实施方式中的、使用不同的示例的视觉化表现，作为BBTI值的函数，对不同的切片各自的每段的灌流曲线的同时一面显示进行视觉化处理的描绘图。
图9B是本实施方式中的、使用不同的示例的视觉化表现，作为BBTI值的函数，对不同的切片各自的每段的灌流曲线的同时一面显示进行视觉化处理的描绘图。
图9C是本实施方式中的、使用不同的示例的视觉化表现，作为BBTI值的函数，对不同的切片各自的每段的灌流曲线的同时一面显示进行视觉化处理的描绘图。
图10A是本实施方式中的、用于对BBTI、与灌流相关联的信号强度(SI：SignalIntensity)、以及在LV的3D图像内构成实质上相邻的切片的LV切片编号的全部或者任一个之间的关系更好地进行视觉化以及理解的新的类型的靶心图。
图10B是本实施方式中的、用于对BBTI、与灌流相关联的信号强度、以及在LV的3D图像内构成实质上相邻的切片的LV切片编号的全部或者任一个之间的关系更好地进行视觉化以及理解的新的类型的靶心图。
图11是本实施方式中的、为了执行对磁共振成像收集到的BBTI标记添加以及标记删除的图像进行分析的系统的示例的实施方式，以流程图的形式表示图1的MRI系统(或者，其他的图像处理装置)所使用的示例的计算机程序代码构造的概略图。
具体实施方式
以下，参照附图，说明实施方式所涉及的磁共振成像装置(以下，适当地称为MRI(Magnetic Resonance Imaging)系统)以及图像处理装置。另外，实施方式并不限定于以下的实施方式。关于本申请人在2012年8月16日申请的美国专利申请第13/587，294号说明书，通过参照在本说明书中引用。
图1是表示构成为提供非造影动态MRI心肌灌流的分析以及视觉化的、示例的实施方式所涉及的MRI系统的高水平概略框图。图1所示的MRI系统包含架台10(概略剖面所示)、相互接口连接的各种关联系统构成要素20。至少架台10通常被配置于屏蔽室中。图1所示的MRI系统的形状实质上包含静磁场B0磁铁12、Gx、Gy、Gz倾斜磁场线圈组14、RF(RadioFrequency)线圈总成16同轴圆筒状的配置。沿着该圆筒状的要素阵列的水平轴线，存在实质上由被检体床，即，台11支承的被检体9(例如，用于心脏MRI的心脏)的、以包围对象的解剖学组织(即，对象部位(ROI：Region Of Interest)的方式示出的成像区域18。
MRI系统控制部22具备与显示部24、键盘/鼠标26、以及打印机28连接的输入/输出端口。对于本领域的技术人员而言不言而喻，显示部24也可以是还具备控制输入的多样性的触摸屏。
MRI系统控制部22与MRI序列控制部30接口连接，MRI序列控制部30控制Gx、Gy、Gz倾斜磁场线圈驱动器32、RF发送部34、以及发送/接收开关36(相同的RF线圈用于发送以及接收的双方的情况)。对本领域的技术人员而言不言而喻，能够使用多种不同种类的RF线圈(例如，全身线圈、表面线圈、鸟笼型线圈、线圈阵列等)，在与成像容积内的ROI之间进行RF信号的发送以及接收的双方或者一方。不言而喻，将一个或者多个合适的生理学变换器8粘贴于被检体的身体，将心电图(ECG：Electrocardiogram)信号、呼吸同步信号、以及脉搏波同步信号的全部或者任一个门信号供给至MRI序列控制部30。MRI序列控制部30还访问用于执行在MRI序列控制部30的防御范围内已经能够利用的MRI数据收集序列的合适的程序代码构造38，该程序代码构造38例如用于使用特定的MRI数据收集序列参数、划定一个或者多个ROI等的系统输入以及操作者的双方或者一方，产生非造影心脏MRI组织图像。
MRI系统的构成要素20包含向数据处理装置42供给输入的RF接收部40，以使得能够制成能够向显示部24发送的处理完成图像数据。MRI数据处理部42还构成为能够向图像重建程序代码构造44和MR图像存储部46(例如，用于保存来自与示例的实施方式以及图像重建程序代码构造44对应的处理的MR图像数据的)访问。
另外，图1是对MRI系统的程序/数据保存部50进行一般化处理的图。在MRI系统的程序/数据保存部50中，保存(例如，用于非造影心脏MRI动态心肌灌流分析以及视觉化的)程序代码构造、和保存于能够访问MRI系统的各种数据处理构成要素的计算机可读保存介质的、相关联的GUI(Graphical User Interface)、对该GUI的操作者输入等。对于本领域的技术人员而言不言而喻，可以将程序/数据保存部50分割成正常运转时迫切需要那样保存的程序代码构造的系统20的处理计算机中的各种计算机，且至少将一部分直接连结(即，代替普通地保存在MRI系统控制部22中或者直接连结)。
实际上，对于本领域的技术人员而言不言而喻，图1是将为了执行在本说明书中后述的示例的实施方式而增加了若干的变更的一般的MRI系统非常高度地简化的图。系统构成要素能够分割成各种逻辑收集的“块”，通常，包含多个数字信号处理装置(DSP：DigitalSignal Processor)、超小型运算处理装置、特殊用途处理电路(例如，高速模拟/数字(A/D：Analog to Digital)变换、高速傅里叶变换、阵列处理用等)。通常，当发生各时钟周期(或者规定数的时钟周期)时，这些处理装置的各个是物理数据处理电路从某一物理状态向另一物理状态前进的时钟动作型的“状态机器”。
在动作中，处理电路(例如，CPU(Central Processing Unit)、寄存器、缓冲器、计算单元等)的物理状态除了从某一时钟周期向另一时钟周期渐进地变化之外，所连结的数据保存介质(例如，磁性存储介质的位保存部)的物理状态也在那样的系统的动作中，从某一状态向另一状态变化。例如，在MR成像重建处理结束时，物理性保存介质的计算机可读可访问的数据值保存位置的阵列(像素值的多数位二进制显示)从几个事先的状态(例如，全部一律为“零”值、或者全部为“1”值)变为新的状态。在该新的状态下，那样的(例如，像素值的)阵列的物理位置的物理状态在最小值与最大值之间变动，能够表现现实世界的物理现象以及状况(例如，被图像化的部位空间一面的被检体的组织)。对于本领域的技术人员而言不言而喻，保存的数据值那样的阵列还表示且构成物理构造。即，当依次读入命令寄存器中并由MRI系统20的一个或者多个CPU执行时，以产生动作状态的特定序列，构成在MRI系统内移动的计算机控制程序代码的特定构造的方式，构成上述阵列。
下述的示例的实施方式提供用于进行MRI数据的收集以及处理的双方或者一方和MR图像的生成以及显示的双方或者一方的改良后的方法。
最初，预先说明以下的示例的实施方式所涉及的MRI系统的概要。示例的实施方式所涉及的MRI系统具备MRI序列控制部30和MRI系统控制部22。MRI序列控制部30关于长度不同的多个规定时间，执行从对标记区域施加标记脉冲开始经过规定时间后收集心脏的切片的数据的数据收集。例如，作为心脏的切片的数据，MRI序列控制部30收集从心尖部到心基底部之间的多个切片的数据。例如，MRI序列控制部30关于多个不同的TI(Time toInversion)，收集从心尖部到心基底部之间的多个切片相应的左心室短轴像的数据。
接着，MRI系统控制部22(还称为显示控制部)将关于多个规定时间收集到的多个数据的分析结果显示于显示部24，以使得能够识别规定时间与切片的关系。另外，在示例的实施方式中，MRI系统控制部22进行分析本身和分析结果的显示控制。
在该分析结果的显示中存在各种方法。例如，MRI系统控制部22将作为各数据的分析结果的分析图像排列成将规定时间的方向和切片的方向作为各轴的矩阵状，并由一面显示(参照图7A～7C)。在分析图像中，例如存在对图7A那样的原始图像(标记删除图像、标记添加图像、或者差分图像)按照像素单位分配了浓度或颜色的图像、对图7B那样的、为了去除非心肌信号而分割的图像按照像素单位分配了浓度或颜色的图像、或图7C那样的、进一步对心肌内进行分割，按照段单位分配了浓度或颜色的图像等。
另外，例如，作为分析结果，MRI系统控制部22显示通过按照每段对切片的数据进行分析而制成的、表示规定时间的方向的信号强度的变化的灌流曲线的曲线图(参照图9A～9C)。此时，当接受了操作者对显示于显示部24的图像(原始图像、或者分析图像)进行的选择操作时，MRI系统控制部22能够进行控制，以使得将与由操作者选择的图像对应的灌流曲线的曲线图显示于显示部24。
另外，例如，作为分析结果，MRI系统控制部22将切片的数据按照每段进行分析，根据各段的信号强度分配不同的浓度或者不同的颜色，显示将同心圆的半径方向作为规定时间的方向而制成的靶心图(图10A～10B)。
另外，在以下的示例的实施方式中，说明为MRI序列控制部30或MRI系统控制部22执行上述的处理，但实施方式并不限定于此，也可以由其他的控制部执行。
与强调灌流剂相关联，可能使用基于钆(Gd：Gadolinium)的造影剂的情况被熟知，但使用那样的造影剂大多数情况下不是检测心肌的梗塞性病变以及缺血性病变所能容许的MRI的技术。LGE(Late Gadolinium contrast Enhancement)是在注入钆造影剂注入后，通过从心肌观察到的MRI信号，观察钆造影剂的T1造影(清洗)，从而测定正常心肌与梗塞心肌之间的差异。为了得到强调下的心肌的现实的观察，为了在暂时强调状况下能够更好地观察血流的异常，能够通过被检体的身体的运动、或者(想要引起与由于身体的运动而引起的情况类似的心拍数的上升、心脏血管的膨胀等)被注入的药剂诱发性强调，诱发强调。
另外，ASL(Arterial Spin Labeling)法还用于基于非造影的心肌缺血评估，但该技术不会产生灌流曲线或者三维(3D)覆盖，对象限定于单一时刻的单一切片。对应的分析法与本发明者们的新的非造影动态MRI灌流分析以及视觉化明显不同。
本发明者们发现了避免那样的造影剂的使用的方法。该方法实施用于生成正常心肌、缺血性心肌、以及梗塞性心肌之间的识别所能够使用的心肌灌流曲线/表的数据的非造影(即，没有注入化学造影剂)MRI技术。实际上，能够只识别(在对血管插入支架的手术、外科的分流血管手术等血管再生技术中被处置的)被血管再生的梗塞性心肌。根据示例的实施方式，那样的心肌识别能够在所希望的任意的对象部位(例如，手术者定义的任意的对象部位(ROI)、标准的AHA的段、单一像素等)内进行。
避免了造影剂(例如，基于钆的)的注入，尽管如此，为了能够只在被检体强调下状态的期间，更好地检测存在或者可能更显著的异常，希望使用与(运动诱发性或者药物性的)被检体强调组合的示例的实施方式。
在示例的实施方式中，在通过最初的空间选择的RF脉冲(例如，一般空间选择性的180°反转脉冲)“添加标记”后，收集一组“标记添加”(2D或者3D)MRI k空间数据。在此，数据收集子序列在施加给反转延迟时间(TI)的时间之后开始。“标记删除”MRI数据类似的集合也使用相同的TI延迟间隔来收集，但不是最初的空间选择性的RF标记添加脉冲。该技术有时在某一技术领域中作为BBTI(Black Blood Time to Inversion)成像而被熟知。相对于多个TI时间中的各个的、一系列的那样的标记添加/标记删除的数据集在k空间中收集。
相对于各个所提供的TI时间，2D/3D傅里叶变换重建后的空间区域标记添加以及标记删除的图像数据集(即，周知的2DFT/3DFT重建处理的结果)(例如，按照像素单位)相减，作为时间的函数提供能够对血液灌流(MR信号强度)进行曲线图化或者表化处理的BBTI血液灌流图像。对于任意的所提供的对象部位(例如，AHA段、任意的手术者定义的ROI、或者单一的像素)，相对于时间轴绘制数据值得到的曲线图作为时间的函数提供血液灌流曲线。该时间的函数提供相对于“正常的”心肌、异常的缺血性心肌、梗塞性心肌、以及血管再生后的(处置后的)缺血性心肌的灌流曲线之间的几个次元微分。如预想的那样，梗塞性心肌没有示出灌流(即，在检测到的MRI信号强度中，没有峰值，即没有增加)。但是，在缺血性心肌中，当与正常的(或者被血管再生的)心肌相比较时，除了更小的检测信号强度之外，存在从发生开始被时间延迟的最大流量时间。从而，缺血性病变的位置能够基于时间以及振幅(或者曲线下的积分面积)的双方或者一方的比较和与规定的阈值等的比较的双方或者一方来识别。
(例如，如向操作者显示的那样，或者如之后为了向其他的操作者显示而保存数据那样)与包含显著的特征的对象部位(例如，规定的AHA心肌段、或者任意的操作者指示的ROI、或者单一的像素)相关联的灌流曲线以及对应的数据表的双方或者一方能够各种各样地描绘ROI，以使得不注入任何造影剂(例如，钆造影剂)而表示正常、缺血性、梗塞性、或者血管再生后的心肌。
图2A是本实施方式所涉及的MRI系统所使用的、标记添加、标记删除MRI数据收集序列的第1个示例的概略描绘图。在图2A所示的代表性的数据收集序列中，为了使“标记删除”以及“标记添加”的任一个子序列均同步，如图2A所示，使用ECG同步。在图2A的左下方所示的标记删除的子序列中，使用短的初始的固定延迟时间TD，一定在RST ECG信号中的所希望的心脏舒张期部分之间，开始数据收集子序列。在延迟时间TD之后，使用非区域选择180°RF脉冲(为了辅助抑制来自后续的不同的图像的背景信号。这是因为来自心肌的信号例如如MRI信号整体的约10％那样，比较小)。如果希望，实际的数据收集子序列能够与预备脉冲(例如，脂肪抑制脉冲)一起开始。之后，接着与切片选择性的(例如，α°、典型的情况为45°或者90°)NMR脉冲一起开始的所希望的数据收集子序列F，例如，开始切片选择性的180°RF再聚焦脉冲的序列那样的所希望的MRI数据收集子序列F，读出其中的RF自旋回波(SE(SpinEcho))响应并在倾斜磁场Gr脉冲继续期间中取出。另外，在各自旋回波之前，出现相位编码Ge倾斜磁场脉冲(为了取出分别对应的k空间的线，其相对于不同的回波进行变更)。对于本领域的技术人员而言不言而喻，该公知的MRI数据收集子序列也可以是bSSFP(balancedSteady-State Free Precession：平衡定常状态自由进动)类型(目前优选今)、或者高速自旋回波(FSE(Fast Spin Echo))类型、或者可能的其他的类型。
图2A的右下方所示的标记添加数据收集子序列与上述的标记删除子序列相同。其中，在延迟时间TD之后，还可以使用区域选择“标记添加”180°RF脉冲(例如，根据情况，由图2A所示的那样的虚线的同时行进倾斜脉冲表示的倾斜角的脉冲、以及Df频移脉冲)。对于本领域的技术人员而言不言而喻，实际上，这使血液的规定的流入容积返回到非反转磁化配向。即，对该血液的流入容积“添加标记”。从而，当该RF标记添加流动血液MR核进入下游的关心区域(ROI)时，将产生与上述标记删除子序列的响应不同的MR信号响应。
图2B是本实施方式所涉及的MRI系统所使用的、标记添加、标记删除MRI数据收集序列的第2个示例的概略描绘图。图2B所示的标记添加/标记删除交替式收集子序列与图2A所示的相同。其中，如观察得知的那样，不存在以背景抑制为目的的最初的非区域选择180°脉冲(标记添加子序列、标记删除子序列的任一个)。
如在图2A以及图2B的双方中由虚线所示的那样，为了相对于区域选择180°标记添加脉冲实现所希望的倾斜角，大概能够同时前进地使用各种所选择的大小的Gs、Gr、以及Ge倾斜磁场脉冲。
以上叙述了示例的实施方式中的MRI数据的收集，但重新对其进行说明。在示例的实施方式中，MRI系统通过对血液进行标记(添加标记)的方法，即，通过非造影的方法，收集描绘向心肌组织供给的血液的动态的图像组(心肌灌流图像组)。并且，MRI系统通过对收集到的图像组进行分析(分析)，显示分析结果来进行视觉化处理。
简单地说明对血液进行标记的方法。例如，当得到描绘出向心肌组织供给的血液的图像时，例如，作为对血液进行标记的标记脉冲，MRI系统例如对血液在心肌组织的上游流动的区域施加使纵向磁化反转180°的反转脉冲。接受了反转脉冲的施加的血液在纵向磁化的恢复过程中，具有能够与其他的血液或背景组织进行识别的信号强度。另外，反转脉冲的施加定时与MR信号的收集定时(例如，k空间的中心线的收集定时)之间的时间的长度典型的情况是与接受反转脉冲的施加的血液的移动距离成比例。因此，通过一边改变反转脉冲的施加定时与MR信号的收集定时之间的时间的长度一边收集MR信号，从而得到描绘流入心肌组织，并流出的血液的动态的图像组。另外，还存在大致与反转脉冲同时施加没有选择区域的非区域选择脉冲，使摄像区域整体的纵向磁化反转180°的方法。此时，心肌组织上游的血液接受2次反转脉冲的施加，因此，还具有能够与其他的血液或背景组织识别的信号强度。图2A表示施加非选择区域脉冲的方法，图2B表示没有施加非区域选择脉冲的方法。另外，标记脉冲并不限定于使纵向磁化反转180°的反转脉冲。
图3是表示本实施方式中的标记脉冲的施加区域的图。血液从心脏的左心室，经由图3所示的大动脉向全身送出。另外，在图3中省略图示，但从大动脉左右分支的动脉是右冠状动脉以及左冠状动脉，左冠状动脉进一步分为向前侧下行的血管和向后侧包围的血管。相对于心肌组织的血液的供给通过这些右冠状动脉(RCA(Right Coronary Artery))、左冠状动脉前下行枝(LAD(Left Anterior Descending coronary artery))、以及左冠状动脉回旋枝(LCX(Left Circumflex coronary artery))的3根血管进行。本实施方式捕捉流入流出心肌组织的血液的动态，并进行分析，因此，如图3所述，标记脉冲的施加区域优选设定于心脏上部的大动脉的附近。
另外，在本实施方式中，除了对血液进行标记的标记添加图像数据之外，收集没有对血液进行标记的标记删除图像数据，将这些差分图像作为分析对象。这是由于当将标记添加图像数据和标记删除图像数据相减时，减除施加了标记脉冲的血液以外的信号的结果，能够更准确地对施加了标记脉冲的血液进行分析。
另外，在本实施方式中，由于是在没有待机时间(没有操作者进行的操作介入)的一系列的脉冲序列中，收集全部标记添加图像数据以及标记删除图像数据的双方，因此，分别称为“标记添加数据收集子序列”、“标记删除数据收集子序列”等。图2A和图2B的不同是使摄像区域整体的纵向磁化通过非区域选择脉冲反转一次(图2A)，或者不反转(图2B)的不同。在前者的情况下，背景信号成为低信号。之后，鉴于通过标记添加图像和标记删除图像生成差分图像，可以说优选标记添加图像与标记删除图像之间的位置偏移的影响少的前者的方法。
“标记添加数据收集子序列”和“标记删除数据收集子序列”的不同是如将图2A的左右的图或者将图2B的左右的图进行比较知道的那样，是否施加区域选择的反转脉冲的不同。另外，如图2A所示，当还施加非区域选择脉冲时，非区域选择脉冲和区域选择脉冲典型的情况是大致同时施加。因此，在图2A或图2B中，“TI”表示从非区域选择脉冲的施加定时到数据的收集开始定时，大概也可以认为是从区域选择脉冲的施加定时到数据的收集开始定时的时间。
另外，在本实施方式中，优选MRI序列控制部30在一次的摒住呼吸(例如，20秒左右)的期间，收集完成多个切片相应的标记添加图像数据以及标记删除图像数据的集合。例如，MRI序列控制部30收集3个切片相应的从心尖部到心基底部的左心室(LV：LeftVentricle)切片图像。之后，MRI系统通过插补处理，得到6个切片相应的标记添加图像数据以及标记删除图像数据的集合。
另外，在本实施方式中，MRI序列控制部30关于长度不同的多个TI，执行上述的交替式的标记添加/标记删除子序列。例如，MRI序列控制部30针对TI＝200ms、TI＝400ms、TI＝600ms、TI＝800ms、TI＝1，000ms、TI＝1，200ms、TI＝1，400ms、以及TI＝1，600ms的各个TI，交替式地执行标记添加/标记删除子序列。不同的TI的子序列可以以没有待机时间(没有基于操作者的操作介入)的一系列的脉冲序列执行，也可以适当地插入待机时间来执行。TI的顺序也能够任意地变更。
另外，由哪一个TI进行收集，或者其间隔也可以适当地变更。这些TI可以通过在用于接受摄像条件的GUI上接受来自操作者的输入来设定，或者也可以在前序列之一中设定。通常，在基于MRI系统的检查中，包含用于收集各种诊断用的图像的成像序列组和在成像序列组之前进行的前序列组，这些一系列的序列组一边在期间插入操作者进行的操作等一边连续地依次执行。用于设定TI的前序列例如是一边变更TI一边在不同的每个TI得到各个的MR信号，将MR图像并列显示在GUI上使操作者进行选择，或者求MR图像上的线ROI(RegionOf Interest)的轮廓，从而设定之后成像序列所使用的合适的TI的前序列。这样的前序列有时被称为BBTI-prep等。另外，在前序列组中，此外，例如，包含用于收集定位图像的序列、用于调整磁场的不均匀性的序列、用于收集线圈的灵敏度图的序列等。
另外，在本实施方式中，MRI序列控制部30在不同的TI的子序列间，使收集的数据的心脏时相一致。例如，MRI序列控制部30关于任一个TI，在舒张期的心动相位收集数据。此时，非区域选择脉冲或区域选择脉冲的施加定时例如能够根据数据的收集定时对TI部分进行反算来求得。
另外，在本实施方式中，说明了分别收集标记添加图像数据组和标记删除图像数据组的方法，但实施方式并不限定于此，只要是通过非造影得到心肌组织的灌流图像的方法即可。即，收集标记删除图像数据组不是必须的构成。另外，在本实施方式中，例如，说明了对1切片的每一个，或者1切片编码的每一个，交替地收集标记添加图像数据组和标记删除图像数据组的方法，但实施方式并不限定于此。例如，在收集了全部的标记添加图像数据组之后，收集标记删除图像数据组，或者也可以相反。另外，在本实施方式中，说明了针对所有的TI收集标记删除图像数据组的方法，但实施方式并不限定于此，也可以在其他的TI中计算差分图像时使用针对某一TI收集到的标记删除图像数据组。
另外，在本实施方式中，说明了在没有待机时间(没有基于操作者的操作介入)的一系列的脉冲序列中连续地收集TI不同的图像数据组的方法，但实施方式并不限定于此。例如，也可以在每当改变TI时，插入待机时间。另外，在本实施方式中，说明了与ECG同步信号同步，以舒张期的心动相位收集图像数据的方法，但实施方式并不限定于此，也可以以其他的心动相位收集图像数据。此外，图2A、图2B、以及图3所示的本实施方式的收集的方法并不限定于此，还能够任意地变更。
以上，说明了MRI系统由非造影收集心肌灌流图像数据组的方法。以下，说明MRI系统、或者与MRI系统不同的图像处理装置进行将通过MRI序列控制部30收集到的心肌灌流图像数据组作为对象的分析的例子。本发明者们以下提出适合本发明者们的非造影动态心肌灌流技术的、可以说可能最优的示例的分析法。例如，新的方法对左心室内的血液灌流进行视觉化处理，能够对动态3D图像数据集进行处理，该数据集具备表示操作者设定的所有的段，即表示ROI的灌流曲线的能力。
对于本领域的技术人员而言不言而喻，本发明者们的灌流分析以及视觉化的方法能够编入图1的MRI系统，或者也可以作为远离收集原来的标记添加以及标记删除的图像数据的图1的MRI系统的不同的图像分析/显示系统(图像处理装置)的一部来执行。
为了得到由本发明者们的非造影灌流技术得到的动态3D图像，本发明者们提出以下的目前所希望的示例的分析步骤(通常不需要或者希望以下的分析步骤的全部)。
1.执行标记添加图像以及标记删除图像间的复数数据的减法。
2.执行图像位置对准，即，执行相对于不同的BBTI的3D图像的、刚性或者非刚性的位置对准。
3.为了进行与由于对被图像化的切片产生影响的磁化率以及带有标记的切片的双方或者一方而造成的心肌信号的缺失相关的确认，制成标记添加图像、标记删除图像的直方图。通过直方图的相减，能够检测相对于被图像化的心肌的不合适的磁化率以及带有标记的切片的影响。
4.使用心肌的分割。
5.将被分割的心肌以彩色图形式显示。
6.涵盖全部切片对每个段，即，对每个ROI制成灌流曲线。
7.同时发生地显示相对于BBTI的3D切片的全部。
8.当选择段以及ROI的双方或者一方时，同时发生地自动显示分别对应的灌流曲线。
9.为了定量化目的进行灌流曲线的拟合。
(1：针对差分图像的生成)
MRI系统控制部22将由同一切片且同一BBTI收集到的、标记删除图像以及标记添加图像的集合作为对象，按照像素单位进行基于复数的差分处理。另外，如上所述，在本实施方式中，标记删除图像和标记添加图像在摒住呼吸下交替地收集，且心动相相位匹配地收集，因此，认为不一定在差分时必须进行位置对准，但实施方式并不限定于此。也可以在差分处理之前，进行标记删除图像与标记添加图像的位置对准。图4是表示本实施方式中的、标记添加以及标记删除的图像的像素单位的减法的图。如图4所示，MRI系统控制部22按照像素单位进行从标记添加图像的实数部中减去标记删除图像的实数部的计算。另外，MRI系统控制部22按照像素单位进行从标记添加图像的虚数部中减去标记删除图像的虚数部的计算。相减得到的差分图像必须是复数形式(实数部以及虚数部、R+jI)。这是由于相减的像素的强度对由于在标记添加以及标记删除的信号间可能发生的相移而造成的信号变化迟钝。从而，原来的标记添加以及标记删除的图像还必须是复数的形式。复数的数据在相减后，用于制成用于灌流分析的强度图像。
在此，说明在本实施方式中进行基于复数的差分处理的理由。原子核的磁化矢量由复合平面上的实数分量(同相位分量)以及虚数分量(正交相位分量)表现。因此，在MRI中，收集实数分量以及虚数分量各自的k空间数据，通过傅里叶变换生成实数图像以及虚数图像，之后，生成作为绝对值图像的振幅图像或相位图像。
另外，例如，在标记删除图像与标记添加图像之间，存在是否施加区域选择脉冲的差异。该差异表现为纵向磁化分量的差，但在施加激发脉冲后的xy平面上，表现为横磁化分量的相位的差。例如，这是由于在朝上的磁化矢量落入xy平面时和朝下的磁化矢量落入xy平面时，其相位不同。因此，在本实施方式中，为了准确地考虑该相位差，MRI系统控制部22进行基于复数的差分处理。
(2：针对图像位置对准)
在本实施方式中，针对同一BBTI，收集多个左心室切片图像(另外，一部分的左心室切片像也可以通过插补处理生成)，另外，各个位置的左心室切片图像针对不同的BBTI收集。即，不同的BBTI的3D图像在不同的收集时间收集，因此，有时需要不同的BBTI图像间的位置对准。因此，MRI系统控制部22在同一BBTI的左心室切片图像间进行位置对准，同时在不同的BBTI的左心室切片图像间进行位置对准。另外，实施方式并不限定于此，也可以省略这些位置对准的双方或者一方。
图5A是表示本实施方式中的不同的BBTI间的位置对准前的左心室心肌的强度图像的图，图5B是表示本实施方式中的不同的BBTI间的位置对准后的左心室心肌的强度图像的图。图5A以及5B的图像表示涵盖所有的BBTI的相对于成像平板的3D位置对准(为了简化说明，只示出1个切片)。
MRI系统控制部22进行刚性位置对准或者非刚性位置对准。即，位置对准并不限定于刚性位置对准。例如，某一BBTI中的左心室切片图像的环状在另一BBTI中有时是稍微不同的形状。此时，希望执行非刚性位置对准。
在提出的方法中，优选左心室的位置对准。在图5A以及图5B中，示出位置对准窗口只配置于左心室的例子。当将位置对准窗口只配置于左心室时，能够进行相对于左心室的局部的位置对准。例如，MRI系统控制部22将不同的BBTI的左心室切片图像组中，至少1个左心室切片图像显示于显示部24，同时将位置对准窗口(例如，箱型)显示于显示部24，接受来自操作者的操作。例如，操作者将位置对准窗口移动到包含左心室切片图像上的左心室心肌的那样的位置。接着，MRI系统控制部22将其他的BBTI的左心室切片图像显示于显示部24，同时，在之前由操作者设定的位置显示位置对准窗口。假设，当在此显示的位置对准窗口不包含左心室心肌时，例如，操作者如今移动左心室切片图像，以使得在位置对准窗口内包含左心室心肌。例如，通过重复该处理，从而MRI系统控制部22进行不同的BBTI的左心室切片图像间的位置对准。在位置对准前的图5A中，知道位于从位置对准窗口稍微向下方偏移的位置的心肌在位置对准后的图5B中，大致在位置对准窗口的中央。另外，该位置对准的处理只不过是一个例子，位置对准例如也可以通过基于图像分析的自动的处理来进行。另外，也可以将这些方法适用于同一BBTI的左心室切片图像间。
在此，当施加了非区域选择脉冲时，心腔与周围的心肌之间的造影在几个BBTI中反转。此时，标记添加以及标记删除的两图像必须用于位置对准处理。例如，能够选择正对比(心肌信号强度SI＞LV腔血液SI)的图像，关于它们执行位置对准。另外，能够关于负对比的图像，执行其他的位置对准处理。在心腔中血液的替换剧烈。因此，如参照图7A得知的那样，心腔内部的信号值根据BBTI有时高，有时低。即，根据BBTI，存在心肌的信号强度SI比心腔血液的信号强度SI高时(正对比)的情况，还存在心腔血液的信号强度SI比心肌的信号强度SI高时(负对比)的情况。当假设选择进行位置对准时成为基准的图像时，希望根据BBTI选择最优的对比度(正对比或者负对比)的图像。另外，假设为了进行复合位置对准处理而记录相对于参照图像的各图像的像素的位置对准移动。有时实施最高的位置对准有时需要基于目视的手动移动。
(3：心肌信号缺失的确认)
在此，针对心肌信号缺失的确认进行说明。另外，不是一定在该阶段进行心肌信号缺失的确认处理，也可以与其他的处理平行，或者在其他的处理的前后等、任意的定时进行即可(例如，也可以是1：差分图像的生成、2：图像位置对准之前的定时)。即，当检测到存在心肌信号的缺失时，MRI系统控制部22通过将该事实与之后向操作者提供的图像或分析结果(分析图像、曲线图、靶心图等)等一起作为信息输出，从而唤起阅览图像或分析结果的工作人员(例如，医师)的注意(例如，不要将心肌信号的缺失与病变混同)。在心肌信号的缺失中，存在摄像区域(成像平板)与施加标记脉冲(在本实施方式中，区域选择脉冲)的标记区域重合时发生伪影的类型的缺失和由于磁化率的差异而发生伪影的类型的缺失。
例如，MRI系统控制部22在标记删除图像、标记添加图像、或者差分图像上，设定横穿左心室心肌的线ROI，得到线ROI上的信号轮廓(在上述中，将该信号轮廓称为直方图)。并且，MRI系统控制部22根据该信号轮廓，判定在该线ROI上，是否存在心肌信号的缺失(例如，局部的信号值的降低)。MRI系统控制部22对图像整体(适当地离散地)设定线ROI，判定有无心肌信号的缺失。另外，例如，MRI系统控制部22也可以使用减除标记添加图像上的信号轮廓和标记删除图像上的信号轮廓得到的信号轮廓，来判定有无心肌信号的缺失。
并且，当判定为存在心肌信号的缺失时，MRI系统控制部22对判定对象的图像附带该信息，之后当显示图像或分析结果时，作为参考信息，显示该缺失信息等。另外，有无缺失的判定并不限定于自动的方法，例如，还能够通过显示信号轮廓，使操作者通过目视来确认那样的手动的方法进行。另外，在此，说明了在标记删除图像、标记添加图像、或者差分图像上设定线ROI的例子，但实施方式并不限定于此，例如，也可以在后述那样的分割后的图像上设定线ROI。或者，在一系列的脉冲序列组中，作为其他的串联，例如，当事先收集了电影图像时，例如，MRI系统控制部22通过将该电影图像在显示部24上进行动画显示，从而使操作者能够通过目视确认缺失。另外，在此，叙述了当判定为存在心肌信号的缺失时，作为参考信息，显示缺失信息那样的方法，但并不限定于此。例如，MRI系统控制部22在使用该缺失信息进行上述的差分图像的生成或图像位置对准时，对对象的图像实施校正。
(4：心肌的分割)
图6A是表示本实施方式中的、对标记添加或者标记删除的图像施加的分割的图，图6B是表示本实施方式中的、与差分图像合成的对象部位的图，图6C是表示本实施方式中的、与除去被分割的容积的外侧的非心肌信号的清晰的图像合成的对象部位的图。
如图6A、6B、以及6C所示，LV分割通过沿着连续的BBTI排列LV，从而在位置对准后实施。在标记删除图像中，能够将心内膜以及心外膜的轮廓按照每个切片(手动或者半手动)描绘并保存。接着，能够将保存的轮廓粘贴在差分图像上。为了只对LV心肌进行视觉化处理，除去其他的信号即可(例如，参照图6C的“清晰的”图像)。另外，为了排除来自伪影(例如，磁化率伪影)以及带有标记的切片的干涉的所有的影响，必须慎重地对LV轮廓进行定位。
即，如图6A所示，MRI系统控制部22在标记添加图像或者标记删除图像上，操作者手动地接受心内膜以及心外膜的轮廓的指定，从而进行分割。或者，MRI系统控制部22通过对标记添加图像或者标记删除图像进行图像处理，从而从标记添加图像或者标记删除图像中提取心内膜以及心外膜的轮廓。
接着，如图6B所示，MRI系统控制部22在差分图像上，嵌入由标记添加图像或者标记删除图像得到的心内膜以及心外膜的轮廓。并且，如图6C所示，MRI系统控制部22使用心内膜以及心外膜的轮廓，从差分图像中除去非心肌信号。即，MRI系统控制部22从差分图像中，除去心内膜的内侧的区域和心外膜的外侧的区域。另外，在该阶段中，例如，当操作者发现心肌信号的缺失时，例如，MRI系统控制部22通过接受操作者进行的调整轮廓(例如，为了除去缺失部分而对轮廓进行调整)的操作，从而能够将心肌信号的缺失从分析对象中除去。另外，在手动进行心内膜以及心外膜的轮廓的指定(或者自动提取)中，MRI系统控制部22能够从标记添加图像或者标记删除图像中，适当地选择易于指定轮廓(或者，易于提取)的图像(例如，正对比或负对比高的图像)。
(5：心肌的彩色图显示)
图7A是本实施方式中的、将分割前的LV切片图像在1个显示面板上一起作为BBTI的时间周期的函数显示的示例的灌流视觉化的描绘图。另外，图7B以及图7C是本实施方式中的、将分割后的LV切片图像在1个显示面板上一起作为BBTI的时间周期的函数显示的示例的灌流视觉化的描绘图。灌流视觉化能够如图7A～7C那样实施。在此，从所有的BBTI中的心尖部到心基底部的切片的LV图像由一面显示。在该一面显示中，能够观察左心室的血液灌流(为了更好地对信号强度变化进行视觉处理，希望利用彩色图)。图7A表示LV分割前的视觉化。图7B以及图7C表示分割后的视觉化。
即，如图7A或图7B所示，MRI系统控制部22分别对分割前或者分割后的LV切片图像，按照像素单位分配与该信号值对应的颜色，生成彩色图图像。并且，MRI系统控制部22例如将横轴方向作为BBTI方向，将纵轴方向作为从心尖部到心基底部的切片方向，矩阵状地将所有的彩色图图像由一面显示。另外，一面显示的方法并不限定于此。例如，横轴和纵轴也可以相反。另外，例如，还能够根据来自操作者的指定，适当地省略显示一部分的行或一部分的列。
另外，在上述的例子中，说明了按照像素单位分配与信号值对应的颜色的方法，但实施方式并不限定于此。如图7C所示，MRI系统控制部22能够进一步对心肌进行分割，按照被分割的段的单位，分配与该段内的信号值(例如，平均值)对应的颜色，生成彩色图图像。另外，当显示图7A～7C的彩色图图像时，MRI系统控制部22能够将表示颜色的分割与信号强度的关系的彩色条一起显示。
图8A～8D是本实施方式中的、左心室的冠状动脉区域的分割的描绘图。在本图中，使用AHA的6分割模型，得到左心室被分割的示例的冠状动脉区域的颜色区分显示。例如，LV的冠状动脉区域的分割也可以是标准的AHA的6段模型，或者也可以根据需要是任意的数量的其他的操作者定义的段。图8A、8B、8C、以及8D所示的AHA6段从LV与RV(右心室)之间的槽开始自动地顺时针地移动。即，槽成为作为用于识别冠状动脉区域的开始点的标记。添加了各连续编号的AHA段如图8A～8D所示，为了在所有的段间识别强度变化并示出，能够进行标记添加、平均值化、以及彩色图化处理。
如上所述，相对于心肌组织的血液的供给通过右冠状动脉(RCA)、左冠状动脉前下行枝(LAD)、以及左冠状动脉回旋枝(LCX)的3根血管来进行。另外，在心肌组织中，段1以及段2通过左冠状动脉前下行枝接受血液的供给，段3以及段4通过右冠状动脉接受血液的供给，段5以及6通过左冠状动脉回旋枝接受血液的供给。
例如，如图8A所示，首先，MRI系统控制部22在标记添加图像或者标记删除图像上，操作者手动或者自动地将心肌整体分割成AHA的6段。此时，当自动地进行分割时，例如，将存在于LV与RV之间的界标作为起点，进行分割。接着，如图8B所示，MRI系统控制部22在差分图像上，嵌入由标记添加图像或者标记删除图像得到的6段的边界，从差分图像中除去非心肌信号。
接着，MRI系统控制部22对该每个段对信号值进行分析。即，例如，MRI系统控制部22首先对每个像素进行信号值的分析，接着，对每个段计算每个像素的分析结果的平均值，从而计算每个段的分析结果。另外，平均值的计算只不过是一个例子。并且，如图8C或8D所示，MRI系统控制部22按照段的单位，分配与该段内的信号值(例如，平均值)对应的颜色，生成彩色图图像。图8C是由彩色图图像，或者灰度(灰色的浓淡)显示的例，图8D是分配不同的颜色的彩色图图像的例子。MRI系统控制部22也可以将这样生成的彩色图图像如上述的图7C那样，由一面显示。另外，如图8C或8D所示，MRI系统控制部22相对于各彩色图图像，能够适当地显示切片编号、BBTI值、或各段的识别编号等。
(6：每个ROI的灌流曲线的制成)
另外，MRI系统控制部22能够制成按照每段表示信号值的时间变化的曲线图，并显示。图9A～9C是本实施方式中的、使用不同的示例的视觉化表现，作为BBTI值的函数，对不同的切片各自的每段的灌流曲线的同时一面显示进行视觉化的描绘图。沿着不同的BBTI图像的各个不同的切片的各段的灌流曲线能够如图9A～9C所示的那样制成。当使用所有的段的被保存的ROI时，能够对所有的切片内的所有的段自动地生成灌流曲线。图9A的例子表示来自一被检体的非平滑化的原始的灌流曲线。灌流曲线还能够由操作者指定的任意的ROI生成。
如图9C所示的那样，为了辅助进一步的定量化分析，也可以适用以往的多项式曲线拟合法，即，曲线平滑化法。两个曲线拟合方程式希望用于使灌流曲线最好地拟合。拟合后的参数必须能够表示灌流峰值强度、峰值的定时、以及峰值下的面积。希望所有的可识别的异常的参数(例如，如通过被曲线拟合的灌流数据的定量分析得知的那样)在被拟合的参数表内、以及直接显示的灌流曲线上的双方或者一方中必须添加标记。同样地，与检测到的异常对应的所有的部位(单个或者多个)希望在进行对应的AHA段的视觉化时必须直接添加标记。
图9B是模拟灌流曲线的3D视觉化的图。在该图中，(使用x，y，z的直角坐标轴的一般的方向)，设x轴为BBTI值，设y轴为相对信号强度(灌流)，设z轴为切片编号进行曲线图化。
灌流曲线在图9A～9C中由一面显示，但当操作者选择了特定的分别对应的切片或者段时(例如，通过特定的BBTI中的该选择的切片或者特定的切片的段上“点击”确定位置的鼠标箭头)，能够取出相对于图7A或者图7B内的重合的显示的任意的一灌流曲线。
即，MRI系统控制部22按照每段对各LV切片图像的信号值进行分析，按照每段制成表示BBTI方向的信号值的变化的曲线。并且，例如，如图9A或9C所示，MRI系统控制部22将6个切片相应的灌流曲线进行排列，由一面显示。另外，此时，MRI系统控制部22能够将各段的灌流曲线的差异以颜色、或线种类的差异来表现。另外，例如，当显示图7A～7C时，如果在图7A～7C的显示上接受规定的LV切片图像的选择，则MRI系统控制部22也可以适当地显示与所接受的LV切片图像对应的曲线图。或者，相反，当显示图9A～9C时，如果在这些曲线图上接受规定的曲线图的选择，则MRI系统控制部22也可以显示与所接受的曲线图对应的LV切片图像(例如，1行相应的)。这样，LV切片图像或曲线图的显示并不限定于上述的例子，还能够显示其一部分，或者使用相互的对应关系成为能够调出的构造等，能够任意地变更。
(其他：靶心图)
图10A～10B是本实施方式中的、用于对BBTI、与灌流相关联的信号强度、以及在LV的3D图像内构成实质上相邻的切片的LV切片编号的全部或者任一个之间的关系更好地视觉化以及理解的新的类型的靶心图。图10A～10B是使用周知的靶心技术的变形表示新的灌流的视觉化的图。在此，在图10A中，对于所提供的切片编号，同心圆将AHA段与沿着半径方向示出的BBTI值一起示出。各段希望相对于平均信号强度(SI：Signal Intensity)或者相对灌流值用颜色区分。在该视觉化中，作为相对于各心脏段的BBTI的函数能够容易地观察SI的倾向并易于理解。在图10B中，对于所提供的BBTI值，同心圆还示出AHA段，目前，切片编号沿着半径方向示出。在此，各段还希望为了表示平均SI或者灌流值而用颜色区分。目前，在该视觉化中，作为所提供的BBTI中的切片编号的函数能够容易地观察SI的倾向并易于理解。
即，作为分析结果的一个例子，例如，MRI系统控制部22制成将同心圆的半径方向作为BBTI方向的靶心图，显示制成的靶心图。另外，例如，MRI系统控制部22制成将同心圆的半径方向作为切片方向的靶心图，显示制成的靶心图。另外，在图10A或图10B中，以AHA的6段为例进行了说明，但分割例如也可以是17段等、其他的分割。这在上述的其他的分析中也相同。
示例的分析方法特别地为了使用本发明者们的非造影动态心肌灌流技术而设计。分析步骤均有益于对心肌内部的血液的灌流进行视觉化处理，根据健康的心肌识别梗塞性部位或者缺血性部位。生成的灌流曲线对于缺血性疾病或者梗塞的定量评估而言很重要。另外，作为分析结果的显示方法，说明了彩色图图像、灌流曲线的曲线图或靶心图等，但MRI系统控制部22例如也可以根据操作者的希望适当地显示所需的图像。
(示例的计算机程序代码构造)
图11是本实施方式中的、为了执行对磁共振成像收集到的BBTI标记添加以及标记删除的图像进行分析的系统的示例的实施方式，以流程图的形式表示图1的MRI系统(或者，其他的图像处理装置)所使用的示例的计算机程序代码构造的概略图。图11的步骤S900表示对动态非造影标记添加以及标记删除的BBTI图像数据的分析的入口。对于本领域的技术人员而言不言而喻，该流程图表示在能够通过高水平的应用程序或者基本软件调出的计算机程序子程序中被发现的那样的、可执行的计算机程序代码构造。
MRI系统控制部22根据需要，显示最新(最初或者最后使用的)设定的初始参数(步骤S902)，生成每个BBTI的强度图像(步骤S904)。MRI系统控制部22确认操作者是否接受这些初始参数(步骤S908)，当没有接受时(步骤S908，“否”)，显示S906所示的参数设定画面，能够接受进一步的基于操作者的调整/输入。接着，这些参数再次由步骤S902显示，在步骤S908中确定受理或者不受理。另外，作为参数，准备以下等选项：A.由刚性或者非刚性的哪一个来进行位置对准，B.是否进行伪影的确认，C.由AHA或者用户定义的哪一个来进行分割，D.是否进行全部切片对BBTI的一面显示，E.是否进行灌流曲线的显示。另外，在图11的例子中，经由步骤S904的处理进行参数的设定，但不一定需要按照该顺序进行处理。例如，也可以不执行步骤S904，而只进行步骤S902、908、906的处理。
一旦知道能够受理子程序的初始设定(实际上，对操作者进一步提供那样的调整的机会)，接着，MRI系统控制部22使用复数计算将预先收集到的标记添加以及标记删除的图像相减，生成3D图像的每个BBTI值以及每个切片的强度图像(步骤S904)。
之后，MRI系统控制部22为了得到不同的BBTI值中的3D图像执行图像位置对准(步骤S910)。为了实施图像位置对准，根据需要也可以允许基于操作者的手动辅助。当然，也有时不需要图像位置对准。
当希望进行伪影检测时(步骤S912，“是”)，接着，MRI系统控制部22制成标记添加以及标记删除的图像的直方图，并显示制成的直方图。进行减法，以使得能够提供表示磁化率的误差与由于在图像数据收集中对心肌产生影响的添加有标记的脉冲而造成的误差的双方或者一方的数据。对本领域的技术人员而言不严而喻，MRI系统控制部22进行检测是否存在超过某一一定的阈值的误差的判定，当存在误差时，进行校正动作与操作者辅助的申请的双方或者一方，或者可能只能结束处理。
MRI系统控制部22也可以根据需要执行心肌的分割(步骤S916)。接着，MRI系统控制部22将被分割的心肌的图像根据信号强度以被彩色图化的像素显示(步骤S918)。另外，MRI系统控制部22对3D图像的每个切片的每个段(或者ROI)制成灌流曲线。MRI系统控制部22也可以将LV切片图像以及灌流曲线的双方或者一方相对于BBTI值进行显示(步骤S922)。如上所述，用于3D图像的彩色化切片图像希望作为BBTI值的函数由一面显示。同样地，3D图像的每段以及每个切片的灌流曲线希望显示于单一面板。
当LV切片图像由单一面板显示时，接着，提供给操作者通过显示部(例如，通过鼠标或触摸等)选择特定的段，即选择ROI的机会(步骤S924)。当进行那样的操作者选择时，接着，MRI系统控制部22显示与该特定的段对应的灌流曲线(步骤S926)。并且，提供给操作者是否结束处理的选项(步骤S928)。当希望处理结束时，接着，MRI系统控制部22结束该子程序，控制返回到呼叫高水平程序或者基本软件。
显示所有的类型，即，除了灌流关联的图像之外，来自沿着BBTI轴的磁共振血管摄影法(MRA：Magnetic Resonance Angiography)等的非造影磁共振的多切片图像被认为新鲜且有利。例如，为了对通过磁共振成像取得的图像进行分析而计算机化的系统包含与相关联的存储部、显示部、以及输入/输出端口连结的至少1台计算机处理装置，可以构成为，(a)在对象部位(ROI)中，取得多个BBTI间隔中的每个间隔的左心室心肌的多切片非造影MR图像，(b)作为用于3D图像的多个切片的、以及用于单一显示面板内的多个BBTI值的BBTI的函数显示从心尖部到心基底部的LV切片图像。作为单一的显示面板内的BBTI的函数，MR切片那样的视觉化(例如，与图7A、7B、9A～9C、以及10A～10B的全部或者任一个所示的那样的视觉化相同)有益于操作者迅速地“看到”一系列的BBTI值与根据变动的BBTI值取得的多种类型的MR图像之间的重要的关系。
(其他的实施方式)
实施方式并不限定于上述的实施方式。
(具体的数值、处理步骤)
在上述的实施方式中示例的具体的数值或处理步骤只不过是一个例子。例如，在上述的实施方式中，使用图11说明了处理步骤的一个例子，但实施方式并不限定于图11所示的处理步骤，也能够根据运用的方式等适当地变更。例如，图11所示的处理的顺序能够任意地变更。
(触发信号)
在上述的实施方式中，说明了将ECG信号作为触发信号一边进行心电同步一边进行数据收集的例子，但实施方式并不限定于此。也可以代替ECG信号，将脉搏波信号、呼吸信号等其他的生物体信号或MRI系统的时钟信号等作为触发信号来使用。
(图像处理装置)
在上述的实施方式中，说明了MRI系统执行数据收集、分析、分析结果的显示的全部的例子，但实施方式并不限定于此。例如，包含MRI系统和图像处理装置的图像处理系统也可以执行上述的各种处理。在此，所谓图像处理装置例如是工作站、PACS(PictureArchiving and Communication System)的图像保管装置(图像服务器)、阅览器、电子病历系统的各种装置等。此时，例如，MRI系统进行基于MRI序列控制部30的数据收集。另一方面，图像处理装置具备相当于上述的MRI系统控制部22的控制部，通过从MRI装置100或者从图像服务器经由网络接收通过MRI系统取得的MR数据或k空间数据，或者通过经由记录介质接受操作者进行的输入，并存储于存储部。并且，图像处理装置将存储于存储部的该MR数据或k空间数据作为对象，执行上述的各种处理(例如，基于MRI系统控制部22的分析处理或显示控制处理)即可。
如上所述，根据至少一个实施方式的磁共振成像装置以及图像处理装置，对通过磁共振成像(MRI：Magnetic Resonance Imaging)收集的BBTI(Black Blood Time toInvers ion)标记添加以及标记删除的图像进行分析，作为表示对象部位(ROI：Region OfInterest)内的血液灌流的时间(BBTI值)的函数，生成不同的强度的3D图像。例如，具有正常以及异常的心肌组织的不同的值的ROI的灌流数据针对3D图像的多个切片，以及针对多个BBTI值，在单一的显示面板内显示。
(1)一种对通过磁共振成像收集到的图像进行分析的计算机化的系统，其特征在于，上述系统与相关联的存储部、显示部、以及输入/出力端口连结，且包含至少1台计算机处理装置，该计算机处理装置构成为：(a)在对象部位(ROI)中，收集每个BBTI的左心室(LV：Left Ventricle)心肌的多切片非造影MR图像，(b)针对3D图像的多个切片，以及针对多个BBTI值，将从心尖部到心基底部的LV切片图像显示在单一的显示面板内。
(2)一种对通过磁共振成像收集到的BBTI标记添加以及标记删除的图像进行分析的计算机化的系统，其特征在于，上述系统与相关联的存储部、显示部、以及输入/输出端口连结，且包含至少1台计算机处理装置，该计算机处理装置构成为：(a)在对象部位(ROI)中，作为表示血液灌流的时间的函数，为了生成不同的强度的3D图像，对多个BBTI中的每个BBTI，使用复数值计算，按照像素单位，预先将收集到的标记删除以及标记添加的图像进行相减，(b)针对3D图像的多个切片，以及针对多个BBTI值，将具有作为相对于正常以及异常的心肌组织的差的值的ROI的灌流数据显示在单一显示面板内。
(3)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为在上述3D图像之间，实施3D图像位置对准。
(4)根据(3)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为执行刚性或者非刚性的3D图像位置对准。
(5)根据(4)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为为了实施上述3D图像位置对准，接受基于操作者的手动的排列配置输入。
(6)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置(a)实施上述ROI被图像化的左心室心肌组织的分割，(b)构成为在一个显示面板上作为BBTI的函数显示LV切片图像，在此，上述显示的LV切片图像包含：(i)分割前的LV切片图像、(ii)分割后的LV切片图像、(iii)AHA(American Heart Association)的分割后的LV切片图像中的至少一个。
(7)根据(2)所述的系统，其特征在于，上述ROI包含左心室心肌，上述至少1台计算机处理装置还构成为由至少一个靶心构造显示上述灌流数据，上述靶心构造表示(a)具有沿着半径方向曲线图化的不同的BBTI值的每个值的沿着同心圆而曲线图化的灌流值的LV心肌的所提供的分割切片、和(b)相对于所提供的BBTI值的、沿着半径方向曲线图化的LV心肌的不同的切片中的每个切片的沿着同心圆曲线图化的LV心肌灌流值中的至少一个。
(8)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为，(a)生成直方图，(b)将标记添加以及标记删除的图像的上述直方图相减，上述减法结果提供表示由于磁化率而造成的伪影、和由在图像数据收集处理中对心肌图像产生影响的添加有标记的脉冲而造成的伪影的双方或者一方的数据。
(9)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为，(a)为了上述ROI的左心室心肌图像数据而执行冠状动脉区域的分割，(b)生成上述左心室心肌被分割的每个部分的平均图像数据值，(c)被分割的部分显示伴随着表示分别相关联的不同的平均图像数据值的视觉差而显示的上述被分割的左心室心肌。
(10)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为(a)作为3D图像的每个切片的BBTI的函数生成至少1个血液灌流曲线，(b)由一个显示面板同时显示所得到的多个灌流曲线，上述灌流曲线表示：(a)作为曲线拟合前的BBTI的函数的原来的灌流数据的阵列、(b)与作为BBTI的函数的原来的灌流数据拟合的灌流曲线的阵列、(c)由与作为BBTI的函数的原来的灌流数据拟合的灌流曲线得到的定量化的灌流曲线数据中的至少一个。
(11)根据(2)所述的系统，其特征在于，上述至少1台计算机处理装置还构成为，当对用于3D图像的各切片的上述ROI的多个被分割的部分中的每个部分作为BBTI的函数生成血液灌流曲线，操作者选择切片的段时，同时显示分别对应的灌流曲线。
(12)一种磁共振成像系统，其特征在于，构成为包含：MRI架台，具备形成能够配置被检体的对象部位(ROI)的成像容积的静磁场磁铁、倾斜磁场线圈、以及至少一个RF(RadioFrequency)线圈；MRI控制电路，构成为与上述MRI架台内的控制构成要素连接，当将ROI配置在上述架台内时，执行施加从被检体组织取出MRI信号的RF脉冲以及倾斜磁场脉冲的MRI数据收集序列，使上述取出的MRI信号(a)在数据收集子序列中与最初的空间选择性的RF标记脉冲的使用连动(标记添加)，且(b)在数据收集子序列中不与最初的空间选择性的RF标记脉冲的使用连动(标记删除)，收集MR图像数据并进行处理，上述MRI控制电路(a)不注入造影剂，而使用多个反转时间(TI：Time to Inversion)间隔中的每个间隔的上述标记添加以及标记删除的数据收集子序列取得用于被检体ROI的多维MR k空间数据，(b)将上述收集到的k空间数据重建为空间区域标记添加以及标记删除的MR图像数据，(c)使用复数值计算将多个TI间隔中的每个间隔的上述收集到的标记添加以及标记删除的MR图像数据相减，作为表示正常、缺血性、以及梗塞性的组织间出现差异的、相对于上述ROI的、对应的流动时间、以及相对的峰值血流强度的双方或者一方、或者存在这些不足的上述ROI的血液灌流的时间的函数生成不同的强度图像数据，(d)针对3D图像的多个切片以及多个BBTI值，将具有正常、缺血性、以及梗塞性的组织的不同的值的上述ROI的灌流数据在单一显示面板上显示。
根据以上所述的至少一个实施方式的磁共振成像装置，能够恰当地显示由非造影得到的心肌灌流图像的分析结果。
虽然说明了本发明的几个实施方式,但这些实施方式是作为例子而提示的,并不意图限定本发明的范围。这些实施方式能够以其他的各种方式进行实施，在不脱离发明的要旨的范围内，能够进行各种的省略、置换、变更。这些实施方式或其变形与包含于发明的范围或要旨中一样，包含于权利要求书记载的发明及其均等的范围中。

Claims (6)

1.一种磁共振成像装置，其中，具备：
序列控制部，关于长度不同的多个规定时间，一边执行数据收集，一边使收集的数据的心脏时相一致，上述数据收集是在从对标记区域施加标记脉冲开始经过上述规定时间后通过非造影的方法一边进行心电同步一边收集心脏的切片的数据的数据收集；和
显示控制部，将是关于多个上述规定时间收集到的多个上述数据的分析结果亦即分析图像且按照段的单位与段内的信号值对应地分配了颜色的彩色图图像以排列成将接受所述标记脉冲的施加的血液的移动距离和上述切片作为各轴的矩阵状的方式在显示部上一面显示，
上述显示控制部在分割后的上述分析图像上设定横穿左心室心肌的线ROI，基于上述线ROI上的信号轮廓，判定在上述线ROI上有无心肌信号的缺失，在判定为有上述心肌信号的缺失的情况下，在显示上述分析图像时，作为参考信息，对有上述缺失的意思的信息进行显示。
2.根据权利要求1所述的磁共振成像装置，其中，
作为上述心脏的切片的数据，上述序列控制部收集从心尖部到心基底部之间的多个切片的数据。
3.根据权利要求1所述的磁共振成像装置，其中，
作为上述分析图像，上述显示控制部将根据各段的信号强度分配了不同的浓度的颜色的图像、或者根据各段的信号强度分配了不同的颜色的图像排列成矩阵状，并由一面显示。
4.根据权利要求1所述的磁共振成像装置，其中，
作为上述数据的分析结果，上述显示控制部将按照每段对切片的数据进行分析而制成的、表示与所述移动距离对应的时间的方向的信号强度的变化的灌流曲线的曲线图显示于显示部。
5.根据权利要求4所述的磁共振成像装置，其中，
当接受了操作者对显示于显示部的图像进行的选择操作时，上述显示控制部将与由上述操作者选择的图像对应的上述灌流曲线的曲线图显示于显示部。
6.根据权利要求1所述的磁共振成像装置，其中，
作为上述数据的分析结果，上述显示控制部显示将切片的数据按照每段进行分析，根据各段的信号强度分配不同的浓度或者不同的颜色，并将同心圆的半径方向作为与所述移动距离对应的时间的方向而制成的靶心图。
CN201380043226.3A 2012-09-14 2013-09-17 磁共振成像装置 Active CN104582567B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US13/618,001 US10368777B2 (en) 2012-09-14 2012-09-14 Non-contrast dynamic MRI myocardial perfusion analysis and visualization
US13/618,001 2012-09-14
PCT/JP2013/075056 WO2014042275A1 (ja) 2012-09-14 2013-09-17 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
CN104582567A CN104582567A (zh) 2015-04-29
CN104582567B true CN104582567B (zh) 2020-10-16

Family

ID=50275170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380043226.3A Active CN104582567B (zh) 2012-09-14 2013-09-17 磁共振成像装置

Country Status (5)

Country Link
US (1) US10368777B2 (zh)
EP (1) EP2886052B1 (zh)
JP (1) JP6381879B2 (zh)
CN (1) CN104582567B (zh)
WO (1) WO2014042275A1 (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6091989B2 (ja) * 2013-05-10 2017-03-08 東芝メディカルシステムズ株式会社 画像処理装置及び磁気共鳴イメージング装置
KR101623834B1 (ko) * 2013-08-09 2016-05-24 삼성전자주식회사 의료 영상 촬영과 관련된 컨텐츠를 제공하기 위한 방법 및 그 장치
US9833165B2 (en) * 2014-04-29 2017-12-05 Biosense Webster (Israel) Ltd. Checking for perforation of the epicardium using magnetic resonance imaging
US9949643B2 (en) * 2014-10-18 2018-04-24 International Business Machines Corporation Automatic visualization of regional functional parameters of left ventricle from cardiac imaging
JP6953121B2 (ja) * 2016-09-23 2021-10-27 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
CN108459286A (zh) * 2017-02-20 2018-08-28 冯原 基于磁共振成像的体内软组织力学特性测试方法及装置
CN107621616B (zh) * 2017-08-01 2019-12-13 国家纳米科学中心 一种磁共振投影成像方法及装置
WO2019090418A1 (en) 2017-11-12 2019-05-16 Synex Medical Inc. Wearable blood analyte measurement device and method for measuring blood analyte concentration
EP3599586B1 (de) * 2018-07-24 2020-10-14 Siemens Healthcare GmbH Verfahren und vorrichtung zur analyse von magnetresonanz-aufnahmen
CN110133556B (zh) * 2019-05-29 2021-01-19 上海联影医疗科技股份有限公司 一种磁共振图像处理方法、装置、设备及存储介质
CN113285876B (zh) * 2020-02-19 2024-04-23 中兴通讯股份有限公司 路由方法、路由装置及计算机可读存储介质

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5377681A (en) 1989-11-13 1995-01-03 University Of Florida Method of diagnosing impaired blood flow
US5431161A (en) * 1993-04-15 1995-07-11 Adac Laboratories Method and apparatus for information acquistion, processing, and display within a medical camera system
US5647360A (en) * 1995-06-30 1997-07-15 Siemens Corporate Research, Inc. Digital subtraction angiography for 3D diagnostic imaging
US5908386A (en) * 1995-12-14 1999-06-01 Regents Of The Universotiy Of Minnesota Fast MRI for assessment of myocardial perfusion with arrythmia insensitive magnetization preparation
US5827187A (en) * 1996-04-23 1998-10-27 Mayo Foundation For Medical Education And Research Dynamic MR digital subtraction angiography with complex subtraction
US6564080B1 (en) * 1999-03-31 2003-05-13 Kabushiki Kaisha Toshiba MR imaging on ASL technique
DE10064768B4 (de) * 2000-12-22 2006-12-07 Siemens Ag Verfahren zur Untersuchung eines Lebewesens mittels eines bildgebenden Verfahrens
US7047060B1 (en) 2001-11-26 2006-05-16 Koninklijke Philips Electronics N.V. Multiple preparatory excitations and readouts distributed over the cardiac cycle
US7289841B2 (en) * 2002-10-25 2007-10-30 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
JP4446049B2 (ja) 2003-11-06 2010-04-07 株式会社三重ティーエルオー 心筋血流の定量化装置
DE102004038670B4 (de) * 2004-08-09 2014-06-26 Siemens Aktiengesellschaft Verfahren zur Segmentierung eines medizinischen Datensatzes
JP4763989B2 (ja) 2004-09-22 2011-08-31 株式会社東芝 磁気共鳴イメージング装置、磁気共鳴データ処理装置及び磁気共鳴データ処理プログラム
DE102005002949A1 (de) 2005-01-21 2006-08-03 Siemens Ag Verfahren zur Visualisierung von Schädigungen im Myokard
WO2007092159A2 (en) 2006-02-02 2007-08-16 Wake Forest University Health Sciences Cardiac visualization systems for displaying 3-d images of cardiac voxel intensity distributions with optional physician interactive boundary tracing tools
US7587232B2 (en) 2006-02-28 2009-09-08 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus, magnetic resonance data processing apparatus, magnetic resonance data processing program and magnetic resonance imaging apparatus control method
DE102006020864A1 (de) 2006-05-04 2007-11-08 Siemens Ag Verfahren zur Ermittlung und Anzeige von wenigstens einer Information über ein Zielvolumen
US8046045B2 (en) * 2006-05-08 2011-10-25 Siemens Medical Solutions Usa, Inc. Anatomical feature condition assessment system using tissue displacement characteristics
US8532741B2 (en) * 2006-09-08 2013-09-10 Medtronic, Inc. Method and apparatus to optimize electrode placement for neurological stimulation
US8160676B2 (en) * 2006-09-08 2012-04-17 Medtronic, Inc. Method for planning a surgical procedure
US10295638B2 (en) 2007-06-28 2019-05-21 Toshiba Medical Systems Corporation Image processing apparatus, image diagnostic apparatus and image processing method
US8571288B2 (en) * 2007-12-07 2013-10-29 Kabushiki Kaisha Toshiba Image display apparatus and magnetic resonance imaging apparatus
CN102231963B (zh) 2008-12-03 2014-09-24 皇家飞利浦电子股份有限公司 重新参数化的牛眼图
US8581582B2 (en) * 2009-09-18 2013-11-12 Kabushiki Kaisha Toshiba MRI non-contrast time-slip angiography using variably positioned cine sub-sequence
CN102860827B (zh) 2009-09-18 2017-05-17 东芝医疗系统株式会社 磁共振成像装置以及磁共振成像方法
US8848990B2 (en) * 2010-09-28 2014-09-30 Siemens Aktiengesellschaft Automatic registration of image series with varying contrast based on synthetic images derived from intensity behavior model
US20120078085A1 (en) 2010-09-29 2012-03-29 Siemens Corporation Method of Analysis for Dynamic Magnetic Resonance Perfusion Imaging
JP5984239B2 (ja) 2010-11-15 2016-09-06 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
US9566014B2 (en) * 2011-11-09 2017-02-14 Siemens Healthcare Gmbh System for cardiac MR and MR cine imaging using parallel image processing

Also Published As

Publication number Publication date
EP2886052A1 (en) 2015-06-24
WO2014042275A1 (ja) 2014-03-20
EP2886052A4 (en) 2016-04-20
JP6381879B2 (ja) 2018-08-29
JP2014057862A (ja) 2014-04-03
US20140081125A1 (en) 2014-03-20
US10368777B2 (en) 2019-08-06
CN104582567A (zh) 2015-04-29
EP2886052B1 (en) 2018-07-18

Similar Documents

Publication Publication Date Title
CN104582567B (zh) 磁共振成像装置
US9474455B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5575722B2 (ja) 磁気共鳴イメージング装置
US10169866B2 (en) Medical image processing and diagnostic image generation device for predetermined types of diagnostic information
US9585575B2 (en) Non-contrast MRI with differentiation of ischemic, infarct and normal tissue
JP6261871B2 (ja) 画像処理装置及び磁気共鳴イメージング装置
JP6218436B2 (ja) 磁気共鳴イメージング装置
WO2012067123A1 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
US8315450B2 (en) Method and system for display of medical image data
JP2017140132A (ja) 画像処理装置およびｍｒｉ装置
JP6991728B2 (ja) 画像処理装置、磁気共鳴イメージング装置及び画像処理方法
JP6517031B2 (ja) 医用画像処理装置および磁気共鳴イメージング装置
JP5984239B2 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2009028511A (ja) 画像処理装置、画像診断装置および画像処理方法
US11112477B2 (en) Magnetic resonance imaging apparatus and image processing apparatus
JP7237612B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP5002214B2 (ja) 磁気共鳴イメージング装置
US11497412B2 (en) Combined oxygen utilization, strain, and anatomic imaging with magnetic resonance imaging
JP2011024722A (ja) 磁気共鳴イメージング装置及びその画像処理方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20160712

Address after: Japan Tochigi Tian Yuan City

Applicant after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Applicant before: Toshiba Corp

Applicant before: Toshiba Medical System Co., Ltd.

GR01 Patent grant
GR01 Patent grant