JP2023049144A - 磁気共鳴イメージング装置、画像生成方法、および画像生成プログラム - Google Patents

Abstract

【課題】急峻な位相の変化に伴って発生するアーチファクトを低減すること。【解決手段】実施形態に係る磁気共鳴イメージング装置は、取得部と、生成部と、を備える。取得部は、複数のコイルに関する第１コイル感度マップと、前記第１コイル感度マップとは位相が異なり前記複数のコイルに関する第２コイル感度マップとを取得する。生成部は、前記第１コイル感度マップと前記複数のコイルに関する磁気共鳴データとに基づいて第１画像を生成し、前記第２コイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて第２画像を生成し、前記第１画像と前記第２画像とに基づいて、最終画像を生成する。【選択図】図３

Description

本明細書及び図面に開示の実施形態は、磁気共鳴イメージング装置、画像生成方法、および画像生成プログラムに関する。

磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ：以下、ＭＲＩと呼ぶ）装置において、撮像時間の短縮のために、間引き収集が行われることがある。このとき、間引き収集によって得られた間引きデータからコイル感度のマップを推定することがある。推定されたコイル感度マップを用いて再構成を行った場合、再構成された画像の位相に、滑らかに変化しない個所が発生することがある。

このため、推定されたコイル感度マップを用いて再構成された画像に対して補間処理が実行されると、補間処理後の画像において、位相が大きく変化している部分においてアーチファクトが発生する。

特開２０１３－２４０５７１号公報

本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、急峻な位相の変化に伴って発生するアーチファクトを低減することにある。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。

実施形態に係る磁気共鳴イメージング装置は、取得部と、生成部と、を備える。取得部は、複数のコイルに関する第１コイル感度マップと、前記第１コイル感度マップとは位相が異なり前記複数のコイルに関する第２コイル感度マップとを取得する。生成部は、前記第１コイル感度マップと前記複数のコイルに関する磁気共鳴データとに基づいて第１画像を生成し、前記第２コイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて第２画像を生成し、前記第１画像と前記第２画像とに基づいて、最終画像を生成する。

図１は、実施形態に係る磁気共鳴イメージング装置の一例を示すブロック図。 図２は、実施形態に係る画像生成処理の手順の一例を示すフローチャート。 図３は、実施形態に係り、画像生成処理における各種処理と当該各種処理に関するデータとの関連の一例を示す図。 図４は、実施形態に係る推定感度強度画像の一例を示す図。 図５は、実施形態に係り、第１コイル感度マップにおける複数の位相画像の一例を示す図。 図６は、実施形態に係り、複数のコイル画像における強度画像の一例を示す図。 図７は、実施形態に係り、複数のコイル画像における位相画像の一例を示す図。 図８は、実施形態に係る第１画像の一例を示す図。 図９は、実施形態に係り、複数のコイルに対応する複数の分解画像における強度画像の一例を示す図。 図１０は、実施形態に係り、複数のコイルに対応する複数の分解画像における位相画像の一例を示す図。 図１１は、実施形態に係り、第２コイル感度マップにおける複数の位相画像の一例を示す図。 図１２は、実施形態に係り、第３コイル感度マップにおける複数の位相画像の一例を示す図。 図１３は、実施形態に係り、コイル合成された第２画像の一例を示す図。 図１４は、実施形態に係り、コイル合成された第３画像の一例を示す図。 図１５は、実施形態に係る第１補間画像の一例を示す図。 図１６は、実施形態に係る第２補間画像の一例を示す図。 図１７は、実施形態に係る第３補間画像の一例を示す図。 図１８は、実施形態に係る最終画像の一例を示す図。

以下、図面を参照しながら、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ：以下、ＭＲＩと呼ぶ）装置、画像生成方法、および画像生成プログラムの実施形態について詳細に説明する。なお、本実施形態の技術的思想は、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ：陽電子放出コンピュータ断層撮像）－ＭＲＩ装置、ＳＰＥＣＴ（ｓｉｎｇｌｅ ｐｈｏｔｏｎ ｅｍｉｓｓｉｏｎ ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ：単一光子放出コンピュータ断層撮像）－ＭＲＩ装置などのＭＲＩ装置と複合的な各種モダリティに適用されてもよい。

（第１の実施形態）
図１は、本実施形態に係るＭＲＩ装置１００の一例を示す図である。図１に示すように、ＭＲＩ装置１００は、静磁場磁石１０１と、傾斜磁場コイル１０３と、傾斜磁場電源１０５と、寝台１０７と、寝台制御回路１０９と、送信回路１１３と、送信コイル１１５と、受信コイル１１７と、受信回路１１９と、撮像制御回路（撮像制御部）１２１と、システム制御回路（システム制御部）１２３と、メモリ１２５と、入力インターフェース１２７と、ディスプレイ１２９と、処理回路１３１と、を備える。

静磁場磁石１０１は、中空の略円筒状に形成された磁石である。静磁場磁石１０１は、内部の空間に略一様な静磁場を発生する。静磁場磁石１０１としては、例えば、超伝導磁石等が使用される。

傾斜磁場コイル１０３は、中空の略円筒形状に形成されたコイルであり、円筒形の冷却容器の内面側に配置される。傾斜磁場コイル１０３は、傾斜磁場電源１０５から個別に電流供給を受けて、互いに直交するＸ、Ｙ、及びＺの各軸に沿って磁場強度が変化する傾斜磁場を発生する。傾斜磁場コイル１０３によって発生されるＸ、Ｙ、Ｚ各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場、位相エンコード用傾斜磁場および周波数エンコード用傾斜磁場を形成する。スライス選択用傾斜磁場は、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場は、空間的位置に応じて磁気共鳴信号（以下、ＭＲ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ）信号と呼ぶ）の位相を変化させるために利用される。周波数エンコード用傾斜磁場は、空間的位置に応じてＭＲ信号の周波数を変化させるために利用される。

傾斜磁場電源１０５は、撮像制御回路１２１の制御により、傾斜磁場コイル１０３に電流を供給する電源装置である。

寝台１０７は、被検体Ｐが載置される天板１０７１を備えた装置である。寝台１０７は、寝台制御回路１０９による制御のもと、被検体Ｐが載置された天板１０７１を、ボア１１１内へ挿入する。

寝台制御回路１０９は、寝台１０７を制御する回路である。寝台制御回路１０９は、入出力インターフェース１７を介した操作者の指示により寝台１０７を駆動することで、天板１０７１を長手方向および上下方向、場合によっては左右方向へ移動させる。

