JP6554717B2 - 磁気共鳴装置およびプログラム - Google Patents

Description

本発明は、拡散情報を含むデータを収集する磁気共鳴装置、およびこの磁気共鳴装置に適用されるプログラムに関する。

近年、拡散強調イメージング（ＤＷＩ：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：励起回数）がｎ回の例（ＮＥＸ＝ｎ）が示されている。スキャンＳＣ１は、ｎ個の期間Ｖ１〜Ｖｎを有している。各期間では、３つのシーケンスＥ１、Ｅ２、およびＥ３が実行される。したがって、スキャンＳＣ１では、３つのシーケンスＥ１、Ｅ２、およびＥ３のセットがｎ回実行される。以下に、スキャンＳＣ１の期間Ｖ１〜Ｖｎで実行されるシーケンスＥ１、Ｅ２、Ｅ３について説明する。

スキャンＳＣ１では、先ず、期間Ｖ１において、シーケンスＥ１が実行される（図７参照）。

図７は、シーケンスＥ１の一例を示す図である。
シーケンスＥ１は、以下のパルス（１）〜（４）を有している。
（１）ＲＦパルス
（１ａ）スライスを励起するための励起パルスα１
（１ｂ）スピンを再収束させるためのリフォーカスパルスα２
（２）スライス選択方向の勾配パルス
（２ａ）スライスを選択するためのスライス選択勾配パルスＧ_ｓ１およびＧ_ｓ２
（３）周波数エンコード方向の勾配パルス
（３ａ）勾配パルスＧ_ｆｓ
（３ｂ）ｂ＝ｂ１のＭＰＧ
（３ｃ）エコーを収集するための読出し勾配パルスＧ_ｒｅ
（４）位相エンコード方向の勾配パルス
（４ａ）勾配パルスＧ_ｐｓ
（４ｂ）位相エンコード勾配パルスＧｐ

シーケンスＥ１は、スライス１を励起し、スライス１の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。本形態では、シーケンスＥ１は、スライス１の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間の全データを１回のショットで収集する（シングルショット法）。画像作成手段９１（図２参照）は、シーケンスＥ１により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＥ１が終了した後、シーケンスＥ２が実行される。

シーケンスＥ２も、シーケンスＥ１と同様に、図７に示すシーケンスチャートで表される。ただし、シーケンスＥ２は、スライス２が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＥ２は、スライス２を励起し、スライス２の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＥ２により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＥ２が終了した後、シーケンスＥ３が実行される。

シーケンスＥ３も、シーケンスＥ１と同様に、図７に示すシーケンスチャートで表される。ただし、シーケンスＥ３は、スライス３が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＥ３は、スライス３を励起し、スライス３の拡散強調画像（ｂ＝ｂ１）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＥ３により収集されたデータに基づいて、拡散強調画像を作成する。

したがって、期間Ｖ１においてシーケンスＥ１、Ｅ２、およびＥ３を実行することにより、スライス１〜３の各々の拡散強調画像を取得することができる。図８に、期間Ｖ１のシーケンスＥ１〜Ｅ３により得られた拡散強調画像ＤＷ_１１、ＤＷ_２１、ＤＷ_３１（ｂ＝ｂ１）を概略的に示す。
期間Ｖ１の後、次の期間Ｖ２においてシーケンスが実行される。

期間Ｖ２でも、期間Ｖ１と同様に、シーケンスＥ１、Ｅ２、およびＥ３が実行される。したがって、期間Ｖ２でも、スライス１〜３の各々の拡散強調画像（ｂ＝ｂ１）を取得することができる。図９に、期間Ｖ２のシーケンスＥ１〜Ｅ３により得られた拡散強調画像を概略的に示す。図９では、期間Ｖ１で取得された拡散強調画像と、期間Ｖ２で取得された拡散強調画像とを区別するために、期間Ｖ２で取得された拡散強調画像を符号「ＤＷ_１２」、「ＤＷ_２２」、「ＤＷ_３２」で示してある。

以下同様に、３つのシーケンスＥ１、Ｅ２、およびＥ３のセットが繰り返し実行される。そして、スキャンＳＣ１の最後の期間Ｖｎにおいて、シーケンスＥ１、Ｅ２、およびＥ３が実行される。図１０に、期間ＶｎにおいてシーケンスＥ１、Ｅ２、およびＥ３を実行することにより取得された拡散強調画像を概略的に示す。図１０では、期間Ｖｎで取得された拡散強調画像を符号「ＤＷ_１ｎ」、「ＤＷ_２ｎ」、「ＤＷ_３ｎ」で示してある。
次に、スキャンＳＣ２について説明する。

図１１はスキャンＳＣ２の説明図である。
図１１は、図６と同様に、ＮＥＸ＝ｎの例を示している。スキャンＳＣ２は、ｎ個の期間Ｗ１〜Ｗｎを有している。各期間では、３つのシーケンスＦ１、Ｆ２、およびＦ３が実行される。したがって、スキャンＳ２では、３つのシーケンスＦ１、Ｆ２、およびＦ３のセットがｎ回実行される。以下に、スキャンＳＣ２の期間Ｗ１〜Ｗｎで実行されるシーケンスＦ１、Ｆ２、Ｆ３について説明する。
スキャンＳＣ２では、先ず、期間Ｗ１において、シーケンスＦ１が実行される（図１２参照）。

