JP2019188205A - 磁気共鳴イメージング装置 - Google Patents

磁気共鳴イメージング装置 Download PDF

Info

Publication number
JP2019188205A
JP2019188205A JP2019124132A JP2019124132A JP2019188205A JP 2019188205 A JP2019188205 A JP 2019188205A JP 2019124132 A JP2019124132 A JP 2019124132A JP 2019124132 A JP2019124132 A JP 2019124132A JP 2019188205 A JP2019188205 A JP 2019188205A
Authority
JP
Japan
Prior art keywords
imaging
sequence
unit
data
image
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
JP2019124132A
Other languages
English (en)
Other versions
JP6786670B2 (ja
Inventor
徳典 木村
Tokunori Kimura
徳典 木村
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
Canon 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 Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Publication of JP2019188205A publication Critical patent/JP2019188205A/ja
Application granted granted Critical
Publication of JP6786670B2 publication Critical patent/JP6786670B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/563Image 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/5635Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/70Means for positioning the patient in relation to the detecting, measuring or recording means
    • A61B5/704Tables
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • A61B5/7425Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
    • 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/5608Data 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
    • 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/561Image 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/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • 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/561Image 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/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • G01R33/5614Generating steady state signals, e.g. low flip angle sequences [FLASH] using a fully balanced steady-state free precession [bSSFP] pulse sequence, e.g. trueFISP
    • 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/563Image 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/56308Characterization of motion or flow; Dynamic imaging
    • 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/563Image 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/56366Perfusion imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • 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/5602Image 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
    • 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/5607Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood imaging
    • 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/563Image 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/56308Characterization of motion or flow; Dynamic imaging
    • G01R33/56333Involving spatial modulation of the magnetization within an imaged region, e.g. spatial modulation of magnetization [SPAMM] tagging

Landscapes

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

Abstract

【課題】撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる磁気共鳴イメージング装置を提供する。【解決手段】実施形態に係る磁気共鳴イメージング装置は、分割部と、収集部と、合成部とを備える。分割部は、被検体の撮像領域に流入する流体を標識化するためのRF波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間に基づいて、前記複数のタイミングを少なくとも2つの時間的な範囲に分割する。収集部は、前記撮像領域の全体について、第1の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が第1の読み出しシーケンスと異なる第2の読み出しシーケンスそれぞれを用いてデータ収集を行う。合成部は、第1の読み出しシーケンスによって収集されたデータから生成された画像と、第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。【選択図】図1

Description

本発明の実施形態は、磁気共鳴イメージング装置に関する。
従来、磁気共鳴イメージング装置で用いられるパルスシーケンスに含まれる読み出しシーケンスの種類には、大きく分類すると、SE(Spin Echo)系とFE(Field Echo)(GRE(Gradient Echo))系とがあり、さらに、FE系には、FLASH(Fast Low-Angle shot)系とSSFP(Steady State Free Precession)系とがある。一般的に、これらの読み出しシーケンスは、前者ほど、単位収集時間あたりのSN比(SNR:Signal-to-Noise Ratio)は低いが磁場不均一性(歪みや信号抜け)に強く、後者ほど、単位収集時間あたりのSN比は高いが磁場不均一性に弱いことが知られている。このため、磁気共鳴イメージング装置の操作者は、各読み出しシーケンスの特性を把握したうえで、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類や撮像条件を適宜に使い分けている。
Edelmann RR et al、Radiology 192:513-519(1994) 木村徳典: Modified STAR using asymmetric inversion slabs(ASTAR)法による非侵襲血流イメージング、日磁医誌 2001;20(8)、374-385 Kwong KK, Chesler DA, koff RM, Donahue KM, et al、MR perfusion studies with T1-weighted echo planar imaging、Magn Reson Med 1995; 34:878-887 Mani S et al. MRM 37:898-905 (1997) Wells JA at al. In vivo hadamard encoded continuous arterial spin labeling(H−CASL). MRM 63:1111-1118 (2010)
本発明が解決しようとする課題は、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる磁気共鳴イメージング装置を提供することである。
実施形態に係る磁気共鳴イメージング(Magnetic Resonance Imaging:MRI)装置は、分割部と、収集部と、合成部とを備える。分割部は、被検体の撮像領域に流入する流体を標識化するためのRF波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間に基づいて、前記複数のタイミングを少なくとも2つの時間的な範囲に分割する。収集部は、前記撮像領域の全体について、第1の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第1の読み出しシーケンスと異なる第2の読み出しシーケンスそれぞれを用いてデータ収集を行う。合成部は、前記第1の読み出しシーケンスによって収集されたデータから生成された画像と、前記第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。
図1は、第1の実施形態に係るMRI装置の構成例を示す図である。 図2は、従来のASL法の一例を説明するための図(1)である。 図3は、従来のASL法の一例を説明するための図(2)である。 図4は、従来のASL法の一例を説明するための図(3)である。 図5は、従来の3DTOF−MRAを説明するための図(1)である。 図6は、従来の3DTOF−MRAを説明するための図(2)である。 図7は、第1の実施形態に係るFE系の読み出しシーケンスの一例を示す図である。 図8は、第1の実施形態に係るSSFP系の読み出しシーケンスの一例を示す図である。 図9は、第1の実施形態に係るMRI装置の詳細な構成例を示す図である。 図10は、従来のST−MI−ASLによるデータ収集を示す図である。 図11は、第1の実施形態に係る収集部によって行われるST−MI−ASLを示す図である。 図12は、第1の実施形態に係る合成部によって行われる画像補正を説明するための図である。 図13は、第1の実施形態に係るMRI装置によって行われる処理の流れを示すフローチャートである。 図14は、第2の実施形態に係るMRI装置によって行われるマルチスラブ法を用いたASL−MRAを説明するための図である。 図15は、第2の実施形態に係るMRI装置によって行われるST−MI−ASLを用いたASL−MRAを説明するための図である。 図16は、第3の実施形態に係るMRI装置によって行われるPASL法を説明するための図(1)である。 図17は、第3の実施形態に係るMRI装置によって行われるPASL法を説明するための図(2)である。 図18は、第3の実施形態に係るMRI装置によって行われるPASL法を説明するための図(3)である。 図19は、第3の実施形態に係るMRI装置によって行われるマルチスラブ法を用いたPASL法を説明するための図である。 図20は、第4の実施形態に係るMRI装置によって行われるHE−MT−PASLを説明するための図である。 図21は、第4の実施形態に係るMRI装置によって行われるST−MR−ASLを説明するための図(1)である。 図22は、第4の実施形態に係るMRI装置によって行われるST−MR−ASLを説明するための図(2)である。 図23は、第5の実施形態に係るMRI装置によって行われるマルチスラブ3DTOF−MRAを説明するための図である。 図24は、第6の実施形態に係るMRI装置によって行われるマルチスラブASLを説明するための図である。
以下では添付図面を参照して、MRI装置の実施形態を詳細に説明する。
(第1の実施形態)
まず、第1の実施形態について説明する。図1は、第1の実施形態に係るMRI装置の構成例を示す図である。図1に示すように、MRI装置100は、静磁場磁石1、傾斜磁場コイル2、傾斜磁場電源3、寝台4、寝台制御部5、送信RFコイル6、送信部7、受信RFコイル8、受信部9、シーケンサ10、ECG(Electrocardiogram)センサ21、ECGユニット22、及び計算機システム30を備える。
静磁場磁石1は、中空の円筒形状に形成された磁石であり、内部の空間に一様な静磁場を発生する。この静磁場磁石1としては、例えば永久磁石、超伝導磁石等が使用される。
傾斜磁場コイル2は、中空の円筒形状に形成されたコイルであり、静磁場磁石1の内側に配置される。この傾斜磁場コイル2は、互いに直交するX,Y,Zの各軸に対応する3つのコイルが組み合わされて形成されており、これら3つのコイルは、後述する傾斜磁場電源3から個別に電流供給を受けて、X,Y,Zの各軸に沿って磁場強度が変化する傾斜磁場を発生させる。なお、Z軸方向は、静磁場と同方向とする。傾斜磁場電源3は、傾斜磁場コイル2に電流を供給する。
ここで、傾斜磁場コイル2によって発生するX,Y,Z各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Ge及びリードアウト用傾斜磁場Grにそれぞれ対応する。スライス選択用傾斜磁場Gsは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Geは、空間的位置に応じて磁気共鳴信号の位相を変化させるために利用される。リードアウト用傾斜磁場Grは、空間的位置に応じて磁気共鳴信号の周波数を変化させるために利用される。
寝台4は、被検体Pが載置される天板4aを備え、後述する寝台制御部5による制御のもと、被検体Pが載置された状態で天板4aを傾斜磁場コイル2の空洞(撮像口)内へ挿入する。通常、この寝台4は、長手方向が静磁場磁石1の中心軸と平行になるように設置される。寝台制御部5は、制御部36による制御のもと、寝台4を制御する装置であり、寝台4を駆動して、天板4aを長手方向及び上下方向へ移動する。
送信RFコイル6は、傾斜磁場コイル2の内側に配置され、送信部7からRF(Radio Frequency)パルスの供給を受けて、被検体PにRF波を印加する。送信部7は、ラーモア周波数に対応するRF波を送信RFコイル6に送信する。
受信RFコイル8は、傾斜磁場コイル2の内側に配置され、RF波の影響によって被検体Pから放射される磁気共鳴信号を受信する。この受信RFコイル8は、磁気共鳴信号を受信すると、その磁気共鳴信号を受信部9へ出力する。
受信部9は、受信RFコイル8から出力される磁気共鳴信号に基づいてMR(Magnetic Resonance:磁気共鳴)データを生成する。具体的には、この受信部9は、受信RFコイル8から出力される磁気共鳴信号をデジタル変換することによってMRデータを生成する。このMRデータは、前述したスライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Ge及びリードアウト用傾斜磁場Grによって、PE(Phase Encode)方向、RO(Read Out)方向、SE(Slice Encode)方向の空間周波数の情報が対応付けられることで、k空間に対応するデータとして生成される。そして、MRデータを生成すると、受信部9は、そのMRデータをシーケンサ10へ送信する。なお、受信部9は、静磁場磁石1や傾斜磁場コイル2などを備える架台装置側に備えられていてもよい。
シーケンサ10は、計算機システム30から送信されるシーケンス実行情報にしたがって傾斜磁場電源3、送信部7及び受信部9を駆動することによって、被検体Pの画像を生成するためのデータ収集を行う。ここで、シーケンス実行情報とは、被検体PからMRデータを収集するためのパルスシーケンスを定義する情報である。具体的には、シーケンス実行情報は、傾斜磁場電源3が傾斜磁場コイル2に供給する電源の強さや電源を供給するタイミング、送信部7が送信RFコイル6に送信するRF信号の強さやRF信号を送信するタイミング、受信部9が磁気共鳴信号を検出するタイミングなど、データ収集を行うための手順を定義した情報である。
なお、シーケンサ10は、傾斜磁場電源3、送信部7及び受信部9を駆動して被検体Pをスキャンした結果、受信部9からMRデータが送信されると、そのMRデータを計算機システム30へ転送する。
ECGセンサ21は、被検体Pの体表に付着され、被検体Pの心拍、脈波、呼吸などを示すECG信号を電気信号として検出する。ECGユニット22は、ECGセンサ21により検出されたECG信号にA/D変換処理やディレー処理を含む各種処理を施してゲート信号を生成し、生成したゲート信号をシーケンサ10に送信する。
計算機システム30は、MRI装置100の全体制御を行う。例えば、計算機システム30は、上述した各部を駆動することで、データ収集や画像再構成などを行う。この計算機システム30は、インタフェース部31、画像再構成部32、記憶部33、入力部34、表示部35、及び制御部36を有する。
インタフェース部31は、計算機システム30とシーケンサ10との間でやり取りされる各種信号の送受信を制御する。例えば、このインタフェース部31は、シーケンサ10に対してシーケンス実行情報を送信し、シーケンサ10からMRデータを受信する。MRデータを受信すると、インタフェース部31は、各MRデータを被検体Pごとに記憶部33に格納する。
画像再構成部32は、記憶部33によって記憶されたMRデータに対して後処理すなわちフーリエ変換等の再構成処理を施すことで、被検体Pの体内が描出された画像データを生成する。
記憶部33は、インタフェース部31により受信されたMRデータや、画像再構成部32により生成された画像データなどを被検体Pごとに記憶する。
入力部34は、操作者からの各種指示や情報入力を受け付ける。この入力部34としては、マウスやトラックボールなどのポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスを適宜に利用可能である。
表示部35は、制御部36による制御のもと、スペクトラムデータあるいは画像データ等の各種の情報を表示する。この表示部35としては、液晶表示器などの表示デバイスを利用可能である。
制御部36は、図示していないCPU(Central Processing Unit)やメモリ等を有し、MRI装置100の全体制御を行う。具体的には、この制御部36は、入力部34を介して操作者から受け付けられた各種指示に基づいてシーケンス実行情報を生成し、生成したシーケンス実行情報をシーケンサ10に送信することによってスキャンを制御したり、スキャンの結果としてシーケンサ10から送られるMRデータに基づいて行われる画像の再構成を制御したりする。
以上、第1の実施形態に係るMRI装置100の構成例について説明した。このような構成のもと、MRI装置100は、各種の読み出しシーケンスを用いて、様々な撮像法で被検体を撮像することが可能である。なお、ここでは、励起用のRFパルスを印加してから磁気共鳴信号を取得するまでのシーケンスを「読み出しシーケンス」と呼び、読み出しシーケンスの前に印加される標識化用のRFパルスなどの各種プリパルスと読み出しシーケンスとを含んだ全体的なシーケンスを「パルスシーケンス」と呼ぶ。
ここで、従来、MRI装置で用いられるパルスシーケンスに含まれる読み出しシーケンスの種類には、大きく分類すると、SE系とFE(GRE)系とがあり、さらに、FE系には、FLASH系とSSFP系とがある。一般的に、これらの読み出しシーケンスは、前者ほど、単位収集時間あたりのSN比は低いが磁場不均一性(歪みや信号抜け)に強く、後者ほど、単位収集時間あたりのSN比は高いが磁場均一性に弱いことが知られている。このため、磁気共鳴イメージング装置の操作者は、各読み出しシーケンスの特性を把握したうえで、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類や撮像条件を適宜に使い分けている。
例えば、MRI装置による撮像法の一つに、ASL(Arterial Spin Labeling)法がある。ASL法は、撮像領域に流入する流体を標識化するためのRF波を印加してから所定の待ち時間が経過した後にデータ収集を行うことによって、非造影で血流を画像化する撮像法である。
図2〜4は、従来のASL法を説明するための図である。一般的に、ASL法では、背景信号を抑制するために、タグ画像と呼ばれる画像と、コントロール画像と呼ばれる画像とがそれぞれ生成され、両画像の差分画像が生成される。
例えば、図2に示すように、ASL法では、撮像領域(imaging slab)の上流に設定されたタグ領域(tag slab)に、撮像領域に流入する血液(arterial blood flow)を標識化するためのRF波が印加される。そして、図3に示すように、RF波(tag)が印加されてから所定の待ち時間TIが経過した後にデータ収集(imaging)が行われ、収集されたデータからタグ画像が生成される。さらに、ASL法では、図2に示すように、撮像領域を挟んでタグ領域の反対側に設定されたコントロール領域(control slab)にもRF波が印加される。そして、図3に示すように、RF波(control)が印加されてから、タグ画像用のデータ収集が行われたときと同じ待ち時間TIが経過した後にデータ収集(imaging)が行われ、収集されたデータからコントロール画像が生成される。その後、図4に示すように、生成されたタグ画像(tag image)とコントロール画像(control image)との差分画像を生成することで、背景信号を抑制したASL画像(ASL image)が得られる。
また、ASL法は、その用途から、大きく、ASL−MRA(MR Angio)とASL−MRP(MR Perfusion)とに分けられる。ASL−MRAは、血管画像を得るための方法であり、ASL−MRPは、灌流画像を得るための方法である。さらに、近年では、ASL−MRA又はASL−MRPの時間的な変化(hemodynamics)を観察するために、標識化用のRF波を印加してからデータ収集までの待ち時間を動的に変えて、複数回撮像を行うことも行われている。
ここで、ASL−MRAとASL−MRPとは、どちらも撮像法の種類はASL法であるが、技術的なポイントが異なる。例えば、ASL−MRAでは、血管のSN比が大きいため加算回数は少なくてよいが、時間分解能及び空間分解能が要求される。また、ASL−MRPでは、ASL−MRAに比べて血流信号のSN比が背景の静止組織よりもはるかに低い(1/1000以下)画像が得られるため、タグ画像とコントロール画像との差分を行う場合に、それぞれの画像でMTC(Magnetization Transfer Contrast)効果を同じにしたり、低マトリクスでの複数回の加算平均が行われたりする。さらに、ASL−MRPの中でも、血流量は腎や肺などは高く、脳や前立腺、筋などは低いため、撮像対象の臓器によって読み出しシーケンスの撮像条件が異なる。
一方、一般的に最も広く利用されている3DTOF(3 Dimensional Time-Of-Flight)−MRAでは、血管流のタイムオブフライト(time-of-flight)効果を用いるため、1回の撮像当たりの撮像領域(imaging slab)の厚さは薄いほどSN比が高くなる。そのため、3DTOF−MRAでは、撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行うマルチスラブ(multi-slab)法が用いられる。
図5及び6は、従来の3DTOF−MRAを説明するための図である。図5は、シングルスラブ(single-slab)法による撮像の一例を示しており、図6は、マルチスラブ(multi-slab)法による撮像の一例を示している。
例えば、図5に示すように、シングルスラブ法では、撮像領域(imaging slab)に照射するRF波に傾斜をつける(流入部から距離的に遠いほどRF波のフリップ角(Flip Angle)を浅くする)TONE(Tilted Optimized Non-saturation Excitation)法(ISCE(Inclined Slab for Contrast Enhancement)法)が併用される場合が多い。また、図6に示すように、マルチスラブ法でTONE法が併用される場合もある。このような3D−MRAで頭部を撮像する場合は、ウィリス動脈輪(Willis ring)部の近傍にある主幹動脈部の抜け(void)が少ない画像が得られるFE(GRE)系の読み出しシーケンスが用いられるのが一般的である。
また、MRI装置による他の撮像法の一つに、拡散イメージング(diffusion imaging)がある。この拡散イメージングでは、シングルショットのSE−EPI(Echo Planar Imaging)系の読み出しシーケンスが主に用いられるが、歪みの大きな部位についてはFSE(Fast Spin Echo)系の読み出しシーケンスも用いられている。また、EPIの中でも、TE(Echo Time)、収集マトリクス(matrix)、SENSE(Sensitivity Encoding)−Factor、ショット数が使い分けられている。さらには、k空間の軌跡(trajectory)の違いであるPROPELLER(Periodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction)、radial系などのサンプリング(sampling)法も使い分けられている。
さらに、最近の傾向として、T1やT2、プロトン(proton)密度、ADC(Apparent Diffusion Coefficient)、血流量、ケミカルシフト、磁化率、温度などのMRIに関する各種パラメータを強調画像ではなく定量値画像(パラメータ画像)に変換する手法が盛んに用いられている。これらの各種パラメータの定量値画像を得るための撮像法では、様々な読み出しシーケンスが利用可能になっており、SN比や歪み、時間などの特性によって、使用する読み出しシーケンスの種類や撮像条件が決められるようになっている。
このように、MRI装置では、各種の撮像法で被検体を撮像することが可能であり、各撮像法において、撮像対象の部位や撮像の目的に応じて読み出しシーケンスの種類や撮像が使い分けられている。
例えば、前述したASL法では、エコー信号を収集するための読み出しシーケンス(imaging sequence)として、FE(GRE)、SSFP、FE−EPIなどのFE(GRE)系の読み出しシーケンスが用いられ、さらには、FSE、SE−EPIなどのSE系の読み出しシーケンスも用いられる。また、k空間を充填する方法も、Cartesian、non−Cartesian(spiral、radial)などの様々な方法が用いられる。これらの中には、対象部位(対象臓器)や使用目的の全ての場合に最適な方法、すなわちオールマイティな方法はなく、MRAかMRPかに応じて、又は、撮像対象の臓器に応じて、適宜使い分けられている。
ここで、前述したように、FE(GRE)系はSE系と比べてSN比は高いが磁場不均一性に弱く、さらにFE系の中でも、SSFP(Fast Imaging with Steady-state Precession:FISPとも呼ばれる)系はFLASH系と比べてSN比は高いが、さらに磁場不均一性に弱い。そのため、ASL−MRAで脳の主幹部を撮像する場合は、副鼻腔や側頭骨などがあるためFE系が主に用いられ、腹部の臓器を撮像する場合は、SSFP系が用いられることが多い。また、脳の中でも、頭頂部の末梢動脈を撮像する場合は磁場不均一性がさほど問題にならないため、SN比が高いSSFP系が優れていることが知られている。
また、Time−Resolve MRA(trMRA)では、標識化された血液は、TIが短いときは主幹動脈にあり、TIが長くなると末梢動脈に移行していく。そのため、TIが短いときはFE(GRE)系が、TIが長いときはSSFP系が適していることになる。また、3DTOF−MRAでは、ウィリス動脈輪(Willis ring)部の近傍はFE系で、頭頂部はSSFP系が適していることになる。
このように、MRI装置による撮像法では、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、それぞれの位置に適した読み出しシーケンスの種類や撮像条件を使い分けるのが望ましい。しかしながら、一般的に、操作者が時間的な位置や空間的な位置に応じて読み出しシーケンスの種類や撮像条件を切り替えることは難しく、装置側で制御する仕組みが求められていた。
上述した従来技術の課題に対し、第1の実施形態に係るMRI装置100は、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも2つの範囲に分割し、第1の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、第2の範囲についてはシーケンスの種類又は撮像条件が第1の読み出しシーケンスと異なる第2の読み出しシーケンスを用いてデータ収集を行う。そして、MRI装置100は、第1の読み出しシーケンスによって収集されたデータから生成された画像と第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。
このような構成によれば、第1の読み出しシーケンスの種類又は撮像条件を第1の範囲を撮像するために適したものとし、第2の読み出しシーケンスの種類又は撮像条件を第2の範囲を撮像するために適したものとすることで、被検体を撮像する際の時間的な位置又は空間的な位置に応じて、適切な読み出しシーケンスを用いて撮像が行われることになる。すなわち、第1の実施形態に係るMRI装置100によれば、被検体を撮像する際の時間的な位置又は空間的な位置に応じて、各位置に適した読み出しシーケンスの種類又は撮像条件を用いて撮像を行うことで、より高画質な画像を容易に得ることができる。
以下では、第1の実施形態に係るMRI装置100の詳細な構成、及び、第1の実施形態に係るMRI装置100によって行われる処理の流れについて説明する。なお、本実施形態では、被検体の撮像領域に流入する流体を標識化するための1回のRF波を印加してから複数の待ち時間が経過するごとに連続的にデータ収集を行うST−MI(Single-Tag Multiple-TI)−ASLで被検体の頭部血管を撮像する場合の例を説明する。また、本実施形態では、第1の読み出しシーケンスとして、一般的なFE系(以下、単にFE系と呼ぶ)の読み出しシーケンスを用い、第2の読み出しシーケンスとして、SSFP系の読み出しシーケンスを用いる場合の例を説明する。
図7は、第1の実施形態に係るFE系の読み出しシーケンスの一例を示す図である。また、図8は、第1の実施形態に係るSSFP系の読み出しシーケンスの一例を示す図である。なお、図7及び8に例示する各読み出しシーケンスは、いずれも3D撮像用の読み出しシーケンスである。具体的には、図7に例示するFE系の読み出しシーケンスは、リードアウト用傾斜磁場Grを反転させることでエコー信号を発生させるものであり、さらに、位相エンコード用傾斜磁場Geについて、印加した傾斜磁場をデータ収集後にキャンセルするためのリワインド用傾斜磁場を印加するものである。また、図8に例示するSSFP系の読み出しシーケンスは、図7に示す読み出しシーケンスにおいて、スライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Ge及びリードアウト用傾斜磁場Grの全てについて、繰り返し時間TRの中で傾斜磁場の積分値がゼロになるようにリワインド用傾斜磁場を印加するものである。
ここで、図8に例示するSSFP系の読み出しシーケンスは、図7に例示するFE系の読み出しシーケンスと比べると、単位収集時間あたりのSN比は高いが、3軸方向にリワインド用傾斜磁場を印加して傾斜磁場の積分値を調整するため、磁場不均一性に弱い。一方、図7に例示するFE系の読み出しシーケンスは、図8に例示するSSFP系の読み出しシーケンスと比べると、磁場不均一性に強く、かつ、空間分解能が高いが、単位収集時間あたりのSN比は低い。
そこで、本実施形態では、MRI装置100は、ST−MI−ASLで被検体頭部血管を撮像する場合に、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも2つの範囲に分割する。そして、MRI装置100は、分割した2つの範囲のうち、待ち時間が短い方の範囲については、磁場不均一性に強いFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については、SN比が高いSSFP系の読み出しシーケンスを用いてデータ収集を行う。
図9は、第1の実施形態に係るMRI装置100の詳細な構成例を示す図である。図9は、図1に示したシーケンサ10及び計算機システム30を示している。また、図9は、計算機システム30が有する機能部のうち、インタフェース部31、画像再構成部32、記憶部33、入力部34、表示部35、及び制御部36を示している。
図9に示すように、記憶部33は、撮像条件記憶部33aと、分割情報記憶部33bと、MRデータ記憶部33cと、画像データ記憶部33dとを有する。
撮像条件記憶部33aは、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を記憶する。なお、ここでいう撮像条件とは、読み出しシーケンスの生成に必要な各種撮像パラメータのパラメータ値であり、例えば、収集マトリクス、NAQ、スラブ厚、1スラブ内のスライス枚数、ショット数、SENSE−Factorなどである。また、撮像条件には、PROPELLER、radial系などのサンプリング法の種類も含まれる。本実施形態では、撮像条件記憶部33aは、少なくとも、FE系の読み出しシーケンスに関する撮像条件と、SSFP系の読み出しシーケンスに関する撮像条件とを記憶する。
分割情報記憶部33bは、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割するための基準となる情報を記憶する。例えば、分割情報記憶部33bは、撮像対象の部位に関する解剖学的な情報を記憶する。本実施形態では、分割情報記憶部33bは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、標識化された血液が主幹動脈を流れる範囲と末梢動脈を流れる範囲とに時間的な範囲を分割するための情報として、MCA(Middle-Cerebral Artery)の左右に分かれる分岐を第1分枝とした場合の第2分枝に至る辺りに標識化された血液が到達するときの待ち時間を記憶する。
なお、分割情報記憶部33bが記憶する情報は、解剖学的な情報に限られない。例えば、過去に行われた撮像において基準として用いられた待ち時間のうち、操作者が最も適していると判断した待ち時間を記憶しておいてもよい。すなわち、待ち時間に関する操作者の経験値を記憶しておいてもよい。
MRデータ記憶部33cは、インタフェース部31を介してシーケンサ10から受信されたMRデータを記憶する。画像データ記憶部33dは、画像再構成部32によってMRデータから再構成された画像を記憶する。
また、制御部36は、撮像条件設定部36aと、分割部36bと、収集部36cと、血流画像生成部36dと、合成部36eとを有する。
撮像条件設定部36aは、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を設定する。具体的には、撮像条件設定部36aは、入力部34を介して、操作者から各種撮像パラメータのパラメータ値の入力を受け付け、受け付けたパラメータ値を撮像条件として撮像条件記憶部33aに登録する。本実施形態では、撮像条件設定部36aは、少なくとも、FE系の読み出しシーケンスに関する撮像条件と、SSFP系の読み出しシーケンスに関する撮像条件とを設定する。
分割部36bは、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも2つの範囲に分割する。具体的には、分割部36bは、分割情報記憶部33bによって記憶された情報に基づいて、時間的な範囲又は空間的な範囲を分割する。本実施形態では、分割部36bは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも2つの範囲に分割する。
例えば、分割部36bは、複数の待ち時間のうち少なくとも1つの待ち時間を閾値として設定することで、撮像時の時間的な範囲を少なくとも2つの範囲に分割する。具体的には、分割部36bは、ST−MI−ASLで用いられる複数の待ち時間のうち、分割情報記憶部33bによって記憶されている待ち時間と一致する待ち時間を閾値として設定する。本実施形態では、分割部36bは、MCAの左右に分かれる分岐を第1分枝とした場合の第2分枝に至る辺りに標識化された血液が到達するときの待ち時間を閾値として設定することで、標識化された血液が主幹動脈を流れる範囲と末梢動脈を流れる範囲とに時間的な範囲を分割する。
収集部36cは、分割部36bによって分割された少なくとも2つの範囲のうち、第1の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、第2の範囲については読み出しシーケンスの種類又は撮像条件が第1の読み出しシーケンスと異なる第2の読み出しシーケンスを用いてデータ収集を行う。
本実施形態では、収集部36cは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、分割された2つの範囲のうち、待ち時間が短い方の範囲についてはFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはSSFP系の読み出しシーケンスを用いてデータ収集を行う。具体的には、収集部36cは、撮像条件記憶部33aによって記憶されている撮像条件に基づいて、FE系のシーケンス情報及びSSFP系のシーケンス情報をそれぞれ生成する。そして、収集部36cは、生成した各シーケンス情報をシーケンサ10に送信し、それらの各シーケンス情報を用いてST−MI−ASLを実行するようにシーケンサ10を制御する。
図10は、従来のST−MI−ASLによるデータ収集を示す図である。また、図11は、第1の実施形態に係る収集部36cによって行われるST−MI−ASLを示す図である。図10に示すように、従来のST−MI−ASLでは、複数の待ち時間TI1〜TI6ごとに、同じ種類の読み出しシーケンス(例えば、FE系の読み出しシーケンス)を用いてデータ収集が行われていた(S1〜S6)。
これに対し、本実施形態では、収集部36cは、分割部36bによって設定された待ち時間の閾値に基づいて、データ収集に用いる読み出しシーケンスを切り替える。例えば、待ち時間の閾値がTI3であった場合には、収集部36cは、図11に示すように、待ち時間がTI3以下であるときは、FE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間がTI3を超えているときは、SSFP系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ10を制御する。
ここで、本実施形態では、分割部36bによって、MCA(Middle-Cerebral Artery)の左右に分かれる分岐を第1分枝とした場合の第2分枝に至る辺りに標識化された血液が到達するときの待ち時間が閾値として設定される。これにより、本実施形態では、標識化された血液が主幹動脈を流れる範囲については、FE系の読み出しシーケンスでデータ収集が行われ、標識化された血液が末梢動脈を流れる範囲については、SSFP系の読み出しシーケンスでデータ収集を行われることになる。
なお、収集部36cは、FE系の読み出しシーケンス及びSSFP系の読み出しシーケンスのどちらを用いてデータ収集を行う場合も、各待ち時間において、タグ画像用のデータ収集とコントロール画像用のデータ収集とをそれぞれ行うようにシーケンサ10を制御する。そして、収集部36cは、インタフェース部31を介してシーケンサ10からMRデータを受信すると、受信したMRデータをMRデータ記憶部33cに格納する。この結果、画像再構成部32によって、ST−MI−ASLで用いられる複数の待ち時間ごとに、シーケンサ10によって収集されたMRデータからタグ画像及びコントロール画像がそれぞれ生成されることになる。
図9に戻って、血流画像生成部36dは、第1の読み出しシーケンスによって収集されたデータ及び第2の読み出しシーケンスによって収集されたデータそれぞれから血流画像を生成する。本実施形態では、血流画像生成部36dは、ST−MI−ASLで用いられる複数の待ち時間ごとに、画像再構成部32によって生成されたタグ画像及びコントロール画像の差分画像を血流画像として生成する。すなわち、血流画像生成部36dは、FE系の読み出しシーケンスによって収集されたデータ及びSSFP系の読み出しシーケンスによって収集されたデータそれぞれについて、待ち時間ごとに血流画像を生成する。
合成部36eは、第1の読み出しシーケンスによって収集されたデータから生成された画像と第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。本実施形態では、合成部36eは、ST−MI−ASLが実行される場合に、血流画像生成部36dによって待ち時間ごとに生成された各血流画像を合成する。すなわち、合成部36eは、FE系の読み出しシーケンスによって収集されたデータから生成された血流画像と、SSFP系の読み出しシーケンスによって収集されたデータから生成された血流画像とを合成する。そして、合成部36eは、合成した血流画像を表示部35に表示させる。
ここで、合成部36eは、第1の読み出しシーケンスによって収集されたデータから生成された血流画像と第2の読み出しシーケンスによって収集されたデータから生成された血流画像の少なくとも一方について、ゲイン補正、T1緩和補正及び画像補間のうち少なくとも1つの補正処理を行った後に、各画像を合成してもよい。本実施形態では、FE系の読み出しシーケンスを用いて得られた血流画像及びSSFP系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方について補正処理を行う。
図12は、第1の実施形態に係る合成部36eによって行われる画像補正を説明するための図である。図12の上段、中段及び下段に示す図において、横軸はST−MI−ASLで用いられる待ち時間TIを示しており、縦軸は血流画像における血流部分の信号強度S.I.(Signal Intensity)を示している。
例えば、図12の上段に示すように、0<TI<TIthの範囲ではFE系の読み出しシーケンスで、TI=TIthではFE系及びSSFP系の両方の読み出しシーケンスで、TIth<TI<TImaxの範囲ではSSPF系の読み出しシーケンスでデータ収集を行うとする。この場合には、FE系及びSSFP系それぞれの読み出しシーケンスでデータ収集を行う際の受信ゲイン及びノイズSD(Noise S.D.)値を同じにしたとすると、血流部分の信号強度は、磁場不均一性が問題とならない範囲(TIが長い範囲)では、FE系の読み出しシーケンスで得られた血流画像よりも、SSFP系の読み出しシーケンスで得られた血流画像のほうが大きい。なお、ここでいうノイズSD値は、画像に含まれる画素ごとの信号強度の標準偏差(Standard Deviation)であり、ノイズの特性を示す指標値である。
そこで、合成部36eは、FE系の読み出しシーケンスを用いて得られた血流画像及びSSFP系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方について、ゲイン補正を行う。例えば、図12の中段に示すように、TI=TIthとしたときにFE系及びSSFP系それぞれの読み出しシーケンスで収集されるデータの信号強度が一致するように、少なくとも一方の読み出しシーケンスでデータ収集を行う際の受信ゲインを調整する。
なお、図12の中段に示す例では、FE系の読み出しシーケンスでデータ収集を行う際の受信ゲインに、SSFP系の読み出しシーケンスでデータ収集を行う際の受信ゲインを合わせた場合を示している。しかし、逆に、SSFP系の読み出しシーケンスでデータ収集を行う際の受信ゲインに、FE系の読み出しシーケンスでデータ収集を行う際の受信ゲインを合わせてもよい。また、FE系の読み出しシーケンスでデータ収集を行う際の受信ゲイン、及び、SSFP系の読み出しシーケンスでデータ収集を行う際の受信ゲインの両方を調整してもよい。
また、図12に示す例では、TI=TIthとしてFE系及びSSFP系の両方の読み出しシーケンスでデータ収集を行うこととしたが、必ずしも、このように異なる同じTIで重複してデータ収集を行わなくてもよい。その場合には、例えば、合成部36eは、FE系の読み出しシーケンスで収集される最後のTIのデータと、SSFP系の読み出しシーケンスで収集される最初のTIのデータとで信号強度が一致するように、少なくとも一方の読み出しシーケンスでデータ収集を行う際の受信ゲインを調整する。
また、合成部36eは、上述したゲイン補正を行った後に、さらにT1緩和補正を行ってもよい。その場合には、例えば、図12の下段に示すように、合成部36eは、TI=0としてFE系の読み出しシーケンスで得られた血流画像における血流部分の信号強度に合うように、他のTIで得られた血流画像における血流部分の信号強度を補正する。例えば、合成部36eは、血液のT1値をT1bとした場合に、各TIで得られた血流画像それぞれの信号値に{1/exp(−TI/T1b)}を乗じることによって、各血流画像における信号強度を補正する。例えば、このT1緩和補正は、ASL−MRPで血液の流量を定量化する際に行われることが望ましい。
また、合成部36eは、上述したゲイン補正を行った後に、さらに画像補間を行ってもよい。例えば、合成部36eは、待ち時間TIごとに生成された血流画像について、連続するTIの2つの血流画像を用いて、重み付き補間により画像を生成し、生成した画像を当該2つの血流画像の間に補間する。
このように、FE系の読み出しシーケンスを用いて得られた血流画像及びSSFP系の読み出しシーケンスを用いて得られた血流画像の少なくとも一方についてゲイン補正を行うことによって、一方の読み出しシーケンスで得られた血流画像の信号強度を他方の読み出しシーケンスで得られた血流画像の信号強度に合わせることができ、合成画像における血管の連続性を向上させることができる。また、ゲイン補正に加えてT1緩和補正及び画像補間を行うことによって、合成画像における血管の連続性をさらに向上させることができる。なお、上述したゲイン補正、T1緩和補正及び画像補間は必ずしも全ての補正処理を行わなくてもよく、いずれか1つ又は2つの補正処理だけを行ってもよい。
図13は、第1の実施形態に係るMRI装置100によって行われる処理の流れを示すフローチャートである。図13に示すように、MRI装置100では、撮像条件設定部36aが、被検体の撮像に用いられる読み出しシーケンスの種類ごとに撮像条件を設定する。本実施形態では、撮像条件設定部36aは、少なくとも、FE系の読み出しシーケンスに関する撮像条件と、SSFP系の読み出しシーケンスに関する撮像条件とを設定する(ステップS101)。
そして、操作者から撮像開始の指示を受け付けると(ステップS102,Yes)、分割部36bが、被検体を撮像する際の時間的な範囲又は空間的な範囲を少なくとも2つの範囲に分割する。本実施形態では、分割部36bは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも2つの範囲に分割する(ステップS103)。
続いて、収集部36cが、分割部36bによって分割された少なくとも2つの範囲のうち、第1の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、第2の範囲については第1の読み出しシーケンスと種類又は撮像条件が異なる第2の読み出しシーケンスを用いてデータ収集を行う。本実施形態では、収集部36cは、分割部36bによって分割された2つの範囲のうち、待ち時間が短い方の範囲についてはFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはSSFP系の読み出しシーケンスを用いてデータ収集を行う(ステップS104)。
血流画像生成部36dは、第1の読み出しシーケンスによって収集されたデータ及び第2の読み出しシーケンスによって収集されたデータそれぞれから血流画像を生成する。本実施形態では、血流画像生成部36dは、ST−MI−ASLで用いられる複数の待ち時間ごとに、画像再構成部32によって生成されたタグ画像及びコントロール画像の差分画像を血流画像として生成する(ステップS105)。
合成部36eは、第1の読み出しシーケンスによって収集されたデータから生成された画像と第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する。本実施形態では、FE系の読み出しシーケンスによって収集されたデータから生成された血流画像と、SSFP系の読み出しシーケンスによって収集されたデータから生成された血流画像とを合成する(ステップS106)。そして、合成部36eは、合成した血流画像を表示部35に表示させる(ステップS107)。
上述したように、第1の実施形態に係るMRI装置100は、ST−MI−ASLで被検体の頭部血管を撮像する場合に、撮像時の時間的な範囲を少なくとも2つの範囲に分割し、待ち時間が短い方の範囲については、磁場均一性に強く、かつ、空間分解能が高いFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については、SN比が高いSSFP系の読み出しシーケンスを用いてデータ収集を行うので、主幹動脈及び末梢動脈の両方を高精度に描出した画像を得ることができる。したがって、第1の実施形態に係るMRI装置100によれば、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる。
なお、第1の実施形態では、ST−MI−ASLで被検体を撮像する場合について説明したが、撮像法がこれに限られない。例えば、本実施形態は、ST−SI(Single-Tag Single-TI)−ASLを用いる場合にも同様に実施が可能である。ST−SI−ASLは、被検体の撮像領域に流入する流体を標識化するためのRF波を1回印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを、待ち時間を変えて複数回実行する撮像法である。
このようなST−SI−ASLで被検体を撮像する場合に、例えば、分割部36bは、ST−MI−ASLを用いる場合と同様に、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間に基づいて、撮像時の時間的な範囲を少なくとも2つの範囲に分割する。そして、収集部36cは、分割部36bによって分割された少なくとも2つの範囲のうち、第1の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、第2の範囲については第1の読み出しシーケンスと種類又は撮像条件が異なる第2の読み出しシーケンスを用いてデータ収集を行う。例えば、収集部36cは、ST−MI−ASLを用いる場合と同様に、分割された2つの範囲のうち、待ち時間が短い方の範囲についてはFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲についてはSSFP系の読み出しシーケンスを用いてデータ収集を行う。
さらに、1回のRF波の印加でk空間全体のデータをサンプリングすることが難しい場合には、k空間を複数の領域に分割し、分割した領域ごとに、標識化用のRF波を印加してデータ収集を行ってもよい。その場合には、分割した領域全てについてデータを収集するまでは、同じ待ち時間でデータ収集を繰り返し、分割した領域全てについてデータを収集した後に、待ち時間を変更して、次のk空間のデータを収集する。このような場合も、同様に、分割部36bが、待ち時間に基づいて、撮像時の時間的な範囲を少なくとも2つの範囲に分割する。また、収集部36cが、第1の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、第2の範囲については第1の読み出しシーケンスと種類又は撮像条件が異なる第2の読み出しシーケンスを用いてデータ収集を行う。
また、第1の実施形態では、収集部36cが、分割部36bによって設定された待ち時間の閾値に基づいて、データ収集に用いる読み出しシーケンスを変更する場合について説明したが、読み出しシーケンスの制御はこれに限られない。例えば、収集部36cは、撮像領域における所定の位置に標識化された流体が到達したことを検出した場合に、データ収集用の読み出しシーケンスを第1の読み出しシーケンスから第2の読み出しシーケンスに変更してもよい。
例えば、収集部36cは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、血流画像生成部36dによって血流画像が生成されるごとに、生成された血流画像における標識化された血液の位置を検出する。そして、収集部36cは、検出した血液がMCAの第2分枝に至る辺りに到達したことを検出した場合に、データ収集に用いる読み出しシーケンスをFE系の読み出しシーケンスからSSFP系の読み出しシーケンスに変更する。
また、第1の実施形態では、収集部36cが、ST−MI−ASLで被検体の頭部血管を撮像する場合に、データ収集用に用いる読み出しシーケンスの種類をFE系からSSFP系に変更する場合について説明したが、実施形態はこれに限られない。例えば、収集部36cは、読み出しシーケンスの種類を変更するのではなく、読み出しシーケンスの撮像条件を変更してもよい。
例えば、収集部36cは、ST−MI−ASLで被検体の頭部血管を撮像する場合に、同じ種類の読み出しシーケンスを用いて、分割部36bによって分割された2つの範囲のうち、待ち時間が短い方の範囲については空間分解重視の撮像条件でデータ収集を行い、待ち時間が長い方の範囲についてはSN比重視の撮像条件でデータ収集を行う。この場合には、撮像条件記憶部33aが、同一種類の読み出しシーケンスについて、空間分解重視の撮像条件とSN比重視の撮像条件とをそれぞれ記憶する。なお、空間分解重視の撮像条件とSN比重視の撮像条件とは、それぞれ、収集マトリクス、NAQ、スラブ厚、1スラブ内のスライス枚数、ショット数などを適宜に調整することによって設定される。
このように読み出しシーケンスの撮像条件を変更することによって、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、撮像法をASL−MRAからASL−MRPへの切り替えることも可能になる。ASL−MRAで空間分解能重視のまま待ち時間TIを延長していくと、灌流はSN比の限界のためノイズレベル以下になって描出されないが、途中でSN比重視モードに変更することによって、標識化された血液が待ち時間TIの延長とともに毛細血管床に移行する様子をシームレスに画像化することができる。
なお、収集部36cは、データ収集に用いるパルシーケンスの種類及び撮像条件のいずれか一方を変更するだけではなく、パルシーケンスの種類及び撮像条件の両方を変更してもよい。これにより、被検体を時系列に撮像する際の時間的な位置や撮像領域の空間的な位置に応じて、より適切な方法でデータを収集することができるようになり、さらに高画質な画像を得ることができる。
(第2の実施形態)
次に、第2の実施形態について説明する。第2の実施形態では、被検体を撮像する際の時間的な範囲又は空間的な範囲の全体について、種類が異なる2つの読み出しシーケンスそれぞれを用いてデータ収集を行う場合の例を説明する。なお、第2の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
本実施形態では、マルチスラブ法を用いたASL−MRAで被検体の頭部血管を撮像する場合の例と、ST−MI−ASLを用いたASL−MRAで被検体の頭部血管を撮像する場合の例とを説明する。ここで、マルチスラブ法は、被検体の撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行う撮像法である。
第2の実施形態では、収集部36cは、被検体を撮像する際の時間的な範囲又は空間的な範囲の全体については第1の読み出しシーケンス及び第2の読み出しシーケンスそれぞれを用いてデータ収集を行う。また、合成部36eは、第1の読み出しシーケンスによって収集されたデータから生成された画像のうち第1の範囲における画像と、第2の読み出しシーケンスによって収集されたデータから生成された画像のうち第2の範囲における画像とを合成する。
図14は、第2の実施形態に係るMRI装置によって行われるマルチスラブ法を用いたASL−MRAを説明するための図である。図14では、マルチスラブ法を用いたASL−MRAで被検体の頭部血管を撮像する場合に、主幹動脈側の分割領域141及び末梢動脈側の分割領域142の2つの領域に撮像領域を分割してデータ収集を行う場合の例を示している。
例えば、図14の左側に示すように、収集部36cは、待ち時間TIごとに、FE系の読み出しシーケンス及びSSFP系の読み出しシーケンスそれぞれでデータ収集を行う。そして、図14の右側に示すように、合成部36eは、FE系の読み出しシーケンスを用いて得られた血流画像のうち主幹動脈側の分割領域141の画像と、SSFP系の読み出しシーケンスを用いて得られた血流画像のうち末梢動脈側の分割領域142の画像とを合成する。これにより、主幹動脈から末梢動脈までが良好に描出された血流画像が得られる。
このとき、例えば、合成部36eは、FE系の読み出しシーケンスで得られた血流画像については、主幹動脈側の分割領域141における信号値の重みを末梢動脈側の分割領域142における信号値の重みより大きくし、SSFP系の読み出しシーケンスで得られた血流画像については、末梢動脈側の分割領域142における信号値の重みを主幹動脈側の分割領域141における信号値の重みより大きくしたうえで、それぞれの血流画像を合成する。
具体的な例として、例えば、合成部36eは、FE系の読み出しシーケンスで得られた血流画像については、分割領域141における信号値の重み:分割領域142における信号値の重み=1:0とし、SSFP系の読み出しシーケンスで得られた血流画像については、分割領域141における信号値の重み:分割領域142における信号値の重み=0:1としたうえで、それぞれの血流画像を合成する。このとき、例えば、合成部36eは、重み付き加算又は最大値投影(Maximum Intensity Projection:MIP)により、血流画像を合成する。
図15は、第2の実施形態に係るMRI装置によって行われるST−MI−ASLを用いたASL−MRAを説明するための図である。図15では、ST−MI−ASLを用いたASL−MRAで、被検体の頭部血管を撮像する場合の例を示している。図15において、左上の図は、FE系の読み出しシーケンスでTI=600msとしたときに得られる画像を示しており、右上の図は、FE系の読み出しシーケンスでTI=1400msとしたときに得られる画像を示している。また、図15において、左下の図は、SSFP系の読み出しシーケンスでTI=600msとしたときに得られる画像を示しており、右下の図は、SSFP系の読み出しシーケンスでTI=1400msとしたときに得られる画像を示している。なお、図15に示す例では、標識化された血液が、TI=600msとしたときに主幹動脈に達し、TI=1400msとしたときに末梢動脈に達することとしている。
例えば、図15に示すように、収集部36cは、複数の待ち時間TIごとに、FE系の読み出しシーケンス及びSSFP系の読み出しシーケンスそれぞれでデータ収集を行う。そして、合成部36eは、TI=600msについてはFE系の読み出しシーケンスで得られた血流画像(図15の左上)を用い、TI=1400msについてはSSFP系の読み出しシーケンスで得られた血流画像(図15の右下)を用いて、複数のTIに渡るダイナミック画像を生成する。これにより、主幹動脈から末梢動脈までが良好に描出されたダイナミック画像が得られる。
このとき、例えば、合成部36eは、TI=600msについては、FE系の読み出しシーケンスで得られた血流画像における信号値の重みをSSFP系の読み出しシーケンスで得られた血流画像における信号値の重みより大きくし、TI=1400msについては、SSFP系の読み出しシーケンスで得られた血流画像における信号値の重みをFE系の読み出しシーケンスで得られた血流画像における信号値の重みより大きくしたうえで、TIごとに、それぞれの血流画像を合成する。
具体的な例として、例えば、合成部36eは、TI=600msについては、FE系の読み出しシーケンスで得られた血流画像における信号値の重み:SSFP系の読み出しシーケンスで得られた血流画像における信号値の重み=1:0とし、TI=1400msについては、FE系の読み出しシーケンスで得られた血流画像における信号値の重み:SSFP系の読み出しシーケンスで得られた血流画像における信号値の重み=0:1としたうえで、TIごとに、それぞれの血流画像を合成する。このとき、例えば、合成部36eは、重み付き加算又は最大値投影(Maximum Intensity Projection:MIP)により、血流画像を合成する。
(第3の実施形態)
次に、第3の実施形態について説明する。第3の実施形態では、標識化用のRF波としてパルス波を用いるPASL(Pulsed ASL)法で被検体の頭部血管を撮像する場合の例を説明する。なお、第3の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
例えば、PASL法で被検体を撮像する場合に、標識化用のRF波を印加する領域の大きさを血管の走行方向に長く設定すると、待ち時間TIを長くしたときに、標識化された血液は、末梢動脈に流入した状態となるとともに、主幹動脈にも流入した状態となる。そのため、待ち時間TIが長いときにSSFP系の読み出しシーケンスを用いると、主幹動脈の部分でフローボイド(flow-void)が出現したり、血管以外の部分にゴーストなどのアーチファクトが出現したりする場合がある。このことから、血管画像を合成する場合に、待ち時間TIを長くしたときの血管画像には、主幹動脈に流入した血液を残さないようにするのが望ましい。
そこで、本実施形態に係るMRI装置は、PASL法で被検体を撮像する際に、標識化用のRF波をタグ領域に印加してから所定の時間が経過した後にタグ領域にサチュレーションパルスを印加することで、RF波によって標識化された血液の後端における血液信号を低減させる。ここで用いられるサチュレーションパルスは、TEC(Tag-End-Cut)パルスと呼ばれる。このように、標識化された血液の後端部分の血液信号を低減させることで、待ち時間TIを長くしたときの血管画像において、主幹動脈に流入した血液を除去することができる。
図16〜18は、第3の実施形態に係るMRI装置によって行われるPASL法を説明するための図である。図16は、第3の実施形態に係るPASL法の読み出しシーケンスを示している。また、図17は、第3の実施形態に係るPASL法によって得られる血流画像を示している。また、図18は、第3の実施形態に係るPASL法によって得られる合成画像を示している。
本実施形態に係るPASL法の読み出しシーケンスでは、標識化用のRF波を印加してから一定の時間Ttecが経過した後に、標識化用のRF波を印加する時間幅と等価なTECパルスを印加する。そして、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間TIを、Ttecの整数倍となるように段階的に延長する。例えば、待ち時間TIを、短い方から、TI(1)=Ttec、TI(2)=2*Ttec、・・・、TI(n−1)=(n−1)*Ttec、TI(n)=n*Ttecと変える(n:自然数)。ここで、例えば、図16に示すように、TIを変える幅を少し短めにして、(Ttec−a)としてもよい(a:定数、a<Ttec)。
これにより、例えば、図17に示すように、所定の幅を有するボーラス(bolus)状の血流が待ち時間TIごとに段階的に進んでいく血流画像が得られる。そして、これらの血流画像を合成することで、例えば、図18に示すように、主幹動脈及び末梢動脈を含んだ頭部全体の血流画像が得られる。なお、図16に示すようにTIを変える幅を(Ttec−a)とした場合には、血流画像を合成する際に、標識化された血液の境界をオーバーラップさせることができる。これにより、TIが異なる血流画像の境界で血管の端部に無信号の部分が生じるのを防ぐことができ、合成画像における血管の連続性を向上させることができる。
そして、本実施形態では、このようなPASL法で被検体を撮像する場合に、収集部36cが、所定の待ち時間TIを境に、データ収集用の読み出しシーケンスを、FE系の読み出しシーケンスからSSFP系の読み出しシーケンスに変更する。例えば、図16に示すように、収集部36cは、待ち時間の閾値がTI3であった場合には、収集部36cは、図11に示すように、待ち時間がTI3以下であるときは、FE系の読み出しシーケンスを用いてデータ収集(S1〜S3)を行い、待ち時間がTI3を超えているときは、SSFP系の読み出しシーケンスを用いてデータ収集(S4〜S6)を行うように、シーケンサ10を制御する。
なお、ここでは、PASL法を用いる場合について説明したが、本実施形態は、例えば、CASL(Continuous ASL)法やPCASL(Pulsed Continuous ASL)法を用いる場合も同様に実施が可能である。CASL法は、標識化用のRF波として連続波を用いる方法であり、PCASL法は、CASL法の実用化を目的としたもので、標識化用のRF波として、短いパルス波を複数用いる方法である。
CASL法又はPCASL法では、Ttecは、標識化用のRF波の印加時間Ttagに相当し、TECパルスを印加してからデータ収集を行うまでの待ち時間は、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間Tpld(post labeling delay)に相当する。そのため、例えば、CASL法又はPCASL法を用いる場合は、TpldをTtagの整数倍となるように段階的に延長する。そして、収集部36cは、PASL法を用いる場合と同様に、所定の待ち時間TIを境に、データ収集用の読み出しシーケンスを、FE系の読み出しシーケンスからSSFP系の読み出しシーケンスに変更する。
さらに、本実施形態は、上述したPASL法とマルチスラブ法とを組み合わせて用いる場合も同様に実施が可能である。この場合には、MRI装置は、被検体の撮像領域を複数の分割領域に分割し、撮像領域に流入する流体を標識化するためのRF波を印加してから複数の分割領域にそれぞれに対応する複数の待ち時間が経過するごとにデータ収集を行う撮像法で被検体を撮像する撮像法を行う。そして、分割部36bは、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間に基づいて、空間的な範囲を少なくとも2つの範囲に分割する。また、収集部36cは、分割された2つの範囲のうち、待ち時間が短い方の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、待ち時間が長い方の範囲については第2の読み出しシーケンスを用いてデータ収集を行う。そして、合成部36eは、第1の読み出しシーケンスによって収集されたデータから生成された画像と第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成することで撮像領域の画像を生成する。つまり、合成部36eは、複数の分割領域の画像を合成することで、分割前の撮像領域の画像を生成する。
図19は、第3の実施形態に係るMRI装置によって行われるマルチスラブ法を用いたPASL法を説明するための図である。図19では、マルチスラブ法を用いたPASL法で被検体の頭部血管を撮像する場合に、主幹動脈側の分割領域191と末梢動脈側の分割領域192とに撮像領域を分割してデータ収集を行う場合の例を示している。この場合に、例えば、分割部36bは、被検体の頭部に設定された撮像領域を、待ち時間TIに基づいて、待ち時間T1(1)及びTI(2)に対応する主幹動脈側の分割領域191と、待ち時間T1(3)及びTI(4)に対応する末梢動脈側の分割領域192とに分割する。
そして、この場合に、例えば、収集部36cは、待ち時間T1(1)及びTI(2)についてはFE系の読み出しシーケンスを用いてデータ収集を行い、待ち時間T1(3)及びTI(4)についてはSSFP系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ10を制御する。この場合も、例えば、図17に示したように、所定の幅を有するボーラス状の血流が待ち時間TIごとに段階的に進んでいく血流画像が得られる。そして、これらの血流画像を合成することで、例えば、図18に示したように、主幹動脈及び末梢動脈を含んだ頭部全体の血流画像が得られる。
このように、短い待ち時間TIに対応する主幹動脈側の分割領域については、磁場均一性に強く、かつ、空間分解能が高いFE系の読み出しシーケンスを用いてデータ収集を行い、長い待ち時間TIに対応する末梢動脈側の分割領域については、SN比が高いSSFP系の読み出しシーケンスを用いてデータ収集を行うことによって、主幹動脈及び末梢動脈の両方を高精度に描出した画像を得ることができる。
なお、ここでは、収集部36cが、FE系の読み出しシーケンスからSSFP系の読み出しシーケンスに読み出しシーケンスの種類を変更する場合の例を説明したが、例えば、読み出しシーケンスの条件を変更してもよい。その場合には、例えば、収集部36cは、同じ種類の読み出しシーケンスを用いて、短い待ち時間TIに対応する主幹動脈側の分割領域については、空間分解重視の撮像条件でデータ収集を行い、長い待ち時間TIに対応する末梢動脈側の分割領域については、SN比重視の撮像条件でデータ収集を行う。
さらに、例えば、収集部36cは、単一の読み出しシーケンスで、血管信号が多い分割領域からデータを収集し、血管信号が少ない分割領域からはデータを収集しないように、シーケンサ10を制御してもよい。なお、この場合には、各分割領域は、データを収集する領域としない領域とを細かく区分けできるように、それぞれ血管の走行方向における幅を狭く設定するのが望ましい。
(第4の実施形態)
次に、第4の実施形態について説明する。第4の実施形態では、被検体の撮像領域に流入する流体を標識化するためのRF波を印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを待ち時間を変えて複数回実行する撮像法で被検体を撮像する場合の例を説明する。なお、第4の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
例えば、このような撮像法の一例として、Hadamard Encoding(HE)法を用いたMulti−Tag(MT)−PASLであるHE−MT−PASLがある。HE法については、例えば、CASL法と組み合わせた方法が、Wells JA et al. In vivo hadamard encoded continuous arterial spin labeling (H-CASL). MRM 63:1111-1118 (2010)により提案されている。HE−MT−PASLでは、タグ画像用のRF波とコントロール画像用のRF波とを所定の時間間隔ごとに効率よく配置してデータを収集し、収集した画像を加減算することで、待ち時間TIが異なる複数のASL画像が得られる。この方法によれば、ST−SI−ASLやST−MI−ASLと比べて、SN比をより向上させることができる。
図20は、第4の実施形態に係るMRI装置によって行われるHE−MT−PASLを説明するための図である。HE−MT−PASLでは、例えば、異なる待ち時間TIの組み合わせ数Nが、3、7、15・・・など、すなわち、N=2n−1(n:自然数)に設定される。図20では、異なる待ち時間TIの組み合わせ数をN=7(n=3)とした場合の例を示している。
なお、図20に示す複数のRF波(tag/control RF)のうち、左下がりの斜線が付されたものはタグ画像用のRF波であり、右下がりの斜線が付されたものはコントロール画像用のRF波である。図20に示すように、例えば、N=7とした場合には、2種類の画像用のRF波の順序が異なる8種類((N+1)種類)の読み出しシーケンスの組み合わせが用いられる。また、各読み出しシーケンスにおける待ち時間TI(i)は、k=1〜8の各パターンで同じであれば、任意の時間間隔でよい。
ここで、k=1〜8の各パターンで収集されるデータの複素信号をSkとし、各TI(i)における4倍の差分信号S{TI(i)}、i=1〜7は、以下の式で求められる。
S{TI(i)}=4[S{TI(i)}−Scont{TI(i)}]
すなわち、i=1〜7とした場合には、S{TI(1)}〜S{TI(7)}は、以下の式により求められる。
S{TI(1)}=S1−S2+S3−S4+S5−S6+S7−S8
S{TI(2)}=S1+S2−S3−S4+S5+S6−S7−S8
S{TI(3)}=S1−S2−S3+S4+S5−S6−S7+S8
S{TI(4)}=S1+S2+S3+S4−S5−S6−S7−S8
S{TI(5)}=S1−S2+S3−S4+S5−S6+S7−S8
S{TI(6)}=S1+S2−S3−S4−S5−S6+S7+S8
S{TI(7)}=S1−S2−S3+S4−S5+S6+S7−S8
そして、この場合のSN比は、1回あたりの差分のNAQ=4回加算に相当するsqrt(4)=2倍となる。
また、1回の繰り返し時間をTrepeatとすると、撮像時間は、N組のTIに対して、(N+1)*Trepeatとなる。一方、標識用のRF波を1回印加するごとにデータ収集を行う場合には、撮像時間は、N組のTIに対して、2*(N+1)/2*N*Trepeat=N(N+1)*Trepeatとなる。したがって、HE−MT−SI法では、標識用のRF波を1回印加するごとにデータ収集を行う場合に比べて、撮像時間比はN+1/N(N+1)=1/Nとなる。つまり、N=7とした場合には、標識用のRF波を1回印加するごとにデータ収集を行う場合と同等のSN比のTIの異なる差分画像が、1/7の時間で得られることになる。
なお、ASL−MRPにおいて一般的なST−PASLでは、十分なSN比を確保するためには数回程度の加算平均が行われる。これに対し、上述したHE−MT−PASLによれば、同一のデータが1回の収集のみで得られ、かつ、数段階のTIの加算平均と等価な画像が得られるので、収集時間を有効に短縮することができる。なお、コントロール画像用のRF波は、MRP用でMTCが問題となる場合のみ印加されればよい。したがって、血管のMRA用などでMTCが問題とならない場合には、コントロール画像用のRF波は省略されてもよい。
また、上述したHE−MT−PASLでは、血液の流速に応じたパルス間隔の制御など、タグ条件に制約がある。さらに、HE法では、血液の流速が時間的に変動する場合、例えば拍動流下で用いられる場合には、同一時刻でのRF波のタイミングで血液の流速が同じでないと、画像を加減算する際、特に減算の際に、不要なタイミングで標識化された流体信号が意図した通りに消えない場合もある。これは、PASL法と組み合わせた場合でも、CASL法と組み合わせた場合でも同様である。そこで、例えば、収集部36cは、同一時刻のRF波が同一の心位相で行われるように、ECGセンサ21及びECGユニット22によって検出される心拍のゲート信号におけるR波をトリガ信号として、最初のRF波を印加するようにしてもよい。
そして、本実施形態では、収集部36cは、上述したHE−MT−PASLを用いて被検体を撮像する場合に、k=1〜8の8種類のパターンのうち、いずれか4種類のパターンについてはFE系の読み出しシーケンスでデータ収集を行い、残りの4種類のパターンについてはSSFP系の読み出しシーケンスでデータ収集を行う。この場合には、平均効果によって、SSFP系の読み出しシーケンスで得られる画像とFE系の読み出しシーケンスで得られる画像とが混じった画像が得られる。
例えば、図20に示すように、収集部36cは、8種類のパターンを順に実行する際に、FE系の読み出しシーケンスとSSFP系の読み出しシーケンスとを交互に用いるように、シーケンサ10を制御する。すなわち、収集部36cは、k=1、3、5、7のパターンについてはSSFP系の読み出しシーケンスでデータ収集を行い、k=2、4、6、8のパターンについてはFE系の読み出しシーケンスでデータ収集を行うように制御する。
または、例えば、収集部36cは、TIが短い側にタグ画像用のRF波を多く配置したパターンについてはFE系の読み出しシーケンスを用い、TIが長い側にタグ画像用のRF波を多く配置したパターンについてはFE系の読み出しシーケンスを用いるようにしてもよい。
さらに、被検体の撮像領域に流入する流体を標識化するためのRF波を印加してから所定の待ち時間が経過した後にデータ収集を行うパターンを待ち時間を変えて複数回実行する撮像法の他の例として、ST−MI−ASLを複数回実行する撮像法がある。この撮像法では、標識化用のRF波を印加してから最初のデータ収集を行うまでの待ち時間を変えながら、ST−MI−ASLのパターンを複数回実行し、収集された複数のデータを、RF波が印加されてからデータ収集を行うまでの待ち時間が同じものごとに加算する例について説明する。
図21及び22は、第4の実施形態に係るMRI装置によって行われるST−MR−ASLを説明するための図である。図21に示すように、例えば、収集部36cは、ST−MI−ASLによるデータ収集をM回実行する。この場合に、例えば、収集部36cは、標識用のRF波を印加してから最初のデータ収集までの待ち時間を所定の時間幅Tintの整数倍となるように変えながら、ST−MI−ASLによるデータ収集をN回実行する。また、例えば、収集部36cは、ST−MI−ASLによるデータ収集を実行するごとに、標識化用のRF波を印加してから所定の時間幅Tintの整数倍の時刻ごとに、N個のMRデータを収集する。なお、本実施形態では、標識用のRF波が印加されてから最初のデータ収集を行うまでの待ち時間と、2回目のデータ収集において標識用のRF波が印加されてからデータ収集を行うまでの経過時間とを、それぞれTIと呼ぶ。
ここで、例えば、収集部36cは、各TIをTImnと表した場合に、以下に示すTI11、TI21、TI31、TI41のように、標識用のRF波が印加されてから最初のデータ収集を行うまでの待ち時間を変化させながら、ST−MI−PASLによるデータ収集を実行する。そして、収集部36cは、1回目については、待ち時間TI11が経過した後にデータ収集を行い、その後、TI12、TI13、TI14が経過するごとにデータ収集を実行する。また、収集部36cは、2回目については、待ち時間TI21が経過した後にデータ収集を行い、その後、TI22、TI23、TI24が経過するごとにデータ収集を実行する。また、収集部36cは、3回目については、待ち時間TI31が経過した後にデータ収集を行い、その後、TI32、TI33、TI34が経過するごとにデータ収集を実行する。また、収集部36cは、4回目については、待ち時間TI41が経過した後にデータ収集を行い、その後、TI42、TI43、TI44が経過するごとにデータ収集を実行する。
TI14=7Tint
TI13=TI24=6Tint
TI12=TI23=TI34=5Tint
TI11=TI22=TI33=TI44=4Tint
TI21=TI32=TI43=3Tint
TI31=TI42=2Tint
TI41=1Tint
そして、収集部36cは、収集された複数のMRデータを補正したうえで、標識用のRF波が印加されてからデータ収集を行うまでの経過時間が同じものごとに加算する。例えば、収集部36cは、加算されるMRデータをS(TI)とし、TIに依存したT1減衰を補正するための補正係数をa(TI)とし、データを収集する順番nに依存したT1減衰を補正するための補正係数をbnとした場合に、以下の式によりS(TI)を求める。
S(7Tint)=a(7Tint)*b4*S14
S(6Tint)=a(6Tint)*{b3*S13+b4*S24}/2
S(5Tint)=a(5Tint)*{b2*S12+b3*S23+b4*S34}/3
S(4Tint)=a(4Tint)*{b1*S11+b2*S22+b3*S33+b4*S44}/4
S(3Tint)=a(5Tint)*{b1*S21+b2*S32+b3*S43}/3
S(2Tint)=a(2Tint)*{b1*S31+b2*S42}/2
S(1Tint)=a(1Tint)*b1*S41
この方法によれば、図21に示すように、M+N−1個の異なるTIのMRデータが得られ、かつ、中間のTIほどSN比が向上する。なお、TIが中間より短い場合にはもともとSN比が大きいので問題ないが、TIが中間より長い場合にはSN比は低下し、最も長いTIのSN比は収集されるデータが1個なので、1組のMT−SI−PASLと同じになる。そのため、Tintごとに標識用のRF波とデータ収集とを配置するのが望ましい。さらに、N回のデータ収集の組と最初のデータ収集までの待ち時間TIをTintずつシフトさせたm組の独立したtagのM組として、m=nとするのがより望ましい。なお、m、nは異なってもよい。
また、TIが最も短くなる最後(m=M)のRF波と1回目(n=1)のデータ収集までの時間はTintである必要はない。また、TIは、Tintの整数倍であればよく、ここでいう整数は必ずしも連続していなくてもよい。抜けたTIについては、合成計算時にゼロとしてないものとすればよい。
なお、ASLでは、TIが長いほどSN比は低下するので加算回数を増やしたいが、ST−MI−PASLでは時間的に後のデータほどSN比は低下するため、あまり多くのdataは収集できない。したがって、nは数個にしておき、独立したRF波で収集したm組のデータから合成して中間以降の長いTIでSN比の小さなデータを捨てたとしても、処理効率を向上できる。
そして、本実施形態では、収集部36cは、上述した撮像法を用いて被検体を撮像する場合に、標識用のRF波が印加されてから最初のデータ収集を行うまでの待ち時間が短いパターンについてはFE系の読み出しシーケンスでデータ収集を行い、標識用のRF波が印加されてから最初のデータ収集を行うまでの待ち時間が長いパターンについてはSSFP系の読み出しシーケンスでデータ収集を行う。
例えば、図21に示すように、収集部36cは、標識用のRF波が印加されてから最初のデータ収集を行うまでの待ち時間をTintずつ短くしながらST−MT−ASLのパターンを4回実行する場合に、1回目及び2回目のパターンについてはSSFP系の読み出しシーケンスでデータ収集を行い、3回目及び4回目のパターンについてはSSFP系の読み出しシーケンスでデータ収集を行うように、シーケンサ10を制御する。
または、収集部36cは、各パターンにおけるデータ収集ごとに、標識用のRF波が印加されてからデータ収集を行うまでの待ち時間が短いものについてはFE系の読み出しシーケンスを用い、標識用のRF波が印加されてからデータ収集を行うまでの待ち時間が長いものについてはSSFP系の読み出しシーケンスを用いるようにしてもよい。
例えば、図22に示すように、収集部36cは、4回のパターンについて、標識用のRF波が印加されてからデータ収集を行うまでの待ち時間が4Tint以下であるデータ収集S11、S21、S22、S31〜S33、S41〜S44についてはFE系の読み出しシーケンスを用い、標識用のRF波が印加されてからデータ収集を行うまでの待ち時間が4Tintより長いデータ収集S12〜S14、S23〜S24、S34についてはSSFP系の読み出しシーケンスを用いるように、シーケンサ10を制御する。
上述した第1〜4の実施形態に係るMRI装置によれば、単一の読み出しシーケンスを用いる場合と比べて、全ての待ち時間TIにおいて、磁化率効果の影響が少ない、かつ、SN比に優れた、臓器全体のMRA画像又はMRP画像を得ることができる。また、第1〜4の実施形態に係るMRI装置によれば、例えば脳の動脈血管を撮像する場合に、主幹動脈だけでなく、末梢動脈における血流の描出能を向上させることができる。また、第1〜4の実施形態に係るMRI装置によれば、標識化用のRF波を印加してから標識化された血液が到達するまでの遅延時間が長い組織の部分について、SN比が高い灌流画像が得られる。
(第5の実施形態)
次に、第5の実施形態について説明する。第5の実施形態では、被検体の撮像領域を複数の分割領域に分割して当該複数の分割領域ごとにデータ収集を行う撮像法で被検体の頭部血管を撮像する場合の例を説明する。なお、第5の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
例えば、このような撮像法の一例として、マルチスラブ法を適用した3DTOF−MRAであるマルチスラブ3DTOF−MRAがある。マルチスラブ3DTOF−MRAでは、標識化用のRF波を印加する代わりに、縦磁化が十分に回復した血液のインフロー(in-flow)効果と、静止部における縦磁化の抑制効果とによるコントラスト差を利用するため、撮像領域(imaging slab)のみが設定される。
このようなマルチスラブ3DTOF−MRAで被検体を撮像する場合に、分割部36bは、分割領域に基づいて、被検体を撮像する際の空間的な範囲を少なくとも2つの範囲に分割する。また、収集部36cは、分割部36bによって分割された2つの範囲のうち、主幹動脈を含む範囲については第1の読み出しシーケンスを用いてデータ収集を行い、末梢動脈を含む範囲については第2の読み出しシーケンスを用いてデータ収集を行う。
図23は、第5の実施形態に係るMRI装置によって行われるマルチスラブ3DTOF−MRAを説明するための図である。図23に示すように、例えば、マルチスラブ3DTOF−MRAでは、被検体の撮像領域を4つの分割領域に分割して、それぞれの分割領域についてデータ収集が行われる。
この場合に、例えば、分割部36bは、標識化された血液が主幹動脈を流れる範囲と、末梢動脈を流れる範囲とに空間的な範囲を分割する。また、例えば、収集部36cは、ウィリス動脈輪部の近傍にある主幹動脈部については、磁場不均一性に強いFE系の読み出しシーケンスを用いてデータ収集を行い、頭頂などの末梢動脈部については、SN比に優れたSSFP系の読み出しシーケンスを用いてデータ収集を行うように、シーケンサ10を制御する。
なお、上述したマルチスラブ3DTOF−MRAでは、撮像領域に照射するRF波のフリップ角に傾斜をつけることで末梢動脈の描出能の向上を図るTONE法が併用されてもよい。また、2Dマルチスライス法でも適用が可能である。その場合、読み出しシーケンスの撮像条件を変えてもよい。TOF法は基本的にはGREであるが、その中の撮像条件うち、TONE法での傾斜は対象スラブ内の血流速に応じて変化させてもよい。すなわち、流速が遅いほど傾斜を急にする。またTEとGMN(Gradient Moment Nulling)の組み合わせを、短いTE+0次GMNと、長いTE+1次GMNとで使い分ける。脳の主幹動脈部など磁化率が大きく流速が速ければ前者、抹消で逆なら後者にすることで、位相分散による信号低下を抑制する。さらに、BW(BandWidth)を変えたものを組み合わせても良い。抹消は狭いBWでSN比を向上させ、主幹部は広めのBWで位相分散による信号低下を抑制する。
(第6の実施形態)
次に、第6の実施形態について説明する。第6の実施形態では、マルチスラブ法による撮像をTOF−MRA及びASL−MRA/MRPと組み合わせたマルチスラブASLで被検体の頭部血管を撮像する場合の例を説明する。なお、第2の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
マルチスラブASLでは、被検体の撮像領域を複数の分割領域に分割し、各分割領域について、当該分割領域から一定の距離だけ流体の上流側に離れた標識化領域に当該分割領域に流入する流体を標識化するためのRF波を印加してから一定の待ち時間が経過した後にデータ収集を行う撮像法で被検体を撮像する。
図24は、第6の実施形態に係るMRI装置によって行われるマルチスラブASLを説明するための図である。図24に示すように、例えば、マルチスラブASLでは、撮像領域(imaging slab)を4つの分割領域121〜124に分割する。そして、4つの分割領域それぞれのデータを収集する際に、各分割領域から一定の距離だけ血流の上流側に離れたタグ領域(tag slab)に標識化用のRF波を印加する。すなわち、マルチスラブASLでは、データ収集を行う分割領域を移動しながら、分割領域の移動にともなって、タグ領域を分割領域から一定の距離を保った状態で追随させて移動する。なお、図24では図示を省略しているが、コントロール領域についても、それぞれの分割領域に対して相対的な位置を同じにして追随させる。
このようなマルチスラブASLによれば、複数の分割領域についてタグ領域の位置を撮像領域全体の上流側に固定する場合と比べて、下流側の分割領域になるほど、標識化された血液が到達するまでの時間が短くなり、SN比を向上させることができる。そのため、標識化用のRF波を印加してからデータ収集を行うまでの待ち時間TIを延長しなくても、タグ領域から遠い分割領域内にある血管(MRA)や灌流(MRP)も良好に描出することができる。
このようなマルチスラブASLで被検体の頭部血管を撮像する場合に、例えば、分割部36bは、分割領域に基づいて、撮像領域を少なくとも2つの範囲に分割する。また、収集部36cは、分割部36bによって分割された範囲のうち、主幹動脈を含む範囲については第1の読み出しシーケンスを用いてデータ収集を行い、末梢動脈を含む範囲については第2の読み出しシーケンスを用いてデータ収集を行う。
例えば、図24に示すように、収集部36cは、4つの分割領域のうち、主幹動脈を含む上流側の分割領域121及び122については、FE系の読み出しシーケンスを用いてデータ収集を行う。また、収集部36cは、末梢動脈を含む下流側の分割領域123及び124については、SSFP系の読み出しシーケンスを用いてデータ収集を行う。これにより、被検体の頭部血管が全体にわたって良好に描出することができる。
上述した第5及び6の実施形態に係るMRI装置によれば、単一の読み出しシーケンスを用いる場合と比べて、撮像時間を延長することなく、磁化率効果の影響の少ない、かつ、SN比が優れた画像を得ることができる。また、第5及び6の実施形態に係るMRI装置によれば、例えば脳の動脈血管を撮像する場合に、主幹動脈だけでなく、末梢動脈における血流の描出能を向上させることができる。
(第7の実施形態)
次に、第7の実施形態について説明する。第7の実施形態では、拡散イメージングで被検体を撮像する場合の例を説明する。なお、第7の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
近年、拡散イメージングは、脳神経領域だけでなく、腫瘍検出用に全身臓器でも使用されてきている。また、拡散イメージングでは、一般的に、SE系のシングルショットSE−EPI(S−EPI)の読み出しシーケンスが用いられる。
ここで、脳神経領域では、小脳、前髄など脳底部や頸部、胸部の脊髄などの部分、腹部領域では、肺や肺に近い肝臓部や直腸や大腸に近い部分は、磁場不均一性が大きくなることが知られている。しかし、S−EPIの読み出しシーケンスは磁場不均一性に弱く、歪みが大きくなったり、別のボクセルに描出されるべき信号が1つのボクセル内で重なって描出されてしまったりする場合もある。
また、拡散イメージングで用いられる他の読み出しシーケンスには、FSE系(例えば、HASTE:Half-Fourier Single shot Turbo spin Echo)や、ナビゲータエコー(navigator echo)を同時に収集して動きによるショットごとの位相差を補正するマルチショットEPI(M−EPI)がある。さらに、k空間の軌跡(trajectory)の違いによって、non−Cartesian系のPROPELLER(JET)、radial系などがある。しかし、これらの読み出しシーケンスは、S−EPIの読み出しシーケンスと比べて、単位時間当たりのSN比が低いなどの理由で、1種類の読み出しシーケンスや撮像条件で全身などの広い領域を撮像するのは難しかった。
このような拡散イメージングで被検体を撮像する場合に、分割部36bは、静磁場強度(B0)の分布に基づいて、被検体を撮像する際の空間的な範囲を少なくとも2つの範囲に分割する。例えば、分割部36bは、deltaB0の標準偏差に関する所定の閾値に基づいて、被検体を撮像する際の撮像領域を2つに分割する。
また、収集部36cは、分割部36bによって分割された2つの範囲のうち、静磁場不均一性が大きい方の範囲については第1の読み出しシーケンスを用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については第2の読み出しシーケンスを用いてデータ収集を行う。例えば、収集部36cは、2つの範囲のうち、deltaB0の標準偏差が大きい方を静磁場不均一性が大きい範囲とし、deltaB0の標準偏差が小さい方を静磁場不均一性が小さい範囲とする。なお、撮像対象が頭部であれば、静磁場不均一性が大きい方の範囲は頸部や脳底部などに相当し、静磁場不均一性が小さい方の範囲は頭頂部などに相当する。また、撮像対象が腹部であれば、静磁場不均一性が大きい方の範囲は肺や肝臓の肺に近い部分などに相当し、静磁場不均一性が小さい方の範囲は腔腹膜臓器などに相当する。
例えば、収集部36cは、静磁場不均一性が大きい方の範囲については、歪みに強いFSE系の読み出しシーケンスを用いてデータ収集を行う。また、収集部36cは、静磁場不均一性が小さい方の範囲については、FSE系と比べて歪みには強くないが、SN比が高いS−EPI系の読み出しシーケンスを用いてデータ収集を行う。
また、例えば、収集部36cは、静磁場不均一性が大きい方の範囲については、FSE系の読み出しシーケンスでPROPELLER(JET)を用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については、S−EPIの読み出しシーケンスでデータ収集を行ってもよい。
なお、例えば、収集部36cは、静磁場不均一性が大きい方の範囲については、同じ読み出しシーケンス又は異なる読み出しシーケンスを用いて、歪みに強い撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、SN比が高い撮像条件でデータ収集を行ってもよい。このとき、例えば、収集部36cは、歪みに強い撮像条件として、SENSE−Factorを高くしたものを用い、SN比が高い撮像条件として、SENSE−Factorを低くしたものを用いる。
また、例えば、収集部36cは、同じEPIの読み出しシーケンスを用いて、静磁場不均一性が大きい方の範囲については、ショット数を多くした撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、ショット数を少なくした撮像条件でデータ収集を行ってもよい。例えば、収集部36cは、静磁場不均一性が大きい方の範囲については、M−EPIの読み出しシーケンスを用いてデータ収集を行い、静磁場不均一性が小さい方の範囲については、S−EPIの読み出しシーケンスでデータ収集を行う。
また、拡散イメージングでは、MPG強度やb−valueが大きくなると、渦電流による歪みも大きくなることが知られている。そこで、例えば、収集部36cは、静磁場不均一性が大きい方の範囲については、MPG強度やb−valueを小さくした撮像条件でデータ収集を行い、静磁場不均一性が小さい方の範囲については、MPG強度やb−valueを大きくした撮像条件でデータ収集を行ってもよい。
(第8の実施形態)
次に、第8の実施形態について説明する。第8の実施形態では、定量値画像(パラメータ画像)を撮像する場合の例を説明する。なお、第8の実施形態に係るMRI装置の構成は、基本的には、図1及び9に示したものと同じである。
近年、検査間での比較が容易なことから、T1やT2、プロトン密度、ADC、血流量、ケミカルシフト、磁化率、温度などのMRIに関する各種パラメータを画像化することが行われている。
このように、所定の定量値を空間マッピングするための撮像法で被検体を撮像する場合に、分割部36bは、対象部位の特性に応じて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割する。そして、収集部36cは、分割部36bによって分割された複数の範囲ごとに、各範囲に含まれる対象部位の特性に応じて、データ収集用の読み出しシーケンスの種類又は撮像条件を変更する。
以上、第1〜8の実施形態に係るMRI装置について説明したが、実施形態はこれに限られない。例えば、上述した実施形態では、分割部36bが、撮像対象の部位に関する解剖学的な情報に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割する場合について説明したが、分割の方法はこれに限られない。
例えば、分割部36bは、シミング用に測定された静磁場強度の分布に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割してもよい。例えば、分割部36bは、deltaB0に関する標準偏差や平均値などの統計値に基づいて、時間的な範囲又は空間的な範囲を少なくとも2つの範囲に分割する。この場合には、分割部36bは、deltaB0に関する統計値について所定の閾値を設定し、その閾値を境に、時間的な範囲又は空間的な範囲を分割する。
また、例えば、分割部36bは、撮像領域の位置決めに用いられる画像に基づいて、被検体を撮像する際の時間的な範囲又は空間的な範囲を分割してもよい。例えば、分割部36bは、位置決め画像として撮像されたアキシャル像、コロナル像又はサジタル像に対して操作者が指定した位置を境に、時間的な範囲又は空間的な範囲を分割する。この場合、例えば、分割部36bは、3DTOFの画像を取得済みであれば、3DTOFの画像を位置決め用の画像として用い、未取得であれば、2DPC(Phase Contrast)法で取得した画像を位置決め用の画像として用いてもよい。
また、上述した実施形態では、FE系の読み出しシーケンスとSSFP系の読み出しシーケンスとを用いた場合の例を中心に説明したが、読み出しシーケンスの種類はこれに限られない。すなわち、収集部36cは、分割部36bによって分割された少なくとも2つの範囲について、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの種類を適宜に使い分ける。例えば、収集部36cは、2つの範囲のうち、一方の範囲でSN比が求められる場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも単位収集時間あたりのSN比が高い読み出しシーケンスを用いる。また、例えば、収集部36cは、2つの範囲のうち、一方の範囲で空間分解能が求められる場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも空間分解能が高い読み出しシーケンスを用いる。また、例えば、収集部36cは、2つの範囲のうち、一方の範囲が静磁場不均一性の影響を受けやすい場合は、その範囲については、他方の範囲で用いる読み出しシーケンスよりも静磁場不均一性に強い読み出しシーケンスを用いる。
また、収集部36cは、分割部36bによって分割された少なくとも2つの範囲について、撮像対象の部位や撮像の目的に応じて、読み出しシーケンスの撮像条件を適宜に使い分けてもよい。例えば、収集部36cは、2つの範囲のうち、一方の範囲でSN比が求められる場合は、その範囲については、他方の範囲で用いる撮像条件よりも単位収集時間あたりのSN比が高い撮像条件を用いる。また、例えば、収集部36cは、2つの範囲のうち、一方の範囲で空間分解能が求められる場合は、その範囲については、他方の範囲で用いる撮像条件よりも空間分解能が高い撮像条件を用いる。また、例えば、収集部36cは、2つの範囲のうち、一方の範囲が静磁場不均一性の影響を受けやすい場合は、その範囲については、他方の範囲で用いる撮像条件よりも静磁場不均一性に強い撮像条件を用いる。
また、上述した実施形態では、被検体の頭部血管を撮像する場合の例を中心に説明したが、撮像対象となる部位は頭部に限られない。例えば、撮像対象の部位は、腹部などの他の部位であってもよい。また、上述した実施形態では、血流画像を生成する場合の例を中心に説明したが、撮像対象となる流体は血液に限らない。例えば、撮像対象の流体は、CSF(cerebrospinal fluid)などの他の流体であってもよい。
前述したように、読み出しシーケンスの撮像条件は、所定の撮像領域を撮像する際に用いられる撮像パラメータのパラメータ値である。つまり、撮像条件が異なるということは、あるスライス位置のデータを収集する際の撮像パラメータのパラメータ値が異なることを意味する。例えば、撮像領域を時間的な範囲に分割する場合には、同じスライス位置について、分割された時間的な範囲ごとに、読み出しシーケンスの撮像パラメータのパラメータ値を変えてデータ収集が行われる。また、例えば、撮像領域を空間的な範囲に分割する場合には、異なるスライス位置について、分割された空間的な範囲ごとに、読み出しシーケンスの撮像パラメータのパラメータ値を変えてデータ収集が行われる。
上述した第1〜8の実施形態に係るMRI装置によれば、時間的な範囲又は空間的な範囲の異なる全ての部分について、単一の読み出しシーケンス種又は単一の撮像条件で撮像する場合と比べて、ほぼ同一の撮像時間で、よりアーチファクトの少ない、空間分解能やSN比に優れたMRI画像を得ることができる。
また、上述した第1〜8の実施形態では、撮像領域を時間的又空間的な範囲に分割し、分割された範囲ごとに、読み出しシーケンスの種類及び撮像条件の少なくとも一方を変える場合の例を説明したが、実施形態はこれに限られない。例えば、分割された範囲ごとに、読み出しシーケンスの前に印加されるプリパルスの条件を変えてもよい。
例えば、分割された範囲ごとに、ASL法で用いられる標識化用のRFパルスのタグ条件を変えてもよい。例えば、第6の実施形態で説明したように、収集部が、マルチスラブASLで被検体を撮像する場合に、時間的又空間的に分割された範囲ごとに、データ収集を行う分割領域の移動に追随させて、標識化用のRFパルスを印加する領域すなわちタグ領域の位置を移動させる。このとき、例えば、収集部は、タグ領域と分割領域との相対位置が変わらないように、タグ領域の絶対位置を変える。
また、例えば、アディアバティックパルスを用いた撮像が行われる場合に、分割された範囲ごとに、アディアバティックパルスである反転パルスのRFパワー(フリップ角)を変えてもよい。アディアバティックパルスは、理論的には、縦磁化が180度倒れる以上のRFパワー(フリップ角)に対しても縦磁化は180度以上回転することはない。しかし、RFパワーが大きくなれば、送信RFコイルに送信されるRF波を増幅するRFアンプの負荷が大きくなるので、短時間に複数のアディアバティックパルスを印加する場合は、デューティサイクルを考慮した最適な大きさにRFパワーを配分する必要がある。例えば、送信RFコイルの感度が不均一であることを考慮して、収集部が、分割された範囲のうち、送信RFパワーが低めになるRFコイルの中心から遠い範囲については、反転パルスのフリップ角を300度とし、送信RFパワーが十分なRFコイルの中心に近い範囲については、反転パルスのフリップ角を180度とするように制御する。
上述した少なくとも1つの実施形態に係るMRI装置によれば、撮像対象の部位や撮像の目的に応じて、より高画質な画像を得ることができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
100 MRI(磁気共鳴イメージング)装置
30 計算機システム
36 制御部
36b 分割部
36c 収集部
36e 合成部

Claims (2)

  1. 被検体の撮像領域に流入する流体を標識化するためのRF波を印加した後に、待ち時間の異なる複数のタイミングでデータ収集を行う場合に、前記待ち時間に基づいて、前記複数のタイミングを少なくとも2つの時間的な範囲に分割する分割部と、
    前記撮像領域の全体について、第1の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第1の読み出しシーケンスと異なる第2の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
    前記第1の読み出しシーケンスによって収集されたデータから生成された画像と、前記第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
    を備える、磁気共鳴イメージング装置。
  2. 被検体の撮像領域に流入する流体を標識化するためのRF波を印加して、待ち時間が経過した後にデータ収集を行うパターンを、前記待ち時間を異ならせて複数回実行する場合に、前記待ち時間に基づいて、前記待ち時間に応じた複数のタイミングを少なくとも2つの時間的な範囲に分割する分割部と、
    前記撮像領域の全体について、第1の読み出しシーケンス、及び、シーケンスの種類及び撮像条件の少なくとも一方が前記第1の読み出しシーケンスと異なる第2の読み出しシーケンスそれぞれを用いてデータ収集を行う収集部と、
    前記第1の読み出しシーケンスによって収集されたデータから生成された画像と、前記第2の読み出しシーケンスによって収集されたデータから生成された画像とを合成する合成部と
    を備える、磁気共鳴イメージング装置。
JP2019124132A 2013-04-04 2019-07-03 磁気共鳴イメージング装置 Active JP6786670B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013078873 2013-04-04
JP2013078873 2013-04-04

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2014077341A Division JP6552789B2 (ja) 2013-04-04 2014-04-03 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
JP2019188205A true JP2019188205A (ja) 2019-10-31
JP6786670B2 JP6786670B2 (ja) 2020-11-18

Family

ID=51658455

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2014077341A Active JP6552789B2 (ja) 2013-04-04 2014-04-03 磁気共鳴イメージング装置
JP2019124132A Active JP6786670B2 (ja) 2013-04-04 2019-07-03 磁気共鳴イメージング装置

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2014077341A Active JP6552789B2 (ja) 2013-04-04 2014-04-03 磁気共鳴イメージング装置

Country Status (4)

Country Link
US (1) US10663552B2 (ja)
JP (2) JP6552789B2 (ja)
CN (1) CN105072991B (ja)
WO (1) WO2014163164A1 (ja)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6552789B2 (ja) * 2013-04-04 2019-07-31 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
WO2014168033A1 (ja) * 2013-04-08 2014-10-16 国立大学法人京都大学 磁気共鳴イメージング装置
CN106415298B (zh) * 2014-04-18 2020-01-31 皇家飞利浦有限公司 涉及对血管造影加权图像和灌注加权图像的采集的mri
US10460452B2 (en) 2014-11-14 2019-10-29 The Regents Of The University Of California Ultrasound-based volumetric particle tracking method
US9962142B2 (en) * 2014-11-14 2018-05-08 The Regents Of The University Of California Ultrasound-based volumetric particle tracking method
DE102014223734B4 (de) * 2014-11-20 2018-08-09 Siemens Healthcare Gmbh Korrektur von bildgebenden Verfahren in einer Magnetresonanzvorrichtung
US10524686B2 (en) * 2014-12-01 2020-01-07 The Regents Of The University Of California Diffusion reproducibility evaluation and measurement (DREAM)-MRI imaging methods
DE102014225846B4 (de) * 2014-12-15 2016-07-28 Siemens Healthcare Gmbh Ermittlung von Magnetresonanz-Angiographiebildern mit Time-of-Flight-Angiographie und Magnetresonanzeinrichtung
CN107949325B (zh) * 2014-12-26 2022-02-15 东芝医疗系统株式会社 磁共振成像装置、扩散加权图像的生成方法以及图像处理装置
JP6483269B2 (ja) * 2015-08-31 2019-03-13 株式会社日立製作所 磁気共鳴イメージング装置および撮像シーケンス生成方法
JP6580938B2 (ja) 2015-10-19 2019-09-25 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
US10058393B2 (en) 2015-10-21 2018-08-28 P Tech, Llc Systems and methods for navigation and visualization
CN109477879B (zh) * 2016-04-21 2021-09-21 皇家飞利浦有限公司 动脉结构的磁共振成像
JP2019514475A (ja) * 2016-04-21 2019-06-06 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 動脈構造の磁気共鳴画像診断
JP6877933B2 (ja) * 2016-09-23 2021-05-26 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び画像解析装置
JP7020930B2 (ja) * 2018-01-24 2022-02-16 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、磁気共鳴イメージングシステム及びパラメータ推定方法
JP7187206B2 (ja) * 2018-08-06 2022-12-12 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
DE102018214325A1 (de) * 2018-08-24 2020-02-27 Siemens Healthcare Gmbh Verfahren und Bereitstellungseinheit zum Bereitstellen eines virtuellen tomographischen Schlaganfall-Nachfolgeuntersuchungsbildes
WO2020235505A1 (ja) * 2019-05-17 2020-11-26 国立大学法人新潟大学 核磁気共鳴イメージング装置、核磁気共鳴イメージング方法、及びプログラム
US10871537B1 (en) * 2019-06-19 2020-12-22 GE Precision Healthcare LLC Systems and methods for background suppression in time-of-flight magnetic resonance angiography
JP7479959B2 (ja) * 2020-06-25 2024-05-09 富士フイルムヘルスケア株式会社 画像処理装置、医用撮像装置および画像処理プログラム
JP2022111589A (ja) * 2021-01-20 2022-08-01 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
CN116990736B (zh) * 2023-09-22 2023-12-12 山东奥新医疗科技有限公司 一种核磁共振成像方法、装置、设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10323336A (ja) * 1997-03-28 1998-12-08 Toshiba Corp Mri装置およびmr撮像方法
JP2000041967A (ja) * 1998-07-27 2000-02-15 Ge Yokogawa Medical Systems Ltd 磁気共鳴撮像方法および装置
JP2004024637A (ja) * 2002-06-27 2004-01-29 Toshiba Corp Mri装置およびmri画像撮影方法
US20080071166A1 (en) * 2006-09-06 2008-03-20 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
US20100030062A1 (en) * 2008-03-31 2010-02-04 Bolar Divya S System and method to analyze blood parameters using magnetic resonance imaging
US20110130668A1 (en) * 2009-11-27 2011-06-02 Shigeharu Ohyu Blood flow perfusion analyzing apparatus, blood flow perfusion analyzing method, fluid analyzing apparatus and fluid analyzing method
WO2012060192A1 (ja) * 2010-11-05 2012-05-10 株式会社 日立メディコ 磁気共鳴撮像装置、照射磁場計測方法

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62253043A (ja) * 1986-04-25 1987-11-04 富士電機株式会社 核磁気共鳴による断層画像診断装置
JPH05123314A (ja) * 1991-10-31 1993-05-21 Hitachi Medical Corp 磁気共鳴イメージング装置におけるマルチスライス撮像方法
JP4297541B2 (ja) * 1999-02-22 2009-07-15 株式会社日立メディコ 磁気共鳴イメージング装置
US6564080B1 (en) 1999-03-31 2003-05-13 Kabushiki Kaisha Toshiba MR imaging on ASL technique
JP4594482B2 (ja) * 2000-03-14 2010-12-08 株式会社東芝 磁気共鳴イメージング方法及び磁気共鳴イメージング装置
JP2002253527A (ja) * 2001-03-02 2002-09-10 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006304818A (ja) * 2005-04-26 2006-11-09 Hitachi Medical Corp 画像合成方法及び磁気共鳴イメージング装置
WO2006127687A2 (en) * 2005-05-20 2006-11-30 Advanced Mri Technologies Llc Efficient visualization of cortical vascular territories and functional mri by cycled arterial spin labeling mri
CN101138497B (zh) * 2006-09-06 2012-07-18 株式会社东芝 磁共振成像装置
US8159221B2 (en) * 2007-11-22 2012-04-17 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and method using SSFP having non-zero first moment gradients within the repetition time
US7898254B2 (en) * 2008-02-19 2011-03-01 Advanced Mri Technologies, Llc Arterial spin labeled, segmented, interleaved 3D GRASE MRI
US8664954B2 (en) * 2009-03-31 2014-03-04 Max-Planck-Gesellschaft Zur Foerderung Der Wissenschaften E.V. Magnetic resonance imaging with improved imaging contrast
DE102009019895B4 (de) * 2009-05-04 2011-05-12 Siemens Aktiengesellschaft Verfahren und Magnetresonanzanlage zur diffusionsgewichteten Aufnahme von MR-Signalen
JP5683987B2 (ja) 2010-02-12 2015-03-11 株式会社東芝 磁気共鳴イメージング装置
JP5689766B2 (ja) * 2011-08-04 2015-03-25 株式会社東芝 磁気共鳴イメージング装置
US9759797B2 (en) * 2012-03-20 2017-09-12 Oxford University Innovation Limited Motion sensitized and motion suppressed imaging using dante prepared pulse trains
JP6552789B2 (ja) * 2013-04-04 2019-07-31 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10323336A (ja) * 1997-03-28 1998-12-08 Toshiba Corp Mri装置およびmr撮像方法
JP2000041967A (ja) * 1998-07-27 2000-02-15 Ge Yokogawa Medical Systems Ltd 磁気共鳴撮像方法および装置
JP2004024637A (ja) * 2002-06-27 2004-01-29 Toshiba Corp Mri装置およびmri画像撮影方法
US20080071166A1 (en) * 2006-09-06 2008-03-20 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2008086748A (ja) * 2006-09-06 2008-04-17 Toshiba Corp 磁気共鳴イメージング装置
US20100030062A1 (en) * 2008-03-31 2010-02-04 Bolar Divya S System and method to analyze blood parameters using magnetic resonance imaging
US20110130668A1 (en) * 2009-11-27 2011-06-02 Shigeharu Ohyu Blood flow perfusion analyzing apparatus, blood flow perfusion analyzing method, fluid analyzing apparatus and fluid analyzing method
JP2011131041A (ja) * 2009-11-27 2011-07-07 Toshiba Corp 血流動態解析装置、血流動態解析プログラム、流体解析装置及び流体解析プログラム
WO2012060192A1 (ja) * 2010-11-05 2012-05-10 株式会社 日立メディコ 磁気共鳴撮像装置、照射磁場計測方法
US20130207653A1 (en) * 2010-11-05 2013-08-15 Kosuke Ito Magnetic resonanace imaging apparatus and method for measuring irradiation magnetic field

Also Published As

Publication number Publication date
JP2014210175A (ja) 2014-11-13
US20160018501A1 (en) 2016-01-21
WO2014163164A1 (ja) 2014-10-09
CN105072991A (zh) 2015-11-18
JP6786670B2 (ja) 2020-11-18
JP6552789B2 (ja) 2019-07-31
US10663552B2 (en) 2020-05-26
CN105072991B (zh) 2019-12-17

Similar Documents

Publication Publication Date Title
JP6786670B2 (ja) 磁気共鳴イメージング装置
US9213075B2 (en) Apparatus and method for magnetic resonance imaging
JP4342809B2 (ja) Mri装置
US9517033B2 (en) Magnetic resonance imaging apparatus
USRE46567E1 (en) Magnetic resonance imaging apparatus
KR101664433B1 (ko) 자기 공명 영상 장치 및 자기 공명 영상 장치의 영상 생성 방법
US9123121B2 (en) Magnetic resonance imaging apparatus and fluid imaging method
Edelman et al. Quiescent‐inflow single‐shot magnetic resonance angiography using a highly undersampled radial k‐space trajectory
US20080064951A1 (en) Magnetic resonance imaging system
CN110095742A (zh) 一种基于神经网络的平面回波成像方法和装置
JP2015196096A (ja) 磁気共鳴イメージング装置、方法及びプログラム
JP5777425B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP5689766B2 (ja) 磁気共鳴イメージング装置
JP2017153747A (ja) 磁気共鳴イメージング装置
JP2013070765A (ja) 磁気共鳴イメージング装置
US10444316B2 (en) Reduction of eddy currents during flow encoded magnetic resonance imaging
JP2012254361A (ja) 磁気共鳴イメージング装置
US10871538B2 (en) Planning support for selective arterial spin labeling MR imaging methods
JP2012096092A (ja) 磁気共鳴イメージング装置
JP6449374B2 (ja) 磁気共鳴イメージング装置
JP2016040024A (ja) 磁気共鳴イメージング装置
Menza Accelerated, high spatial and temporal resolution phase contrast techniques for functional analysis of the myocardium
Thomas et al. MR angiography and arterial spin labelling
JP5658312B2 (ja) 磁気共鳴イメージング装置
Mora Álvarez Arterial Spin Labeling (ASL) Imaging of Cerebral Perfusion at High Magnetic Fields

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190703

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200609

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200728

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200929

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201028

R150 Certificate of patent or registration of utility model

Ref document number: 6786670

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150