JPWO2010064572A1 - 磁気共鳴イメージング装置および同期撮像方法 - Google Patents

磁気共鳴イメージング装置および同期撮像方法 Download PDF

Info

Publication number
JPWO2010064572A1
JPWO2010064572A1 JP2010541301A JP2010541301A JPWO2010064572A1 JP WO2010064572 A1 JPWO2010064572 A1 JP WO2010064572A1 JP 2010541301 A JP2010541301 A JP 2010541301A JP 2010541301 A JP2010541301 A JP 2010541301A JP WO2010064572 A1 JPWO2010064572 A1 JP WO2010064572A1
Authority
JP
Japan
Prior art keywords
echo signal
measurement
period
magnetic resonance
imaging apparatus
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
JP2010541301A
Other languages
English (en)
Other versions
JP5372015B2 (ja
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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical 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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2010541301A priority Critical patent/JP5372015B2/ja
Publication of JPWO2010064572A1 publication Critical patent/JPWO2010064572A1/ja
Application granted granted Critical
Publication of JP5372015B2 publication Critical patent/JP5372015B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7285—Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
    • 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
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

Landscapes

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

Abstract

被検者の周期的体動情報に同期させた撮像において、所望の画像コントラストを維持して、撮像時間を短縮する。そのために、周期的体動を有する被検者の該周期的体動情報から検出したトリガー情報に同期したエコー信号の同期計測において、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、K空間を複数の部分領域に分割し、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。

Description

本発明は、被検者中の水素や燐等からの核磁気共鳴(以下、「NMR」という)信号を測定し、核の密度分布や緩和時間分布等を画像化する磁気共鳴イメージング装置および同期計測方法に関し、特に同期計測の撮像効率向上に関するものである。
核磁気共鳴イメージング(以下、「MRI」という)装置は、被検者の組織を構成する原子核スピンが発生するNMR信号(エコー信号)を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮像においては、エコー信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたエコー信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。
上記MRI装置における撮像では、一般的に数分から数十分もの撮像時間が必要となるため、撮像中に被検者の心拍や呼吸などの体動を避けることができない。このため、画像上に体動に由来するアーチファクトが発生し画質が劣化することが知られている。
MRI装置における撮像では、上記問題点を回避する方法として、特許文献1〜3に開示されているような、被検者に心電極や脈波センサーなどを装着して生体信号を検出し、検出された生体信号をトリガー信号として、エコー信号を収集するタイミングを心臓の動きなどに同期させて撮像する方法が用いられる。つまり、特許文献1〜3では、トリガー信号に同期して、被検者の比較的小さな動きの時相のみ選択的にエコー信号の計測を行うことにより体動に由来する画像上のアーチファクトを良好に抑制している。
特開2008-125986号公報 特開2008-136851号公報 特許第4090619号公報
特許文献1や特許文献2では、被検者の比較的小さな動きの時相のみ選択的にエコー信号の計測を行うために、撮像時間が延長するという課題が残されている。特に特許文献3で述べられている心電同期計測方法では、拡張期に所定スライスエンコード量分のエコー信号を収集する動作を一定の複数心拍毎に繰り返す3次元スキャンを実行している。これにより、乱流的な血流の影響を排除して、安定した血流状態でのエコー信号を計測し、再構成画像のコントラストを高めているが、一定の複数心拍毎にパルスシーケンスを繰り返すために、撮像時間の延長という課題は未解決のまま残ってしまう。
そこで本発明の目的は、MRI装置を用いた、被検者の周期的体動情報に同期させた撮像において、所望の画像コントラストを維持して、撮像時間を短縮することにある。
上記の目的を達成するために、本発明は、周期的体動を有する被検者の該周期的体動情報から検出したトリガー情報に同期したエコー信号の同期計測において、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、K空間を複数の部分領域に分割し、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
具体的には、本発明のMRI装置は、被検体の周期的体動情報からトリガー情報を検出する検出部と、トリガー情報に同期して被検体からのエコー信号を計測する同期計測を制御する計測制御部と、エコー信号のデータがK空間に配置されて成るK空間データに基づいて、被検体の画像を取得する演算処理部と、を備え、同期計測は、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、演算処理部は、K空間を複数の部分領域に分割し、計測制御部は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
また、本発明の同期撮像方法は、被検体の周期的体動情報からトリガー情報を検出する検出工程と、トリガー情報に同期して被検体からのエコー信号を計測する同期計測を制御する計測制工程と、エコー信号のデータが配置されるK空間を複数の部分領域に分割する工程と、エコー信号を用いて、被検体の画像を取得する演算処理工程と、を備え、同期計測は、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、計測制御工程は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
本発明のMRI装置及び同期撮像方法によれば、被検者の周期的体動情報に同期させた撮像において、所望の画像コントラストを維持して、撮像時間を短縮することが可能になる。
本発明に係るMRI装置の一実施例の全体基本構成を示すブロック図。 本発明に係る3次元K空間の分割の一例を示す図。 本発明に係る3次元K空間の分割設定用GUIの一例を示す図。 本発明に係る3次元K空間の分割設定用GUIの他の一例を示す図。 本発明の第1の実施形態に係る心電同期計測の一例を示す図。(a)は、トリガー信号としての心電波形のR波と、エコー信号計測のタイムテーブルをそれぞれ示す図。(b)は、エコー信号計測のためのパルスシーケンスのシーケンスチャートを示す図。 本発明の第1の実施形態の処理フォローを示すフォローチャート。 1心拍動周期に計測するエコーデータ数(N)の決定処理フローを示すフローチャート。 本発明の第2の実施形態に係る心電同期計測の一例を示す図。(a)は、T2強調画像を取得する場合のトリガー信号とエコー信号計測のタイムテーブルをそれぞれ示す図(b)は、エコー信号計測のためのパルスシーケンスのシーケンスチャートを示す図。 本発明の第2の実施形態の処理フォローを示すフォローチャート。 部分領域に応じたトリガー信号待機回数(NT)の設定処理フローを示すフォローチャート。
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検者の断層画像を得るもので、図1に示すように、MRI装置は静磁場発生部2と、傾斜磁場発生部3と、送信部5と、受信部6と、情報処理部7と、計測制御部4と、を備えて構成される。
静磁場発生部2は、垂直磁場方式であれば、被検者1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させるもので、被検者1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
傾斜磁場発生部3は、MRI装置の座標系(静止座標系)であるX，Y，Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成る。後述の計測制御部4からの命令に従ってそれぞれのコイルの傾斜磁場電源10が駆動されることにより、X，Y，Zの3軸方向の傾斜磁場Gx，Gy，Gzが、被検者1が横たわる静磁場空間に印加される。撮像時には、スライス面(撮像断面)に直交する方向にスライス方向傾斜磁場パルス(Gs)が印加されて被検者1に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード方向傾斜磁場パルス(Gp)と周波数エンコード方向傾斜磁場パルス(Gf)とが印加されて、エコー信号にそれぞれの方向の位置情報がエンコードされる。
送信部5は、被検者1の生体組織を構成する原子の原子核スピンにNMR現象を誘起するために、被検者1に高周波磁場パルス(以下、「RFパルス」という)を照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとから成る。高周波発振器11から出力された高周波パルスが、計測制御部4からの指令によるタイミングで変調器12により振幅変調され、この振幅変調された高周波パルスが高周波増幅器13で増幅された後に、被検者1に近接して配置された高周波コイル14aに供給されることにより、RFパルスが被検者1に照射される。
受信部6は、被検者1の生体組織を構成する原子核スピンのNMR現象により放出されるエコー信号を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射されたRFパルスによって誘起された被検者1の応答のエコー信号が被検者1に近接して配置された高周波コイル14bで検出され、信号増幅器15で増幅された後、計測制御部4からの指令によるタイミングで、直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、エコー信号として演算処理部7に送られる。以下デジタル量に変換されたエコー信号をエコーデータという。
計測制御部4は、ある所定のパルスシーケンスに基づいて、傾斜磁場発生部3、送信部5、および受信部6を制御して、RFパルスと傾斜磁場パルスの印加と、エコー信号の計測とを、繰り返す制御手段である。計測制御部4は、CPU8の制御で動作し、被検者1の断層画像の再構成に必要なエコー信号収集に必要な種々の命令を傾斜磁場発生部3、送信部5、および受信部6に送って、これらを制御する。
情報処理部7は、各種データ処理と処理結果の表示及び保存等を行うもので、CPU(演算処理部)8と、光ディスク19や磁気ディスク18等の外部記憶装置と、ディスプレイ20と、から成る。受信部6からのエコーデータがCPU8に入力されると、CPU8内のK空間に対応するメモリに、このエコーデータが記憶される(以下、エコー信号またはエコーデータをK空間に配置する旨の記載は、エコーデータがこのメモリに書き込まれて記憶されることを意味する。そして、K空間に配置されたエコーデータをK空間データという)。そして、CPU8はこのK空間データに対して、信号処理、画像再構成等の演算処理を実行し、その結果である被検者1の断層画像をディスプレイ20に表示すると共に、外部記憶装置に記録する。
操作部25は、操作者からの、MRI装置の各種制御情報や上記演算処理部7で行う処理の制御情報の入力を受け付け、トラックボール又はマウス23、及び、キーボード24から成る。この操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。
本発明に係るMRI装置は、さらに、被検者に装着されて被検者から心電波形信号を取得する心電極31と、心電極からの信号が入力されて、被検者の心電波形及びそのR波(トリガー信号)を検出する心電波形モニタ32とを備える。心電波形モニタ32が検出した心電波形情報(周期的体動情報の一例)は、計測制御部4を介してCPU8に入力され、計測制御部4がトリガー信号(トリガー情報)に同期して、所定のパルスシーケンスに基づいて上記各部を制御することにより同期計測を制御する。
なお、図1において、送信側の高周波コイル14aと傾斜磁場コイル9は、被検者1が挿入される静磁場発生部2の静磁場空間内に、垂直磁場方式であれば被検者1に対向して、水平磁場方式であれば被検者1を取り囲むようにして設置されている。また、受信側の高周波コイル14bは、被検者1に対向して、或いは取り囲むように設置されている。
現在のMRI装置の撮像対象核種は、臨床で普及しているものとしては、被検者の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
(K空間の分割)
最初に、本発明に係るK空間の分割について説明する。本発明のMRI装置及び同期撮像方法は、被検体の周期的体動情報から検出したトリガー情報に同期した同期計測において、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設ける。そして、K空間を複数の部分領域に分割し、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
K空間が2次元である場合には、K空間の分割を位相エンコード方向に行い、K空間が3次元である場合には、K空間の分割をスライスエンコード方向と位相エンコード方向の内の少なくとも一方向に行う。K空間の分割の一例を図2に示す。図2は、3次元K空間200を、スライスエンコード方向(kz)に、原点を通るSlice(kz)軸に垂直な平面に関して正負対称に、正負それぞれの側を3つの領域(領域201、202、203)に分割する例を示す。具体的には、CPU8は、3次元K空間を、原点を含む低空間周波数側の部分領域201と、その両側の高空間周波数側をそれぞれ2つの部分領域202-1及び203-1と、202-2及び203-2 (ハイフン"-"以降の数字は、1がK空間スライスエンコード方向(Kz)の正側、2が負側を意味する)とに分割する。さらに、分割の割合は操作者の設定に応じて変更可能とする。図2は、部分領域201が3次元K空間全体の50%と設定され、部分領域202と部分領域203はそれぞれ3次元K空間全体の25%ずつとなるように設定されて3次元K空間が分割された例を示す。なお、さらに位相エンコード方向にもK空間を分割してもよい。
各分割領域の3次元K空間全体に対する割合は、操作者が変更可能な撮像パラメータとする。この各分割領域のK空間全体に対する割合の設定は、例えば図3に示すようなGUI300を介して操作者が行うことを可能とする。図3に示す例は、操作者がGUI300上の3次元K空間の分割数Segment#(301)を入力設定した後、CPU8は各部分領域の全体に対する割合(Data Rate)をそれぞれ決定するためのタグをディスプレイ20に表示させ、操作者が部分領域毎の3次元K空間全体に対する割合タグ(302,303,304)に所望の数値を入力設定する。これらの操作者の設定入力に基づいて、CPU8はそれぞれの部分領域の割合を決定すると共に、この設定情報を計測制御部4に通知する。
或いは、図4に示すように、3次元K空間の分割数Segment#301と、K空間の低空間周波数側の部分領域の、K空間全体に対する割合Data Rate302のみを操作者に入力設定してもらい、CPU8は、残りの部分領域のK空間全体に対する割合を均等に分割するものとして、残りの部分領域のK空間全体に対する割合を算出して決定するとともに、各部分領域の割合を計測制御部4に通知することも可能である。
以下、体動として心臓の拍動を、周期的体動情報として心電波形を想定し、心電波形のR波をトリガー信号(トリガー情報の一例)として、このトリガー信号に同期してエコー信号の計測を制御する心電同期計測を例にして、本発明の各実施形態を説明する。しかし、本発明は、心電同期計測に限らず、血流の脈動や関節の周期的運動等の他の周期的体動にも適用可能である。
(第1の実施形態)
次に、本発明のMRI装置及び同期撮像方法の第1の実施形態を説明する。本実施形態は、第1の期間として、トリガー情報とエコー信号の計測期間との間にディレイ時間を設定し、高空間周波数側の部分領域に対応するエコー信号の計測におけるディレイ時間を、低空間周波側の部分領域に対応するエコー信号の計測におけるディレイ時間よりも短くする。本実施形態は、例えばＴ1強調画像の取得に好適である。以下、図5〜図7に基づいて本実施形態を説明する。
最初に本実施形態の概要を、図5を用いて説明する。図5は、被検者の心電波形を検出しそのR波をトリガー信号としてこれに同期して行う、本実施形態に係る心電同期計測の一例を示す。図5(a)は、トリガー信号としての心電波形のR波と、エコー信号計測のタイムテーブルをそれぞれ示している。また、図5(b)は、エコー信号計測のためのパルスシーケンスの一例のシーケンスチャートを示している。
操作者により図5に示す心電同期計測が起動されると、計測制御部4は、トリガー信号501に同期してパルスシーケンスを駆動する。計測制御部4は、このパルスシーケンスを短い繰り返し時間(TR)で繰り返してエコー信号を計測(503)するが、部分領域201のエコー信号計測中では、トリガー信号から一定時間(Delay time)経過まではエコー信号収集は行わず空打ち502を行う。この空打ち502は、同じパルスシーケンスであるがエコー信号を計測しないか、計測しても画像再構成に用いない。このようにトリガー信号から一定のDelay timeを設けて空打ちを行う理由は、Delay timeの時間帯は心臓の収縮期で心臓の動きが速いためであり、心臓の動きが比較的ゆっくりしている拡張期にエコー信号の計測を行うためである。これにより、心臓の動きに起因して画像上に発生するアーチファクトを低減して高画質の画像を取得することが可能になる。
計測制御部4は、トリガー信号から一定時間(Delay time)経過後、被検者の体動が少ない時間帯に、短い繰り返し時間(TR)で位相エンコードを変えてパルスシーケンスを少なくとも1回繰り返して、部分領域201に対応するエコー信号を少なくとも1以上計測する。そして、計測制御部4は、このようなDelay time後に短い繰り返し時間(TR)でパルスシーケンスを繰り返して行うエコー信号の計測を、1心拍動周期毎に複数心拍繰り返し実行して、部分領域201に対応するエコー信号を全て計測する。
次に、部分領域202のエコー信号の計測へ移行すると、計測制御部4はDelay Timeを短く設定し、1心拍動周期に計測するエコー信号の数を増加させる。図5では、部分領域201のエコー信号計測時と比較して、部分領域202のエコー信号計測時のDelay Timeが、パルシーケンスの繰り返し時間(TR)の一繰り返し時間分(1TR)短くなっているが、Delay Timeの短縮は、1TRに限定されるものではなく、また、パルシーケンスの繰り返し時間(TR)単位に限定されるものでもなく、任意に設定することが可能である。
同様にして、部分領域203のエコー信号計測へ移行した段階で、計測制御部4は部分領域202のエコー信号計測時と比較して、部分領域203のエコー信号計測時のDelay Timeをさらに短く設定し、1心拍動周期に計測するエコー信号の数を増加させる。図5の例は、Delay Time=0(ゼロ)と設定した場合を示す。すなわち、計測制御部4はトリガー信号直後から空打ち無しにエコー信号の計測を行う。
一般的に、低空間周波数領域のエコー信号は、その強度が大きく、画像のコントラストを含む画質を支配する重要な信号である。このような低空間周波数領域のエコー信号を、心臓の動きが比較的ゆっくりしている拡張期に計測できるようDelay Timeを設定することにより、動きに起因するアーチファクトを低減して、所望のコントラストを有する画像を取得できるようになる。
一方、高空間周波数領域のエコー信号は、低空間周波数領域のエコー信号と異なり、画像のコントラストに影響するものではく、画像の分解能に寄与する。また、高空間周波数領域のエコー信号は強度も小さいので画像への影響が少ない。そのため、本実施形態のように、部分領域202と203のエコー信号を、Delay Timeを少なくして収縮期にかかる時間帯で計測して、これらのエコー信号を画像再構成に用いても画像への影響が軽微となる。そのため、本実施形態の様に、高空間周波数領域のエコー信号ほど、Delay Timeを少なくして計測して良いことになる。逆に、高空間周波数領域のエコー信号計測の際のDelay Timeを短くすることによって、計測できるエコー信号数を増加できるので、全体の撮像時間を短縮できる。
また、本実施形態は、前述の通り、1心拍動周期内で短い繰り返し時間(TR)でパルスシーケンスを繰り返してエコー信号を計測する。そのため、本実施形態はT1強調画像の取得に好適となる。そこで、T1強調画像の取得に好適なパルスシーケンスとして、例えば図5(b)に示すスピンエコー(SE)シーケンスを用いる。図5(b)に示すシーケンスチャートのRF、Gs、Gp、Gf及びEchoは、それぞれRFパルス、スライス傾斜磁場、位相エンコード傾斜磁場、周波数エンコード傾斜磁場、及びエコー信号を表す。スライス選択傾斜磁場512を被検者に印加した状態で、90度RFパルス511を被検者に印加して所望の領域の磁化を90度励起して横磁化を発生させる。その直後に、スライスリフェーズ傾斜磁場513を印加して、所望領域の励起に伴う横磁化の位相分散を再収束させる。次に、スライスエンコード傾斜磁場514を印加して、エコー信号の位相にスライスエンコード方向の空間情報をエンコードする。また、対称なエコー信号を発生させるために周波数ディフェーズ傾斜磁場515を印加して、横磁化の位相を分散させる。その後に、スライス傾斜磁場517と180度RFパルス516を印加して、横磁化を180度反転させて再収束させる。そして、位相エンコード傾斜磁場518を印加して、エコー信号の位相に位相エンコード方向の空間情報をエンコードし、周波数エンコード傾斜磁場519を印加しながらエコー信号520を計測する。これによりエコー信号520に周波数エンコード方向の空間情報がエンコードされる。以上の各パルスを、スライスエンコード傾斜磁場514と位相エンコード傾斜磁場518の少なくとも一方を変えながら短い繰り返し時間(TR)で繰り返して、各部分領域のエコー信号を計測する。
また、本実施形態は、前述のDelay Timeを利用して空打ちを行うことにより、エコー信号計測前に励起領域からのエコー信号の強度を安定化させる。特に低空間周波数領域に対応するエコー信号の計測の際には、積極的にDelay Timeを利用して空打ちを行うことにより、強度の安定したエコー信号を計測して高画質の画像を取得する。
以上迄が本実施形態の概要である。次に本実施形態の動作フローを、図6,図7に示すフローチャートを用いて説明する。図6は、本実施形態の動作フローの概要を、図7は特にステップ602の処理の詳細を示す。これらの処理フローは、プログラムとして外部記憶装置に記憶されており、必要に応じてCPU8がそのメモリに読み込んで実行することにより、実行される。なお、図6,図7は、例えば図3，図4に示したGUIを介して事前に操作者が3次元K空間の分割数を入力設定し、CPU8がその入力された設定値に応じて、K空間の分割数を決定した後の処理フォローを表している。
ステップ601で、操作者により心電同期計測が起動されると、計測制御部4は、被検体の心電波形におけるR波を、同期計測を行うためのトリガー信号として、その読み込みを開始する。
ステップ602で、3次元K空間の部分領域に応じた、1心拍動周期で計測するエコー信号数(N)やDelay time(Td)等が決定される。この決定処理の詳細は図7に基づいて後述する。なお、ステップ601と602の実行順序はどちらが先でも良い。
ステップ603で、計測制御部4は、トリガー信号を待つ。
ステップ604で、トリガー信号が受信されると、計測制御部4は、ステップ602で決定されたDelay time(Td)だけ待機する処理を行う。このDelay time(Td)期間中には、計測制御部4は、前述の通り空打ち605を行う。
ステップ606で、Delay time(Td)が経過した時点が最も被検者の動きが少ないタイミングであるとして、図5に示すような、1心拍動周期のエコー信号計測が実行される。つまり、計測制御部4は、スライスエンコード傾斜磁場もしくは位相エンコード傾斜磁場の印加量を変えて、ステップ602で設定された、1心拍動周期に計測するエコー信号数分の計測を行う。
ステップ607で、1心拍動周期のエコー信号計測が終了し、且つ心電同期計測が終了でない場合はステップ601に戻って、再度トリガー信号の読み込みの開始と、K空間の部分領域に応じた1心拍動周期に収集するエコー信号数の決定を行い、撮像終了となるまでステップ601〜606を繰り返し実行する。
以上迄が、本実施形態の処理フローの概要である。次にステップ602の1心拍動周期に計測するエコー信号数の決定処理フローを、図7に示すフローチャートに基づいて説明する。
ステップ701で、CPU8は、部分領域201のエコー信号計測時におけるDelay Time(Td1)と計測するエコー信号数(N1)を基準として、部分領域202のエコー信号計測時における、Delay Time(Td2)、計測エコー信号数(N2)と、部分領域203のエコー信号計測時における、Delay Time(Td3)、計測エコー信号数(N3)と、を求める。なお、計測エコー信号数（N1)とDelay Time(Td1)は、例えば、操作者が計測エコー信号数（N1)とDelay time(Td1)を設定する、操作者が計測エコー信号数（N1)を設定した後、CPU8がDelay Time(Td1)を算出する、操作者がDelay time(Td1)を設定した後、CPU8が計測エコー信号数（N1)を算出する、のいずれでも良い。分割領域202、分割領域203でのDelay Time、計測エコー信号数は、例えば、撮像パラメータとして操作者により設定される撮像高速化率Rapid Rateを用いて下式で求めることが出来る。
N2=N1×(1＋Rapid Rate/100), Td2=Td1＋TR×(N2-N1) (1)
N3=N2×(1＋Rapid Rate/100), Td3=Td2＋TR×(N3-N2)
撮像高速化率Rapid Rateは、例えば、最終的にDelay Timeが最小値となるように、つまり、最後の部分領域(この場合では部分領域203)のエコー信号計測時のDelay Timeが最小値となるように、自動的に算出することも可能である。
ステップ702で、CPU8は、部分領域202のエコー信号の計測であるか否かを判断し、そうであれば、ステップ703で、Delay Timeと計測エコー信号数とを、(1)式で求めた部分領域202の値Td2,N2にそれぞれ設定し、計測制御部4に通知する。
ステップ704で、CPU8は、部分領域203のエコー信号の計測であるか否かを判断し、そうであれば、ステップ705で、Delay Timeと計測エコー信号数とを、(1)式で求めた部分領域203の値Td3,N3にそれぞれ設定し、計測制御部4に通知する。
以上迄が、ステップ602の1心拍動周期に計測するエコー信号数の決定処理フローの説明である。
以上説明したように、本実施形態のMRI装置及び同期撮像方法によれば、3次元K空間をスライスエンコード方向に複数分割し、分割された部分領域のスライスエンコード方向の位置に応じて、エコー信号を計測する際のDelay Timeを異ならせることにより、1心拍動周期に計測するエコー信号数を増加させて撮像時間の短縮化を行うことが可能となる。具体的には、低空間周波数側の部分領域のエコー信号の計測の際のDelay Timeを所望の画像コントラストとなるように設定することにより3次元画像のコントラストを所望の状態に維持できる。そして、高空間周波数側の部分領域のエコー信号の計測の際のDelay Timeを低空間周波数側の部分領域のエコー信号の計測の際のDelay Timeよりも短くすることにより、高空間周波数側の部分領域において1心拍動周期に計測するエコー信号数を増加させて撮像時間の短縮化を行うことができる。
なお、以上の本実施形態の説明では、3次元計測を例にしたが、2次元の心電同期計測においても、2次元K空間を位相エンコード方向に複数分割し、分割された部分領域の位相エンコード方向の位置に応じて、エコー信号を計測する際のDelay Timeを異ならせることにより、2次元画像のコントラストを所望の状態に維持しつつ、1心拍動周期に計測するエコー信号の数を増加させて撮像時間の短縮化を行うことが可能である。
(第2の実施形態)
次に、本発明のMRI装置及び同期撮像方法の第2の実施形態を説明する。本実施形態は、第2の期間として、一つのエコー信号計測期間から、次のエコー信号計測期間のトリガーとなるトリガー情報までの待ち時間を設定し、高空間周波数側の部分領域に対応するエコー信号の計測における待ち時間を、低空間周波数側の部分領域に対応するエコー信号の計測における待ち時間よりも短くする。好ましくは、周期的体動の一周期の整数倍だけ、待ち時間を異ならせる。すなわち、エコー信号の計測を行うためのパルスシーケンスの繰り返し時間について、高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くする。本実施形態は、例えばＴ2強調画像の取得に好適である。以下、前述の第1の実施形態と異なる点のみを図8〜図10に基づいて説明し、同じ点の説明は省略する。
最初に本実施形態の概要を、図8を用いて説明する。図8は、被検者の心電波形を検出しそのR波をトリガー信号としてこれに同期して行う、本実施形態に係る心電同期計測の一例を示す。図8(a)は、T2強調画像を取得する場合のトリガー信号とエコー信号計測のタイムテーブルをそれぞれ示している。図8(b)は、エコー信号計測のためのパルスシーケンスの一例のシーケンスチャートを示している。
計測制御部4は、トリガー信号801-1aの検出から一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号計測802-1aを行う。このパルスシーケンスには、図8(b)に示すように、複数のエコー信号を計測可能なFSEシーケンスやEPIシーケンスなどのマルチエコータイプのパルスシーケンスが用いられる。詳細は後述する。計測制御部4は、このようなパルスシーケンスに基づいて所定数のエコー信号の計測を制御する。この心拍動周期でのエコー信号計測後、T2強調画像を得るために、数(例えば2,3)心拍を待ち時間として、この待ち時間経過後、再度トリガー信号801-1bの検出後にパルスシーケンスを起動してエコー信号の計測802-1bを行う。
計測制御部4は、このような、Delay Time後のエコー信号の計測と、数心拍動周期分の待ち時間の待機と、を一単位とする処理を、複数心拍動周期毎(ここでは3心拍動周期毎)に複数回繰り返して、つまりパルスシーケンスの繰り返し時間(TR1)=3心拍動周期として、部分領域201のエコー信号を計測する。
次に、分割領域202のエコー信号計測へ移行すると、計測制御部4は、トリガー信号801-2aの検出から、部分領域201の計測時と同じ一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号の計測802-2aを行う。その後の待ち時間は、分割領域201のエコー信号計測時と比較して短く設定される。図8(a)の例では、パルスシーケンス終了後の待ち時間が、分割領域201の場合と比較して1心拍動周期分短くされている。ただし、本実施形態は、待ち時間の短縮度合いとして1心拍動周期に規定されるものではなく、任意に設定可能である。そして、この待ち時間後に、再度トリガー信号801-2bの検出後にパルスシーケンスを起動してエコー信号の計測802-2bを行う。
計測制御部4は、このような、Delay Time後のエコー信号の計測と、部分領域201のエコー信号計測時よりも短い1以上の心拍動周期分の待ち時間の待機と、を一単位とする処理を、複数心拍動周期毎(ここでは2心拍動周期毎)に複数回繰り返して、つまりパルスシーケンスの繰り返し時間(TR2)=2心拍動周期として、部分領域202のエコー信号を計測する。
同様にして、分割領域203のエコー信号計測へ移行した段階で、計測制御部4は、パルスシーケンス動作後の待ち時間を分割領域202のエコー信号計測時よりも短く設定して、エコー信号の計測周期を短縮させる。つまり、分割領域203のエコー信号計測へ移行すると、計測制御部4は、トリガー信号801-3aの検出から、部分領域201の計測時と同じ一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号の計測802-3aを行う。その後の待ち時間は、分割領域202のエコー信号計測時と比較して短く設定される。図8(a)の例では、パルスシーケンス終了後の待ち時間が、分割領域202の場合と比較して1心拍動周期分短くされて、結果として1心拍動周期以下の時間とされている。そして、この待ち時間後に、再度トリガー信号801-3bの検出後にパルスシーケンスを起動してエコー信号の計測802-3bを行う。つまり、計測制御部4は、パルスシーケンスの繰り返し時間(TR3)=1心拍動周期として、部分領域203のエコー信号を計測する。
以上の説明を纏めると、本実施形態は、エコー信号の計測を行うためのパルスシーケンスの繰り返し時間(TR)について、待ち時間を異ならせることにより、高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くする。
なお、高空間周波数側の部分領域に対応するエコー信号の計測においては、Delay time＝0(ゼロ)としてもよい。
また、本実施形態は、前述の通り、部分領域に応じて1以上の心拍動周期毎にマルチエコータイプのパルスシーケンスを繰り返してエコー信号を計測する。特に低空間周波数側の部分領域では、複数心拍動周期毎にパルスシーケンスを繰り返す。そのため、本実施形態はT2強調画像の取得に好適となる。そこで、T2強調画像の取得に好適なパルスシーケンスとして、例えば図8(b)に示すFSEシーケンスを用いる。図8(b)に示すシーケンスチャートのRF、Gs、Gp、Gf及びEchoは、図5(b)と同様である。FSEシーケンスは、図5(b)に示すスピンエコーシーケンスの180度RFパルス(516-1〜516-6)とスライス傾斜磁場(517-1〜517-5)と周波数エンコード傾斜磁場(519-1〜519-5)とを繰り返して複数のエコー信号(520-1〜520-5)を計測するパルスシーケンスであって、各エコー信号の計測の前後で、絶対量が同じで極性の異なる印加量となる一対の位相エンコード傾斜磁場(801-1a、801-1b〜801-5a、ただし801-5b、801-3は印加量ゼロ)を、エコー信号の計測毎にその印加量を変えて印加する。そのため位相エンコード傾斜磁場518は印加しない。これにより、位相エンコードの異なる複数のエコー信号(520-1〜520-5)を一回の90度RFパルス511で計測することが可能になる。
以上迄が本実施形態の概要である。次に本実施形態の動作フローを、図9, 図10に示すフローチャートを用いて説明する。図9は、本実施形態の動作フローの概要を、図10は特にステップ901の処理の詳細を示す。これらの処理フローは、プログラムとして外部記憶装置に記憶されており、必要に応じてCPU8がそのメモリに読み込んで実行することにより、実行される。なお、前述の第1の実施形態の図6,図7の説明と同様に、操作者が3次元K空間の分割数を入力設定し、CPU8がその入力された設定値に応じて、K空間の分割数を決定したものとする。
ステップ901で、操作者により心電同期計測が起動されると、撮像を行うK空間の部分領域に応じてトリガー信号待機回数(NT)の設定が行われる。この設定処理の詳細は図10に基づいて後述する。
ステップ902で、上記トリガー信号待機回数(NT)の設定完了後、計測制御部4は、トリガー信号の読み込みを開始する。
ステップ903で、計測制御部4は、トリガー信号を受信すると、ステップ904で、設定されたたトリガー信号待機回数(NT)から1を減算する。
ステップ905で、計測制御部4は、事前に設定されたDelay time(Td)だけ待機する。このDelay Ｔｉｍｅ期間中に前述の空打ちを行えば、その後に計測するエコー信号の強度が安定して高画質の画像を得ることができる。
ステップ906で、Delay time(Td)が経過した時点が最も被検者の動きが少ないタイミングであるとして、計測制御部4は、例えば図8(b)に示すようなパルスシーケンスを実行する。
ステップ907で、このイメージングスキャンが終了した後、計測制御部4は、再度トリガー信号の読み込みを行い、トリガー信号を受信すると、設定したトリガー信号待機回数(NT)から1を減算する処理ステップ908を行い、NT=0となるまで実行(909)する。
ステップ910で、上記待機時間が終了すると、撮像終了でない場合、ステップ901に戻って、再度トリガー信号待機回数(NT)の設定が行われ、撮像終了となるまで上記ステップ901〜909が繰り返される。
以上迄が、本実施形態の処理フローの概要である。次に、ステップ901のK空間の部分領域に応じたトリガー信号待機回数(NT)の設定処理フローの概要を図10に示すフローチャートに基づいて説明する。
ステップ1001で、CPU8は、分割領域201計測時におけるトリガー信号待機回数(NT1)を基準として、部分領域202のトリガー信号待機回数(NT2)と、部分領域203のトリガー信号待機回数(NT3)の設定を行う。なお、トリガー信号待機回数(NT1)は、例えば、操作者がトリガー信号待機回数(NT1)を設定する、予め所定の値を記憶しておく、のいずれでも良い。分割領域202、分割領域203でのトリガー信号待機回数は、例えば、撮像パラメータとして操作者により設定される撮像高速化率Rapid Rateを用いて下式で求めることが出来る。
NT2=NT1×(1−Rapid Rate/100) (2)
NT3=NT2×(1−Rapid Rate/100)
撮像高速化率Rapid Rateは、例えば、前述の第1の実施形態と同様に、最終的にDelay Timeが最小値となるように自動的に算出することも可能である。
ステップ1002で、CPU8は、部分領域202のエコー信号の計測であるか否かを判断し、そうであれば、ステップ1003で、トリガー信号待機回数を、(2)式で求めた部分領域202の値NT2に設定する。
ステップ1004で、CPU8は、部分領域203のエコー信号の計測であるか否かを判断し、そうであれば、ステップ1005で、トリガー信号待機回数を、(2)式で求めた部分領域203の値NT2に設定する。
以上迄が、ステップ901のK空間の部分領域に応じたトリガー信号待機回数(NT)の設定処理フローの説明である。
以上説明したように、本実施形態のMRI装置及び同期撮像方法によれば、3次元K空間をスライスエンコード方向に複数分割し、分割された部分領域のスライスエンコード方向の位置に応じて、エコー信号を計測した後の待ち時間を異ならせることにより、撮像時間の短縮化を行うことが可能となる。具体的には、低空間周波数側の部分領域のエコー信号の計測の際の待ち時間を所望の画像コントラストとなるように設定することにより3次元画像のコントラストを所望の状態に維持できる。そして、高空間周波数側の部分領域のエコー信号の計測の際の待ち時間を低空間周波数側の部分領域のエコー信号の計測の際の待ち時間よりも短くすることにより、撮像時間の短縮化を行うことが可能となる。
以上までが、本発明のMRI装置及び同期撮像方法に係る心電同期計測の各実施形態の説明である。しかし、本発明のMRI装置及び同期撮像方法は、上記各実施形態の説明で開示された内容にとどまらず、本発明の趣旨を踏まえた上で他の形態を取り得る。例えば、第1の実施形態の部分領域毎のDelay Timeの制御と第2の実施形態の部分領域毎の待ち時間の制御を組み合わせて実施してもよい。また、第1の実施形態でFSEシーケンスを用いてT2強調画像を取得してもよいし、第2の実施形態で、プロトン強調画像の取得や、スピンエコー(SE)シーケンスを用いてもよい。
1 被検者、2 静磁場発生系、3 傾斜磁場発生系、4 計測制御部、5 送信系、6 受信系、7 信号処理系、8 中央処理装置(CPU)、9 傾斜磁場コイル、10 傾斜磁場電源、11 高周波発振器、12 変調器、13 高周波増幅器、14a 高周波コイル(送信コイル)、14b 高周波コイル(受信コイル)、15 信号増幅器、16 直交位相検波器、17 A/D変換器、18 磁気ディスク、19 光ディスク、20 ディスプレイ、21 ROM、22 RAM、23 トラックボール又はマウス、24 キーボード、51 ガントリ、52 テーブル、53 筐体、54 処理装置
本発明は、被検体中の水素や燐等からの核磁気共鳴(以下、「NMR」という)信号を測定し、核の密度分布や緩和時間分布等を画像化する磁気共鳴イメージング装置および同期計測方法に関し、特に同期計測の撮像効率向上に関するものである。
核磁気共鳴イメージング(以下、「MRI」という)装置は、被検体の組織を構成する原子核スピンが発生するNMR信号(エコー信号)を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮像においては、エコー信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたエコー信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。
上記MRI装置における撮像では、一般的に数分から数十分もの撮像時間が必要となるため、撮像中に被検体の心拍や呼吸などの体動を避けることができない。このため、画像上に体動に由来するアーチファクトが発生し画質が劣化することが知られている。
MRI装置における撮像では、上記問題点を回避する方法として、特許文献1〜3に開示されているような、被検体に心電極や脈波センサーなどを装着して生体信号を検出し、検出された生体信号をトリガー信号として、エコー信号を収集するタイミングを心臓の動きなどに同期させて撮像する方法が用いられる。つまり、特許文献1〜3では、トリガー信号に同期して、被検体の比較的小さな動きの時相のみ選択的にエコー信号の計測を行うことにより体動に由来する画像上のアーチファクトを良好に抑制している。
特開2008-125986号公報 特開2008-136851号公報 特許第4090619号公報
特許文献1や特許文献2では、被検体の比較的小さな動きの時相のみ選択的にエコー信号の計測を行うために、撮像時間が延長するという課題が残されている。特に特許文献3で述べられている心電同期計測方法では、拡張期に所定スライスエンコード量分のエコー信号を収集する動作を一定の複数心拍毎に繰り返す3次元スキャンを実行している。これにより、乱流的な血流の影響を排除して、安定した血流状態でのエコー信号を計測し、再構成画像のコントラストを高めているが、一定の複数心拍毎にパルスシーケンスを繰り返すために、撮像時間の延長という課題は未解決のまま残ってしまう。
そこで本発明の目的は、MRI装置を用いた、被検体の周期的体動情報に同期させた撮像において、所望の画像コントラストを維持して、撮像時間を短縮することにある。
上記の目的を達成するために、本発明は、周期的体動を有する被検体の該周期的体動情報から検出したトリガー情報に同期したエコー信号の同期計測において、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、K空間を複数の部分領域に分割し、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
具体的には、本発明のMRI装置は、被検体の周期的体動情報からトリガー情報を検出する検出部と、トリガー情報に同期して被検体からのエコー信号を計測する同期計測を制御する計測制御部と、エコー信号のデータがK空間に配置されて成るK空間データに基づいて、被検体の画像を取得する演算処理部と、を備え、同期計測は、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、演算処理部は、K空間を複数の部分領域に分割し、計測制御部は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
また、本発明の同期撮像方法は、被検体の周期的体動情報からトリガー情報を検出する検出工程と、トリガー情報に同期して被検体からのエコー信号を計測する同期計測を制御する計測制工程と、エコー信号のデータが配置されるK空間を複数の部分領域に分割する工程と、エコー信号を用いて、被検体の画像を取得する演算処理工程と、を備え、同期計測は、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設け、計測制御工程は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
本発明のMRI装置及び同期撮像方法によれば、被検体の周期的体動情報に同期させた撮像において、所望の画像コントラストを維持して、撮像時間を短縮することが可能になる。
本発明に係るMRI装置の一実施例の全体基本構成を示すブロック図。 本発明に係る3次元K空間の分割の一例を示す図。 本発明に係る3次元K空間の分割設定用GUIの一例を示す図。 本発明に係る3次元K空間の分割設定用GUIの他の一例を示す図。 本発明の第1の実施形態に係る心電同期計測の一例を示す図。(a)は、トリガー信号としての心電波形のR波と、エコー信号計測のタイムテーブルをそれぞれ示す図。(b)は、エコー信号計測のためのパルスシーケンスのシーケンスチャートを示す図。 本発明の第1の実施形態の処理フローを示すフローチャート。 1心拍動周期に計測するエコーデータ数(N)の決定処理フローを示すフローチャート。 本発明の第2の実施形態に係る心電同期計測の一例を示す図。(a)は、T2強調画像を取得する場合のトリガー信号とエコー信号計測のタイムテーブルをそれぞれ示す図(b)は、エコー信号計測のためのパルスシーケンスのシーケンスチャートを示す図。 本発明の第2の実施形態の処理フローを示すフローチャート。 部分領域に応じたトリガー信号待機回数(NT)の設定処理フローを示すフローチャート。
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体の断層画像を得るもので、図1に示すように、MRI装置は静磁場発生部2と、傾斜磁場発生部3と、送信部5と、受信部6と、情報処理部7と、計測制御部4と、を備えて構成される。
静磁場発生部2は、垂直磁場方式であれば、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
傾斜磁場発生部3は、MRI装置の座標系(静止座標系)であるX，Y，Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成る。後述の計測制御部4からの命令に従ってそれぞれのコイルの傾斜磁場電源10が駆動されることにより、X，Y，Zの3軸方向の傾斜磁場Gx，Gy，Gzが、被検体1が横たわる静磁場空間に印加される。撮像時には、スライス面(撮像断面)に直交する方向にスライス方向傾斜磁場パルス(Gs)が印加されて被検体1に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード方向傾斜磁場パルス(Gp)と周波数エンコード方向傾斜磁場パルス(Gf)とが印加されて、エコー信号にそれぞれの方向の位置情報がエンコードされる。
送信部5は、被検体1の生体組織を構成する原子の原子核スピンにNMR現象を誘起するために、被検体1に高周波磁場パルス(以下、「RFパルス」という)を照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとから成る。高周波発振器11から出力された高周波パルスが、計測制御部4からの指令によるタイミングで変調器12により振幅変調され、この振幅変調された高周波パルスが高周波増幅器13で増幅された後に、被検体1に近接して配置された高周波コイル14aに供給されることにより、RFパルスが被検体1に照射される。
受信部6は、被検体1の生体組織を構成する原子核スピンのNMR現象により放出されるエコー信号を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射されたRFパルスによって誘起された被検体1の応答のエコー信号が被検体1に近接して配置された高周波コイル14bで検出され、信号増幅器15で増幅された後、計測制御部4からの指令によるタイミングで、直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、エコー信号として演算処理部7に送られる。以下デジタル量に変換されたエコー信号をエコーデータという。
計測制御部4は、ある所定のパルスシーケンスに基づいて、傾斜磁場発生部3、送信部5、および受信部6を制御して、RFパルスと傾斜磁場パルスの印加と、エコー信号の計測とを、繰り返す制御手段である。計測制御部4は、CPU8の制御で動作し、被検体1の断層画像の再構成に必要なエコー信号収集に必要な種々の命令を傾斜磁場発生部3、送信部5、および受信部6に送って、これらを制御する。
情報処理部7は、各種データ処理と処理結果の表示及び保存等を行うもので、CPU(演算処理部)8と、光ディスク19や磁気ディスク18等の外部記憶装置と、ディスプレイ20と、から成る。受信部6からのエコーデータがCPU8に入力されると、CPU8内のK空間に対応するメモリに、このエコーデータが記憶される(以下、エコー信号またはエコーデータをK空間に配置する旨の記載は、エコーデータがこのメモリに書き込まれて記憶されることを意味する。そして、K空間に配置されたエコーデータをK空間データという)。そして、CPU8はこのK空間データに対して、信号処理、画像再構成等の演算処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置に記録する。
操作部25は、操作者からの、MRI装置の各種制御情報や上記演算処理部7で行う処理の制御情報の入力を受け付け、トラックボール又はマウス23、及び、キーボード24から成る。この操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。
本発明に係るMRI装置は、さらに、被検体に装着されて被検体から心電波形信号を取得する心電極31と、心電極からの信号が入力されて、被検体の心電波形及びそのR波(トリガー信号)を検出する心電波形モニタ32とを備える。心電波形モニタ32が検出した心電波形情報(周期的体動情報の一例)は、計測制御部4を介してCPU8に入力され、計測制御部4がトリガー信号(トリガー情報)に同期して、所定のパルスシーケンスに基づいて上記各部を制御することにより同期計測を制御する。
なお、図1において、送信側の高周波コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生部2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置されている。また、受信側の高周波コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。
現在のMRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
(K空間の分割)
最初に、本発明に係るK空間の分割について説明する。本発明のMRI装置及び同期撮像方法は、被検体の周期的体動情報から検出したトリガー情報に同期した同期計測において、エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設ける。そして、K空間を複数の部分領域に分割し、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、第1の期間と第2の期間の少なくとも一方を異ならせる。
K空間が2次元である場合には、K空間の分割を位相エンコード方向に行い、K空間が3次元である場合には、K空間の分割をスライスエンコード方向と位相エンコード方向の内の少なくとも一方向に行う。K空間の分割の一例を図2に示す。図2は、3次元K空間200を、スライスエンコード方向(kz)に、原点を通るSlice(kz)軸に垂直な平面に関して正負対称に、正負それぞれの側を3つの領域(領域201、202、203)に分割する例を示す。具体的には、CPU8は、3次元K空間を、原点を含む低空間周波数側の部分領域201と、その両側の高空間周波数側をそれぞれ2つの部分領域202-1及び203-1と、202-2及び203-2 (ハイフン"-"以降の数字は、1がK空間スライスエンコード方向(Kz)の正側、2が負側を意味する)とに分割する。さらに、分割の割合は操作者の設定に応じて変更可能とする。図2は、部分領域201が3次元K空間全体の50%と設定され、部分領域202と部分領域203はそれぞれ3次元K空間全体の25%ずつとなるように設定されて3次元K空間が分割された例を示す。なお、さらに位相エンコード方向にもK空間を分割してもよい。
各分割領域の3次元K空間全体に対する割合は、操作者が変更可能な撮像パラメータとする。この各分割領域のK空間全体に対する割合の設定は、例えば図3に示すようなGUI300を介して操作者が行うことを可能とする。図3に示す例は、操作者がGUI300上の3次元K空間の分割数Segment#(301)を入力設定した後、CPU8は各部分領域の全体に対する割合(Data Rate)をそれぞれ決定するためのタグをディスプレイ20に表示させ、操作者が部分領域毎の3次元K空間全体に対する割合タグ(302,303,304)に所望の数値を入力設定する。これらの操作者の設定入力に基づいて、CPU8はそれぞれの部分領域の割合を決定すると共に、この設定情報を計測制御部4に通知する。
或いは、図4に示すように、3次元K空間の分割数Segment#301と、K空間の低空間周波数側の部分領域の、K空間全体に対する割合Data Rate302のみを操作者に入力設定してもらい、CPU8は、残りの部分領域のK空間全体に対する割合を均等に分割するものとして、残りの部分領域のK空間全体に対する割合を算出して決定するとともに、各部分領域の割合を計測制御部4に通知することも可能である。
以下、体動として心臓の拍動を、周期的体動情報として心電波形を想定し、心電波形のR波をトリガー信号(トリガー情報の一例)として、このトリガー信号に同期してエコー信号の計測を制御する心電同期計測を例にして、本発明の各実施形態を説明する。しかし、本発明は、心電同期計測に限らず、血流の脈動や関節の周期的運動等の他の周期的体動にも適用可能である。
(第1の実施形態)
次に、本発明のMRI装置及び同期撮像方法の第1の実施形態を説明する。本実施形態は、第1の期間として、トリガー情報とエコー信号の計測期間との間にディレイ時間を設定し、高空間周波数側の部分領域に対応するエコー信号の計測におけるディレイ時間を、低空間周波側の部分領域に対応するエコー信号の計測におけるディレイ時間よりも短くする。本実施形態は、例えばＴ1強調画像の取得に好適である。以下、図5〜図7に基づいて本実施形態を説明する。
最初に本実施形態の概要を、図5を用いて説明する。図5は、被検体の心電波形を検出しそのR波をトリガー信号としてこれに同期して行う、本実施形態に係る心電同期計測の一例を示す。図5(a)は、トリガー信号としての心電波形のR波と、エコー信号計測のタイムテーブルをそれぞれ示している。また、図5(b)は、エコー信号計測のためのパルスシーケンスの一例のシーケンスチャートを示している。
操作者により図5に示す心電同期計測が起動されると、計測制御部4は、トリガー信号501に同期してパルスシーケンスを駆動する。計測制御部4は、このパルスシーケンスを短い繰り返し時間(TR)で繰り返してエコー信号を計測(503)するが、部分領域201のエコー信号計測中では、トリガー信号から一定時間(Delay time)経過まではエコー信号収集は行わず空打ち502を行う。この空打ち502は、同じパルスシーケンスであるがエコー信号を計測しないか、計測しても画像再構成に用いない。このようにトリガー信号から一定のDelay timeを設けて空打ちを行う理由は、Delay timeの時間帯は心臓の収縮期で心臓の動きが速いためであり、心臓の動きが比較的ゆっくりしている拡張期にエコー信号の計測を行うためである。これにより、心臓の動きに起因して画像上に発生するアーチファクトを低減して高画質の画像を取得することが可能になる。
計測制御部4は、トリガー信号から一定時間(Delay time)経過後、被検体の体動が少ない時間帯に、短い繰り返し時間(TR)で位相エンコードを変えてパルスシーケンスを少なくとも1回繰り返して、部分領域201に対応するエコー信号を少なくとも1以上計測する。そして、計測制御部4は、このようなDelay time後に短い繰り返し時間(TR)でパルスシーケンスを繰り返して行うエコー信号の計測を、1心拍動周期毎に複数心拍繰り返し実行して、部分領域201に対応するエコー信号を全て計測する。
次に、部分領域202のエコー信号の計測へ移行すると、計測制御部4はDelay Timeを短く設定し、1心拍動周期に計測するエコー信号の数を増加させる。図5では、部分領域201のエコー信号計測時と比較して、部分領域202のエコー信号計測時のDelay Timeが、パルシーケンスの繰り返し時間(TR)の一繰り返し時間分(1TR)短くなっているが、Delay Timeの短縮は、1TRに限定されるものではなく、また、パルシーケンスの繰り返し時間(TR)単位に限定されるものでもなく、任意に設定することが可能である。
同様にして、部分領域203のエコー信号計測へ移行した段階で、計測制御部4は部分領域202のエコー信号計測時と比較して、部分領域203のエコー信号計測時のDelay Timeをさらに短く設定し、1心拍動周期に計測するエコー信号の数を増加させる。図5の例は、Delay Time=0(ゼロ)と設定した場合を示す。すなわち、計測制御部4はトリガー信号直後から空打ち無しにエコー信号の計測を行う。
一般的に、低空間周波数領域のエコー信号は、その強度が大きく、画像のコントラストを含む画質を支配する重要な信号である。このような低空間周波数領域のエコー信号を、心臓の動きが比較的ゆっくりしている拡張期に計測できるようDelay Timeを設定することにより、動きに起因するアーチファクトを低減して、所望のコントラストを有する画像を取得できるようになる。
一方、高空間周波数領域のエコー信号は、低空間周波数領域のエコー信号と異なり、画像のコントラストに影響するものではく、画像の分解能に寄与する。また、高空間周波数領域のエコー信号は強度も小さいので画像への影響が少ない。そのため、本実施形態のように、部分領域202と203のエコー信号を、Delay Timeを少なくして収縮期にかかる時間帯で計測して、これらのエコー信号を画像再構成に用いても画像への影響が軽微となる。そのため、本実施形態の様に、高空間周波数領域のエコー信号ほど、Delay Timeを少なくして計測して良いことになる。逆に、高空間周波数領域のエコー信号計測の際のDelay Timeを短くすることによって、計測できるエコー信号数を増加できるので、全体の撮像時間を短縮できる。
また、本実施形態は、前述の通り、1心拍動周期内で短い繰り返し時間(TR)でパルスシーケンスを繰り返してエコー信号を計測する。そのため、本実施形態はT1強調画像の取得に好適となる。そこで、T1強調画像の取得に好適なパルスシーケンスとして、例えば図5(b)に示すスピンエコー(SE)シーケンスを用いる。図5(b)に示すシーケンスチャートのRF、Gs、Gp、Gf及びEchoは、それぞれRFパルス、スライス傾斜磁場、位相エンコード傾斜磁場、周波数エンコード傾斜磁場、及びエコー信号を表す。スライス選択傾斜磁場512を被検体に印加した状態で、90度RFパルス511を被検体に印加して所望の領域の磁化を90度励起して横磁化を発生させる。その直後に、スライスリフェーズ傾斜磁場513を印加して、所望領域の励起に伴う横磁化の位相分散を再収束させる。次に、スライスエンコード傾斜磁場514を印加して、エコー信号の位相にスライスエンコード方向の空間情報をエンコードする。また、対称なエコー信号を発生させるために周波数ディフェーズ傾斜磁場515を印加して、横磁化の位相を分散させる。その後に、スライス傾斜磁場517と180度RFパルス516を印加して、横磁化を180度反転させて再収束させる。そして、位相エンコード傾斜磁場518を印加して、エコー信号の位相に位相エンコード方向の空間情報をエンコードし、周波数エンコード傾斜磁場519を印加しながらエコー信号520を計測する。これによりエコー信号520に周波数エンコード方向の空間情報がエンコードされる。以上の各パルスを、スライスエンコード傾斜磁場514と位相エンコード傾斜磁場518の少なくとも一方を変えながら短い繰り返し時間(TR)で繰り返して、各部分領域のエコー信号を計測する。
また、本実施形態は、前述のDelay Timeを利用して空打ちを行うことにより、エコー信号計測前に励起領域からのエコー信号の強度を安定化させる。特に低空間周波数領域に対応するエコー信号の計測の際には、積極的にDelay Timeを利用して空打ちを行うことにより、強度の安定したエコー信号を計測して高画質の画像を取得する。
以上迄が本実施形態の概要である。次に本実施形態の動作フローを、図6,図7に示すフローチャートを用いて説明する。図6は、本実施形態の動作フローの概要を、図7は特にステップ602の処理の詳細を示す。これらの処理フローは、プログラムとして外部記憶装置に記憶されており、必要に応じてCPU8がそのメモリに読み込んで実行することにより、実行される。なお、図6,図7は、例えば図3，図4に示したGUIを介して事前に操作者が3次元K空間の分割数を入力設定し、CPU8がその入力された設定値に応じて、K空間の分割数を決定した後の処理フローを表している。
ステップ601で、操作者により心電同期計測が起動されると、計測制御部4は、被検体の心電波形におけるR波を、同期計測を行うためのトリガー信号として、その読み込みを開始する。
ステップ602で、3次元K空間の部分領域に応じた、1心拍動周期で計測するエコー信号数(N)やDelay time(Td)等が決定される。この決定処理の詳細は図7に基づいて後述する。なお、ステップ601と602の実行順序はどちらが先でも良い。
ステップ603で、計測制御部4は、トリガー信号を待つ。
ステップ604で、トリガー信号が受信されると、計測制御部4は、ステップ602で決定されたDelay time(Td)だけ待機する処理を行う。このDelay time(Td)期間中には、計測制御部4は、前述の通り空打ち605を行う。
ステップ606で、Delay time(Td)が経過した時点が最も被検体の動きが少ないタイミングであるとして、図5に示すような、1心拍動周期のエコー信号計測が実行される。つまり、計測制御部4は、スライスエンコード傾斜磁場もしくは位相エンコード傾斜磁場の印加量を変えて、ステップ602で設定された、1心拍動周期に計測するエコー信号数分の計測を行う。
ステップ607で、1心拍動周期のエコー信号計測が終了し、且つ心電同期計測が終了でない場合はステップ601に戻って、再度トリガー信号の読み込みの開始と、K空間の部分領域に応じた1心拍動周期に収集するエコー信号数の決定を行い、撮像終了となるまでステップ601〜606を繰り返し実行する。
以上迄が、本実施形態の処理フローの概要である。次にステップ602の1心拍動周期に計測するエコー信号数の決定処理フローを、図7に示すフローチャートに基づいて説明する。
ステップ701で、CPU8は、部分領域201のエコー信号計測時におけるDelay Time(Td1)と計測するエコー信号数(N1)を基準として、部分領域202のエコー信号計測時における、Delay Time(Td2)、計測エコー信号数(N2)と、部分領域203のエコー信号計測時における、Delay Time(Td3)、計測エコー信号数(N3)と、を求める。なお、計測エコー信号数（N1)とDelay Time(Td1)は、例えば、操作者が計測エコー信号数（N1)とDelay time(Td1)を設定する、操作者が計測エコー信号数（N1)を設定した後、CPU8がDelay Time(Td1)を算出する、操作者がDelay time(Td1)を設定した後、CPU8が計測エコー信号数（N1)を算出する、のいずれでも良い。分割領域202、分割領域203でのDelay Time、計測エコー信号数は、例えば、撮像パラメータとして操作者により設定される撮像高速化率Rapid Rateを用いて下式で求めることが出来る。
N2=N1×(1＋Rapid Rate/100), Td2=Td1＋TR×(N2-N1) (1)
N3=N2×(1＋Rapid Rate/100), Td3=Td2＋TR×(N3-N2)
撮像高速化率Rapid Rateは、例えば、最終的にDelay Timeが最小値となるように、つまり、最後の部分領域(この場合では部分領域203)のエコー信号計測時のDelay Timeが最小値となるように、自動的に算出することも可能である。
ステップ702で、CPU8は、部分領域202のエコー信号の計測であるか否かを判断し、そうであれば、ステップ703で、Delay Timeと計測エコー信号数とを、(1)式で求めた部分領域202の値Td2,N2にそれぞれ設定し、計測制御部4に通知する。
ステップ704で、CPU8は、部分領域203のエコー信号の計測であるか否かを判断し、そうであれば、ステップ705で、Delay Timeと計測エコー信号数とを、(1)式で求めた部分領域203の値Td3,N3にそれぞれ設定し、計測制御部4に通知する。
以上迄が、ステップ602の1心拍動周期に計測するエコー信号数の決定処理フローの説明である。
以上説明したように、本実施形態のMRI装置及び同期撮像方法によれば、3次元K空間をスライスエンコード方向に複数分割し、分割された部分領域のスライスエンコード方向の位置に応じて、エコー信号を計測する際のDelay Timeを異ならせることにより、1心拍動周期に計測するエコー信号数を増加させて撮像時間の短縮化を行うことが可能となる。具体的には、低空間周波数側の部分領域のエコー信号の計測の際のDelay Timeを所望の画像コントラストとなるように設定することにより3次元画像のコントラストを所望の状態に維持できる。そして、高空間周波数側の部分領域のエコー信号の計測の際のDelay Timeを低空間周波数側の部分領域のエコー信号の計測の際のDelay Timeよりも短くすることにより、高空間周波数側の部分領域において1心拍動周期に計測するエコー信号数を増加させて撮像時間の短縮化を行うことができる。
なお、以上の本実施形態の説明では、3次元計測を例にしたが、2次元の心電同期計測においても、2次元K空間を位相エンコード方向に複数分割し、分割された部分領域の位相エンコード方向の位置に応じて、エコー信号を計測する際のDelay Timeを異ならせることにより、2次元画像のコントラストを所望の状態に維持しつつ、1心拍動周期に計測するエコー信号の数を増加させて撮像時間の短縮化を行うことが可能である。
(第2の実施形態)
次に、本発明のMRI装置及び同期撮像方法の第2の実施形態を説明する。本実施形態は、第2の期間として、一つのエコー信号計測期間から、次のエコー信号計測期間のトリガーとなるトリガー情報までの待ち時間を設定し、高空間周波数側の部分領域に対応するエコー信号の計測における待ち時間を、低空間周波数側の部分領域に対応するエコー信号の計測における待ち時間よりも短くする。好ましくは、周期的体動の一周期の整数倍だけ、待ち時間を異ならせる。すなわち、エコー信号の計測を行うためのパルスシーケンスの繰り返し時間について、高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くする。本実施形態は、例えばＴ2強調画像の取得に好適である。以下、前述の第1の実施形態と異なる点のみを図8〜図10に基づいて説明し、同じ点の説明は省略する。
最初に本実施形態の概要を、図8を用いて説明する。図8は、被検体の心電波形を検出しそのR波をトリガー信号としてこれに同期して行う、本実施形態に係る心電同期計測の一例を示す。図8(a)は、T2強調画像を取得する場合のトリガー信号とエコー信号計測のタイムテーブルをそれぞれ示している。図8(b)は、エコー信号計測のためのパルスシーケンスの一例のシーケンスチャートを示している。
計測制御部4は、トリガー信号801-1aの検出から一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号計測802-1aを行う。このパルスシーケンスには、図8(b)に示すように、複数のエコー信号を計測可能なFSEシーケンスやEPIシーケンスなどのマルチエコータイプのパルスシーケンスが用いられる。詳細は後述する。計測制御部4は、このようなパルスシーケンスに基づいて所定数のエコー信号の計測を制御する。この心拍動周期でのエコー信号計測後、T2強調画像を得るために、数(例えば2,3)心拍を待ち時間として、この待ち時間経過後、再度トリガー信号801-1bの検出後にパルスシーケンスを起動してエコー信号の計測802-1bを行う。
計測制御部4は、このような、Delay Time後のエコー信号の計測と、数心拍動周期分の待ち時間の待機と、を一単位とする処理を、複数心拍動周期毎(ここでは3心拍動周期毎)に複数回繰り返して、つまりパルスシーケンスの繰り返し時間(TR1)=3心拍動周期として、部分領域201のエコー信号を計測する。
次に、分割領域202のエコー信号計測へ移行すると、計測制御部4は、トリガー信号801-2aの検出から、部分領域201の計測時と同じ一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号の計測802-2aを行う。その後の待ち時間は、分割領域201のエコー信号計測時と比較して短く設定される。図8(a)の例では、パルスシーケンス終了後の待ち時間が、分割領域201の場合と比較して1心拍動周期分短くされている。ただし、本実施形態は、待ち時間の短縮度合いとして1心拍動周期に規定されるものではなく、任意に設定可能である。そして、この待ち時間後に、再度トリガー信号801-2bの検出後にパルスシーケンスを起動してエコー信号の計測802-2bを行う。
計測制御部4は、このような、Delay Time後のエコー信号の計測と、部分領域201のエコー信号計測時よりも短い1以上の心拍動周期分の待ち時間の待機と、を一単位とする処理を、複数心拍動周期毎(ここでは2心拍動周期毎)に複数回繰り返して、つまりパルスシーケンスの繰り返し時間(TR2)=2心拍動周期として、部分領域202のエコー信号を計測する。
同様にして、分割領域203のエコー信号計測へ移行した段階で、計測制御部4は、パルスシーケンス動作後の待ち時間を分割領域202のエコー信号計測時よりも短く設定して、エコー信号の計測周期を短縮させる。つまり、分割領域203のエコー信号計測へ移行すると、計測制御部4は、トリガー信号801-3aの検出から、部分領域201の計測時と同じ一定時間(Delay time)経過後、パルスシーケンスを起動してエコー信号の計測802-3aを行う。その後の待ち時間は、分割領域202のエコー信号計測時と比較して短く設定される。図8(a)の例では、パルスシーケンス終了後の待ち時間が、分割領域202の場合と比較して1心拍動周期分短くされて、結果として1心拍動周期以下の時間とされている。そして、この待ち時間後に、再度トリガー信号801-3bの検出後にパルスシーケンスを起動してエコー信号の計測802-3bを行う。つまり、計測制御部4は、パルスシーケンスの繰り返し時間(TR3)=1心拍動周期として、部分領域203のエコー信号を計測する。
以上の説明を纏めると、本実施形態は、エコー信号の計測を行うためのパルスシーケンスの繰り返し時間(TR)について、待ち時間を異ならせることにより、高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くする。
なお、高空間周波数側の部分領域に対応するエコー信号の計測においては、Delay time＝0(ゼロ)としてもよい。
また、本実施形態は、前述の通り、部分領域に応じて1以上の心拍動周期毎にマルチエコータイプのパルスシーケンスを繰り返してエコー信号を計測する。特に低空間周波数側の部分領域では、複数心拍動周期毎にパルスシーケンスを繰り返す。そのため、本実施形態はT2強調画像の取得に好適となる。そこで、T2強調画像の取得に好適なパルスシーケンスとして、例えば図8(b)に示すFSEシーケンスを用いる。図8(b)に示すシーケンスチャートのRF、Gs、Gp、Gf及びEchoは、図5(b)と同様である。FSEシーケンスは、図5(b)に示すスピンエコーシーケンスの180度RFパルス(516-1〜516-6)とスライス傾斜磁場(517-1〜517-5)と周波数エンコード傾斜磁場(519-1〜519-5)とを繰り返して複数のエコー信号(520-1〜520-5)を計測するパルスシーケンスであって、各エコー信号の計測の前後で、絶対量が同じで極性の異なる印加量となる一対の位相エンコード傾斜磁場(801-1a、801-1b〜801-5a、ただし801-5b、801-3は印加量ゼロ)を、エコー信号の計測毎にその印加量を変えて印加する。そのため位相エンコード傾斜磁場518は印加しない。これにより、位相エンコードの異なる複数のエコー信号(520-1〜520-5)を一回の90度RFパルス511で計測することが可能になる。
以上迄が本実施形態の概要である。次に本実施形態の動作フローを、図9, 図10に示すフローチャートを用いて説明する。図9は、本実施形態の動作フローの概要を、図10は特にステップ901の処理の詳細を示す。これらの処理フローは、プログラムとして外部記憶装置に記憶されており、必要に応じてCPU8がそのメモリに読み込んで実行することにより、実行される。なお、前述の第1の実施形態の図6,図7の説明と同様に、操作者が3次元K空間の分割数を入力設定し、CPU8がその入力された設定値に応じて、K空間の分割数を決定したものとする。
ステップ901で、操作者により心電同期計測が起動されると、撮像を行うK空間の部分領域に応じてトリガー信号待機回数(NT)の設定が行われる。この設定処理の詳細は図10に基づいて後述する。
ステップ902で、上記トリガー信号待機回数(NT)の設定完了後、計測制御部4は、トリガー信号の読み込みを開始する。
ステップ903で、計測制御部4は、トリガー信号を受信すると、ステップ904で、設定されたたトリガー信号待機回数(NT)から1を減算する。
ステップ905で、計測制御部4は、事前に設定されたDelay time(Td)だけ待機する。このDelay Ｔｉｍｅ期間中に前述の空打ちを行えば、その後に計測するエコー信号の強度が安定して高画質の画像を得ることができる。
ステップ906で、Delay time(Td)が経過した時点が最も被検体の動きが少ないタイミングであるとして、計測制御部4は、例えば図8(b)に示すようなパルスシーケンスを実行する。
ステップ907で、このイメージングスキャンが終了した後、計測制御部4は、再度トリガー信号の読み込みを行い、トリガー信号を受信すると、設定したトリガー信号待機回数(NT)から1を減算する処理ステップ908を行い、NT=0となるまで実行(909)する。
ステップ910で、上記待機時間が終了すると、撮像終了でない場合、ステップ901に戻って、再度トリガー信号待機回数(NT)の設定が行われ、撮像終了となるまで上記ステップ901〜909が繰り返される。
以上迄が、本実施形態の処理フローの概要である。次に、ステップ901のK空間の部分領域に応じたトリガー信号待機回数(NT)の設定処理フローの概要を図10に示すフローチャートに基づいて説明する。
ステップ1001で、CPU8は、分割領域201計測時におけるトリガー信号待機回数(NT1)を基準として、部分領域202のトリガー信号待機回数(NT2)と、部分領域203のトリガー信号待機回数(NT3)の設定を行う。なお、トリガー信号待機回数(NT1)は、例えば、操作者がトリガー信号待機回数(NT1)を設定する、予め所定の値を記憶しておく、のいずれでも良い。分割領域202、分割領域203でのトリガー信号待機回数は、例えば、撮像パラメータとして操作者により設定される撮像高速化率Rapid Rateを用いて下式で求めることが出来る。
NT2=NT1×(1−Rapid Rate/100) (2)
NT3=NT2×(1−Rapid Rate/100)
撮像高速化率Rapid Rateは、例えば、前述の第1の実施形態と同様に、最終的にDelay Timeが最小値となるように自動的に算出することも可能である。
ステップ1002で、CPU8は、部分領域202のエコー信号の計測であるか否かを判断し、そうであれば、ステップ1003で、トリガー信号待機回数を、(2)式で求めた部分領域202の値NT2に設定する。
ステップ1004で、CPU8は、部分領域203のエコー信号の計測であるか否かを判断し、そうであれば、ステップ1005で、トリガー信号待機回数を、(2)式で求めた部分領域203の値NT2に設定する。
以上迄が、ステップ901のK空間の部分領域に応じたトリガー信号待機回数(NT)の設定処理フローの説明である。
以上説明したように、本実施形態のMRI装置及び同期撮像方法によれば、3次元K空間をスライスエンコード方向に複数分割し、分割された部分領域のスライスエンコード方向の位置に応じて、エコー信号を計測した後の待ち時間を異ならせることにより、撮像時間の短縮化を行うことが可能となる。具体的には、低空間周波数側の部分領域のエコー信号の計測の際の待ち時間を所望の画像コントラストとなるように設定することにより3次元画像のコントラストを所望の状態に維持できる。そして、高空間周波数側の部分領域のエコー信号の計測の際の待ち時間を低空間周波数側の部分領域のエコー信号の計測の際の待ち時間よりも短くすることにより、撮像時間の短縮化を行うことが可能となる。
以上までが、本発明のMRI装置及び同期撮像方法に係る心電同期計測の各実施形態の説明である。しかし、本発明のMRI装置及び同期撮像方法は、上記各実施形態の説明で開示された内容にとどまらず、本発明の趣旨を踏まえた上で他の形態を取り得る。例えば、第1の実施形態の部分領域毎のDelay Timeの制御と第2の実施形態の部分領域毎の待ち時間の制御を組み合わせて実施してもよい。また、第1の実施形態でFSEシーケンスを用いてT2強調画像を取得してもよいし、第2の実施形態で、プロトン強調画像の取得や、スピンエコー(SE)シーケンスを用いてもよい。
1 被検体、2 静磁場発生系、3 傾斜磁場発生系、4 計測制御部、5 送信系、6 受信系、7 信号処理系、8 中央処理装置(CPU)、9 傾斜磁場コイル、10 傾斜磁場電源、11 高周波発振器、12 変調器、13 高周波増幅器、14a 高周波コイル(送信コイル)、14b 高周波コイル(受信コイル)、15 信号増幅器、16 直交位相検波器、17 A/D変換器、18 磁気ディスク、19 光ディスク、20 ディスプレイ、21 ROM、22 RAM、23 トラックボール又はマウス、24 キーボード、51 ガントリ、52 テーブル、53 筐体、54 処理装置