送信回路１１３は、撮像制御回路１２１の制御により、ラーモア周波数で変調された高周波パルスを送信コイル１１５に供給する。例えば、送信回路１１３は、発振部や位相選択部、周波数変換部、振幅変調部、ＲＦ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ）アンプなどを有する。発振部は、静磁場中における対象原子核に固有の共鳴周波数のＲＦパルスを発生する。位相選択部は、発振部によって発生したＲＦパルスの位相を選択する。周波数変換部は、位相選択部から出力されたＲＦパルスの周波数を変換する。振幅変調部は、周波数変換部から出力されたＲＦパルスの振幅を例えばｓｉｎｃ関数に従って変調する。ＲＦアンプは、振幅変調部から出力されたＲＦパルスを増幅して送信コイル１１５に供給する。

送信コイル１１５は、傾斜磁場コイル１０３の内側に配置されたＲＦコイルである。送信コイル１１５は、送信回路１１３からの出力に応じて、高周波磁場に相当するＲＦパルスを発生する。

受信コイル１１７は、傾斜磁場コイル１０３の内側に配置されたＲＦコイルである。受信コイル１１７は、高周波磁場によって被検体Ｐから放射されるＭＲ信号を受信する。受信コイル１１７は、受信されたＭＲ信号を受信回路１１９へ出力する。受信コイル１１７は、例えば、１以上、典型的には複数のコイルエレメント（以下、複数のコイルと呼ぶ）を有するコイルアレイである。以下、説明を具体的にするために、受信コイル１１７は、複数のコイルを有するコイルアレイとして説明する。

なお、図１において送信コイル１１５と受信コイル１１７とは別個のＲＦコイルとして記載されているが、送信コイル１１５と受信コイル１１７とは、一体化された送受信コイルとして実施されてもよい。送受信コイルは、被検体Ｐの撮像部位に対応し、例えば、頭部コイルのような局所的な送受信ＲＦコイルである。

受信回路１１９は、撮像制御回路１２１の制御により、受信コイル１１７から出力されたＭＲ信号に基づいて、デジタルのＭＲ信号（以下、ＭＲデータと呼ぶ）を生成する。具体的には、受信回路１１９は、受信コイル１１７から出力されたＭＲ信号に対して、検波、フィルタリングなどの信号処理を施した後、当該信号処理が施されたデータに対してアナログ／デジタル（Ａ／Ｄ（Ａｎａｌｏｇ ｔｏ Ｄｉｇｉｔａｌ））変換（以下、Ａ／Ｄ変換と呼ぶ）して、ＭＲデータを生成する。受信回路１１９は、生成されたＭＲデータを、撮像制御回路１２１に出力する。例えば、ＭＲデータは、複数のコイル各々において生成され、複数のコイル各々を識別するタグとともに、撮像制御回路１２１に出力される。

撮像制御回路１２１は、処理回路１５から出力された撮像プロトコルに従って、傾斜磁場電源１０５、送信回路１１３及び受信回路１１９等を制御し、被検体Ｐに対する撮像を行う。撮像プロトコルは、検査の種類に応じたパルスシーケンスを有する。撮像プロトコルには、傾斜磁場電源１０５により傾斜磁場コイル１０３に供給される電流の大きさ、傾斜磁場電源１０５により電流が傾斜磁場コイル１０３に供給されるタイミング、送信回路１１３により送信コイル１１５に供給される高周波パルスの大きさや時間幅、送信回路１１３により送信コイル１１５に高周波パルスが供給されるタイミング、受信コイル１１７によりＭＲ信号が受信されるタイミング等が定義されている。撮像制御回路１２１は、傾斜磁場電源１０５、送信回路１１３及び受信回路１１９等を駆動して被検体Ｐを撮像した結果、受信回路１１９からＭＲデータを受信すると、受信したＭＲデータを処理回路１３１へ転送する。

撮像制御回路１２１は、例えば、間引き（アンダーサンプリング（ｕｎｄｅｒ－ｓａｍｐｌｉｎｇ））収集での撮像（以下、間引き撮像と呼ぶ）に関するパルスシーケンスを実行するスキャンにより、ＭＲデータを収集する。以下、説明を具体的にするために、被検体Ｐに対するスキャンは、ｋ空間において自動校正信号（オートキャリブレーションシグナル（Ａｕｔｏ－Ｃａｌｉｂｒａｔｅｄ Ｓｉｇｎａｌ）：以下、ＡＣＳと呼ぶ））に対応するトラジェクトリーに関するＭＲ信号を収集する撮像プロトコル（以下、ＡＣＳプロトコルと呼ぶ）を有するパラレルイメージングであるものとする。例えば、本実施形態におけるスキャンに関するパラレルイメージングは、ＧＲＡＰＰＡ（ｇｅｎｅｒａｌｉｚｅｄ ａｕｔｏｃａｌｉｂｒａｔｉｎｇ ｐａｒｔｉａｌｌｙ ｐａｒａｌｌｅｌ ａｃｑｕｉｓｉｔｉｏｎｓ）に類する撮像プロトコルを有する。

なお、当該パラレルイメージングにおけるＡＣＳに対応するトラジェクトリーは、ＧＲＡＰＰＡのようにｋ空間の中心付近に限定されず、任意に設定可能である。また、被検体Ｐに対するスキャンは、ＡＣＳプロトコルを有すれば、上記パラレルイメージングに限定されず、圧縮センシング（ｃｏｍｐｒｅｓｓｅｄ ｓｅｎｓｉｎｇ）に類する撮像プロトコルにより実現されてもよい。また、被検体Ｐに対するスキャンは、ＡＣＳプロトコルを有すれば、圧縮センシングとパラレルイメージングとを組み合わせた撮像プロトコルにより実現されてもよい。

複数のコイル各々においてパラレルイメージングにより収集されたＭＲデータは、ユーザの指示等により予め設定された間引き率（Ｒｅｄａｃｔｉｏｎ ｆａｃｔｏｒ：Ｒ、加速率（ａｃｃｅｌｅｒａｔｉｏｎ ｆａｃｔｏｒ）ともいう）に従ってｋ空間において位相エンコードステップを間引いて収集されたデータと、ＡＣＳと、を有する。ＡＣＳは、データの間引きの位置に関連するターゲット点（ｔａｒｇｅｔ ｐｏｉｎｔ）におけるターゲットデータ（ｔａｒｇｅｔ ｄａｔａ）とターゲット点の近隣のソース点（ｓｏｕｒｃｅ ｐｏｉｎｔ）におけるソースデータ（ｓｏｕｒｃｅ ｄａｔａ）と、を有する。

