[go: up one dir, main page]

CN103371819A - 磁共振成像装置以及图像处理装置 - Google Patents

磁共振成像装置以及图像处理装置 Download PDF

Info

Publication number
CN103371819A
CN103371819A CN2012105019838A CN201210501983A CN103371819A CN 103371819 A CN103371819 A CN 103371819A CN 2012105019838 A CN2012105019838 A CN 2012105019838A CN 201210501983 A CN201210501983 A CN 201210501983A CN 103371819 A CN103371819 A CN 103371819A
Authority
CN
China
Prior art keywords
interval
space data
magnetic resonance
space
channels
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
CN2012105019838A
Other languages
English (en)
Other versions
CN103371819B (zh
Inventor
竹岛秀则
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN103371819A publication Critical patent/CN103371819A/zh
Application granted granted Critical
Publication of CN103371819B publication Critical patent/CN103371819B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • 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/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • 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/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
    • 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/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/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • G01R33/5614—Generating steady state signals, e.g. low flip angle sequences [FLASH] using a fully balanced steady-state free precession [bSSFP] pulse sequence, e.g. trueFISP
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308—Characterization of motion or flow; Dynamic imaging
    • G01R33/56325—Cine imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5673—Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本实施方式涉及磁共振成像装置以及图像处理装置。提供一种能够提高画质的磁共振成像装置以及图像处理装置。实施方式所涉及的磁共振成像装置具备配置部、生成部以及重建部。上述配置部将通过并行成像收集到的多个通道量的磁共振信号在k空间的第1区域中以第1间隔进行配置，在比上述第1区域宽范围的第2区域中以比上述第1间隔大的第2间隔进行配置。上述生成部根据被配置于上述k空间的k空间数据中的第2间隔的k空间数据，生成多个通道量的第1间隔的k空间数据。上述重建部根据被配置于上述第1区域的第1间隔的k空间数据、根据第2间隔的k空间数据而生成的第1间隔的k空间数据、以及多个通道量的灵敏度分布，重建磁共振图像。

Description

磁共振成像装置以及图像处理装置
本申请主张2012年4月27日申请的日本专利申请号2012-103918及2012年9月13日申请的日本专利申请号2012-201752的优先权,并在本申请中引用上述日本专利申请的全部内容。
技术领域
实施方式涉及磁共振成像（imaging）装置以及图像处理装置。
背景技术
磁共振成像装置（以下，适当地称为“MRI（Magnetic ResonanceImaging）装置”）通过线圈（coil）来检测从被检体的原子（例如，氢原子）发出的核磁共振信号（以下，适当地称为“NMR（NuclearMagnetic Resonance）信号”）。另外，MRI装置将所检测到的NMR信号变换成被称为k空间数据（data）的原始数据，通过变换并重建该k空间数据，从而得到表示成为对象的原子的分布的图像（以下，适当地称为“MR图像”）。MRI装置对k空间数据应用傅立叶（Fourier）逆变换，而得到与通过线圈检测到的信号对应的MR图像。
k空间数据是从三维的对象物拍摄的数据。MRI装置为了收集二维的剖面图像或者三维的体（volume）图像所需要的数据，重复进行一维的拍摄。MRI装置通过被称为频率编码（encode）的方法来进行一维的拍摄。关于其他维的拍摄，MRI装置一边通过被称为相位编码的方法使相位发生变化一边重复进行拍摄。另外，当得到二维的剖面图像时，MRI装置通过对三维的对象物进行切片（slice）选择，来选择作为拍摄对象的剖面。
在此，例如，当收集与k空间中的所希望的分辨率“256×256”对应的k空间数据的情况下，MRI装置通过相位编码将进行拍摄的k空间线（line）的数量设定为“256”，因为重复进行256次一维的拍摄，所以拍摄时间变长。因此，以往，作为通过相位编码来减少进行拍摄的k空间线的数量的方法，讨论了并行成像（parallel imaging）技术。并行成像技术是着眼于在通过多个线圈同时收集k空间数据时，根据线圈的配置而灵敏度不同的技术。在并行成像技术中，存在被称为SENSE系统（sensitivity encoding：灵敏度编码）的技术、或被称为SMASH（simultaneous acquisition of spatial harmonics：空间谐波并行采集）系统的技术（例如，GRAPPA（generaized autocalibratingpartially parallel acquisition：全局自动校准部分并行采集)），被搭载在大量的MRI装置中而被广泛地使用。
发明内容
本发明要解决的问题在于提供一种能够提高画质的磁共振成像装置以及图像处理装置。
实施方式所涉及的磁共振成像装置具备配置部、生成部、以及重建部。上述配置部将通过并行成像收集到的多个通道量的磁共振信号在k空间的第1区域中以第1间隔进行配置，在比上述第1区域宽范围的第2区域中以比上述第1间隔大的第2间隔进行配置。上述生成部根据被配置在上述k空间的k空间数据中的第2间隔的k空间数据，生成多个通道量的第1间隔的k空间数据。上述重建部根据被配置在上述第1区域的第1间隔的k空间数据、基于第2间隔的k空间数据生成的第1间隔的k空间数据、以及多个通道量的灵敏度分布，重建磁共振图像。
根据实施方式的磁共振成像装置以及图像处理装置，能够提高画质。
附图说明
图1是表示第1实施方式所涉及的MRI装置的框（block）图。
图2是用于说明第1实施方式中的k空间数据的收集模式（pattern）的图。
图3是表示第1实施方式中的k空间数据的收集模式的图。
图4是表示第1实施方式中的k空间数据的收集模式的图。
图5是表示第1实施方式中的处理步骤的流程图（flowchart）。
图6是表示第1实施方式中的处理步骤的流程图。
图7是表示第1实施方式中的MR图像计算处理的图。
图8是表示第1实施方式的变形例中的处理步骤的流程图。
图9是表示第2实施方式中的MR图像计算处理的图。
图10是表示其他的实施方式中的k空间数据的采样（sampling）间隔的图。
具体实施方式
以下，详细地说明各实施方式所涉及的磁共振成像装置以及图像处理装置。另外，除非特别注明，以二维的k空间数据为例进行了说明，但即使k空间数据的维数不是2也能够利用以下所述的方法。另外，在各实施方式的说明中，假设在k空间的kx方向（频率编码方向）不进行间拔来进行说明，但该说明并不限制对在kx方向进行间拔的收集方法的应用。另外，在以下的说明中，示出了格子状地收集k空间的例子，但只要是能够利用并行成像技术的收集方法，也可以不限制使用格子状以外的扫描方法来收集k空间数据。另外，除非特别注明，设k空间数据以及MR图像通过复数来表现。在医疗诊断中所使用的MR图像例如通过对在以下的实施方式中说明的MR图像，针对各像素求复数的大小来得到。从k空间向图像空间的变换方法可以是傅立叶变换（Fourier transform）或傅立叶逆变换的某一个。从k空间向图像空间的变换方法，按照k空间数据的采用方法来决定。
（第1实施方式）
图1是表示第1实施方式所涉及的MRI装置100的框图。如图1所示，MRI装置100具备静磁场磁铁101、倾斜磁场线圈102、倾斜磁场电源103、床104、床控制部105、发送线圈106、发送部107、接收线圈阵列（coil array）108、接收部109、序列（sequence）控制部110、计算机系统120。另外，MRI装置100中不包含被检体P（例如，人体）。
静磁场磁铁101是形成为中空的圆筒形状的磁铁，在内部的空间中产生一致的静磁场。静磁场磁铁101例如是永久磁铁、超导磁铁等。倾斜磁场线圈102是形成为中空的圆筒状的线圈，被配置在静磁场磁铁101的内侧。倾斜磁场线圈102由与相互正交的X、Y、Z的各轴对应的3个线圈组合形成，这3个线圈从倾斜磁场电源103独立地接受电流，产生磁场强度沿着X、Y、Z的各轴发生变化的倾斜磁场。另外，Z轴方向与静磁场是同方向。
倾斜磁场电源103对倾斜磁场线圈102供给电流。在此，由倾斜磁场线圈102产生的X、Y、Z各轴的倾斜磁场例如分别与切片选择用倾斜磁场Gs、相位编码用倾斜磁场Ge、以及读出（readout）用倾斜磁场Gr对应。切片选择用倾斜磁场Gs用于任意地确定拍摄剖面。相位编码用倾斜磁场Ge用于根据空间位置使NMR信号的相位发生变化。读出用倾斜磁场Gr用于根据空间位置使NMR信号的频率发生变化。
床104具备载置被检体P的顶板104a，在基于床控制部105的控制下，将顶板104a以载置有被检体P的状态插入倾斜磁场线圈102的空洞（拍摄口）内。通常，床104被设置成长度方向与静磁场磁铁101的中心轴平行。床控制部105在基于计算机系统（system）120的控制下，驱动床104而使顶板104a向长度方向以及上下方向移动。
发送线圈106被配置在倾斜磁场线圈102的内侧，从发送部107接受RF脉冲的供给，产生高频磁场。发送部107对发送线圈106供给与拉莫尔（Larmor）频率对应的RF脉冲（pulse），该拉莫尔频率由作为对象的原子的种类以及磁场的强度决定。
接收线圈阵列108被配置在倾斜磁场线圈102的内侧，通过高频磁场的影响而接收从被检体P发出的NMR信号。如果接收线圈阵列108接收NMR信号，则将所接收到的NMR信号向接收部109输出。另外，在第1实施方式中，接收线圈阵列108是具有多个接收线圈的线圈阵列。
接收部109根据从接收线圈阵列108输出的NMR信号生成NMR数据。具体而言，接收部109通过对从接收线圈阵列108输出的NMR信号进行数字变换来生成NMR数据。另外，接收部109向序列控制部110发送所生成的NMR数据。另外，接收部109也可以被安装在具备静磁场磁铁101或倾斜磁场线圈102等的架台装置侧。在此，在第1实施方式中，从接收线圈阵列108的各线圈元件（coil element）（各接收线圈）输出的NMR信号通过适当地分配合成，以被称为通道（channel）等的单位输出至接收部109。因此，在接收部109以后的后级的处理中针对每个通道来处理NMR数据。线圈元件的总数与通道的总数的关系存在相同的情况，也存在通道的总数比线圈元件的总数少的情况或者相反地通道的总数比线圈元件的总数多的情况。在以下中，如记述为“每个接收线圈（每个通道）”的情况下，表示该处理可以针对每个线圈元件进行，或者也可以针对分配合成线圈元件后的每个通道进行。另外，分配合成的定时并不限定于上述的定时（timing）。在直到基于后述的图像重建部122的重建处理前，也可以按照通道单位来分配合成NMR信号或者NMR数据。
序列控制部110根据从计算机系统120发送的序列信息来驱动倾斜磁场电源103、发送部107以及接收部109，从而进行被检体P的拍摄。在此，序列信息是对用于进行拍摄的步骤进行定义的信息。在序列信息中，对倾斜磁场电源103向倾斜磁场线圈102供给的电源的强度或供给电源的定时、发送部107向发送线圈106发送的RF脉冲的强度或施加RF脉冲的定时、接收部109检测NMR信号的定时等进行定义。
另外，序列控制部110驱动倾斜磁场电源103、发送部107以及接收部109而对被检体P进行拍摄的结果是，如果从接收部109接收NMR数据，则将所接收到的NMR数据向计算机系统120传送。
计算机系统120进行MRI装置100的整体控制、或数据收集、图像重建等，具有接口（interface）部121、图像重建部122、存储部123、输入部124、显示部125以及控制部126。
接口部121向序列控制部110发送序列信息，从序列控制部110接收NMR数据。另外，如果接口部121接收NMR数据，则将接收到的NMR数据保存到存储部123中。保存于存储部123中的NMR数据由后述的收集部126a配置在k空间中。其结果，存储部123存储多个通道量的k空间数据。
图像重建部122通过对存储于存储部123的k空间数据实施傅立叶变换处理等重建处理，来生成频谱（spectrum）数据或图像数据。存储部123存储通过接口部121接收到的NMR数据、通过收集部126a配置在k空间中的k空间数据、以及由图像重建部122生成的图像数据等。例如，存储部123是RAM（Random Access Memory：随机存取存储器）、闪存器（flash memory）等半导体存储器元件、硬盘（harddisk）、光盘（disk）等。输入部124接受来自操作者的各种指示或信息输入。输入部124例如是鼠标（mouse）或轨迹球（track ball）等定位设备（pointing device）、模式（mode）切换开关（switch）等选择设备、或者键盘（keyboard）等输入设备。显示部125在基于控制部126的控制下，显示频谱数据或图像数据等各种信息。显示部125例如是液晶显示器等显示设备。
控制部126进行MRI装置100的整体控制。具体而言，控制部126根据经由输入部124从操作者输入的拍摄条件生成序列信息，并将所生成的序列信息发送至序列控制部110，从而进行拍摄。另外，控制部126控制根据作为拍摄的结果从序列控制部110发送的NMR数据而进行的图像的重建，或者控制显示部125进行的显示。例如，控制部126是ASIC（Application Specific Integrated Circuit：可编程专用集成电路）、FPGA（Field Programmable Gate Array：现场可编程门阵列）等集成电路、CPU（Central Processing Unit：中央处理单元）、MPU（Micro Processing Unit：微处理单元）等电子电路。
在此，第1实施方式所涉及的MRI装置100通过规则地使用基于第1间隔的间拔、和基于比第1间隔大的第2间隔的间拔来收集k空间数据，从而在实现缩短拍摄时间的同时提高被重建的MR图像的画质。该功能例如通过基于被安装在上述的控制部126内的各部的控制来实现。如图1所示，控制部126具备收集部126a和计算部126b。
收集部126a对与所希望的分辨率对应的k空间线，规则地使用基于第1间隔的间拔、和基于比第1间隔大的第2间隔的间拔，针对每个接收线圈（每个通道）收集k空间数据。例如，收集部126a控制序列控制部110等，针对每个接收线圈（每个通道）收集NMR数据，并将所收集到的NMR数据配置到k空间中，并保存到存储部123中。此时，收集部126a将所收集到的NMR数据在k空间的第1区域中以第1间隔进行配置，在比第1区域宽范围的第2区域中，以比第1间隔大的第2间隔来配置。基于收集部126a的处理的细节之后进行叙述。
计算部126b使用由收集部126a针对每个接收线圈（每个通道）集、并保存于存储部123的k空间数据计算MR图像。例如，计算部126b控制图像重建部122等，另外，使用存储部123所保存的k空间数据来计算MR图像。另外，例如，计算部126b进行控制，以使得计算出的MR图像显示在显示部125上。如图1所示，计算部126b具备k空间数据生成部126c和MR图像计算部126d。
k空间数据生成部126c使用由收集部126a以基于第2间隔的间拔针对每个接收线圈（每个通道）收集到的k空间数据，针对每个接收线圈（每个通道）生成作为填充以基于第2间隔的间拔而间拔了的k空间线的位置的k空间数据的第1间隔的k空间数据。另外，k空间数据生成部126c将针对每个接收线圈（每个通道）生成的第1间隔的k空间数据发送至MR图像计算部126d。另外，基于k空间数据生成部126c的处理的细节之后进行叙述。
MR图像计算部126d针对每个接收线圈（每个通道）组合通过k空间数据生成部126c针对每个接收线圈（每个通道）生成的第1间隔的k空间数据、与以基于第1间隔的间拔针对每个接收线圈（每个通道）收集到的k空间数据，计算MR图像。另外，基于MR图像计算部126d的处理细节之后进行叙述。
在此，对SENSE以及GRAPPA进行说明。一般而言，在SENSE的情况下，MRI装置为了预先训练（training）用而收集k空间数据，对所收集到的k空间数据重建，预先求出各接收线圈（各通道）的灵敏度分布（也被称为灵敏度图）。例如，作为事前拍摄（在对诊断图像进行拍摄的主拍摄之前进行的拍摄），MRI装置进行灵敏度图拍摄（也被称为校准（calibration）扫描（scan）等）来收集训练用的k空间数据（也被称为校准数据等），求出灵敏度分布。当进行训练后的拍摄时，MRI装置等间隔地间拔收集k空间数据。例如，MRI装置只收集第偶数个k空间线（line）。接着，如果MRI装置通过收集到的k空间数据对各通道（每个接收线圈的k空间数据）进行重建，则得到折叠（fold）了的MRI图像。之后，MRI装置利用各通道的MRI图像和通过训练求得的灵敏度分布来分离（unfolding）折叠分量，求出高分辨率的MRI图像。由于进行间拔采样，因此，能够缩短拍摄时间。例如，如果只进行第偶数个采样（sample），则其拍摄时间变为1/2。在SENSE中的unfold的过程中生成的线性系统根据未知的像素数与线圈数的关系，有时成为被称作ill-posed或者ill-conditioned的线性系统，此时，为了使线性系统稳定，也可以导入与解相关的已知知识进行归一化（regularize）处理。
另一方面，在GRAPPA的情况下，MRI装置等间隔地间拔收集k空间数据，在k空间上根据附近的数据来推定被间拔了的数据。该推定通过利用全部通道的k空间数据通过加权线性和来计算各通道的欠缺数据的方法来进行。该加权是在为了训练用而对k空间的一部分（例如，k空间的中心部分）进行了全（full）采样（不间拔地收集）的基础上，根据训练用的数据而计算出的。
众所周知，如果将k空间数据看作拍摄作为对象的原子分布的傅立叶变换而得到的数据，则相当于其低频区域的部分（称为“k空间的中心部分”）和相当于高频区域的部分（称为“k空间的周边部分”）相比较所收集到的值的大小变大。因为在SENSE或GRAPPA中（除了GRAPPA的训练用数据以外）等间隔地间拔，所以以同一密度来收集k空间的中心部分以及周边部分。但是，考虑通过使对MRI图像的画质做出更大贡献的k空间的中心部分的采样密度比周边部分高，从而提高相同的k空间线数下的重建画质。
也考虑对以可变密度随机地收集到的k空间数据应用并行成像技术。但是，例如，当生成所希望的采样密度函数、并按照该采样密度函数使采样密度发生变化时，如果是SENSE，则由unfolding进行处理的像素数增加。如果是GRAPPA，则推定式的数量爆发式地增加，另外，稀疏地收集到的部分的推定也需要大量的训练用线。
第1实施方式所涉及的MRI装置100想要通过对以基于“密集的”第1间隔的间拔而收集到的k空间数据应用并行成像技术来得到高画质的MR图像。当然，根据与拍摄时间的关系，以基于第1间隔的间拔来收集全部k空间线未必是优选的。因此，第1实施方式所涉及的MRI装置100针对k空间整体，一边以基于“稀疏的”第2间隔的间拔收集k空间数据，一边对k空间的中心部分以基于“密集的”第1间隔的间拔收集k空间数据，使用以基于该第2间隔的间拔收集到的k空间数据，生成填充被间拔了的k空间线的第1间隔的k空间数据。
在此，第1实施方式所涉及的MRI装置100在使用“稀疏的”第2间隔的k空间数据来生成“密集的”第1间隔的k空间数据时也应用并行成像技术。这样，第1实施方式所涉及的MRI装置100是通过多次应用并行成像技术来得到高画质的MR图像的装置，为了该目的，不是随机（random）地、而是规则地收集以不同的间隔间拔的数据。即，当至少以基于第2间隔的间拔来收集k空间数据时，MRI装置100通过规则地进行收集以使得遍及k空间整体而等间隔地收集数据，从而能够应用并行成像技术。另外，这也可以说是以能够利用既存的并行成像技术的形式对k空间数据的收集增加了制约。另外，在第1实施方式中，说明了当至少以基于第2间隔的间拔收集k空间数据时，遍及k空间整体而等间隔地收集数据的例子，但实施方式并不限定于此。例如，当以基于第2间隔的间拔收集k空间数据时，MRI装置100也可以对k空间的周边部分的一部分不收集k空间数据，而执行填零（zero fill）处理。
以下，详细地说明基于收集部126a的处理。图2是用于说明第1实施方式中的k空间数据的采样间隔的图，表示在等间隔的并行成像技术中所使用的收集模式（pattern）的一个例子。另外，在图2中，kx表示频率编码方向，ky表示相位编码方向。另外，在图2中，表示填充在kx方向的数据点组的箭头示出1k空间线。另外，实线的箭头表示未间拔而收集的k空间线，虚线的箭头表示通过间拔而未收集的k空间线。
图2的（A）所示的收集模式201是间拔（只收集1k空间线）2k空间线中的1k空间线的收集模式（以下，称为“2采样间隔”）。另外，图2的（B）所示的收集模式202是间拔（只收集1k空间线）4k空间线中的3k空间线的收集模式（以下，称为“4采样间隔”）。
在此说明的收集方法基于以下的考虑。以下，设A、B是比1大的实数。另外，以下，以A、B是整数时为例进行说明，但例如由于作为并行成像技术的SENSE（如果考虑被折叠的像素数根据位置而变化的情况）对于不是整数的间拔间隔也能够应用，因此，当使用SENSE来重建时A、B不需要是整数。
对于以长的间隔A的间拔得到的k空间数据，如果能够进行与间隔A对应的图像重建，则也能够根据重建后的MR图像，重新生成以更短的间隔B间拔的k空间。因此，分别关于在更短的间隔B中收集到的数据使用真值、关于未收集的数据使用重新生成的k空间数据，从而能够进行与间隔B对应的图像重建。因此，考虑以间隔A收集全部的k空间线，但以间隔B只收集一部分的k空间线。如果考虑收集的k空间线的数量，则能够认为通过由这样的收集方法（以下，称为第1实施方式的收集方法）得到的k空间数据的重建而获得的MR图像的画质比由对间隔A下的等间隔收集的重建而得到的MR图像的画质高，比由对间隔B下的等间隔收集的重建而得到的MR图像的画质低。考虑为了一边抑制收集的k空间线数的增加量，一边使使用第1实施方式的收集方法时的画质接近间隔B下的等间隔收集，优选确定收集模式，以使得大多数情况下在k空间中在信号强度高的中心部分中以间隔B、在信号强度低的周边部分中间隔A来进行收集。
图3是表示第1实施方式中的k空间数据的收集模式的图，表示适合第1实施方式的收集模式的一个例子。图3所示的收集模式301表示设A=4、B=2时的（对于二维k空间的）收集模式。设间隔A（=4）下的k空间线304为全部采样对象，而设间隔B（=2）下的k空间线302以及303中的、靠近k空间的中心部分的k空间线302为收集的对象，但没有将作为k空间的周边部分的k空间线303设为收集的对象。
作为另一个例子，图4是表示第1实施方式中的k空间数据的收集模式的图，表示对于三维k空间的收集模式的一个例子。另外，在图4中，只图示出ky以及kz方向（相位编码方向），省略了kx方向（频率编码方向）的图示。因此，各圆相当于“1”k空间线。另外，黑色圆表示不间拔而收集的k空间线，白色圆表示通过间拔而没有收集的k空间线。
图4所示的收集模式表示分别针对ky方向、kz方向，设A=4、B=2时的收集模式。该收集模式是通过将对图3的ky方向的采样位置应用于图4中的ky方向、kz方向的各个轴，且只将在ky方向以及kz方向的两轴中设为收集对象的位置作为收集对象，从而三维地扩张了的例子。例如，将位置401设为收集对象，但没有将位置402设为收集对象。另外，在图4中，示例出了在ky方向以及kz方向的两轴进行间拔的收集模式，但实施方式并不限定于此。例如，也可以是只在某一轴进行间拔、在另一轴上进行全采样的采样模式。例如，也可以是在ky方向进行间拔、但在kz方向进行全采样的采样模式。另外，例如，也可以是在kz方向进行间拔、但ky方向进行全采样的采样模式。
另外，上述示出了A=4、B=2的例子。即，第1实施方式所涉及的收集部126a以不间拔而收集k空间线的间隔在整数次中是1次（例如，在4次中是1次，在2次中是1次），且间隔A是间隔B的整数倍的方式，设定了间隔A以及间隔B，但例如也可以如A=6、B=2那样使用2以外的比率，也可以如A=3、B=1.5那样作为A或者B使用整数以外的值，也可以是如A=5、B=2那样A不是B的整数倍。其中，如果将A设定为B的整数倍，则在前级的并行成像中使用的k空间数据能够在全部后级的并行成像中利用，由于减少了利用以间隔B重新生成的包含噪声（noise）的推定采样值的次数，因此，与没有将A设定为B的整数倍时相比较能够期待提高画质。另外，根据所使用的图像重建方法（例如，如果是GRAPPA）而也可以不进行2阶段下的重建。针对使用了GRAPPA的例子之后进行叙述。
[使用了实施方式中的采样方法的整体的流程]
图5是表示第1实施方式中的处理步骤的流程图，表示在以限制采样位置的形式进行可变密度采样的情况下，由MRI装置100执行的图像重建的流程。
首先，收集部126a以基于间隔A的间拔来收集k空间整体，同时以基于间隔B的间拔来收集k空间的中心部分（步骤（step）S501）。换而言之，收集部126a收集包含间隔A的k空间线全部、和间隔B的k空间线的一部分的k空间线的集合（k空间线组）。该过程如之前所说明的那样，但在形式上可以说如以下那样。在步骤S501中，当将A以及B设为整数（A＞B≥2）时，收集满足以下的3个条件的k空间线。（1）包含全部的在k空间内以基于间隔A的间拔而收集到的k空间线即第1k空间线组（line set）。（2）包含一个以上的在k空间内以基于间隔B的间拔而收集到的k空间线即第2k空间线组中的、与第1k空间线组内的k空间线不同的k空间线。（3）使第2k空间线组内的k空间线的总数比间隔B下的全部k空间线数小。例如，收集部126a通过控制序列控制部110，从而能够进行上述的收集。
接着，计算部126b使用在步骤S501中收集到的k空间数据来计算MR图像（步骤S502）。例如，计算部126b通过控制图像重建部122，对存储于存储部123的k空间数据进行重建而得到MR图像。另外，设计算部126b通过某种方法已了解了通过收集部126a进行收集的位置。例如，计算部126b预先确定采样位置。另外，例如，设为在收集部126a对计算部126b发送k空间数据时，发送关于采样位置的信息。
接着，计算部126b进行控制，使得在步骤S502中重建的MR图像显示在显示部125上（步骤S503）。
图6是表示第1实施方式中的处理步骤的流程图，表示基于计算部126b的MR图像的重建方法的流程的一个例子。另外，作为说明图参照图7而对MR图像的重建方法进行说明。
首先，k空间数据生成部126c接受由收集部126a收集到的k空间数据（例如，图7所示的k空间数据701）的输入（步骤S601）。另外，当非既知k空间数据的采样位置的情况下，k空间数据生成部126c也可以在该步骤S601中进一步从收集部126a接受k空间数据的采样位置的信息输入。
接着，k空间数据生成部126c使用通过基于间隔A的间拔针对每个接收线圈（每个通道）收集到的k空间数据，来重建中间MR图像（步骤S602）。另外，在此生成的MR图像不是作为最终输出而使用的MR图像，因此，为了便于说明，称为中间MR图像。例如，如果使用SENSE（之后对SENSE中的求解的式子的一个例子进行叙述），则k空间数据生成部126c针对各通道（针对每个接收线圈（每个通道））重建接受了输入的k空间数据701，并在此基础上使用灵敏度分布703来执行与间隔A对应的展开。另外，在图7中，示出了以间拔前的分辨率重建了的折叠图像702，但也可以以间拔后的分辨率重建并执行SENSE。
如果对该步骤S602换而言之，则首先k空间数据生成部126c分别重建以基于间隔A的间拔针对每个接收线圈（每个通道）收集到的k空间数据701而分别得到折叠图像702。接着，k空间数据生成部126c针对折叠图像702的每一个分别应用各接收线圈（各通道）的灵敏度分布703进行并行成像技术的unfolding（展开），从而计算中间MR图像704。
接着，k空间数据生成部126c使用灵敏度分布703来逆重建在步骤S602中重建了的中间MR图像（步骤S603）。具体而言，通过以下步骤进行逆重建。中间MR图像是1个图像。首先，k空间数据生成部126c对该中间MR图像内的各像素应用各灵敏度分布703，生成加权了的每个通道（每个接收线圈）的中间MR图像。由于所生成的每个通道的中间MR图像是没有folding的图像（没有折叠的图像），因此，k空间数据生成部126c执行folding过程（按照依存于作为收集对象的位置的相位而确定的权重来进行加权的相加。权重由离散傅立叶变换的变换式可知）。
如果对该步骤S603换而言之，则首先k空间数据生成部126c对计算出的中间MR图像704分别应用各接收线圈（各通道）的灵敏度分布703，分别计算出反映各接收线圈（各通道）的灵敏度的中间MR图像704′，之后分别逆重建计算出的中间MR图像704′，针对每个接收线圈（每个通道）生成k空间数据。另外，k空间数据生成部126c通过对针对每个接收线圈（每个通道）生成的k空间数据分别进行间拔，从而针对每个接收线圈（每个通道）生成间隔B的k空间数据705（推定）。
并且，MR图像计算部126d利用在步骤S603中通过逆重建得到的推定k空间数据705、以及在步骤S601中作为输入而接受的k空间数据701，（针对每个通道）生成间隔B的k空间数据，进行间隔B下的重建，得到最终的MR图像709（步骤S604）。间隔B下的重建例如能够使用SENSE、使用灵敏度分布703以及在之前的操作中生成的k空间数据进行unfolding（展开），得到MR图像709这样的操作来执行。关于间隔B的k空间数据，例如对于在步骤S601中作为输入而提供的k空间线使用其k空间数据，对于在步骤S601的输入中欠缺但间隔B下的重建需要的k空间线，能够使用通过步骤S603的逆重建得到的k空间数据。另外，例如，针对在步骤S601中作为输入而提供的k空间数据，也可以使用与通过逆重建得到的k空间数据的加权平均。
如果对该步骤S604换而言之，则首先MR图像计算部126d组合由k空间数据生成部126c针对每个接收线圈（每个通道）生成的间隔B的k空间数据705、和通过收集部126a实际收集到的k空间数据701中的间隔B的k空间数据，分别生成每个接收线圈（每个通道）的间隔B的k空间数据706。另外，在图7中，与间隔B的k空间数据705进行组合的、实际收集到的间隔B的k空间数据是被椭圆707包围的中心部分的实线的k空间线。另外，组合可以是基于实际收集到的k空间数据的覆盖处理，或者也可以将两者加权组合。并且，MR图像计算部126d通过分别重建间隔B的k空间数据706而分别得到折叠图像708，分别对折叠图像708应用各接收线圈（各通道）的灵敏度分布703而进行并行成像技术的unfolding（展开），计算作为最终输出的MR图像709。
在以上的说明中，对将图像重建设为SENSE的2级结构来执行的情况进行了说明，但也可以是3级以上的结构。例如，如果是3级结构，则计算部126b如以下那样地执行重建。将C以及D设为满足C＞D的A的约数，计算部126b首先对A以及C（也就是将C作为之前叙述的B）进行2阶段的重建（步骤S601～604）。之后，计算部126b对C以及D（也就是将D作为之前叙述的B）执行使用了灵敏度分布的逆重建（步骤S603）以及重建（步骤S604）。另外，与之前说明的例子相同，C和D未必需要是A的约数。另外，针对3级以上的结构中的最后的1阶段，也可以是没有间拔的全采样。
[SENSE中的先验知识（Prior）的利用]
然而，作为通过高速拍摄而利用欠缺了k空间的一部分的k空间数据来重建MR图像的方法，知道有除了拍摄的k空间数据之外对先验知识（Prior）进行评价，并综合这些来进行MR图像的方法。先验知识也被称为已知知识。在此，当计算部126b使用已知知识来计算中间MR图像、并使用已知知识来计算最终输出的MR图像的情况下，中间MR图像的计算所应用的已知知识的强度（以下，适当地称为“先验知识强度”）也可以设定为比最终输出的MR图像的计算所应用的先验知识强度高。另外，关于先验知识，例如，参照“Tsao J.et al.,“Unifying Linear Prior-Information-Driven Methods for AcceleratedImage Acquisition,”Mag.Reson.In Med.46:652-660,2001.”（特别地，在p.653内的Table 1.Summary of Reconstruction Methods中，能够利用被称为Static refererence image的归一化项）。如果设先验知识强度为λ、folding后的输入图像为ρalias、推定MR图像为ρ、静态参照MR图像为ρstatic（也被称为“已知知识”、“既定图像”）、噪声协方差行列式为Ψ、S为灵敏度行列式、R为归一化所使用的行列式，则例如使用了归一化的SENSE能够以以下的形式来提供。
【式1】
ρ＝(S*Ψ-1S+λR*R)-1(S*Ψ-1ρa1ias+λR*Rρstatic)
针对SENSE的观点或其细节，记载在“Pruessmann K.P.etal.,“SENSE:Sensitivity Encoding for Fast MRI,”Mag.Reson.InMed.42:952-962,1999.”中。另外，如果λ=0，则变为不使用归一化的SENSE的式子。先验知识与SENSE例如能够通过对由SENSE求解的方程式追加先验知识项来组合利用。此时，作为参数（parameter）设定先验知识强度。当使用第1实施方式中的2级结构时，例如，设前级的先验知识强度比后级的先验知识强度强，从而能够在易于生成噪声的间隔A中，抑制SENSE的噪声。另外，例如也可以设在前级中λ＞0，在后级中λ=0（即也可以只在前级中使用先验知识）。
另外，在此所述的先验知识的利用方法并不特别地限定于SENSE，例如也可以是k-tSENSE。另外，当以时间序列对k空间数据进行拍摄时，也可以通过在进行收集的各采样位置在时间方向取平均来计算平均k空间数据，并在重建了该平均k空间数据的基础上，将重建后的MR图像作为静态参照MR图像来使用（参照“Tsao J.et al.，“k-tBLAST and k-t SENSE:Dynamic MRI With High Frame RateExploiting Spatiotemporal Correlations,”Mag.Reson.In Med.50:1031-1042,2003.”）。
[滤波器（filter）]
利用2级以上的结构时的优点在于，能够对前级的输出应用高画质化操作，例如，能够应用非线性噪声除去滤波器。图8是表示第1实施方式的变形例中的处理步骤的流程图。例如，k空间数据生成部126c能够在步骤S802的重建处理与步骤S804的逆重建处理之间，执行滤波器应用步骤S803。
作为在此利用的滤波器，例如，能够列举出ε滤波器。所谓ε滤波器是指一边将各像素切换为逐次关注的像素一边进行以下的操作的滤波器，也能够对复数数据利用。首先，k空间数据生成部126c针对在步骤S802中重建了的中间MR图像，提取以关注像素为中心的另行确定的范围的像素值。接着，如果关注像素与周边像素的像素值的差的大小为一定以上，则k空间数据生成部126c将周边像素的像素值看作另行确定的值（例如关注像素自身的值），否则，在将周边像素的像素值看作其自身的像素值的基础上，取周边像素的加权平均。如果使用ε滤波器，则不怎么增加处理量就能够一边保持轮廓一边选择性地减少噪声。另外，在滤波器应用步骤S803中应用的滤波器不一定是ε滤波器，也可以使用其他的滤波器。
如上述那样，根据第1实施方式，能够适当地收集k空间数据，进而能够提高MR图像的画质。
（第2实施方式）
在第1实施方式中，说明了作为并行成像技术使用了SENSE的方法，但实施方式并不限定于此。例如，作为并行成像技术，计算部126b也可以使用SMASH系统的GRAPPA。
[利用了GRAPPA的重建]
当利用GRAPPA进行MR图像的重建时，收集部126a除了与间隔A以及间隔B对应的采样位置的k空间数据之外，作为GRAPPA的训练用的数据还预先不间拔地收集k空间的中心部分的k空间数据。此时，k空间数据生成部126c能够根据间隔B的k空间数据，计算用于以间隔B插补间隔A的k空间数据的插补系数。并且，k空间数据生成部126c通过将该插补系数应用于间隔A的k空间数据，来推定被间拔的位置的间隔B的k空间数据。另外，k空间数据生成部126c能够根据训练用数据，计算用于插补间隔B的k空间数据的插补系数。并且，k空间数据生成部126c通过将该插补系数应用于间隔B的k空间数据，来推定其余的位置的k空间线。这样，k空间数据生成部126c能够推定k空间内的所有k空间线。另外，有时将不该间拔而收集到的训练用的数据称为“灵敏度分布”、“校准数据”等。
即，k空间数据生成部126c通过对以基于间隔A的间拔针对每个接收线圈（每个通道）收集到的k空间数据，应用由以基于间隔B的间拔针对每个接收线圈（每个通道）收集到的k空间数据导出的插补系数，而针对每个接收线圈（每个通道）推定间隔B的k空间数据。该间隔B的k空间数据（推定）是填充以基于间隔A的间拔而被间拔的k空间线的位置的k空间数据。另外，MR图像计算部126d通过对组合了推定出的间隔B的k空间数据、与实际上收集到的间隔B的k空间数据的每个接收线圈（每个通道）的k空间数据应用插补系数，来推定k空间内的所有k空间线，计算MR图像。
这样，当利用了GRAPPA时，具有以下效果，即，k空间数据生成部126c将间隔不是任意地而是例如限定为间隔A以及间隔B这2个，因此，仅仅通过进行限定数量（例如，2个）的插补系数集合的训练，就能够实现可变密度采样。
[基于GRAPPA以及SENSE的多级的重建]
另外，MRI装置100也可以设为通过基于GRAPPA以及SENSE的多级的重建来进行图像重建，并首先利用间隔大的GRAPPA，计算更小的间隔下的SENSE所需的k空间线，接着执行SENSE进行重建。
图9是表示第2实施方式中的MR图像计算处理的图。首先，与第1实施方式相同，k空间数据生成部126c接受由收集部126a收集到的k空间数据（例如，图9所示的k空间数据901）的输入。接着，k空间数据生成部126c对以基于间隔A的间拔针对每个接收线圈（每个通道）收集到的k空间数据，应用由以基于间隔B的间拔针对每个接收线圈（每个通道）收集到的k空间数据导出的插补系数，而针对每个接收线圈（每个通道）生成间隔B的k空间数据902（推定）。该间隔B的k空间数据902（推定）是填充以基于间隔A的间拔而间拔了的k空间线的位置的k空间数据。
之后，与第1实施方式相同，MR图像计算部126d组合由k空间数据生成部126c针对每个接收线圈（每个通道）生成的间隔B的k空间数据902（推定）、和由收集部126a实际收集到的k空间数据901中的间隔B的k空间数据，分别生成每个接收线圈（每个通道）的间隔B的k空间数据903。另外，在图9中，与间隔B的k空间数据902（推定）进行组合的、实际收集到的间隔B的k空间数据是被椭圆904包围的中心部分的实线的k空间线。并且，MR图像计算部126d分别重建间隔B的k空间数据903而分别得到折叠图像905，对各个折叠图像905分别应用各接收线圈（各通道）的灵敏度分布906进行并行成像技术的unfolding（展开），从而计算作为最终输出的MR图像907。
（其他的实施方式）
实施方式并不限定于上述的第1实施方式以及第2实施方式。
在上述的实施方式中，例如假设收集部126a规则地使用基于第1间隔的间拔、和基于比第1间隔大的第2间隔的间拔，针对每个接收线圈（每个通道）收集k空间数据的情况。然而，实施方式并不限定于此。即，实施方式并不限定于如何收集k空间数据。只要存在“在第1区域以第1间隔进行配置，在比第1区域大的范围的第2区域中，以比第1间隔大的第2间隔来进行配置的k空间数据”即可。即，在上述的实施方式以及后述的实施方式中，以“收集”为前提进行了说明，但实施方式并不限定于如何收集，只要存在以规定的配置模式配置在k空间中的k空间数据即可。
此时，例如，图3或图4等能够解释为表示k空间数据的配置模式（也被称为采样模式等）的一个例子的图。例如，图3所示的配置模式301表示设第1间隔=2，第2间隔=4时的（对于二维k空间的）配置模式。将第2间隔（=4）下的k空间线304全部设为配置的对象，将第1间隔下的k空间线302以及303中的靠近k空间的中心部分的k空间线302设为配置的对象，但没有将作为k空间的周边部分的k空间线303设为配置的对象。
例如，在第1实施方式中说明了的例子的情况下，收集部126a将通过并行成像收集到的多个通道量的NMR数据在存储部123的k空间的第1区域中以第1间隔进行配置，在比第1区域大的范围的第2区域中以比第1间隔大的第2间隔进行配置。例如，收集部126a将NRM数据在与k空间的中心部分对应的第1区域中以第1间隔进行配置，在包含k空间的周边部分的宽范围的第2区域中以第2间隔进行配置。k空间数据生成部126c利用第2间隔的k空间数据重建多个通道量的第1折叠图像组（例如，图7的702），并对第1折叠图像组应用多个通道量的灵敏度分布（例如，图7的703），从而重建中间磁共振图像（例如，图7的704），通过根据多个通道量的灵敏度分布（例如，图7的703）逆重建中间磁共振图像，从而生成多个通道量的第1间隔的k空间数据（例如，图7的705）。MR图像计算部126d根据配置在第1区域中的第1间隔的k空间数据、以及根据第2间隔的k空间数据生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组（例如，图7的708），通过对第2折叠图像组应用多个通道量的灵敏度分布（例如，图7的703），来重建磁共振图像（例如，图7的709）。
另外，例如，在第2实施方式中说明的例子（利用了GRAPPA的重建）的情况下，收集部126a将通过并行成像收集到的多个通道量的NMR数据在存储部123的k空间的第1区域中以第1间隔进行配置，在比第1区域宽范围的第2区域中以比第1间隔大的第2间隔进行配置。例如，收集部126a将NRM数据在与k空间的中心部分对应的第1区域中以第1间隔进行配置，在包含k空间的周边部分的宽范围的第2区域中以第2间隔进行配置。k空间数据生成部126c通过对第2间隔的k空间数据应用根据第1间隔的k空间数据导出的插补系数，来生成多个通道量的第1间隔的k空间数据。MR图像计算部126d对在第1区域中配置的第1间隔的k空间数据、以及根据第2间隔的k空间数据生成的第1间隔的k空间数据中的至少一部分的k空间数据，应用根据灵敏度分布导出的插补系数，来重建磁共振图像。
另外，例如，在第2实施方式中说明了的例子（基于GRAPPA以及SENSE的多级的重建）的情况下，收集部126a将通过并行成像收集到的多个通道量的NMR数据在存储部123的k空间的第1区域中以第1间隔进行配置，在比第1区域宽范围的第2区域中，以比第1间隔大大第2间隔进行配置。例如，收集部126a将NRM数据在与k空间的中心部分对应的第1区域中以第1间隔进行配置，在包含k空间的周边部分的宽范围的第2区域中以第2间隔进行配置。k空间数据生成部126c通过对第2间隔的k空间数据应用根据第1间隔的k空间数据导出的插补系数，来生成多个通道量的第1间隔的k空间数据（例如，图9的902）。MR图像计算部126d根据被配置在第1区域中的第1间隔的k空间数据、以及由第2间隔的k空间数据生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组（例如，图9的905），通过对第2折叠图像组应用多个通道量的灵敏度分布（例如，图9的906）来重建磁共振图像（例如，图9的907）。
另外，在k空间的第1区域中以第1间隔进行配置、在第2区域中以第2间隔进行配置的定时未必需要是相同的定时。例如，在某一阶段的k空间中，只要在其第2区域中以第2间隔配置k空间数据即可。换而言之，在该阶段中，也可以不在k空间的第1区域中配置第1间隔的k空间数据。例如，在该阶段中k空间数据生成部126c能够执行根据该第2间隔的k空间数据生成第1间隔的k空间数据的处理。在第1区域中以第1间隔来配置k空间数据的处理既可以在基于k空间数据生成部126c的处理的前后进行，也可以并行进行。当并行进行的情况下，能够缩短从配置到MR图像重建的整体的处理时间。另外，实施方式包含以第1间隔来配置k空间数据的k空间、和以第2间隔来配置k空间数据的k空间不是同一k空间的情况。
[时间序列数据的采样]
在上述的实施方式中，以收集二维或者三维的k空间数据而不是时间序列的情况为例进行了说明，但实施方式并不限定于此。实施方式所涉及的MRI装置100即使在收集按照时间序列排列的规定帧（frame）数的k空间数据（以下，适当地称为“时间序列数据”）的情况下，也能够执行上述的各种处理。例如，当收集部126a收集时间序列数据的情况下，在相位编码方向规则地使用基于间隔A的间拔和基于间隔B的间拔来收集k空间数据，并且在时间方向规则地使用基于间隔A的间拔和基于间隔B的间拔来收集k空间数据。另外，在收集部126a收集时间序列数据的情况下，针对每帧变更不间拔而被收集的k空间线的位置。
图10是表示其他的实施方式中的k空间数据的采样间隔的图，表示例如k-tSENSE中使用的那样的采样位置的例子。另外，在图10中，只图示出ky方向（相位编码方向）以及时间方向，省略了kx方向（频率编码方向）的图示。因此，各圆相当于“1”k空间线。另外，黑色圆表示不间拔而收集的采样位置1001。另外，虚线的框1002～1005分别表示在时刻t=1～4时收集到的k空间数据。关于k-tSENSE，参照“Tsao J.et al.,“k-t BLAST and k-t SENSE:Dynamic MRI With HighFrame Rate Exploiting Spatiotemporal Correlations,”Mag.Reson.In Med.50:1031-1042,2003.”。
在图10所示的例子的情况下，收集部126a对各时刻t，设间隔A=4（4采样间隔）、间隔B=2（2采样间隔）来确定相位编码方向的采样位置。另外，收集部126a对以4采样间隔进行收集的部分（k空间的周边部分）以周期4（以4时刻单位）、针对以2采样间隔进行收集的部分（k空间的中心部分）以周期2（以2时刻单位）来依次使采样位置发生变化。即，即使对时间方向，收集部126a也规则地使用基于间隔A的间拔和基于间隔B的间拔来收集k空间数据。这样，由于收集部126a能够针对采样间隔A以及采样间隔B的双方确保时间方向的周期性的采样位置，因此能够原样地使用k-tSENSE。
另外，k空间数据的充填方法未必需要如图10所示的那样地从左到右的顺序，例如，也可以是从右到左的顺序。更一般而言，作为将各时刻的k空间数据称为帧的例子，在使按照时间序列排列的连续的A帧量的采样位置重叠了的情况下，收集部126a只要将采样位置确定为与k空间内的全部采样位置一致即可。这样，收集部126a能够一边按照时间序列切换相位编码方向的采样位置，一边收集k空间数据。
[其他的结构]
另外，在上述的实施方式中，说明了作为医用图像诊断装置的MRI装置100执行各种处理的情况，但实施方式并不限定于此。例如，图像处理装置或包含MRI装置100与图像处理装置的图像处理系统也可以执行上述的各种处理，代替MRI装置100。在此，所谓图像处理装置，例如，是工作站、PACS（Picture Archiving and CommunicationSystem：图片存档和通信系统）的图像保管装置（图像服务器）或浏览器（viewer）、电子病历系统的各种装置等。此时，例如，图像处理装置通过从MRI装置100、或者从图像服务器经由网络接收、或者通过经由记录介质由操作者输入等而接受由MRI装置100收集到的k空间数据并保存到存储部中。并且，图像处理装置也可以将存储于存储部的该k空间数据作为对象，执行上述的各种处理（例如，基于计算部126b的处理）。另外，在上述的实施方式中，说明了将MR图像显示在显示部125的例子，但实施方式并不限定于此，MRI装置100也可以在计算了MR图像之后不进行显示，而是例如对存储部123输出图像数据等，并结束处理。此时，例如，通过MRI装置100或图像处理装置计算出的MR图像由其他的装置等利用。
[程序（program）]
在上述的实施方式中示出的处理步骤所示的指示能够根据作为软件（software）的程序来执行。通用的计算机系统通过预先存储该程序，并读入该程序，从而也能够得到与基于上述的实施方式的MRI装置或图像处理装置的效果相同的效果。在上述的实施方式中叙述的指示作为能够使计算机执行的程序，记录在磁盘（软（flexible）盘、硬盘等）、光盘（CD-ROM、CD-R、CD-RW、DVD-ROM、DVD±R、DVD±RW等）、半导体存储器、或者与这些类似的记录介质中。如果是计算机或者固件（firmware）系统可读的存储介质，其存储形式也可以是任何方式。如果计算机从该记录介质中读入程序，并根据该程序在CPU中执行程序所叙述的指示，则能够实现与上述的实施方式的MRI装置或图像处理装置相同的动作。当然，计算机取得或者读入程序时，也可以通过网络来取得或者读入。
另外，根据从存储介质安装到（install）计算机或编入系统的程序的指示在计算机上运行的OS（操作（operating）系统）、或数据库（database）管理软件、网络（network）等MW（中间软件（middleware））等也可以执行用于实现上述的实施方式的各处理的一部分。
另外，存储介质并不限定于与计算机或者编入系统独立的介质，还包含下载（download）并存储或者暂时存储通过LAN（Local AreaNetwork：局域网）或因特网（internet）等来传送的程序的存储介质。
另外，存储介质并不限定为1个，在由多个介质执行上述的实施方式中的处理时也包含在实施方式中的存储介质中，介质的结构也可以是任何结构。
另外，实施方式中的计算机或编入系统用于根据存储于存储介质的程序，来执行上述的实施方式中的各处理，也可以是由个人电脑、微型计算机（micro computer）等之一构成的装置、多个装置网络连接而成的系统等任意的结构。
另外，实施方式中的计算机并不限定于个人电脑，还包含信息处理设备所包含的运算处理装置、微型计算机等，是能够通过程序实现实施方式中的功能的设备、装置的总称。
根据以上所述的至少一个实施方式的磁共振成像装置、图像处理装置、系统以及方法，能够适当地进行k空间数据的收集。
虽然说明了本发明的几个实施方式,但这些实施方式是作为例子而提示的,并不意图限定本发明的范围。这些实施方式能够以其他的各种方式进行实施，在不脱离发明的要旨的范围内，能够进行各种的省略、置换、变更。这些实施方式或其变形与包含于发明的范围或要旨中，并且包含于权利要求书记载的发明及其均等的范围中。

