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

磁共振成像装置 Download PDF

Info

Publication number
CN107072586B
CN107072586B CN201580011025.4A CN201580011025A CN107072586B CN 107072586 B CN107072586 B CN 107072586B CN 201580011025 A CN201580011025 A CN 201580011025A CN 107072586 B CN107072586 B CN 107072586B
Authority
CN
China
Prior art keywords
parameter
resolution
estimation
image
parameters
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
CN201580011025.4A
Other languages
English (en)
Other versions
CN107072586A (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.)
Fujifilm Healthcare Corp
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Publication of CN107072586A publication Critical patent/CN107072586A/zh
Application granted granted Critical
Publication of CN107072586B publication Critical patent/CN107072586B/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
    • 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/387—Compensation of inhomogeneities
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56536—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to magnetic susceptibility variations
    • 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246—Spatial mapping of the RF magnetic field B1
    • 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
    • G01R33/4835—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 of multiple 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI

Abstract

本发明提供一种技术，在使用同时生成多个参数的计算图像的方法的情况下，也提高计算图像的画质、参数推定精度。因此，本发明使用对各个推定对象的参数最优的分辨率的重构图像，推定该参数的值，得到该参数的值的分布即计算图像。通过调整在摄影时以推定对象的参数中的、最优分辨率最高的推定参数的最优分辨率取得的重构图像的分辨率，而得到最优分辨率的重构图像。或者，在摄影时，只对计算预定的推定对象参数时使用的重构图像，以对该推定对象参数最优的分辨率而取得。

Description