なお、撮像制御回路１２１は、被検体Ｐの撮像に用いられる受信コイル１１７の感度の分布を示す画像の生成に関するＭＲデータを、任意の撮像手法により収集してもよい。コイルの感度を示す画像は、複素数のデータで表現される。受信コイル１１７の感度の分布を示す画像の生成に関するＭＲデータの収集は、例えば、被検体Ｐに対するスキャンに先立って、ロケータスキャンなどを含むプリスキャンにおいて撮像制御回路１２１により実行される。撮像制御回路１２１は、例えばプロセッサにより実現される。

「プロセッサ」という文言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。

システム制御回路１２３は、ハードウェア資源として図示していないプロセッサ、ＲＯＭ（Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）やＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）等のメモリ等を有し、システム制御機能によりＭＲＩ装置１００を制御する。具体的には、システム制御回路１２３は、メモリに記憶されたシステム制御プログラムを読み出してメモリ上に展開し、展開されたシステム制御プログラムに従って本ＭＲＩ装置１００の各回路を制御する。

例えば、システム制御回路１２３は、入力インターフェース１２７を介して操作者から入力された撮像条件に基づいて、撮像プロトコルをメモリ１２５から読み出す。システム制御回路１２３は、撮像プロトコルを撮像制御回路１２１に送信し、被検体Ｐに対する撮像を制御する。システム制御回路１２３は、例えばプロセッサにより実現される。なお、システム制御回路１２３は、処理回路１５に組み込まれてもよい。このとき、システム制御機能は処理回路１５により実行され、処理回路１５は、システム制御回路１２３の代替として機能する。システム制御回路１２３を実現するプロセッサは、上述と同様な内容なため、説明は省略する。

メモリ１２５は、システム制御回路１２３において実行されるシステム制御機能に関する各種プログラム、各種撮像プロトコル、撮像プロトコルを規定する複数の撮像パラメータを含む撮像条件等を記憶する。また、メモリ１２５は、処理回路１５により実現される取得機能３３、特定機能３５、および生成機能３７を、コンピュータによって実行可能なプログラムの形態で記憶する。

また、メモリ１２５は、取得機能３３により取得された各種データ、特定機能３５および生成機能３７により実施される処理に用いられる各種データ、生成機能３７により生成されたＭＲ画像、当該ＭＲ画像の生成過程において生成される各種データなどを記憶する。メモリ１２５に記憶される各種データのついては、後ほど説明する。また、メモリ１２５は、被検体Ｐに対するスキャンにより取得されたＭＲデータおよび当該ＭＲデータに基づいてＭＲ画像を再構成するアルゴリズムを記憶する。

なお、メモリ１２５は、不図示の通信インターフェースを介して受信された各種データを記憶してもよい。例えば、メモリ１２５は、放射線情報システム（ＲＩＳ：Ｒａｄｉｏｌｏｇｙ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｙｓｔｅｍ）等の医療機関内の情報処理システムから受信した被検体Ｐの検査オーダに関する情報（撮像対象部位、検査目的等）を記憶する。

メモリ１２５は、例えば、ＲＯＭ、ＲＡＭ、フラッシュメモリ等の半導体メモリ素子、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉｖｅ）、光ディスク等により実現される。また、メモリ１２５は、ＣＤ（Ｃｏｍｐａｃｔ Ｄｉｓｃ）－ＲＯＭドライブやＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ）ドライブ、フラッシュメモリ等の可搬型記憶媒体との間で種々の情報を読み書きする駆動装置等で実現されてもよい。

入力インターフェース１２７は、操作者からの各種指示（例えば、電源投入指示）や情報入力を受け付ける。入力インターフェース１２７は、例えば、トラックボール、スイッチボタン、マウス、キーボード、操作面へ触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一体化されたタッチスクリーン、光学センサを用いた非接触入力回路、及び音声入力回路等によって実現される。入力インターフェース１２７は、処理回路１３１に接続されており、操作者から受け取った入力操作を電気信号へ変換し処理回路１３１へと出力する。なお、本明細書において入力インターフェース１２７は、マウス、キーボードなどの物理的な操作部品を備えるものだけに限られない。例えば、ＭＲＩ装置１００とは別体に設けられた外部の入力機器から入力操作に対応する電気信号を受け取り、この電気信号を制御回路へ出力する電気信号の処理回路も入力インターフェース１２７の例に含まれる。

入力インターフェース１２７は、ディスプレイ１２９に表示されたプリスキャン画像に対して、ユーザの指示によりＦＯＶを入力する。具体的には、入力インターフェース１２７は、ディスプレイ１２９に表示されたロケータ画像において、ユーザによる範囲の設定指示によりＦＯＶを入力する。また、入力インターフェース１２７は、検査オーダに基づくユーザの指示により、スキャンに関する各種撮像パラメータを入力する。

ディスプレイ１２９は、処理回路１３１またはシステム制御回路１２３による制御のもとで、各種のＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）や、処理回路１３１によって生成されたＭＲ画像、生成機能３７により生成された最終画像等を表示する。また、ディスプレイ１２９は、スキャンに関する撮像パラメータ撮像、および画像処理に関する各種情報などを表示する。ディスプレイ１２９は、例えば、ＣＲＴディスプレイや液晶ディスプレイ、有機ＥＬディスプレイ、ＬＥＤディスプレイ、プラズマディスプレイ、又は当技術分野で知られている他の任意のディスプレイ、モニタ等の表示デバイスにより実現される。

処理回路１３１は、例えば、上述のプロセッサなどにより実現される。処理回路１３１は、取得機能３３、特定機能３５、および生成機能３７などを備える。取得機能３３、特定機能３５、および生成機能３７をそれぞれ実現する処理回路１３１は、取得部、特定部、および生成部に相当する。取得機能３３、特定機能３５、および生成機能３７などの各機能は、コンピュータによって実行可能なプログラムの形態でメモリ１２５に記憶されている。例えば、処理回路１３１は、プログラムをメモリ１２５から読み出し、実行することで各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の処理回路１３１は、取得機能３３、特定機能３５、および生成機能３７などの各機能を有することとなる。

上記説明では、「プロセッサ」が各機能に対応するプログラムをメモリ１２５から読み出して実行する例を説明したが、実施形態はこれに限定されない。プロセッサが例えばＣＰＵである場合、プロセッサはメモリ１２５に保存されたプログラムを読み出して実行することで機能を実現する。一方、プロセッサがＡＳＩＣである場合、メモリ１２５にプログラムを保存する代わりに、当該機能がプロセッサの回路内に論理回路として直接組み込まれる。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。また、単一の記憶回路が各処理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散して配置して、処理回路１３１は個別の記憶回路から対応するプログラムを読み出す構成としても構わない。