図１２は、シーケンスＦ１の一例を示す図である。
シーケンスＦ１は、シーケンスＥ１（図７参照）と比較すると、ｂ＝ｂ２に設定されたＭＰＧを有しているが、その他のパルスは、シーケンスＥ１と同じである。したがって、シーケンスＦ１は、スライス１を励起し、スライス１の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ１により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＦ１が終了した後、シーケンスＦ２が実行される。

シーケンスＦ２も、シーケンスＦ１と同様に、図１２に示すシーケンスチャートで表される。ただし、シーケンスＦ２は、スライス２が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＦ２は、スライス２を励起し、スライス２の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ２により収集されたデータに基づいて、拡散強調画像を作成する。
シーケンスＦ２が終了した後、シーケンスＦ３が実行される。

シーケンスＦ３も、シーケンスＦ１と同様に、図１２に示すシーケンスチャートで表される。ただし、シーケンスＦ３は、スライス３が励起されるようにＲＦパルスの周波数が設定されている。シーケンスＦ３は、スライス３を励起し、スライス３の拡散強調画像（ｂ＝ｂ２）の再構成に必要なｋ空間のデータを収集する。画像作成手段９１は、シーケンスＦ３により収集されたデータに基づいて、拡散強調画像を作成する。

以下同様に、期間Ｗ２〜Ｗｎについても、３つのシーケンスＦ１、Ｆ２、およびＦ３のセットが繰り返し実行される。したがって、スキャンＳＣ２を実行することにより、１枚のスライスにつき、ｎ枚の拡散強調画像を取得することができる。図１３に、期間Ｗ１〜Ｗｎにおいて取得された拡散強調画像を概略的に示す。

このようにして、スキャンＳＣ１およびＳＣ２が実行される。スキャンＳＣ１およびＳＣ２は被検体が息止めしている間に実行してもよいし、呼吸同期法を用いて実行してもよい。

ただし、ＭＰＧの印加中に、拍動、蠕動などが原因で器官が動いてしまうと、器官の動きの影響を受けて、拡散強調画像に画像強度の損失（信号損失）が現れることが知られている。図１４に、拡散強調画像ＤＷ_１１の一例を示す。図１４では、拡散強調画像ＤＷ_１１の中の一部（例えば、肝臓の左葉、胃、および脾臓）に、画像強度の損失ＬＯ１、ＬＯ２、およびＬＯ３が現れている。このような画像強度の損失は、後述するステップＳＴ８においてＡＤＣマップを作成する場合、ＡＤＣ値の誤差の原因になるので、拡散強調画像ＤＷ_１１の中で画像強度の損失が生じている部分の画像強度はできるだけ復元する必要がある。そこで、本形態では、スキャンＳＣ１およびＳＣ２を実行した後、画像強度を復元するためにステップＳＴ２を実行する。以下に、ステップＳＴ２について説明する。尚、ステップＳＴ２は、ステップＳＴ２１およびＳＴ２２を有しているので、ステップＳＴ２１およびＳＴ２２について順に説明する。

ステップＳＴ２１では、マスク作成手段９２（図２参照）が、後述するステップＳＴ２２で使用するマスク（図２２参照）を作成する。以下に、マスクの作成方法について説明する。

先ず、マスク作成手段９２は、ｂ＝ｂ１のスキャンＳＣ１により得られた拡散強調画像の中から、スライス１の拡散強調画像ＤＷ_１１〜ＤＷ_１ｎを選択する。図１５に、選択された拡散強調画像ＤＷ_１１〜ＤＷ_１ｎを示す。

次に、マスク作成手段９２は、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎのピクセル値の標準偏差を表すＳＤマップを作成する。図１６に、ＳＤマップの作成方法を概略的に示す。

マスク作成手段９２は、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎから同じ位置のピクセルのピクセル値を選択し、選択したピクセル値の標準偏差を求める。例えば、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎのピクセルＰａのピクセル値の標準偏差を求める場合、マスク作成手段９２は、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎからピクセルＰａのピクセル値Ｘａ１〜Ｘａｎを選択する。そして、ピクセル値Ｘａ１〜Ｘａｎの標準偏差σａを計算する。このようにして計算された標準偏差σａが、ＳＤマップ１のピクセルＰａにおける値として用いられる。

上記の説明では、ＳＤマップのピクセルＰａにおける標準偏差σａを求める方法について説明したが、ＳＤマップの他のピクセルにおける標準偏差も同じ方法で計算される。例えば、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎのピクセルＰｂにおけるピクセル値Ｘｂ１〜Ｘｂｎの標準偏差σｂを計算することにより、ＳＤマップ１のピクセルＰｂにおける値を得ることができる。同様に、拡散強調画像ＤＷ_１１〜ＤＷ_１ｎのピクセルＰｃにおけるピクセル値Ｘｃ１〜Ｘｃｎの標準偏差σｃを計算することにより、ＳＤマップ１のピクセルＰｃにおける値を得ることができる。

したがって、ピクセルごとに標準偏差を計算することにより、ｂ＝ｂ１におけるスライス１のＳＤマップ１を得ることができる。

図１５および図１６では、スライス１のＳＤマップ１を作成する手順について示されているが、他のスライス２およびスライス３についても同様の方法でＳＤマップを得ることができる。図１７に、スライス２のＳＤマップ２およびスライス３のＳＤマップ３を概略的に示す。このように、スライスごとにＳＤマップを得ることができる。

尚、図１５〜図１７では、ｂ＝ｂ１のスキャンＳＣ１により取得された拡散強調画像を用いてＳＤマップを作成した例について示されているが、ｂ＝ｂ２のスキャンＳＣ２により取得された拡散強調画像に対しても、同様の方法でスライスごとにＳＤマップが作成される。図１８に、ｂ＝ｂ２において、スライスごとに取得されたＳＤマップを概略的に示す。
したがって、全部で６枚のＳＤマップが作成される。図１９に、ＳＤマップの一例を示す。

