(57)【要約】  関心のある領域の磁気共鳴（ＭＲ）画像を作成する方法が提供される。この方法は、関心のある領域の初期ＭＲ画像を取得するステップであって、初期ＭＲ画像は、領域にわたる位置での感ＭＲ性の物理的性質の値をマッピングしている、取得するステップと、領域にわたる感ＭＲ性の物理的性質の値における推定された不確実性の対応するマップを決定するステップと、領域の強調ＭＲ画像を計算するステップであって、強調ＭＲ画像は、初期画像の各位置において、感ＭＲ性の物理的性質のそれぞれの値とそれぞれの推定された不確実性とを組み合わせる関数の値をマッピングしており、関数は、相対的に高く推定された不確実性を有する位置よりも、相対的に低く推定された不確実性を有する位置に、より高い加重を適用している、計算するステップとを含む。
2019513515.tif 000002
【特許請求の範囲】
【請求項１】
関心のある領域の磁気共鳴（ＭＲ）画像を作成する方法であって、
前記関心のある領域の初期ＭＲ画像を取得するステップであって、前記初期ＭＲ画像は、前記領域にわたる位置で、感ＭＲ性の物理的性質の値をマッピングしている、取得するステップと、
前記領域にわたる前記感ＭＲ性の物理的性質の値において推定された不確実性の対応するマップを決定するステップと、
前記領域の強調ＭＲ画像を計算するステップであって、前記強調ＭＲ画像は、前記初期画像の各位置において、前記感ＭＲ性の物理的性質のそれぞれの値とそれぞれの前記推定された不確実性とを組み合わせる関数の値をマッピングし、前記関数は、より高い強調を、相対的に高く推定された不確実性を有する位置よりも、相対的に低く推定された不確実性を有する位置へ適用する、計算するステップとを含む、方法。
【請求項２】
前記関数は、逆指数強調を前記推定された不確実性に適用する、請求項１に記載の方法。
【請求項３】
前記推定された不確実性は、前記領域にわたる前記感ＭＲ性の物理的性質の値の逆モデルから決定される、請求項１または２に記載の方法。
【請求項４】
前記決定するステップにおいて、前記逆モデルは、最小二乗法を用いて前記感ＭＲ性の物理的性質の前記値に適合される、請求項３に記載の方法。
【請求項５】
前記推定された不確実性は、前記感ＭＲ性の物理的性質の前記値の分散である、請求項１から４のいずれか一項に記載の方法。
【請求項６】
前記感ＭＲ性の物理的性質が、見かけの拡散係数、Ｔ１、Ｔ２、またはプロトン密度である、請求項１から５のいずれか一項に記載の方法。
【請求項７】
前記感ＭＲ性の物理的性質が見かけの拡散係数であり、前記関数は、より高い強調を、相対的に高い見かけの拡散係数を有する位置よりも、相対的に低い見かけの拡散係数を有する位置に適用する、請求項１から６のいずれか一項に記載の方法。
【請求項８】
前記関数は、前記見かけの拡散係数に逆指数強調を適用する、請求項７に記載の方法。
【請求項９】
前記領域にわたる位置での見かけの拡散係数の値をマッピングしている前記初期ＭＲ画像は、異なるｂ値で前記領域にわたる複数の拡散強調画像を取得し、そこから前記初期ＭＲ画像を導出することによって取得される、請求項７または８に記載の方法。
【請求項１０】
関心のある立体を通る、それぞれ平行なスライスである複数の関心のある領域の各々について請求項１から９のいずれか一項に記載の方法を実行して、それにより対応する強調ＭＲ画像を各スライスに対して取得するステップと、
前記強調ＭＲ画像を投影三次元画像として表示するステップとを含む、３Ｄ  ＭＲ画像データを表示する方法。
【請求項１１】
請求項１から１０のいずれか一項に記載の方法を実行するように動作可能に構成されたコンピュータシステム。
【請求項１２】
請求項１１に記載の前記コンピュータシステムを有する磁気共鳴撮像装置。
【請求項１３】
請求項１から１０のいずれか一項に記載の方法を実行するためのコンピュータプログラム。
【請求項１４】
請求項１３に記載のプログラムを実行するコンピュータプログラム製品。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴（ＭＲ）撮像法およびシステムに関し、特に、ＭＲ拡散強調画像に関する。
【背景技術】
【０００２】
全身拡散強調ＭＲ撮像（ＷＢＤＷＩ）は、前立腺［１、２］、乳房［３］、リンパ腫［４］および骨髄腫［５］を含む原発性悪性腫瘍の範囲からの転移性骨疾患の病期診断および応答評価にとって魅力的なツールになりつつある。ＷＢＤＷＩは、外因性造影剤を必要とせずに、病変組織と健常組織との間に非侵襲性のコントラストを提供し、臨床医に単一の放射線ビューにおいて骨格全体の転移性の程度を検査する能力を提供する。異なる拡散強調（ｂ値）で同じ解剖学的位置での複数の画像を取得することにより、ＷＢＤＷＩは見かけの拡散係数（ＡＤＣ）、腫瘍の細胞充実度の間接的定量を提供するバイオマーカーの測定も提供する［６］。ＡＤＣの定量化は、治療後のＡＤＣの増加が腫瘍の壊死を示すと考えられる治療応答の評価に大きな有用性を示している［７］。
【０００３】
ＷＢＤＷＩの比較的最近の適用は、ｃｏｍｐｕｔｅｄ  ＤＷＩ（ｃＤＷＩ）［８ａ、８ｂ］である。この技術は、直接測定から得られるものよりも高いｂ値で画像を合成することにより、良性組織と病変組織との間のコントラストを改善する。これは、病変の誤解の一般的な原因であるＴ２  ｓｈｉｎｅ−ｔｈｒｏｕｇｈ効果を排除し、高いｂ値画像の信号対ノイズ比を改善するのに役立つことができる。この方法論はまた、全身の容量的腫瘍量（ｔＤＶ）および広範囲の病変ＡＤＣ（ｇＡＤＣ）の定量化を提供する、ＷＢＤＷＩの研究における疾患の半自動セグメンテーションに対しても利用されており、この両者は治療応答を評価するためのバイオマーカーとして使用してもよい［９］。ｃＤＷＩは、高ｂ値画像を合成的に作成するためにＡＤＣ測定と一緒に使用される、モデルフィッティングから取得または推定されるｂ＝０画像（Ｓ０）の使用に依存している。Ｓ０画像は、画像化された組織のＴ１、Ｔ２およびプロトン密度の影響を受け、結果としてのｃＤＷＩのコントラストは必ずしも組織間のＡＤＣの差に起因するとは限らない。さらに、Ｓ０画像は、一般的にコイルの感度に大きく影響され、視界全体で厳密な信号の不均一性をもたらし、これは、収集ステーション間の「シグナル−ステップ」アーティファクトとして一般的に明示される［１］。そのような画像の不均一性は、施設内および施設間の絶対的な信号強度の違いに起因したＷＢＤＷＩにおけるセグメンテーション技術の標準化を妨げる。その他、指数関数的強調ＤＷ撮像（ｅＤＷＩ）を介して、これらの困難性に対する解決策を提供しようと試みている［１０］。この方法論は、Ｓ０を全視界にわたって一定の値に設定することで、合成高ｂ値画像を純粋にＡＤＣで生成したコントラストで得ることができる。しかしながら、この技術は、ＡＤＣ計算が撮像ノイズの影響を強く受ける分野には一般的には鋭さに欠けるので、本質的に低コントラスト−ノイズ比（ＣＮＲ）を有している。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明は、上述した方法の欠点に対処することを目的とする。
【課題を解決するための手段】
【０００５】
一般的には、本発明は、ＡＤＣなどの感ＭＲ性の物理的性質の値における不確実性を推定し、推定された不確実性を使用して、強調を適用し、ＭＲ画像におけるノイズ低減を生じさせることを提案する。そのようなアプローチは、例えば、ｅＤＷＩ画像において観察されたＣＮＲを改善するために使用することができる。
【０００６】
従って、第１の態様では、本発明は、
関心のある領域の初期ＭＲ画像を取得するステップであって、初期ＭＲ画像は、領域にわたる位置で、感ＭＲ性の物理的性質の値をマッピングしている、取得するステップと、
領域にわたる感ＭＲ性の物理的性質の値において、推定された不確実性の対応するマップを決定するステップと、
領域の強調ＭＲ画像を計算するステップであって、強調ＭＲ画像は、初期画像の各位置において、感ＭＲ性の物理的性質のそれぞれの値とそれぞれの推定された不確実性とを組み合わせる関数の値をマッピングしており、関数は、相対的に高く推定された不確実性を有する位置よりも相対的に低く推定された不確実性を有する位置に対してより高い強調を適用する、計算するステップとを含む関心のある領域の磁気共鳴（ＭＲ）画像を生成する方法を提供する。
【０００７】
上述したように、強調ＭＲ画像は、改善されたＣＮＲ特性を有することができる。加えて、強調ＭＲ画像は、ｃＤＷＩにおいて観察され得る信号の不均一性のような、信号の不均一性の影響を受けにくい。
【０００８】
本発明の第２の態様は、３ＤのＭＲ画像データを表示する方法であって、この方法は、
関心のある立体（ｖｏｌｕｍｅ）を通る、それぞれ平行なスライスである複数の関心のある領域の各々について第１の態様の方法を実行し、それにより各スライスに対して対応する強調ＭＲ画像を得るステップと、
投影された立体画像として、例えば、関心のある立体の最大強度投影として強調ＭＲ画像を表示するステップとを含む、３ＤのＭＲ画像データを表示する方法を提供する。
【０００９】
有利には、この方法は、同じ患者に対して複数の訪問にわたって取得される投影された立体画像の視覚的な比較を改善することができ、それ故、施設間で全身疾患のセグメント化を統一するために使用することができる。
【００１０】
本発明の第３の態様は、第１または第２の態様の方法を実行するように動作可能に構成されたコンピュータシステムを提供する。例えば、関心のある領域の強調ＭＲ画像を生成するためのコンピュータシステムは、
関心のある領域の初期ＭＲ画像を保管するためのコンピュータ可読媒体であって、初期ＭＲ画像は領域にわたる位置で、感ＭＲ性の物理的性質の値をマッピングしている、コンピュータ可読媒体と、
領域にわたる感ＭＲ性の物理的性質の値における推定された不確実性の対応するマップを決定し、領域の強調ＭＲ画像を計算するための１つまたは複数のプロセッサであって、この強調ＭＲ画像は、初期画像の各位置で、感ＭＲ性の物理的性質のそれぞれの値とそれぞれの推定された不確実性とを結びつける関数の値をマッピングしており、この関数は、相対的に高く推定された不確実性を有する位置よりも相対的に低く推定された不確実性を有する位置により高い強調を適用している、１つまたは複数のプロセッサとを含んでいてもよい。
【００１１】
このコンピュータシステムは、初期ＭＲ画像、推定された不確実性のマップ、強調ＭＲ画像、および（このシステムが第２の態様の方法を実行するように動作可能に構成されている場合）投影立体画像のいずれか１つまたは複数を表示するためのディスプレイを含むことができる。
【００１２】
本発明のさらなる態様は、第３の態様のコンピュータシステムを有する磁気共鳴撮像装置を提供する。
【００１３】
本発明のさらなる態様は、第１または第２の態様の方法を実行するためのコンピュータプログラム、および第１または第２態様の方法を実行するためのコンピュータプログラムを備えるコンピュータプログラム製品を提供する。
【００１４】
本発明の任意の特徴をここで説明する。これらは単独でも、または本発明の任意の態様の任意の組み合わせにおいて適用可能である。
【００１５】
この関数は、推定された不確実性に逆指数加重を適用し得る。
【００１６】
都合のよいことに、推定された不確実性は、領域にわたる感ＭＲ性の物理的性質の値の逆モデルから決定することができる。例えば、決定において、逆モデルは、最小二乗法を用いて感ＭＲ性の物理的性質の値にフィッティングされてもよい。
【００１７】
推定される不確実性は、感ＭＲ性の物理的性質の値の分散であってもよい。
【００１８】
感ＭＲ性の物理的性質は、見かけの拡散係数、Ｔ１、Ｔ２、またはプロトン密度であり得る。感ＭＲ性の物理的性質が見かけの拡散係数である場合、関数は、相対的に高い見かけの拡散係数を有する位置よりも、相対的に低い見かけの拡散係数を有する位置に、より高い強調を適用してもよく、例えばこの関数は、見かけの拡散係数に逆指数の強調を適用することができる。従って、この機能は、特に、転移性骨疾患の領域を同定するのに役立ち得る。あるいは、感ＭＲ性の物理的性質が見かけの拡散係数である場合、関数は、相対的に高い見かけの拡散係数を有する位置よりも、相対的に低い見かけの拡散係数を有する位置に、より低い強調を適用してもよく、例えばこの関数は、見かけの拡散係数に指数関数的な強調を適用することができる。このようにして、この機能は壊死の領域を特定するのに役立つ。領域にわたる位置での見かけの拡散係数の値をマッピングしている初期ＭＲ画像は、異なるｂ値で領域にわたる複数の拡散強調画像を取得し、そこから初期ＭＲ画像を導出することによって取得することができる。
【図面の簡単な説明】
【００１９】
本発明の実施形態を、添付の図面を参照して例として以下に説明する。
【００２０】
【図１】（ａ）は所与の断面積および所与のｂ値に対する拡散強調ＭＲ画像示しており、（ｂ）は複数のそのような画像についてのｂ値に対するＩｎＳ（ｂ）のプロットを示しており、プロット内のＩｎＳ（ｂ）に対する値は（ａ）に示された画素で得られ、その画素でのＡＤＣおよびＡＤＣにおける標準偏差σＡＤＣについて計算された値を示しており、（ｃ）は他の画素に対するＡＤＣ計算を繰り返すことにより得られた断面に対するＡＤＣ画像を示しており、（ｄ）は他の画素に対するσＡＤＣ計算を繰り返すことによって得られた断面に対するσＡＤＣのマップを示している。
【図２】（ａ）は図１（ｃ）のＡＤＣ画像を示しており、（ｂ）は図１（ｄ）のσＡＤＣマップを示しており、（ｃ）はＡＤＣ画像から導かれたｅｘｐ｛−ｂｃ・ＡＤＣ｝のマップを示しており、（ｄ）はσＡＤＣマップから導かれたｅｘｐ｛−ａ・σＡＤＣ｝のマップを示しており、（ｅ）はＳｎｃ（ａ，ｂｃ）＝ｅｘｐ｛−ｂｃ・ＡＤＣ｝・ｅｘｐ｛−ａ・σＡＤＣ｝の強調ＭＲ画像を示している。
【図３】第１の患者について、（ａ）は基準軸断面に対するσＡＤＣマップを示しており、（ｂ）は断面に対する対応するＡＤＣ画像を示しており、（ｃ）は断面に対するノイズ補正されて、指数関数的に加重された対応するＤＷＩ（ｎｉｃｅＤＷＩ）画像を示しており、画像を計算するために使用された（ａ，ｂｃ）の値が示されており、（ｄ）は断面に対する対応する高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を示しており、（ｅ）は病変組織とバックグラウンド組織との間のコントラストを視覚的に最大化するｂｃおよびａ値でのｎｉｃｅＤＷＩ画像から得られた最大強度投影（ＭＩＰ）立体を示しており、（ｆ）は通常取得される高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を用いて生成された対応するＭＩＰ立体を示している。
【図４】第２の患者について、（ａ）は基準軸断面に対するσＡＤＣマップを示しており、（ｂ）は断面に対する対応するＡＤＣ画像を示しており、（ｃ）は断面に対する対応するｎｉｃｅＤＷＩ画像を示しており、画像を計算するために使用された（ａ，ｂｃ）の値が示されており、（ｄ）は断面に対する対応する高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を示しており、（ｅ）は病変組織とバックグラウンド組織との間のコントラストを視覚的に最大化するｂｃおよびａ値でのｎｉｃｅＤＷＩ画像から得られた最大強度投影（ＭＩＰ）立体を示しており、（ｆ）は通常取得される高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を用いて生成された対応するＭＩＰ立体を示している。
【図５】第３の患者について、（ａ）は基準軸断面に対するσＡＤＣマップを示しており、（ｂ）は断面に対する対応するＡＤＣ画像を示しており、（ｃ）は断面に対する対応するｎｉｃｅＤＷＩ画像を示しており、画像を計算するために使用された（ａ、ｂｃ）の値が示されており、（ｄ）は断面に対する対応する高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を示しており、（ｅ）は病変組織とバックグラウンド組織との間のコントラストを視覚的に最大化するｂｃおよびａ値でのｎｉｃｅＤＷＩ画像から得られた最大強度投影（ＭＩＰ）立体を示しており、（ｆ）は通常取得される高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を用いて生成された対応するＭＩＰ立体を示している。
【図６】バックグラウンドノイズ分散の計算の例を示している。最も低いｂ値画像の対数からのデータは、ガウスの混合としてモデル化される（ａ）。最も低い平均値を有するガウス成分は、最も低いｂ値のｍａｇｎｉｔｕｄｅ画像上のバックグラウンドの領域（斜線領域の外側）を分類するために使用される（ｂ）。バックグラウンド領域内のデータは、画像に対するσを推定するために使用される（ｃ）。
【図７】複数の取得を通した画像ノイズ分散の直接計算（「真の」と表示され、各ケースにおける左の画像）および平均画像からの推定（「推定された」と表示され、各ケースにおける右の画像）を介したσＡＤＣの計算の４つの例を示している。散布図は、真のσＡＤＣ計算と推定されたσＡＤＣ計算との間のピクセル単位の相関を示しており、破線は、均等の線を表す。一般的に、両方の技術の間で、特に骨のような動きに影響されない領域（矢印Ｂ）において良好な一致がある。動きに影響を受けるこれらの領域は、特に腸、腎臓および脾臓において、計算上の相違をもたらし得る。これらの領域では、推定値は真の値を過小評価する傾向がある（矢印Ｓ）。動きに対する非感受性は有利となり得る。右上の例では、真の計算は、肝病変に対するσＡＤＣの高い値を明らかにしており、一方、推定された値は、周囲の健康な肝臓（矢印Ｌ）とのコントラストをはるかに低く改善している。
【図８】（ｉ）従来の高ｂ値画像と、（ｉｉ）繰り返しのｂ値取得が利用可能なｎｉｃｅＤＷＩと、（ｉｉｉ）繰り返しのｂ値取得が不可能なσＡＤＣの推定値とで取得されたｎｉｃｅＤＷＩ画像を示している。手法（ｉｉｉ）は、手法（ｉｉ）と比較して、精巣（矢印Ｅ）のような健康な組織と腫瘍（矢印Ｔ）との間で類似のコントラストを提供することができる。手法（ｉｉ）および（ｉｉｉ）の両方は、手法（ｉ）に対して、健康組織と病変組織との間のコントラストを改善している。
【発明を実施するための形態】
【００２１】
以下の説明は、好ましい例示的な実施形態のみを提供するものであり、本発明の範囲、適用性または構成を限定するものではない。むしろ、好ましい例示的な実施形態の以下の説明は、当業者に発明の好ましい例示的な実施形態の実現を可能にする説明を提供するが、本発明の範囲から逸脱することなく要素の機能および配置において様々な変更を加えることができると理解される。
【００２２】
本明細書で開示するように、「コンピュータ可読媒体」という用語は、リードオンリーメモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）、磁気ＲＡＭ、コアメモリ、磁気ディスク記憶媒体、光学記憶媒体、フラッシュメモリ装置、および／または情報を記憶するための他の機械可読媒体を含む、データを記憶するための一つまたは複数の装置を表す。「コンピュータ可読媒体」という用語は、携帯用または固定された記憶装置に限定されず、光学記憶装置、ワイヤレスチャネル、および命令および／またはデータを記憶、格納または搬送することができる様々な他の媒体を含む。
【００２３】
さらに、実施形態は、ハードウェア、ソフトウェア、ファームウェア、ミドルウェア、マイクロコード、ハードウェア記述言語、またはそれらの任意の組み合わせによって実施されてもよい。ソフトウェア、ファームウェア、ミドルウェアまたはマイクロコードで実施される場合、必要なタスクを実行するプログラムコードまたはコードセグメントは、記憶媒体などの機械可読媒体に格納されてもよい。プロセッサは、必要なタスクを実施することができる。コードセグメントは、プロシージャ、ファンクション、サブプログラム、プログラム、ルーチン、サブルーチン、モジュール、ソフトウェアパッケージ、クラス、または命令、データ構造、またはプログラム文の任意の組み合わせを表すことができる。コードセグメントは、情報、データ、引数、パラメータ、またはメモリ内容を渡す、および／または受け取ることによって、別のコードセグメントまたはハードウェア回路に結合されてもよい。情報、引数、パラメータ、データなどは、メモリ共有、メッセージパッシング、トークンパッシング、ネットワーク送信などの任意の適切な手段を介して渡され、転送され、または送信され得る。
【００２４】
理論
磁気共鳴画像における信号強度
【数１】
2019513515.tif 000003
に対するモデルは、
【数２】
2019513515.tif 000004
によって与えることができ、ここで、
【数３】
2019513515.tif 000005
は、感ＭＲ性の、組織の物理的性質（例えば、ＡＤＣ、スピン−格子緩和Ｔ１、スピン−スピン緩和Ｔ２、プロトン密度）を意味し、βはＭＲＩシーケンスパラメータ（例えば、ｂ値、エコー時間、フリップ角）を表し、
【数４】
2019513515.tif 000006
は撮像システムに対するモデルであり、εはいくつかの確率分布（例えば、ライス分布［１１］）からの統計ノイズである。一般的に、αはすべて不明であるが、使用者は、ＭＲ撮像システムの所与の制約内で、βを選択することができる。従って、組織パラメータの定量化は、逆撮像モデルｍβ−１の同定に依存しており、それにより、αのマップは、別個のスキャナ設定βｉでそれぞれ取得された一連のＭ個の画像Ｓ＝｛Ｓｉ：ｉ＝１．．．Ｍ｝から推定することができ、
【数５】
2019513515.tif 000007
画像ノイズεからもたらされる推定の不確実性はηαで表される。
【００２５】
計算された画像化は、想定される撮像モデル、
【数６】
2019513515.tif 000008
を使用して、シーケンスパラメータβＣのいくつかの任意のセットで画像ＳＣを合成するために、推定された組織パラメータを使用するポストプロセッシング技術である。
【００２６】
そのようなプロセスは、臨床スキャナで物理的に利用可能でないβの値で画像コントラストを生成することができ、信号−ノイズ比および結果の画像の幾何学的忠実性を改善することができる［８ａ、８ｂ］。このアプローチは、撮像モデルｍを、出力画像で望ましい性質を有するいくつかの任意の関数
【数７】
2019513515.tif 000009
で置き換えることにより、さらに一般化することができる［１２］：
【数８】
2019513515.tif 000010
【００２７】
しかしながら、本発明はまた、推定されたパラメータの不確実性ηαを計算された画像強調に組み込む。より具体的には、推定された不確実性の分散Σα’を使用してもよく、それによりパラメータ推定に大きな不確実性がある計算された画素は、推定の信頼性が高いこれらの画素と比較して大きく減衰される。そのような関数の例は、逆指数であり、それにより計算された信号強度は、
【数９】
2019513515.tif 000011
となり、ここで、
【数１０】
2019513515.tif 000012
はノイズ抑制の度合いに影響を与えるユーザ制御可能なパラメータである。
【００２８】
とりわけ（この技術のより一般的な適用性に留意しながら）ｅＤＷＩに焦点を当てると、ＤＷＩにおける信号強度は、以下の
【数１１】
2019513515.tif 000013
に従ってモデル化される。ここで、シーケンスパラメータは、必ずしも固有ではないＭ個のｂ値のベクトルβＤＷＩ＝（ｂ１，ｂ２，．．．，ｂＭ）からなる。結果の画像の対数を取ることにより、これは線形方程式：
【数１２】
2019513515.tif 000014
に変換される。
【００２９】
逆モデルに適合させるために、標準の最小二乗ルーチンを採用することができる。すなわち、
【数１３】
2019513515.tif 000015
であり、ここで、
【数１４】
2019513515.tif 000016
であり、α’＝（ＡＤＣ’，ＩｎＳ０’）である。ＡＤＣの推定における分散は、
【数１５】
2019513515.tif 000017
によって与えられることが示され得る。ここで、σ２は、各画素について計算された、線形フィットからの残差平方和を表す。分散マップに沿ったＡＤＣマップの生成が図１に示されており、（ａ）は所与の断面および所与のｂ値に対する拡散強調ＭＲ画像を示しており、（ｂ）は複数のそのような画像についてのｂ値に対するＩｎＳ（ｂ）のプロットを示しており、プロット内のＩｎＳ（ｂ）に対する値は（ａ）において示された画素でとられ、およびその画素でのＡＤＣに対する計算された値およびＡＤＣにおける標準偏差σＡＤＣを示しており、（ｃ）は他の画素に対してＡＤＣ計算を繰り返すことによって得られた断面に対するＡＤＣ画像を示しており、および（ｄ）は他の画素に対するσＡＤＣ計算を繰り返すことにより得られた断面に対するσＡＤＣマップを示している。
【００３０】
σＡＤＣ２の上記計算は、少なくとも３つのｂ値（Ｍ≧３）を取得する必要がある。しかしながら、２つのｂ値のみが利用可能である場合、単一の画像からの画像ノイズを計算するための多数の方法（例えば［１３］参照）のいずれかが、感ＭＲ性の物理的性質における不確実性を推定するために使用されてもよく、その一例が付録に記載されている。データサポートがある場合、ＡＤＣおよびσＡＤＣ２の推定を改善するために、加重最小二乗推定を使用することができる。σＡＤＣ２のノンパラメトリックな推定は、当業者に周知なブートストラップ手法によっても達成することもできる。また、前述のように、そのような手法はＤＷＩの範囲外で適用することができる。従って、それらは例えば、Ｔ１ｗ及びＴ２ｗの撮像に適用することができる。
【００３１】
ノイズ補正された指数強調ＤＷＩ（ｎｉｃｅＤＷＩ）のコントラストは、次のように定義することができる。
【数１６】
2019513515.tif 000018
ここで、ｂｃとａは、合成画像における拡散およびノイズ強調の量をそれぞれ制御するユーザ定義可能なパラメータである。逆指数関数の使用は、最終的な画像で望ましい特性を有している。特に、高いｂｃ値は、ＡＤＣが低い領域において高い信号強度を有する画像をもたらし、高いａ値は、適合されたＡＤＣ推定値の高い信頼度がある領域において高い信号強度を有する画像をもたらす。この手法は、図２に示されており、（ａ）は図１（ｃ）のＡＤＣ画像を示しており、（ｂ）は図１（ｄ）のσＡＤＣマップを示しており、（ｃ）はＡＤＣ画像から導かれたｅｘｐ｛−ｂｃ・ＡＤＣ｝のマップを示しており、（ｄ）はσＡＤＣマップから導かれたｅｘｐ｛−ａ・σＡＤＣ｝のマップを示しており、（ｅ）は強調ＭＲ画像  Ｓｎｃ（ａ，ｂｃ）＝ｅｘｐ｛−ｂｃ・ＡＤＣ｝・ｅｘｐ｛−ａ・σＡＤＣ｝を示している。
【００３２】
結果
全身のＤＷ画像は、以下のパラメータ：ｂ値＝５０／９００または５０／６００／９００ｓ／ｍｍ２、直行拡散符号化方向、ＳＴＩＲ脂肪抑制（Ｔｌ＝１８０ｍｓ）、スライス厚さ＝５ｍｍ、ＧＲＡＰＰＡ画像取得（減衰因子＝２）を有する１．５Ｔ  ＭＲＩスキャナ（ＭＡＧＮＥＴＯＭ  Ａｅｒａ，Ｓｉｅｍｅｎｓ  ＡＧ、医療部門，エアランゲン，ドイツ）で得られた。画像は、それぞれが全胴体をカバーするように４０スライスからなる４ステーションにおける移動テーブルで取得された。各ステーションでの画像は、それぞれ１信号平均（ＮｅＸ）で３回取得され、個々の拡散符号化方向からのデータが個別にエクスポートされた。従って、最終のデータセットは、各スライス位置で、各ｂ値で取得された９個の画像からなっていた。
【００３３】
前立腺癌由来の転移性骨疾患を有する３人の患者に対して、このようにして全身のＭＲ撮像を行った。
【００３４】
各患者について、以下の画像（ａ）から（ｆ）が生成された：
（ａ）基準軸断面に対するσＡＤＣマップ
（ｂ）断面に対する対応するＡＤＣ画像
（ｃ）断面に対する対応するｎｉｃｅＤＷＩ画像。画像を計算するために使用された（ａ，ｂｃ）の値が示される。
（ｄ）断面に対する対応する高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）画像；
（ｅ）病変組織およびバックグラウンド組織との間のコントラストを視覚的に最大化したｂｃ値およびａ値でのｎｉｃｅＤＷＩ画像から得られた最大強度投影（ＭＩＰ）立体；
（ｆ）従来法で取得された高ｂ値画像（ｂ＝９００ｓ／ｍｍ２）を用いて生成された対応するＭＩＰ立体。
【００３５】
図３（ａ）から（ｆ）は、第１の患者のこれらの画像を示す。ｎｉｃｅＤＷＩは、ステーション間インターフェースに存在する信号強度のスプリアス不連続性を除去することができ（（ｅ）および（ｆ）の上矢印）、病変とバックグラウンドのコントラストを改善することができる（（ｅ）および（ｆ）の下矢印）。
【００３６】
図４（ａ）から（ｆ）は、第２の患者のこれらの画像を示す。再び、ｎｉｃｅＤＷＩは、ステーション間インターフェースに存在する信号強度のスプリアス不連続性を除去することができる（（ｅ）および（ｆ）の上矢印）。また、腎臓や脾臓などの正常組織の存在を取り除くことができる（（ｅ）と（ｆ）の下矢印）。
【００３７】
図５（ａ）から（ｆ）は、第３の患者のこれらの画像を示す。再び、ｎｉｃｅＤＷＩは、腎臓および脾臓などの正常組織の存在を除去することができる（（ｅ）および（ｆ）の矢印）。
【００３８】
付録−複数取得のないＡＤＣ不確定性の推定
１．理論
先の述べたように、ＡＤＣ計算は、Ｎ個の異なるｂ値から得られる拡散強調画像Ｓｉからの線形最小二乗近似を通して達成され得る。
【数１７】
2019513515.tif 000019
【００３９】
誤差伝搬により、このようにＡＤＣ推定値の分散を
【数１８】
2019513515.tif 000020
として推定することが可能である。ここで、
【数１９】
2019513515.tif 000021
式（２）から、
【数２０】
2019513515.tif 000022
および各ｂ値での対数信号（ｌｏｇ−ｓｉｇｎａｌ）の分散の推定が必要であることは明らかである。一つの方法は、
【数２１】
2019513515.tif 000023
による評価に至る各ｂ値でのＹｉのＭｉ回の繰り返し測定を実施することである。
【００４０】
そのようなデータが不均一分散性であるので定数σｉ２を仮定していないことは注目に値する。しかしながら、直接測定が可能でない場合は、以下の近似を使用し得る。非対数ドメインにおける画像ノイズがほぼガウス分布（すなわち、近似Ｓｉ〜Ｎ（μｉ，σ））であると仮定すると、誤差伝搬により、
【数２２】
2019513515.tif 000024
を有しており、それにより
【数２３】
2019513515.tif 000025
を有する［１４］。
【００４１】
式（４）を介して生成されたσＡＤＣ２のボクセル式のマップのコントラストは、ｂ値のｂｉ、および各ｂ値での測定された信号値、Ｓｉの選択にのみ依存することは明らかである。絶対測定を決定するには、下記で議論するようにσ２の推定値が必要である。
【００４２】
２．σ２の推定値
ｍａｇｎｉｔｕｄｅ  ＤＷ画像のノイズ分散を（例えば、ウェーブレット分解またはハイパスフィルタリングを介して）推定するために、多数のオプションが利用可能である。ここでは、最も低いｂ値画像からの対数−空間データの混合ガウスモデル（ＧＭＭ）を使用することの可能性を調査し、バックグラウンドノイズをセグメント化し、元のデータセットにおけるこの領域からσ２を近似する。低ｂ値の対数−画像Ｙ１のデータは、Ｐ個のガウス分布
【数２４】
2019513515.tif 000026
からなると考慮している。
【００４３】
ＧＭＭに続いて、最も低い平均を有するガウス成分がバックグラウンドノイズのクラス、すなわち
【数２５】
2019513515.tif 000027
を表すと仮定する。バックグラウンド画素の分類は次いで、従来の最大事後推定によって簡単となり、そこからノイズ分散を計算することができる（図６に示される）。
【００４４】
３．実施例
以下のｂ値：｛５０×９；６００×９；９００×９｝で軸方向ＤＷ画像が患者から取得された。式（２）および（３）を用いてσＡＤＣの黄金律（「真の」）計算が計算され、一方で式（４）および（５）を用いて、各異なるｂ値での平均画像を含む第２の画像データセットを用いて推定値が計算された。両方の推定方法の実施例は、図７において単一の患者の４つの軸方向位置に対して実証される。σＡＤＣの真の値と推定された値との間に良好な一般的対応関係があることが明らかにされる。式（４）を使用するときに、組織内でわずかな過小推定があり、バックグラウンドノイズのわずかな過大推定がある。良好な推定が特に骨などの静的な場所で見られるが、動きにより影響を受ける領域（例えば、腎臓、腸および脾臓）は一般的に過小推定される。この過小推定は有利となることができ、肝臓では、黄金律の技術を使用するときに、転移組織と健全組織との間で極めて少ないコントラストがある。式（４）の運動に対する無感覚は、病変のコントラストを改善する。
【００４５】
図８で証明されたものは、（ｉ）従来法で得られた高ｂ値（ｂ９００）画像、（ｉｉ）各ｂ値での９回の繰り返し取得から計算されたσＡＤＣを用いて生成されたｎｉｃｅＤＷＩ画像、および（ｉｉｉ）各ｂ値での繰り返し取得が可能でない上で概説された方法論から計算されたσＡＤＣを使用して生成されたｎｉｃｅＤＷＩ画像、に対する最大強度投影（ＭＩＰｓ）の比較である。これは、異なるｂ値での繰り返し取得が可能でないところで、上で概説したような手法は、異なるｂ値での取得が可能である手法と比較して、腫瘍組織（黒い矢印）と精巣などの健全組織（白い矢印）との間に類似の画像コントラストを生成することができる。
【００４６】
以下の全ての参考文献は参照により本明細書に組み込まれる。
［１］Koh DM, Blackledge M, Padhani AR, Takahara T, Kwee TC, et al. (2012) Whole-body diffusion-weighted MRI: tips, tricks, and pitfalls. AJR Am J Roentgenol 199: 252-262. 
［２］Eiber M, Holzapfel K, Ganter C, Epple K, Metz S, et al. (2011) Whole-body MRI including diffusion-weighted imaging (DWI) for patients with recurring prostate cancer: technical feasibility and assessment of lesion conspicuity in DWI. J Magn Reson Imaging 33: 1160-1170.
［３］Padhani AR, Koh OM, Collins DJ (2011) Whole-body diffusion-weighted MR imaging in cancer: current status and research directions. Radiology 261: 700-718.
［４］van Ufford HM, Kwee TC, Beek FJ, van Leeuwen MS, Takahara T, et al. (2011) Newly diagnosed lymphoma: initial results with whole-body T1-weighted, STIR, and diffusion weighted MRI compared with 18F-FDG PET/CT. AJR Am J Roentgenol 196: 662-669.
［５］Giles SL, Messiou C, Collins DJ, Morgan VA, Simpkin CJ, et al. (2014) Whole-body diffusion-weighted MR imaging for assessment of treatment response in myeloma. Radiology 271: 785-794.
［６］Koh OM, Collins DJ (2007) Diffusion-weighted MRI in the body: applications and challenges in oncology. AJR Am J Roentgenol 188: 1622-1635.
［７］Padhani AR, Makris A, Gall P, Collins DJ, Tunariu N, et al. (2014) Therapy monitoring of skeletal metastases with whole-body diffusion MRI. J Magn Reson Imaging 39: 1049-1078.
［８ａ］Blackledge MD, Leach MO, Collins DJ, Koh OM (2011) Computed diffusion-weighted MR imaging may improve tumor detection. Radiology 261: 573-581.
［８ｂ］W02010/116124.
［９］Blackledge MD, Collins DJ, Tunariu N, Orton MR, Padhani AR, et al. (2014) Assessment of treatment response by total tumor volume and global apparent diffusion coefficient using diffusion-weighted MRI in patients with metastatic bone disease: a feasibility study. PLoS One 9: e91779.
［１０］ Provenzale JM, Engelter ST, Petrella JR, Smith JS, MacFall JR (1999) Use of MR exponential diffusion-weighted images to eradicate T2 "shine-through" effect. AJR Am J Roentgenol 172: 537-539.
［１１］Gudbjartsson H, Patz S (1995) The Rician distribution of noisy MRI data. Magn Reson Med 34: 910-914.
［１２］Gatidis S, Schmidt H, Martirosian P, Nikolaou K, Schwenzer NF (2015) Apparent diffusion coefficient-dependent voxelwise computed diffusion-weighted imaging: An approach for improving SNR and reducing T shine-through effects. J Magn Reson Imaging.
［１３］Coupe P, Manj6n, JV, Gedamu E, Arnold D, Robles M, Collins L (2010) Robust Rician noise estimation for MR images. Med. Image Analysis 14(4): 483-493.
［１４］Jelle Veraart, Jan Sijbers, Stefan Sunaert, Alexander Leemans, and Ben Jeurissen. Weighted linear least squares estimation of di_usion mri parameters: strengths, limitations, and pitfalls. Neuroimage, 81:335-46, Nov 2013.
【図１】
2019513515.tif 000028
【図２】
2019513515.tif 000029
【図３】
2019513515.tif 000030
【図４】
2019513515.tif 000031
【図５】
2019513515.tif 000032
【図６】
2019513515.tif 000033
【図７】
2019513515.tif 000034
【図８】
2019513515.tif 000035
【手続補正書】
【提出日】平成30年12月11日(2018.12.11)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
関心のある領域の磁気共鳴（ＭＲ）画像を作成する方法であって、
前記関心のある領域にわたる位置で、感ＭＲ性の物理的性質の値をマッピングするステップであって、前記マッピングされた値は前記領域の初期ＭＲ画像から取得されている、マッピングするステップと、
前記領域にわたる前記位置での前記感ＭＲ性の物理的性質の値において推定された不確実性の対応するマップを決定するステップと、
前記領域の強調ＭＲ画像を計算するステップであって、前記強調ＭＲ画像は、前記初期画像の各位置において、前記感ＭＲ性の物理的性質のそれぞれの値とそれぞれの前記推定された不確実性とを組み合わせる関数の値をマッピングし、前記関数は、より高い強調を、相対的に高く推定された不確実性を有する位置よりも、相対的に低く推定された不確実性を有する位置へ適用する、計算するステップとを含み、
前記推定された不確実性は、前記感ＭＲ性の物理的性質の前記値の分散であることを特徴とする、方法。
【請求項２】
前記関数は、逆指数強調を前記推定された不確実性に適用する、請求項１に記載の方法。
【請求項３】
前記推定された不確実性は、前記領域にわたる前記感ＭＲ性の物理的性質の値の逆モデルから決定される、請求項１または２に記載の方法。
【請求項４】
前記決定するステップにおいて、前記逆モデルは、最小二乗法を用いて前記感ＭＲ性の物理的性質の前記値に適合される、請求項３に記載の方法。
【請求項５】
前記感ＭＲ性の物理的性質が、見かけの拡散係数、Ｔ１、Ｔ２、またはプロトン密度である、請求項１から４のいずれか一項に記載の方法。
【請求項６】
前記感ＭＲ性の物理的性質が見かけの拡散係数であり、前記関数は、より高い強調を、相対的に高い見かけの拡散係数を有する位置よりも、相対的に低い見かけの拡散係数を有する位置に適用する、請求項１から５のいずれか一項に記載の方法。
【請求項７】
前記関数は、前記見かけの拡散係数に逆指数強調を適用する、請求項６に記載の方法。
【請求項８】
前記領域にわたる位置での見かけの拡散係数の値をマッピングしている前記初期ＭＲ画像は、異なるｂ値で前記領域にわたる複数の拡散強調画像を取得し、そこから前記初期ＭＲ画像を導出することによって取得される、請求項６または７に記載の方法。
【請求項９】
関心のある立体を通る、それぞれ平行なスライスである複数の関心のある領域の各々について請求項１から８のいずれか一項に記載の方法を実行して、それにより対応する強調ＭＲ画像を各スライスに対して取得するステップと、
前記強調ＭＲ画像を投影三次元画像として表示するステップとを含む、３Ｄ  ＭＲ画像データを表示する方法。
【請求項１０】
請求項１から９のいずれか一項に記載の方法を実行するように動作可能に構成されたコンピュータシステム。
【請求項１１】
請求項１０に記載の前記コンピュータシステムを有する磁気共鳴撮像装置。
【請求項１２】
請求項１から９のいずれか一項に記載の方法を実行するためのコンピュータプログラム。
【請求項１３】
請求項１２に記載のプログラムを実行するコンピュータプログラム製品。
【国際調査報告】
2019513515.tif 000036
2019513515.tif 000037
2019513515.tif 000038