処理回路１３１は、取得機能３３により、複数のコイルに関する第１コイル感度マップを取得する。具体的には、取得機能３３は、受信コイル１１７における複数のコイル各々により取得されたＭＲデータに基づいて、複数のコイル各々の感度マップを推定する。より詳細には、取得機能３３は、例えば、複数のコイル各々に関して、ＭＲデータにおけるＡＣＳを用いたＥＳＰＩＲｉＴ（Ｕｅｃｋｅｒ Ｍ， Ｌａｉ Ｐ， Ｍｕｒｐｈｙ ＭＪ， ｅｔ ａｌ． ： ＥＳＰＩＲｉＴ －ａｎ ｅｉｇｅｎｖａｌｕｅ ａｐｐｒｏａｃｈ ｔｏ ａｕｔｏｃａｌｉｂｒａｔｉｎｇ ｐａｒａｌｌｅｌ ＭＲＩ ： ｗｈｅｒｅ ＳＥＮＳＥ ｍｅｅｔｓ ＧＲＡＰＰＡ．Ｍａｇｎ Ｒｅｓｏｎ Ｍｅｄ ２０１４ ； ７１ ： ９９０ １００１）により、感度マップを推定する。これにより、取得機能３３は、間引き収集されたＭＲデータに基づいて、複数のコイルに関する感度マップ（以下、推定感度マップと呼ぶ）を推定する。推定感度マップは、複数のコイルにおいて感度の分布を濃淡値で示す強度画像と当該強度画像に対応する位相画像とを有する。なお、感度マップの推定方法は、上記ＥＳＰＩＲｉＴに限定されず、既知の各種技術が適用可能である。

次いで、処理回路１３１は、取得機能３３により、複数のコイルのうち特定機能３５により特定されたコイル（以下、第１特定コイルと呼ぶ）の推定感度マップにおける位相画像を基準として、複数のコイルに関する推定感度マップの位相をシフトさせる。これにより、取得機能３３は、複数のコイルに関する第１コイル感度マップを取得する。すなわち、第１コイル感度マップは、複数のコイルに関し、第１特定コイルの推定感度マップにおける位相画像（以下、第１特定位相画像と呼ぶ）の位相を一様にするように、第１特定コイルとは異なる他のコイルに関する推定感度マップにおける位相画像の位相をシフトさせた位相画像を有する。換言すれば、第１コイル感度マップにおいて複数のコイルの感度を濃淡値で示す強度画像は、推定感度マップにおける強度画像と同様であるが、第１コイル感度マップにおける強度画像に対応する位相画像は、第１特定位相画像における位相を一様にするように推定感度マップにおける位相画像の位相をシフトさせた位相画像に相当する。特定機能３５による第１特定コイルの特定については、後ほど説明する。

処理回路１３１は、取得機能３３により、第１コイル感度マップにおける位相画像とは位相が異なり複数のコイルに関する第２コイル感度マップを取得する。具体的には、取得機能３３は、第１特定コイルとは異なり、特定機能３５により特定されたコイル（以下、第２特定コイルと呼ぶ）の推定感度マップにおける位相画像を基準として、複数のコイルに関する推定感度マップにおける位相画像の位相をシフトさせる。これにより、取得機能３３は、複数のコイルに亘る第２コイル感度マップを取得する。特定機能３５による第２特定コイルの特定については、後ほど説明する。以上により、取得機能３１は、複数のコイルのうち１つのコイルの感度をそれぞれ基準として複数のコイルの感度の位相画像における位相を変化させて、第１コイル感度マップと第２コイル感度マップとを取得する。すなわち、第１コイル感度マップと第２コイル感度マップとの相違は、基準となる位相画像が異なることにある。換言すれば、第１コイル感度マップと第２コイル感度マップとは、複数のコイル各々において、強度画像は同一であるが位相画像における位相値が異なる。

なお、処理回路１３１は、取得機能３３により、第１コイル感度マップおよび第２コイル感度マップとは位相が異なり複数のコイルに関する奇数個のコイル感度マップをさらに取得してもよい。例えば、取得機能３３は、特定機能３５により特定されたコイル（以下、第３特定コイルと呼ぶ）の推定感度マップにおける位相画像を基準として、複数のコイルに関する推定感度マップの位相をシフトさせる。これにより、取得機能３３は、複数のコイルに亘る第３コイル感度マップを取得する。以下同様にして、取得機能３３は、特定機能３５により特定されたコイルを基準に応じて、複数のコイル感度マップを取得する。例えば、複数のコイルが８つの場合、取得機能３３は、８つのコイル各々を基準とした８つのコイル感度マップを取得してもよい。

処理回路１３１は、特定機能３５により、複数のコイルに対応する複数の感度マップ（推定感度マップ）の強度画像またはＭＲデータに基づいて再構成された強度画像に基づいて、位相の変化の基準となる１つのコイルを特定する。具体的には、特定機能３５は、複数のコイルに対応する複数の推定感度マップにおける複数の強度画像（以下、推定感度強度画像と呼ぶ）において、例えば、特異値分解により、信号値を多く有している順（以下、多信号値順と呼ぶ）に複数のコイルを特定する。多信号順は、特異値の高い順に相当する。特異値分解は、コイルコンプレッションなどの既知の技術が適用可能であるため、説明は省略する。なお、特異値分解の対象画像は、生成機能３７によりＭＲデータを再構成したコイルごとの強度画像（コントラスト画像）であってもよい。

具体的には、処理回路１３１は、特定機能３５により、特異値分解の結果に基づいて、複数の推定感度強度画像のうち、最も多く信号値を有している推定感度強度画像に対応するコイルを、第１特定コイルとして特定する。特定機能３５は、複数の推定感度強度画像のうち、２番目に多く信号値を有している推定感度強度画像に対応するコイルを、第２特定コイルとして特定する。特定機能３５は、複数の推定感度強度画像のうち、３番目に多く信号値を有している推定感度強度画像に対応するコイルを、第３特定コイルとして特定する。以下、同様にして、特定機能３５は、コイルを特定する。

処理回路１３１は、生成機能３７により、第１コイル感度マップと複数のコイルに関するＭＲデータとに基づいて第１画像を生成する。具体的には、生成機能３７は、複数のコイル各々に関するＭＲデータに基づいて、複数のコイル各々に対応する画像（以下、コイル画像と呼ぶ）を再構成する。ＭＲデータは被検体Ｐに対する間引き収集により取得されているため、コイル画像は、折り返しを有する画像となる。生成機能３７は、複数のコイルに対応する複数のコイル画像と、第１コイル感度マップとに基づく複素展開処理により、折り返しの無い１枚の第１画像を生成する。複素展開処理は、例えば、ＦＩＳＴＡ（Ｆａｓｔ ＩＳＴＡ（Ｉｔｅｒａｔｉｖｅ Ｓｈｒｉｎｋａｇｅ Ｓｏｆｔ－ｔｈｒｅｓｈｏｌｄｉｎｇ Ａｌｇｏｒｉｔｈｍ）などの既知の技術が適用可能であるため、説明は省略する。