Claims (16)

1.一种磁共振成像装置，其特征在于，具备：
配置部，将通过并行成像收集到的多个通道量的磁共振信号在k空间的第1区域中以第1间隔进行配置，在比上述第1区域宽范围的第2区域中以比上述第1间隔大的第2间隔进行配置；
生成部，根据被配置在上述k空间中的k空间数据中的第2间隔的k空间数据，生成多个通道量的第1间隔的k空间数据；以及
重建部，根据被配置于上述第1区域的第1间隔的k空间数据、基于第2间隔的k空间数据而生成的第1间隔的k空间数据、以及多个通道量的灵敏度分布，重建磁共振图像。
2.根据权利要求1所述的磁共振成像装置，其特征在于，
上述配置部将上述磁共振信号在与k空间的中心部分对应的上述第1区域中以上述第1间隔进行配置，在包含k空间的周边部分的宽范围的上述第2区域中以上述第2间隔进行配置。
3.根据权利要求1所述的磁共振成像装置，其特征在于，
上述生成部根据上述第2间隔的k空间数据重建中间磁共振图像，并根据多个通道量的灵敏度分布逆重建上述中间磁共振图像，从而生成上述多个通道量的第1间隔的k空间数据。
4.根据权利要求2所述的磁共振成像装置，其特征在于，
上述生成部根据上述第2间隔的k空间数据重建中间磁共振图像，并根据多个通道量的灵敏度分布逆重建上述中间磁共振图像，从而生成上述多个通道量的第1间隔的k空间数据。
5.根据权利要求1所述的磁共振成像装置，其特征在于，
上述生成部根据上述第2间隔的k空间数据重建多个通道量的第1折叠图像组，并对上述第1折叠图像组应用多个通道量的灵敏度分布，从而重建中间磁共振图像，通过根据多个通道量的灵敏度分布逆重建上述中间磁共振图像，从而生成上述多个通道量的第1间隔的k空间数据，
上述重建部根据被配置于上述第1区域的第1间隔的k空间数据、以及基于第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组，对上述第2折叠图像组应用多个通道量的灵敏度分布，从而重建上述磁共振图像。
6.根据权利要求2所述的磁共振成像装置，其特征在于，
上述生成部根据上述第2间隔的k空间数据重建多个通道量的第1折叠图像组，并对上述第1折叠图像组应用多个通道量的灵敏度分布，从而重建中间磁共振图像，通过根据多个通道量的灵敏度分布逆重建上述中间磁共振图像，生成上述多个通道量的第1间隔的k空间数据，
上述重建部根据被配置于上述第1区域的第1间隔的k空间数据、以及基于第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组，通过对上述第2折叠图像组应用多个通道量的灵敏度分布，重建上述磁共振图像。
7.根据权利要求5所述的磁共振成像装置，其特征在于，
在通过上述生成部使用已知知识来重建上述中间磁共振图像，并且通过上述重建部使用已知知识来重建上述磁共振图像的情况下，上述中间磁共振图像的重建所应用的已知知识的强度被设定为比上述磁共振图像的重建所应用的已知知识的强度高。
8.根据权利要求6所述的磁共振成像装置，其特征在于，
在通过上述生成部使用已知知识来重建上述中间磁共振图像，并且通过上述重建部使用已知知识来重建上述磁共振图像的情况下，上述中间磁共振图像的重建所应用的已知知识的强度被设定为比上述磁共振图像的重建所应用的已知知识的强度高。
9.根据权利要求1所述的磁共振成像装置，其特征在于，
上述生成部对上述第2间隔的k空间数据应用根据上述第1间隔的k空间数据导出的插补系数，从而生成上述多个通道量的第1间隔的k空间数据，
上述重建部对被配置于上述第1区域的第1间隔的k空间数据、以及根据第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，应用根据上述灵敏度分布而导出的插补系数，从而重建磁共振图像。
10.根据权利要求2所述的磁共振成像装置，其特征在于，
上述生成部对上述第2间隔的k空间数据应用根据上述第1间隔的k空间数据导出的插补系数，从而生成上述多个通道量的第1间隔的k空间数据，
上述重建部对被配置于上述第1区域的第1间隔的k空间数据、以及基于第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，应用根据上述灵敏度分布导出的插补系数，从而重建磁共振图像。
11.根据权利要求1所述的磁共振成像装置，其特征在于，
上述生成部对上述第2间隔的k空间数据应用根据上述第1间隔的k空间数据导出的插补系数，从而生成上述多个通道量的第1间隔的k空间数据，
上述重建部利用被配置于上述第1区域的第1间隔的k空间数据、以及根据第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组，并且对上述第2折叠图像组应用多个通道量的灵敏度分布，从而重建上述磁共振图像。
12.根据权利要求2所述的磁共振成像装置，其特征在于，
上述生成部对上述第2间隔的k空间数据应用根据上述第1间隔的k空间数据导出的插补系数，从而生成上述多个通道量的第1间隔的k空间数据，
上述重建部利用被配置于上述第1区域的第1间隔的k空间数据、以及根据第2间隔的k空间数据而生成的第1间隔的k空间数据中的至少一部分的k空间数据，重建多个通道量的第2折叠图像组，并且对上述第2折叠图像组应用多个通道量的灵敏度分布，从而重建上述磁共振图像。
13.根据权利要求1所述的磁共振成像装置，其特征在于，
上述配置部以上述第2间隔为上述第1间隔的整数倍的方式，将上述磁共振信号配置在k空间中。
14.根据权利要求1所述的磁共振成像装置，其特征在于，
上述配置部将沿着时间序列排列的规定帧数量的磁共振信号在时间方向也以上述第1间隔和上述第2间隔进行配置。
15.根据权利要求1所述的磁共振成像装置，其特征在于，
在将按照时间序列排列的规定帧数量的磁共振信号配置到k空间时，上述配置部将配置k空间数据的位置、和不配置k空间数据的位置配置成针对每帧而不同。
16.一种图像处理装置，其特征在于，具备：
存储部，存储将通过并行成像收集到的多个通道量的磁共振信号在k空间的第1区域中以第1间隔进行配置，在比上述第1区域宽范围的第2区域中以比上述第1间隔大的第2间隔进行配置的k空间数据；
生成部，根据被配置于上述k空间的k空间数据中的第2间隔的k空间数据，生成多个通道量的第1间隔的k空间数据；以及
重建部，根据被配置于上述第1区域的第1间隔的k空间数据、基于第2间隔的k空间数据而生成的第1间隔的k空间数据、以及多个通道量的灵敏度分布，重建磁共振图像。
CN201210501983.8A 2012-04-27 2012-11-30 磁共振成像装置以及图像处理装置 Active CN103371819B (zh)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2012103918 2012-04-27
JP2012-103918 2012-04-27
JP2012201752A JP5984258B2 (ja) 2012-04-27 2012-09-13 磁気共鳴イメージング装置及び画像処理装置
JP2012-201752 2012-09-13