Claims (18)

  1. 被検体の周期的体動情報からトリガー情報を検出する検出部と、
    前記トリガー情報に同期して前記被検体からのエコー信号を計測する同期計測を制御する計測制御部と、
    前記エコー信号のデータがK空間に配置されて成るK空間データに基づいて、前記被検体の画像を取得する演算処理部と、
    を備え、
    前記同期計測は、前記エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設けた磁気共鳴イメージング装置であって、
    前記演算処理部は、前記K空間を複数の部分領域に分割し、
    前記計測制御部は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、前記第1の期間と前記第2の期間の少なくとも一方を異ならせることを特徴とする磁気共鳴イメージング装置。
  2. 請求項1記載の磁気共鳴イメージング装置において、
    前記K空間はスライスエンコード方向と位相エンコード方向を含む3次元空間であり、前記K空間の分割はスライスエンコード方向に行われることを特徴とする磁気共鳴イメージング装置。
  3. 請求項1記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記第1の期間として、前記トリガー情報と前記エコー信号の計測期間との間にディレイ時間を設けることを特徴とする磁気共鳴イメージング装置。
  4. 請求項3記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記高空間周波数側の部分領域に対応するエコー信号の計測における前記ディレイ時間を、前記低空間周波側の部分領域に対応するエコー信号の計測における前記ディレイ時間よりも短くすることを特徴とする磁気共鳴イメージング装置。
  5. 請求項3記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記部分領域に対応するエコー信号の計測を、1心拍動周期毎に繰り返すことを特徴とする磁気共鳴イメージング装置。
  6. 請求項1記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記第2の期間として、一つのエコー信号計測期間から、次のエコー信号計測期間のトリガーとなるトリガー情報までの待ち時間を設けることを特徴とする磁気共鳴イメージング装置。
  7. 請求項6記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記高空間周波数側の部分領域に対応するエコー信号の計測における前記待ち時間を、前記低空間周波数側の部分領域に対応するエコー信号の計測における前記待ち時間よりも短くすることを特徴とする磁気共鳴イメージング装置。
  8. 請求項6記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記周期的体動の一周期の整数倍だけ、前記待ち時間を異ならせることを特徴とする磁気共鳴イメージング装置。
  9. 請求項6記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記低空間周波数側の部分領域に対応するエコー信号の計測を、複数心拍動周期毎に繰り返すことを特徴とする磁気共鳴イメージング装置。
  10. 請求項6記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記第1の期間として、前記トリガー情報と前記エコー信号の計測期間との間にディレイ時間を設け、前記部分領域に依らずに、前記ディレイ時間を同じにすることを特徴とする磁気共鳴イメージング装置。
  11. 請求項1記載の磁気共鳴イメージング装置において、
    前記計測制御部は、前記エコー信号の計測を行うためのパルスシーケンスの繰り返し時間について、前記高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、前記低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くすることを特徴とする磁気共鳴イメージング装置。
  12. 請求項3記載の磁気共鳴イメージング装置において、
    前記計測制御部は、スピンエコーシーケンスを用いて前記エコー信号の計測を行うことを特徴とする磁気共鳴イメージング装置。
  13. 請求項6記載の磁気共鳴イメージング装置において、
    前記計測制御部は、マルチエコータイプのパルスシーケンスを用いて前記エコー信号の計測を行うことを特徴とする磁気共鳴イメージング装置。
  14. 請求項1記載の磁気共鳴イメージング装置において、
    前記K空間を複数の部分領域に分割するための分割数と、少なくとも一つの部分領域の前記K空間に対する割合の設定を受け付ける入力設定部を備えたことを特徴とする磁気共鳴イメージング装置。
  15. 被検体の周期的体動情報からトリガー情報を検出する検出工程と、
    前記トリガー情報に同期して前記被検体からのエコー信号を計測する同期計測を制御する計測制工程と、
    前記エコー信号を用いて、前記被検体の画像を取得する演算処理工程と、
    を備え、
    前記同期計測は、前記エコー信号の計測期間の、前に第1の期間と後に第2の期間の少なくとも一方を設けた同期撮像方法であって、
    前記エコー信号のデータが配置されるK空間を複数の部分領域に分割する工程を更に備え、
    前記計測制御工程は、低空間周波数側の部分領域に対応するエコー信号の計測と、高空間周波数側の部分領域に対応するエコー信号の計測とで、前記第1の期間と前記第2の期間の少なくとも一方を異ならせることを特徴とする同期撮像方法。
  16. 請求項15記載の同期撮像方法において、
    前記計測制御工程は、前記第1の期間として、前記トリガー情報と前記エコー信号の計測期間との間にディレイ時間を設けることを特徴とする同期撮像方法。
  17. 請求項15記載の同期撮像方法において、
    前記計測制御工程は、前記第2の期間として、一つのエコー信号計測期間から、次のエコー信号計測期間のトリガーとなるトリガー情報までの待ち時間を設けることを特徴とする同期撮像方法。
  18. 請求項15記載の同期撮像方法において、
    前記計測制御工程は、前記エコー信号の計測を行うためのパルスシーケンスの繰り返し時間について、前記高空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間を、前記低空間周波数側の部分領域に対応するエコー信号の計測の場合の繰り返し時間よりも短くすることを特徴とする同期撮像方法。