処理回路１３１は、生成機能３７により、第２コイル感度マップと第１コイル感度マップと第１画像とに基づいて第２画像を生成する。例えば、生成機能３７は、第１画像と第１コイル感度マップとに基づいて、複数のコイルに応じて第１画像を分解し、分解画像を生成する。分解画像は、複数のコイル各々に対応し、複数のコイル各々の感度を反映した折り返しの無いＭＲ画像である。第１画像を分解した分解画像を生成する分解処理は、既知の技術が適用可能であるため、説明は省略する。生成機能３７は、第２コイル感度マップを用いて複数のコイルに対応する複数の分解画像を合成すること（以下、コイル合成と呼ぶ）により、第２画像を生成する。

処理回路１３１は、生成機能３７により、第１画像と第２画像とに基づいて、最終画像を生成する。例えば、生成機能３７は、第１画像に対して補間処理を実行して、第１補間画像を生成する。生成機能３７は、第２画像に対して補間処理を実行して、第２補間画像を生成する。補間処理は、例えば、複素補間などによる高解像度化を実現する処理（Ｚｏｏｍ）であって、例えば超解像処理に相当する。超解像処理は、隣接する画素における画素値の線形補間、またはｋ空間におけるゼロフィルなど、既知の技術が適用可能であるため、説明は省略する。

処理回路１３１は、生成機能３７により、第１補間画像と第２補間画像とに基づいて最終画像を生成する。具体的には、処理回路１３１は、生成機能３７により、第１補間画像と第２補間画像とに対して統計処理を実行して、最終画像を生成する。より詳細には、生成機能３７は、統計処理として、第１補間画像と第２補間画像とにおいて対応する位置の２つの画素における２つの画素値の平均を計算すること、または当該２つの画素値のうち最大値または中央値を選択することにより、最終画像を生成する。

また、取得機能３３により、第１コイル感度マップおよび第２コイル感度マップとは位相が異なり複数のコイルに関する奇数個のコイル感度マップがさらに取得された場合、処理回路１３１は、生成機能３７により、奇数個のコイル感度マップと第１コイル感度マップと第１画像とに基づいて奇数個の画像を生成する。次いで、生成機能３７は、第１画像と第２画像と奇数個の画像とにおいて対応する奇数個の画素における奇数個の画素値において中央値を選択することにより、第１画像と第２画像と奇数個の画像とに基づく最終画像を生成する。

以上のように構成された本実施形態のＭＲＩ装置１００により実行される画像生成処理について、図２乃至図１８を用いて説明する。画像生成処理は、急峻な位相の変化に伴って発生するアーチファクトを低減して、例えば、超解像画像を最終画像として生成する処理である。以下、説明を具体的にするために、受信コイル１１７における複数のコイルは８つであって、特定機能３５により特定されるコイルは、３つであるものとする。

図２は、画像生成処理の手順の一例を示すフローチャートである。図３は、画像生成処理における各種処理と当該各種処理に関するデータとの関連の一例を示す図である。

（画像生成処理）
（ステップＳ２０１）
撮像制御回路１２１は、被検体Ｐに対して間引き撮像を実行する（Ｐｒ１）。これにより、処理回路１３１は、取得機能３３により、複数のコイル各々に関するＭＲデータＤＡ１を取得する。

（ステップＳ２０２）
処理回路１３１は、取得機能３３により、間引き収集されたＭＲデータＤＡ１に基づいて、例えば、ＥＳＰＩＲｉＴにより、複数のコイル各々において推定感度マップを推定する（Ｐｒ２）。これにより、取得機能３３は、複数のコイルに対応する複数の推定感度マップＤＡ２を生成する。このとき、特定機能３５は、推定感度マップＤＡ２における推定感度強度画像に対する特異値分解により、多信号順に３つのコイル、すなわち第１特定コイルと第２特定コイルと第３特定コイルとを特定する。なお、特定機能３５は、コイル画像ＤＡ４に対する特異値分解により、多信号順に３つのコイルを特定してもよい。このとき、ステップＳ２０４は、ステップＳ２０３より前に実行されることとなる。

図４は、推定感度強度画像ＥＳＩの一例を示す図である。特定機能３５は、推定感度強度画像ＥＳＩ各々において特異値分解を行い、当該特異値分解の結果に基づいて、多信号順に、第１特定コイルと第２特定コイルと第３特定コイルとを特定する。

（ステップＳ２０３）
処理回路１３１は、取得機能３３により、複数の推定感度マップＤＡ２に基づいて、所定のコイルに対応する推定感度マップを基準とした第１コイル感度マップＤＡ３を生成する。具体的には、取得機能３３は、第１特定コイルに対応する推定感度マップの位相画像（以下、第１位相画像と呼ぶ）を基準として、複数のコイルにおける位相画像の位相をシフトさせる（Ｐｒ３）。これにより、取得機能３３は、第１コイル感度マップＤＡ３を生成する。

図５は、第１コイル感度マップＤＡ３における複数の位相画像ＰＩ１の一例を示す図である。図５に示すように、第１位相画像ＰＩ１１は、位相シフトにより、均一な値となる。換言すれば、位相シフトは、特定された第１特定コイルに対応する第１位相画像ＰＩ１１を均一にするように、推定感度マップにおける複数の位相画像における位相をシフトさせることに対応する。

（ステップＳ２０４）
処理回路１３１は、生成機能３７により、ＭＲデータＤＡ１に基づいて、複数のコイルに対応する複数のコイル画像ＤＡ４を再構成する（Ｐｒ４）。図６は、複数のコイル画像ＤＡ４における強度画像ＣＩＭの一例を示す図である。図６に示すように、複数のコイル画像ＤＡ４における強度画像は、間引き撮像により、折り返しアーチファクトを有する。図７は、複数のコイル画像ＤＡ４における位相画像ＣＩＰの一例を示す図である。図６および図７に示す複数のコイル画像ＤＡ４は、第１画像ＤＡ５の生成に用いられる。なお、本ステップにおける処理は、ステップＳ２０１の後の任意の段階で実行可能である。

（ステップＳ２０５）
処理回路１３１は、生成機能３７により、複数のコイル画像ＤＡ４と第１コイル感度マップＤＡ３とに基づいて、第１画像ＤＡ５を生成する。具体的には、生成機能３７は、複数のコイル画像ＤＡ４と第１コイル感度マップＤＡ３とを用いたＦＩＳＴＡなどの複素展開処理を実行する（Ｐｒ５）。これにより、生成機能３７は、第１画像を生成する（ＤＡ５）。図８は、第１画像ＤＡ５の一例を示す図である。図８に示すように、第１画像ＤＡ５において、折り返しは解消されている。