Publications (2)

Publication Number Publication Date
CN103371819A true CN103371819A (zh) 2013-10-30
CN103371819B CN103371819B (zh) 2016-06-08

Family

ID=49458106

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210501983.8A Active CN103371819B (zh) 2012-04-27 2012-11-30 磁共振成像装置以及图像处理装置

Country Status (3)

Country Link
US (2) US9188655B2 (zh)
JP (1) JP5984258B2 (zh)
CN (1) CN103371819B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106028927A (zh) * 2013-12-18 2016-10-12 三星电子株式会社 磁共振成像设备及操作磁共振成像设备的方法
CN106659418A (zh) * 2014-07-29 2017-05-10 株式会社日立制作所 磁共振成像装置及图像重构方法
CN107072586A (zh) * 2015-08-06 2017-08-18 株式会社日立制作所 磁共振成像装置
CN110114685A (zh) * 2016-12-22 2019-08-09 皇家飞利浦有限公司 使用最大似然重建的对心脏的t1映射
CN111257809A (zh) * 2020-01-23 2020-06-09 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备
CN112782629A (zh) * 2020-12-30 2021-05-11 深圳市联影高端医疗装备创新研究院 一种磁共振扫描控制方法、存储介质及系统
CN117607770A (zh) * 2024-01-22 2024-02-27 中国科学技术大学先进技术研究院 磁共振图像重建方法、装置、电子设备及存储介质

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102011086369B4 (de) * 2011-11-15 2022-08-11 Siemens Healthcare Gmbh Verfahren zur Erstellung von MR-Angiographiebildern und entsprechende Magnetresonanzanlage
JP5984258B2 (ja) * 2012-04-27 2016-09-06 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置
JP6073627B2 (ja) 2012-10-01 2017-02-01 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置
JP6482842B2 (ja) 2013-12-27 2019-03-13 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6494986B2 (ja) 2014-01-16 2019-04-03 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置
US10330760B2 (en) * 2014-03-21 2019-06-25 Beth Israel Deaconess Medical Center, Inc. System and method for assessing T2-relaxation times with improved accuracy
KR101630762B1 (ko) 2014-07-30 2016-06-15 삼성전자주식회사 자기 공명 영상 생성 장치 및 방법
WO2016069602A1 (en) * 2014-10-31 2016-05-06 The Trustees Of The University Of Pennsylvania A method and device for magnetic resonance imaging data acquisition guided by physiologic feedback
KR102008498B1 (ko) * 2016-05-19 2019-08-07 삼성전자주식회사 자기 공명 영상 장치 및 그 동작방법
KR101836235B1 (ko) * 2016-05-27 2018-03-08 한국과학기술원 자기공명영상 생성 방법 및 장치
US10732246B2 (en) 2016-08-30 2020-08-04 Canon Medical Systems Corporation Magnetic resonance imaging apparatus
US11255942B2 (en) 2016-08-30 2022-02-22 Canon Medical Systems Corporation Magnetic resonance imaging apparatus
WO2018085969A1 (zh) * 2016-11-08 2018-05-17 中国科学院深圳先进技术研究院 磁共振血管壁成像方法和设备
GB2565786B (en) * 2017-08-22 2020-01-15 Siemens Healthcare Ltd A method of performing magnetic resonance imaging and a magnetic resonance apparatus
US10670676B2 (en) 2018-01-31 2020-06-02 Canon Medical Systems Corporation Image reconstructing method and reconstructing apparatus
JP7164320B2 (ja) * 2018-05-11 2022-11-01 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置、医用画像処理装置、及び画像再構成方法
JP7355532B2 (ja) * 2018-06-29 2023-10-03 キヤノンメディカルシステムズ株式会社 医用情報処理装置、方法及びプログラム
US11493585B2 (en) * 2018-06-29 2022-11-08 Canon Medical Systems Corporation Medical information processing apparatus and medical information processing method
US10527699B1 (en) * 2018-08-01 2020-01-07 The Board Of Trustees Of The Leland Stanford Junior University Unsupervised deep learning for multi-channel MRI model estimation
JP7206073B2 (ja) * 2018-08-21 2023-01-17 キヤノンメディカルシステムズ株式会社 医用画像処理装置、医用画像生成装置、医用画像処理方法、およびプログラム
JP6988001B2 (ja) * 2018-08-30 2022-01-05 オリンパス株式会社 記録装置、画像観察装置、観察システム、観察システムの制御方法、及び観察システムの作動プログラム
US10768252B2 (en) * 2018-10-24 2020-09-08 General Electric Company Methods and systems for sampling k-space data in magnetic resonance imaging
CN112494029B (zh) * 2019-11-29 2025-04-11 上海联影智能医疗科技有限公司 实时mr电影数据重建方法和系统
CN113075603B (zh) * 2021-03-25 2022-07-29 上海交通大学 磁共振成像脉冲序列设计方法
JP7746039B2 (ja) * 2021-06-10 2025-09-30 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び方法
JP7749394B2 (ja) 2021-09-29 2025-10-06 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置、画像生成方法、および画像生成プログラム
DE102022214456A1 (de) * 2022-12-30 2024-07-11 Siemens Healthineers Ag Korrektur von Bewegungseffekten in mittels einer Magnetresonanzanlage aufgenommenen Magnetresonanz(MR)-Daten

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070255130A1 (en) * 2006-04-13 2007-11-01 Jiang Du Magnetic Resonance Imaging with High Spatial and Temporal Resolution
CN101627910A (zh) * 2008-07-17 2010-01-20 株式会社东芝 磁共振成像装置及磁共振成像方法
US20100053677A1 (en) * 2008-08-31 2010-03-04 Brother Kogyo Kabushiki Kaisha Printing system, printer and terminal device therefor, and recording medium containing program for printer or terminal device

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6801800B2 (en) 1999-11-29 2004-10-05 Kabushiki Kaisha Toshiba MR imaging using ECG-prep scan
JP3952247B2 (ja) * 1999-12-08 2007-08-01 株式会社日立メディコ 核磁気共鳴撮影装置
CN100396240C (zh) * 2005-02-28 2008-06-25 西门子（中国）有限公司 频域灵敏度编码磁共振并行成像方法
US9700220B2 (en) 2006-04-25 2017-07-11 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5366370B2 (ja) * 2006-09-06 2013-12-11 株式会社東芝 磁気共鳴イメージング装置
US7382127B2 (en) * 2006-09-15 2008-06-03 General Electric Company System and method of accelerated MR propeller imaging
WO2008111416A1 (ja) * 2007-03-09 2008-09-18 Hitachi Medical Corporation 磁気共鳴イメージング装置および磁気共鳴イメージング方法
KR100902518B1 (ko) 2007-06-15 2009-06-15 한국과학기술원 일반급수 병렬영상법을 이용한 고해상도 자기공명영상생성방법 및 그 기록매체
US7688068B2 (en) 2008-05-06 2010-03-30 General Electric Company System and method for using parallel imaging with compressed sensing
US8736265B2 (en) * 2008-09-17 2014-05-27 Koninklijke Philips N.V. B1-mapping and B1L-shimming for MRI
KR101030676B1 (ko) 2009-01-05 2011-04-22 한국과학기술원 높은 시공간 해상도의 기능적 자기 공명 영상을 위한 고차 일반 급수 병렬 영상법 및 샘플링 법
US8581582B2 (en) 2009-09-18 2013-11-12 Kabushiki Kaisha Toshiba MRI non-contrast time-slip angiography using variably positioned cine sub-sequence
US8543187B2 (en) 2010-10-13 2013-09-24 Kabushiki Kaisha Toshiba Non-contrast MRI signal acquisition
DE102012205811B4 (de) * 2012-04-10 2014-12-11 Siemens Aktiengesellschaft Bestimmen von Abtastungspunkten eines zufälligen Unterabtastungsschemas in der MR-Bildgebung
JP5984258B2 (ja) * 2012-04-27 2016-09-06 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置
JP6073627B2 (ja) * 2012-10-01 2017-02-01 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像処理装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070255130A1 (en) * 2006-04-13 2007-11-01 Jiang Du Magnetic Resonance Imaging with High Spatial and Temporal Resolution
CN101627910A (zh) * 2008-07-17 2010-01-20 株式会社东芝 磁共振成像装置及磁共振成像方法
US20100053677A1 (en) * 2008-08-31 2010-03-04 Brother Kogyo Kabushiki Kaisha Printing system, printer and terminal device therefor, and recording medium containing program for printer or terminal device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
翁卓 等: "基于k空间加速采集的磁共振成像技术", 《中国生物医学工程学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106028927A (zh) * 2013-12-18 2016-10-12 三星电子株式会社 磁共振成像设备及操作磁共振成像设备的方法
CN106028927B (zh) * 2013-12-18 2019-04-05 三星电子株式会社 磁共振成像设备及操作磁共振成像设备的方法
CN106659418A (zh) * 2014-07-29 2017-05-10 株式会社日立制作所 磁共振成像装置及图像重构方法
CN107072586A (zh) * 2015-08-06 2017-08-18 株式会社日立制作所 磁共振成像装置
CN107072586B (zh) * 2015-08-06 2020-08-18 株式会社日立制作所 磁共振成像装置
CN110114685A (zh) * 2016-12-22 2019-08-09 皇家飞利浦有限公司 使用最大似然重建的对心脏的t1映射
CN110114685B (zh) * 2016-12-22 2022-03-08 皇家飞利浦有限公司 使用最大似然重建的对心脏的t1映射
CN111257809A (zh) * 2020-01-23 2020-06-09 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备
CN112782629A (zh) * 2020-12-30 2021-05-11 深圳市联影高端医疗装备创新研究院 一种磁共振扫描控制方法、存储介质及系统
CN112782629B (zh) * 2020-12-30 2022-05-03 深圳市联影高端医疗装备创新研究院 一种磁共振扫描控制方法、存储介质及系统
CN117607770A (zh) * 2024-01-22 2024-02-27 中国科学技术大学先进技术研究院 磁共振图像重建方法、装置、电子设备及存储介质
CN117607770B (zh) * 2024-01-22 2024-04-16 中国科学技术大学先进技术研究院 磁共振图像重建方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
US20130285655A1 (en) 2013-10-31
US9213075B2 (en) 2015-12-15
US9188655B2 (en) 2015-11-17
JP2013240571A (ja) 2013-12-05
CN103371819B (zh) 2016-06-08
JP5984258B2 (ja) 2016-09-06
US20130285662A1 (en) 2013-10-31