肝臓の左葉、胃、脾臓など、心臓の近くに位置している部位は、拍動の影響を受けて動いてしまうので、肝臓の左葉、胃、脾臓などの部位のピクセル値は時間とともに変化する傾向がある。一方、肝臓の右葉は、心臓から遠いので、心臓の拍動の影響をあまり受けない。したがって、肝臓の右葉のピクセル値の時間変化はそれほど大きくならない傾向がある。このため、ＳＤマップから、心臓の拍動の影響を受けやすい部位と、心臓の拍動の影響をあまり受けない部位とをおおまかに区別することができる。

次に、マスク作成手段９２は、ＳＤマップに基づいてマスクを作成する。図２０および図２１に、ＳＤマップに基づいてマスクを作成する方法を概略的に示す。マスク作成手段９２は、ＳＤマップ１〜６の各々に対してしきい値処理ＴＨを行う。しきい値ＴＨは、ＳＤマップのピクセル値のうち、しきい値より小さいピクセル値を「０」に設定し、一方、しきい値以上のピクセル値を「１」に設定する。したがって、各ＳＤマップは、「０」の領域と、「１」の領域とに分けられる。このように、各ＳＤマップを「０」の領域と「１」の領域とに分けることにより、マスクＭＫ１〜ＭＫ６を作成することができる。図２２に、マスクＭＫ１〜ＭＫ６のうち、代表して、マスクＭＫ１の具体的な一例を示す。図２２では、ピクセル値１の領域が白で示されており、ピクセル値０の領域が黒で示されている。

上記のように、ＳＤマップは、心臓の拍動などの影響を受けやすい部位のピクセル値と、心臓の拍動などの影響をあまり受けない部位のピクセル値が異なっている。したがって、ＳＤマップをしきい値処理することにより、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）と、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）とを区別することができる。マスクＭＫ１を参照すると、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）にはピクセル値１が割り当てられており、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）には主にピクセル値０が割り当てられていることがわかる。また、心臓の拍動の影響を受けやすい部位（肝臓の左葉、胃、脾臓など）は画像強度の損失が生じやすく、一方、心臓の拍動の影響をあまり受けない部位（肝臓の右葉など）は画像強度の損失は生じにくい。したがって、マスクＭＫ１のピクセル値１の領域は、画像強度の損失が生じる可能性が高い部位に対応しており、一方、マスクＭＫ１のピクセル値０の領域は、画像強度の損失が生じにくい部位に対応していることがわかる。
マスクＭＫ１〜ＭＫ６を作成した後、ステップＳＴ２２に進む。

ステップＳＴ２２では、拡散強調画像の準同型フィルタリングを実行する。以下にステップＳＴ２２について説明する。尚、どの拡散強調画像に対しても準同型フィルタリングの処理手順は同じである。したがって、以下の説明では、スキャンＳＣ１及びＳＣ２により取得された拡散強調画像のうち、代表して拡散強調画像ＤＷ_１１を取り上げ、拡散強調画像ＤＷ_１１に対して準同型フィルタリングを実行するときの手順について説明する。

ステップＳＴ２２では、先ず、モジュレーションマップ作成手段９３（図２参照）が、準同型フィルタリングに用いられるモジュレーションマップを作成する（図２３参照）。

図２３はモジュレーションマップの作成方法の説明図である。
モジュレーションマップ作成手段９３は、先ず、拡散強調画像ＤＷ_１１のローパスフィルタリングを実行する。図２３では、拡散強調画像ＤＷ_１１をローパスフィルタリングすることにより得られた画像（以下、「ローパスフィルタ画像」と呼ぶ）を、符号「ＤＷ_ｌｐｆ」で示してある。

次に、モジュレーションマップ作成手段９３は、ローパスフィルタ画像ＤＷ_ｌｐｆを反転し、反転画像を作成する。これにより、モジュレーションマップＤ_ｍｏｄ１が作成される。

モジュレーションマップＤ_ｍｏｄ１を作成した後、準同型フィルタ手段９４（図２参照）が、モジュレーションマップＤ_ｍｏｄ１を用いて、拡散強調画像ＤＷ_１１に対して準同型フィルタリングを実行する。具体的には、以下の式に基づいて準同型フィルタリングを実行する。
ＨＭ_１１＝ＤＷ_１１＊Ｄ_ｍｏｄ１ ・・・（１）

式（１）のＨＭ_１１は、準同型フィルタリングにより得られた画像（以下、「準同型フィルタ画像」と呼ぶ）を表している。

以下に、準同型フィルタリングについて説明する。尚、本形態では、準同型フィルタリングを実行する場合、マスクＭＫ１に基づいて準同型フィルタリングが実行される領域を限定している。以下に、この理由を明確にするために、先に、準同型フィルタリングが実行される領域を限定せずに式（１）を用いて準同型フィルタリングを実行する例を説明する。その後、マスクＭＫ１に基づいて準同型フィルタリングが実行される領域を限定し、式（１）を用いて準同型フィルタリングを実行する本形態の方法について説明する。

図２４は、準同型フィルタリングが実行される領域を限定せずに式（１）を用いて準同型フィルタリングを実行する例の説明図である。
拡散強調画像ＤＷ_１１にモジュレーションマップＤ_ｍｏｄ１を乗算することにより、拡散強調画像ＤＷ_１１の準同型フィルタリングを実行する。これにより、準同型フィルタ画像ＨＭが得られる。