（ステップＳ２０６）
処理回路１３１は、生成機能３７により、第１コイル感度マップＤＡ３を用いて第１画像ＤＡ５を分解する（Ｐｒ６）。これにより、生成機能３７は、複数のコイルに対応する複数の分解画像ＤＡ６を生成する。図９は、複数のコイルに対応する複数の分解画像ＤＡ６における強度画像ＤＩＭの一例を示す図である。図９に示す強度画像ＤＩＭは、図６に示す強度画像ＣＩＭに比べて、折り返しアーチファクトが低減されていることにある。図１０は、複数のコイルに対応する複数の分解画像ＤＡ６における位相画像ＤＩＰの一例を示す図である。図９および図１０に示す複数の分解画像ＤＡ６は、第２画像の生成に用いられる。

（ステップＳ２０７）
処理回路１３１は、生成機能３７により、第１コイル感度マップＤＡ３に基づいて、第１コイル感度マップＤＡ３とは位相が異なる第２コイル感度マップＤＡ７を生成する。具体的には、生成機能３７は、第１コイル感度マップＤＡ３において、第２特定コイルに対応する位相画像（以下、第２位相画像と呼ぶ）ＰＩ２２を基準として、複数のコイルにおける位相画像の位相をシフトさせる（Ｐｒ７）。これにより、取得機能３３は、第２コイル感度マップＤＡ７を生成する。

図１１は、第２コイル感度マップＤＡ７における複数の位相画像ＰＩ２の一例を示す図である。図１１に示すように、図５に示す第２位相画像ＰＩ２２は、位相シフトにより、均一な値となる。換言すれば、位相シフトは、第１コイル感度マップＤＡ３における第２位相画像ＰＩ２２を均一にするように、第１コイル感度マップＤＡ３における複数の位相画像の位相をシフトさせることに対応する。

処理回路１３１は、生成機能３７により、第２コイル感度マップＤＡ７に基づいて、第１コイル感度マップＤＡ３および第２コイル感度マップＤＡ７とは位相が異なる第３コイル感度マップＤＡ８を生成する。具体的には、生成機能３７は、第２コイル感度マップＤＡ７において、第３特定コイルに対応する位相画像（以下、第３位相画像と呼ぶ）ＰＩ３３を基準として、複数のコイルにおける位相画像の位相をシフトさせる（Ｐｒ８）。これにより、取得機能３３は、第３コイル感度マップＤＡ８を生成する。

図１２は、第３コイル感度マップＤＡ８における複数の位相画像ＰＩ３の一例を示す図である。図１２に示すように、図１１に示す第３位相画像ＰＩ３３は、位相シフトにより、均一な値となる。換言すれば、位相シフトは、第２コイル感度マップＤＡ７における第３位相画像ＰＩ３３を均一にするように、第２コイル感度マップＤＡ７における複数の位相画像の位相をシフトさせることに対応する。

（ステップＳ２０８）
処理回路１３１は、生成機能３７により、分解画像ＤＡ６と第２コイル感度マップＤＡ７とに基づいて、第２画像ＤＡ９を生成する（Ｐｒ９）。具体的には、生成機能３７は、第２コイル感度マップＤＡ７を用いて、分解画像ＤＡ６に対してコイル合成を実行し、第２画像ＤＡ９を生成する。

図１３は、コイル合成された第２画像ＤＡ９の一例を示す図である。図１３に示す第２画像ＤＡ９と、図８に示す第１画像ＤＡ５との差異は、第２画像ＤＡ９における位相画像ＰＩ９と、第１画像ＤＡ５における位相画像ＰＩ５とが、位相シフトの影響により異なることにある。

処理回路１３１は、生成機能３７により、分解画像ＤＡ６と第３コイル感度マップＤＡ８とに基づいて、第３画像ＤＡ１０を生成する（Ｐｒ１０）。具体的には、生成機能３７は、第３コイル感度マップＤＡ８を用いて、分解画像ＤＡ６に対してコイル合成を実行し、第３画像ＤＡ１０を生成する。

図１４は、コイル合成された第３画像ＤＡ１０の一例を示す図である。図１４に示す第３画像ＤＡ１０と、図８に示す第１画像ＤＡ５と、図１３に示す第２画像ＤＡ９との差異は、第３画像ＤＡ１０における位相画像ＰＩ１０と、第２画像ＤＡ９における位相画像ＰＩ９と、第１画像ＤＡ５における位相画像ＰＩ５とが、位相シフトの影響により異なることにある。

（ステップＳ２０９）
処理回路１３１は、生成機能３７により、第１画像ＤＡ５に対して補間処理Ｐｒ１１を実行し、第１補間画像ＤＡ１１を生成する。第１補間画像の生成は、第１画像ＤＡ５の生成の後であれば、任意の段階で実施可能である。

図１５は、第１補間画像ＤＡ１１の一例を示す図である。図１５に示すように、第１補間画像ＤＡ１１における強度画像ＭＩ１１には、アーチファクトＡが発生している。図１５に示すように、アーチファクトＡの発生位置は、第１補間画像ＤＡ１１における位相画像ＰＩ１１において、急峻な位相の変化が発生している位置（以下、位相急峻変化位置と呼ぶ）に対応する。

処理回路１３１は、生成機能３７により、第２画像ＤＡ９に対して補間処理Ｐｒ１２を実行し、第２補間画像ＤＡ１２を生成する。第２補間画像の生成は、第２画像ＤＡ９の生成の後であれば、任意の段階で実施可能である。

図１６は、第２補間画像ＤＡ１２の一例を示す図である。図１６に示すように、第２補間画像ＤＡ１２における強度画像ＭＩ１２には、アーチファクトＡが発生している。図１６に示すように、アーチファクトＡの発生位置は、第２補間画像ＤＡ１２における位相画像ＰＩ１２において、急峻な位相の変化が発生している位置に対応する。

処理回路１３１は、生成機能３７により、第３画像ＤＡ１０に対して補間処理Ｐｒ１３を実行し、第３補間画像ＤＡ１３を生成する。第３補間画像の生成は、第３画像ＤＡ１０の生成の後であれば、任意の段階で実施可能である。

図１７は、第３補間画像ＤＡ１３の一例を示す図である。図１７に示すように、第３補間画像ＤＡ１３における強度画像ＭＩ１３には、アーチファクトＡが発生している。図１７に示すように、アーチファクトＡの発生位置は、第３補間画像ＤＡ１３における位相画像ＰＩ１３において、急峻な位相の変化が発生している位置に対応する。