Similar Documents

Publication Publication Date Title
CN103371819B (zh) 磁共振成像装置以及图像处理装置
JP6073627B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
US10302725B2 (en) MRI apparatus that continuously generates and collects multiple echo K-space data from which odd and/or even numbered collected groups of K-space data are extracted and separately processed
US9396562B2 (en) MRI reconstruction with incoherent sampling and redundant haar wavelets
JP6400375B2 (ja) 画像処理装置及び磁気共鳴イメージング装置
JP4152381B2 (ja) 磁気共鳴イメージング装置
US10281549B2 (en) Magnetic resonance imaging apparatus and image processing apparatus
JP6494986B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JPWO2008111416A1 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
CN105103001A (zh) 使用高snr同相图像和较低snr至少部分地异相图像的dixon型水/脂肪分离mri
US11016154B2 (en) Magnetic resonance imaging apparatus and image reconstruction program
JP6548531B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
US9535148B2 (en) Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution
JP7140551B2 (ja) 磁気共鳴イメージング装置、処理装置、および医用画像処理方法
JP2017136170A (ja) 磁気共鳴イメージング装置
JP4675936B2 (ja) 核磁気共鳴撮影装置
JP2017100014A (ja) 磁気共鳴イメージング装置及び画像処理装置
JP4906952B2 (ja) 核磁気共鳴撮影装置
JP2004528099A (ja) Ｓｍａｓｈを用いた合成平均化
JP2007050278A (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
C14 Grant of patent or utility model
C41 Transfer of patent application or patent right or utility model
GR01 Patent grant
TA01 Transfer of patent application right

Effective date of registration: 20160517

Address after: Japan Tochigi

Applicant after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Applicant before: Toshiba Corp

Applicant before: Toshiba Medical System Co., Ltd.