磁共振成像装置
技术领域
本发明涉及一种磁共振成像(MRI：Magnetic Resonance Imaging)技术。特别涉及通过计算推定被检测体参数的技术。
背景技术
MRI装置是一种医用图像诊断装置，其使横断被检测体的任意的平面内的氢原子核引起核磁共振，根据所产生的核磁共振信号(NMR信号：回波信号)摄影核磁共振平面内的断层像。一般，与施加用于确定摄影面的切片倾斜磁场同时地，赋予激励该面内的磁化的激励脉冲，由此得到在所激励的磁化收敛的阶段产生的回波信号。这时，对磁化赋予位置信息，因此在从激励到得到回波信号的期间，在断层面内施加相互垂直的方向的相位编码倾斜磁场和读出倾斜磁场。
根据预先设定的脉冲时序，施加用于产生回波信号的激励脉冲和赋予位置信息的各倾斜磁场。根据目的已知有各种该脉冲时序。例如，在梯度回波(GE) 类型的高速摄影法中，使该脉冲时序重复工作，在每次重复时顺序地改变相位编码倾斜磁场，由此顺序地测量得到一张断层像所需个数的回波。
在MRI中，一般通过重构所得到的回波信号而得到图像。除此以外，有根据使用不同的摄影参数执行脉冲时序所得到的多个图像，通过计算对每个像素求出希望的定量值而得的定量值图像(计算图像、参数分布、参数图表)。在计算的定量值中有被检测体参数、装置参数等。
摄影参数是重复时间、高频磁场的设定强度、高频磁场的相位等在摄影时设定的参数，被检测体参数是纵缓和时间T1、横缓和时间T2、自旋密度ρ、共振频率f0、扩散系数D、高频磁场的照射强度分布(B1分布)等被检测体的物性值、依存于它的值，装置参数是磁场强度B0、接收线圈的灵敏度分布 Sc等依存于摄影所使用的MRI装置的参数。
一般，解析地求出摄影参数、被检测体参数或装置参数与像素值之间的关系作为信号函数，使用该信号函数，得到上述定量值。但是，根据摄影时序的不同，有时没有解析地求出信号函数。在这样的摄影时序的情况下，存在以下的方法，即通过数值模拟生成信号函数，由此生成计算图像(例如参照专利文献1)。
现有技术文献
专利文献
专利文献1：日本特开2011-024926号公报
发明内容
发明要解决的问题
根据推定对象的参数计算图像的最优的分辨率不同。例如，作为被检测体参数的B1分布在空间上平滑地变化，因此与T1、T2分布相比，不需要高分辨率。在分别测量计算各个参数的现有的方法的情况下，能够对各参数以最优的分辨率计算出计算图像。
但是，在专利文献1所公开的方法中，同时生成多个被检测体参数和装置参数的计算图像。因此，最优的分辨率被固定为最高的参数的分辨率，而计算全部的计算图像。例如，在同时计算B1、T1、T2、ρ的分布时，也以与其他分布相同的高分辨率求出B1分布的分辨率。但是，以高分辨率计算出不需要高分辨率的参数的计算图像(例如B1分布)，由此SN比降低，其结果造成计算图像的画质、精度的降低。
本发明就是鉴于上述情况而提出的，其目的在于：提供一种技术，其在使用同时生成多个参数的计算图像的方法的情况下，也提高计算图像的画质、参数推定精度。
用于解决问题的手段
本发明使用对各个推定对象的参数最优的分辨率的重构图像，推定该参数的值，得到该参数的值的分布即计算图像。通过调整在摄影时以推定对象的参数中的最优分辨率最高的推定参数的最优分辨率取得的重构图像的分辨率，而得到最优分辨率的重构图像。或者，在摄影时，只对计算预定的推定对象参数时使用的重构图像，以对该推定对象参数最优的分辨率而取得。
发明效果
根据本发明，在同时生成多个参数的计算图像的情况下，也能够提高计算图像的画质、参数推定精度。
附图说明
图1是表示第一实施方式的MRI装置的概要结构的框图。
图2是第一实施方式的计算机的功能框图。
图3是第一实施方式的计算图像生成处理的流程图。
图4的(a)是用于说明第一实施方式的摄影时序的说明图，(b)是用于说明第一实施方式的摄影时序的k空间的说明图。
图5是第一实施方式的信号函数生成处理的流程图。
图6是用于说明第一实施方式的信号函数的强度分布例子的说明图。
图7是第一实施方式的参数推定处理的流程图。
图8是用于说明第一实施方式的参数推定处理的流程的具体例子的说明图。
图9的(a)是用于说明在图8所示的处理中最初得到的重构图像例子的说明图，(b)～(e)是用于说明在图8所示的处理的中途得到的重构图像或参数分布例子的说明图，(f)是说明通过现有方法得到的参数分布例子的说明图。
图10(a)～(d)是用于分别说明通过第一实施方式的方法得到的参数分布、通过现有方法得到的参数分布的说明图。
图11的(a)是用于说明第二实施方式的单回波时序的说明图，(b)是用于说明第二实施方式的多回波时序的说明图。
图12(a)～(c)是用于说明第二实施方式的参数推定处理的说明图。
图13是用于说明本发明的各实施方式的接受画面例子的说明图。
图14的(a)是用于说明本发明的变形例子的摄影时序的说明图，(b)是用于说明(a)的摄影时序的k空间的说明图。
具体实施方式
<第一实施方式>
以下，说明应用本发明的第一实施方式。以下，在用于说明本发明的实施方式的全部图中，基本上具有相同功能的部分附加相同的符号，省略其重复的说明。
首先，说明本实施方式的MRI装置。图1是表示本实施方式的MRI装置 100的概要结构的框图。MRI装置100具备产生静磁场的磁铁101、产生倾斜磁场的倾斜磁场线圈102、时序发生器104、倾斜磁场电源105、高频磁场发生器106、照射高频磁场并且检测核磁共振信号的收发线圈107、接收器108、计算机109、显示装置110、存储装置111。收发线圈107在图中表示出一个，但也可以分别地具备发送线圈和接收线圈。
被检测体(例如生物体)103被载置在卧台(台(table))等上，配置在磁铁101所产生的静磁场空间内。
时序发生器104向倾斜磁场电源105和高频磁场发生器106发送指令，使得分别产生倾斜磁场和高频磁场。通过收发线圈107向被检测体103施加高频磁场。通过收发线圈107接收从被检测体103产生的核磁共振信号(NMR信号：回波信号)，通过接收器108进行检波。
通过时序发生器104设置成为检波的基准的核磁共振频率(检波基准频率；共振频率f0)。检波出的回波信号被发送到计算机109，在此进行图像重构等信号处理。其结果被显示在显示装置110。也可以根据需要，将检波出的信号、测定条件存储到存储装置111。
时序发生器104通常控制各装置以预先编程的定时、强度动作。将程序中的特别是记述了高频磁场、倾斜磁场、信号接收的定时、强度的部分称为脉冲时序(摄影时序)。
计算机109依照脉冲时序，使各部动作，测量回波信号。另外，针对得到的回波信号，实施各种信号处理，得到希望的图像。此外，通过时序发生器 104进行回波信号的测量。
计算机109具备CPU和存储器。另外，由计算机109的CPU将存储在存储装置111中的程序装载到存储器并执行，来实现计算机109实现的各功能。另外，也可以通过ASIC(Application Specific Integrated Circuit：专用集成电路)、 FPGA(field-programmable gate array：现场可编程门阵列)等硬件来实现全部或一部分的功能。另外，将各功能的处理所使用的各种数据、处理中生成的各种数据存储在存储装置111中。
此外，磁铁101作为静磁场产生部发挥功能，倾斜磁场线圈102和倾斜磁场电源105作为倾斜磁场产生部发挥功能，收发线圈107和高频磁场发生器 106作为高频发送部发挥功能，收发线圈107和接收器108作为接收部发挥功能，时序发生器104和计算机109作为控制部发挥功能，计算机109还作为图像重构部发挥功能。
在本实施方式中，根据改变摄影参数而取得的多个重构图像，分别以最优分辨率得到被检测体参数和/或装置参数的计算图像。为了实现它，如图2所示，本实施方式的计算机109具备接受部210、重构图像取得部220、参数推定部230、显示图像生成部240。另外，如后述那样，在需要生成信号函数的情况下，还具备信号函数生成部250。
以下，在本实施方式中，在摄影时序执行时，用户将能够在任意设定的参数称为摄影参数，将依存于被检测体103的参数称为被检测体参数，将MRI 装置固有的参数称为装置参数。
摄影参数例如有重复时间(TR)、回波时间(TE)、高频磁场的强度(倾斜角(FlipAngle：FA))、高频磁场的相位(θ)等。另外，被检测体参数有纵缓和时间(T1)、横缓和时间(T2)、自旋密度(ρ)、共振频率差(Δf0)、扩散系数(D)、高频磁场强度(高频磁场的照射强度分布：B1)等。进而，装置参数有静磁场强度(B0)、接收线圈的灵敏度分布(Sc)等。此外，共振频率差Δf0是各像素的共振频率和基准频率f0之间的差。
[接受部]
接受部210接受来自用户的指示。接受部210在重构图像取得部220执行摄影时，接受必要的信息。进而，在本实施方式中，接受用于生成计算图像(参数分布)的对象的被检测体参数和/或装置参数(推定参数)的指定、摄影所使用的摄影参数组的指定。经由显示在显示装置110的接受画面，接受这些信息。另外，预先将能够指定的参数的信息保存在存储装置111中。
[重构图像取得部]
本实施方式的重构图像取得部220依照预先确定的摄影时序，使用预先设定的摄影参数，执行摄影，得到重构图像。通过向时序发生器104发出指示以便控制各部依照摄影参数和摄影时序而执行摄影。
在本实施方式中，使用一个以上的摄影参数的值不同的多个摄影参数组，依照预先确定的摄影时序，进行摄影，针对每个摄影参数组取得重构图像。
这时，本实施方式的重构图像取得部220依照与由接受部210接受的各推定参数对应的摄影条件，进行摄影，取得重构图像。此外，与各推定参数对应的摄影条件包括对该推定参数最优的分辨率，预先保存在存储装置111等中。
[参数推定部]
参数推定部230使用多个重构图像、依照取得该重构图像时的摄影时序的信号函数，推定依存于被检测体103的被检测体参数、以及依存于装置的装置参数中的至少一个参数即推定对象的参数(推定参数)的值。在本实施方式中，使用对每个推定参数最优的分辨率的重构图像，对推定参数的值进行推定。
对每个像素进行推定参数的值的推定。因此，作为结果，本实施方式的参数推定部230输出每个像素的参数值即参数分布。将在后面说明本实施方式的参数推定部230的参数推定处理的细节。
[显示图像生成部]
显示图像生成部240根据由参数推定部230生成的参数分布生成显示图像，显示在显示装置110。此外，也可以不具备显示图像生成部240。
首先，在图3中表示这些各部的计算图像生成处理的流程的概要。
首先，重构图像取得部220依照预先确定的摄影时序，进行摄影，得到多个重构图像320(步骤S1101)。这时，在本实施方式中，如上述那样，使用一个以上的摄影参数的值不同的多个摄影参数组310，进行摄影。然后，参数推定部230使用信号函数330和多个重构图像320，进行参数推定处理(步骤 S1102)，得到希望的推定参数的每个像素的参数值(参数分布)340。然后，显示图像生成部240生成显示图像(步骤S1103)，显示在显示装置110(步骤S1104)，结束处理。
[参数推定处理]
以下，说明本实施方式的参数推定处理。
如上述那样，参数推定部230在参数推定时使用信号函数330。在说明参数推定处理之前，说明信号函数330。
信号函数330是在对每个摄影时序确定摄影参数、被检测体参数以及装置参数与重构图像的各像素值之间的关系的函数。一般，对每个摄影时序解析地计算并公式化。在参数推定时使用具有公式化了的信号函数的摄影时序的情况下，使用该信号函数330。
另一方面，在使用信号函数没有公式化的摄影时序的情况下，在本实施方式中，信号函数生成部250预先生成信号函数330。即，本实施方式的信号函数生成部250在生成计算图像时使用的摄影时序是无法解析地求出信号函数的时序、或不是已知的时序的情况下，通过数值模拟生成信号函数330。此外，在使用解析地求出了信号函数的摄影时序、或已知的时序的情况下，也可以不具备信号函数生成部250。
[摄影时序]
在说明信号函数生成部250的信号函数生成处理的细节之前，说明没有解析地求出信号函数的摄影时序的例子。在此，作为该摄影时序，列举 RF-Sppoiled GE时序为例子进行说明。
图4(a)是RF-Spoiled GE时序610的时序图。RF、Gs、Gp、Gr分别表示高频磁场(RF)脉冲、切片选择倾斜磁场脉冲、相位编码倾斜磁场脉冲、读出倾斜磁场脉冲的施加定时。以下，在本说明书中相同。
在该RF-Spoiled GE时序610中，首先施加切片倾斜磁场脉冲611，并且照射高频磁场(RF)脉冲612，激励对象物体内的预定的切片的磁化。接着，施加切片相位倾斜磁场脉冲613、用于对磁化的相位附加相位编码方向的位置信息的相位编码倾斜磁场脉冲614、失相(dephase)用读出倾斜磁场615。然后，一边施加用于附加读出方向的位置信息的读出倾斜磁场脉冲616，一边测量磁共振信号(回波信号)617。然后，最后施加失相用相位编码倾斜磁场脉冲619。
在使用RF-Spoiled GE时序610进行摄影的情况下，重构图像取得部220 使相位编码倾斜磁场脉冲614和失相用相位编码倾斜磁场脉冲619的强度(相位编码量kp)改变，并且使RF脉冲612的相位的增加值每次变化117度，同时以重复时间TR重复以上的步骤，测量得到一张图像所需要的回波。此外，这时，第n个施加的RF脉冲612的相位θ(n)为θ(n-1)+117n。
由此，如图4(b)所示，将各回波信号配置在k空间上。然后，重构图像取得部220通过二维傅立叶逆变换而重构图像。此外，如果使用RF-Spoiled GE时序610，则得到强调了T1(纵缓和时间)的重构图像。
[信号函数生成处理]
以下，说明信号函数生成部250生成RF-Spoiled GE时序610的信号函数的处理。
如图5所示，信号函数生成部250分别使用一个以上的摄影参数的值不同的多个摄影参数组410，进行数值模拟(步骤S1201)，生成信号函数420。多个摄影参数组410是多个不同的上述摄影参数和上述被检测体参数的组合。
[信号函数]
使用被检测体参数的纵缓和时间T1、横缓和时间T2、自旋密度ρ、高频磁场强度B1、装置参数的接收线圈灵敏度Sc、作为能够在RF-Spoiled GE时序 610中变更的摄影参数的倾斜角FA、重复时间TR、回波时间TE、RF相位增加值θ，用以下的公式(1)表示RF-Spoiled GE时序610的信号函数fs。
[公式1]
I＝fs(ρ，T1，T2，B1，FA，TR，TE，θ，Sc)＝ρScf(T1，T2，B1×FA，θ，TR，TE)…(1)
此外，B1在摄影时为FA的系数，因此为与FA的积的形式。另外，ρ和 Sc作为比例系数对信号强度起作用，因此处于函数的外侧。
[摄影参数组]
本实施方式的信号函数生成部250对被检测体参数的T1、T2的任意的值涵盖地改变摄影参数FA、TR、TE，通过数值模拟得到信号，对所得的信号进行插补，由此生成信号函数。
此外，能够在RF-Spoiled GE时序610中变更的参数如上述那样是FA、 TR、TE、θ。但是，其中，RF相位增加值θ被固定为117度。θ被固定是为了得到与作为高速摄影法之一的FLASH同等的T2依存性少的图像对比度。如果改变θ，则图像对比度的T2依存性有很大变化。
另外，这时，摄影对象的自旋密度ρ和B1、Sc是固定的(例如为1)。
以下表示生成信号函数时使用的被检测体参数T1、T2的值、改变的摄影参数FA、TR、θ的值的一个例子。分别包含实际摄影时使用的摄影参数的范围、被检测体的T1、T2的范围。
T1(15个)[s]：0.05，0.07，0.1，0.14，0.19，0.27，0.38，0.53，0.74， 1.0，1.5，2.0，2.8，4.0，5.6。
T2(17个)[s]：0.01，0.02，0.03，0.04，0.05，0.07，0.1，0.14，0.19， 0.27，0.38，0.53，0.74，1.0，1.4，2.0，2.8。
TR(4个)[ms]：10，20，30，40。
FA(10个)度：5，10，15，20，25，30，35，40，50，60
θ(17个)度：0，1，2，3，4，5，6，7，8，9，10，12，14，16，18， 20，22。
在该情况下，信号函数生成部250生成由以上的摄影参数和被检测体参数的全部组合构成的173400个摄影参数组410，通过数值模拟计算出基于各个摄影参数组410的信号值。
[数值模拟]
数值模拟以在格子点上配置自旋所得的被检测体模型、摄影时序、摄影参数、装置参数作为输入，求解作为磁共振现象的基础方程式的Bloch的公式，输出磁共振信号。
作为自旋的空间分布(γ，M0，T1，T2，Cs)而赋予被检测体模型。此处，γ是磁旋比，M0是热平衡磁化(自旋密度)，T1和T2分别是纵缓和时间和横缓和时间。
Bloch的公式是1阶线性常微分方程式，用下式表示。
[公式2]
Figure BDA0001094301550000091
H＝Bo+GxX+Gyy+Gzz+2πfo/γ
在此，(x，y，z)表示三维的正交坐标系，z与静磁场(强度为B0)的方向相等。另外，(Mx，My，Mz)是自旋，Gx、Gy、Gz分别是下标方向的倾斜磁场强度，H1是高频磁场强度，f0是旋转坐标系的频率。
[插补]
如果通过上述公式(2)得到每个摄影参数组410的信号值，则信号函数生成部250对它们进行插补，生成信号函数fs(420)。对插补，可以使用1次到3次左右的线性插补、样条插补。
在图6中表示本实施方式的信号函数生成部250通过上述方法生成的信号函数420的强度的一部分。在此，表示在设T1为900ms，T2为100ms，θ为 5度的情况下的FA和TR的每个值的强度。横轴是FA(deg.)，纵轴是TR(s)。信号强度越大，则用越浓的颜色表示。
[参数推定处理]
接着，说明本实施方式的参数推定部230的参数推定处理的细节。如上述那样，在本实施方式中，对每个推定参数，使用最优分辨率的重构图像，对推定参数的值进行推定。如上述那样，预先将每个推定参数的最优分辨率保存在存储装置111等中。
在推定时，参数推定部230使用多个重构图像320、摄影时序的信号函数330。重构图像取得部220使用多个摄影参数组310，依照预先确定的摄影时序进行摄像，分别取得该多个重构图像320。在多个摄影参数组中，也可以包含摄影参数相同的参数组。
这时，取得的重构图像320的张数为推定对象的推定参数个数以上。另外，以成为推定对象的推定参数中的最优分辨率最高的推定参数的分辨率(最高分辨率)取得要取得的多个重构图像320。
例如，在推定被检测体参数的T1、T2、B1、ρ和装置参数Sc的积即a(＝ρSc)的情况下，对于B1，如上述那样，在空间上平滑地变化，因此与T1、 T2相比，不需要高分辨率。因此，在该情况下，取得推定T1、T2所需要的分辨率的图像。
此外，根据摄像视野FOV和编码步进数确定分辨率。经由接受部210设定这些摄影参数。
另外，参数推定部230针对最优分辨率比最高分辨率低的推定参数，将重构图像取得部220取得的重构图像的分辨率降低为该推定参数的最优分辨率，对推定参数的值进行推定。
这时，参数推定部230从最优分辨率低的推定参数开始顺序地对推定参数的值进行推定。在进行推定时，将已经推定的推定参数的值调整为推定对象的推定参数的最优分辨率，与信号函数330一起使用。
以下，依照处理的流程说明本实施方式的参数推定部230的参数推定处理的细节。图7是本实施方式的参数推定处理的处理流程。在此，将推定参数个数设为M(M为2以上的整数)。
在该情况下，预先准备多个摄影参数组(FA、TR、θ)310，在上述步骤 S1101中，重构图像取得部220取得所准备的摄影参数组310的个数的重构图像320。
准备的摄影参数组310的个数一般需要作为推定参数个数的M以上。这是因为在求解方程式的解时，方程式(摄影参数组和根据它摄影的图像)的个数必须与未知数(推定参数)的个数相等或比其更多。在此，作为一个例子，假设准备N(N为M以上的整数)组。因此，重构图像320的个数也是N张。
例如，在推定参数个数是4个(T1、T2、B1、a)的情况下，摄影参数组 310的个数为比4组多的例如6组。此外，摄影参数组个数越多，则越是提高推定精度，但摄影时间变长。另外，在此，以推定对象的推定参数中的要求最高分辨率的推定参数的分辨率，取得多个重构图像320。
参数推定部230推定M个推定参数。在本实施方式中，从最优分辨率小的推定参数开始顺序地进行推定(步骤S1301)。另外，同时对最优分辨率相同的推定参数进行推定。以下，将最优分辨率为第m个低的推定参数个数设为k(m)个。
首先，将N张重构图像320调整为最优分辨率为第m个低的k(m)个推定参数的最优分辨率(步骤S1302)。在此，降低N张重构图像320的分辨率。例如通过将多个像素的值相加等，来降低分辨率。
此外，在降低分辨率时，根据重构图像320的SN比和对处理对象的推定参数最优的分辨率，决定设为怎样程度的分辨率。即，在重构图像320的SN 比充分高的情况下，可以任意地决定分辨率。但是，在SN比低的情况下，与推定参数的分布对应地设为尽量低的分辨率。由此，即使SN比低，也能够尽量保持参数的推定精度高。
将已经推定的推定参数的参数分布的分辨率调整为与在步骤S1302中调整后的图像的分辨率一致(步骤S1303)。在此，先推定出的推定参数的参数分布的分辨率为当前正在推定的推定参数以下，因此进行插补扩大处理，调整分辨率。在插补放大时，使用以下的公式(3)所示的预先确定的次数n以下的同次函数的和。例如，理想的是次数n为6左右。
[公式3]
Figure BDA0001094301550000111
另外，使用在步骤S1302中得到的分辨率的图像、已经推定出的推定参数的已经调整分辨率后的参数分布、先生成出的信号函数330，对推定对象的参数进行推定(步骤S1304)。其结果是得到该推定参数的参数分布340。
具体地说，通过使每个像素的信号值I与对公式(1)进行变形所得的以下的公式(4)的函数f拟合，由此来推定推定参数。
[公式4]
Figure BDA0001094301550000121
此外，通过用以下的公式(5)表示的最小二乘法来进行函数拟合。
[公式5]
Figure BDA0001094301550000122
在此，χ是信号函数和幻像(phantom)的像素值的残差的总和，I是(FA，θ，TR)处的像素值。
对全部推定参数重复进行以上的处理(步骤S1305、S1306)，结束处理。
<实施例>
以下，使用图8、图9(a)～图9(f)通过具体例子说明本实施方式的参数以上的参数推定部230的参数推定处理。
在此，将推定参数设为(B1、a、T1、T2)这4个。因此，准备的摄影参数组310的个数为比其多的6组。将各摄影参数组310分别设为P1～P6。
P1～P6的摄影参数组310为分别通过FA为10度和30度、θ为2度、6 度、20度、TR为10ms、40ms的组合所得到的12组摄像参数组中的6组。此外，TE固定为5ms。此外，摄影参数组310例如依照误差传输的法则等，选择能够高精度地推定的组。
另外，推定对象的推定参数中只有B1是最优分辨率比其他推定参数低的推定参数。即，在本例子中，最优分辨率为第一个低的推定参数只有B1一个，最优分辨率为第二个低的推定参数是a、T1、T2这3个。
因此，首先，重构图像取得部220与a、T1、T2的最优分辨率对应地取得6张重构图像320。在图9(a)中表示重构图像取得部220使用各摄影参数组310得到的6张重构图像320的例子。
参数推定部230首先推定B1的参数值(B1分布)(步骤S1410)。这时，与B1的最优分辨率对应地使图像低分辨率化并推定。
在B1分布的推定中，首先，将6张重构图像320的分辨率降低为对推定参数B1最优的分辨率(步骤S1401)。在图9(b)中表示分别使图9(a)所示的重构图像320低分辨率化后的图像321。
图9(b)所示的图像321将原重构图像320的4×4像素相加，成为一个像素。例如在原重构图像320的分辨率是1.6mm×1.6mm的情况下，将4×4 像素相加成为一个像素，为6.4mm×6.4mm，得到图像321。
然后，使用信号函数330、低分辨率化后的6张图像321，推定B1的每个像素的参数值(步骤S1402)，得到参数分布(B1分布)341。此外，这时，还以与B1相同的分辨率推定a、T1、T2。在图9(c)中表示它们的推定结果341。
接着，参数推定部230推定剩余的推定参数(T1、T2、a)(步骤S1420)。在此，使用插补扩大了的B1分布和取得时的分辨率的多个重构图像320。
因此，首先对B1分布341进行插补扩大，将其分辨率设为与原6张重构图像320的分辨率相等(步骤S1421)。在图9(d)中表示插补扩大了的B1分布342。
另外，使用图9(a)所示的原重构图像320、插补扩大后的B1分布342、信号函数330，以原重构图像320的分辨率推定a、T1、T2、a的3个参数(步骤S1402)，得到参数分布。在图9(e)中表示插补扩大后的B1分布342、推定出的a、T1、T2的参数分布343。
此外，在步骤S1402中的推定处理中，对于B1，使用插补放大后的B1分布342作为正解值，因此，成为已知，未知数是a、T1、T2的3个。未知数的个数减少，因此处理速度变快，推定精度提高。
在此，在图9(f)中表示现有的方法，即根据以与最优分辨率最高的推定参数对应的分辨率得到的重构图像推定全部推定参数(B1、a、T1、T2)的分布的情况下的结果(参数分布)344。
如果对图9(e)所示的参数分布343和图9(f)所示的参数分布344进行比较，则可知本实施方式的方法的结果343与现有技术相比，噪声减少了。
为了更明确地表示效果，在图10(a)和图10(b)中分别表示图9(e) 的B1分布和a分布343的纵向的一行(横向的大致中心位置的行)的轮廓513、 523，另外，在图10(c)和图10(d)中分别表示图9(f)的B1分布和a分布344的同一轮廓514、524。在本图中，横轴是一行上的位置，纵轴是信号强度。
如该图所示，对于B1分布，如果对本实施方式的方法的轮廓513和现有方法的轮廓514进行比较，则被插补扩大了，因此可知在本实施方式的方法的分布中几乎没有噪声。对于a分布，如果对本实施方式的方法的轮廓523和现有方法的轮廓524进行比较，则也可知通过本实施方式的方法得到的分布的噪声少。
如以上说明的那样，本实施方式的MRI装置100具备：参数推定部230，其使用多个重构图像320和取得该重构图像320时所遵从的摄影时序的信号函数330，推定依存于被检测体103的被检测体参数和依存于装置的装置参数的至少一个推定参数的值，其中，使用多个摄影参数组310，依照上述摄影时序进行摄像，分别取得上述多个重构图像320，上述信号函数330是对每个上述摄影时序确定上述摄影参数、上述被检测体参数以及上述装置参数与上述重构图像的各像素值之间的关系的函数，上述参数推定部230使用对每个上述推定参数最优的分辨率的上述重构图像来推定上述值。
这样，根据本实施方式，针对推定对象的参数，使用对它们分别最优的分辨率的图像，推定其值。因此，能够高速、高精度、SN比不降低地得到计算图像。
另外，也可以是上述多个重构图像320具有上述推定参数中的上述最优分辨率最高的推定参数的该最优分辨率即最高分辨率，上述参数推定部230针对上述最优分辨率比上述最高分辨率低的推定参数，将上述重构图像320的分辨率降低为对该推定参数最优的分辨率，来推定上述值。
这时，上述参数推定部230从上述最优分辨率低的上述推定参数开始顺序地推定上述值，并且在推定时，将已经推定出的推定参数的值调整为推定对象的上述推定参数的最优分辨率，与上述信号函数330一起使用。
因此，根据本实施方式，根据计算量少的最优分辨率低的推定参数进行推定，并且随着最优分辨率变高，未知的推定对象参数个数变少。因此，处理中的发散等降低，能够得到高精度的推定结果。
<第二实施方式>
接着，说明本发明的第二实施方式。在第一实施方式中，以与最优分辨率最高的推定参数的计算图像一致的分辨率取得图像，在推定参数值时，分别降低为最优分辨率。另一方面，在本实施方式中，在对每个推定参数有帮助的重构图像清楚的情况下，在取得重构图像时，以最优分辨率取得该重构图像。
本实施方式的MRI装置基本上具有与第一实施方式的MRI装置100相同的结构。但是，参数推定时的处理不同，因此计算机109的各部的处理内容不同。以下，着眼于与第一实施方式不同的处理来说明本实施方式。
此外，能够使用独立的图像计算定量值的推定参数例如有共振频率差Δ f0。在计算磁化率图像时使用Δf0。
即，Δf0与相位为比例关系。因此，通过将预定的系数乘以Δf0分布的各像素值，能够得到相位图像。另外，相位图像反映了组织间的磁化率。因此，通过基于相位图像来求解逆问题，能够得到磁化率图表。
如果有TE不同的2张以上的图像，则能够计算出该Δf0分布。在计算Δ f0分布时，与T1、T2、a等参数相比，要求高分辨率。
因此，在本实施方式中，首先得到TE不同的2张第一重构图像，根据该 2张图像计算出Δf0分布。然后，以与第一重构图像不同的分辨率，改变摄影参数地得到多个第二重构图像，通过与第一实施方式相同的方法，推定其他推定参数。这时，将第一重构图像变换为第二重构图像的分辨率地使用第一重构图像。另外，在本实施方式中，通过多回波时序取得TE不同的2张第一重构图像。
即，本实施方式的重构图像取得部220取得推定参数中最优分辨率最高的推定参数的具有作为最优分辨率的第一分辨率的第一重构图像、具有与第一分辨率不同的第二分辨率的第二重构图像。该第一重构图像是通过多回波时序取得的多回波图像。此外，第二分辨率为在对其他推定参数最优的分辨率中的最高的分辨率。
然后，本实施方式的参数推定部230根据多回波图像，对推定参数中的Δ f0的值进行推定。然后，参数推定部230将多回波图像的分辨率调整为第二分辨率，与信号函数一起使用调整后的多回波图像和第二重构图像，推定Δf0以外的推定参数的值。
以下，在详细说明本实施方式的参数推定部230的处理之前，说明取得 TE不同的图像的摄像时序。在本实施方式中，也与第一实施方式同样地，基本上使用RF-Spoiled GE时序610。在取得TE不同的图像的情况下，使用在一次TR中测量2个以上的回波信号的多回波时序。
图11(a)是单回波时序的RF-Spoiled GE时序(以下称为单回波时序。) 610，图11(b)是多回波时序的RF-Spoiled GE时序(以下称为多回波时序。) 630。此外，在此，作为多回波时序630，示例在一次TR内取得2个回波信号的情况。此外，在本图中，A/D表示回波信号取得定时。
图11(a)所示的单回波时序610如图4(a)所说明的那样。其中，在图 11(a)中表示三维摄像的情况下的时序。因此，613是切片相位兼切片相位编码倾斜磁场脉冲。
重构图像取得部220改变相位编码倾斜磁场脉冲613和614的强度(相位编码量ks、kp)，并且使RF脉冲612的相位的增加值每次变化117度(第n 个RF脉冲的相位为θ(n)＝θ(n-1)+117n)，同时以重复时间TR重复进行以上的步骤，测量得到图像所需要的回波。各回波被配置在三维的k空间 (kr-kp-ks空间)上。
然后，重构图像取得部220通过三维傅立叶逆变换而重构图像。通过该脉冲时序610得到的图像是强调了T1(纵缓和时间)的图像。
对于图11(b)所示的多回波时序630，所施加的脉冲基本上与单回波时序610的情况相同。即，首先施加切片倾斜磁场脉冲611并且照射高频磁场 (RF)脉冲612，激励对象物体内的某切片的磁化。接着，在施加切片相位兼切片相位编码脉冲613、相位编码倾斜磁场脉冲614、失相用读出倾斜磁场脉冲615后，一边施加用于附加读出方向的位置信息的读出倾斜磁场脉冲616，一边测量第一磁共振信号(回波信号)617。
在测量回波信号617后，接着，一边施加读出倾斜磁场脉冲618，一边测量第二个回波信号622。最后，施加失相用相位编码倾斜磁场脉冲619、620、以及破碎脉冲621。
在本时序630中，回波信号617和回波信号622相对于RF脉冲612的时间(TE)不同。
重构图像取得部220与单回波时序610同样地控制动作，分别将TE不同的2个回波信号617、622配置在不同的三维的k空间。然后，重构图像取得部220通过三维傅立叶逆变换，从各k空间分别重构图像。重构出的2张图像的TE不同，因此T2强调度不同。
[参数推定处理]
参数推定部230首先根据该TE不同的2张图像，计算Δf0分布。
被激励的自旋的横磁化的旋转量(相位)是将TE乘以该横磁化的位置处的频率和成为基准的共振频率f0之间的差(Δf0)所得的值。因此，如果用2 张图像的相位差除以TE的差，则得到Δf0。对每个像素求出Δf0，得到Δf0分布。
此外，TE不同的图像的个数并不限于2个。例如，也可以使用进一步追加读出脉冲而测量3个以上的TE不同的回波的多回波时序，进而增加TE不同的图像而用于Δf0分布的计算。如果增加图像个数，则Δf0分布的计算精度提高。因此，一般使用在一次TR中测量3个以上(通常是5个左右)的回波信号的时序。
本实施方式的参数推定部230推定包含Δf0的预先确定的推定参数。
如上述那样，在第一实施方式中，参数推定部230进行参数推定处理，因此分别改变一个以上的摄影参数的值地取得推定参数个数以上的图像。但是，在本实施方式中，为了计算Δf0分布而执行多回波时序630。通过该多回波时序630，得到摄影参数相同但TE不同的2张以上的图像。如果为了计算Δf0分布而设执行多回波时序630所得到的图像(称为多回波图像)的张数为J张，推定参数个数为M，则在本实施方式中，依照单回波时序610取得M-J张以上的图像(单图像)。
这时，与Δf0以外的推定参数中其参数分布(计算图像)的最优分辨率最高的推定参数的最优分辨率一致地取得单图像。
此后的处理与第一实施方式相同。即，从最优分辨率低的推定参数开始，顺序地对推定参数的值进行推定，得到参数分布。在推定时，使用预先准备的信号函数fs、J张多回波图像、(M-J)张单图像。这时，对于J张多回波图像，与单图像的分辨率一致地调整分辨率(在此使其降低)。
此外，例如在本实施方式中，在对Δf0以外的推定参数最优的分辨率是相同水平的情况，或多回波图像、单图像都具有充分的SN比的情况下，也可以不使单图像的分辨率变化地推定Δf0以外的推定参数。
[参数推定处理的流程]
使用图12(a)～图12(c)说明本实施方式的参数推定处理的流程。在此，推定参数为T1、T2、B1、ρSc(＝a)、Δf0的5个。另外，列举以下的情况为例子进行说明，即使用单回波时序610和多回波时序630，取得各推定参数的三维分布。在多回波时序630中，得到TE不同的2张多回波图像。
此外，通过在第一实施方式中说明的方法等，预先生成信号函数fs。
在该情况下，也与第一实施方式同样地，准备比推定的推定参数个数多的个数的FA、TR、θ的任意一个的值不同的摄影参数组310。即，从根据FA 为10度和30度、θ为2度、6度、20度、TR为10ms、40ms的组合得到的 12组摄像参数组中，与第一实施方式同样地进行选择。但是，在本实施方式中，还通过多回波时序630取得2张多回波图像。因此，在第一实施方式中，准备了6组摄影参数组(P1～P6)310，但在此准备5组摄影参数组(P1～P5)。由此，得到6张图像。
另外，在多回波时序630中，如果TR不适度地变长，则无法测量TE长的回波。因此，多回波时序630所使用的摄影参数组310为长的TR的组。例如，在准备了TR为10ms和40ms的摄影参数组310的情况下，应用于多回波时序630的摄影参数组310是TR为40ms的参数组。
另外，在第一实施方式的例子中，TE固定为5ms。在本实施方式中，在多回波时序630中，将第一回波的TE(TE1)设为与单回波时序610相同的 5ms，将第二回波的TE(TE2)设为20ms。
进而，如上述那样，Δf0在此后用于磁化率分布的计算。在该磁化率分布的计算中，需要比通常的形态图像高的分辨率。因此，设定摄影参数使得通过多回波时序630得到的图像(多回波图像)的分辨率比通过单回波时序610 得到的图像(单回波图像)高。
例如，为了得到根据Δf0计算磁化率分布的精度，理想的是以表示多回波图像的分辨率的体素大小为从0.5mm到1mm左右的方式来设定摄影参数。在此，为0.8mm。另一方面，单回波图像的分辨率为1.6mm。这是其他推定参数T1、T2、B1、a的最优分辨率中最高的。
此外，将执行多回波时序630所使用的摄影参数组设为P1，将执行单回波时序610所使用的摄影参数组设为P2～P5。
在图12(a)中表示使用这些摄影参数组得到的图像。在该图中，图像711 和712是通过多回波时序630得到的多回波图像，分别为在TE1、TE2中得到的多回波图像。
另外，图像713～716分别是通过以摄影参数组P2～P5执行单回波时序610 而得到的单回波图像。这些图像713～716是在TE1中得到的图像。
本实施方式的参数推定部230首先使用多回波图像711、712推定Δf0。
接着，参数推定部230如图12(b)所示那样，将多回波图像711、712 的分辨率变换为与单回波图像713～716相同的分辨率。在此，降低分辨率。设降低后的多回波图像为721、722。
在此，将多回波图像711、712的体素大小设定为单回波图像713～716的 2倍，因此例如通过简单地将多回波图像711、712的2×2×2体素相加而成为1个体素，来实现分辨率的降低。除此以外，也可以使用线性插补、2次到 3次左右的插补。
然后，参数推定部230使用预先准备的信号函数fs、多回波图像721、722、单回波图像713～716，对推定参数B1、a、T1、T2进行推定。这些推定参数的推定方法与第一实施方式相同。例如，使用上述公式(4)和公式(5)，通过最小二乘拟合，推定这些推定参数。这时，也可以与第一实施方式同样地，首先以低分辨率求出B1图表后，以原分辨率求出T1、T2、a。
通过以上，如图12(c)所示，得到Δf0、B1、a、T1、T2各自的参数分布(计算图像)740。
如以上说明的那样，本实施方式的MRI装置100具备与第一实施方式相同的参数推定部230。另外，在对推定参数的值进行推定时使用的上述多个重构图像由具有上述推定参数中的上述最优分辨率最高的推定参数的该最优分辨率即第一分辨率的第一重构图像(711、712)、具有与上述第一分辨率不同的第二分辨率的第二重构图像(713～716)构成。
另外，上述第一重构图像(711、712)也可以是通过多回波时序取得的多回波图像，参数推定部230根据上述多回波图像，推定上述推定参数中的共振频率差Δf0的值。另外，上述参数推定部230将上述多回波图像的分辨率调整为上述第二分辨率，与上述信号函数一起使用调整后的多回波图像和上述第二重构图像(713～716)，推定上述共振频率差以外的上述推定参数的值。
这时，上述参数推定部230也可以针对上述最优分辨率比上述第二分辨率低的推定参数，将上述重构图像的分辨率降低为对该推定参数最优的分辨率来推定上述值。进而，也可以针对Δf0以外的推定参数，从上述最优分辨率低的上述推定参数开始顺序地推定上述值，并且在推定时，将已经推定出的推定参数的值调整为推定对象的上述推定参数的最优分辨率，与上述信号函数一起使用。
这样，根据本实施方式，与第一实施方式同样地，能够以必要充分的分辨率得到希望的定量值(被检测体参数、装置参数的值)。因此，能够得到高质量的计算图像。另外，能够抑制推定时的发散产生次数，能够高效地推定值。
另外，在本实施方式中，根据通过多回波时序取得的高分辨率的图像计算Δf0，另外将高分辨率的图像设为低分辨率，与其他低分辨率图像一起使用，推定T1、T2、B1、a这样的推定参数。通过这样使用多回波时序，能够减少取得的参数推定用的图像个数。与分别个别地进行摄影、推定相比，能够提高摄影效率。
此外，上述各实施方式也可以具备从用户接受推定对象的参数的指定的接受画面。由接受部210生成接受画面，显示在显示装置110。在图13中表示该接受画面800的例子。
如该图所示，接受画面800例如具备接受推定对象的推定参数的指定的推定参数接受区域810、接受所设定的摄像参数组的参数组接受区域820。
接受部210接受来自用户的指示，向重构图像取得部220和参数推定部 230通知该指示。重构图像取得部220依照所接受的指示取得图像，参数推定部230依照该指示推定预定的推定参数。
这时，例如也可以针对Δf0等预先决定了对推定最优的摄影时序的推定参数，对应地保存摄影时序。由此，在用户选择了该推定参数的情况下，即使不设定摄影时序，也自动地选择对应地保存的摄影时序。
<变形例>
上述各实施方式的显示图像生成部240能够使用所得到的计算图像和信号函数fs，任意地生成变更了推定出的被检测体参数和装置参数、摄影时序、摄影参数的一部分的图像(例如T1强调像、T2强调像)。
其中，能够变更的摄影时序是信号函数fs已知、或已经通过数值模拟生成出的任意的摄影时序。
以下，作为摄影时序，以使用最普通的自旋回波(SE)时序的情况为例子进行说明。首先，说明SE时序。图14(a)是SE时序900的例子。
在SE时序900中，首先施加切片选择倾斜磁场脉冲901，并且照射高频磁场(RF)脉冲902，激励对象物体内的某切片的磁化。接着，施加切片相位倾斜磁场脉冲903、用于对磁化的相位附加相位编码方向的位置信息的相位编码倾斜磁场脉冲904、失相用读出倾斜磁场905。然后，在对各轴施加用于抑制不需要的信号的破碎脉冲906、907、908后，与切片选择倾斜磁场脉冲909 一起照射重聚焦脉冲910。再次施加破碎脉冲911、912、913。然后，一边施加用于附加读出方向的位置信息的读出倾斜磁场脉冲914，一边通过A/D915 测量回波信号918，最后施加破碎脉冲916和917。将从RF脉冲902的照射到回波峰值的时间称为回波时间TE。
一边改变相位编码倾斜磁场脉冲904的强度(相位编码量kp)和切片位置，一边以重复时间TR重复进行以上的步骤，测量与必要的切片张数对应的量的回波信号。根据RF脉冲902的频率，改变切片位置。
对每个切片，如图14(b)所示那样，将各回波信号918配置在k空间上，通过傅立叶逆变换重构图像。
在该SE时序900中，通过适当地改变TR和TE，得到强调了T1的对比度的图像(T1强调像)、强调了T2的对比度的图像(T2强调像)。例如，如果设TR为数百毫秒，TE为十毫秒左右，则得到T1强调像，如果设TR为数秒，TE为一百毫秒左右，则得到T2强调像。
解析地求出SE时序900的信号函数，用以下的公式(6)表示以SE时序 900摄影的图像的亮度值ISE。
[公式6]
ISE＝a(1-exp(-TR/T1))exp(-TE/T2)×sin(B1×90°)…(6)
在公式(6)中，使用推定出的被检测体参数T1、T2，并且指定摄影参数的TR和TE的值，由此能够以任意的强调度，自由地生成T1强调像和T2强调像。另外，通过将被检测体参数之一的B1设为均匀(例如设定为B1＝1)，能够得到B1分布均匀的情况下的图像。
另外，例如如果使用在上述各实施方式中推定出的被检测体参数T1、T2，将TR设定为500ms，将TE设定为15ms，则得到T1强调像。在通过SE时序900得到的T1强调像中，如果T1长，则亮度值小。例如，如果将TE直接设为15ms，将TR设为100ms那样的短，则得到强调了对比度、进而提高了 T1强调度的T1强调像。
另一方面，如果使TR和TE更长(例如(TR、TE)＝(4000ms、100ms))，则得到T2强调像。在T2强调像中，如果T2长，则亮度值大。
此外，成为对象的脉冲时序并不限于SE时序900。只要是梯度回波时序、 RSSG时序等如上述那样解析地得到信号函数fs的时序、或能够通过数值模拟生成的时序即可。
此外，在信号函数fs不是已知、或没有通过数值模拟生成的情况下，也可以使用所设定的摄影时序、摄影参数、被检测体参数、装置参数，通过数值模拟来生成图像。
这样，根据本变形例，能够使用推定出的推定参数的值、预先确定的摄影时序的信号函数，得到该摄影时序的预先确定的画质(任意的对比度)的图像。
此外，在上述各实施方式中，MRI装置100具备参数推定部230、显示图像生成部240、信号函数生成部250，但并不限于此。也可以将其中的至少一个结构构筑在例如能够与MRI装置100进行数据的收发的与MRI装置100独立的信息处理装置上。
附图标记说明
100：MRI装置；101：磁铁；102：倾斜磁场线圈；103：被检测体；104：时序发生器；105：倾斜磁场电源；106：高频磁场发生器；107：收发线圈； 108：接收器；109：计算机；110：显示装置；111：存储装置；210：接受部； 220：重构图像取得部；230：参数推定部；240显示图像生成部；250：信号函数生成部；310：摄影参数组；320：重构图像；321：低分辨率化后的图像； 330：信号函数；340：参数分布；341：参数分布；342：插补扩大后的参数分布；343：参数分布；344：现有方法的参数分布；410：摄影参数组；420：信号函数；513：B1分布的轮廓；514：B1分布的轮廓；523：a分布的轮廓；524： a分布的轮廓；610：RF-Spoiled GE时序；611：切片倾斜磁场脉冲；612：RF 脉冲；613：切片相位(兼切片相位编码)倾斜磁场脉冲；614：相位编码倾斜磁场脉冲；615：失相用读出倾斜磁场脉冲；616：读出倾斜磁场脉冲；617：回波信号；618：读出倾斜磁场脉冲；619：失相用相位编码倾斜磁场脉冲；620：失相用相位编码倾斜磁场脉冲；621：破碎脉冲；622：回波信号；630：多回波时序；711：多回波图像；712：多回波图像；713：单回波图像；714：单回波图像；715：单回波图像；716：单回波图像；721：多回波图像；722：多回波图像；740：参数分布；800：接受画面；810：推定参数接受区域；820：参数组接受区域；900：SE时序；901：切片选择倾斜磁场脉冲；902：RF脉冲； 903：切片相位倾斜磁场脉冲；904：相位编码倾斜磁场脉冲；905：失相用读出倾斜磁场脉冲；906：破碎脉冲；907：破碎脉冲；908：破碎脉冲；909：切片选择倾斜磁场脉冲；910：重聚焦脉冲；911：破碎脉冲；912：破碎脉冲； 913：破碎脉冲；914：读出倾斜磁场脉冲；915：A/D；916：破碎脉冲；918：回波信号。