図１５乃至図１７に示すように、位相シフトの影響により、アーチファクトＡの発生位置はそれぞれ異なる。また、図１５乃至図１７に示すように、アーチファクトＡの画素値は、アーチファクトＡの周辺の画素値に比べて黒く、すなわち小さい画素値となる。

（ステップＳ２１０）
処理回路１３１は、生成機能３７により、第１補間画像ＤＡ１１と第２補間画像ＤＡ１２と第３補間画像ＤＡ１３とに対して統計処理を実行する（Ｐｒ１４）。これにより、生成機能３７は、最終画像ＤＡ１４を生成する。具体的には、生成機能３７は、第１補間画像ＤＡ１１における強度画像ＭＩ１１と第２補間画像ＤＡ１２における強度画像ＭＩ１２と第３補間画像ＤＡ１３における強度画像ＭＩ１３とにおいて、同一の位置の画素における３つの画素値に対して最大値または中央値（Ｍｅｄｉａｎ）を選択する。生成機能３７は、当該画素値の選択を強度画像の全画素について実行し、選択された画素値を用いて最終画像ＤＡ１４を生成する。このような画素値の選択により、アーチファクトＡに対応する画素値は、非選択となる。

図１８は、最終画像ＤＡ１４の一例を示す図である。図１８に示すように、最終画像ＤＡ１４では、図１５に示す強度画像ＭＩ１１におけるアーチファクトＡと、図１６に示す強度画像ＭＩ１２におけるアーチファクトＡと、図１７に示す強度画像ＭＩ１３におけるアーチファクトＡとが除去されている。

システム制御回路１２３は、生成された最終画像ＤＡ１４をディスプレイ１２９に表示する。また、システム制御回路１２３は、生成された最終画像ＤＡ１４を、メモリ１２５に記憶する。以上により、画像生成処理は終了する。

以上に述べた実施形態に係るＭＲＩ装置１００は、複数のコイルに関する第１コイル感度マップＤＡ３と、第１コイル感度マップＤＡ３とは位相が異なり複数のコイルに関する第２コイル感度マップＤＡ７とを取得し、第１コイル感度マップＤＡ３と複数のコイルに関するＭＲデータＤＡ１とに基づいて第１画像ＤＡ５を生成し、第２コイル感度マップＤＡ７と第１コイル感度マップＤＡ３と第１画像ＤＡ５とに基づいて第２画像ＤＡ９を生成し、第１画像ＤＡ５と第２画像ＤＡ９とに基づいて、最終画像ＤＡ１４を生成する。

具体的には、本ＭＲＩ装置１００は、第１画像ＤＡ５に対して補間処理Ｐｒ１１を実行して、第１補間画像ＤＡ１１を生成し、第２画像ＤＡ９に対して補間処理Ｐｒ１２を実行して、第２補間画像ＤＡ１２を生成し、第１補間画像ＤＡ１１と第２補間画像ＤＡ１２とに基づいて最終画像ＤＡ１４を生成する。より詳細には、本ＭＲＩ装置１００は、第１補間画像ＤＡ１１と第２補間画像ＤＡ１２とに対して統計処理Ｐｒ１４を実行して、最終画像ＤＡ１４を生成する。例えば、本ＭＲＩ装置１００は、統計処理Ｐｒ１４として、第１補間画像ＤＡ１１と第２補間画像ＤＡ１２とにおいて対応する位置の２つの画素における２つの画素値の平均を計算すること、または２つの画素値のうち最大値または中央値を選択することにより、最終画像ＤＡ１４を生成する。

また、本ＭＲＩ装置１００は、第１コイル感度マップＤＡ３および第２コイル感度マップＤＡ７とは位相が異なり複数のコイルに関する奇数個のコイル感度マップ（例えば、第３コイル感度マップＤＡ８）をさらに取得し、奇数個のコイル感度マップと第１コイル感度マップＤＡ３と第１画像ＤＡ５とに基づいて奇数個の画像（例えば、第３画像ＤＡ１０）を生成し、第１画像ＤＡ５と第２画像ＤＡ９と奇数個の画像とに関して対応する奇数個の画素における奇数個の画素値において中央値を選択することにより、第１画像ＤＡ５と第２画像ＤＡ９と奇数個の画像とに基づく最終画像ＤＡ１４を生成する。

また、本ＭＲＩ装置１００は、間引き収集されたＭＲデータＤＡ１に基づいて、複数のコイルに関する推定感度マップＤＡ２を推定し、推定された推定感度マップＤＡ２に基づいて、第１コイル感度マップＤＡ３を取得する。

これらのことから、実施形態に係るＭＲＩ装置１００によれば、感度マップにおける位相をシフトさせることで、異なる位相シフトに応じたコイル合成後の複数の位相画像におおいて、アーチファクトＡの発生位置を移動させることができる。このため、本ＭＲＩ装置１００によれば、アーチファクトＡの発生位置となる位相急峻変化位置を避けて、最終画像ＤＡ１４を生成することができる。

以上のことから、本ＭＲＩ装置１００によれば、アーチファクトＡを低減させて、または、アーチファクトＡを効果的に除去して、例えば、高解像度化された最終画像ＤＡ１４を生成することができる。したがって、本ＭＲＩ装置１００によれば、最終画像ＤＡ１４における画質を改善することができ、被検体Ｐに対する検査効率および診断効率を向上させることができる。なお、本実施形態の変形例として、特定機能８５による処理を省略してもよい。この場合においても、実施形態と同様な効果を実現することができる。

また、本ＭＲＩ装置１００は、複数のコイルのうち１つのコイルの感度を基準として複数のコイルの感度の位相画像における位相を変化させて、第１コイル感度マップＤＡ３と第２コイル感度マップＤＡ７とを取得する。このとき、本ＭＲＩ装置１００は、複数のコイルに対応する複数の推定感度マップの強度画像またはＭＲデータに基づいて再構成されたコイル画像の強度画像に基づいて、位相の変化の基準となる１つのコイルを特定する。

これらのことから、実施形態に係るＭＲＩ装置１００によれば、全てのコイルに亘って位相シフトを行ってコイル感度マップを生成する必要はなく、位相シフトによる影響が大きい順に、例えば第１画像ＤＡ５、第２画像ＤＡ９、第３画像ＤＡ１０などを用いて、最終画像ＤＡ１４を生成することができる。これにより、本ＭＲＩ装置１００によれば、画像生成処理における計算コストを抑えて、位相急峻変化位置を効率的に移動させることができる。このため、本ＭＲＩ装置１００によれば、計算コストに関する各種費用を低減し、かつ計算時間を短縮することができる。したがって、本ＭＲＩ装置１００によれば、最終画像ＤＡ１４における画質を改善しつつ、被検体Ｐに対する検査効率および診断効率をさらに向上させることができる。

