(57)【要約】      （修正有）【課題】積算処理を行うパラレルイメージング法において展開誤差を低減すること。【解決手段】磁気共鳴イメージング装置は、記憶部、再構成部、積算部及び展開部を有する。記憶部は、複数の受信チャネルを用いたパラレルイメージングにより収集された、積算回数に応じた複数のデータセットを複数の受信チャネル毎に記憶する。再構成部は、複数のデータセットに対して複数の受信チャネル毎に再構成処理を施して、複数の折り返し画像を生成する。積算部は、複数の折り返し画像を複数の受信チャネル毎に積算して、複数の受信チャネルと同数の複数の積算画像を生成する。展開部は、複数の受信チャネルに関する複数の感度分布画像を利用して複数の積算画像に展開処理を施して、最終画像を生成する。【選択図】図４
2019005136.tif 000002
【特許請求の範囲】
【請求項１】
複数の受信チャネルを用いたパラレルイメージングにより収集された、積算回数に応じた複数のデータセットを前記複数の受信チャネル毎に記憶する記憶部と、
前記複数のデータセットに対して前記複数の受信チャネル毎に再構成処理を施して、複数の折り返し画像を生成する再構成部と、
前記複数の折り返し画像を前記複数の受信チャネル毎に積算して、前記複数の受信チャネルと同数の複数の積算画像を生成する積算部と、
前記複数の受信チャネルに関する複数の感度分布画像を利用して前記複数の積算画像に展開処理を施して、最終画像を生成する展開部と、
を具備する磁気共鳴イメージング装置。
【請求項２】
前記複数の受信チャネルを用い、前記積算回数と同数の複数のＤＷＩ収集を実行する撮像制御部を更に備える、請求項１記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
パラレルイメージング法において複数のショットで得られた複数の画像に積算処理を行う場合、再構成処理及び展開処理を実行してショット毎の展開画像を生成した後、これら展開画像に積算処理を行い最終画像を生成している。この方法では、積算回数を増やしても、展開処理の入力画像のＳＮＲ（Signal to Noise Ratio）が上がらないので、展開処理の誤差を低減することができない。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】米国特許出願公開第２０１２／３１９６８６号明細書
【特許文献２】米国特許第８０２６７２０号明細書
【特許文献３】米国特許出願公開第２０１５／０４２３４１号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
発明が解決しようとする課題は、積算処理を行うパラレルイメージング法において展開誤差を低減することである。
【課題を解決するための手段】
【０００５】
本実施形態に係る磁気共鳴イメージング装置は、記憶部、再構成部、積算部及び展開部を有する。前記記憶部は、複数の受信チャネルを用いたパラレルイメージングにより収集された、積算回数に応じた複数のデータセットを前記複数の受信チャネル毎に記憶する。前記再構成部は、前記複数のデータセットに対して前記複数の受信チャネル毎に再構成処理を施して、複数の折り返し画像を生成する。前記積算部は、前記複数の折り返し画像を前記複数の受信チャネル毎に積算して、前記複数の受信チャネルと同数の複数の積算画像を生成する。前記展開部は、前記複数の受信チャネルに関する複数の感度分布画像を利用して前記複数の積算画像に展開処理を施して、最終画像を生成する。
【図面の簡単な説明】
【０００６】
【図１】図１は、本実施形態に係る磁気共鳴イメージング装置１の構成を示す図である。
【図２】図２は、本実施形態に係る簡易的なパルスシーケンスを示す図である。
【図３】図３は、図２の１ショット分のＤＷＩ収集のパルスシーケンスの一例を示す図である。
【図４】図４は、本実施形態に係る、処理回路により行われる最終画像の生成処理の過程を模式的に示す図である。
【図５】図５は、比較例による最終画像の生成処理の過程を模式的に示す図である。
【発明を実施するための形態】
【０００７】
以下、図面を参照しながら本実施形態に係わる磁気共鳴イメージング装置を説明する。
【０００８】
図１は、本実施形態に係る磁気共鳴イメージング装置１の構成を示す図である。図１に示すように、磁気共鳴イメージング装置１は、架台１１、寝台１３、傾斜磁場電源２１、送信回路２３、受信回路２５、寝台駆動装置２７、シーケンス制御回路２９及びホストＰＣ５０を有する。
【０００９】
架台１１は、静磁場磁石４１と傾斜磁場コイル４３とを有する。静磁場磁石４１と傾斜磁場コイル４３とは架台１１の筐体に収容されている。架台１１の筐体には中空形状を有するボアが形成されている。架台１１のボア内には送信コイル４５と受信コイル４７とが配置される。
【００１０】
静磁場磁石４１は、中空の略円筒形状を有し、略円筒内部に静磁場を発生する。静磁場磁石４１としては、例えば、永久磁石、超伝導磁石または常伝導磁石等が使用される。ここで、静磁場磁石４１の中心軸をＺ軸に規定し、Ｚ軸に対して鉛直に直交する軸をＹ軸と呼び、Ｚ軸に水平に直交する軸をＸ軸と呼ぶことにする。Ｘ軸、Ｙ軸及びＺ軸は、直交３次元座標系を構成する。
【００１１】
傾斜磁場コイル４３は、静磁場磁石４１の内側に取り付けられ、中空の略円筒形状に形成されたコイルユニットである。傾斜磁場コイル４３は、傾斜磁場電源２１からの電流の供給を受けて傾斜磁場を発生する。より詳細には、傾斜磁場コイル４３は、互いに直交するＸ軸、Ｙ軸、Ｚ軸に対応する３つのコイルを有する。当該三つのコイルは、Ｘ軸、Ｙ軸、Ｚ軸の各軸に沿って磁場強度が変化する傾斜磁場を形成する。Ｘ軸、Ｙ軸、Ｚ軸の各軸に沿う傾斜磁場は合成されて互いに直交するスライス選択傾斜磁場Ｇｓ、位相エンコード傾斜磁場Ｇｐ及びリードアウト傾斜磁場Ｇｒが所望の方向に形成される。これら傾斜磁場は、静磁場に重畳されて被検体Ｐに印加される。スライス選択傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード傾斜磁場Ｇｐは、空間的位置に応じてＭＲ信号の位相を変化させるために利用される。リードアウト傾斜磁場Ｇｒは、空間的位置に応じてＭＲ信号の周波数を変化させるために利用される。なお、以下の説明においてスライス選択傾斜磁場Ｇｓの傾斜方向はＺ軸、位相エンコード傾斜磁場Ｇｐの傾斜方向はＹ軸、リードアウト傾斜磁場Ｇｒの傾斜方向はＸ軸であるとする。
【００１２】
傾斜磁場電源２１は、シーケンス制御回路２９からのシーケンス制御信号に従い傾斜磁場コイル４３に電流を供給する。傾斜磁場電源２１は、傾斜磁場コイル４３に電流を供給する。傾斜磁場電源２１による電流の供給を受けて、Ｘ軸、Ｙ軸及びＺ軸の各軸に沿う傾斜磁場が傾斜磁場コイル４３により発生される。当該傾斜磁場は、静磁場磁石４１により形成された静磁場に重畳されて被検体Ｐに印加される。
【００１３】
送信コイル４５は、例えば、傾斜磁場コイル４３の内側に配置され、送信回路２３から高周波電流パルス（以下、ＲＦ電流パルスと呼ぶ）の供給を受けて高周波磁場パルス（以下、ＲＦ磁場パルスと呼ぶ）を発生する。
【００１４】
送信回路２３は、被検体Ｐ内に存在する対象プロトンを励起するためのＲＦ磁場パルスを、送信コイル４５を介して被検体Ｐに印加する。対象プロトンとしては、典型的には、プロトンが用いられる。ＲＦ磁場パルスは、対象プロトンに固有の共鳴周波数で振動し、対象プロトンを励起させる。励起された対象プロトンから磁気共鳴信号（以下、ＭＲ信号と呼ぶ）が発生され、受信コイル４７により検出される。送信コイル２３は、例えば、全身用コイル（ＷＢコイル）である。後述のように、全身用コイルは、送受信コイルとして使用されても良い。
【００１５】
受信コイル４７は、ＲＦ磁場パルスの作用を受けて被検体Ｐ内に存在する対象プロトンから発せられるＭＲ信号を受信する。受信コイル４７は、ＭＲ信号を受信可能な複数の受信コイルエレメントを有する。受信されたＭＲ信号は、有線又は無線を介して受信回路２５に供給される。図１に図示しないが、受信コイル４７は、並列的に実装された複数の受信チャネルを有している。受信チャネルは、ＭＲ信号を受信する受信コイルエレメント及びＭＲ信号を増幅する増幅器等を有している。ＭＲ信号は、受信チャネル毎に出力される。受信チャネルの総数と受信コイルエレメントの総数とは同一であっても良いし、受信チャネルの総数が受信コイルエレメントの総数に比して多くても良いし、少なくても良い。
【００１６】
受信回路２５は、励起された対象プロトンから発生されるＭＲ信号を受信コイル４７を介して受信する。受信回路２５は、受信されたＭＲ信号を信号処理してデジタルのＭＲ信号を発生する。デジタルのＭＲ信号を生データと呼ぶことにする。生データは、有線又は無線を介してホストＰＣ５０に供給される。
【００１７】
なお、上記の送信コイル４５と受信コイル４７とは一例に過ぎない。送信コイル４５と受信コイル４７との代わりに、送信機能と受信機能とを備えた送受信コイルが用いられても良い。また、送信コイル４５、受信コイル４７及び送受信コイルが組み合わされても良い。
【００１８】
架台１１に隣接して寝台１３が設置される。寝台１３は、天板１３１と基台１３３とを有する。天板１３１には被検体Ｐが載置される。基台１３３は、天板１３１をＸ軸、Ｙ軸、Ｚ軸各々に沿ってスライド可能に支持する。基台１３３には寝台駆動装置２７が収容される。寝台駆動装置２７は、シーケンス制御回路２９からの制御を受けて天板１３１を移動する。寝台駆動装置２７としては、例えば、サーボモータやステッピングモータ等の如何なるモータが用いられても良い。
【００１９】
シーケンス制御回路２９は、ハードウェア資源として、ＣＰＵ（Central Processing Unit）あるいはＭＰＵ（Micro Processing Unit）のプロセッサとＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）等のメモリとを有する。シーケンス制御回路２９は、通信ＩＦ６１を介してホストＰＣ５０から供給されるパルスシーケンス情報に基づいて傾斜磁場電源２１、送信回路２３及び受信回路２５を同期的に制御し、当該パルスシーケンス情報に応じたパルスシーケンスで被検体Ｐに関するデータ収集を実行する。例えば、本実施形態に係るシーケンス制御回路２９は、複数の受信チャネルを利用して、位相エンコードステップを所定の倍速率に応じて間引いて生データを収集する、パラレルイメージング法のためのデータ収集を実行する。本実施形態に係るパラレルイメージングにおいては、同一位相エンコードにおいて複数の積算回数分のデータ収集が行われる。積算回数は、データ収集回数（Number Of Acquisition）や励起回数（Number Of Excitation）とも呼ばれる。本実施形態に係るパラレルイメージング法のデータ収集により、複数の積算回数にそれぞれ対応する生データセットが、複数の受信チャネル毎に収集されることとなる。これら複数の積算回数且つ複数の受信チャネル毎の生データセットは、記憶回路５３に記憶される。なお、生データセットは、当該所定の倍速率に応じて間引いてｋ空間に充填される一揃いの生データを意味する。
【００２０】
図１に示すように、ホストＰＣ５０は、処理回路５１、記憶回路５３、表示回路５５、入力回路５７、ネットワークＩＦ５９及び通信ＩＦ６１を有するコンピュータ装置である。
【００２１】
処理回路５１は、ハードウェア資源として、ＣＰＵやＧＰＵ（Graphical processing unit）、ＭＰＵ等のプロセッサとＲＯＭやＲＡＭ等のメモリとを有する。処理回路５１は、各種プログラムの実行により再構成機能５１１、積算機能５１３、展開機能５１５、画像処理機能５１７及びシステム制御機能５１９を有する。なお処理回路５１は、再構成機能５１１、積算機能５１３、展開機能５１５、画像処理機能５１７及びシステム制御機能５１９を実現する特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）やフィールド・プログラマブル・ゲート・アレイ（Field Programmable Gate Array：ＦＰＧＡ）、他の複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）により実現されても良い。再構成機能５１１、積算機能５１３、展開機能５１５、画像処理機能５１７及びシステム制御機能５１９は、単一の基板に実装されても良いし、複数の基板に分散されても良い。
【００２２】
再構成機能５１１において処理回路５１は、生データセットに再構成処理を実行して画像を再構成する。再構成処理としては、具体的には、ＦＦＴ（Fast Fourier Transform）等が用いられる。
【００２３】
積算機能５１３において処理回路５１は、複数の受信チャネル各々について積算回数分の複数の画像に積算処理を施して、当該複数の受信チャネル毎に画像を生成する。以下、生成機能５１３において生成された画像を積算画像と呼ぶことにする。
【００２４】
展開機能５１５において処理回路５１は、複数の受信チャネルに関する複数の感度分布画像を利用して、複数の積算画像にパラレルイメージング法の展開（Unfolding）処理を施して画像を生成する。感度分布画像は、各受信チャネルの感度の空間分布を示す画像であり、各画素の画素値は、全身用コイル画像の画素の画素値に対する当該受信コイルエレメント画像の画素の画素値の比を示す。全身用コイル画像は、全身用コイルを用いた感度分布測定スキャンにより収集された生データに基づいて再構成機能５１１により生成される。同様に、各受信コイルエレメント画像は、当該受信チャネルを用いた感度分布測定スキャンにより収集された生データに基づいて再構成機能５１１により生成される。以下、展開機能５１５により生成された画像を最終画像と呼ぶことにする。
【００２５】
画像処理機能５１７において処理回路５１は、最終画像等の画像に種々の画像処理を施す。例えば、処理回路５１は、ボリュームレンダリングや、サーフェスレンダリング、画素値投影処理、ＭＰＲ（Multi-Planer Reconstruction）処理、ＣＰＲ（Curved MPR）処理等の画像処理を施す。また、処理回路５１は、最終画像にテンソル解析を施して拡散強調パラメータの空間分布を示す機能画像を生成しても良い。
【００２６】
システム制御機能５１９において処理回路５１は、本実施形態に係る磁気共鳴イメージング装置１の全体を制御する。
【００２７】
記憶回路５３は、種々の情報を記憶するＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）、集積回路記憶装置等の記憶装置である。また、記憶回路５３は、ＣＤ−ＲＯＭドライブやＤＶＤドライブ、フラッシュメモリ等の可搬性記憶媒体との間で種々の情報を読み書きする駆動装置等であっても良い。例えば、記憶回路５３は、複数の受信コイルエレメントに関する複数の感度分布画像を記憶する。また、記憶回路５３は、各種スキャンにより収集された生データセットやプログラム等を記憶する。
【００２８】
表示回路５５は、種々の情報を表示する。例えば、表示回路５５は、展開機能５１５により生成された最終画像や、画像処理機能５１７により画像処理された画像等を表示する。表示回路５５は、表示インタフェースと表示機器とを有する。表示インタフェースは、表示対象を表すデータを映像信号に変換する。映像信号は、表示機器に供給される。表示機器は、表示対象を表す映像信号を表示する。表示機器としては、例えば、ＣＲＴディスプレイや液晶ディスプレイ、有機ＥＬディスプレイ、ＬＥＤディスプレイ、プラズマディスプレイ、又は当技術分野で知られている他の任意のディスプレイが適宜利用可能である。
【００２９】
入力回路５７は、具体的には、入力機器と入力インタフェースとを有する。入力機器は、ユーザからの各種指令を受け付ける。入力機器としては、キーボードやマウス、各種スイッチ、タッチスクリーン、タッチパッド等が利用可能である。入力インタフェースは、入力機器からの出力信号をバスを介して処理回路５１に供給する。なお、入力回路５７は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限らない。例えば、磁気共鳴イメージング装置１とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、受け取った電気信号を種々の回路へ出力するような電気信号の処理回路も入力回路の例に含まれる。
【００３０】
ネットワークＩＦ５９は、ＬＡＮ（Local Area Network）等を介して磁気共鳴イメージング装置１と、ワークステーションやＰＡＣＳ（Picture Archiving and Communication System）、ＨＩＳ（Hospital Information System）、ＲＩＳ（Radiology Information System）等とを接続するインタフェースである。ネットワークＩＦは、各種情報を接続先のワークステーション、ＰＡＣＳ、ＨＩＳ及びＲＩＳとの間で送受信する。
【００３１】
通信ＩＦ６１は、有線又は無線を介してホストＰＣ５０にシーケンス制御回路２９と受信回路２５とを接続するインタフェースである。例えば、通信ＩＦ６１は、シーケンス制御回路２９にパルスシーケンス情報を送信する。また、通信ＩＦ６１は、受信回路２５から生データを受信する。
【００３２】
なお、上記の構成は一例であって、これに限定されない。例えば、シーケンス制御回路２９は、ホストＰＣ５０に組み込まれても良い。また、シーケンス制御回路２９と処理回路５１とが同一の基板に実装されても良い。シーケンス制御回路２９、傾斜磁場電源２１、送信回路２３及び受信回路２５は、ホストＰＣ５０とは異なる単一の制御装置に実装されても良いし、複数の装置に分散して実装されても良い。
【００３３】
次に、本実施形態に係る磁気共鳴イメージング装置の動作例について説明する。
【００３４】
図２は、本実施形態に係る簡易的なパルスシーケンスを示す図である。図２に示すように、本実施形態に係るシーケンス制御回路２９は、傾斜磁場電源２１、送信回路２３及び受信回路２５を同期的に制御し、被検体Ｐの同一撮像部位を対象として複数のショット（データ収集）を実行する。ショットの回数Ｎは、積算回数に一致し、２以上であれば幾つでも良い。各ショットにおいてシーケンス制御回路２９は、傾斜磁場電源２１、送信回路２３及び受信回路２５を同期的に制御し、複数の受信チャネルを利用して、位相エンコードステップを所定の倍速率に応じて間引いて生データを収集する。これにより、積算回数Ｎ（Ｎは整数）と受信チャネル数Ｍ（Ｍは整数）との組合せ毎の生データセットが収集される。なお、積算回数と倍速率とは、例えば、撮像プロトコルに応じて自動的に又はユーザにより入力回路５７を介した指示に従い手動的に設定される。
【００３５】
本実施形態に係るショットの撮像法は、特に限定されないが、以下の説明を具体的に行うため、最終画像の生成処理に有用な拡散強調撮像、すなわち、ＤＷＩ（Diffusion Weighted Imaging）収集であるものとする。ＤＷＩ収集は、組織内の水の拡散を測定するために使用される。ＤＷＩ収集においては、ＭＰＧパルスによる位相分散の影響を除去するため複数のショットが実行される。
【００３６】
図３は、１ショット分のＤＷＩ収集のパルスシーケンスの一例を示す図である。図３のＤＷＩ収集のパルスシーケンスは、エコープラナー（ＥＰＩ：echo-planar imaging）法を基礎にする。図３に示すＲＦのタイムラインは、送信回路２３により印加されるＲＦ磁場パルスの磁場強度の時系列を示し、Ｇｓのタイムラインは、傾斜磁場電源２１により印加されるスライス選択傾斜磁場の磁場強度の時系列を示し、Ｇｐのタイムラインは、傾斜磁場電源２１により印加される位相エンコード傾斜磁場の磁場強度の時系列を示し、Ｇｒのタイムラインは、傾斜磁場電源２１により印加されるリードアウト用傾斜磁場の磁場強度の時系列を示す。
【００３７】
図３に示すように、送信回路２３により９０°パルスと１８０°パルスとが所定の時間間隔をあけて印加される。９０°パルスの印加時において傾斜磁場電源２１によりスライス選択傾斜磁場が印加され、１８０°パルスの印加時においても傾斜磁場電源２１によりスライス選択傾斜磁場が印加される。ＤＷＩ収集において傾斜磁場電源２１は、２つの運動検出勾配磁場（ＭＰＧ：motion-probing gradient）パルスを、所定時間間隔を空けて印加し、拡散された流体の信号を選択的に減衰させる。傾斜磁場電源２１は、９０°パルスと１８０°パルスとの間に第１のＭＰＧパルスを印加する。第１のＭＰＧパルスの印加により、プロトンのスピンの位相分散が生じる。傾斜磁場電源２１は、１８０°パルスの印加後、第２のＭＰＧパルスを印加する。第２のＭＰＧパルスは、第１のＭＰＧパルスと同じ強度を有する。第１のＭＰＧパルスと第２のＭＰＧパルスとは、スライス選択傾斜磁場方向、位相エンコード傾斜磁場方向及びリードアウト傾斜磁場方向に印加される。なお、第１のＭＰＧパルスと第２のＭＰＧパルスとは、スライス選択傾斜磁場方向、位相エンコード傾斜磁場方向及びリードアウト傾斜磁場方向の全方向に印加される必要はなく、１方向又は２方向に印加されても良い。
【００３８】
第１のＭＰＧパルスの印加と第２のＭＰＧパルスの印加との間において水分子が動かなかった場合、第２のパルスは、第１のＭＰＧパルスにより生じた位相分散を補償する。水分子が拡散した場合、位相分散を完全に補償することはできない。ＭＰＧパルスの印加間において拡散された水は、拡散しない水よりも弱いＭＲ信号を発生する。ＭＲ信号の減衰は拡散の量に比例する。拡散強調の程度はｂ値により表される。ｂ値は、ｂ（s/mm2）＝γ２Ｇｘ２δ２（Δ−δ／３）により規定され、ＭＰＧパルスの振幅Ｇｘ、ＭＰＧパルスの印加時間δ、２つのＭＰＧパルス間の時間間隔Δに依存する。なお、γは磁気回転比を示す。
【００３９】
図３に示すように、１８０°パルスの印加から所定時間経過後、受信回路２５は、受信コイル４７を介してＭＲ信号を受信する。この際、傾斜磁場電源２１は、正極性のリードアウト傾斜磁場と負極性のリードアウト傾斜磁場とを傾斜磁場コイル４３から交互に高速に印加し、正極性のリードアウト傾斜磁場と負極性のリードアウト傾斜磁場との間に位相エンコード傾斜磁場を高速に印加する。この際、倍速率に応じて位相エンコードステップが間引かれる。このようにして、一回の励起によりｋ空間を充填するのに必要な生データセットが収集される。
【００４０】
シーケンス制御回路２９は、上記１ショット分のＤＷＩ収集を積算回数分だけ繰り返すことにより、積算回数Ｎと受信チャネル数Ｍとの組合せ毎の生データセットを収集する。すなわち、Ｎ×Ｍの生データセットが収集される。なお、傾斜磁場電源２１は、ショット毎に、ｂ値及び印加方向の少なくとも一方を変えてＤＷＩパルスを印加すると良い。ＤＷＩパルスの印加方向は、Ｘ軸、Ｙ軸、及びＺ軸の各々のＭＰＧパルスの強度のベクトル和により決定される。
【００４１】
上記の説明において図３のＤＷＩ収集のパルスシーケンスは、エコープラナー（ＥＰＩ：echo-planar imaging）法を基礎にするものとしたが、本実施形態に係るＤＷＩ収集のパルスシーケンスは、これに限定されず、ＳＥ（Spin Echo）法やＧＲＥ（Gradient Echo）法等の既存の如何なるパルスシーケンスを基礎にしても良い。
【００４２】
次に、処理回路５１により行われる最終画像の生成処理について説明する。
【００４３】
図４は、処理回路５１により行われる最終画像の生成処理の過程を模式的に示す図である。図４に示すように、上記のＤＷＩ収集により、Ｍ個の受信チャネル毎に、Ｎショット分の生データセットが収集され、Ｎ×Ｍ個の生データセットが記憶回路５３に記憶される。なお、最終画像の生成処理は、シーケンス制御回路２９によるＤＷＩ収集の終了を契機として自動的に開始されても良いし、ユーザにより入力回路５７を介して開始指示がなされた事を契機として開始されても良い。
【００４４】
まず、処理回路５１は、再構成機能５１１を実行する（ステップＳ１）。ステップＳ１において処理回路５１は、Ｎ×Ｍ個の生データセット各々にＦＦＴを施して折り返し画像を生成する。折り返し画像は、パラレルイメージングのＦＯＶ（以下、撮像ＦＯＶと呼ぶ）に比して、当該パラレルイメージング法の倍速率に応じて縮小されたＦＯＶを有する。当該画像は、倍速率に応じた位相エンコードステップの間引きに起因する折り返しアーチファクトを有している。
【００４５】
ステップＳ１が行われると処理回路５１は、積算機能５１３を実行する（ステップＳ２）。ステップＳ２において処理回路５１は、Ｎ枚の折り返し画像をＭ個の受信チャネル毎に積算して、受信チャネルと同数であるＭ個の積算画像を生成する。すなわち、処理回路５１は、各ショットについて、Ｎ枚の折り返し画像に積算処理を施して積算画像を生成する。処理回路５１は、例えば、Ｎ枚の折り返し画像の同一座標の画素値を加算することにより積算画像を生成する。なお、積算画像の画素値は、Ｎ枚の折り返し画像の同一座標の画素値の合計値であることに限定されず、当該合計値に基づく統計値であっても良い。統計値としては、Ｎ枚の折り返し画像の同一座標の画素値の平均値、中間値、最小値及び最大値等が計算される。積算画像のＦＯＶは、折り返し画像のＦＯＶと同一である。積算画像には、折り返し画像に含まれていた折り返しアーチファクトが残存する。
【００４６】
ステップＳ２が行われると処理回路５１は、展開機能５１５を実行する（ステップＳ３）。ステップＳ３において処理回路５１は、Ｍ個の受信チャネルに関するＭ枚の感度分布画像を利用して、ステップＳ２において生成されたＭ枚の積算画像に展開処理を施して１枚の最終画像を生成する。最終画像のＦＯＶは、縮小されておらず、撮像ＦＯＶに略一致する。最終画像の折り返しアーチファクトは低減されている。生成された最終画像は、表示回路５５に表示される。また、画像処理機機能５１７により、処理回路５１は、最終画像にテンソル解析を施して拡散強調パラメータの空間分布を示す機能画像を生成しても良い。拡散強調パラメータとしては、平均ＡＤＣ（apparent diffusion coefficient）や各種の異方性指標、例えば、ＦＡ（fractional anisotropy）、ＲＡ（relative anisotropy）、ＶＲ（volume ratio）等が知られている。機能画像は表示回路５５に表示される。
【００４７】
次に、本実施形態に係る最終画像の生成処理の効果について説明する。
【００４８】
図５は、以上で説明した実施形態とは異なる、最終画像の生成処理の過程の比較例を模式的に示す図である。図５に示すように、比較例においては、各生データセットにＦＦＴ等の再構成処理を施してＮ×Ｍ毎の折り返し画像を生成し（ステップＳＺ１）、Ｎ回のショット各々についてＭ枚の折り返し画像に展開処理を施してＮ枚の展開画像を生成し（ステップＳＺ２）、その後、Ｎ枚の展開画像に積算処理を施して単一の積算画像（最終画像）を生成している（ステップＳＺ３）。展開処理においてはＳＮＲが高い事が望まれる。比較例においては、展開処理の入力画像が積算前の折り返し画像であるので、積算回数（ショット回数）を増やしても、展開画像及び最終画像のＳＮＲは上がらない。
【００４９】
本実施形態によれば、図４に示すように、積算処理（Ｓ２）を施した後に展開処理（Ｓ３）が行われる。従って、展開処理の入力画像が積算後の折り返し画像であるので、積算回数（ショット回数）を増やすことにより、最終画像のＳＮＲを上げることができる。従って本実施形態によれば、ＤＷＩ収集において積算回数を増やすことにより、パラレルイメージング法の展開処理の展開誤差の少ない画像を生成することができる。
【００５０】
以上、上記の少なくとも一つの実施形態によれば、積算処理を行うパラレルイメージング法において展開誤差を低減することが可能になる。
【００５１】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００５２】
１…磁気共鳴イメージング装置、１１…架台、１３…寝台、２１…傾斜磁場電源、２３…送信回路、２５…受信回路、２７…寝台駆動装置、２９…シーケンス制御回路、４１…静磁場磁石、４３…傾斜磁場コイル、４５…送信コイル、４７…受信コイル、５０…ホストＰＣ、５１…処理回路、５３…記憶回路、５５表示回路、５７…入力回路、１３１…天板、１３３…基台、５１１…再構成機能、５１３…積算機能、５１５…展開機能、５１７…画像処理機能、５１９…システム制御機能。
【図１】
2019005136.tif 000003
【図２】
2019005136.tif 000004
【図３】
2019005136.tif 000005
【図４】
2019005136.tif 000006
【図５】
2019005136.tif 000007