Claims (14)

1.一种磁共振成像装置，其特征在于，具备：
参数推定部，其使用多个重构图像以及取得该重构图像时所依照的摄影时序的信号函数，来推定依存于被检测体的一个以上的被检测体参数以及依存于装置的一个以上的装置参数中的至少一个推定参数的值，其中，
上述多个重构图像是使用多个摄影参数组，依照上述摄影时序进行摄像，分别取得的重构图像，
上述信号函数是按照每个上述摄影时序确定摄影参数、上述被检测体参数以及上述装置参数与上述重构图像的各像素值之间的关系的函数，
上述参数推定部使用对每个上述推定参数最优的分辨率的上述重构图像来推定上述值。
2.根据权利要求1所述的磁共振成像装置，其特征在于，
上述多个重构图像分别具有每个上述推定参数的最优分辨率中的最高的分辨率即最高分辨率，
上述参数推定部针对上述最优分辨率比上述最高分辨率低的推定参数，将上述重构图像的分辨率降低为对该推定参数最优的分辨率，推定上述值。
3.根据权利要求2所述的磁共振成像装置，其特征在于，
上述参数推定部从上述最优分辨率低的上述推定参数开始顺序地推定上述值，并且在推定时，将已经推定出的推定参数的值调整为对推定对象的上述推定参数最优的分辨率，与上述信号函数一起使用。
4.根据权利要求1所述的磁共振成像装置，其特征在于，
上述推定参数包含接收线圈灵敏度、高频磁场强度以及共振频率差的至少一个。
5.根据权利要求1所述的磁共振成像装置，其特征在于，
上述推定参数包含纵缓和时间、横缓和时间以及质子密度的至少一个。
6.根据权利要求1所述的磁共振成像装置，其特征在于，
上述多个重构图像由具有第一分辨率的第一重构图像、具有与上述第一分辨率不同的第二分辨率的第二重构图像构成，其中，上述第一分辨率是上述推定参数中上述最优的分辨率最高的推定参数的该最优的分辨率。
7.根据权利要求6所述的磁共振成像装置，其特征在于，
上述第一重构图像是通过多回波时序取得的多回波图像，
上述参数推定部根据上述多回波图像，推定上述推定参数中的共振频率差的值。
8.根据权利要求7所述的磁共振成像装置，其特征在于，
上述参数推定部将上述多回波图像的分辨率调整为上述第二分辨率，将调整后的多回波图像和上述第二重构图像与上述信号函数一起使用，推定上述共振频率差以外的上述推定参数的值。
9.根据权利要求1所述的磁共振成像装置，其特征在于，
还具备：重构图像取得部，其取得上述多个重构图像。
10.根据权利要求1所述的磁共振成像装置，其特征在于，
对于上述多个摄影参数组，上述摄影参数中的重复时间、高频磁场的强度、高频磁场的相位的至少一个值分别不同。
11.根据权利要求1所述的磁共振成像装置，其特征在于，
上述重构图像的张数是上述推定参数个数以上。
12.根据权利要求1所述的磁共振成像装置，其特征在于，
还具备：显示图像生成部，其使用推定出的上述推定参数的值、预先确定的摄影时序的信号函数，得到该摄影时序的预先确定的画质的图像。
13.根据权利要求1所述的磁共振成像装置，其特征在于，
还具备：接受部，其接受上述推定参数的指示，
上述接受部经由推定参数接受画面，接受上述推定参数的设定，
重构图像取得部依照包含与接受的推定参数对应地预先保存的上述最优的分辨率的摄影条件，取得上述重构图像。
14.根据权利要求1所述的磁共振成像装置，其特征在于，
针对多个不同的上述摄影参数和上述被检测体参数的组合，分别执行数值模拟，对得到的结果进行插补，由此生成上述信号函数。
CN201580011025.4A 2015-08-06 2015-08-06 磁共振成像装置 Active CN107072586B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/072424 WO2017022136A1 (ja) 2015-08-06 2015-08-06 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
CN107072586A CN107072586A (zh) 2017-08-18
CN107072586B true CN107072586B (zh) 2020-08-18