図２４では、画像強度の損失ＬＯ１、ＬＯ２、ＬＯ３のうち、ＬＯ１およびＬＯ２において、画像強度が復元された例が示されている。したがって、準同型フィルタリングにより、脾臓における画像強度の損失ＬＯ３は残っているが、肝臓の左葉および胃は画像強度が復元された準同型フィルタ画像ＨＭを取得することができる。

しかし、準同型フィルタリングは、画像強度を復元できる一方で、コントラストを増強させる作用を有している。したがって、準同型フィルタリングによりコントラストが変化してしまうという問題がある。特に、コントラストが変化すると、後述するステップＳＴ８において算出するＡＤＣの誤差が大きくなりやすいので、コントラストはできるだけ変化しないようにすることが望まれる。

そこで、本形態では、準同型フィルタ手段９４は、マスクＭＫ１（図２２参照）に基づいて、準同型フィルタリングを実行する領域を限定している。マスクＭＫ１に基づいて準同型フィルタリングを実行する領域を限定することにより、拡散強調画像ＤＷ_１１の中で画像強度の損失があまり生じていない部分では、コントラストの変化を防止することができるという効果がある。以下に、この効果が得られる理由について、図２５を参照しながら説明する。

マスクＭＫ１は、ピクセル値１の領域と、ピクセル値０の領域とを有している。マスクＭＫ１のピクセル値１の領域は、画像強度の損失が生じている部位（例えば、肝臓の左葉および胃）に対応しているが、マスクＭＫ１のピクセル値０の領域は、画像強度の損失があまり生じていない部位（例えば、肝臓の右葉）に対応している。したがって、拡散強調画像ＤＷ_１１の中で、マスクＭＫ１のピクセル値１の領域に対応する部分では、画像強度を復元するための準同型フィルタリングを実行する必要があるが、拡散強調画像ＤＷ_１１の中で、マスクＭＫ１のピクセル値０の領域に対応する部分では、画像強度の損失があまり生じていないので準同型フィルタリングを行う必要がない。そこで、本形態では、拡散強調画像ＤＷ_１１の中で、マスクＭＫ１のピクセル値１の領域に対応する部分では準同型フィルタリングを実行するが、マスクＭＫ１のピクセル値０の領域に対応する部分では、準同型フィルタリングを行わないようにしている。これにより、肝臓の右葉の大部分では準同型フィルタリングは実行されないので、コントラストの変化を防止することができる。したがって、準同型フィルタ画像ＨＭ_１１は、肝臓の左葉および胃においては画像強度が復元されており、肝臓の右葉ではコントラストの変化が防止されている。このような理由から、本形態では、マスクＭＫ１（図２２参照）に基づいて、準同型フィルタリングを実行する領域を限定している。

上記の説明では、拡散強調画像ＤＷ_１１の準同型フィルタ画像ＨＭ_１１を得る方法について述べたが、他の拡散強調画像についても、同様の方法で準同型フィルタ画像を得ることができる。したがって、拡散強調画像ごとに準同型フィルタ画像を得ることができる。図２６および図２７に準同型フィルタ画像を概略的に示す。図２６では、ｂ＝ｂ１における準同型フィルタ画像が示されており、図２７では、ｂ＝ｂ２における準同型フィルタ画像が示されている。尚、図２６および図２７では、図面のスペースの制約上、スライス１の拡散強調画像から得られた準同型フィルタ画像のみが示されているが、他のスライスの拡散強調画像に対しても同様に準同型フィルタ画像が得られる。

上記のように、マスクＭＫ１を用いて準同型フィルタリングを実行することにより、肝臓の右葉においてコントラストの変化を防止しながら、肝臓の左葉や胃において画像強度を復元することができる。しかし、マスクＭＫ１のピクセル値１の領域においては準同型フィルタリングが実行されているので、肝臓の左葉や胃では、画像強度が復元される一方で、コントラストが増強する。したがって、図２５に示す準同型フィルタ画像ＨＭ_１１を、画像強度が復元された復元画像として使用した場合、肝臓の左葉や胃においてＡＤＣの誤差が大きくなると考えられる。そこで、本形態では、ステップＳＴ２２において準同型フィルタ画像を作成した後、ステップＳＴ３〜ＳＴ６が実行される。ステップＳＴ３〜ＳＴ６では、準同型フィルタ画像を用いて、良好なコントラストを有するとともに画像強度が復元された復元画像が作成される。以下に、ステップＳＴ３〜ＳＴ６について説明する。尚、どの準同型フィルタ画像に対してもステップＳＴ３〜ＳＴ６の処理は同じである。したがって、以下の説明では、ステップＳＴ２２で得られた準同型フィルタ画像のうち、代表して準同型フィルタ画像ＨＭ_１１を取り上げ、準同型フィルタ画像ＨＭ_１１に対してステップＳＴ３〜ＳＴ６を実行するときの方法について説明する。

ステップＳＴ３では、差分画像作成手段９５（図２参照）が、以下の式（２）に基づいて、準同型フィルタ画像ＨＭ_１１と拡散強調画像ＤＷ_１１との差分画像ＤＩＦを求める。
ＤＩＦ＝ＨＭ_１１−ＤＷ_１１ ・・・（２）