JP2010541301A 2008-12-04 2009-11-26 磁気共鳴イメージング装置および同期撮像方法 Expired - Fee Related JP5372015B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2010541301A JP5372015B2 (ja) 2008-12-04 2009-11-26 磁気共鳴イメージング装置および同期撮像方法

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2008309192 2008-12-04
JP2008309192 2008-12-04
PCT/JP2009/069945 WO2010064572A1 (ja) 2008-12-04 2009-11-26 磁気共鳴イメージング装置および同期撮像方法
JP2010541301A JP5372015B2 (ja) 2008-12-04 2009-11-26 磁気共鳴イメージング装置および同期撮像方法

Publications (2)

Publication Number Publication Date
JPWO2010064572A1 true JPWO2010064572A1 (ja) 2012-05-10
JP5372015B2 JP5372015B2 (ja) 2013-12-18

Family

ID=42233223

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2010541301A Expired - Fee Related JP5372015B2 (ja) 2008-12-04 2009-11-26 磁気共鳴イメージング装置および同期撮像方法

Country Status (4)

Country Link
US (1) US20110237931A1 (ja)
JP (1) JP5372015B2 (ja)
CN (1) CN102238909B (ja)
WO (1) WO2010064572A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102238909B (zh) * 2008-12-04 2013-09-04 株式会社日立医疗器械 磁共振成像装置以及同步摄像方法
CN103185876B (zh) * 2011-12-30 2015-05-13 西门子(深圳)磁共振有限公司 磁共振成像方法及磁共振成像装置
CN104427934B (zh) * 2012-08-13 2017-03-01 株式会社日立制作所 磁共振成像装置以及磁共振成像方法
DE102012216773B4 (de) * 2012-09-19 2018-11-08 Siemens Healthcare Gmbh Aufnahme von mindestens zwei Messdatensätzen eines Untersuchungsobjekts mittels einer Magnetresonanzanlage
CN108567422B (zh) * 2018-03-02 2021-03-05 上海东软医疗科技有限公司 一种血管成像方法和装置
CN113644920B (zh) * 2021-06-16 2022-10-14 北京协同创新研究院 一种用于脑部成像的微波多路收发系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19901171C2 (de) * 1999-01-14 2001-12-13 Axel Haase Verfahren und Vorrichtung zum Gewinnen von Daten für Magnetresonanz-Bildgebung
JP4632535B2 (ja) * 2000-12-27 2011-02-16 株式会社東芝 Ｍｒｉ装置
US6801800B2 (en) * 1999-11-29 2004-10-05 Kabushiki Kaisha Toshiba MR imaging using ECG-prep scan
US6611701B2 (en) * 2000-12-30 2003-08-26 Ge Medical Systems Global Technology Company, Llc Method and apparatus for fast breath-held 3D MR data acquisition using variable sampling
EP1764620A1 (en) * 2005-09-19 2007-03-21 Kabushiki Kaisha Toshiba MRI involving retrospective data extraction
US7840045B2 (en) * 2006-04-21 2010-11-23 The University Of Utah Research Foundation Method and system for parallel reconstruction in the K-space domain for application in imaging systems
US9020575B2 (en) * 2006-11-10 2015-04-28 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5269381B2 (ja) * 2006-11-10 2013-08-21 株式会社東芝 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP5019576B2 (ja) * 2006-11-24 2012-09-05 株式会社東芝 磁気共鳴イメージング装置
CN102238909B (zh) * 2008-12-04 2013-09-04 株式会社日立医疗器械 磁共振成像装置以及同步摄像方法