Family

ID=57572457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580011025.4A Active CN107072586B (zh) 2015-08-06 2015-08-06 磁共振成像装置

Country Status (4)

Country Link
US (1) US11122988B2 (zh)
JP (1) JP6049954B1 (zh)
CN (1) CN107072586B (zh)
WO (1) WO2017022136A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11228758B2 (en) 2016-01-22 2022-01-18 Peking University Imaging method and device
EP4063628A1 (en) 2017-04-24 2022-09-28 General Electric Company Adaptive linear linked piston electric power generator
US10451697B2 (en) * 2017-07-25 2019-10-22 Spintech, Inc. Systems and methods for strategically acquired gradient echo imaging
JP6782681B2 (ja) * 2017-11-06 2020-11-11 株式会社日立製作所 磁気共鳴イメージング装置、撮像パラメータセット生成用演算装置及び撮像パラメータセット生成プログラム
EP3511725A1 (en) * 2018-01-11 2019-07-17 Koninklijke Philips N.V. Dual resolution dixon magnetic resonance imaging
JP7020930B2 (ja) * 2018-01-24 2022-02-16 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、磁気共鳴イメージングシステム及びパラメータ推定方法
JP7280997B2 (ja) * 2018-08-29 2023-05-24 富士フイルムヘルスケア株式会社 画像処理装置及び磁気共鳴イメージング装置
CN110495887B (zh) * 2019-07-30 2022-02-15 哈尔滨医科大学 一种针对磁共振双信号纳米探针使用rarevtr序列同时获得t1加权成像方法
JP7448418B2 (ja) * 2020-05-29 2024-03-12 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、及び、磁気共鳴イメージング方法
JP7377231B2 (ja) 2021-02-02 2023-11-09 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、及び、画像処理方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007020829A (ja) * 2005-07-15 2007-02-01 Hitachi Ltd 画像データ解析方法およびシステム
CN101658420A (zh) * 2008-08-26 2010-03-03 株式会社东芝 磁共振成像装置以及摄像条件设定方法
JP2011024926A (ja) * 2009-07-29 2011-02-10 Hitachi Medical Corp 磁気共鳴イメージング装置
CN103371819A (zh) * 2012-04-27 2013-10-30 株式会社东芝 磁共振成像装置以及图像处理装置
CN104323775A (zh) * 2013-07-22 2015-02-04 株式会社日立制作所 磁共振成像装置、图像处理装置以及图像处理方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL8801731A (nl) * 1988-07-08 1990-02-01 Philips Nv Werkwijze en kernspinresonantie-inrichting om de dwarsrelaxatietijdkonstante t2 snel te bepalen.
JP4924965B2 (ja) 2005-07-15 2012-04-25 マスセット株式会社 子供椅子用構造体
EP1991887B1 (en) * 2006-02-17 2018-10-17 Regents of the University of Minnesota High field magnetic resonance

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007020829A (ja) * 2005-07-15 2007-02-01 Hitachi Ltd 画像データ解析方法およびシステム
CN101658420A (zh) * 2008-08-26 2010-03-03 株式会社东芝 磁共振成像装置以及摄像条件设定方法
JP2011024926A (ja) * 2009-07-29 2011-02-10 Hitachi Medical Corp 磁気共鳴イメージング装置
CN103371819A (zh) * 2012-04-27 2013-10-30 株式会社东芝 磁共振成像装置以及图像处理装置
CN104323775A (zh) * 2013-07-22 2015-02-04 株式会社日立制作所 磁共振成像装置、图像处理装置以及图像处理方法