図２８に、差分画像ＤＩＦの一例を示す。差分画像ＤＩＦは、準同型フィルタ画像ＨＭ_１１と拡散強調画像ＤＷ_１１とのピクセル値の差を表している。図２８から、差分画像ＤＩＦは、肝臓の左葉や胃において、大きいピクセル値を有していることがわかる。したがって、準同型フィルタリングを実行することにより、肝臓の左葉や胃では画像強度が復元されていることがわかる。差分画像ＤＩＦを求めた後、ステップＳＴ４に進む。

ステップＳＴ４では、ローパスフィルタリング手段９６（図２参照）が、マスクＭＫ１に対してローパスフィルタリングを実行する。図２９は、マスクＭＫ１をローパルフィルタリングするときの説明図である。マスクＭＫ１をローパスフィルタリングすることにより、マスクＭＫ１のピクセル値１と０との境界において滑らかな信号変化を実現することができる。以下では、ローパスフィルタリングが実行された後のマスクＭＫ１を、符号「ＭＫ１_ｌｐｆ」で表すことにする。マスクＭＫ１_ｌｐｆを作成した後、ステップＳＴ５に進む。

ステップＳＴ５では、ローパスフィルタリング手段９６が、差分画像ＤＩＦ（図２８参照）に対してローパスフィルタリングを実行する。図３０は、差分画像ＤＩＦをローパスフィルタリングするときの説明図である。図３０では、ローパルフィルタリングが実行された後の差分画像ＤＩＦを、符号「ＤＩＦ_ｌｐｆ」で示す。差分画像ＤＩＦをローパスフィルタリングすることにより、ピクセル値を滑らかに変化させることができるので、準同型フィルタリングにより増強されてしまったコントラストを低減することができる。差分画像ＤＩＦ_ｌｐｆは、差分画像ＤＩＦと比較すると、ピクセル値が滑らかに変化しており、コントラストが低減されている。差分画像ＤＩＦ_ｌｐｆを作成した後、ステップＳＴ６に進む。

ステップＳＴ６では、復元画像作成手段９７（図２参照）が、差分画像ＤＩＦ_ｌｐｆを用いて、画像強度が復元されるとともに良好なコントラストを有する復元画像を作成する。以下に、復元画像の作成方法について説明する。本形態では、復元画像は、以下の式（３）に従って作成される。
ＲＳ_１１＝（１−ＭＫ１_ｌｐｆ）＊ＤＷ_１１＋ＭＫ１_ｌｐｆ＊（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）
・・・（３）
ここで、ＲＳ_１１：復元画像
ＭＫ１_ｌｐｆ：ローパスフィルタリングされた後のマスク
ＤＷ_１１：拡散強調画像
ＤＩＦ_ｌｐｆ：ローパルフィルタリングされた後の差分画像

以下に、式（３）について説明する。尚、式（３）の説明については、先ず、右辺第２項について先に説明し、次に、右辺第１項について説明する。

図３１は第２項の説明図である。
第２項の計算を行う場合、復元画像作成手段９７は、拡散強調画像ＤＷ_１１に差分画像ＤＩＦ_ｌｐｆを加算する。先に説明したように、差分画像ＤＩＦ_ｌｐｆは、差分画像ＤＩＦと比較すると、準同型フィルタリングにより増強されたコントラストが低減されている。したがって、拡散強調画像ＤＷ_１１に差分画像ＤＩＦ_ｌｐｆを加算することにより、肝臓の左葉および胃の画像強度が復元され、且つコントラストも低減された加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）を得ることができる。

そして、復元画像作成手段９７は、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）に対してマスクＭＫ１_ｌｐｆを乗算する。マスクＭＫ１_ｌｐｆを乗算することにより、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）の中から、マスクＭＫ１_ｌｐｆのピクセル値１の領域に対応する画像部分Ｄ２を取り出すことができる。尚、マスクＭＫ１_ｌｐｆのピクセル値１の領域は、肝臓の左葉および胃を含んでいる。したがって、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）の中から取り出された画像部分Ｄ２には、画像強度が復元された肝臓の左葉および胃が含まれている。

図３２は第１項の説明図である。
第１項の（１−ＭＫ１_ｌｐｆ）は、マスクＭＫ１_ｌｐｆのピクセル値を反転させることにより得られる反転マスクを表している。マスクＭＫ１_ｌｐｆのピクセル値１は、反転マスク（１−ＭＫ１_ｌｐｆ）ではピクセル値０であり、マスクＭＫ１_ｌｐｆのピクセル値０は、反転マスク（１−ＭＫ１_ｌｐｆ）ではピクセル値１である。

第１項の計算を行う場合、復元画像作成手段９７は、拡散強調画像ＤＷ_１１に反転マスク（１−ＭＫ１_ｌｐｆ）を乗算する。反転マスク（１−ＭＫ１_ｌｐｆ）を乗算することにより、拡散強調画像ＤＷ_１１の中から、反転マスク（１−ＭＫ１_ｌｐｆ）のピクセル値１の領域に対応する画像部分Ｄ１を取り出すことができる。

次に、復元画像作成手段９７は、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）から取り出した画像部分Ｄ２と、拡散強調画像ＤＷ_１１から取り出した画像部分Ｄ１とを加算する。図３３は、画像部分Ｄ１およびＤ２を加算するときの説明図である。加算により復元画像ＲＳ_１１を得ることができる。

復元画像ＲＳ_１１は、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）から取り出された画像部分Ｄ２を含んでいるので、肝臓の左葉および胃では画像強度が復元され、且つコントラストも低減されている。

また、本形態では、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）から画像部分Ｄ２を取り出す場合、ローパスフィルタリングされたマスクＭＫ１_ｌｐｆが使用されている。したがって、加算画像（ＤＷ_１１＋ＤＩＦ_ｌｐｆ）から取り出された画像部分Ｄ２のエッジにおいて、ピクセル値を滑らかに変化させることができる。更に、拡散強調画像ＤＷ_１１から画像部分Ｄ１を取り出す場合、ローパスフィルタリングされたマスクＭＫ１_ｌｐｆの反転マスク（１−ＭＫ１_ｌｐｆ）が使用されている。したがって、拡散強調画像ＤＷ_１１から取り出された画像部分Ｄ１のエッジにおいて、ピクセル値を滑らかに変化させることができる。このため、画像部分Ｄ１とＤ２とを加算しても、画像部分Ｄ１とＤ２との境界における位置ずれが原因で生じるアーチファクトを低減することができる。

上記の説明では、準同型フィルタ画像ＨＭ_１１に基づいて復元画像ＲＳ_１１を得る方法について述べられているが、他の準同型フィルタ画像についても、同様の方法で復元画像を得ることができる。したがって、準同型フィルタ画像ごとに復元画像を得ることができる。図３４および図３５に復元画像を概略的に示す。図３４では、ｂ＝ｂ１における復元画像が示されており、図３５では、ｂ＝ｂ２における復元画像が示されている。尚、図３４および図３５では、図面のスペースの制約上、スライス１の準同型フィルタ画像から得られた復元画像ＲＳ_１１〜ＲＳ_１ｎおよびＲＳ_４１〜ＲＳ_４ｎのみが示されているが、他のスライスの準同型フィルタ画像からも同様に復元画像が得られる。

本形態では、準同型フィルタリング（ステップＳＴ２２）の手法を用いることにより、肝臓の左葉および胃の画像強度が復元された復元画像（図３３参照）を得ることができる。

しかし、準同型フィルタリングを行っただけでは画像強度を十分に復元できない部位がある。図３３を参照すると、肝臓の左葉や胃の画像強度は復元されているが、脾臓では準同型フィルタリングを実行しても画像強度は十分に復元できていないことがわかる。そこで、本形態では、脾臓の画像強度を復元するため、ステップＳＴ７が実行される。

ステップＳＴ７では、合成画像作成手段９８（図２参照）が、復元画像を合成する。
図３６は、復元画像の合成方法の一例を示す図である。
先ず、合成画像作成手段９８は、ステップＳＴ６で得られた復元画像の中から、スライス１の復元画像ＲＳ_１１〜ＲＳ_１ｎを選択する。そして、合成画像作成手段９８は、以下の式を用いて復元画像ＲＳ_１１〜ＲＳ_１ｎを合成する。

式（４）は、振幅二乗重み付け平均法により復元画像ＲＳ_１１〜ＲＳ_１ｎを重み付けし、重み付けされた復元画像ＲＳ_１１〜ＲＳ_１ｎの平均を計算している。

合成画像作成手段９８は、式（４）に従ってスライス１の復元画像ＲＳ_１１〜ＲＳ_１ｎを合成する。したがって、合成画像ＤＣ_１を得ることができる。図３７に、合成画像ＤＣ_１の一例を示す。

上記のように、準同型フィルタリングを実行しても脾臓の画像強度は十分に復元することは難しいので、復元画像ＲＳ_１１〜ＲＳ_１ｎは脾臓に画像強度の損失を持つことが多い（図３３参照）。しかし、復元画像ＲＳ_１１〜ＲＳ_１ｎを比較すると、一般的に、脾臓の中に現れる画像強度の損失の位置は、脾臓の動きなどが原因で異なっている。したがって、復元画像ＲＳ_１１〜ＲＳ_１ｎを合成することにより、或る復元画像の脾臓に現れる画像強度の損失を、別の復元画像によって復元することができる。したがって、図３７に示すように、脾臓の画像強度が復元された合成画像ＤＣ_１を得ることができる。尚、本形態では、振幅二乗重み付け平均法を用いているが、画像強度を復元できるのであれば、振幅二乗重み付け平均法とは別の方法で復元画像を合成してもよい。例えば、復元画像ＲＳ_１１〜ＲＳ_１ｎのうちの一部の復元画像には重み付けを行うが、残りの復元画像には重み付けを行わず、重み付けされた復元画像と重み付けされていない復元画像とを合成してもよい。また、復元画像ＲＳ_１１〜ＲＳ_１ｎの相加平均法で復元画像ＲＳ_１１〜ＲＳ_１ｎを合成してもよい。

図３６では、ｂ＝ｂ１におけるスライス１の合成画像ＤＣ_１を作成する例について説明したが、ｂ＝ｂ２におけるスライス１の合成画像についても、振幅二乗重み付け平均法により合成画像を取得することができる。図３８に、ｂ＝ｂ２におけるスライス１の合成画像ＤＣ_４を概略的に示す。

尚、上記の説明では、スライス１の合成画像の作成方法について説明したが、他のスライスについても同様に合成画像が作成される。したがって、ｂ＝ｂ１における各スライスの合成画像と、ｂ＝ｂ２における各スライスの合成画像とが作成される。
合成画像を作成した後、ステップＳＴ８に進む。

ステップＳＴ８では、ＡＤＣマップ作成手段９９（図２参照）が、ｂ＝ｂ１における合成画像と、ｂ＝ｂ２における合成画像に基づいて、スライスごとにＡＤＣマップを作成する。例えば、スライス１のＡＤＣマップは、合成画像ＤＣ_１およびＤＣ_４を用いて作成することができる。各スライスのＡＤＣマップを作成した後、フローを終了する。