実施形態における技術的思想を画像生成方法で実現する場合、当該画像生成方法は、複数のコイルに関する第１コイル感度マップＤＡ３と、第１コイル感度マップＤＡ３とは位相が異なり複数のコイルに関する第２コイル感度マップＤＡ７とを取得し、第１コイル感度マップＤＡ３と複数のコイルに関するＭＲデータＤＡ１とに基づいて第１画像ＤＡ５を生成し、第２コイル感度マップＤＡ７と第１コイル感度マップＤＡ３と第１画像ＤＡ５とに基づいて第２画像ＤＡ９を生成し、第１画像ＤＡ５と第２画像ＤＡ９とに基づいて、最終画像ＤＡ１４を生成する。本画像生成方法に関する画生成処理の手順および効果は、実施形態における記載と同様なため、説明は省略する。

実施形態における技術的思想を画像生成プログラムで実現する場合、当該画像生成プログラムは、コンピュータに、複数のコイルに関する第１コイル感度マップＤＡ３と、第１コイル感度マップＤＡ３とは位相が異なり複数のコイルに関する第２コイル感度マップＤＡ７とを取得し、第１コイル感度マップＤＡ３と複数のコイルに関するＭＲデータＤＡ１とに基づいて第１画像ＤＡ５を生成し、第２コイル感度マップＤＡ７と第１コイル感度マップＤＡ３と第１画像ＤＡ５とに基づいて第２画像ＤＡ９を生成し、第１画像ＤＡ５と第２画像ＤＡ９とに基づいて、最終画像ＤＡ１４を生成すること、を実現させる。

例えば、ＭＲＩ装置１００などのモダリティやＰＡＣＳサーバなどにおけるコンピュータに画像生成プログラムをインストールし、これらをメモリ上で展開することによっても、画像生成処理を実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（ハードディスクなど）、光ディスク（ＣＤ－ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記憶媒体に格納して頒布することも可能である。画像生成プログラムによる画像生成処理の手順および効果は、実施形態と同様なため、説明は省略する。

以上説明した少なくとも１つの実施形態等によれば、急峻な位相の変化に伴って発生するアーチファクトを低減することができる。

いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

３３ 取得機能
３５ 特定機能
３７ 生成機能
１００ 磁気共鳴イメージング装置
１０１ 静磁場磁石
１０３ 傾斜磁場コイル
１０５ 傾斜磁場電源
１０７ 寝台
１０９ 寝台制御回路
１１１ ボア
１１３ 送信回路
１１５ 送信コイル
１１７ 受信コイル
１１９ 受信回路
１２１ 撮像制御回路
１２３ システム制御回路
１２５ メモリ
１２７ 入力インターフェース
１２９ ディスプレイ
１３１ 処理回路

Claims (10)

1. 複数のコイルに関する第１コイル感度マップと、前記第１コイル感度マップとは位相が異なり前記複数のコイルに関する第２コイル感度マップとを取得する取得部と、
   前記第１コイル感度マップと前記複数のコイルに関する磁気共鳴データとに基づいて第１画像を生成し、前記第２コイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて第２画像を生成し、前記第１画像と前記第２画像とに基づいて、最終画像を生成する生成部と、
   を備える磁気共鳴イメージング装置。
2. 前記生成部は、
   前記第１画像に対して補間処理を実行して、第１補間画像を生成し、
   前記第２画像に対して補間処理を実行して、第２補間画像を生成し、
   前記第１補間画像と前記第２補間画像とに基づいて前記最終画像を生成する、
   請求項１に記載の磁気共鳴イメージング装置。
3. 前記生成部は、前記第１補間画像と前記第２補間画像とに対して統計処理を実行して、前記最終画像を生成する、
   請求項２に記載の磁気共鳴イメージング装置。
4. 前記生成部は、前記統計処理として、前記第１補間画像と前記第２補間画像とにおいて対応する位置の２つの画素における２つの画素値の平均を計算すること、または前記２つの画素値のうち最大値または中央値を選択することにより、前記最終画像を生成する、
   請求項３に記載の磁気共鳴イメージング装置。
5. 前記取得部は、前記第１コイル感度マップおよび前記第２コイル感度マップとは位相が異なり前記複数のコイルに関する奇数個のコイル感度マップをさらに取得し、
   前記生成部は、
   前記奇数個のコイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて奇数個の画像を生成し、
   前記第１画像と前記第２画像と前記奇数個の画像とに関して対応する奇数個の画素における奇数個の画素値において中央値を選択することにより、前記第１画像と前記第２画像と前記奇数個の画像とに基づく前記最終画像を生成する、
   請求項１に記載の磁気共鳴イメージング装置。
6. 前記取得部は、間引き収集された前記磁気共鳴データに基づいて、複数のコイルに関する感度マップを推定し、前記推定された感度マップに基づいて、前記第１コイル感度マップを取得する、
   請求項１乃至５のうちいずれか一項に記載の磁気共鳴イメージング装置。
7. 前記取得部は、前記複数のコイルのうち１つのコイルの感度を基準として前記複数のコイルの感度の位相画像における前記位相を変化させて、前記第１コイル感度マップと前記第２コイル感度マップとを取得する、
   請求項１乃至６のうちいずれか一項に記載の磁気共鳴イメージング装置。
8. 前記複数のコイルに対応する複数の感度マップの強度画像または前記磁気共鳴データに基づいて再構成された強度画像に基づいて、前記位相の変化の基準となる前記１つのコイルを特定する特定部をさらに備える、
   請求項７に記載の磁気共鳴イメージング装置。
9. 複数のコイルに関する第１コイル感度マップと、前記第１コイル感度マップとは位相が異なり前記複数のコイルに関する第２コイル感度マップとを取得し、
   前記第１コイル感度マップと前記複数のコイルに関する磁気共鳴データとに基づいて第１画像を生成し、
   前記第２コイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて第２画像を生成し、
   前記第１画像と前記第２画像とに基づいて、最終画像を生成すること、
   を備える画像生成方法。
10. コンピュータに、
    複数のコイルに関する第１コイル感度マップと、前記第１コイル感度マップとは位相が異なり前記複数のコイルに関する第２コイル感度マップとを取得し、
    前記第１コイル感度マップと前記複数のコイルに関する磁気共鳴データとに基づいて第１画像を生成し、
    前記第２コイル感度マップと前記第１コイル感度マップと前記第１画像とに基づいて第２画像を生成し、
    前記第１画像と前記第２画像とに基づいて、最終画像を生成すること、
    を実現させる画像生成プログラム。

Images