Also Published As

Publication number Publication date
CN107072586A (zh) 2017-08-18
US20170209067A1 (en) 2017-07-27
JP6049954B1 (ja) 2016-12-21
WO2017022136A1 (ja) 2017-02-09
JPWO2017022136A1 (ja) 2017-08-03
US11122988B2 (en) 2021-09-21

Similar Documents

Publication Publication Date Title
CN107072586B (zh) 磁共振成像装置
US8736265B2 (en) B1-mapping and B1L-shimming for MRI
US8754645B2 (en) Method for spatially resolved determination of an MR parameter
US9575154B2 (en) MR imaging using a multi-point dixon technique
US8710840B2 (en) Method and system for phase-sensitive magnetic resonance imaging
JP5449903B2 (ja) 磁気共鳴イメージング装置
US7701214B2 (en) Magnetic resonance imaging apparatus
JP6371554B2 (ja) 磁気共鳴イメージング装置
US10551467B2 (en) Method and apparatus for movement compensation during magnetic resonance imaging
US20150362576A1 (en) Metal resistant mr imaging
EP3044604B1 (en) Metal resistant mr imaging
EP4071494A1 (en) Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset
US11428769B2 (en) Magnetic resonance imaging device, calculation device for generation of imaging parameter set, and imaging parameter set generation program
US10866296B2 (en) Magnetic resonance imaging apparatus and calculation image generation method using pulse sequence
JP5684888B2 (ja) 磁気共鳴イメージング装置
CN113729684B (zh) 磁共振成像装置以及磁共振成像方法
JP6487554B2 (ja) 磁気共鳴イメージング装置
US10823803B2 (en) Method and apparatus for reconstructing contrast levels from magnetic resonance acquisitions
JP6114846B2 (ja) 磁気共鳴イメージング装置
JP5881869B2 (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211122

Address after: Chiba County, Japan

Patentee after: Fujifilm medical health Co.,Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi, Ltd.