本形態では、ステップＳＴ７で復元画像を合成しているので、準同型フィルタリングでは十分に復元することができなかった脾臓の画像強度を復元することができる。

また、本形態では、マスクＭＫ１のピクセル値１の領域では、準同型フィルタリングが行われるが、マスクＭＫ１のピクセル値０の領域では、準同型フィルタリングが行われないようにしている。したがって、準同型フィルタリングが実行されない領域ではコントラストの変化を防止することができる。

また、本形態のシーケンスＳＣでは、２つのｂ値（ｂ＝ｂ１およびｂ＝ｂ２）が用いられているが、２つのｂ値のうちの一方は、ｂ＝０であってもよい。ｂ＝０に設定する場合は、図７（又は図１２）に示すシーケンスのＭＰＧが印加されないようにすればよい。また、ｂ値の数は２つに限定されることはなく、１つのｂ値のみを用いた例を実行してもよいし、３つ以上のｂ値を用いた例を実行してもよい。

次に、本形態の効果を検証するために、シミュレーションを実行した。以下に、シミュレーション結果について説明する。

図３９は、シミュレーション条件の説明図である。
シミュレーションでは、スピンエコーＥＰＩ画像（レファレンス画像）に対してランダムに特定された２５箇所（図３９に「＋」で示されている箇所）における信号値を表すシミュレーションデータを得た。

次に、シミュレーションデータを合成した。具体的には、以下の３つの方法で合成を行った。
（ａ）シミュレーションデータを相加平均法で合成した。
（ｂ）シミュレーションデータを振幅二乗重み付き平均法で合成した。
（ｃ）シミュレーションデータを準同型フィルタリングし、準同型フィルタリングされたデータを振幅二乗重み付き平均法で合成した。

図４０はシミュレーション結果を示す図である。
図４０は、上記の方法（ａ）、（ｂ）、および（ｃ）で合成されたデータの正規化二乗平均誤差（Normalized Root Mean Squared Error）を表すグラフである。横軸は、レファンレス画像に特定された２５個の各箇所を表している。縦軸は正規化二乗平均誤差である。図４０を参照すると、（ｃ）の方法が最も小さい誤差を示していることがわかる。したがって、準同型フィルタリングと振幅二乗重み付き平均法の組合せが画像強度の復元に有効であることがわかる。

また、ｂ＝０のスキャンＡ（ＮＥＸ＝４）と、ｂ＝１０００のスキャンＢ（ＮＥＸ＝４）とを実行し、両方のスキャンＡおよびＢにより得られたデータに基づいてＡＤＣマップを作成した。具体的には、以下の２つの方法でＡＤＣマップを作成した。

（ｍ１）スキャンＡにより得られた４枚の画像を相加平均法で合成し、合成画像を得た。同様に、スキャンＢにより得られた４枚の画像を相加平均法で合成し、合成画像を得た。そして、この２つの合成画像に基づいてＡＤＣマップを作成した。
（ｍ２）スキャンＡにより得られた４枚の画像を準同型フィルタリングし、４枚の準同型フィルタ画像ＨＭ_ａ１〜ＨＭ_ａ４を得た。そして、４枚の準同型フィルタ画像ＨＭ_ａ１〜ＨＭ_ａ４に基づいて４枚の復元画像ＲＳ_ａ１〜ＲＳ_ａ４を作成した。次に、４枚の復元画像ＲＳ_ａ１〜ＲＳ_ａ４を振幅二乗重み付け平均法で合成し、合成画像ＤＣａを得た。

同様に、スキャンＢにより得られた４枚の画像を準同型フィルタリングし、４枚の準同型フィルタ画像ＨＭ_ｂ１〜ＨＭ_ｂ４を得た。そして、４枚の準同型フィルタ画像ＨＭ_ｂ１〜ＨＭ_ｂ４に基づいて４枚の復元画像ＲＳ_ｂ１〜ＲＳ_ｂ４を作成した。次に、４枚の復元画像ＲＳ_ｂ１〜ＲＳ_ｂ４を振幅二乗重み付け平均法で合成し、合成画像ＤＣｂを得た。そして、この２つの合成画像に基づいてＡＤＣマップを作成した。

図４１はＡＤＣマップを示す図である。
図４１（ａ）は、方法（ｍ１）を用いて作成されたＡＤＣマップを示し、図４１（ｂ）は、方法（ｍ２）を用いて作成されたＡＤＣマップを示す。

図４１（ａ）では、肝臓の左葉や脾臓におけるＡＤＣ値は、拍動の影響を受けて、実際のＡＤＣ値よりも大きい値になる。これに対し、図４１（ｂ）では、肝臓の左葉や脾臓のＡＤＣ値は、図４１（ａ）のＡＤＣ値よりも小さくなっている。したがって、準同型フィルタリングと振幅二乗重み付け平均法との両方を用いることにより、高品質のＡＤＣマップが得られることがわかる。

２ マグネット
３ テーブル
３ａ クレードル
４ コイル
５ 送信器
６ 勾配磁場電源
７ 受信器
８ コンピュータ
９ プロセッサ
１０ メモリ
１１ 操作部
１２ 表示部
１３ 被検体
２１ ボア
９１ 画像作成手段
９２ マスク作成手段
９３ モジュレーションマップ作成手段
９４ 準同型フィルタ手段
９５ 差分画像作成手段
９６ ローパスフィルタ手段
９７ 復元画像作成手段
９８ 合成画像作成手段
９９ ＡＤＣマップ作成手段
１００ ＭＲ装置