Also Published As

Publication number Publication date
WO2010064572A1 (ja) 2010-06-10
CN102238909A (zh) 2011-11-09
US20110237931A1 (en) 2011-09-29
CN102238909B (zh) 2013-09-04
JP5372015B2 (ja) 2013-12-18

Similar Documents

Publication Publication Date Title
JP5002222B2 (ja) 磁気共鳴イメージング装置
US9301704B2 (en) Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same
JP6467341B2 (ja) 磁気共鳴イメージング装置、画像処理装置、画像診断装置、画像解析装置、ｍｒｉ画像作成方法およびプログラム
EP2335572A1 (en) Magnetic resonance imaging apparatus, and breath-holding imaging method
JP5372015B2 (ja) 磁気共鳴イメージング装置および同期撮像方法
US20120281901A1 (en) Magnetic resonance imaging apparatus and blood vessel image capturing method
JP5536665B2 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2007190114A (ja) 磁気共鳴イメージング装置
US8594767B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2018033691A (ja) 磁気共鳴測定装置および画像処理方法
JP4230875B2 (ja) 磁気共鳴イメージング装置
JP5336731B2 (ja) 磁気共鳴イメージング装置
KR101541290B1 (ko) 자기 공명 신호 측정 방법 및 장치
JP4018107B2 (ja) 磁気共鳴イメージング装置
JP5371620B2 (ja) 核磁気共鳴イメージング装置
JP5319835B2 (ja) 磁気共鳴イメージング装置
JP4994786B2 (ja) 磁気共鳴イメージング装置
JP4454268B2 (ja) 磁気共鳴イメージング装置
JP5421600B2 (ja) 核磁気共鳴イメージング装置および核磁気共鳴イメージング装置の作動方法
JP2015002834A (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2016096989A (ja) 磁気共鳴イメージング装置
JP6157976B2 (ja) 磁気共鳴イメージング装置、及び方法
JP2012254189A (ja) 磁気共鳴イメージング装置
JP2014171833A (ja) 磁気共鳴イメージング装置および非造影ｍｒａ方法
JP5514293B2 (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121115

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20121115

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130716

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130809

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130917

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees