(57)【要約】      （修正有）【課題】拡散情報を含むデータ収集において、できるだけコントラストの変化を生じさせない磁気共鳴装置を提供する。【解決手段】シーケンスＳＣ１およびＳＣ２により収集されたデータに基づいて、複数の拡散強調画像を作成する第１の画像作成手段９１と、複数の拡散強調画像に基づいてマスクを作成するマスク作成手段９２と、拡散強調画像に基づいてモジュレーションマップを作成するモジュレーションマップ作成手段９３と、マスクとモジュレーションマップとを用いて拡散強調画像の準同型フィルタリングを実行する準同型フィルタ手段９４とを有するＭＲ装置。【選択図】図２
2015208420.tif 000002
【特許請求の範囲】
【請求項１】
複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行するスキャン手段と、
前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成手段と、
前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための第１のフィルタリングを実行する第１のフィルタ手段であって、前記第１の画像の中の第２の部分において前記第１のフィルタリングを実行しない第１のフィルタ手段と、
を有する磁気共鳴装置。
【請求項２】
前記第１の画像の前記第１の部分に対応する第１の領域と、前記第１の画像の前記第２の部分に対応する第２の領域とを含むマスクを作成するマスク作成手段を有し、
前記フィルタ手段は、
前記第１の画像の中で、前記マスクの第１の領域に対応する部分において、前記第１のフィルタリングを実行するが、前記第１の画像の中で、前記マスクの第１の領域に対応する部分において、前記第１のフィルタリングを実行しない、請求項１に記載の磁気共鳴装置。
【請求項３】
前記マスク作成手段は、
前記複数の第１の画像のピクセル値の標準偏差を表す標準偏差マップを作成し、前記標準偏差マップに基づいて前記マスクを作成する、請求項２に記載の磁気共鳴装置。
【請求項４】
前記第１のフィルタリングは準同型フィルタリングである、請求項２又は３に記載の磁気共鳴装置。
【請求項５】
前記第１の画像を準同型フィルタリングするためのモジュレーションマップを作成するモジュレーションマップ作成手段を有し、
前記フィルタ手段は、
前記モジュレーションマップを用いて、前記第１の画像の準同型フィルタリングを実行する、請求項４に記載の磁気共鳴装置。
【請求項６】
前記モジュレーションマップ作成手段は、
前記第１の画像に基づいて前記モジュレーションマップを作成する、請求項５に記載の磁気共鳴装置。
【請求項７】
前記第１の画像を準同型フィルタリングすることにより得られた第２の画像と、前記第１の画像とを差分することにより、差分画像を作成する差分画像作成手段と、
前記差分画像をローパスフィルタリングするための第２のフィルタ手段と、
ローパスフィルタリングされた差分画像と前記第１の画像とを加算することによりに得られる第１の加算画像を作成する第２の画像作成手段と、
を有する、請求項４〜６のうちのいずれか一項に記載の磁気共鳴装置。
【請求項８】
前記第２のフィルタ手段は、
前記マスクの前記第１の領域と前記第２の領域との境界におけるピクセル値が滑らかに変化するように、前記マスクをローパスフィルタリングし、
前記第２の画像作成手段は、
ローパスフィルタリングされたマスクを用いて、前記第１の加算画像から第１の画像部分を取り出し、
前記ローパスフィルタリングされたマスクの反転マスクを用いて、前記第１の画像から第２の画像部分を取り出し、
前記第１の画像部分と前記第２の画像部分とを加算することにより、第３の画像を作成する、請求項７に記載の磁気共鳴装置。
【請求項９】
前記第２の画像作成手段は複数の第３の画像を作成し、
前記複数の第３の画像を合成することにより合成画像を作成する第３の画像作成手段を有する、請求項８に記載の磁気共鳴装置。
【請求項１０】
前記第３の画像作成手段は、
前記複数の第３の画像の各々を重み付けし、重み付けされた複数の第３の画像に基づいて前記合成画像を作成する、請求項９に記載の磁気共鳴装置。
【請求項１１】
前記第３の画像作成手段は、
振幅二乗重み付き平均法を用いて前記複数の第２の画像の各々を重み付けする、請求項１０に記載の磁気共鳴装置。
【請求項１２】
前記合成画像と、前記合成画像とは別の合成画像とに基づいて、ＡＤＣマップを作成するＡＤＣマップ作成手段を有する、請求項９〜１１のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１３】
前記スキャン手段は、
前記スライスからデータを収集するための第２のシーケンスを複数回実行し、
前記別の合成画像は、
前記第２のシーケンスを複数回実行することにより得られたデータに基づいて作成されている、請求項１２に記載の磁気共鳴装置。
【請求項１４】
前記第１のシーケンスは、拡散強調を行うための第１の勾配パルスであって、拡散強調を強さを表すｂ値が第１の値に設定された第１の勾配パルスを含み、
前記第２のシーケンスは、拡散強調を行うための第２の勾配パルスであって、拡散強調を強さを表すｂ値が第２の値に設定された第２の勾配パルスを含む、請求項１３に記載の磁気共鳴装置。
【請求項１５】
前記第１のシーケンスは、拡散強調を行うための第１の勾配パルスであって、拡散強調を強さを表すｂ値が第１の値に設定された第１の勾配パルスを含み、
前記第２のシーケンスは、拡散強調を行うための勾配パルスを含んでいない、請求項１３に記載の磁気共鳴装置。
【請求項１６】
前記第１の画像の前記第１の部分は肝臓の左葉および胃を含んでおり、前記第１の画像の前記第２の部分は肝臓の右葉を含んでいる、請求項１〜１５のうちのいずれか一項に記載の磁気共鳴装置。
【請求項１７】
複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成処理と、
前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための第１のフィルタリングを実行する第１のフィルタ処理であって、前記第１の画像の中の第２の部分において前記第１のフィルタリングを実行しない第１のフィルタ処理と、
をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、拡散情報を含むデータを収集する磁気共鳴装置、およびこの磁気共鳴装置に適用されるプログラムに関する。
【背景技術】
【０００２】
近年、拡散強調イメージング（ＤＷＩ：Diffusion Weighted Imaging）が普及している。ＤＷＩは、ＭＰＧ（Motion Probing Gradient）と呼ばれる磁場強度の大きい勾配磁場を印加することによって、撮影部位の拡散情報を含む拡散強調画像を取得する方法である。拡散強調画像を用いることにより、見かけの拡散係数（ＡＤＣ：Apparent Diffusion Coefficient）など、種々の拡散情報を取得することができる。しかし、ＭＰＧの印加中に、拍動、蠕動により組織（tissue）が動いてしまうと、画像強度の損失（image intensity loss）が生じることが知られている。このような損失が生じると、見かけの拡散係数（以下、「ＡＤＣ」と呼ぶ）の値の推定誤差が大きくなるという問題がある。そこで、この問題に対処するために、モジュレーションマップ（Modulation Map）を用いて準同型フィルタリング（Homomorphic Filtering）を実行する方法が知られている（非特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１５７６８７号公報
【非特許文献】
【０００４】
【非特許文献１】Lim, Two-dimensional signal and image processing, Prentice Hall, 1990
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、準同型フィルタリングを実行した場合、画像のコントラストが変化するという問題がある。
したがって、できるだけコントラストの変化を生じさせない技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行するスキャン手段と、
前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成手段と、
前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための第１のフィルタリングを実行する第１のフィルタ手段であって、前記第１の画像の中の第２の部分において前記第１のフィルタリングを実行しない第１のフィルタ手段とを有する磁気共鳴装置である。
【０００７】
本発明の第２の観点は、複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行する磁気共鳴装置に適用されるプログラムであって、
前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成処理と、
前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための第１のフィルタリングを実行する第１のフィルタ処理であって、前記第１の画像の中の第２の部分において前記第１のフィルタリングを実行しない第１のフィルタ処理と、
をコンピュータに実行させるためのプログラムである。
【発明の効果】
【０００８】
第１の画像の中の第２の部分において第１のフィルタリングが実行されないので、コントラストの変化を防止することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の一形態の磁気共鳴装置の概略図である。
【図２】プロセッサ９が実行する処理を示す図である。
【図３】ＭＲ装置の動作フローを示す図である。
【図４】スキャンを実行するときに設定されるスライスを示す図である。
【図５】スライスの拡散強調画像を取得するためのスキャンＳＣの説明図である。
【図６】スキャンＳＣ１の説明図である。
【図７】シーケンスＥ１の一例を示す図である。
【図８】期間Ｖ１のシーケンスＥ１〜Ｅ３により得られた拡散強調画像ＤＷ１１、ＤＷ２１、ＤＷ３１（ｂ＝ｂ１）を概略的に示す図である。
【図９】期間Ｖ２のシーケンスＥ１〜Ｅ３により得られた拡散強調画像を概略的に示す図である。
【図１０】期間ＶｎにおいてシーケンスＥ１、Ｅ２、およびＥ３を実行することにより取得された拡散強調画像を概略的に示す図である。
【図１１】スキャンＳＣ２の説明図である。
【図１２】シーケンスＦ１の一例を示す図である。
【図１３】期間Ｗ１〜Ｗｎにおいて取得された拡散強調画像を概略的に示す図である。
【図１４】拡散強調画像ＤＷ１１の一例を示す図である。
【図１５】拡散強調画像ＤＷ１１〜ＤＷ１ｎを示す図である。
【図１６】ＳＤマップの作成方法を概略的に示す図である。
【図１７】スライス２のＳＤマップ２およびスライス３のＳＤマップ３を概略的に示す図である。
【図１８】ｂ＝ｂ２において、スライスごとに取得されたＳＤマップを概略的に示す図である。
【図１９】ＳＤマップの一例を示す図である。
【図２０】マスクＭＫ１〜ＭＫ３を作成する方法の説明図である。
【図２１】マスクＭＫ４〜ＭＫ６を作成する方法の説明図である。
【図２２】マスクＭＫ１の具体的な一例を示す図である。
【図２３】モジュレーションマップの作成方法の説明図である。
【図２４】準同型フィルタリングが実行される領域を限定せずに式（１）を用いて準同型フィルタリングを実行する例の説明図である。
【図２５】準同型フィルタリングが実行される領域を限定し、式（１）を用いて準同型フィルタリングを実行する例の説明図である。
【図２６】ｂ＝ｂ１における準同型フィルタ画像を示す図である。
【図２７】ｂ＝ｂ２における準同型フィルタ画像を示す図である。
【図２８】差分画像ＤＩＦの一例を示す図である。
【図２９】マスクＭＫ１をローパルフィルタリングするときの説明図である。
【図３０】差分画像ＤＩＦをローパスフィルタリングするときの説明図である。
【図３１】第２項の説明図である。
【図３２】第１項の説明図である。
【図３３】画像部分Ｄ１およびＤ２を加算するときの説明図である。
【図３４】ｂ＝ｂ１における復元画像を示す図である。
【図３５】ｂ＝ｂ２における復元画像を示す図である。
【図３６】復元画像の合成方法の一例を示す図である。
【図３７】合成画像ＤＣ１の一例を示す図である。
【図３８】ｂ＝ｂ２におけるスライス１の合成画像ＤＣ４を概略的に示す図である。
【図３９】シミュレーション条件の説明図である。
【図４０】シミュレーション結果を示す図である。
【図４１】ＡＤＣマップを示す図である。
【発明を実施するための形態】
【００１０】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１１】
図１は、本発明の一形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ）１００は、マグネット２、テーブル３、受信ＲＦコイル４などを有している。
【００１２】
マグネット２は、被検体１３が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、ＲＦコイルなど（図示せず）を有している。超伝導コイルは静磁場を印加し、勾配コイルは勾配磁場を印加し、ＲＦコイルはＲＦパルスを印加する。
【００１３】
テーブル３は、被検体１３を支持するためのクレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１３はボア２１に搬送される。
【００１４】
受信ＲＦコイル（以下、単に「コイル」と呼ぶ）４は、被検体１３の胴部に取り付けられている。
【００１５】
ＭＲ装置１００は、更に、送信器５、勾配磁場電源６、受信器７、コンピュータ８、操作部１１、および表示部１２などを有している。
【００１６】
送信器５はＲＦコイルに電流を供給し、勾配磁場電源６は勾配コイルに電流を供給する。受信器７は、コイル４が受信した信号に対して、検波などの信号処理を行う。尚、マグネット２、送信器５、勾配磁場電源６、受信器７を合わせたものがスキャン手段に相当する。
【００１７】
コンピュータ８はＭＲ装置１００の各部の動作を制御する。コンピュータ８はプロセッサ９およびメモリ１０などを有している。
【００１８】
メモリ１０には、プロセッサ９により実行されるプログラムが記憶されている。プロセッサ９は、メモリ１０に記憶されているプログラムを読み出し、プログラムに記述されている処理を実行する。図２に、プロセッサ９が実行する処理を示す。
【００１９】
画像作成手段９１は、スキャンにより収集されたデータに基づいて画像を作成する。
マスク作成手段９２は、準同型フィルタリングを実行するときに使用するマスク（図２２参照）を作成する。マスクの詳細については後述する。
モジュレーションマップ作成手段９３は、準同型フィルタリングを実行するときに使用するモジュレーションマップ（図２３参照）を作成する。モジュレーションマップの詳細については後述する。
準同型フィルタ手段９４は、マスクとモジュレーションマップを用いて、画像作成手段９１により作成された画像に対して準同型フィルタリングを実行する。準同型フィルタ手段９４は第１のフィルタ手段に相当する。
【００２０】
差分画像作成手段９５は、後述する差分画像（図２８参照）を作成する。
ローパスフィルタ手段９６は、マスクや差分画像に対してローパルフィルタリングを実行する。ローパルフィルタ手段９６は第２のフィルタ手段に相当する。
復元画像作成手段９７は、後述する式（３）に基づいて復元画像を求める。復元画像作成手段９７は第２の画像作成手段に相当する。
合成画像作成手段９８は、復元画像を合成することにより合成画像を作成する。
ＡＤＣマップ作成手段９９は、合成画像に基づいてＡＤＣマップを作成する。
【００２１】
プロセッサ９は、画像作成手段９１〜ＡＤＣマップ作成手段９９を構成する一例であり、メモリ１０に記憶されているプログラムを実行することにより、これらの手段として機能する。
【００２２】
操作部１１は、オペレータにより操作され、種々の情報をコンピュータ８に入力する。表示部１２は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２３】
本形態では、ＭＲ装置１００を用いて被検体をスキャンし、スキャンにより得られたデータに基づいてＡＤＣマップを作成する。以下に、ＭＲ装置１００の動作フローについて説明する。
【００２４】
図３は、ＭＲ装置の動作フローを示す図である。
ステップＳＴ１では、拡散強調画像を取得するためのスキャンが実行される（図４および図５参照）。
【００２５】
図４はスキャンを実行するときに設定されるスライスを示す図、図５はスライスの拡散強調画像を取得するためのスキャンＳＣの説明図である。
本形態では、肝臓を横切る位置にスライスが設定される。図４では、説明の便宜上、３枚のスライス１〜３が設定された例が示されている。
【００２６】
スキャンＳＣは、２つのスキャンＳＣ１およびＳＣ２に分けられる。スキャンＳＣ１は、拡散強調の強さを表すｂ値がｂ＝ｂ１に設定されたＭＰＧ（Motion Probing Gradient）を有するシーケンスによりスライス１〜３の拡散情報を取得するためのスキャンである。一方、スキャンＳＣ２は、ｂ＝ｂ２に設定されたＭＰＧを有するシーケンスによりスライス１〜３の拡散情報を取得するためのスキャンである。
【００２７】
図６はスキャンＳＣ１の説明図である。
図６には、各スライスのＮＥＸ（Number of Excitation：励起回数）がｎ回の例（ＮＥＸ＝ｎ）が示されている。スキャンＳＣ１は、ｎ個の期間Ｖ１〜Ｖｎを有している。各期間では、３つのシーケンスＥ１、Ｅ２、およびＥ３が実行される。したがって、スキャンＳＣ１では、３つのシーケンスＥ１、Ｅ２、およびＥ３のセットがｎ回実行される。以下に、スキャンＳＣ１の期間Ｖ１〜Ｖｎで実行されるシーケンスＥ１、Ｅ２、Ｅ３について説明する。
【００２８】
スキャンＳＣ１では、先ず、期間Ｖ１において、シーケンスＥ１が実行される（図７参照）。
【００２９】
図７は、シーケンスＥ１の一例を示す図である。
シーケンスＥ１は、以下のパルス（１）〜（４）を有している。
（１）ＲＦパルス
（１ａ）スライスを励起するための励起パルスα１
（１ｂ）スピンを再収束させるためのリフォーカスパルスα２
（２）スライス選択方向の勾配パルス
（２ａ）スライスを選択するためのスライス選択勾配パルスＧｓ１およびＧｓ２
（３）周波数エンコード方向の勾配パルス
（３ａ）勾配パルスＧｆｓ
（３ｂ）ｂ＝ｂ１のＭＰＧ
（３ｃ）エコーを収集するための読出し勾配パルスＧｒｅ
（４）位相エンコード方向の勾配パルス
（４ａ）勾配パルスＧｐｓ
（４ｂ）位相エンコード勾配パルスＧｐ
【００３０】
シーケンスＥ１は、スライス１を励起し、スライス１の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。本形態では、シーケンスＥ１は、スライス１の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間の全データを１回のショットで収集する（シングルショット法）。画像作成手段９１（図２参照）は、シーケンスＥ１により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＥ１が終了した後、シーケンスＥ２が実行される。
【００３１】
シーケンスＥ２も、シーケンスＥ１と同様に、図７に示すシーケンスチャートで表される。ただし、シーケンスＥ２は、スライス２が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＥ２は、スライス２を励起し、スライス２の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＥ２により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＥ２が終了した後、シーケンスＥ３が実行される。
【００３２】
シーケンスＥ３も、シーケンスＥ１と同様に、図７に示すシーケンスチャートで表される。ただし、シーケンスＥ３は、スライス３が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＥ３は、スライス３を励起し、スライス３の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＥ３により収集されたデータに基づいて、拡散強調画像を作成する。
【００３３】
したがって、期間Ｖ１においてシーケンスＥ１、Ｅ２、およびＥ３を実行することにより、スライス１〜３の各々の拡散強調画像を取得することができる。図８に、期間Ｖ１のシーケンスＥ１〜Ｅ３により得られた拡散強調画像ＤＷ１１、ＤＷ２１、ＤＷ３１（ｂ＝ｂ１）を概略的に示す。
期間Ｖ１の後、次の期間Ｖ２においてシーケンスが実行される。
【００３４】
期間Ｖ２でも、期間Ｖ１と同様に、シーケンスＥ１、Ｅ２、およびＥ３が実行される。したがって、期間Ｖ２でも、スライス１〜３の各々の拡散強調画像（ｂ＝ｂ１）を取得することができる。図９に、期間Ｖ２のシーケンスＥ１〜Ｅ３により得られた拡散強調画像を概略的に示す。図９では、期間Ｖ１で取得された拡散強調画像と、期間Ｖ２で取得された拡散強調画像とを区別するために、期間Ｖ２で取得された拡散強調画像を符号「ＤＷ１２」、「ＤＷ２２」、「ＤＷ３２」で示してある。
【００３５】
以下同様に、３つのシーケンスＥ１、Ｅ２、およびＥ３のセットが繰り返し実行される。そして、スキャンＳＣ１の最後の期間Ｖｎにおいて、シーケンスＥ１、Ｅ２、およびＥ３が実行される。図１０に、期間ＶｎにおいてシーケンスＥ１、Ｅ２、およびＥ３を実行することにより取得された拡散強調画像を概略的に示す。図１０では、期間Ｖｎで取得された拡散強調画像を符号「ＤＷ１ｎ」、「ＤＷ２ｎ」、「ＤＷ３ｎ」で示してある。
次に、スキャンＳＣ２について説明する。
【００３６】
図１１はスキャンＳＣ２の説明図である。
図１１は、図６と同様に、ＮＥＸ＝ｎの例を示している。スキャンＳＣ２は、ｎ個の期間Ｗ１〜Ｗｎを有している。各期間では、３つのシーケンスＦ１、Ｆ２、およびＦ３が実行される。したがって、スキャンＳ２では、３つのシーケンスＦ１、Ｆ２、およびＦ３のセットがｎ回実行される。以下に、スキャンＳＣ２の期間Ｗ１〜Ｗｎで実行されるシーケンスＦ１、Ｆ２、Ｆ３について説明する。
スキャンＳＣ２では、先ず、期間Ｗ１において、シーケンスＦ１が実行される（図１２参照）。
【００３７】
図１２は、シーケンスＦ１の一例を示す図である。
シーケンスＦ１は、シーケンスＥ１（図７参照）と比較すると、ｂ＝ｂ２に設定されたＭＰＧを有しているが、その他のパルスは、シーケンスＥ１と同じである。したがって、シーケンスＦ１は、スライス１を励起し、スライス１の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ１により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＦ１が終了した後、シーケンスＦ２が実行される。
【００３８】
シーケンスＦ２も、シーケンスＦ１と同様に、図１２に示すシーケンスチャートで表される。ただし、シーケンスＦ２は、スライス２が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＦ２は、スライス２を励起し、スライス２の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ２により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＦ２が終了した後、シーケンスＦ３が実行される。
【００３９】
シーケンスＦ３も、シーケンスＦ１と同様に、図１２に示すシーケンスチャートで表される。ただし、シーケンスＦ３は、スライス３が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＦ３は、スライス３を励起し、スライス３の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ３により収集されたデータに基づいて、拡散強調画像を作成する。
【００４０】
以下同様に、期間Ｗ２〜Ｗｎについても、３つのシーケンスＦ１、Ｆ２、およびＦ３のセットが繰り返し実行される。したがって、スキャンＳＣ２を実行することにより、１枚のスライスにつき、ｎ枚の拡散強調画像を取得することができる。図１３に、期間Ｗ１〜Ｗｎにおいて取得された拡散強調画像を概略的に示す。
【００４１】
このようにして、スキャンＳＣ１およびＳＣ２が実行される。スキャンＳＣ１およびＳＣ２は被検体が息止めしている間に実行してもよいし、呼吸同期法を用いて実行してもよい。
【００４２】
ただし、ＭＰＧの印加中に、拍動、蠕動などが原因で器官が動いてしまうと、器官の動きの影響を受けて、拡散強調画像に画像強度の損失（信号損失）が現れることが知られている。図１４に、拡散強調画像ＤＷ１１の一例を示す。図１４では、拡散強調画像ＤＷ１１の中の一部（例えば、肝臓の左葉、胃、および脾臓）に、画像強度の損失ＬＯ１、ＬＯ２、およびＬＯ３が現れている。このような画像強度の損失は、後述するステップＳＴ８においてＡＤＣマップを作成する場合、ＡＤＣ値の誤差の原因になるので、拡散強調画像ＤＷ１１の中で画像強度の損失が生じている部分の画像強度はできるだけ復元する必要がある。そこで、本形態では、スキャンＳＣ１およびＳＣ２を実行した後、画像強度を復元するためにステップＳＴ２を実行する。以下に、ステップＳＴ２について説明する。尚、ステップＳＴ２は、ステップＳＴ２１およびＳＴ２２を有しているので、ステップＳＴ２１およびＳＴ２２について順に説明する。
【００４３】
ステップＳＴ２１では、マスク作成手段９２（図２参照）が、後述するステップＳＴ２２で使用するマスク（図２２参照）を作成する。以下に、マスクの作成方法について説明する。
【００４４】
先ず、マスク作成手段９２は、ｂ＝ｂ１のスキャンＳＣ１により得られた拡散強調画像の中から、スライス１の拡散強調画像ＤＷ１１〜ＤＷ１ｎを選択する。図１５に、選択された拡散強調画像ＤＷ１１〜ＤＷ１ｎを示す。
【００４５】
次に、マスク作成手段９２は、拡散強調画像ＤＷ１１〜ＤＷ１ｎのピクセル値の標準偏差を表すＳＤマップを作成する。図１６に、ＳＤマップの作成方法を概略的に示す。
【００４６】
マスク作成手段９２は、拡散強調画像ＤＷ１１〜ＤＷ１ｎから同じ位置のピクセルのピクセル値を選択し、選択したピクセル値の標準偏差を求める。例えば、拡散強調画像ＤＷ１１〜ＤＷ１ｎのピクセルＰａのピクセル値の標準偏差を求める場合、マスク作成手段９２は、拡散強調画像ＤＷ１１〜ＤＷ１ｎからピクセルＰａのピクセル値Ｘａ１〜Ｘａｎを選択する。そして、ピクセル値Ｘａ１〜Ｘａｎの標準偏差σａを計算する。このようにして計算された標準偏差σａが、ＳＤマップ１のピクセルＰａにおける値として用いられる。
【００４７】
上記の説明では、ＳＤマップのピクセルＰａにおける標準偏差σａを求める方法について説明したが、ＳＤマップの他のピクセルにおける標準偏差も同じ方法で計算される。例えば、拡散強調画像ＤＷ１１〜ＤＷ１ｎのピクセルＰｂにおけるピクセル値Ｘｂ１〜Ｘｂｎの標準偏差σｂを計算することにより、ＳＤマップ１のピクセルＰｂにおける値を得ることができる。同様に、拡散強調画像ＤＷ１１〜ＤＷ１ｎのピクセルＰｃにおけるピクセル値Ｘｃ１〜Ｘｃｎの標準偏差σｃを計算することにより、ＳＤマップ１のピクセルＰｃにおける値を得ることができる。
【００４８】
したがって、ピクセルごとに標準偏差を計算することにより、ｂ＝ｂ１におけるスライス１のＳＤマップ１を得ることができる。
【００４９】
図１５および図１６では、スライス１のＳＤマップ１を作成する手順について示されているが、他のスライス２およびスライス３についても同様の方法でＳＤマップを得ることができる。図１７に、スライス２のＳＤマップ２およびスライス３のＳＤマップ３を概略的に示す。このように、スライスごとにＳＤマップを得ることができる。
【００５０】
尚、図１５〜図１７では、ｂ＝ｂ１のスキャンＳＣ１により取得された拡散強調画像を用いてＳＤマップを作成した例について示されているが、ｂ＝ｂ２のスキャンＳＣ２により取得された拡散強調画像に対しても、同様の方法でスライスごとにＳＤマップが作成される。図１８に、ｂ＝ｂ２において、スライスごとに取得されたＳＤマップを概略的に示す。
したがって、全部で６枚のＳＤマップが作成される。図１９に、ＳＤマップの一例を示す。
【００５１】
肝臓の左葉、胃、脾臓など、心臓の近くに位置している部位は、拍動の影響を受けて動いてしまうので、肝臓の左葉、胃、脾臓などの部位のピクセル値は時間とともに変化する傾向がある。一方、肝臓の右葉は、心臓から遠いので、心臓の拍動の影響をあまり受けない。したがって、肝臓の右葉のピクセル値の時間変化はそれほど大きくならない傾向がある。このため、ＳＤマップから、心臓の拍動の影響を受けやすい部位と、心臓の拍動の影響をあまり受けない部位とをおおまかに区別することができる。
【００５２】
次に、マスク作成手段９２は、ＳＤマップに基づいてマスクを作成する。図２０および図２１に、ＳＤマップに基づいてマスクを作成する方法を概略的に示す。マスク作成手段９２は、ＳＤマップ１〜６の各々に対してしきい値処理ＴＨを行う。しきい値ＴＨは、ＳＤマップのピクセル値のうち、しきい値より小さいピクセル値を「０」に設定し、一方、しきい値以上のピクセル値を「１」に設定する。したがって、各ＳＤマップは、「０」の領域と、「１」の領域とに分けられる。このように、各ＳＤマップを「０」の領域と「１」の領域とに分けることにより、マスクＭＫ１〜ＭＫ６を作成することができる。図２２に、マスクＭＫ１〜ＭＫ６のうち、代表して、マスクＭＫ１の具体的な一例を示す。図２２では、ピクセル値１の領域が白で示されており、ピクセル値０の領域が黒で示されている。
【００５３】
上記のように、ＳＤマップは、心臓の拍動などの影響を受けやすい部位のピクセル値と、心臓の拍動などの影響をあまり受けない部位のピクセル値が異なっている。したがって、ＳＤマップをしきい値処理することにより、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）と、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）とを区別することができる。マスクＭＫ１を参照すると、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）にはピクセル値１が割り当てられており、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）には主にピクセル値０が割り当てられていることがわかる。また、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）は画像強度の損失が生じやすく、一方、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）は画像強度の損失は生じにくい。したがって、マスクＭＫ１のピクセル値１の領域は、画像強度の損失が生じる可能性が高い部位に対応しており、一方、マスクＭＫ１のピクセル値０の領域は、画像強度の損失が生じにくい部位に対応していることがわかる。
マスクＭＫ１〜ＭＫ６を作成した後、ステップＳＴ２２に進む。
【００５４】
ステップＳＴ２２では、拡散強調画像の準同型フィルタリングを実行する。以下にステップＳＴ２２について説明する。尚、どの拡散強調画像に対しても準同型フィルタリングの処理手順は同じである。したがって、以下の説明では、スキャンＳＣ１及びＳＣ２により取得された拡散強調画像のうち、代表して拡散強調画像ＤＷ１１を取り上げ、拡散強調画像ＤＷ１１に対して準同型フィルタリングを実行するときの手順について説明する。
【００５５】
ステップＳＴ２２では、先ず、モジュレーションマップ作成手段９３（図２参照）が、準同型フィルタリングに用いられるモジュレーションマップを作成する（図２３参照）。
【００５６】
図２３はモジュレーションマップの作成方法の説明図である。
モジュレーションマップ作成手段９３は、先ず、拡散強調画像ＤＷ１１のローパスフィルタリングを実行する。図２３では、拡散強調画像ＤＷ１１をローパスフィルタリングすることにより得られた画像（以下、「ローパスフィルタ画像」と呼ぶ）を、符号「ＤＷｌｐｆ」で示してある。
【００５７】
次に、モジュレーションマップ作成手段９３は、ローパスフィルタ画像ＤＷｌｐｆを反転し、反転画像を作成する。これにより、モジュレーションマップＤｍｏｄ１が作成される。
【００５８】
モジュレーションマップＤｍｏｄ１を作成した後、準同型フィルタ手段９４（図２参照）が、モジュレーションマップＤｍｏｄ１を用いて、拡散強調画像ＤＷ１１に対して準同型フィルタリングを実行する。具体的には、以下の式に基づいて準同型フィルタリングを実行する。
ＨＭ１１＝ＤＷ１１＊Ｄｍｏｄ１    ・・・（１）
【００５９】
式（１）のＨＭ１１は、準同型フィルタリングにより得られた画像（以下、「準同型フィルタ画像」と呼ぶ）を表している。
【００６０】
以下に、準同型フィルタリングについて説明する。尚、本形態では、準同型フィルタリングを実行する場合、マスクＭＫ１に基づいて準同型フィルタリングが実行される領域を限定している。以下に、この理由を明確にするために、先に、準同型フィルタリングが実行される領域を限定せずに式（１）を用いて準同型フィルタリングを実行する例を説明する。その後、マスクＭＫ１に基づいて準同型フィルタリングが実行される領域を限定し、式（１）を用いて準同型フィルタリングを実行する本形態の方法について説明する。
【００６１】
図２４は、準同型フィルタリングが実行される領域を限定せずに式（１）を用いて準同型フィルタリングを実行する例の説明図である。
拡散強調画像ＤＷ１１にモジュレーションマップＤｍｏｄ１を乗算することにより、拡散強調画像ＤＷ１１の準同型フィルタリングを実行する。これにより、準同型フィルタ画像ＨＭが得られる。
【００６２】
図２４では、画像強度の損失ＬＯ１、ＬＯ２、ＬＯ３のうち、ＬＯ１およびＬＯ２において、画像強度が復元された例が示されている。したがって、準同型フィルタリングにより、脾臓における画像強度の損失ＬＯ３は残っているが、肝臓の左葉および胃は画像強度が復元された準同型フィルタ画像ＨＭを取得することができる。
【００６３】
しかし、準同型フィルタリングは、画像強度を復元できる一方で、コントラストを増強させる作用を有している。したがって、準同型フィルタリングによりコントラストが変化してしまうという問題がある。特に、コントラストが変化すると、後述するステップＳＴ８において算出するＡＤＣの誤差が大きくなりやすいので、コントラストはできるだけ変化しないようにすることが望まれる。
【００６４】
そこで、本形態では、準同型フィルタ手段９４は、マスクＭＫ１（図２２参照）に基づいて、準同型フィルタリングを実行する領域を限定している。マスクＭＫ１に基づいて準同型フィルタリングを実行する領域を限定することにより、拡散強調画像ＤＷ１１の中で画像強度の損失があまり生じていない部分では、コントラストの変化を防止することができるという効果がある。以下に、この効果が得られる理由について、図２５を参照しながら説明する。
【００６５】
マスクＭＫ１は、ピクセル値１の領域と、ピクセル値０の領域とを有している。マスクＭＫ１のピクセル値１の領域は、画像強度の損失が生じている部位（例えば、肝臓の左葉および胃）に対応しているが、マスクＭＫ１のピクセル値０の領域は、画像強度の損失があまり生じていない部位（例えば、肝臓の右葉）に対応している。したがって、拡散強調画像ＤＷ１１の中で、マスクＭＫ１のピクセル値１の領域に対応する部分では、画像強度を復元するための準同型フィルタリングを実行する必要があるが、拡散強調画像ＤＷ１１の中で、マスクＭＫ１のピクセル値０の領域に対応する部分では、画像強度の損失があまり生じていないので準同型フィルタリングを行う必要がない。そこで、本形態では、拡散強調画像ＤＷ１１の中で、マスクＭＫ１のピクセル値１の領域に対応する部分では準同型フィルタリングを実行するが、マスクＭＫ１のピクセル値０の領域に対応する部分では、準同型フィルタリングを行わないようにしている。これにより、肝臓の右葉の大部分では準同型フィルタリングは実行されないので、コントラストの変化を防止することができる。したがって、準同型フィルタ画像ＨＭ１１は、肝臓の左葉および胃においては画像強度が復元されており、肝臓の右葉ではコントラストの変化が防止されている。このような理由から、本形態では、マスクＭＫ１（図２２参照）に基づいて、準同型フィルタリングを実行する領域を限定している。
【００６６】
上記の説明では、拡散強調画像ＤＷ１１の準同型フィルタ画像ＨＭ１１を得る方法について述べたが、他の拡散強調画像についても、同様の方法で準同型フィルタ画像を得ることができる。したがって、拡散強調画像ごとに準同型フィルタ画像を得ることができる。図２６および図２７に準同型フィルタ画像を概略的に示す。図２６では、ｂ＝ｂ１における準同型フィルタ画像が示されており、図２７では、ｂ＝ｂ２における準同型フィルタ画像が示されている。尚、図２６および図２７では、図面のスペースの制約上、スライス１の拡散強調画像から得られた準同型フィルタ画像のみが示されているが、他のスライスの拡散強調画像に対しても同様に準同型フィルタ画像が得られる。
【００６７】
上記のように、マスクＭＫ１を用いて準同型フィルタリングを実行することにより、肝臓の右葉においてコントラストの変化を防止しながら、肝臓の左葉や胃において画像強度を復元することができる。しかし、マスクＭＫ１のピクセル値１の領域においては準同型フィルタリングが実行されているので、肝臓の左葉や胃では、画像強度が復元される一方で、コントラストが増強する。したがって、図２５に示す準同型フィルタ画像ＨＭ１１を、画像強度が復元された復元画像として使用した場合、肝臓の左葉や胃においてＡＤＣの誤差が大きくなると考えられる。そこで、本形態では、ステップＳＴ２２において準同型フィルタ画像を作成した後、ステップＳＴ３〜ＳＴ６が実行される。ステップＳＴ３〜ＳＴ６では、準同型フィルタ画像を用いて、良好なコントラストを有するとともに画像強度が復元された復元画像が作成される。以下に、ステップＳＴ３〜ＳＴ６について説明する。尚、どの準同型フィルタ画像に対してもステップＳＴ３〜ＳＴ６の処理は同じである。したがって、以下の説明では、ステップＳＴ２２で得られた準同型フィルタ画像のうち、代表して準同型フィルタ画像ＨＭ１１を取り上げ、準同型フィルタ画像ＨＭ１１に対してステップＳＴ３〜ＳＴ６を実行するときの方法について説明する。
【００６８】
ステップＳＴ３では、差分画像作成手段９５（図２参照）が、以下の式（２）に基づいて、準同型フィルタ画像ＨＭ１１と拡散強調画像ＤＷ１１との差分画像ＤＩＦを求める。
ＤＩＦ＝ＨＭ１１−ＤＷ１１  ・・・（２）
【００６９】
図２８に、差分画像ＤＩＦの一例を示す。差分画像ＤＩＦは、準同型フィルタ画像ＨＭ１１と拡散強調画像ＤＷ１１とのピクセル値の差を表している。図２８から、差分画像ＤＩＦは、肝臓の左葉や胃において、大きいピクセル値を有していることがわかる。したがって、準同型フィルタリングを実行することにより、肝臓の左葉や胃では画像強度が復元されていることがわかる。差分画像ＤＩＦを求めた後、ステップＳＴ４に進む。
【００７０】
ステップＳＴ４では、ローパスフィルタリング手段９６（図２参照）が、マスクＭＫ１に対してローパスフィルタリングを実行する。図２９は、マスクＭＫ１をローパルフィルタリングするときの説明図である。マスクＭＫ１をローパスフィルタリングすることにより、マスクＭＫ１のピクセル値１と０との境界において滑らかな信号変化を実現することができる。以下では、ローパスフィルタリングが実行された後のマスクＭＫ１を、符号「ＭＫ１ｌｐｆ」で表すことにする。マスクＭＫ１ｌｐｆを作成した後、ステップＳＴ５に進む。
【００７１】
ステップＳＴ５では、ローパスフィルタリング手段９６が、差分画像ＤＩＦ（図２８参照）に対してローパスフィルタリングを実行する。図３０は、差分画像ＤＩＦをローパスフィルタリングするときの説明図である。図３０では、ローパルフィルタリングが実行された後の差分画像ＤＩＦを、符号「ＤＩＦｌｐｆ」で示す。差分画像ＤＩＦをローパスフィルタリングすることにより、ピクセル値を滑らかに変化させることができるので、準同型フィルタリングにより増強されてしまったコントラストを低減することができる。差分画像ＤＩＦｌｐｆは、差分画像ＤＩＦと比較すると、ピクセル値が滑らかに変化しており、コントラストが低減されている。差分画像ＤＩＦｌｐｆを作成した後、ステップＳＴ６に進む。
【００７２】
ステップＳＴ６では、復元画像作成手段９７（図２参照）が、差分画像ＤＩＦｌｐｆを用いて、画像強度が復元されるとともに良好なコントラストを有する復元画像を作成する。以下に、復元画像の作成方法について説明する。本形態では、復元画像は、以下の式（３）に従って作成される。
ＲＳ１１＝（１−ＭＫ１ｌｐｆ）＊ＤＷ１１＋ＭＫ１ｌｐｆ＊（ＤＷ１１＋ＤＩＦｌｐｆ）
・・・（３）
ここで、ＲＳ１１：復元画像
ＭＫ１ｌｐｆ：ローパスフィルタリングされた後のマスク
ＤＷ１１：拡散強調画像
ＤＩＦｌｐｆ：ローパルフィルタリングされた後の差分画像
【００７３】
以下に、式（３）について説明する。尚、式（３）の説明については、先ず、右辺第２項について先に説明し、次に、右辺第１項について説明する。
【００７４】
図３１は第２項の説明図である。
第２項の計算を行う場合、復元画像作成手段９７は、拡散強調画像ＤＷ１１に差分画像ＤＩＦｌｐｆを加算する。先に説明したように、差分画像ＤＩＦｌｐｆは、差分画像ＤＩＦと比較すると、準同型フィルタリングにより増強されたコントラストが低減されている。したがって、拡散強調画像ＤＷ１１に差分画像ＤＩＦｌｐｆを加算することにより、肝臓の左葉および胃の画像強度が復元され、且つコントラストも低減された加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）を得ることができる。
【００７５】
そして、復元画像作成手段９７は、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）に対してマスクＭＫ１ｌｐｆを乗算する。マスクＭＫ１ｌｐｆを乗算することにより、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）の中から、マスクＭＫ１ｌｐｆのピクセル値１の領域に対応する画像部分Ｄ２を取り出すことができる。尚、マスクＭＫ１ｌｐｆのピクセル値１の領域は、肝臓の左葉および胃を含んでいる。したがって、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）の中から取り出された画像部分Ｄ２には、画像強度が復元された肝臓の左葉および胃が含まれている。
【００７６】
図３２は第１項の説明図である。
第１項の（１−ＭＫ１ｌｐｆ）は、マスクＭＫ１ｌｐｆのピクセル値を反転させることにより得られる反転マスクを表している。マスクＭＫ１ｌｐｆのピクセル値１は、反転マスク（１−ＭＫ１ｌｐｆ）ではピクセル値０であり、マスクＭＫ１ｌｐｆのピクセル値０は、反転マスク（１−ＭＫ１ｌｐｆ）ではピクセル値１である。
【００７７】
第１項の計算を行う場合、復元画像作成手段９７は、拡散強調画像ＤＷ１１に反転マスク（１−ＭＫ１ｌｐｆ）を乗算する。反転マスク（１−ＭＫ１ｌｐｆ）を乗算することにより、拡散強調画像ＤＷ１１の中から、反転マスク（１−ＭＫ１ｌｐｆ）のピクセル値１の領域に対応する画像部分Ｄ１を取り出すことができる。
【００７８】
次に、復元画像作成手段９７は、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）から取り出した画像部分Ｄ２と、拡散強調画像ＤＷ１１から取り出した画像部分Ｄ１とを加算する。図３３は、画像部分Ｄ１およびＤ２を加算するときの説明図である。加算により復元画像ＲＳ１１を得ることができる。
【００７９】
復元画像ＲＳ１１は、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）から取り出された画像部分Ｄ２を含んでいるので、肝臓の左葉および胃では画像強度が復元され、且つコントラストも低減されている。
【００８０】
また、本形態では、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）から画像部分Ｄ２を取り出す場合、ローパスフィルタリングされたマスクＭＫ１ｌｐｆが使用されている。したがって、加算画像（ＤＷ１１＋ＤＩＦｌｐｆ）から取り出された画像部分Ｄ２のエッジにおいて、ピクセル値を滑らかに変化させることができる。更に、拡散強調画像ＤＷ１１から画像部分Ｄ１を取り出す場合、ローパスフィルタリングされたマスクＭＫ１ｌｐｆの反転マスク（１−ＭＫ１ｌｐｆ）が使用されている。したがって、拡散強調画像ＤＷ１１から取り出された画像部分Ｄ１のエッジにおいて、ピクセル値を滑らかに変化させることができる。このため、画像部分Ｄ１とＤ２とを加算しても、画像部分Ｄ１とＤ２との境界における位置ずれが原因で生じるアーチファクトを低減することができる。
【００８１】
上記の説明では、準同型フィルタ画像ＨＭ１１に基づいて復元画像ＲＳ１１を得る方法について述べられているが、他の準同型フィルタ画像についても、同様の方法で復元画像を得ることができる。したがって、準同型フィルタ画像ごとに復元画像を得ることができる。図３４および図３５に復元画像を概略的に示す。図３４では、ｂ＝ｂ１における復元画像が示されており、図３５では、ｂ＝ｂ２における復元画像が示されている。尚、図３４および図３５では、図面のスペースの制約上、スライス１の準同型フィルタ画像から得られた復元画像ＲＳ１１〜ＲＳ１ｎおよびＲＳ４１〜ＲＳ４ｎのみが示されているが、他のスライスの準同型フィルタ画像からも同様に復元画像が得られる。
【００８２】
本形態では、準同型フィルタリング（ステップＳＴ２２）の手法を用いることにより、肝臓の左葉および胃の画像強度が復元された復元画像（図３３参照）を得ることができる。
【００８３】
しかし、準同型フィルタリングを行っただけでは画像強度を十分に復元できない部位がある。図３３を参照すると、肝臓の左葉や胃の画像強度は復元されているが、脾臓では準同型フィルタリングを実行しても画像強度は十分に復元できていないことがわかる。そこで、本形態では、脾臓の画像強度を復元するため、ステップＳＴ７が実行される。
【００８４】
ステップＳＴ７では、合成画像作成手段９８（図２参照）が、復元画像を合成する。
図３６は、復元画像の合成方法の一例を示す図である。
先ず、合成画像作成手段９８は、ステップＳＴ６で得られた復元画像の中から、スライス１の復元画像ＲＳ１１〜ＲＳ１ｎを選択する。そして、合成画像作成手段９８は、以下の式を用いて復元画像ＲＳ１１〜ＲＳ１ｎを合成する。
【数１】
2015208420.tif 000003
【００８５】
式（４）は、振幅二乗重み付け平均法により復元画像ＲＳ１１〜ＲＳ１ｎを重み付けし、重み付けされた復元画像ＲＳ１１〜ＲＳ１ｎの平均を計算している。
【００８６】
合成画像作成手段９８は、式（４）に従ってスライス１の復元画像ＲＳ１１〜ＲＳ１ｎを合成する。したがって、合成画像ＤＣ１を得ることができる。図３７に、合成画像ＤＣ１の一例を示す。
【００８７】
上記のように、準同型フィルタリングを実行しても脾臓の画像強度は十分に復元することは難しいので、復元画像ＲＳ１１〜ＲＳ１ｎは脾臓に画像強度の損失を持つことが多い（図３３参照）。しかし、復元画像ＲＳ１１〜ＲＳ１ｎを比較すると、一般的に、脾臓の中に現れる画像強度の損失の位置は、脾臓の動きなどが原因で異なっている。したがって、復元画像ＲＳ１１〜ＲＳ１ｎを合成することにより、或る復元画像の脾臓に現れる画像強度の損失を、別の復元画像によって復元することができる。したがって、図３７に示すように、脾臓の画像強度が復元された合成画像ＤＣ１を得ることができる。尚、本形態では、振幅二乗重み付け平均法を用いているが、画像強度を復元できるのであれば、振幅二乗重み付け平均法とは別の方法で復元画像を合成してもよい。例えば、復元画像ＲＳ１１〜ＲＳ１ｎのうちの一部の復元画像には重み付けを行うが、残りの復元画像には重み付けを行わず、重み付けされた復元画像と重み付けされていない復元画像とを合成してもよい。また、復元画像ＲＳ１１〜ＲＳ１ｎの相加平均法で復元画像ＲＳ１１〜ＲＳ１ｎを合成してもよい。
【００８８】
図３６では、ｂ＝ｂ１におけるスライス１の合成画像ＤＣ１を作成する例について説明したが、ｂ＝ｂ２におけるスライス１の合成画像についても、振幅二乗重み付け平均法により合成画像を取得することができる。図３８に、ｂ＝ｂ２におけるスライス１の合成画像ＤＣ４を概略的に示す。
【００８９】
尚、上記の説明では、スライス１の合成画像の作成方法について説明したが、他のスライスについても同様に合成画像が作成される。したがって、ｂ＝ｂ１における各スライスの合成画像と、ｂ＝ｂ２における各スライスの合成画像とが作成される。
合成画像を作成した後、ステップＳＴ８に進む。
【００９０】
ステップＳＴ８では、ＡＤＣマップ作成手段９９（図２参照）が、ｂ＝ｂ１における合成画像と、ｂ＝ｂ２における合成画像に基づいて、スライスごとにＡＤＣマップを作成する。例えば、スライス１のＡＤＣマップは、合成画像ＤＣ１およびＤＣ４を用いて作成することができる。各スライスのＡＤＣマップを作成した後、フローを終了する。
【００９１】
本形態では、ステップＳＴ７で復元画像を合成しているので、準同型フィルタリングでは十分に復元することができなかった脾臓の画像強度を復元することができる。
【００９２】
また、本形態では、マスクＭＫ１のピクセル値１の領域では、準同型フィルタリングが行われるが、マスクＭＫ１のピクセル値０の領域では、準同型フィルタリングが行われないようにしている。したがって、準同型フィルタリングが実行されない領域ではコントラストの変化を防止することができる。
【００９３】
また、本形態のシーケンスＳＣでは、２つのｂ値（ｂ＝ｂ１およびｂ＝ｂ２）が用いられているが、２つのｂ値のうちの一方は、ｂ＝０であってもよい。ｂ＝０に設定する場合は、図７（又は図１２）に示すシーケンスのＭＰＧが印加されないようにすればよい。また、ｂ値の数は２つに限定されることはなく、１つのｂ値のみを用いた例を実行してもよいし、３つ以上のｂ値を用いた例を実行してもよい。
【００９４】
次に、本形態の効果を検証するために、シミュレーションを実行した。以下に、シミュレーション結果について説明する。
【００９５】
図３９は、シミュレーション条件の説明図である。
シミュレーションでは、スピンエコーＥＰＩ画像（レファレンス画像）に対してランダムに特定された２５箇所（図３９に「＋」で示されている箇所）における信号値を表すシミュレーションデータを得た。
【００９６】
次に、シミュレーションデータを合成した。具体的には、以下の３つの方法で合成を行った。
（ａ）シミュレーションデータを相加平均法で合成した。
（ｂ）シミュレーションデータを振幅二乗重み付き平均法で合成した。
（ｃ）シミュレーションデータを準同型フィルタリングし、準同型フィルタリングされたデータを振幅二乗重み付き平均法で合成した。
【００９７】
図４０はシミュレーション結果を示す図である。
図４０は、上記の方法（ａ）、（ｂ）、および（ｃ）で合成されたデータの正規化二乗平均誤差（Normalized Root Mean Squared Error）を表すグラフである。横軸は、レファンレス画像に特定された２５個の各箇所を表している。縦軸は正規化二乗平均誤差である。図４０を参照すると、（ｃ）の方法が最も小さい誤差を示していることがわかる。したがって、準同型フィルタリングと振幅二乗重み付き平均法の組合せが画像強度の復元に有効であることがわかる。
【００９８】
また、ｂ＝０のスキャンＡ（ＮＥＸ＝４）と、ｂ＝１０００のスキャンＢ（ＮＥＸ＝４）とを実行し、両方のスキャンＡおよびＢにより得られたデータに基づいてＡＤＣマップを作成した。具体的には、以下の２つの方法でＡＤＣマップを作成した。
【００９９】
（ｍ１）スキャンＡにより得られた４枚の画像を相加平均法で合成し、合成画像を得た。同様に、スキャンＢにより得られた４枚の画像を相加平均法で合成し、合成画像を得た。そして、この２つの合成画像に基づいてＡＤＣマップを作成した。
（ｍ２）スキャンＡにより得られた４枚の画像を準同型フィルタリングし、４枚の準同型フィルタ画像ＨＭａ１〜ＨＭａ４を得た。そして、４枚の準同型フィルタ画像ＨＭａ１〜ＨＭａ４に基づいて４枚の復元画像ＲＳａ１〜ＲＳａ４を作成した。次に、４枚の復元画像ＲＳａ１〜ＲＳａ４を振幅二乗重み付け平均法で合成し、合成画像ＤＣａを得た。
【０１００】
同様に、スキャンＢにより得られた４枚の画像を準同型フィルタリングし、４枚の準同型フィルタ画像ＨＭｂ１〜ＨＭｂ４を得た。そして、４枚の準同型フィルタ画像ＨＭｂ１〜ＨＭｂ４に基づいて４枚の復元画像ＲＳｂ１〜ＲＳｂ４を作成した。次に、４枚の復元画像ＲＳｂ１〜ＲＳｂ４を振幅二乗重み付け平均法で合成し、合成画像ＤＣｂを得た。そして、この２つの合成画像に基づいてＡＤＣマップを作成した。
【０１０１】
図４１はＡＤＣマップを示す図である。
図４１（ａ）は、方法（ｍ１）を用いて作成されたＡＤＣマップを示し、図４１（ｂ）は、方法（ｍ２）を用いて作成されたＡＤＣマップを示す。
【０１０２】
図４１（ａ）では、肝臓の左葉や脾臓におけるＡＤＣ値は、拍動の影響を受けて、実際のＡＤＣ値よりも大きい値になる。これに対し、図４１（ｂ）では、肝臓の左葉や脾臓のＡＤＣ値は、図４１（ａ）のＡＤＣ値よりも小さくなっている。したがって、準同型フィルタリングと振幅二乗重み付け平均法との両方を用いることにより、高品質のＡＤＣマップが得られることがわかる。
【符号の説明】
【０１０３】
２  マグネット
３  テーブル
３ａ  クレードル
４  コイル
５  送信器
６  勾配磁場電源
７  受信器
８  コンピュータ
９  プロセッサ
１０  メモリ
１１  操作部
１２  表示部
１３  被検体
２１  ボア
９１  画像作成手段
９２  マスク作成手段
９３  モジュレーションマップ作成手段
９４  準同型フィルタ手段
９５  差分画像作成手段
９６  ローパスフィルタ手段
９７  復元画像作成手段
９８  合成画像作成手段
９９  ＡＤＣマップ作成手段
１００  ＭＲ装置
【図１】
2015208420.tif 000004
【図２】
2015208420.tif 000005
【図３】
2015208420.tif 000006
【図４】
2015208420.tif 000007
【図５】
2015208420.tif 000008
【図６】
2015208420.tif 000009
【図７】
2015208420.tif 000010
【図８】
2015208420.tif 000011
【図９】
2015208420.tif 000012
【図１０】
2015208420.tif 000013
【図１１】
2015208420.tif 000014
【図１２】
2015208420.tif 000015
【図１３】
2015208420.tif 000016
【図１４】
2015208420.tif 000017
【図１５】
2015208420.tif 000018
【図１６】
2015208420.tif 000019
【図１７】
2015208420.tif 000020
【図１８】
2015208420.tif 000021
【図１９】
2015208420.tif 000022
【図２０】
2015208420.tif 000023
【図２１】
2015208420.tif 000024
【図２２】
2015208420.tif 000025
【図２３】
2015208420.tif 000026
【図２４】
2015208420.tif 000027
【図２５】
2015208420.tif 000028
【図２６】
2015208420.tif 000029
【図２７】
2015208420.tif 000030
【図２８】
2015208420.tif 000031
【図２９】
2015208420.tif 000032
【図３０】
2015208420.tif 000033
【図３１】
2015208420.tif 000034
【図３２】
2015208420.tif 000035
【図３３】
2015208420.tif 000036
【図３４】
2015208420.tif 000037
【図３５】
2015208420.tif 000038
【図３６】
2015208420.tif 000039
【図３７】
2015208420.tif 000040
【図３８】
2015208420.tif 000041
【図３９】
2015208420.tif 000042
【図４０】
2015208420.tif 000043
【図４１】
2015208420.tif 000044