Claims (14)

1. 複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行するスキャン手段と、
   前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成手段と、
   前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための準同型フィルタリングを実行する第１のフィルタ手段であって、前記第１の画像の中の第２の部分において前記準同型フィルタリングを実行しない第１のフィルタ手段と、
   前記第１の画像に基づいて、前記準同型フィルタリングするためのモジュレーションマップを作成するモジュレーションマップ作成手段と
   を有し、
   前記第１のフィルタ手段は、
   前記モジュレーションマップを用いて、前記準同型フィルタリングを実行する、磁気共鳴装置。
2. 前記第１の画像の前記第１の部分に対応する第１の領域と、前記第１の画像の前記第２の部分に対応する第２の領域とを含むマスクを作成するマスク作成手段を有し、
   前記第１のフィルタ手段は、
   前記第１の画像の中で、前記マスクの第１の領域に対応する部分において、前記準同型フィルタリングを実行するが、前記第１の画像の中で、前記マスクの第２の領域に対応する部分において、前記準同型フィルタリングを実行しない、請求項１に記載の磁気共鳴装置。
3. 前記マスク作成手段は、
   前記複数の第１の画像のピクセル値の標準偏差を表す標準偏差マップを作成し、前記標準偏差マップに基づいて前記マスクを作成する、請求項２に記載の磁気共鳴装置。
4. 前記第１の画像を準同型フィルタリングすることにより得られた第２の画像と、前記第１の画像とを差分することにより、差分画像を作成する差分画像作成手段と、
   前記差分画像をローパスフィルタリングするための第２のフィルタ手段と、
   ローパスフィルタリングされた差分画像と前記第１の画像とを加算することによりに得られる第１の加算画像を作成する第２の画像作成手段と、
   を有する、請求項２又は３に記載の磁気共鳴装置。
5. 前記第２のフィルタ手段は、
   前記マスクの前記第１の領域と前記第２の領域との境界におけるピクセル値が滑らかに変化するように、前記マスクをローパスフィルタリングし、
   前記第２の画像作成手段は、
   ローパスフィルタリングされたマスクを用いて、前記第１の加算画像から第１の画像部分を取り出し、
   前記ローパスフィルタリングされたマスクの反転マスクを用いて、前記第１の画像から第２の画像部分を取り出し、
   前記第１の画像部分と前記第２の画像部分とを加算することにより、第３の画像を作成する、請求項４に記載の磁気共鳴装置。
6. 前記第２の画像作成手段は複数の第３の画像を作成し、
   前記複数の第３の画像を合成することにより合成画像を作成する第３の画像作成手段を有する、請求項５に記載の磁気共鳴装置。
7. 前記第３の画像作成手段は、
   前記複数の第３の画像の各々を重み付けし、重み付けされた複数の第３の画像に基づいて前記合成画像を作成する、請求項６に記載の磁気共鳴装置。
8. 前記第３の画像作成手段は、
   振幅二乗重み付き平均法を用いて前記複数の第２の画像の各々を重み付けする、請求項７に記載の磁気共鳴装置。
9. 前記合成画像と、前記合成画像とは別の合成画像とに基づいて、ＡＤＣマップを作成するＡＤＣマップ作成手段を有する、請求項６〜８のうちのいずれか一項に記載の磁気共鳴装置。
10. 前記スキャン手段は、
    前記スライスからデータを収集するための第２のシーケンスを複数回実行し、
    前記別の合成画像は、
    前記第２のシーケンスを複数回実行することにより得られたデータに基づいて作成されている、請求項９に記載の磁気共鳴装置。
11. 前記第１のシーケンスは、拡散強調を行うための第１の勾配パルスであって、拡散強調の強さを表すｂ値が第１の値に設定された第１の勾配パルスを含み、
    前記第２のシーケンスは、拡散強調を行うための第２の勾配パルスであって、拡散強調の強さを表すｂ値が第２の値に設定された第２の勾配パルスを含む、請求項１０に記載の磁気共鳴装置。
12. 前記第１のシーケンスは、拡散強調を行うための第１の勾配パルスであって、拡散強調の強さを表すｂ値が第１の値に設定された第１の勾配パルスを含み、
    前記第２のシーケンスは、拡散強調を行うための勾配パルスを含んでいない、請求項１０に記載の磁気共鳴装置。
13. 前記第１の画像の前記第１の部分は肝臓の左葉および胃を含んでおり、前記第１の画像の前記第２の部分は肝臓の右葉を含んでいる、請求項１〜１２のうちのいずれか一項に記載の磁気共鳴装置。
14. 複数の部位を含む撮影部位を横切るスライスから拡散情報を含むデータを収集するための第１のシーケンスを複数回実行する磁気共鳴装置に適用されるプログラムであって、
    前記第１のシーケンスを複数回実行することにより収集されたデータに基づいて、拡散情報を含む複数の第１の画像を作成する第１の画像作成処理と、
    前記第１の画像ごとに、前記第１の画像のうちの画像強度の損失が生じた第１の部分において画像強度を復元するための準同型フィルタリングを実行する第１のフィルタ処理であって、前記第１の画像の中の第２の部分において前記準同型フィルタリングを実行しない第１のフィルタ処理と、
    前記第１の画像に基づいて、前記準同型フィルタリングするためのモジュレーションマップを作成するモジュレーションマップ作成処理と
    をコンピュータに実行させるためのプログラムであり、
    前記第１のフィルタ処理は、
    前記モジュレーションマップを用いて、前記準同型フィルタリングを実行する、プログラム。