CN1676097A - 磁共振成象装置及其磁共振成象装置的数据处理方法 - Google Patents

磁共振成象装置及其磁共振成象装置的数据处理方法 Download PDF

Info

Publication number
CN1676097A
CN1676097A CNA2005100717180A CN200510071718A CN1676097A CN 1676097 A CN1676097 A CN 1676097A CN A2005100717180 A CNA2005100717180 A CN A2005100717180A CN 200510071718 A CN200510071718 A CN 200510071718A CN 1676097 A CN1676097 A CN 1676097A
Authority
CN
China
Prior art keywords
data
sensitivity
signal
aforementioned
map diagram
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
CNA2005100717180A
Other languages
English (en)
Other versions
CN100444791C (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
Priority claimed from JP2004052312A external-priority patent/JP4723814B2/ja
Priority claimed from JP2004052313A external-priority patent/JP4664605B2/ja
Application filed by Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN1676097A publication Critical patent/CN1676097A/zh
Application granted granted Critical
Publication of CN100444791C publication Critical patent/CN100444791C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • 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/24Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246Spatial mapping of the RF magnetic field B1
    • 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/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field

Abstract

磁共振成象装置(20)具有执行用于生成RF线圈(24)的敏感度映象图数据的扫描作业的扫描实施组件;对扫描作业获得的图象数据的位于无信号区域附近位置处的有信号区域,实施区域抑制处理的区域抑制处理组件(44b);使用实施区域抑制处理后的图象数据,生成敏感度映象图数据的敏感度映象图数据生成组件(44);以及时敏感度映象图数据实施三维平滑滤波处理的平滑处理组件(44i)。

Description

磁共振成象装置及其 磁共振成象装置的数据处理方法
技术领域
本发明涉及利用核磁共振信号对被检测体图象实施摄象的磁共振成象装置以及磁共振成象装置的数据处理方法,特别涉及能够对信号接收用线圈的敏感度分布产生的图象数据亮度和信号值的不均匀性实施修正的磁共振成象装置以及磁共振成象装置的数据处理方法。
背景技术
在先技术中作为医疗现场使用的检测装置,包括如图24所示的磁共振成象(MRI:Magnetic Resonance Imaging)装置1(比如说可以参见日本专利第3135592号公报)。
磁共振成象装置1是以下这样的装置:可以通过在位于形成静磁场用的筒形静磁场用磁铁2的内部处的被检测体P的摄象区域处,通过倾斜磁场线圈组件3的各倾斜磁场线圈3x、3y、3z形成沿X轴、Y轴、Z轴方向的倾斜磁场,并且通过RF(Radio Frequency)线圈4对拉莫尔(Larmor)频率的高频(RF)信号实施传送,来使位于被检测体P内部处的原子核自旋产生磁共振,并且利用激励所产生的核磁共振(NMR:Nuclear Magnetic Resonance)信号对被检测体P的图象实施再构成处理。
换句话说就是,可以预先使用静磁场用电源5在静磁场用磁铁2的内部形成静磁场。根据由输入装置6接收到的指令信号,由序列信号控制器控制组件7a将作为信号控制信息的序列信号传送至序列信号控制器8处,进而由序列信号控制器8依据序列信号,对将高频信号传送至与各倾斜磁场线圈3x、3y、3z相连接的倾斜磁场用电源9和RF线圈4处的信号发送器10实施控制。采用这种构成形式,可以在摄象区域处形成倾斜磁场,并且可以向被检测体P发送出高频信号。
由倾斜磁场线圈3x、3y、3z在这时产生的X轴倾斜磁场、Y轴倾斜磁场、Z轴倾斜磁场,主要是分别作为相位编码(PE:phase encoding)用倾斜磁场、读取(RO:readout)用倾斜磁场、层面编码(SE:slice encoding)用倾斜磁场使用的。因此,作为原子核位置信息的X坐标、Y坐标、Z坐标,分别被变换为原子核自旋的相位、频率、层面的位置,并且一边改变相位编码量,一边重复进行序列作业。
然后,伴随着对位于被检测体P内部处的原子核自旋实施激励所产生的NMR信号,利用RF线圈4接收信号,同时传送至信号接收器11处并变换为数字化的原始数据(raw data)。该原始数据可以通过序列信号控制器8读取至序列信号控制器控制组件7a处,进而由序列信号控制器控制组件7a将原始数据配置在形成在原始数据用数据库7b处的K空间(傅立叶空间)处。图象再构成组件7c可以通过对配置在K空间处的原始数据实施傅立叶变换,来获得被检测体P的再构成图象数据,并将其储存在图象数据用数据库7d处。而且,还可以通过显示处理组件7e将图象数据传送至适当的显示装置7f处,以对其实施显示。
这种磁共振成象装置1为了能够提高摄象速度,采用着RF线圈4由发送信号用的整体(WB:whole-body)线圈和作为接收信号用的主线圈的相控矩阵线圈(PAC:phased-array coil)构成的构成形式(比如可以参见Roemer PB,et al,The MNR Phased Array,MRM 16,192-225(1990))。相控矩阵线圈具有多个表面线圈,所以可以通过使用各表面线圈同时接收NMR信号以在短时间里收集更多原始数据的方式,缩短摄象所需要的时间。
当采用这种磁共振成象装置1实施图象诊断时,不希望在最终获得的图象数据中出现亮度起伏(unevenness)(信号强度起伏)。然而,当RF线圈4是由具有多个表面线圈的相控矩阵线圈构成时,由于构成RF线圈4用的各表面线圈的敏感度不均匀性,会使NMR信号的信号强度和简单对原始数据实施傅立叶变换等再构成处理获得的图象数据的信号强度产生不均匀性,所以会使图象数据产生亮度起伏。
因此,在先技术中是在实施生成被检测体P图象用的主扫描之前,先进行敏感度预扫描的。所以,通过敏感度预扫描利用WB线圈和相控矩阵线圈获得图象数据,并通过图25所示的流程步骤,依据作为各图象数据的信号强度SPAC、SWB的除法运算值的信号强度比(SPAC/SWB),对作为三维敏感度映象图数据的相控矩阵线圈的敏感度分布实施推定,进而通过所获得的三维敏感度映象图数据对图象数据的亮度实施修正。
首先,通过敏感度预扫描运行组件7g将敏感度推定用序列信号传送至序列信号控制器控制组件7a处,实施敏感度预扫描作业。将通过WB线圈获得的WB再构成图象和通过相控矩阵线圈获得的PAC再构成图象,分别储存在WB再构成图象用数据库7h和PAC再构成图象用数据库7i处。
随后,由敏感度分布推定组件7j依据WB再构成图象和PAC再构成图象,求解相控矩阵线圈的敏感度分布的推定值。
换句话说就是,在步骤S1中,由阈值处理部7k对WB再构成图象和PAC再构成图象实施阈值处理。对WB再构成图象和PAC再构成图象的各信号强度小于等于阈值(包含该阈值)的区域实施屏蔽处理(mask),并且生成出WB绝对值图象数据和PAC绝对值图象数据。
随后在步骤S2中,由区域抑制处理(region reduction)部71对实施PAC再构成图象和WB再构成图象的阈值处理后的区域实施区域抑制处理,将位于缩小后的屏蔽区域附近处的信号强度比较小的部分,从三维敏感度映象图数据制作时使用的数据中除去。
随后在步骤S3中,由除法运算处理部7m通过将阈值处理和区域抑制处理后的各层面中的PAC绝对值图象(PAC再构成图象的信号绝对值),除以WB绝对值图象(WB再构成图象的信号绝对值)的方式,对作为三维敏感度映象图数据的PAC绝对值图象和WB绝对值图象间的信号强度比实施求解。
随后在步骤S4中,通过标准化处理部7n对作为PAC绝对值图象和WB绝对值图象间的信号强度比求解出的三维敏感度映象图数据,按每个层面实施标准化处理。
随后在步骤S5中,由被检测体内区域内插处理部7o对实施标准化处理后的三维敏感度映象图数据中位于被检测体P内部处的无信号区域,实施线型内插处理。
随后在步骤S6中,为了对位于被检测体外部处无信号区域中的敏感度分布实施推定,通过被检测体外区域内插处理部7p实施区域增大(regiongrowing)处理,由此对被检测体外部处无信号区域中的敏感度分布实施内插。
随后在步骤S7中,为了使沿切割层面方向的敏感度分布一样,而由层面方向加权处理部7q对三维敏感度映象图数据沿切割层面方向加权,由此对其实施修正处理。
随后在步骤S8中,由平滑处理部7r实施诸如数据拟合处理和平滑(smoothing)处理等各种处理,以制作出作为整个三维区域中的体积数据的敏感度映象图,并储存在敏感度映象图用数据库7s处。
随后,可以将图象获得用序列信号由主扫描运行组件7t传送至序列信号控制器控制组件7a处,运行主扫描。由图象再构成组件7c实施图象再构成处理,并将所获得的图象数据储存在图象数据用数据库7d处。图象数据修正组件7u利用储存在敏感度映象图用数据库7s处的敏感度映象图,对储存在图象数据用数据库7d处的图象数据的亮度实施修正,并且通过显示处理组件7e将实施亮度修正后的图象数据传送至显示装置7f处,以实施显示。
一般说来,通过敏感度预扫描和主扫描获得的图象数据中包含有无信号区域。在作为摄象区域的被检测体处存在有诸如肺等的部位,所以存在有不会产生NMR信号的无信号区域。而且,位于被检测体内的无信号区域与有信号区域间的边界附近位置处的有信号区域,通常会出现信号强度比较小的现象。而且,这会对相控矩阵线圈的敏感度分布推定值产生影响。
然而,属于在先技术的磁共振成象装置1,在制作相控矩阵线圈的敏感度映象图时的无信号区域内插方法,与位于被检测体内外无关地仅仅对无信号区域实施区域增大处理的。
而且,使用相控矩阵线圈实施的敏感度预扫描与使用WB线圈实施的敏感度预扫描是分别实施的,所以在使用相控矩阵线圈获得的图象数据中的被检测体P的形状,与使用WB线圈获得的图象数据中的被检测体P的形状间往往会出现偏差。
而且,不只是相控矩阵线圈处的各表面线圈的敏感度分布,还会由于配置特性而产生与装置坐标系统中Z轴方向相关的信号强度不均匀性。
然而,属于在先技术的磁共振成象装置1,并没有考虑到实施敏感度预扫描时的图象数据位置偏差和表面线圈的配置特性产生的信号强度不均匀性,仅仅是对WB再构成图象数据和PAC再构成图象数据的各信号强度实施阈值处理,并依据由此分别获得的PAC绝对值图象数据和WB绝对值图象数据间的信号强度比,对相控矩阵线圈的敏感度分布实施推定并生成敏感度映象图。
由于采用这种构成形式,属于在先技术的磁共振成象装置1,难以使通过敏感度预扫描推定出的相控矩阵线圈的敏感度分布具有良好精度,所以存在有难以对实施主扫描获得的再构成图象的亮度实施具有良好精度的修正的问题。
而且,在先技术中作为医疗现场使用的检测装置,还包括如图26所示的磁共振成象装置1(比如说可以参见日本专利第3135592号公报)。
如上所述,当RF线圈4是由相控矩阵线圈和WB线圈构成时,与相控矩阵线圈和WB线圈的敏感度不均匀性相关的NMR信号,以及通过再构成处理方式获得的图象数据的信号强度也会出现不均匀性。一般说来,WB线圈的敏感度不均匀性小到可以被忽略的程度,然而作为为了其它技术目的而设置的线圈的相控矩阵线圈中的表面线圈的敏感度不均匀性则比较大,会对图象数据产生影响。
因此,有必要对由于相控矩阵线圈的敏感度不均匀性而产生的图象数据信号强度中的不均匀性实施修正。
如图26所示的磁共振成象装置1,可以按照如图27的示意性流程图所示,在步骤S1中通过敏感度预扫描运行组件7g,将敏感度推定用序列信号传送至序列信号控制器控制组件7a处,将相控矩阵线圈和WB线圈作为信号接收用线圈并实施敏感度预扫描作业。将通过WB线圈获得的WB线圈图象数据和通过相控矩阵线圈获得的主线圈图象数据,作为相控矩阵线圈的敏感度分布推定用的图象数据实施获取,并分别储存在WB线圈图象用数据库7h和主线圈图象用数据库7v处。采用这种构成形式,可以按两次方式对作为三维图象数据的体积数据实施摄象。
随后在步骤S2中,通过敏感度分布推定组件7j对相控矩阵线圈的敏感度分布推定值实施求解。换句话说就是,通过使用除法运算处理部7w,利用如图28(b)所示的WB线圈图象数据的信号强度SWB,对如图28(a)所示的主线圈图象数据的信号强度SPAC实施除运算的方式,求出作为相控矩阵线圈的敏感度分布推定值的、如图28(c)所示的主线圈图象数据与WB线圈图象数据间的信号强度比(SPAC/SWB)。
在这时,为了不对主线圈图象数据和WB线圈图象数据的信号强度SPAC、SWB小于阈值(不包含该阈值)、比如说小于最大值的10%(不包括10%)的区域,实施除法运算处理,可以使用阈值处理部7x根据需要实施作为除法运算前处理的各信号强度SPAC、SWB的阈值处理,对位于小于阈值(不包含该阈值)的区域处的各信号强度SPAC、SWB实施屏蔽。
通过实施将这种阈值处理作为前处理的除法运算处理,可以消除诸如图象对比度等的、相控矩阵线圈的敏感度分布不均匀性之外的主要因素对图象数据的信号强度产生的影响,从而可以按照具有良好精度的方式对敏感度分布实施推定。
随后,对由于肺叶等的存在而在实施阈值处理时产生的数据失落部分的无信号区域,使用内插处理部7y实施内插处理或外插处理,对敏感度分布实施推定,进而通过平滑处理部7z对整个二维区域实施数据拟合处理和平滑处理,获得如图28(d)所示的敏感度分布推定值曲线。
然后,对整个三维区域的各个剖面实施同样的图象数据处理,求解出作为体积数据的敏感度分布推定值。
随后在步骤S3中,将作为三维敏感度映象图数据的相控矩阵线圈的敏感度分布推定值,储存在敏感度映象图用数据库7s处。
随后在步骤S4中,将图象获得用序列信号由主扫描运行组件7t传送至序列信号控制器控制组件7a处,在将相控矩阵线圈作为信号接收用线圈的条件下运行主扫描。对原始数据实施收集并通过图象再构成组件7c实施的图象再构成处理,获得图象数据。
随后在步骤S5中,按照与实施主扫描时诸如摄象剖面方向、空间分辨率等摄象条件、数据收集条件、图象再构成条件等的各条件相对应的方式,通过图象数据修正组件7u从敏感度映象图用数据库7s中获取出相对应的三维敏感度映象图数据。
随后在步骤S6中,图象数据修正组件7u利用所获取出的三维敏感度映象图数据,对图象数据实施修正处理。采用这种方式,可以改善图象数据的信号强度不均匀性。
在另一方面,还采用以下的方法:通过对使用相控矩阵线圈接收到的NMR信号获得的图象数据自身实施后处理,来对作为敏感度映象图数据的相控矩阵线圈的敏感度分布实施推定,进而利用所获得的相控矩阵线圈的敏感度映象图数据,对图象数据的信号强度起伏实施修正。如果举例来说,可以采用通过对使用相控矩阵线圈获得的图象数据实施平滑处理的方式,制作出具有极低频率成分的图象数据,并将其作为敏感度分布使用的方法。
而且,还提出以下方案:以从WB线圈实施信号发送的高频信号的信号强度作为参考,对由通过相控矩阵线圈接收到的NMR信号获得的图象数据的信号强度起伏实施修正的技术(比如说可以参见日本特开昭63-132645号公报),以及利用预先储存的相控矩阵线圈的敏感度分布和从图象数据实施推定所获得的相控矩阵线圈的位置信息,对图象数据的信号强度起伏实施修正的技术(比如说可以参见日本特开平7-59750号公报)。
然而,如果采用在先技术中依据实施敏感度预扫描时分别通过相控矩阵线圈和WB线圈获得的图象数据的信号强度值SigPAC、SigWB的除法运算值(SigPAC/SigWB),对相控矩阵线圈的敏感度分布实施推定,由此对图象数据的信号强度实施修正的方法时,会出现实施敏感度预扫描所需要的时间比较长的问题。因此,对于对例如被检测体P的腹部实施摄象的场合,与敏感度预扫描的运行时间相对应地暂停呼吸时间也会比较长。
而且,在使用相控矩阵线圈实施图象数据收集与使用WB线圈实施图象数据收集之间,有可能由于被检测体P移动等原因而导致被检测体P的位置偏离(错误记录)。而且,为了能够通过相控矩阵线圈和WB线圈两者获取图象数据,还需要使相控矩阵线圈和WB线圈间形成完全去耦(decoupling)的构成形式。
在另一方面,在使用对由相控矩阵线圈给出的图象数据自身实施后处理推定出的相控矩阵线圈的敏感度分布,对图象数据的信号强度起伏实施修正的方法中,由于对相控矩阵线圈的敏感度分布的推定精度比较低,所以难以对图象数据的信号强度起伏实施充分修正,从而难以使最终获得的修正后的图象数据具有充分的均匀性。
而且,还存在有通常难以对各种各样的图象种类实施精度良好的修正处理的问题。如果举例来说,在通过对图象数据进行T1强调和T2强调的方式,使图象数据为具有所希望的对比度的图象数据的场合,即使能够通过平滑处理而对敏感度分布实施推定,由于表示敏感度分布的图象数据具有对比度,所以也难以使用实施平滑处理后的图象数据作为敏感度分布使用。
而且,还一直存在有通过多重层面方式实施摄象时,难以对全部层面实施图象数据修正的问题。
发明内容
本发明就是解决上述技术问题用的发明,本发明的目的就是提供一种能够依据通过实施敏感度预扫描获得的图象数据,对RF线圈的敏感度分布实施更高精度的推定,从而可以依据所获得的RF线圈的敏感度分布,对实施主扫描而获得的图象数据的亮度实施良好修正的磁共振成象装置,以及磁共振成象装置的数据处理方法。
本发明的另一目的就是提供一种能够在更短的时间里,按照与图象种类等摄象条件无关的方式,对信号接收用线圈的敏感度分布的不均匀性造成的图象数据的信号强度起伏实施精度良好的修正的磁共振成象装置,以及磁共振成象装置的数据处理方法。
为了能够实现上述目的,本发明提供的一种磁共振成象装置具有:执行用于生成RF线圈的敏感度映象图数据用的扫描作业的扫描实施组件;对前述扫描作业获得的图象数据的无信号区域附近位置处的有信号区域,实施区域抑制处理的区域抑制处理组件;使用实施区域抑制处理后的前述图象数据,生成敏感度映象图数据的敏感度映象图数据生成组件;以及对前述敏感度映象图数据实施三维平滑滤波处理的平滑处理组件。
而且,为了能够实现上述目的,本发明提供的另一种磁共振成象装置具有:执行用于生成RF线圈的敏感度映象图数据用的扫描作业的扫描实施组件;使用通过前述扫描作业获得的图象数据,生成敏感度映象图数据的敏感度映象图数据生成组件;对前述敏感度映象图数据的被检测体内部处的无信号区域实施线性内插处理的线性内插处理组件;以及对前述敏感度映象图数据实施三维平滑滤波处理的平滑处理组件。
而且,为了能够实现上述目的,本发明提供的再一种磁共振成象装置具有:将实施上述图象摄象时的信号接收用线圈作为信号接收用线圈,执行用于生成前述图象摄象的信号接收用线圈的敏感度修正数据的扫描作业的扫描实施组件;以及仅仅将前述扫描作业中通过前述图象摄象时的信号接收用线圈获得的数据作为原始数据,生成前述敏感度修正用数据的敏感度修正用数据生成组件。
而且,为了能够实现上述目的,本发明提供的一种磁共振成象装置的数据处理方法具有:对实施用于生成RF线圈的敏感度映象图数据的扫描而获得的图象数据的位于无信号区域附近的有信号区域,实施区域抑制处理的步骤;使用实施区域抑制处理后的前述图象数据,生成敏感度映象图数据的步骤;以及对前述敏感度映象图数据实施三维平滑滤波处理的步骤。
而且,为了能够实现上述目的,本发明提供的另一种磁共振成象装置的数据处理方法具有:使用实施用于生成RF线圈的敏感度映象图数据的扫描而获得的图象数据,生成敏感度映象图数据的步骤;对前述敏感度映象图数据的位于被检测体内部的无信号区域实施线性内插处理的步骤;以及对前述敏感度映象图数据实施三维平滑滤波处理的步骤。
而且,为了能够实现上述目的,本发明提供的再一种磁共振成象装置的数据处理方法具有:将实施图象摄象时的信号接收用线圈作为信号接收用线圈,执行用于生成前述图象摄象的信号接收用线圈的敏感度修正用数据的扫描的步骤;以及仅仅将前述扫描作业中通过前述图象摄象的信号接收用线圈获得的数据作为原始数据,生成前述敏感度修正用数据的步骤。
根据本发明的磁共振成象装置和磁共振成象装置的数据处理方法,能够依据实施敏感度预扫描而获得的图象数据,以更高精度对RF线圈的敏感度分布实施推定,进而能够依据所获得的RF线圈的敏感度分布,对通过实施主扫描而获得的图象数据的亮度实施更为良好的修正处理。
而且,还可以在更短的时间里,与图象种类等摄象条件无关地,对信号接收用线圈的敏感度分布的不均匀性造成的图象数据的信号强度起伏实施精度良好的修正。
附图说明
图1为表示本发明的磁共振成象装置的一种实施形式的功能方框图。
图2为表示图1中所示的RF线圈和信号接收器的详细构成图。
图3为表示图2中所示的WB线圈和相控矩阵线圈的一个配置实例的剖面模式图。
图4为表示使用图1中所示的磁共振成象装置对被检测体实施断层图象摄象时的步骤的流程图。
图5为表示使用图1中所示的磁共振成象装置生成敏感度映象图数据时的详细步骤的流程图。
图6为表示通过图1中所示的磁共振成象装置实施区域抑制处理的一个实例的图。
图7为表示通过图1中所示的磁共振成象装置对三维敏感度映象图数据的位于被检测体内部中的无信号区域实施线性内插处理的一个实例用的图。
图8为说明通过26点方法对三维敏感度映象图数据的位于被检测体外部的无信号区域实施区域增大处理的场合的说明图。
图9为说明通过6点方法对三维敏感度映象图数据的位于被检测体外部的无信号区域实施区域增大处理的场合的说明图。
图10为表示通过图1中所示的磁共振成象装置对三维敏感度映象图数据的位于被检测体外部的无信号区域实施区域增大处理的一个实例的图。
图11为表示使用图1中所示的磁共振成象装置对三维敏感度映象图数据沿切割层面方向实施加权处理时的修正系数的曲线图。
图12为表示使用图1中所示的磁共振成象装置获得的亮度修正后的被检测体的断层图象的图。
图13为表示使用在先技术中的磁共振成象装置获得的亮度修正后的被检测体的断层图象的图。
图14为表示本发明的磁共振成象装置的第二实施形式的功能方框图。
图15为表示图14中所示的RF线圈和信号接收器的一个实例的详细构成图。
图16为表示使用图14中所示的磁共振成象装置对被检测体实施断层图象摄象时的步骤的流程图。
图17为表示在图16中所示的流程图中,对相控矩阵线圈的敏感度分布实施推定时的详细步骤的一个实例的流程图。
图18为表示在低对比度摄象条件下获得的主线圈图象数据的信号强度分布的图。
图19为表示在常规对比度摄象条件下获得的主线圈图象数据的信号强度分布的图。
图20为表示使用在先技术中的磁共振成象装置实施三维敏感度映象图数据生成时所生成出的数据的图。
图21为说明使用图14中所示的磁共振成象装置实施三维敏感度映象图数据生成时所存在问题的说明图。
图22为表示通过能够避免如图21所示问题的步骤,使用磁共振成象装置作为三维敏感度映象图数据时所生成出的数据的图。
图23为表示本发明的磁共振成象装置的第三实施形式的功能方框图。
图24为表示在先技术中的一种磁共振成象装置的功能方框图。
图25为表示通过图24中所示的在先技术的磁共振成象装置生成敏感度映象图的步骤的流程图。
图26为表示在先技术中的一种磁共振成象装置的功能方框图。
图27为表示通过图26中所示的在先技术的磁共振成象装置对图象数据信号强度起伏实施修正的步骤的流程图。
图28为说明通过图26中所示的在先技术的磁共振成象装置推定敏感度分布的步骤的说明图。
具体实施方式
下面参考附图,对本发明的磁共振成象装置以及磁共振成象装置的数据处理方法的实施形式进行说明。
图1为表示本发明的磁共振成象装置的第一实施形式的功能方框图。
磁共振成象装置20构成为:将形成静磁场用的筒形静磁场用磁铁21、设置在该静磁场用磁铁21内部处的均场线圈(shim coil)22、以及倾斜磁场线圈组件23和RF线圈24,安装在图中未示出的机架中。
而且,在磁共振成象装置20中还设置有控制系统25。控制系统25具有静磁场用电源26、倾斜磁场用电源27、均场线圈用电源28、信号发送器29、信号接收器30、序列信号控制器31和计算机32。控制系统25中的倾斜磁场用电源27由X轴倾斜磁场用电源27x、Y轴倾斜磁场用电源27y和Z轴倾斜磁场用电源27z构成。而且,计算机32具有图中未示出的运算装置和储存装置,并且设置有输入装置33和显示装置34。
静磁场用磁铁21与静磁场用电源26相连接,并且具有能够通过由静磁场用电源26供给的电流在摄象区域形成静磁场的功能。而且,在静磁场用磁铁21的内侧处,还同轴地设置着筒型的均场线圈22。均场线圈22与均场线圈用电源28相连接,并且构成为能够由均场线圈用电源28向均场线圈22供给电流而使静磁场均匀化。
倾斜磁场线圈组件23由X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z构成,并且在静磁场用磁铁21的内部形成为筒状。在倾斜磁场线圈组件23的内部设置有躺卧平台35并成为摄象区域,并且可以将被检测体P设置在躺卧平台35上。而且,也可以采用不将RF线圈24设置在机架中,而是设置在躺卧平台35和被检测体P附近的构成形式。
倾斜磁场线圈组件23与倾斜磁场用电源27相连接。倾斜磁场线圈组件23中的X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z,分别与倾斜磁场用电源27中的X轴倾斜磁场用电源27x、Y轴倾斜磁场用电源27y和Z轴倾斜磁场用电源27z相连接。
而且,通过由X轴倾斜磁场用电源27x、Y轴倾斜磁场用电源27y和Z轴倾斜磁场用电源27z分别向X轴倾斜磁场线圈23x、Y轴倾斜磁场线圈23y和Z轴倾斜磁场线圈23z供给电流,可以在摄象区域处分别形成X轴方向的倾斜磁场Gx、Y轴方向的倾斜磁场Gy和Z轴方向的倾斜磁场Gz。
RF线圈24与信号发送器29和信号接收器30相连接。RF线圈24具有从信号发送器29接收高频信号并传送至被检测体P处的功能,并且具有接收随着被检测体P内部的原子核自旋的高频信号形成的激励所产生的NMR信号,并传送至信号接收器30的功能。
图2为表示图1中所示的RF线圈24和信号接收器30的详细构成图。
RF线圈24由例如高频信号发送用的WB线圈24a、NMR信号接收用的相控矩阵线圈24b构成。相控矩阵线圈24b具有多个表面线圈24c,信号接收器30由多个信号接收系统回路30a构成。而且,各表面线圈24c分别独立地与信号接收器30的各信号接收系统回路30a相连接,WB线圈与信号发送器29和信号接收系统回路30a相连接。
图3为表示图2中所示的WB线圈24a和相控矩阵线圈24b的一个配置实例用的剖面模式图。
相控矩阵线圈24b的各表面线圈24c,围绕Z轴对称地配置在包含例如被检测体P的特定关注部位的剖面L的周围。而且,将WB线圈24a设置在相控矩降线圈24b的外侧处。RF线圈24可以通过WB线圈24a向被检测体P发送出高频信号,并且可以利用WB线圈24a或相控矩阵线圈24b的各表面线圈24c,通过多信道接收由包含特定关注部位的剖面L处给出的NMR信号,进而传送至各信号接收器30的各信号接收系统回路30a。
在另一方面,控制系统25的序列信号控制器31,与倾斜磁场用电源27、信号发送器29和信号接收器30相连接。序列信号控制器31具有以下功能:对驱动倾斜磁场用电源27、信号发送器29和信号接收器30所需要的控制信息、描述诸如施加至倾斜磁场用电源27处的脉冲电流的强度和施加时间、施加定时时间等动作控制信息的序列信号信息实施储存的功能;依据所储存的规定的序列信号对倾斜磁场用电源27、信号发送器29和信号接收器30实施驱动,以生成出X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz和高频信号的功能。
而且,序列信号控制器31还可以从信号接收器30接收作为数字化后的NMR信号的原始数据(raw data),并传送至计算机32。
因此,信号发送器29具有依据从序列信号控制器31接收到的控制信息,将高频信号传送至RF线圈24处的功能,信号接收器30具有:对由RF线圈24接收到的NMR信号实施所需要的信号处理并实施A/D变换,由此生成作为数字化后的NMR信号的原始数据的功能;将所生成出的原始数据传送至序列信号控制器31的功能。
而且,计算机32可以通过对程序实施读取和运行的方式,而作为敏感度预扫描运行组件36、主扫描运行组件37、序列信号控制器控制组件38、原始数据用数据库39、图象再构成组件40、图象数据用数据库41、PAC再构成图象用数据库42、WB再构成图象用数据库43、敏感度分布推定组件44、敏感度映象图用数据库45、图象数据修正组件46、显示处理组件47发挥功能。当然,不只是程序,还可以通过设置特定的电路而构成计算机32。
敏感度预扫描运行组件36具有:对为了求出作为相控矩阵线圈24b的敏感度分布的三维敏感度映象图数据,而实施敏感度预扫描时的序列信号(敏感度推定用序列信号)实施生成的功能;将所生成出的敏感度推定用序列信号,传送至序列信号控制器控制组件38以实施敏感度预扫描的功能。
主扫描运行组件37具有:将实施用于获取图象数据的主扫描时使用的各种序列信号,传送至序列信号控制器控制组件38以实施主扫描的功能。
序列信号控制器控制组件38具有:依据由输入装置33或其它构成要素给出的信息,将从敏感度预扫描运行组件36和主扫描运行组件37接收到的序列信号中的所需要的序列信号,传送至序列信号控制器31处以实施敏感度预扫描或主扫描的功能。而且,序列信号控制器控制组件38还具有接收从序列信号控制器31通过实施敏感度预扫描或主扫描所收集到的WB线圈24a和相控矩阵线圈24b的各表面线圈24c的原始数据,并将其配置在形成在原始数据用数据库39处的K空间(傅立叶空间)中的功能。
采用这种构成形式,原始数据用数据库39可以对信号接收器30中生成出的各个WB线圈24a和表面线圈24c的各原始数据实施储存。换句话说就是,将原始数据配置在形成在原始数据用数据库39处的K空间中。
图象再构成组件40具有:对通过运行主扫描而配置在原始数据用数据库39的K空间处的原始数据,实施诸如傅立叶变换(FT)等图象再构成处理以再构成出被检测体P的图象数据的功能;将再构成出的图象数据写入至图象数据用数据库41的功能。
而且,图象再构成组件40还具有:对通过运行敏感度预扫描而配置在原始数据用数据库39的K空间处的原始数据,按照对运行主扫描而获得的原始数据实施再构成处理相同的方式实施再构成处理,从而根据由相控矩阵线圈24b和WB线圈24a分别获得的被检测体P的图象数据,生成出PAC再构成图象和WB再构成图象的功能;将所生成出的PAC再构成图象和WB再构成图象,分别写入至PAC再构成图象用数据库42和WB再构成图象用数据库43处的功能。
敏感度分布推定组件44具有:使用分别储存在PAC再构成图象用数据库42和WB再构成图象用数据库43中的PAC再构成图象和WB再构成图象,制作出相控矩阵线圈24b的三维敏感度映象图数据的功能;将所制作出的三维敏感度映象图数据写入至敏感度映象图用数据库45处的功能。换句话说就是,敏感度分布推定组件44具有作为生成RF线圈24的敏感度映象图数据的敏感度映象图数据生成组件使用的功能。
因此,敏感度分布推定组件44具有阈值处理部44a、区域抑制处理部44b、除法运算处理部44c、标准化处理部44d、数据平坦化处理部44e、被检测体内区域内插处理部44f、被检测体外区域内插处理部44g、层面方向加权处理部44h、三维平滑处理部44i。
阈值处理部44a具有对PAC再构成图象和WB再构成图象实施阈值处理的功能,即具有对PAC再构成图象和WB再构成图象中的各信号强度小于等于分别预先设定的阈值(包括该阈值)的部分的数据实施屏蔽的功能。
区域抑制处理部44b具有通过对敏感度分布推定时使用的PAC再构成图象和WB再构成图象的阈值处理后的区域,实施区域抑制处理,将位于缩小后屏蔽区域附近处的信号强度比较小的部分,从三维敏感度映象图数据制作时使用的数据中除去的功能。
除法运算处理部44c具有通过将作为阈值处理和区域抑制处理后的PAC再构成图象的信号绝对值的PAC绝对值图象,除以作为WB再构成图象的信号绝对值的WB绝对值图象,求出作为三维敏感度映象图数据的PAC绝对值图象和WB绝对值图象间的信号强度的功能。
标准化处理部44d具有对三维敏感度映象图数据实施标准化处理的功能。
数据平坦化处理部44e具有利用变换函数对三维敏感度映象图数据实施数据平坦化处理,并且利用逆变换函数对数据平坦化处理后的三维敏感度映象图数据实施处理,以求解出实施数据平坦化处理前的三维敏感度映象图数据的功能。换句话说就是,数据平坦化处理部44e具有作为暂时将三维敏感度映象图数据变换为适用线性内插的平坦化分布,并且将线性内插处理后的三维敏感度映象图数据恢复至原始分布的数据平坦化处理组件使用的功能。因此,可以与目的相对应地将任意的函数作为这种变换函数。
被检测体内区域内插处理部44f具有作为对三维敏感度映象图数据的被检测体P内部的无信号区域,实施线性内插处理的线性内插处理组件使用的功能。
被检测体外区域内插处理部44g具有作为对三维敏感度映象图数据的被检测体P外部的无信号区域,通过区域增大处理而实施内插的区域增大处理组件使用的功能。
层面方向加权处理部44h具有能够沿切割层面方向对三维敏感度映象图数据实施加权处理的功能。
三维平滑处理部44i具有作为对三维敏感度映象图数据实施三维平滑滤波处理的平滑处理组件使用的功能。
图象数据修正组件46具有:从储存在敏感度映象图用数据库45处的三维敏感度映象图数据中,分离并抽出与实施主扫描时的摄象条件、数据收集条件、图象再构成条件等的图象数据获取条件相对应的三维敏感度映象图数据的功能;使用所抽取出的三维敏感度映象图数据,对通过运行主扫描而储存在图象数据用数据库41处的图象数据的亮度实施修正的功能。
显示处理组件47具有将储存在图象数据用数据库41处的图象数据,传送至显示装置34处以实施显示的功能。
具有上述构成形式的磁共振成象装置20,作为整体,由各个构成要素作为实施诸如主扫描和敏感度预扫描等扫描的扫描运行组件、对通过敏感度预扫描获得的图象数据的位于无信号区域附近处的信号区域实施区域抑制处理的区域抑制处理组件、使用通过敏感度预扫描获得的图象数据生成出敏感度映象图数据的敏感度映象图数据生成组件、以及通过对敏感度映象图数据沿切割层面方向进行加权而对其实施修正的层面方向加权处理组件而发挥功能。
下面,对磁共振成象装置20的作用进行说明。
图4为表示使用图1中所示的磁共振成象装置20对被检测体P实施断层图象摄象时的步骤的流程图,图中在字母S后附加有数字的参考标号,表示流程中的各步骤。
首先在步骤S10中,通过敏感度预扫描运行组件36将敏感度推定用序列信号传送至序列信号控制器控制组件7a处,在用于获取图象数据的主扫描之前,先进行获取相控矩阵线圈24b的敏感度映象图数据的敏感度预扫描作业。
换句话说就是,将被检测体P预先设置在躺卧平台35上,由静磁场用电源26向静磁场用磁铁21供给电流以在摄象区域形成静磁场。而且,由均场线圈用电源28向均场线圈22供给电流以使形成在摄象区域处的静磁场均匀化。
随后,从输入装置33向序列信号控制器控制组件38发出动作指令。序列信号控制器控制组件38将敏感度推定用序列信号传送至序列信号控制器31。序列信号控制器31依据该敏感度推定用序列信号对倾斜磁场用电源27、信号发送器29和信号接收器30实施驱动,从而在设置有被检测体P的摄象区域处形成X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz,而且同时产生出高频信号。
在这时,由倾斜磁场线圈产生的X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz,主要是分别作为相位编码(PE)用倾斜磁场、读取(RO)用倾斜磁场、层面编码(SE)用倾斜磁场使用。因此,在被检测体P的内部的原子核自旋的转动方向上表现出规则性的、作为利用SE用倾斜磁场沿Z轴方向形成的层面中的二维位置信息的X坐标和Y坐标,被分别变换为通过PE用倾斜磁场和RO用倾斜磁场形成的被检测体P内部的原子核自旋的相位变化量和频率变化量。
与敏感度推定用序列信号相对应地将高频信号,从信号发送器29施加至RF线圈24的WB线圈24a处,并且从WB线圈24a将高频信号传送至被检测体P处。而且,在被检测体P的内部与高频信号的频率相对应的包含在切割层面内的原子核核磁共振产生的NMR信号,通过RF线圈24的WB线圈24a和相控矩阵线圈24b的各表面线圈24c,通过多信道被接收并分别传送至各信号接收器30处。
各信号接收器30从WB线圈24a和相控矩阵线圈24b的各表面线圈24c接收NMR信号,并且对其进行诸如前置放大、中间频率变换、相位检波、低频放大、滤波等各种信号处理。各信号接收器30还通过对NMR信号实施A/D变换的方式,生成出作为数字数据的NMR信号的原始数据。信号接收器30还将所生成出的原始数据传送至序列信号控制器31处。
序列信号控制器31将从信号接收器30接收到的原始数据传送至序列信号控制器控制组件38处,由序列信号控制器控制组件38将原始数据配置在形成在原始数据用数据库39处的K空间中。图象再构成组件40通过分别对由WB线圈24a和相控矩阵线圈24b获得的原始数据实施傅立叶变换(FT),而生成出作为被检测体P的三维图象数据的WB再构成图象和PAC再构成图象,并且分别将其写入至WB再构成图象用数据库43和PAC再构成图象用数据库42中。
在这时,对通过敏感度预扫描收集到的原始数据实施的再构成处理方法,可以是与主扫描中的再构成处理方法相同的方法。如果举例来说,在主扫描中使用相控矩阵线圈24b进行的再构成处理方法,可以是对通过各表面线圈24c获得的图象数据求解平方和的平方根的Sum of Square处理(SoS处理)方法,或对通过各表面线圈24c获得的图象数据求解信号强度绝对值的和的再构成处理方法。
随后在步骤S11中,敏感度分布推定组件44通过利用储存在WB再构成图象用数据库43和PAC再构成图象用数据库42处的各切割层面的WB再构成图象和PAC再构成图象,来生成出三维敏感度映象图数据。
图5为表示使用图1中所示的磁共振成象装置20,生成敏感度映象图数据时的详细步骤的流程图,图中在字母S后附加有数字的参考标号,表示流程中的各步骤。
首先在步骤S20中,由阈值处理部44a对PAC再构成图象和WB再构成图象实施阈值处理,对PAC再构成图象和WB再构成图象的各信号强度小于等于分别预先设定的阈值(包含该阈值)的部分的数据实施屏蔽。因此,由于位于肺叶和被检测体外部区域的信号强度会小于等于阈值(包含该阈值),所以将被认为是无信号区域的区域处的PAC再构成图象和WB再构成图象,从制作三维敏感度映象图数据用的数据中排除。
随后在步骤S21中,由区域抑制处理部44b对敏感度分布推定时使用的PAC再构成图象和WB再构成图象的阈值处理后的区域实施区域抑制处理,将位于缩小后屏蔽区域附近处的信号强度比较小的部分从三维敏感度映象图数据制作时使用的数据中除去。
图6为表示通过图1中所示的磁共振成象装置20实施区域抑制处理的一个实例的图。
实施区域抑制处理之前的PAC再构成图象和WB再构成图象,具有如图6(a)所示的、通过阈值处理实施屏蔽后的无信号区域D1和有信号区域D2。位于无信号区域D1附近位置处的有信号区域D2,通常会出现信号强度比较小的现象。由于使用相控矩阵线圈24b实施的敏感度预扫描和使用WB再构成图象实施的敏感度预扫描是分别实施的,所以对于被检测体P的内脏等产生了位置偏离的场合,如果原封不动地使用PAC再构成图象和WB再构成图象的信号强度比,则在位于无信号区域D1附近位置处的有信号区域D2中会出现信号强度不连续的现象。
因此,将如图6(b)所示的有信号区域D2中的位于无信号区域D1附近位置处的部分D2’置换为无信号区域,从而缩小有信号区域D2。采用这种方式,可以将PAC再构成图象和WB再构成图象中的位于无信号区域D1附近位置处的信号强度比较小的有信号区域D2,从三维敏感度映象图数据制作时使用的数据中除去。
随后在步骤S22中,由除法运算处理部44c通过将实施阈值处理和区域抑制处理后的各切割层面的作为PAC再构成图象的信号绝对值的PAC绝对值图象,除以作为WB再构成图象的信号绝对值的WB绝对值图象,从而求出作为三维敏感度映象图数据的PAC绝对值图象和WB绝对值图象间的信号强度比。
随后在步骤S23中,通过标准化处理部44d对作为PAC绝对值图象和WB绝对值图象间的信号强度比求解出的三维敏感度映象图数据,按每个切割层面实施标准化处理。
随后在步骤S24中,由数据平坦化处理部44e利用变换函数,对实施标准化处理后的三维敏感度映象图数据实施数据平坦化处理,将三维敏感度映象图数据变换为适于进行线性内插的平坦化数据。如果举例来说,可以采用诸如n次函数、指数函数、对数函数等任意函数,对三维敏感度映象图数据实施数据拟合,并按照能够降低对线性内插产生影响的局部数据起伏的方式实施适当处理。
随后在步骤S25中,由被检测体内区域内插处理部44f对实施数据平坦化处理后的三维敏感度映象图数据中的位于被检测体P内部处的无信号区域,实施线性内插处理。
图7为表示通过图1中所示的磁共振减象装置20,对三维敏感度映象图数据的位于被检测体P内部的无信号区域实施线性内插处理的一个实例的图。
图7(a)为表示沿切割层面方向看三维敏感度映象图数据时的一个实例的图,图7(b)为表示沿PE层面方向看三维敏感度映象图数据时的一个实例的图。实施标准化处理和数据平坦化处理后的三维敏感度映象图数据,具有无信号区域D1和有信号区域D2。而且,无信号区域D1由位于被检测体P内部的无信号区域D1a和位于被检测体P外部的无信号区域D1b构成。
被检测体内区域内插处理部44f如图7(b)中的箭头所示,对位于被检测体P内部的无信号区域D1a,通过使用直线对诸如沿RO方向剖面的有信号区域D2中的信号强度值实施连接的方式,实施线性内插处理。采用这种方式,可以将位于被检测体P内部的无信号区域D1a置换为有信号区域D2。
而且,线性内插处理并不仅限于沿RO方向,还可以沿诸如PE方向、SL方向等任意方向进行。
随后在步骤S26中,由数据平坦化处理部44e对实施线性内插处理后的三维敏感度映象图数据,利用逆转换函数将其再次转换为实施数据平坦化处理前的状态。
随后在步骤S27中,由被检测体外区域内插处理部44g对三维敏感度映象图数据中的位于被检测体P外部的无信号区域,通过实施区域增大处理的方式实施内插处理。
图8为说明通过26点法对三维敏感度映象图数据的位于被检测体P外部的无信号区域实施区域增大处理的场合下的方法的说明图,图9为说明通过6点方法对三维敏感度映象图数据的位于被检测体P外部的无信号区域实施区域增大处理的场合下的方法的说明图。
区域增大处理是一种用有信号区域处的值原封不动地对无信号区域处的值实施置换的处理方法。如果举例来说,采用26点方法的区域增大处理如图8所示,对于格子中点A的信号强度大于等于阈值(包含该阈值)的有信号区域,且与该中点A邻接的26个格点的各信号强度小于等于阈值(包含该阈值)的无信号区域的场合,使用中点A处的信号强度对无信号区域的26个格点处的信号强度实施置换处理。
采用6点方法的区域增大处理如图9所示,对于格子中点A的信号强度大于等于阈值(包含该阈值)的有信号区域,且与该中点A邻接的6个格点的各信号强度小于等于阈值(包含该阈值)的无信号区域的场合,使用中点A处的信号强度对无信号区域的6个格点处的信号强度实施置换处理。
图10为表示通过图1中所示的磁共振成象装置20,对三维敏感度映象图数据的位于被检测体P外部的无信号区域实施区域增大处理的一个实例的图。
实施线性内插处理后的三维敏感度映象图数据如图10(a)所示,具有位于被检测体P外部的无信号区域D1和位于被检测体P内部的有信号区域D2。对位于被检测体P外部的无信号区域D1实施区域增大处理的结果是,如图10(b)所示,将位于被检测体P外部的无信号区域D1置换为有信号区域D2而实施内插处理,使全部区域均为有信号区域D2。而且,对各个切割层面的所有区域均制作出三维敏感度映象图数据。
根据相控矩阵线圈24b的信道配置形式,相控矩阵线圈24b的敏感度分布在装置座标系统的Z轴方向(层面方向)上有可能不一样。
因此在步骤S28中,由层面方向加权处理部44h通过对三维敏感度映象图数据沿切割层面方向进行加权而对其实施修正。层面方向加权处理部44h可以通过诸如公式(1),对修正系数的倒数Y实施求解,进而通过将修正系数1/Y乘以三维敏感度映象图数据,来对其沿切割层面方向加权。
Y=1                            if Z<B  (1)
Y=A×(A-B)2+1                 if Z≥B
其中,Z:沿Z方向的层面位置
A:系数
B:偏移量
图11为表示使用图1中所示的磁共振成象装置20,对三维敏感度映象图数据沿切割层面方向实施加权处理时的修正系数1/Y的曲线图。
图11中的纵轴表示修正系数1/Y,横轴表示沿Z方向的层面位置Z。而且,图11中的实线表示当公式(1)中系数A=-20、偏移量B=0时的修正系数1/Y的曲线。
换句话说就是,对于沿Z方向的层面位置Z比预先确定的偏移量B=0更小的场合,设置为修正系数1/Y=1,而不对三维敏感度映象图数据实施修正。
在另一方面,对于沿Z方向的层面位置Z大于等于预先确定的偏移量B=0(包括0)的场合,利用二次函数的计算公式时修正系数1/Y实施计算,进而通过将所获得修正系数1/Y乘以三维敏感度映象图数据而对其实施修正。
然而,对修正系数的倒数Y实施计算的计算公式也可以不采用二次函数,而是采用其它任意函数实施近似处理。
随后在步骤S29中,由三维平滑处理部44i对三维敏感度映象图数据实施适当的、具有所需要的强度的3D平滑滤波处理。采用这种方式,可以去除局部位置处的值极端大的部位,从而可以沿RO方向、PE方向、SL方向的各方向上提高连续性,高精度地生成出最终的三维敏感度映象图数据。
随后在步骤S30中,将通过敏感度分布推定组件44生成出的三维敏感度映象图数据,写入并储存在敏感度映象图用数据库45处。
随后在图4的步骤S12中,由主扫描运行组件37将图象获得用序列信号传送至序列信号控制器控制组件38处,执行主扫描。而且,通过由图象再构成组件40对所收集到的原始数据实施图象再构成处理的方式,获取图象数据。
随后在步骤S13中,使用三维敏感度映象图数据对通过主扫描获得的图象数据的亮度实施修正。采用这种方式,可以按照与主扫描中的诸如摄象剖面方向、空间分辨率等摄象条件、数据收集条件、图象再构成条件等的各项条件相对应地,由图象数据修正组件46从敏感度映象图用数据库45中获取出对应的三维敏感度映象图数据。
而且,图象数据修正组件46使用获取出的三维敏感度映象图数据,对图象数据的亮度实施修正。在这时还可以根据需要,对所获取出的三维敏感度映象图数据实施标准化处理。
采用这种构成形式,可以抑制由于相控矩阵线圈24b的敏感度离散产生的信号强度不均匀性的影响,从而可以获得改善亮度后的图象数据。
如果采用具有上述构成形式的磁共振成象装置20,可以依据实施敏感度预扫描所获得的图象数据,更高精度地对相控矩阵线圈24b的敏感度分布实施推定,进而可以依据所获得的相控矩阵线圈24b的敏感度分布,对通过实施主扫描所获得的图象数据的亮度实施更加良好的修正。
对于使用属于在先技术的磁共振成象装置1的图象数据亮度修正方法,在无法对图象数据的亮度实施足够修正的各种场合,比如说在敏感度预扫描摄象区域中包含有诸如肺等的无信号区域的场合,由于与被检测体P的思维无关地内脏形状变化,而在实施敏感度预扫描摄象时通过相控矩阵线圈24b实施摄象获得的被检测体P的形状与通过WB线圈24a实施摄象获得的被检测体P的形状间产生有偏差的场合,以及由于相控矩阵线圈24b的各信道配置特性,会沿装置座标系的Z轴方向产生信号强度不均匀性的场合,如果采用磁共振成象装置20则可以对图象数据的亮度实施良好的修正。
图12为表示使用图1中所示的磁共振成象装置20获得的亮度修正后的被检测体P的一个断层图象的示意图,图13为表示使用在先技术中的磁共振成象装置20获得的亮度修正后的被检测体P的一个断层图象的示意图。
正如图13所示,在实施敏感度预扫描摄象时通过相控矩阵线圈24b实施摄象获得的被检测体P的形状与通过WB线圈24a实施摄象获得的被检测体P的形状间产生有偏差时,无法对图象数据的亮度实施良好的修正。
在另一方面,正如图12所示,通过对WB再构成图象和PAC再构成图象实施区域抑制处理的方式,即使在通过相控矩阵线圈24b实施摄象获得的被检测体P的形状与通过WB线圈24a实施摄象获得的被检测体P的形状间产生有偏差时,也可以高精度地生成提高连续性后的三维敏感度映象图数据,从而可以对图象数据的亮度实施良好的修正。
而且,在具有如上所述构成形式的磁共振成象装置20中,也可以省略数据处理过程中的一部分,并且省略与该部分相关的构成要素中的一部分。而且,还可以代替相控矩阵线圈24b,而是采用单一线圈的构成形式。
图14为表示本发明的磁共振成象装置的第二实施形式的功能方框图。
如图14所示的磁共振成象装置20A,其RF线圈24和信号接收器30的具体构成形式和计算机32所具有的功能,与如图1所示的磁共振成象装置20不同。其它构成要素和作用均与如图1所示的磁共振成象装置20实质上相同,因此采用相同的参考标号表示相同的构成部分,并且省略了对这些部分的详细说明。
图15为表示图14中所示的RF线圈24和信号接收器30的一个实例的详细构成图。
RF线圈24由对高频信号实施发送用的WB线圈24a、作为对主扫描时的NMR信号实施接收的线圈用的、成为主线圈的相控矩阵线圈24b构成。相控矩阵线圈24b具有多个表面线圈24c。
信号接收器30由多个信号接收系统回路30a构成。而且,各表面线圈24c分别与信号接收器30的信号接收系统回路30a相连接,WB线圈24a与信号发送器29相连接。另外,WB线圈24a还可以与信号接收器30的信号接收系统回路30a相连接。
WB线圈24a和相控矩阵线圈24b的具体配置形式与如图3所示的构成形式相同。
换句话说就是,相控矩阵线圈24b处的各表面线圈24c,例如围绕Z轴对称地配置在包含被检测体P的特定关注部位的剖面L周围。而且,将WB线圈24a设置在相控矩阵线圈24b的外侧。RF线圈24构成为通过WB线圈24a向被检测体P发送出高频信号,并且利用相控矩阵线圈24b的各表面线圈24c,通过多信道对从包含特定关注部位的剖面L处给出的NMR信号实施信号接收,进而传送至各信号接收器30。
另外,计算机32通过对程序实施读取和运行的方式,作为敏感度预扫描运行组件36、敏感度预扫描条件设定组件50、主扫描运行组件37、序列信号控制器控制组件38、原始数据用数据库39、图象再构成组件40、图象数据用数据库41、主线圈图象用数据库51、敏感度分布推定组件44、敏感度映象图用数据库45、图象数据修正组件46、显示处理组件47而发挥功能。当然,也可以不采用程序方式,而是通过特定的回路设计构成这种计算机32。
敏感度预扫描运行组件36具有:依据从敏感度预扫描条件设定组件50接收到的摄象条件,对为了求出作为相控矩阵线圈24b的敏感度分布的三维敏感度映象图数据而实施敏感度预扫描时的序列信号(敏感度推定用序列信号)实施生成的功能;将所生成出的敏感度推定用序列信号传送至序列信号控制器控制组件38以实施敏感度预扫描的功能。
敏感度预扫描条件设定组件50具有:按照为了求解三维敏感度映象图数据而使对比度充分低的方式,对敏感度预扫描时的摄象条件实施设定的功能;将所设定的敏感度预扫描摄象条件传送至敏感度预扫描运行组件36处的功能。
主扫描运行组件37具有将执行用于获取图象数据的主扫描时使用的各种序列信号,传送至序列信号控制器控制组件38以实施主扫描的功能。
序列信号控制器控制组件38具有依据由输入装置33或其它构成要素给出的信息,将从敏感度预扫描运行组件36和主扫描运行组件37接收到的序列信号中所需要的序列信号,传送至序列信号控制器31处以实施敏感度预扫描或主扫描的功能。而且,序列信号控制器控制组件38还具有接收由序列信号控制器31通过实施敏感度预扫描或主扫描所收集到的相控矩阵线圈24b的各表面线圈24c的原始数据,并将其配置在形成在原始数据用数据库39处的K空间(傅立叶空间)处的功能。
因此,原始数据用数据库39对在信号接收器30中生成出的每个表面线圈24c的各原始数据实施储存。换句话说就是,可以将原始数据配置在形成在原始数据用数据库39中的K空间处。
图象再构成组件40具有:对通过运行主扫描而配置在原始数据用数据库39的K空间处的原始数据,实施诸如傅立叶变换(FT)等图象再构成处理以再构成出被检测体P的图象数据的功能;将再构成出的图象数据写入至图象数据用数据库41处的功能。
而且,图象再构成组件40还具有:对通过运行敏感度预扫描而配置在原始数据用数据库39的K空间处的原始数据,按照对运行主扫描而获得的原始数据实施再构成处理相同的方式实施再构成处理,将被检测体P的图象数据作为主线圈图象数据实施再构成处理的功能;将再构成出的主线圈图象数据写入至主线圈图象用数据库51的功能。
敏感度分布推定组件44具有:通过将储存在主线圈图象用数据库51处的主线圈图象数据,作为构成敏感度修正用原始数据的敏感度推定用数据使用,从而制作出作为敏感度修正用数据的相控矩阵线圈24b的三维敏感度映象图数据的功能;将所制作出的三维敏感度映象图数据写入至敏感度映象图用数据库45处的功能。敏感度分布推定组件44具有阈值处理部44j、区域缩小部44k、内插处理部441、平滑处理部44m。
阈值处理部44j具有对主线圈图象数据实施阈值处理的功能,即具有对主线圈图象数据的信号强度小于等于预先设定的阈值(包含该阈值)的部分处的数据实施屏蔽的功能。
区域缩小部44k具有对实施敏感度分布推定时使用的主线圈图象数据的区域实施缩小处理,将屏蔽区域附近处的信号强度比较小的部分从三维敏感度映象图数据制作时使用的数据中除去的功能。
内插处理部441具有通过外插处理或内插处理,对作为主线圈图象数据实施区域缩小处理后被屏蔽了的无信号区域中的三维敏感度映象图数据实施推定的方式,对敏感度分布推定用的主线圈图象数据实施内插的功能。
平滑处理部44m具有对敏感度分布推定用的主线圈图象数据,实施平滑处理以制作出最终的三维敏感度映象图数据的功能。
图象数据修正组件46具有:从储存在敏感度映象图用数据库45处的三维敏感度映象图数据中,分离抽取与实施主扫描时的摄象条件、数据收集条件、图象再构成条件等的图象数据获取条件相对应的三维敏感度映象图数据的功能;使用所抽取出的三维敏感度映象图数据,对通过运行主扫描而获得的储存在图象数据用数据库41处的图象数据的信号强度实施修正的功能。
显示处理组件47具有将储存在图象数据用数据库41处的图象数据,传送至显示装置34处以实施显示的功能。
具有上述构成形式的磁共振成象装置20A,作为整体,各个构成要素具有作为将图象摄象时的信号接收线圈作为信号接收线圈,而执行用于生成实施主扫描图象摄象时的信号接收线圈的敏感度修正用数据的敏感度预扫描的扫描运行组件、仅仅将实施敏感度预扫描的图象摄象中的由信号接收线圈获得的数据作为原始数据,而生成敏感度修正数据的敏感度修正用数据生成组件而使用的功能。
图16为表示使用图14中所示的磁共振成象装置20A,对被检测体P实施断层图象摄象时的步骤的流程图,图中在字母S后附加有数字的参考标号,表示流程中的各步骤。
首先在步骤S40中,运行敏感度预扫描作业。敏感度预扫描条件设定组件50对敏感度预扫描的摄象条件进行设置,并将敏感度推定用序列信号传送至敏感度预扫描运行组件36处。在此,将敏感度预扫描的摄象条件设置为,使为了制作三维敏感度映象图数据而执行再构成处理获得的图象具有充分低的对比度。
如果举例来说,低对比度的摄象条件,可以设定为通过延长重复时间(TR:repetition time)且缩短回波时间(TE:echo time)的方式,使对T1(纵缓和时间)和T2(横缓和时间)中的任何一个的影响均比较小的质子密度强调图象摄象时的条件,或是与该条件相近似的摄象条件。
而且如果举例来说,实施敏感度预扫描用的敏感度推定用序列信号,可以是在高速傅立叶回波(FFE)型序列信号时将TE设定为比较短的1-5毫秒(ms)左右,将反转(flip)角设定为比较小的5~10度左右时的信号。而且,当TR为200毫秒(ms)左右时,可以对大于等于20幅(包括20幅)的切割层面的主线圈图象数据实施收集,进而可以获得反映了制作三维敏感度映象图数据所需要的表面线圈24c的敏感度的整个体积的主线圈图象数据。
在另一方面,敏感度预扫描的摄象可以是对2D的多重切割层面进行的摄象,也可以为3D摄象。实施3D摄象时的敏感度推定用序列信号,例如可以采用在FFE型序列信号中将TE设定为比较短的1-5毫秒(ms)左右,将反转(flip)角设定为比较小的小于等于5度(包含5度)的信号。而且,当TR为10毫秒(ms)左右时,可以在与2D摄象同样的时间内获得所需要的主线圈图象数据。
通过敏感度预扫描运行组件36将敏感度推定用序列信号传送至序列信号控制器控制组件38处,并仅仅将实施主扫描时的信号接收用线圈(主线圈)的相控矩阵线圈24b作为信号接收用线圈实施敏感度预扫描。换句话说就是,在进行获取图象数据用的主扫描之前,先进行用于获取相控矩阵线圈24b的敏感度映象图数据的敏感度预扫描作业。
将被检测体P预先设置在躺卧平台35上,由静磁场用电源26向静磁场用磁铁21供给电流以在摄象区域形成静磁场。而且,由均场线圈用电源28向均场线圈22供给电流,以使形成在摄象区域处的静磁场均匀化。
随后,由输入装置33向序列信号控制器控制组件38发出动作指令。在接收到该指令时,序列信号控制器控制组件38将敏感度推定用序列信号传送至序列信号控制器31。序列信号控制器31依据该敏感度推定用序列信号对倾斜磁场用电源27、信号发送器29和信号接收器30实施驱动,从而在设置有被检测体P的摄象区域处形成X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz,并且同时产生出高频信号。
在这时,由倾斜磁场线圈产生的X轴倾斜磁场Gx、Y轴倾斜磁场Gy、Z轴倾斜磁场Gz,主要是分别作为相位编码(PE)用倾斜磁场、读取(RO)用倾斜磁场、层面编码(SE)用倾斜磁场使用。因此,被检测体P的内部的原子核自旋沿转动方向表现出规则性,作为利用SE用倾斜磁场沿Z轴方向形成的层面中的二维位置信息的X坐标和Y坐标,通过PE用倾斜磁场和RO用倾斜磁场分别被变换为被检测体P内部的原子核自旋的相位变化量和频率变化量。
从信号发送器29与敏感度推定用序列信号对应地将高频信号施加至RF线圈24的WB线圈24a处,并且从WB线圈24a将高频信号传送至被检测体P处。而且,在被检测体P的内部与高频信号的频率相对应的包含在切割层面中的原子核磁共振产生的NMR信号,通过RF线圈24的、作为主线圈的相控矩阵线圈24b的各表面线圈24c,通过多信道被接收并分别传送至各信号接收器30处。
而且,WB线圈24a可以仅仅用于对高频信号实施发送,而不用于信号的接收作业。
各信号接收器30从相控矩阵线圈24b的各表面线圈24c处接收NMR信号,并且进行诸如前置放大、中间频率变换、相位检波、低频放大、滤波等的各种信号处理。各信号接收器30还通过对NMR信号实施A/D变换的方式,生成出作为数字数据的NMR信号的原始数据。信号接收器30还将所生成出的原始数据传送至序列信号控制器31处。
序列信号控制器31将从信号接收器30接收到的原始数据传送至序列信号控制器控制组件38处,由序列信号控制器控制组件38将原始数据配置在形成在原始数据用数据库39中的K空间处。图象再构成组件40通过对配置在形成在原始数据用数据库39中的K空间处的原始数据实施傅立叶变换(FT),来对作为主线圈图象数据的、构成被检测体P的三维图象数据的体积数据实施再构成处理,并且将获得相控矩阵线圈24b的三维敏感度映象图数据用的敏感度推定用序列信号,写入至主线圈图象用数据库51处。
在这时,对通过敏感度预扫描收集到的原始数据实施的再构成处理方法,是与主扫描中的再构成处理方法相同的方法。如果举例来说,在主扫描中使用相控矩阵线圈24b进行的再构成处理方法,可以为对通过各表面线圈24c获得的图象数据求解平方和的平方根的Sum of Square处理(SoS处理)方法,和对通过各表面线圈24c获得的图象数据求解信号强度绝对值的和的再构成处理方法等。
随后在步骤S41中,敏感度分布推定组件44通过将储存在主线圈图象用数据库51处的成为体积数据的主线圈图象数据作为敏感度推定用数据使用的,对敏感度分布实施推定。
图17为表示在图16中所示的示意性流程图中,对相控矩阵线圈24b的敏感度分布实施推定时的步骤的流程图,图中在字母S后附加有数字的参考标号,表示流程中的各步骤。
首先在步骤S50中,将主线圈图象数据从主线圈图象用数据库51读取至敏感度分布推定组件44处。
图18为表示在低对比度摄象条件下获得的主线圈图象数据的信号强度分布的示意图,图19为表示在常规对比度摄象条件下获得的主线圈图象数据的信号强度分布的示意图。
在图18和图19中,纵轴表示主线圈图象数据的信号值,横轴表示包含如图3所示的关注区域(ROI:region of interest)的剖面L方向的位置。
正如图19所示,对于在具有常规对比度的摄象条件下获得的主线圈图象数据,由于对比度不同产生的影响并不是充分小,所以难以原封不动地作为相控矩阵线圈24b的敏感度推定用数据使用。
在另一方面,正如图18所示,对于在低对比度摄象条件下获得的主线圈图象数据,由于对比度不同产生的影响充分小,所以即使原封不动地作为相控矩阵线圈24b的敏感度推定用数据使用,所产生的误差也相当小。
而且在此,将主线圈图象数据作为包含ROI的直线L上的一维数据进行说明的,然而实际上也可以构成为采用实施二维或三维摄象,将二维或三维主线圈图象数据作为敏感度推定用数据使用,制作三维敏感度映象图数据用的各种处理对象。
因此,将在低对比度的摄象条件下获得的主线圈图象数据作为敏感度推定用数据使用,实施用于制作三维敏感度映象图数据所需要的各种处理作业。下面,对在先技术中的敏感度映象图数据的制作方法进行比较说明。
图20为表示使用在先技术中的磁共振成象装置1生成三维敏感度映象图数据时所生成出的数据的示意图。
属于在先技术的磁共振成象装置1,为了将相控矩阵线圈24b和WB线圈24a均作为信号接收用线圈实施敏感度预扫描,而将如图20(a)所示的通过相控矩阵线圈24b获得的主线圈图象数据、如图20(b)所示的通过WB线圈24a获得的WB线圈图象数据这两个数据均作为敏感度推定用数据使用。而且,通过将主线圈图象数据和WB线圈图象数据双方与预先设定的阈值ε、ε’进行比较来实施阈值处理,将无信号区域部分由敏感度推定用数据中去除。
通过将实施阈值处理后的敏感度推定用数据区域D1中的主线圈图象数据除以WB线圈图象数据,来制作出如图20(c)所示的无维化后的三维敏感度映象图数据。而且,通过对整个区域实施外插处理或内插处理等的内插处理,来对三维敏感度映象图数据实施推定,进而制作出如图20(d)所示的三维敏感度映象图数据。
在另一方面,如图14所示的磁共振成象装置20A,可以仅采用相控矩阵线圈24b作为信号接收用线圈实施敏感度预扫描。
图21为说明使用图14中所示的磁共振成象装置20A实施三维敏感度映象图数据生成时所存在的问题的说明图,图22为表示通过能够避免如图21所示的问题的步骤,使用磁共振成象装置20A实施三维敏感度映象图数据生成时所生成出的数据的示意图。
换句话说就是,在如图17所示的步骤S51中,通过阈值处理部44j对主线圈图象数据实施阈值处理。即,如图21(a)所示和如图22(a)所示,对主线圈图象数据中信号强度小于等于预先设定的阈值ε(直至与该阈值ε相等)的部分的数据实施屏蔽,以将位于被检测体P的区域之外和肺部等的无信号区域的部分,从敏感度推定用数据中去除。
对于诸如质子密度强调图象等的低对比度的主线圈图象数据,在实施阈值处理后的敏感度推定用数据区域D1中可以不对WB线圈图象数据实施除法运算,而将其原封不动地作为敏感度推定用数据使用,但如果仅仅实施阈值处理,如图21(b)所示不消除位于被检测体P的区域之外和肺部等给出的无信号区域附近位置处的信号强度比较低部分的影响,则敏感度推定用数据的值比较小,将难以对三维敏感度映象图数据实施正确制作。
因此在步骤S52中,由区域缩小部44k对作为敏感度推定用数据使用的区域实施缩小处理。换句话说就是,如图22(b)所示在敏感度推定用数据区域D1的屏蔽区域的边界部附近位置处,通常会出现信号强度比其它部分小的现象,所以通过对敏感度推定用数据区域D1实施缩小处理,来去除信号强度小的边缘部分的区域。
随后在步骤S53中,如图22(c)所示,将实施区域缩小处理后的新的敏感度推定用数据区域D2中的主线圈图象数据的信号强度,看作是三维敏感度映象图数据。而且,通过内插处理部441实施外插处理或内插处理等的内插处理,来对实施区域缩小处理后的被屏蔽了的无信号区域中的三维敏感度映象图数据实施推定,制作出如图22(d)所示的整个区域中的三维敏感度映象图数据。
随后在步骤S54中,通过对内插处理后的整个区域中的三维敏感度映象图数据,进行诸如正交函数展开等的数据拟合处理,来实施平滑处理。采用这种方式,最终制作出连续性好的三维敏感度映象图数据。
而且,在图16所示的步骤S42中,将相控矩阵线圈24b的三维敏感度映象图数据储存在敏感度映象图用数据库45处。
随后在步骤S43中,将图象获取用序列信号由主扫描运行组件37传送至序列信号控制器控制组件38处,将相控矩阵线圈24b作为信号接收用线圈实施主扫描作业。然后,收集原始数据并由图象再构成组件40实施图象再构成处理,获得图象数据。
随后在步骤S44中,按照与主扫描中的诸如摄象剖面方向、空间分辨率等摄象条件、数据收集条件、图象再构成条件等的各项条件相对应的方式,由图象数据修正组件46从敏感度映象图用数据库45中取得相对应的三维敏感度映象图数据。
随后在步骤S45中,由图象数据修正组件46使用取得的三维敏感度映象图数据,对图象数据实施修正。换句话说就是,通过将三维敏感度映象图数据的倒数乘以图象数据的各信号强度的方式,实施修正处理。在这时,还可以适当地进行使三维敏感度映象图数据非零化的处理,或是对于三维敏感度映象图数据为零的场合,进行诸如场合分类处理等的图象数据修正的各种常规的错误处理。
采用这种构成形式,可以对由于相控矩阵线圈24b的敏感度离散产生的信号强度不均匀性的影响实施抑制,可以获得对图象质量实施改善后的图象数据。
如果采用具有上述构成形式的磁共振成象装置20A,则可以采用由多个表面线圈构成的相控矩阵线圈24b进行图象摄象,即使信号接收用线圈的敏感度产生离散,由于并不将WB线圈24a作为实施敏感度预扫描时的信号接收用线圈,所以可以在更短的时间里按照与图象种类等摄象条件无关的方式,对图象数据的信号强度起伏实施修正。而且,与在先技术相比,可以容易地在具有良好精度的条件下获得在整个摄象视野范围内均具有比较高诊断能力的图象。
如果举例来说,对于在如上所述的TR为200毫秒(ms)的摄象条件下,对48×48的矩阵实施摄象时,在先技术中为了能够在实施敏感度预扫描时获得WB线圈24a和相控矩阵线圈24b两者的数据,摄象时间为19.2秒,然而如果采用磁共振成象装置20A,通过在一半的时间、即9.6秒中进行憋气摄象,就可以获得制作三维敏感度映象图数据所需要的全部必需的主线圈图象数据,因此可以将摄象时间减少一半,从而可以降低患者的负担。
而且,采用本方法时由于不需要进行使用WB线圈24a的摄象,所以可以在WB线圈24a和相控矩阵线圈24b的去耦性并不十分充分的状态下运行,并且可以防止出现WB线圈24a与相控矩阵线圈24b间的数据位置偏差等的错误。
图23为表示本发明的磁共振成象装置的第三实施形式的功能方框图。
如图23所示的磁共振成象装置20B,在计算机32作为均场(shimming)用摄影条件设定组件60使用的功能方面,与如图14所示的磁共振成象装置20A不同。其它构成要素和作用均与如图14所示的磁共振减象装置20A实质上相同,因此采用相同的参考标号表示相同的构成部分,并且省略了对这些部分的详细说明。
磁共振成象装置20B的计算机32,还具有能够作为均场用摄影条件设定组件60使用的功能。均场用摄影条件设定组件60具有对为了对静磁场空间不均匀性实施修正而进行的均场处理的摄象条件实施设定,并且将其传送至敏感度预扫描运行组件36处的功能。采用这种构成形式,敏感度预扫描运行组件36构成为生成能够同时实施敏感度预扫描和均场作业的摄象条件下的序列信号,并且将其传送至序列信号控制器控制组件38处。
磁共振成象装置20B可以在实施敏感度预扫描作业的同时实施均场作业。这时的摄象序列信号,例如在FFE中可以将TE设定为4.5毫秒(ms)/9.0毫秒(ms)的两个回波,将反转角设定为5~10度左右。而且,可以根据两个回波间的信号相位差求解磁场分布并进行均场作业,将4.6毫秒(ms)的数据作为敏感度推定用数据。
因此,如果采用磁共振成象装置20B,除了可以获得磁共振成象装置20A所具有的技术效果之外,还可以通过更高的效率实施摄象作业。
而且在磁共振成象装置20A、20B中,不只是相控矩阵线圈24b,还可以将头部用线圈、各种阵列线圈、表面线圈等目的不同的各种线圈,作为实施主扫描时的信号接收用线圈,即作为RF线圈24的主线圈。而且,RF线圈24或主线圈也可以采用单一的线圈构成。
如果采用这种构成形式,也可以使用WB线圈24a自身制作出WB线圈24a的敏感度映象图数据。WB线圈24a的敏感度起伏比相控矩阵线圈24b小,但即使今后进一步使装置小型化,基本上也无法对WB线圈24a产生的敏感度起伏忽略不计。所以,如果仅将WB线圈24a作为实施敏感度预扫描时的信号接收用线圈,制作WB线圈24a的敏感度映象图数据,则可以容易地使装置进一步小型化。
在另一方面,如果仅将实施主扫描时由信号接收用线圈获得的图象数据作为原始数据,来生成敏感度映象图数据,则也可以采用实施主扫描时的信号接收用线圈之外的其它线圈,作为实施敏感度预扫描时的信号接收用线圈。
而且,还可以对上述实施形式的磁共振成象装置20、20A、20B实施相互组合,或是省略其中的一部分构成要素和功能。

Claims (14)

1.一种磁共振成象装置,具有:
扫描实施组件,实施用于生成RF线圈的敏感度映象图数据的扫描作业;
区域抑制处理组件,对由前述扫描作业获得的图象数据的位于无信号区域附近位置处的有信号区域实施区域抑制处理;
使用实施区域抑制处理后的前述图象数据,生成敏感度映象图数据的敏感度映象图数据生成组件;
以及对前述敏感度映象图数据实施三维平滑滤波处理的平滑处理组件。
2.一种磁共振成象装置,具有:
扫描实施组件,实施用于生成RF线圈的敏感度映象图数据的扫描作业;
使用通过前述扫描作业获得的图象数据,生成敏感度映象图数据的敏感度映象图数据生成组件;
对前述敏感度映象图数据的位于被检测体内的无信号区域实施线性内插处理的线性内插处理组件;
以及对前述敏感度映象图数据实施三维平滑滤波处理的平滑处理组件。
3.一种磁共振成象装置,具有:
将实施图象摄象时的信号接收用线圈作为信号接收用线圈,而执行用于生成实施前述图象摄象时的信号接收用线圈的敏感度修正用数据的扫描作业的扫描实施组件;
以及仅仅将在前述扫描作业中由实施前述图象摄象时的信号接收用线圈获得的数据作为原始数据,而生成前述敏感度修正用数据的敏感度修正用数据生成组件。
4.如权利要求1所述的磁共振成象装置,其特征在于:进一步设置有对前述敏感度映象图数据的位于被检测体外的无信号区域实施区域增大处理的区域增大处理组件。
5.如权利要求1所述的磁共振成象装置,其特征在于:进一步设置有通过沿切割层面方向进行加权,来对前述敏感度映象图数据实施修正的层面方向加权处理组件。
6.如权利要求2所述的磁共振成象装置,其特征在于:进一步设置有对前述敏感度映象图数据的位于被检测体外的无信号区域实施区域增大处理的区域增大处理组件。
7.如权利要求2所述的磁共振成象装置,其特征在于:进一步设置有通过沿切割层面方向进行加权的方式,对前述敏感度映象图数据实施修正的层面方向加权处理组件。
8.如权利要求2所述的磁共振成象装置,其特征在于:进一步设置有实施使用任意变换函数将前述敏感度映象图数据变换为适合于线性内插的平坦化分布的数据平坦化处理,并且通过使用逆变换函数从而使实施线性内插处理后的前述敏感度映象图数据恢复为前述数据平坦化处理之前的分布的数据平坦化处理组件,其中
前述线性内插处理组件对前述数据平坦化处理后的前述敏感度映象图数据的无信号区域实施线性内插处理,前述平滑处理组件对恢复为通过前述数据平坦化处理组件实施前述数据平坦化处理之前的分布的前述敏感度映象图数据,实施三维平滑滤波处理。
9.如权利要求3所述的磁共振成象装置,其特征在于:将用于生成前述敏感度修正用数据的扫描时的摄象条件,设置为获得为了能将仅仅使用前述图象摄象时的信号接收用线圈获得的数据作为前述敏感度修正用数据的原始数据而充分低的对比度的图象的条件。
10.如权利要求3所述的磁共振成象装置,其特征在于:使用对前述敏感度修正用数据的原始数据实施阈值处理所抽取出的数据的区域缩小处理后的数据,生成前述敏感度修正用数据。
11.如权利要求3所述的磁共振成象装置,其特征在于:前述扫描实施组件将表面线圈作为信号接收用线圈,而执行用于生成前述敏感度修正用数据的扫描。
12.一种磁共振成象装置的数据处理方法,具有:
对实施用于生成RF线圈的敏感度映象图数据的扫描而获得的图象数据的位于无信号区域附近位置处的有信号区域,实施区域抑制处理的步骤;
使用实施区域抑制处理后的前述图象数据,生成敏感度映象图数据的步骤;
以及对前述敏感度映象图数据实施三维平滑滤波处理的步骤。
13.一种磁共振成象装置的数据处理方法,具有:
使用实施用于生成RF线圈的敏感度映象图数据的扫描而获得的图象数据,生成敏感度映象图数据的步骤;
对前述敏感度映象图数据的位于被检测体内的无信号区域实施线性内插处理的步骤;
以及对前述敏感度映象图数据实施三维平滑滤波处理的步骤。
14.一种磁共振成象装置的数据处理方法,具有:
将实施图象摄象时的信号接收用线圈作为信号接收用线圈,执行用于生成实施前述图象摄象时的信号接收用线圈的敏感度修正用数据的扫描作业的步骤;
以及仅仅将前述扫描作业中由实施前述图象摄象时的信号接收用线圈获得的数据作为原始数据,生成前述敏感度修正用数据的步骤。
CNB2005100717180A 2004-02-26 2005-02-25 磁共振成象装置及其磁共振成象装置的数据处理方法 Expired - Fee Related CN100444791C (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2004052312 2004-02-26
JP2004052313 2004-02-26
JP2004052312A JP4723814B2 (ja) 2004-02-26 2004-02-26 磁気共鳴イメージング装置
JP2004052313A JP4664605B2 (ja) 2004-02-26 2004-02-26 磁気共鳴イメージング装置

Related Child Applications (2)

Application Number Title Priority Date Filing Date
CN 200810098596 Division CN101305909B (zh) 2004-02-26 2005-02-25 磁共振成象装置及其磁共振成象装置的数据处理方法
CN201110310625.4A Division CN102385044B (zh) 2004-02-26 2005-02-25 磁共振成象装置及磁共振成象装置的数据处理方法

Publications (2)

Publication Number Publication Date
CN1676097A true CN1676097A (zh) 2005-10-05
CN100444791C CN100444791C (zh) 2008-12-24

Family

ID=34914445

Family Applications (2)

Application Number Title Priority Date Filing Date
CNB2005100717180A Expired - Fee Related CN100444791C (zh) 2004-02-26 2005-02-25 磁共振成象装置及其磁共振成象装置的数据处理方法
CN201110310625.4A Expired - Fee Related CN102385044B (zh) 2004-02-26 2005-02-25 磁共振成象装置及磁共振成象装置的数据处理方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201110310625.4A Expired - Fee Related CN102385044B (zh) 2004-02-26 2005-02-25 磁共振成象装置及磁共振成象装置的数据处理方法

Country Status (4)

Country Link
US (1) US7385397B2 (zh)
EP (3) EP1723906B1 (zh)
CN (2) CN100444791C (zh)
WO (1) WO2005082245A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102670201A (zh) * 2011-03-11 2012-09-19 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN103222868A (zh) * 2012-01-25 2013-07-31 株式会社东芝 磁共振成像装置、以及磁共振成像方法
CN103874457A (zh) * 2011-09-29 2014-06-18 株式会社日立医疗器械 磁共振摄像装置、高频磁场照射方法以及程序

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7446529B2 (en) * 2005-07-21 2008-11-04 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and radio frequency coil unit
JP4785566B2 (ja) * 2006-03-08 2011-10-05 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置と磁気共鳴イメージング方法
CN101484821A (zh) * 2006-04-13 2009-07-15 皇家飞利浦电子股份有限公司 基于波传播的线圈灵敏度估计
US7622920B2 (en) * 2006-07-06 2009-11-24 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus capable of automatically determining RF coil positions
DE102007037657B4 (de) * 2007-08-09 2009-05-14 Siemens Ag Verfahren zur k-Raum-Aufteilung bei der Aufzeichnung von Messdaten mittels Magnet- Resonanz- Technik und Magnet- Resonanz- Gerät
JP5570740B2 (ja) * 2008-03-03 2014-08-13 株式会社東芝 磁気共鳴イメージング装置および撮像条件生成方法
JP5443695B2 (ja) * 2008-03-05 2014-03-19 株式会社東芝 磁気共鳴イメージング装置
US9036884B2 (en) 2009-09-17 2015-05-19 Koninklijke Philips N.V. Image intensity correction for magnetic resonance imaging
CN103380384B (zh) * 2011-02-15 2016-04-06 皇家飞利浦有限公司 使用交错的水参考扫描的具有自动相位和b0校正的磁共振波谱成像方法和装置
JP5931406B2 (ja) 2011-11-09 2016-06-08 株式会社東芝 磁気共鳴イメージング装置
JP5904815B2 (ja) * 2012-02-15 2016-04-20 株式会社東芝 磁気共鳴イメージング装置
KR101310825B1 (ko) * 2012-05-10 2013-10-14 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
JP5897415B2 (ja) 2012-06-29 2016-03-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
WO2014042678A1 (en) * 2012-09-13 2014-03-20 The Regents Of The University Of California System and method for automated detection of lung nodules in medical images
US10191126B2 (en) * 2013-04-04 2019-01-29 Regents Of The University Of Minnesota Systems and methods for spatial gradient-based electrical property properties tomography using magnetic resonance imaging
DE102014226034B4 (de) * 2014-12-16 2017-01-19 Siemens Healthcare Gmbh Bildkorrektur bei der MR-Bildgebung unter Berücksichtigung des Empfangsprofils
CN108013876B (zh) * 2016-10-31 2023-05-23 通用电气公司 磁共振图像均匀度校正方法及装置
US11579219B2 (en) 2017-10-06 2023-02-14 University Of Cincinnati Systems and methods for estimating complex radiofrequency fields in a magnetic resonance imaging
EP3865890B1 (de) * 2020-02-13 2023-05-24 Siemens Healthcare GmbH Magnetresonanztomograph mit b0-modulation und verfahren zum betrieb

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL8602821A (nl) 1986-11-07 1988-06-01 Philips Nv Werkwijze en inrichting voor het bepalen van gecorrigeerd mri oppervlaktespoelbeeld.
US5208534A (en) * 1989-08-09 1993-05-04 Kabushiki Kaisha Toshiba Magnetic resonance imaging system
JP2808782B2 (ja) * 1990-01-31 1998-10-08 株式会社島津製作所 Mr装置
JP3135592B2 (ja) 1990-02-28 2001-02-19 株式会社東芝 磁気共鳴映像装置
JP3498980B2 (ja) * 1992-09-16 2004-02-23 東芝医用システムエンジニアリング株式会社 磁気共鳴イメージング用スキャン位置決め方法及び磁気共鳴イメージング装置
JPH0759750A (ja) 1993-08-30 1995-03-07 Hitachi Medical Corp 核磁気共鳴イメージング装置
JP3600351B2 (ja) * 1996-03-12 2004-12-15 株式会社東芝 磁気共鳴イメージング装置
US5943433A (en) * 1996-12-30 1999-08-24 General Electric Company Method for correcting inhomogeneity of spatial intensity in an aquired MR image
JP3865887B2 (ja) * 1997-08-27 2007-01-10 株式会社日立メディコ 画像補正方法
JP4008193B2 (ja) * 2000-08-29 2007-11-14 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像処理方法、画像処理装置、mri装置および記録媒体
CN1491095A (zh) * 2001-01-19 2004-04-21 株式会社东芝 使用高精度线圈感度图的并行mr成象
JP3896007B2 (ja) 2002-03-07 2007-03-22 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴撮像装置
JP4030805B2 (ja) * 2002-06-07 2008-01-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー コイル感度マップ作成方法およびmri装置
US6836116B2 (en) * 2002-07-01 2004-12-28 Ge Medical Systems Global Technology Company, Llc. Sensitivity encoding magnetic resonance imaging method and system
US6995559B2 (en) * 2003-10-30 2006-02-07 Ge Medical Systems Global Technology Company, Llc Method and system for optimized pre-saturation in MR with corrected transmitter frequency of pre-pulses

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102670201A (zh) * 2011-03-11 2012-09-19 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN102670201B (zh) * 2011-03-11 2016-08-03 西门子公司 用于产生磁共振图像数据的方法和磁共振断层造影系统
CN103874457A (zh) * 2011-09-29 2014-06-18 株式会社日立医疗器械 磁共振摄像装置、高频磁场照射方法以及程序
CN103874457B (zh) * 2011-09-29 2016-02-03 株式会社日立医疗器械 磁共振摄像装置以及高频磁场照射方法
CN103222868A (zh) * 2012-01-25 2013-07-31 株式会社东芝 磁共振成像装置、以及磁共振成像方法
CN103222868B (zh) * 2012-01-25 2016-08-17 东芝医疗系统株式会社 磁共振成像装置、以及磁共振成像方法
US10393844B2 (en) 2012-01-25 2019-08-27 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method

Also Published As

Publication number Publication date
EP2210556A1 (en) 2010-07-28
EP1723906A4 (en) 2009-07-01
US20060106299A1 (en) 2006-05-18
CN102385044B (zh) 2014-12-10
WO2005082245A1 (ja) 2005-09-09
EP1723906A1 (en) 2006-11-22
CN102385044A (zh) 2012-03-21
EP1723906B1 (en) 2019-05-15
EP2409642A1 (en) 2012-01-25
EP2210556B1 (en) 2019-04-10
CN100444791C (zh) 2008-12-24
US7385397B2 (en) 2008-06-10

Similar Documents

Publication Publication Date Title
CN1676097A (zh) 磁共振成象装置及其磁共振成象装置的数据处理方法
CN1905836A (zh) 磁共振成像装置、图像数据修正装置和图像数据修正方法
CN1732847A (zh) 磁共振成像装置和磁共振成像数据处理方法
CN101051030A (zh) 磁共振成像装置以及磁共振成像装置的摄像条件设定方法
CN101077301A (zh) 图像处理装置和磁共振成像装置
CN1618399A (zh) 磁共振成像方法和磁共振成像装置
CN1491095A (zh) 使用高精度线圈感度图的并行mr成象
CN1204412C (zh) 磁共振成象装置
CN1817304A (zh) 高频线圈组件及磁共振成像装置
CN1264480C (zh) 三维反投影方法和一种x射线计算机层析成像装置
CN100337586C (zh) 磁共振成像仪、流量量化仪及流量量化方法
CN1219715A (zh) 用于医学图象的迭代滤波器结构
CN101051388A (zh) 基于复二维奇异谱分析的磁共振部分k数据图像重建方法
CN100339047C (zh) 磁共振成像设备和磁共振成像方法
CN1214622A (zh) 磁共振成像装置和方法
CN1637781A (zh) 图像处理装置、图像处理方法以及图像处理系统
CN101078753A (zh) 采用由多个单元线圈构成的多线圈的并行磁共振成像
CN1491612A (zh) 获取生物电流源的方法和装置
CN1914515A (zh) 磁共振成像设备
CN1225226C (zh) 自旋激励方法和磁共振成像系统
CN1682660A (zh) 计算关于局部血液流量循环指数的方法和装置
CN101040781A (zh) X射线衰减校正方法、ct装置、图像产生装置及方法
CN1635855A (zh) 靶向血管系统的磁共振成像系统和方法
CN1809838A (zh) 用于视频质量评估的方法和系统
CN1636516A (zh) 散射测量方法、散射校正方法、和x射线ct设备

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: 20160714

Address after: Japan Tochigi

Patentee after: TOSHIBA MEDICAL SYSTEMS Corp.

Address before: Tokyo, Japan

Patentee before: Toshiba Corp.

Patentee before: TOSHIBA MEDICAL SYSTEMS Corp.

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: 20081224