(57)【要約】  本発明は、灌流イメージングに確定的でない確率的方法を適用することによって血行動態パラメータを推定するためのシステムおよびプロセスに関する。そのようなプロセスは、動脈入力関数または相補累積分布関数の推定も可能にし、したがってより一般的には、任意の関心対象量の推定も可能にする。本発明は、恣意的または望ましくない仮説によって所望の推定を強いることを強要せず、生理学的または血行動態的な性質の完全不確実情報の演繹的な導入を必要とする点で、とりわけ公知のプロセスから卓越する。
2013539704.tif 000002
【特許請求の範囲】
【請求項１】
器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量（１４）を推定（５７）するために、灌流強調イメージング解析システムの処理ユニット（４）によって実行され、実験による灌流データ（１５）から前記関心対象量を推定するステップを含んでいる方法であって、前記ステップが、
前記ボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定に関係するパラメータに鑑みた実験による灌流データの直接確率分布の割り当て（５４）を行うこと、および
前記動態系のインパルス応答に関する完全不確実情報の導入（５１ｂ）を行なうことによって、前記量の結合事前確率分布の割り当て（５３）を行なうことによって、前記関心対象量の事後周辺確率分布（５６）をベイズ法に従って評価することを含むことを特徴とする方法。
【請求項２】
器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量（１４）を推定（５７）するために、灌流強調イメージング解析システムの処理ユニット（４）によって実行され、実験による灌流データ（１５）から前記関心対象量を推定するステップを含んでおり、前記動態系は、線形であり、時不変であり、関係Ｃ（ｔ）＝ＢＦ・Ｃａ（ｔ）＊Ｒ（ｔ）によって形式的に決定され、ここでＣ（ｔ）はボクセルを循環する造影剤の濃度であり、Ｃａ（ｔ）は前記ボクセルに流れ込む動脈における前記造影剤の濃度であり、ＢＦは前記ボクセルにおける血流量であり、＊は畳み込み積を表わし、Ｒ（ｔ）は前記ボクセルにおける移動時間の相補累積分布関数である方法であって、前記ステップが、
前記ボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定に関係するパラメータに鑑みた実験による灌流データの直接確率分布の割り当て（５４）を行うこと、および
前記ボクセルにおける相補累積分布関数Ｒ（ｔ）に関する完全不確実情報の導入（５１ｂ）を行なうことによって、前記量の結合事前確率分布の割り当て（５３）を行なうことによって、前記関心対象量の事後周辺確率分布（５６）をベイズ法に従って評価することを含むことを特徴とする方法。
【請求項３】
前記量の結合事前確率分布の割り当てが、ボクセルに流れ込む動脈における前記造影剤の濃度−時間曲線に関する完全不確実情報の導入（５１ｂ）によってさらに達成されることを特徴とする請求項１または２に記載の方法。
【請求項４】
該当の複数のボクセルについての逐次の反復によって実行されることを特徴とする請求項１〜３のいずれか一項に記載の方法。
【請求項５】
推定した関心対象量に関する信頼区間によって表わされる補足の情報を計算するステップ（５８）を含むことを特徴とする請求項１〜４のいずれか一項に記載の方法。
【請求項６】
推定した関心対象量に関する信頼区間への賭け率によって表わされる補足の情報を計算するステップ（５９）を含むことを特徴とする請求項１〜５のいずれか一項に記載の方法。
【請求項７】
該当のボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定の問題に関係するパラメータに鑑みた実験による灌流データの直接確率分布の割り当て（５４）、および
前記量の結合事前確率分布の割り当て（５３）の産物の適切さの指標によって表わされる補足の情報を計算するステップ（６０）を含むことを特徴とする請求項４に記載の方法。
【請求項８】
推定した関心対象量（１４）を、該関心対象量をユーザ（６）へと提供することができるヒューマン・マシン・インターフェイス（５）へともたらすステップを含むことを特徴とする請求項４〜７のいずれか一項に記載の方法。
【請求項９】
前記推定した関心対象量に関する補足の情報を、該情報をユーザ（６）へと提供することができるヒューマン・マシン・インターフェイス（５）へともたらすステップを含むことを特徴とする請求項５〜８のいずれか一項に記載の方法。
【請求項１０】
実験による灌流データが、実験による灌流信号の値のベクトルまたは該ベクトルの濃度−時間曲線の値のベクトルへの変換を含むことを特徴とする請求項１〜９のいずれか一項に記載の方法。
【請求項１１】
記憶手段と、外界との通信手段と、処理手段とを備える処理ユニット（４）であって、
通信手段が、外界から実験による灌流データ（１５）を受け取ることができ、
処理手段が、器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量（１４）を推定するための請求項１〜１０のいずれか一項に記載の方法を実行するように構成されていることを特徴とする処理ユニット。
【請求項１２】
通信手段が、推定された関心対象量（１４）を、該関心対象量をユーザ（６）へと提供することができるヒューマン・マシン・インターフェイス（５）に適したフォーマットでもたらす（５７）ことを特徴とする請求項１１に記載の処理ユニット。
【請求項１３】
通信手段が、推定された関心対象量（１４）に関する補足の情報を、該情報をユーザ（６）へと提供することができるヒューマン・マシン・インターフェイス（５）に適したフォーマットでもたらす（５８、５９）ことを特徴とする請求項１２に記載の処理ユニット。
【請求項１４】
請求項１１〜１３のいずれか一項に記載の処理ユニット（４）と、前記処理ユニット（４）によって実行される請求項１〜１０のいずれか一項に記載の方法に従って、推定された量（１４）をユーザ（６）へと提供することができるヒューマン・マシン・インターフェイス（５）とを備える灌流強調イメージング解析システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、灌流強調イメージング（ｐｅｒｆｕｓｉｏｎ−ｗｅｉｇｈｔｅｄ  ｉｍａｇｉｎｇ）へと確定的でない確率論的手法（ｓｏｆｔ  ｐｒｏｂａｂｉｌｉｓｔｉｃ  ｍｅｔｈｏｄｓ）を適用することによって血行動態パラメータを推定するためのシステムおよび方法に関する。さらに、そのような方法は、相補累積分布関数または動脈入力関数の推定を可能にし、その後に、より一般的に、任意の関心対象量の推定を可能にする。本発明は、とりわけ、必要とされる推定を恣意的な望ましくない仮定によって制約または行使することがなく、不確実な事前の生理学的または血行動態的な情報の導入を必要とする点で、公知の方法から相違する。
【背景技術】
【０００２】
本発明は、詳細には、灌流強調磁気共鳴イメージング（ＰＷ−ＭＲＩ）または灌流コンピュータ断層撮影（ＣＴ）に依存する。それらの技術は、脳または心臓などの器官の血行動態についての有用な情報を迅速に取得することを可能にする。この情報は、急性脳梗塞などの病変の救急処置において診断および治療上の決定を行おうとする医師を助けるためにとくに重要である。
【０００３】
それらの技術は、核磁気共鳴装置またはコンピュータ断層撮影装置に依存する。この装置は、脳などの体の一部分の複数のデジタル画像系列をもたらす。この目的のため、前記装置は、高周波数の電磁波の組み合わせを体の前記部分へと照射し、特定の原子によって再放出される信号を測定する。このやり方で、この装置は、画像化された塊の各点（すなわち、ボクセル）の化学組成を明らかにし、したがって生物学的組織の種類を明らかにすることを可能にする。
【０００４】
画像系列は、専用処理ユニットによって解析される。この処理ユニットは、最終的には、適切なヒューマン・マシン・インターフェイスによって、灌流強調画像からの血行動態パラメータの推定を医師へともたらす。このやり方で、医師は、診断を行い、どの治療決定が適切であるかを判断することができる。
【０００５】
核磁気共鳴またはコンピュータ断層撮影灌流強調画像は、造影剤（例えば、磁気共鳴イメージングにおけるガドリニウムキレート）を静脈に注入し、画像の各ボクセルにおいて長時間にわたってこの造影剤のボーラスを記録することによって取得される。簡潔にするため、ボクセルを識別するための添え字ｘ，ｙ，ｚを省略する。例えば、座標ｘ，ｙ，ｚのボクセルの信号を、Ｓｘ，ｙ，ｚ（ｔ）と表わす代わりに、単にＳ（ｔ）と表わす。後述される操作および演算が、推定すべき血行動態パラメータを表わす画像またはマップを最終的に取得するために、着目される各ボクセルについて一般的に実行されることを、理解すべきである。
【０００６】
標準モデルが、或る時間ｔにわたって測定された信号強度Ｓ（ｔ）を、前記造影剤の濃度Ｃ（ｔ）に関連付けることを可能にする。
【０００７】
例えば、灌流コンピュータ断層撮影では、各ボクセルの信号が濃度に正比例し、すなわちＳ（ｔ）＝ｋ・Ｃ（ｔ）＋Ｓ０である。核磁気共鳴による灌流イメージングにおいては、指数関数の関係Ｓ（ｔ）＝Ｓ０・ｅ−ｋ・ＴＥ・Ｃ（ｔ）が存在する。どちらの場合も、Ｓ０が、造影剤の到達前の平均の信号強度を表わしている。核磁気共鳴イメージングに関しては、ｋが、組織における常磁性磁化率と造影剤の濃度との間の関係に依存する定数であり、ＴＥが、エコー時間である。各ボクセルの定数ｋの値は、未知であるため、着目されるすべてのボクセルについて任意の値に設定される。したがって、絶対的な推定ではなく、相対的な推定が得られる。しかしながら、この相対的な情報は、空間におけるそれらの値の相対的な変化（詳細には、正常な組織と病んだ組織との間での違い）が主たる関心の対象であるため、依然として有意義である。
【０００８】
一般に、理論上の信号Ｓ（ｔ）を造影剤の理論上の濃度Ｃ（ｔ）へと結び付けるモデルが、Ｓ（ｔ）＝Ψ（Ｃ（ｔ），ΘＳ）と称され、ΘＳが、このモデルの自由パラメータのベクトルである。例えば、磁気共鳴またはコンピュータ断層撮影による灌流強調イメージングにおいては、ΘＳ＝（Ｓ０）である。
【０００９】
各々の時点における各ボクセルに囲まれた組織領域の造影剤の質量の保存が、
【００１０】
【数１】
2013539704.tif 000003
【００１１】
と表わされる。Ｃａ（ｔ）が、この組織領域へと送り込まれる動脈の造影剤の濃度であり、動脈入力関数またはＡＩＦとして知られている。ＢＦが、この組織領域における血流量であり、Ｃｖ（ｔ）が、この組織領域から流れ出る静脈の造影剤の濃度であり、静脈出力関数またはＶＯＦとして知られている。
【００１２】
動脈／組織／静脈の動態系が、線形であり、時間変化しないと仮定すると、
【００１３】
【数２】
2013539704.tif 000004
【００１４】
であり、ここでｈ（ｔ）は、系のインパルス応答または組織における造影剤の移動時間の確率密度関数であり、「ｘ」を○で囲んだ符号（以下、「＊」で表すことがある）は、畳み込み積を表わす。したがって、初期条件Ｃ（ｔ＝０）＝０における上述の微分方程式の形式解は、Ｃ（ｔ）＝ＢＦ・Ｃａ（ｔ）＊Ｒ（ｔ）と表わされ、ここでＲ（ｔ）は、
【００１５】
【数３】
2013539704.tif 000005
【００１６】
と定義される相補累積分布関数または剰余関数であり、ここでＨは、ヘビサイド単位ステップ一般化関数である。インパルス応答および相補累積分布関数から、別の血行動態パラメータ、すなわち組織における平均通過時間またはＭＴＴが定められる。
【００１７】
【数４】
2013539704.tif 000006
【００１８】
また、血液容量ＢＶを関係ＢＶ＝ＢＦ・ＭＴＴによって定めることもできる。
【００１９】
核磁気共鳴による灌流強調イメージングにおいて、ＢＦ、ＭＴＴ、またはＢＶなどの血行動態パラメータ、ならびに相補累積分布関数は、現在のところ以下のように推定されている。
【００２０】
各々のボクセルについて、時点ｔｉ，ｉ＝１，Ｎにおいてサンプリングされた実験による灌流信号Ｓｅｘｐ（ｔ）が、関係
【００２１】
【数５】
2013539704.tif 000007
【００２２】
を使用することによって濃度−時間曲線Ｃ（ｔ）へと変換される。定数ｋが、すべてのボクセルについて、ゼロでない任意の値（例えば、ｋ・ＴＥ＝１）に設定される。定数Ｓ０が、例えば造影剤の到着前の平均を取ることによって推定される。これが、灌流信号の取得が造影剤のボーラス到着時間（すなわち、ＢＡＴ）と比べて充分に早期に開始される場合にかぎって可能であることに、注意すべきである。濃度Ｃ（ｔ）から、関連の理論上の動脈入力関数Ｃａ（ｔ）が既知であると仮定し、積ＢＦ・Ｒ（ｔ）が、数値逆畳み込みによって推定される。
【００２３】
濃度−時間曲線Ｃ（ｔ）の逆畳み込みのための理論上の動脈入力関数Ｃａ（ｔ）を得るために、いくつかの手法が提案されている。
【００２４】
第１の手法においては、医師が、実験による大域的な動脈入力関数を手作業で選択する。これを、脳の灌流イメージングにおいて、例えば対側性シルビウス動脈（ｃｏｎｔｒａｌａｔｅｒａｌ  ｓｙｌｖｉａｎ  ａｒｔｅｒｙ）または内頸動脈において測定することができ、あるいは例えば光学的な測定など、追加の測定から得ることができる。高い信号対雑音比の信号を得ることができるにせよ、この手法は、多くの欠点を有している。第１に、人間の介入および／または追加の測定が必要である。これは、臨床の緊急の状況において望ましくなく、手順および最終的な結果の再現を、より困難にする。第２に、とりわけ、この大域的動脈入力関数は、各ボクセルの局所的動脈入力関数に一致しない。局所的動脈入力関数からの相違が、遅延（局所的動脈入力関数は、一般に血管系の上流において取得される大域的動脈入力関数と比べて遅れるため）および分散（造影剤の伝播が、血管系の上流よりも下流において遅いため）に関して存在する。ここで、畳み込み積の対称性により、これらの欠点が相補累積分布関数の推定に直接影響するため、これらの現象が、最終的に、血行動態パラメータの推定にかなりの影響を有することが既知である。したがって、例えば、局所的動脈入力関数と局所的静脈出力関数との間の真の平均通過時間（ＭＴＴ）の推定を最終的に得ることができず、大域的動脈入力関数と静脈出力関数との間の平均通過時間しか得ることができない。この乖離を克服するために、一部の者は、大域的動脈入力関数と局所的動脈入力関数との間の遅延を定量化する
【００２５】
【数６】
2013539704.tif 000008
【００２６】
というパラメータなどの新たな記述パラメータを、たとえそれらが（動脈入力関数が各ボクセルの真の局所的動脈入力関数である）元の標準的な灌流モデルに属さないにせよ、導入している。他の方法は、局所的動脈入力関数のこれらの乖離について、血行動態パラメータの推定への影響を最小限にしようと試みている。しかしながら、それらは大域的問題に新たな未知の要素を導入し、それを回避しているにすぎない。
【００２７】
第２の手法によれば、大域的動脈入力関数が、データクラスタリングまたは独立成分分析（ＩＣＡ）などの信号処理技術によって灌流強調画像から自動的に取得される。この手法は、人間の介在を回避できるようにするかもしれないが、大域的動脈入力関数につきものの遅延および分散の問題を解決しておらず、新たな未知の要素を導入している（例えば、動脈入力関数の代わりに静脈出力関数を得ることができる）。
【００２８】
第３の手法によれば、局所的動脈入力関数が、信号処理技術および選択基準によって灌流強調画像から自動的に取得される。例えば、＜＜最良の＞＞関数が、血行動態パラメータまたは相補累積分布関数を推定すべき現在の組織ボクセルのすぐ近傍において捜索される。この第３の手法の目的は、遅延および分散の問題を少なくとも或る程度は克服することによって、より偏りが少なく、より正確な推定を、最終的に得ることである。しかしながら、このようにして得られた局所的動脈入力関数が対象のボクセルの＜＜真の局所的関数＞＞の適切な近似であることが、演繹的にも帰納的にも保証されない。例えば、この＜＜真の＞＞関数が、該当の近傍の外側に位置する（該当の近傍が小さすぎる場合）可能性があり、反対に、別の動脈入力関数と混同される（該当の近傍が大きすぎる場合）可能性がある。さらに、この＜＜最良の＞＞局所的動脈入力関数は、＜＜正常な＞＞動脈入力関数（すなわち、短い／高い早さの造影剤到着時間を有する、大きな振幅を有する、など、）の中から求められる。しかしながら、目的は、例えば虚血性の動脈入力関数など、病気の動脈入力関数から、正常な動脈入力関数を正確に区別することである。結果として、たとえ最終的な結果が大域的な手法よりも良好になりうる場合でも、局所的な動脈入力関数についての不確かさが依然として大きく、血行動態パラメータまたは相補累積分布関数についての不確かさはなおさらである。
【００２９】
上述の方法から得られるとおりの理論上の動脈入力関数Ｃａ（ｔ）によって実験による濃度−時間曲線Ｃ（ｔ）の逆畳み込みを行なうために、標準的な畳み込みモデルＣ（ｔ）＝ＢＦ・Ｃａ（ｔ）＊Ｒ（ｔ）が、最初に例えば矩形法の近似に従って時間について離散化される。
【００３０】
【数７】
2013539704.tif 000009
【００３１】
ここで、Δｔはサンプリング周期である。このやり方で、下記のように線形システムＡｄ＝ｃに帰着する。
【００３２】
【数８】
2013539704.tif 000010
【００３３】
実際には、行列Ａが、深刻に悪い状態であり、ほとんど特異であるため、無意味な解および異常な推定を得ることなく線形システムを数値的に解くことはできない。したがって、例えば行列Ａの疑似逆行列Ａ〜−１を得、次いで
【００３４】
【数９】
2013539704.tif 000011
【００３５】
によってｄの推定ｄ＾を得るために、種々の方法に頼らなければならない。逆行列を得るためのそれらの方法として、ｓＳＶＤ法（単純特異値分解）、ｃＳＶＤ法（円形特異値分解）、およびｏＳＶＤ法（振動特異値分解）などのＡの特異値の打ち切りにもとづく方法（打ち切り特異値分解または（Ｔ）ＳＶＤ）や、周波数域におけるＨｕｎｔ逆畳み込みが挙げられる。
【００３６】
より一般的には、
【００３７】
【数１０】
2013539704.tif 000012
【００３８】
などの基準を最小にすることができ、ここで
【００３９】
【数１１】
2013539704.tif 000013
【００４０】
は、特定の解を有利にし、
【００４１】
【数１２】
2013539704.tif 000014
【００４２】
によってｄの推定を得ることを可能にする正則化項である。そのような方法として、Ｔｉｋｈｏｎｏｖ正則化およびウェーブレット変換にもとづく方法などが挙げられる。
【００４３】
ひとたびｄ＾が得られると、定義により
【００４４】
【数１３】
2013539704.tif 000015
【００４５】
であるため、ＢＦ＾＝ｄ＾（ｔ１）＝ｄ＾（０）によってＢＦの推定ＢＦ＾を得ることができる。しかしながら、ＢＦは、多くの場合、例えば特異値分解にもとづく方法において、それらの方法に固有のｄ＾（０）の系統的な過小評価を補償するために、
【００４６】
【数１４】
2013539704.tif 000016
【００４７】
として推定される。次いで、
【００４８】
【数１５】
2013539704.tif 000017
【００４９】
によってｂの推定ｂ＾が得られ、例えば矩形法の近似に従うことにより、
【００５０】
【数１６】
2013539704.tif 000018
【００５１】
によって
【００５２】
【数１７】
2013539704.tif 000019
【００５３】
の推定が得られる。最後に、典型的には、たとえ積の推定量が推定量の積でなくても、
【００５４】
【数１８】
2013539704.tif 000020
【００５５】
によってＢＶの推定が得られる。
【００５６】
これらの動脈入力関数にもとづく方法について、多数の変種を見つけることができ、例えば信号対雑音比を人工的に高めるために、実験による動脈入力関数を、最初にパラメトリックまたはセミパラメトリック理論モデルＣａ（ｔ，Θａ）（ここで、Θａはパラメータのベクトルである）へと適合させることができる。同様に、数値逆畳み込みをより安定にし、あるいは造影剤の循環信号（第１の通過）と再循環信号（第２の通過）との間に時間的な重なり合いが存在する場合に、再循環に起因して生じうる問題を克服するために、信号を人工的にオーバーサンプリングすることができる。しかしながら、それらの方法も、依然として打ち切り特異値分解などの数値逆畳み込み法にもとづいている。
【００５７】
典型的な実データを模擬するように意図された合成データについて試験した種々の方法についてのいくつかの比較研究によれば、特異値分解にもとづく逆畳み込みおよびその変種（線形打ち切り特異値分解（ｓＳＶＤ）、打ち切り円形特異値分解（ｃＳＶＤ）、または平滑化打ち切り円形特異値分解（ｏＳＶＤ））が、実務において患者、組織の種類、病状、などに応じて生じうる種々の相補累積分布関数Ｒ（ｔ）および動脈入力関数Ｃａ（ｔ）に対して、偏り（真の値に対する系統誤差）、精度（入力信号の信号対雑音比の関数としての真の値に対する推定の標準偏差）、およびロバストさに関してＢＦおよびＭＴＴパラメータの最良の推定を最終的にもたらす。
【００５８】
しかしながら、上述した動脈入力関数の選択に関する問題に加えて、この種の数値法は、激しい固有の問題を抱えている。
【００５９】
第一に、ＢＦ・ｂの推定ｄ＾が、時間とともに減少するのではなく、ときには負の値をとりうるような程度にまで振動する。しかしながら、時刻ｔにおいてボクセルに残る造影剤の量Ｒ（ｔ）は、当然ながら正であり、減少関数である。ｏＳＶＤなどの特別な方法は、それらの異常な振動を軽減することができるが、それらは特異値分解に固有のものであるがゆえに依然として残る。これが、この種の方法において、血流量ＢＦが
【００６０】
【数１９】
2013539704.tif 000021
【００６１】
によって推定される一方で、標準的な灌流モデルにおいては
【００６２】
【数２０】
2013539704.tif 000022
【００６３】
であるがゆえに、
【００６４】
【数２１】
2013539704.tif 000023
【００６５】
として推定されるべき理由である。始点における値の代わりに最大値をとることによって、それらの振動のＢＦの推定に対する影響だけを消去することが期待される。したがって、これらの方法は、完璧に満足できるものにはなり得ない。詳細には、相補累積分布関数、したがって血行動態パラメータのより厳密な推定方法によって到達できるＢＦの推定の精度は、未知のままである。したがって、これらの数値逆畳み込み法は、モデルの特性を満たさない解をもたらすため、標準的な灌流モデルと矛盾する。
【００６６】
この問題を克服し、相補累積分布関数の生理学的に容認できる推定を得るために、それらの相補累積分布関数のパラメトリックモデルＲ（ｔ，ΘＲ）が導入されており、ここでΘＲは、このモデルのパラメータのベクトルである。これらのモデルが、例えばベイズ法を使用することによって実験による信号に合わせられる。しかしながら、この手法は、現時点においては時期尚早であるように見受けられる。実際、相補累積分布関数を記述するために適したパラメータまたはセミパラメトリックモデルを決定するために、相補累積分布関数の事前の非パラメータ的推定が思いのままでなければならない。なぜならば、それらのモデルが、実務において生じうるあらゆる種類の相補累積分布関数を正しく表わすことに完璧に適していないと、得られるＭＴＴまたはＢＦなどの血行動態パラメータの推定が異常になることが、モンテカルロシミュレーションによって示されているからである。したがって、相補累積分布関数の理論上のモデルの選択が重要であり、モデルを実験データへと適用することによってのみ正しく行なうことが可能である。したがって、前記相補累積分布関数を推定するための非パラメータ法が、おそらくは次の段階における古典的なパラメータまたはセミパラメトリックモデルによる置き換えによって相補累積分布関数の生理学的に容認できる推定を得ることを可能にするために、必要である。
【００６７】
しかしながら、上述のように、特異値分解にもとづく方法などの方法によって得られる推定は、標準的な灌流モデルのもとでの
【００６８】
【数２２】
2013539704.tif 000024
【００６９】
の定義そのものに矛盾する。それらは、生理学的にも、物理学的にも、容認することができない。したがって、パラメータまたはセミパラメトリック理論モデルをそれらの推定へと当てはめることが不可能である。さらに、いくつかのモデルを互いに比較し、実験による相補累積分布関数を表わすために最も適したモデルを選択することも、不可能である。したがって、特異値分解法などの数値法につきものの欠陥ゆえに、灌流現象のモデル化および理解における進歩はもはや不可能である。
【００７０】
さらに、相補累積分布関数および血行動態パラメータの非パラメータ的推定の根底にある問題は、実験に基づく動脈入力関数による、実験に基づく濃度−時間曲線の＜＜単純な＞＞逆畳み込みの問題ではないように思われる。実際、絶対的な確かさおよび無限の精度で知られた動脈入力関数が問題において実際に与えられるならば、そのようであるかもしれない。しかしながら、最良でも、最大でも測定雑音までしか知ることができず、前もって濃度−時間曲線へと変換しなければならない実験的に測定された動脈信号が与えられるにすぎない。換言すると、実験的に測定された動脈入力関数が理論上の動脈入力関数に等しいと仮定することによって、実験による信号における測定雑音ならびにそれらの信号からの濃度−時間曲線の推定または変換に起因する不確かさは無視される。
【００７１】
標準的な畳み込みモデルＣ（ｔ）＝ＢＦ・Ｃａ（ｔ）＊Ｒ（ｔ）は、直接的に測定することが一般に不可能である理論上の信号だけを含む。実際、測定される信号は、典型的には、理論上の信号と測定雑音との合計である。したがって、実験による組織および動脈の灌流信号は、それぞれＳｅｘｐ（ｔ）＝Ｓｔｈ（ｔ）＋ξｔおよびＳａｅｘｐ（ｔ）＝Ｓａｔｈ（ｔ）＋ξｔａと表わされ、ここでξｔおよびξｔａは、測定雑音をモデル化する平均ゼロの確率過程である。さらに
【００７２】
【数２３】
2013539704.tif 000025
【００７３】
および
【００７４】
【数２４】
2013539704.tif 000026
【００７５】
であり、ここでＣｔｈ（ｔ）は、磁気共鳴灌流強調イメージングにおける理論上の濃度−時間曲線であり、Ｃａｔｈ（ｔ）は、理論上の動脈入力関数であり、あるいは、より一般的には、上述のようにＳｔｈ（ｔ）＝Ψ（Ｃｔｈ（ｔ），ΘＳ）である。したがって、実験による信号へと適用される理論上の標準的な灌流モデルは、Ｃｔｈ（ｔ）＝ＢＦ・Ｃａｔｈ（ｔ）＊Ｒ（ｔ）と表わされなければならず、すなわち核磁気共鳴灌流強調イメージングにおいては以下のように表わされなければならず、
【００７６】
【数２５】
2013539704.tif 000027
【００７７】
あるいは同等には、下記のように表わされなければならず、
【００７８】
【数２６】
2013539704.tif 000028
【００７９】
下記のようにではない。
【００８０】
【数２７】
2013539704.tif 000029
【００８１】
ここで、Ｃｅｘｐ（ｔ）は、実験による濃度−時間曲線であり、Ｃａｅｘｐ（ｔ）は、実験による動脈入力関数であり、逆畳み込み法の大部分は、この誤りのある間接的な数学的畳み込みモデルにもとづいている。この時点で、必ずしも正でないランダム変数の対数をとることを避けるために、標準的な灌流モデルを
【００８２】
【数２８】
2013539704.tif 000030
【００８３】
と表わすことが好ましいことを、理解できるであろう。
【００８４】
逆畳み込みされる信号の測定の不確かさが、逆畳み込みプロセスの最終的な結果に大きく影響することが知られており、それらの不確かさに起因する入力信号の非常に小さい変動が、最終的な結果のかなりの変動につながりかねない。まさにこれらの問題を克服し、このような不安定さを軽減するために、Ｔｉｋｈｏｎｏｖ正則化などの逆畳み込み法または特異値分解にもとづく方法が導入されている。さらに、現在のところ完全に無視されている測定の雑音および不確かさの動脈入力関数への影響が、さらにより顕著であり、動脈の測定雑音ξｔａおよびＳ０ａへの不確かさが、今や畳み込み行列Ａにおいて生じ、結果として増幅され、伝播する。動脈入力関数への誤差および不確かさを無視することは、血行動態パラメータの推定に重大な誤差を引き起こし、それらの推定に精度の錯覚も引き起こす。いくつかの方法は、現在まで避けられたままであるこの問題を最小限にするために、動脈入力関数における測定雑音の消去を目的としている。推定の誤差を支配し、定量化できるよう、動脈入力関数の不確かさを血行動態パラメータおよび相補累積分布関数の推定に伝播させることができる方法を思いのままにできることが、好ましいと考えられる。
【００８５】
さらに、標準的な灌流モデルをＣｔｈ（ｔ）＝ＢＦ・Ｃａｔｈ（ｔ）＊Ｒ（ｔ）と表わすことを、避けることが好ましいと考えられる。実際、標準的な灌流モデルが、最初に動脈／組織／静脈の動態系のインパルス応答ｈ（ｔ）を定め、そこから相補累積分布関数Ｒ（ｔ）が、
【００８６】
【数２９】
2013539704.tif 000031
【００８７】
と計算される。したがって、標準的な灌流モデルをインパルス応答ｈ（ｔ）の関数として
【００８８】
【数３０】
2013539704.tif 000032
【００８９】
と表わすことが、例えば、矩形法による近似
【００９０】
【数３１】
2013539704.tif 000033
【００９１】
のように相補累積分布関数の推定Ｒ＾を最終的に得るために、ｈ＾によって測定時点ｔｊ，ｊ＝１，Ｎにおけるインパルス応答ｈ（ｔ）を推定すべくこのモデルを当てはめるために、好都合である。詳細には、インパルス応答ｈ（ｔ）を相補累積分布関数の推定Ｒ＾から推定しない方がよい。数値的な観点から、最初に導関数（ｈ（ｔ））を推定し、第２の段階において不定積分（Ｒ（ｔ））を推定することが好ましく、逆ではないことを、容易に理解できるであろう。しかしながら、相補累積分布関数Ｒ（ｔ）がインパルス応答ｈ（ｔ）の代わりに直接推定されることを、確認することができる。
【００９２】
さらに、逆畳み込みの問題は、不良設定であり、無限に多くの可能解の余地がある。さらに、血行動態パラメータ、インパルス応答、相補累積分布関数、または動脈入力関数の推定に関して灌流強調イメージングにおいて直面される不良設定問題は、測定雑音の伝播の問題および実験による信号の濃度−時間曲線への変換の問題によって３重にされる逆畳み込みの問題である。
【００９３】
おそらくは１つの解だけしか認めない良設定問題へと帰着させるために、先験的に可能な無限の解の中から求められる解への事前の情報、制約を加える必要がある。これが、特定の種類の事前情報を程度の差こそあるが明示的に、あるいは程度の差こそあるが直接的に注入する多数の逆畳み込みおよび推定の方法を、技術水準において見つけることができる理由である。
【００９４】
例として、上述のとおりの古典的なＴｉｋｈｏｎｏｖ正則化を挙げることが可能である。その最も一般的なバージョンにおいては、行列ΓがΓ＝αＩＮに等しい。これは、大きなユークリッドノルムの解を不利にし、スカラーαが、この制約におけるこの不利の重みの大きさを表わしている。
【００９５】
他方で、灌流強調イメージングにおいて古典的である打ち切り特異値分解にもとづく方法は、正則化パラメータの役割を果たす所与のしきい値を下回る畳み込み行列の特異値をキャンセルすることからなる。それらの小さな特異値は、信号の高周波成分に関係しており、したがって特異値の打ち切りは、低域通過フィルタとして機能する。しかしながら、他方では、相補累積分布関数は、ｔ＝０での不連続性ゆえに、無限の帯域幅の信号である（相補累積分布関数が
【００９６】
【数３２】
2013539704.tif 000034
【００９７】
について定義されていることを思い出されたい）。他方で、小さな特異値の打ち切りは、追加の生理学的情報に直接対応するだけでなく、とりわけ追加の代数的情報、すなわち所与の線形部分空間の帰属関係に対応する。高すぎる周波数の振動が残ることが確認される。詳細には、動脈入力関数による相補累積分布関数の逆畳み込みによって推定された理論上の信号は、測定雑音に追従する傾向を有する（過剰適合）。さらに、打ち切りのしきい値が、明示的な生理学的制約の重みを意味せず（解を何らかの線形部分空間に属するように強いることを意味するにすぎない）、その効果が複雑で直接的でなく、さまざまな信号の種類に応じて最も適切な値を決定するための一般的基準を与えることが難しい。反対に、基準に対する特異値の切り捨ての作用が、きわめて間接的であるため、例えばｏＳＶＤ法における相補累積分布関数の粗さなど、そのような基準の適切な最適化の可能性が保証されない。
【００９８】
さらに、現在の血行動態パラメータまたは相補累積分布関数の推定方法は、それらの推定の精度（すなわち、推定量の標準偏差）の評価を可能にしておらず、それらの推定およびそれらの推定の精度の評価における信頼度の評価も可能にしていない。詳細には、逆畳み込みおよび推定の方法による標準的な灌流モデルの適合度の定量化が困難である。したがって、打ち切り特異値分解にもとづく方法が、実験による信号の過剰適合の傾向を有することに、注目することができる。推定による信号が滑らかでなく、むしろ反対に、測定雑音に従う傾向を有する。結果として、古典的なχ２の統計値（すなわち、二乗誤差の合計）などのモデルの適合度の指標に関して、低い値が得られる可能性がある。そのような低い値は、高品質の適合を示しているが、過剰適合ゆえに事実ではない。したがって、それらの統計値は、この場合には不適当であって誤解につながりかねず、過剰適合の検出および考慮を望む場合には避けるべきである。したがって、技術水準による逆畳み込みおよび血行動態パラメータの推定の方法に適合度の＜＜良好な＞＞指標を導入することは、自明ではない。したがって、そのような適合度の指標が一意に定められ、過剰適合を考慮に入れ、標準的な灌流モデルの適合度だけをデータへと定量化することができる方法を思いのままにすることが望まれる。
【００９９】
要約すると、血行動態パラメータ、インパルス応答、または相補累積分布関数を推定するための現在の方法は、多数の方法論的な欠陥および制御できない多数の近似を免れない。したがって、結果の解釈が難しくなり、灌流強調イメージングが充分に活用されていない。
【発明の概要】
【０１００】
本発明の目的は、公知の方法の使用から生じるすべての欠点に回答することにある。本発明の主たる目的は、必ずしも満足されず、とりわけ実験的に検証することが不可能であり、あるいは生理学的または血行動態的な種類ではないその場しのぎの制約を導入することなく、生理学的または血行動態的な種類の特定の制約を満たす解の中からインパルス応答ｈの推定ｈ＾を探し出すことを可能にする新規な方法を提供することにある。さらに、それらの方法は、特別な方法に頼ることを必要とせずに、そのような制約の重みを自動的かつ独自のやり方で決定することを可能にする。
【０１０１】
したがって、本発明は、血行動態パラメータ、相補累積分布関数、または動脈入力関数の推定の問題が最終的に良設定であることを主張することを可能にする。本発明は、その独特な解決策（おそらくは複数）を提供する。
【０１０２】
本発明によってもたらされる主要な利点の中でも、以下が可能であることを挙げることができる（ただし、これらに限られるわけではない）。
【０１０３】
・血行動態パラメータ、インパルス応答、相補累積分布関数、または動脈入力関数についての定性的または半定性的で不確実な生理学的情報を定量的なやり方で変換することができる。
【０１０４】
・動脈入力関数および実験による信号における不確かさおよび誤差を明示的に考慮に入れ、関心対象量の推定に伝播させることができる。
【０１０５】
・それらのパラメータ、インパルス応答、相補累積分布関数、または動脈入力関数の推定を、偏り（系統誤差）、精度（統計誤差）、線形性、および実務において生じうる種々の状況に対するロバストさに関して、改善することができる。
【０１０６】
・血行動態パラメータの推定、インパルス応答、相補累積分布関数、または動脈入力関数について、それらの推定において有することができる信頼を改善および明確化するために、信頼区間を得ることができ、それらの信頼区間についての賭け率さえ得ることができる。
【０１０７】
・インパルス応答、相補累積分布関数、または動脈入力関数について、生理学的および血行動態的な観点から容認することができ、標準的な灌流モデルに準拠している非パラメータ的な推定を得ることができ、したがって次の段階において、これらのインパルス応答、相補累積分布関数、または動脈入力関数について、パラメトリックまたはセミパラメトリックな理論上のモデルの適合、比較、および最終的には選択を行なうことができる。
【０１０８】
・動脈入力関数が対象の各ボクセルの組織に関して実際にどの程度まで容認できる動脈入力関数であるかを定量化することができる。
【０１０９】
・大域的灌流モデルの適合度の客観的かつ定量的な指標をもたらすことができ、最終的には比較および最も適した大域的灌流モデルの選択を可能にすることができる。
【０１１０】
この目的のため、本発明は、器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量を推定するための方法を具現化する。そのような方法は、灌流強調イメージング解析システムの処理ユニットによって実行されるように意図され、実験による灌流データから前記関心対象量を推定するステップを含んでいる。必ずしも満たされず、とりわけ実験的に検証することが不可能であり、あるいは生理学的または血行動態的な種類のものではない特別な制約を導入することなく、特定の生理学的または血行動態的な制約を満たす解の中から前記関心対象量を探し求めるために、本発明による推定のステップは、
・該当のボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定に関係するパラメータに鑑みた灌流データの直接確率分布を割り当て、
・前記動態系のインパルス応答に関する完全不確実情報を導入することによって、前記量の結合事前確率分布を割り当てることによって前記関心対象量の事後周辺分布をベイズ法に従って評価することを含む。
【０１１１】
さらに本発明は、器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量を推定するための方法であって、前記動態系が、線形であり、時不変であり、関係Ｃ（ｔ）＝ＢＦ・Ｃａ（ｔ）＊Ｒ（ｔ）によって形式的に決定され、ここでＣ（ｔ）はボクセルを循環する造影剤の濃度であり、Ｃａ（ｔ）は前記ボクセルに流れ込む動脈における前記造影剤の濃度であり、ＢＦは前記ボクセルにおける血流であり、＊は畳み込み積を表わし、Ｒ（ｔ）は前記ボクセルにおける移動時間の相補累積分布関数である、方法を具現化する。そのような方法は、灌流強調イメージング解析システムの処理ユニットによって実行されるように意図され、実験による灌流データから前記関心対象量を推定するステップを含んでいる。本発明によれば、上述したように、前記推定のステップが、
・該当のボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定に関係するパラメータに鑑みた灌流データの直接確率分布を割り当て、
・前記ボクセルにおける相補累積分布関数Ｒ（ｔ）に関する完全不確実情報を導入することによって、前記量の結合事前確率分布を割り当てることによって前記関心対象量の事後周辺分布をベイズ法に従って評価することを含む。
【０１１２】
好ましい実施の形態によれば、両方の場合において、本発明は、ボクセルに流れ込む動脈における前記造影剤の濃度−時間曲線に関する完全不確実情報の導入による前記量の結合事前確率分布の割り当てをさらに企画する。
【０１１３】
好都合には、本発明による方法を、該当の複数のボクセルについての逐次の反復によって実行することができる。
【０１１４】
関心対象量の推定の精度を評価するために、本発明による方法は、推定した関心対象量に関する信頼区間によって表わされる補足の情報を計算するステップを含むことができる。
【０１１５】
好ましい実施の形態によれば、そのような方法が、推定した関心対象量に関する信頼区間への賭け率によって表わされる補足の情報を計算するステップをさらに含むことができる。
【０１１６】
さらに、本発明による方法は、
・該当のボクセルにおける動脈／組織／静脈の動態系の関心対象量の推定の問題に関係するパラメータに鑑みた実験による灌流データの直接確率分布の割り当て、および
・前記量の結合事前確率分布の割り当ての産物の適切さの指標によって表される補足の情報を計算するステップを含むことができる。
【０１１７】
推定を提供するために、本発明による方法は、推定した関心対象量を、該関心対象量をユーザへと提供することができるヒューマン・マシン・インターフェイスへともたらすステップを含むことができる。
【０１１８】
補足の情報が計算される場合には、そのようなプロセスが、前記関心対象量に関する補足の情報を、該情報をユーザへと提供することができるヒューマン・マシン・インターフェイスへともたらすステップをさらに含むことができる。
【０１１９】
好ましい実施の態様によれば、本発明において、実験による灌流データが、実験による灌流信号の値のベクトルまたは該ベクトルの濃度−時間曲線の値のベクトルへの変換を含む。
【０１２０】
さらに本発明は、第２の目的に従い、記憶手段と、外界との通信手段と、処理手段とを備える処理ユニットに関する。通信手段が、外界から実験による灌流データを受け取ることができる。処理手段が、器官のボクセルと呼ばれる単位容積において複数の動脈／組織／静脈の動態系の中から関心対象量を推定する本発明による方法を実行するように構成されている。
【０１２１】
本発明による方法によって生成された推定をユーザに提供するために、そのような処理ユニットの通信手段は、推定された関心対象量を、該関心対象量をユーザへと提供することができるヒューマン・マシン・インターフェイスに適したフォーマットでもたらすことができる。
【０１２２】
この実施の態様によれば、さらに通信手段が、推定された関心対象量に関する補足の情報を、該情報をユーザへと提供することができるヒューマン・マシン・インターフェイスに適したフォーマットでもたらすことができる。
【０１２３】
第３の目的によれば、さらに本発明は、本発明による処理ユニットと、前記処理ユニットによって実行される本発明による方法に従って推定された量をユーザへと提供することができるヒューマン・マシン・インターフェイスとを備える灌流強調イメージング解析システムに関する。
【０１２４】
他の特徴および利点は、以下の説明を検討し、添付の図面を点検することによって、さらに明らかになる。
【図面の簡単な説明】
【０１２５】
【図１】図１は、灌流強調画像を解析するためのシステムの実施の形態を示している。
【図２】図２は、灌流強調画像を解析するためのシステムの実施の形態を示している。
【図３】図３は、造影剤の注入前における、核磁気共鳴イメージング装置によって得られた人間の脳のスライスの灌流画像を示している。
【図４】図４は、脳組織内での造影剤の循環の最中における、核磁気共鳴イメージング装置によって得られた人間の脳のスライスの灌流画像を示している。
【図５】図５は、人間の脳のボクセルに関する核磁気共鳴による灌流強調信号Ｓ（ｔ）を示している。
【図６】図６は、人間の脳のボクセルを循環する造影剤の典型的な濃度−時間曲線Ｃ（ｔ）を示している。
【図７】図７は、典型的な動脈入力関数Ｃａ（ｔ）を示している。
【図８】図８は、本発明による方法を示している。
【図９】図９は、本発明に従って推定された脳血液量のマップを示している。
【図１０】図１０は、本発明に従って推定された脳虚血の場合における脳血流のマップを示している
【図１１】図１１は、本発明に従って推定された平均通過時間のマップを示している。
【図１２】図１２は、脳血流が信頼区間に属する確率のマップを示している。
【発明を実施するための形態】
【０１２６】
図１が、灌流強調画像の解析のためのシステムを示している。核磁気共鳴イメージング装置またはコンピュータ断層撮影装置１が、コンソール２によって制御される。ユーザが、装置１を制御するためのパラメータ１１を選択することができる。装置１によって生成される情報１０から、人間または動物の体の一部分の複数の一連のデジタル画像１２が得られる。好ましい例として、人間の脳の観察からのデジタル画像を使用して、先行技術からの解決策および本発明を例証する。他の器官も考えられる。
【０１２７】
一連の画像１２を、随意によりサーバ３に保存し、患者の医療記録１３を構成してもよい。そのような記録１３は、灌流強調または拡散強調画像など、さまざまな種類の画像を含むことができる。一連の画像１２を、専用の処理ユニット４を使用して解析することができる。そのような処理ユニットは、画像を集めるべく外界と通信するための手段を備える。さらに、この通信のための手段は、処理ユニットが専用のヒューマン・マシン・インターフェイス５によって医師６または研究者に灌流強調画像１２からの血行動態パラメータ１４の推定を最終的にもたらすことを可能にする。分析システムのユーザ６は、診断の承認または却下、自身が適切であると考える治療行為の決定、さらなる検査の実行、などを行なうことができる。随意により、ユーザは、設定１６を通じて処理ユニット４の動作を設定することができる。例えば、表示のしきい値を定めることができ、あるいは閲覧を希望する推定パラメータを選択することができる。
【０１２８】
図２が、前処理ユニット７が一連の画像１２を解析して各ボクセルの灌流データ１５を取り出す解析システムの実施の形態を示している。このようにして、血行動態パラメータ１４の推定を担当する処理ユニット４が、この動作から解放され、外界との通信のための手段によって受信される灌流データ１５からの推定方法を実行する。
【０１２９】
図３が、人間の脳の厚さ５ｍｍのスライスの画像１２の典型例を示している。この画像は、核磁気共鳴によって得られている。この技術を使用して、各々のスライスについて、１．５×１．５×５ｍｍの寸法を有する１２８×１２８のボクセルの行列を得ることができる。双一次補間を使用して、画像２０などの４５８×４５８ピクセルのフラット画像を生成することができる。
【０１３０】
図４が、図３に関して示した画像と同様の画像２０を示している。しかしながら、この画像は、造影剤の注入後に得られている。この画像は、灌流強調脳画像の典型例である。図３に示した同じ画像と対照的に、動脈がはっきりと現れている。公知の技術に従い、血行動態パラメータを推定するために、異常のある半球の対側にある半球内の１つ以上の動脈入力関数２１を選択することが可能である。
【０１３１】
図５ｂが、図２に関連して説明した前処理ユニット７によってもたらされるデータ１５などの核磁気共鳴による灌流強調信号Ｓ（ｔ）の例を示している。したがって、灌流強調信号は、造影剤の注入後のボクセルの時間ｔにつれての変化を表わしている。例えば、図５ｂは、そのような信号を５０秒の期間にわたって示している。縦座標は、信号強度を任意単位にて示している。そのような信号を得るために、図１による処理ユニット４（あるいは、図２による前処理ユニット７）が、例えば図５ａに示されるような時点ｔ１，ｔ２，…，ｔｉ，…，ｔｎにおける核磁気共鳴による一連のｎ枚の灌流強調画像Ｉ１，Ｉ２，…，Ｉｉ，…，Ｉｎを解析する。所与のボクセル（例えば、ボクセルＶ）について、造影剤の注入後のボクセルの時間ｔにつれての変化を表わす灌流強調信号Ｓ（ｔ）が決定される。
【０１３２】
図６が、図５ｂに示したとおりの灌流強調信号から導出された濃度−時間曲線を示している。すでに述べたように、灌流強調信号と関連の濃度−時間曲線との間には、関係が存在する。すなわち、核磁気共鳴による灌流強調イメージングにおいては、指数関数の関係Ｓ（ｔ）＝Ｓ０・ｅ−ｋ・ＴＥ・Ｃ（ｔ）が存在し、ここでＳ０は、造影剤の到着前の平均の信号強度であり、ＴＥは、エコー時間であり、ｋは、常磁性磁化率と組織における造影剤の濃度との間の関係に依存する定数である。
【０１３３】
図６が、ボクセル内の造影剤の濃度の時間変化の観察を可能にする。ボクセルにおける造影剤の１回目の通過において高い振幅のピークが存在し、その後に造影剤の再循環（２回目の通過）の現象に関係したより低い振幅のピークが存在することを、見て取ることができる。
【０１３４】
図７が、図４に関連して示したボクセル２１などの動脈ボクセルにおける造影剤の循環を表わす典型的な動脈入力関数Ｃａ（ｔ）を示している。図７は、詳細には、造影剤の１回目の通過の後の再循環の現象がきわめて弱いことを示している。
【０１３５】
図８が、器官内のボクセルの複数の動脈／組織／静脈の動態系において関心対象量を推定するための本発明による方法の例を示している。そのような方法を、図１および２に関連して説明した相応に構成された灌流強調イメージング解析システムの処理ユニットによって実行することができる。
【０１３６】
本発明による方法は、主として、血行動態パラメータ、サンプリング時点における理論上のインパルス応答の値、または相補累積分布関数の値などの推定すべき種々の関心対象量について、１つ以上の事後周辺分布を割り当てるステップ５６を含む。さらに、前記推定を計算するステップ５７も含む。
【０１３７】
そのような事後周辺分布を割り当てるために、処理ユニットの設定プロセス５０を行なうことが必要である。処理ユニット自身が、好ましくは、１つ以上の設定セッティングによって、この設定を実行することができる。設定は、１つ以上の事後周辺分布からなるライブラリの構築ももたらすことができ、このライブラリは、あらかじめ確立されて前記ユニットのプログラムメモリに保存される。本発明によれば、前記ライブラリを、外部の処理ユニットによって使用または提供されるときに充実したものにすることができる。この外部の処理ユニットは、前記設定セッティングから前記設定を実行することができ、前記ライブラリを出力するために処理ユニットと協働することができる。
【０１３８】
したがって、本発明による方法は、割り当て５６の前に実行される設定ステップを含むことができ、その中でも、以下が必要かつ充分である。
【０１３９】
・対象ボクセルの動脈／組織／静脈の動態系の関心対象量の推定の問題に関係するすべてのパラメータを鑑みた実験データへの直接確率分布の割り当て５４。
【０１４０】
・それらすべてのパラメータについての結合事前確率分布の割り当て５３。
【０１４１】
本発明は、さまざまな応用の事例において１つ以上の関心対象量の推定を可能にする。
【０１４２】
・理論上の動脈入力関数が、絶対的な確かさおよび無限の精度で知られると仮定される。
【０１４３】
・局所的動脈入力関数は、やはり推定されるべき未知の遅延にて大域的動脈入力関数から相違する。
【０１４４】
・動脈入力信号だけが、おそらくは遅延まで測定される。
【０１４５】
・動脈入力関数が与えられず、動脈入力信号が測定されない。これが最も現実的な事例である。
【０１４６】
設定ステップは、該当の応用の事例に依存することができる。
【０１４７】
図８が、動脈入力信号がおそらくは遅延τまで測定されるときの第１の応用の例における本発明による方法を示している。
【０１４８】
実験の灌流モデルＭを、以下のように表わすことができる。
【０１４９】
【数３３】
2013539704.tif 000035
【０１５０】
ｂ＝［Ｒ（ｔ１），…，Ｒ（ｔＮ）］Ｔは、未知の相補累積分布関数の値のベクトルである。すでに説明したように、このベクトルを、インパルス応答ｈ＝［ｈ（ｔ１）＝０，…，ｈ（ｔＮ）］Ｔの値のベクトルから表現することができる。例えば、左矩形近似法によれば、
【０１５１】
【数３４】
2013539704.tif 000036
【０１５２】
であり、右矩形近似法によれば、
【０１５３】
【数３５】
2013539704.tif 000037
【０１５４】
であり、
ｃ＝［ｃ（ｔ１），…，ｃ（ｔＮ）］Ｔは、ボクセルにおける造影剤の未知の理論上の濃度の値のベクトルであり、
ｓ＝［Ｓ（ｔ１），…，Ｓ（ｔＮ）］Ｔは、灌流強調イメージングによる信号強度の実験による測定値の実または複素ベクトルであり、
ΘＳは、Ｓ（ｔ）をＣ（ｔ）に結び付けるパラメータのベクトルであり、
ξ＝［ξ（ｔ１），…，ξ（ｔＮ）］Ｔは、前記測定信号における測定雑音のベクトルであり、
ｓａ＝［Ｓａ（ｔ１），…，Ｓａ（ｔＮ）］Ｔは、実験による動脈入力信号Ｓａ（ｔ）の測定値のベクトルであり、
ΘＳａは、Ｓａ（ｔ）をＣａ（ｔ）に結び付けるパラメータのベクトルであり、
ξａ＝［ξａ（ｔ１），…，ξａ（ｔＮ）］Ｔは、前記動脈信号における測定雑音のベクトルであり、
ａ＝［Ｃａ（ｔ１），…，Ｃａ（ｔＮ）］Ｔは、未知の理論上の動脈入力関数Ｃａ（ｔ）の値のベクトルであり、
Ａは、矩形近似法、台形規則、あるいはＳｉｍｐｓｏｎ法、Ｂｏｏｌｅ法、Ｇａｕｓｓ−Ｌｅｇｅｎｄｒｅ法などのより高次の方法などによって畳み込み積分を数値的に近似することによって得られる未知の理論上の動脈入力関数ａの値のベクトルに関する畳み込み行列である。Ａは、ｃＳＶＤおよびｏＳＶＤなどのような特異値の打ち切りにもとづく逆畳み込み法において使用されるような循環畳み込み行列であってもよい。
【０１５５】
モデルＭの血行動態パラメータのベクトルをΘと称し、例えばここではΘ＝（ＢＦ，τ）である。
【０１５６】
図１、２、および８に関連し、処理ユニット４によって実行される方法は、実験による信号ｓの測定値のベクトルについての情報ＩＳおよび実験による動脈入力信号ｓａの測定値のベクトルについての情報ＩＳａをそれぞれ導入することからなる２つの初期設定ステップ５１ｄおよび５１ｅを含むことができる。
【０１５７】
例えば、実数値の測定雑音ベクトルの組の最初の２つのモーメント
【０１５８】
【数３６】
2013539704.tif 000038
【０１５９】
および
【０１６０】
【数３７】
2013539704.tif 000039
【０１６１】
に専念すると、最大エントロピー（Ｌｅｂｅｓｇｕｅ基準測度のもとでの微分Ｓｈａｎｎｏｎエントロピー）の原理により、ベクトルξおよびξａが、標準偏差σＳおよびσＳａをそれぞれ有する相互に独立な白色の定常ガウス確率過程と見なされる必要がある。
【０１６２】
より一般的には、測定雑音ベクトルの組（ξ，ξａ）を特徴付けるパラメータを（ＥＳ，ＥＳａ）と称する。例えば、（ＥＳ，ＥＳａ）＝（σＳ，σＳａ）である。
【０１６３】
本発明による方法は、ベクトルａおよびｂ、ベクトルΘ、パラメータのベクトルΘＳおよびΘＳａ、ならびにベクトルの組（ＥＳ，ＥＳａ）に鑑みて測定ベクトルの組（ｓ，ｓａ）に結合直接確率分布を割り当てる設定ステップ５４を含む。したがって、この結合直接確率分布は、ｐ（ｓ，ｓａ｜ａ，ｂ，Θ，ΘＳ，ΘＳａ，ＥＳ，ＥＳａ，ＩＳ，ＩＳａ，Ｍ）と表わされる。
【０１６４】
例えば、
【０１６５】
【数３８】
2013539704.tif 000040
【０１６６】
である（ξおよびξａを、標準偏差σＳおよびσＳａをそれぞれ有する相互に独立な白色の定常ガウス確率過程と考えた場合）。あるいは、本発明は、所望であれば、前記結合直接確率分布を別なやり方で、未知の理論上の相補累積分布関数の値のベクトルｂにて直接にではなく、インパルス応答の値のベクトルｈにて表現することを可能にできる。この目的のために、ｂをｈに関して表現するだけでよい。したがって、結合直接確率分布５４が、例えば下記のように表わされる。
【０１６７】
【数３９】
2013539704.tif 000041
【０１６８】
一般性を失うことなく、上述した理由でベクトルｈが優先的に注目されることを後に考慮すべきである。
【０１６９】
同じやり方で、複素信号の測定値のベクトルの組についての結合直接確率分布を、それらの実部および虚部の直接確率分布を乗算することによって割り当てることもできる。
【０１７０】
さらに、本発明は、例えば核磁気共鳴による灌流強調イメージングにおいて
【０１７１】
【数４０】
2013539704.tif 000042
【０１７２】
を使用することによって灌流信号ｓおよびｓａの値のベクトルを濃度ｃｅｘｐの値のベクトルへと正確に変換できるときに変種を提供する。例えば、造影剤の到着前の平均の灌流信号強度Ｓ０を正確かつ他のパラメータとは別個独立に推定できるとき、そのようにすることが可能である。この実施の形態によれば、他のすべてのパラメータに鑑みたベクトルの組（ｃｅｘｐ，ｃｅｘｐａ）についてのガウス確率分布を下記のように割り当てることが依然として正当であることを示すことができ、
【０１７３】
【数４１】
2013539704.tif 000043
【０１７４】
σおよびσａは、今やｃｅｘｐおよびｃｅｘｐａのそれぞれの測定雑音の標準偏差である。
【０１７５】
より一般的には、実験による灌流信号ｓの値のベクトルならびにその濃度−時間曲線ｃｅｘｐの値のベクトルへの変換を指して、用語＜＜実験灌流データ＞＞が使用される。その後に、一般性を失うことなく、実験灌流データをｓおよびｓａと称する。
【０１７６】
さらに、この方法は、モデルＭの血行動態パラメータΘのついての情報ＩΘ、インパルス応答ｈについての情報Ｉｈ、および動脈入力関数ａについての情報Ｉａをそれぞれ導入することからなる３つの設定ステップ５１ａ、５１ｂ、および５１ｃを含む。設定ステップ５１ａ〜５１ｅにおいて導入される情報は、本発明によれば、処理ユニットを設定し、すなわち処理ユニットが事後周辺分布の割り当て５６を行なうことができ、次いで関心対象量の推定５７を行なうことができるようにするための設定パラメータを構成する。それらの情報または設定パラメータから、本発明による方法は、ｐ（ａ，ｈ，Θ，ΘＳ，ΘＳａ，ＥＳ，ＥＳａ｜ＩΘ，Ｉｈ，Ｉａ，ＩＳ，ＩＳａ，Ｍ）と表わすことができる結合事前確率分布を割り当てるステップ５３を含むことができる。
【０１７７】
最大エントロピーの原理を適用し、この分布を、下記のように典型的に因数分解することができる。
【０１７８】
【数４２】
2013539704.tif 000044
【０１７９】
したがって、そのような方法は、事前確率分布ｐ（Θ｜ＩΘ，Ｍ）を割り当てることからなるステップ５２ａを含む。
【０１８０】
例えば、無情報事前分布を割り当てることができる。例えば、情報ＩΘが、ＢＦおよびτがそれぞれ区間［ＢＦｍｉｎ，ＢＦｍａｘ］および［τｍｉｎ，τｍａｘ］に属することを知ることだけからなる場合、事前確率分布ｐ（Θ｜ＩΘ，Ｍ）を、下記のように表わすことができる。
【０１８１】
【数４３】
2013539704.tif 000045
【０１８２】
他端において、過去の実験から得られ、例えば陽電子放出断層撮影（ＰＥＴ）または動脈スピンラベリング（ＡＳＬ）などの定量的な灌流強調イメージング技術から得られる相対頻度サンプリング分布または事後周辺確率分布などの情報確率分布を割り当てることも可能である。
【０１８３】
本発明による方法は、情報Ｉｈ（あるいは、相補累積分布関数だけが関心の対象である場合には、Ｉｂ）からの事前確率分布
【０１８４】
【数４４】
2013539704.tif 000046
【０１８５】
を割り当てることからなる設定ステップ５２ｂをさらに含む。この情報は、確実および不確実情報で作られる。
【０１８６】
本発明によれば、確実情報が、不確実情報から区別される。確実情報は、確かであると考えられ、すなわち確率が１に等しいあらゆるブール命題（Ｂｏｏｌｅａｎ  ｐｒｏｐｏｓｉｔｉｏｎ）に相当する。例えば、＜＜この曲線は滑らかである＞＞または＜＜この信号はこのモデルに従う＞＞などのブール命題は、確実情報を構成する。対照的に、不確実情報は、ブール命題を示すことからなるが、そのようなブール命題を特定の確率で示すにすぎないあらゆるブール命題に関する。最後に、これは、＜＜この曲線はほぼ滑らかである＞＞または＜＜この信号はこのモデルにほぼ従う＞＞などのブール命題を導入することを意味する。以後、確実でないあらゆる不確実情報を＜＜完全不確実情報＞＞と呼ぶ。
【０１８７】
確実情報と不確実情報との間の相違を強調するために、最も単純かつ最も極端な事例を検討する。ガウス分布Ｎ（ｘ，σ２）に従って独立かつ同一に分布したサンプルｘ１，…，ｘｎから実際の量ｘを推定したいと仮定する。実数ａを設定する。事前の確実な定量的情報＜＜ｘ＝ａ＞＞を考える。この確実情報は、ディラック事前確率分布ｐ（ｘ｜ａ）＝δ（ｘ−ａ）へと変換され、すなわちｘ＝ａであればｐ（ｘ｜ａ）＝１であり、ｘ≠ａであればｐ（ｘ｜ａ）＝０である。
【０１８８】
次に、対応する不確実な定量的情報、すなわち＜＜ｘがａにほぼ近い＞＞、つまりは
【０１８９】
【数４５】
2013539704.tif 000047
【０１９０】
を考える。
【０１９１】
最大エントロピーの原理を適用し、この不確実な定量的情報は、ガウス事前確率分布へと変換される。
【０１９２】
【数４６】
2013539704.tif 000048
【０１９３】
確実情報および対応する不確実情報が、異なる事前確率分布へと変換されることを理解できるであろう。詳細には、不確実情報についての事前分布において、事前確率分布ｐ（ε）を割り当てることによって＜＜ほぼ＞＞を定量的なやり方で先験的に変換することを可能にする新たな（ハイパー）パラメータεが現れる。詳細には、確実情報は、εを０へと向けることによる不確実情報の限定の事例と理解される。しかしながら、本発明は、ε＞０である完全不確実情報だけを取り扱う。
【０１９４】
確実および対応する完全不確実な事前確率分布は常に異なっており、通常はそれ自身が異なる関心対象量の推定をもたらす。上述の極端な例では、実際に、確実な場合にはｘ＝ａが既知であるためｘを推定する必要がない。他方で、不確実な場合には、ガウス分布の数学的期待値を推定するという古典的な問題に帰着する。
【０１９５】
次に、第２の例として、＜＜ｙが（セミ）パラメトリックモデルＭ（ｘ，Θ）に従う＞＞という古典的な確実な定量的情報を考える（ここで、Θは、Ｍのパラメータのベクトルであり、ｙは実数である）。本質的に、この定量的情報は、ディラック事前確率分布へと変換される。
【０１９６】
【数４７】
2013539704.tif 000049
【０１９７】
ここで、ｚがｚ〜Ｎ（ｙ，σ２）などの雑音の多い実験データである場合、尤度関数または直接確率分布
【０１９８】
【数４８】
2013539704.tif 000050
【０１９９】
が得られ、ここから例えばベイズ規則を適用することによって関心対象量Θおよびσを推定することができる。
【０２００】
【数４９】
2013539704.tif 000051
【０２０１】
次に、＜＜ｙが（セミ）パラメトリックモデルＭ（ｘ，Θ）にほぼ従う＞＞という対応する不確実な定量的情報を検討する。この不確実情報は、以下のように変換される。
【０２０２】
【数５０】
2013539704.tif 000052
【０２０３】
最大エントロピーの原理を適用し、この不確実情報は、下記の事前確率分布へと変換される。
【０２０４】
【数５１】
2013539704.tif 000053
【０２０５】
上述のようにｚ〜Ｎ（ｙ，σ２）である場合、尤度関数または直接確率分布
【０２０６】
【数５２】
2013539704.tif 000054
【０２０７】
が得られ、ここから例えばベイズ規則を適用することによって関心対象量を推定することができる。
【０２０８】
【数５３】
2013539704.tif 000055
【０２０９】
この第２の例から、確実情報および対応する完全不確実情報が、第１の例のように事前確率分布のレベルにおいてのみならず、尤度関数のレベルにおいても相違しうることを、見て取ることができる。やはり、一般に、関心対象量の異なる推定が最終的に得られる。この例は、本発明に従って（セミ）パラメトリックモデルおよび確定的な血行動態パラメータ推定法を＜＜不確実に＞＞するやり方を大まかに示している。
【０２１０】
定性的情報を説明するための第３の例を検討する。Ｓ（ｔ）が実際の信号であり、ｓ＝［Ｓ（ｔ１），…，Ｓ（ｔＮ）］がその値のベクトルであるとする。＜＜Ｓ（ｔ）が滑らかである＞＞という確実な定性的情報を検討する。この定性的情報が、
【０２１１】
【数５４】
2013539704.tif 000056
【０２１２】
という定量的情報へと正規に（ｃａｎｏｎｉｃａｌｌｙ）変換され、ここで
【０２１３】
【数５５】
2013539704.tif 000057
【０２１４】
は、関数のＬ２ノルムを表わす。したがって、Ｓ（ｔ）は線形であり、Ｓ（ｔ）＝ａｔ＋ｂなどの（ａ，ｂ）が存在する。実際、直線よりも滑らかなものは存在しない。この確実情報が、事前確率分布ｐ（Ｓ（ｔ）｜ｔ，ａ，ｂ）＝δ［Ｓ（ｔ）−ａｔ−ｂ］へと最終的に変換される。したがって、実務においては、組（ａ，ｂ）の推定に帰着する。
【０２１５】
ここで、＜＜Ｓ（ｔ）がほぼ滑らかである＞＞という対応する完全不確実な定性的情報を考える。これが、
【０２１６】
【数５６】
2013539704.tif 000058
【０２１７】
という完全不確実な定量的情報へと変換される。関数空間を扱っているため、この不確実情報を変換する事前確率プロセスを得るために最大エントロピーの原理を直接適用することができないように思われる。しかしながら、
【０２１８】
【数５７】
2013539704.tif 000059
【０２１９】
という定量的情報の離散時間バージョンを考えることができ、ここでｄ２ｓ／ｄｔ２は、ｓからｄ２ｓ／ｄｔ２＝Ｄｓなどの二次導関数の有限差分数値近似を使用することによって得られる。
【０２２０】
したがって、以下のようになる。
【０２２１】
【数５８】
2013539704.tif 000060
【０２２２】
次いで、最大エントロピーの原理の適用により、ベクトルｓについての多変量ガウス事前確率分布がもたらされる。
【０２２３】
【数５９】
2013539704.tif 000061
【０２２４】
やはり、確実情報と完全不確実情報との間の相違を評価することができる。確実情報の場合には、ｓｉが直線ｙ＝ａｔ＋ｂ上に位置するように制約される一方で、完全不確実情報の場合には、推定しなければならないε＞０の値に応じて直線から離れるように多少なりとも移動する。第１の場合が、２つのパラメータ（ａ，ｂ）の推定に帰着する一方で、第２の場合は、Ｎ＋１個のパラメータの推定に帰着する（ベクトルｓについてのＮ個のパラメータおよびεハイパーパラメータ）。
【０２２５】
本発明による不確実および完全不確実情報の考え方を明確に定義したので、インパルス応答ｈ（ｔ）についての確実および不確実情報に戻る。
【０２２６】
確実な定量的情報は、ｈ（ｔ１）＝ｈ（０）＝０を含み、例えば右矩形法による
【０２２７】
【数６０】
2013539704.tif 000062
【０２２８】
を含む。両方の情報により、推定すべき値／パラメータの数をＮ−２に減らすことができる。理解できるとおり、賢明な選択は、値ｈ’＝［ｈ（ｔ２），…，ｈ（ｔＮ−１）］Ｔを保ち、ｈ（ｔＮ）を
【０２２９】
【数６１】
2013539704.tif 000063
【０２３０】
によって表現することからなる。さらに、確実な定量的情報
【０２３１】
【数６２】
2013539704.tif 000064
【０２３２】
も考慮することができる。したがって、それらの確実な定量的情報を純粋に定性的な完全不確実情報と組み合わせることによってベクトルｈ’の結合事前確率分布を割り当てることが残る。例えば、上述のように、＜＜ｈ（ｔ）がほぼ滑らかである＞＞という完全不確実な定性的情報Ｉｈ１を考える。すでに述べたように、この定性的情報Ｉｈ１を、最初に
【０２３３】
【数６３】
2013539704.tif 000065
【０２３４】
という完全不確実な定量的情報へと変換することができる。
【０２３５】
サンプリング時点ｔｉ＝１，Ｎにおける時間離散化の後で、ｈの二次導関数の二次数値近似が、例えば下記によって与えられる。
【０２３６】
【数６４】
2013539704.tif 000066
【０２３７】
例えば下記の四次の式など、さらに高次の数値近似も使用することが可能である。
【０２３８】
【数６５】
2013539704.tif 000067
【０２３９】
それらの数値近似を、行列表示にて
【０２４０】
【数６６】
2013539704.tif 000068
【０２４１】
と表わすことができ、ここでＤは、例えば二次の近似の場合に次元Ｎの正方行列
【０２４２】
【数６７】
2013539704.tif 000069
【０２４３】
であり、ｈ’は上記にて定めたとおりである。
【０２４４】
したがって、ｈ’の二次導関数のユークリッドノルムの平方を、
【０２４５】
【数６８】
2013539704.tif 000070
【０２４６】
と表わすことができる。
【０２４７】
したがって、インパルス応答ｈ（ｔ）が多少なりとも滑らかになると定性的に仮定することを、
【０２４８】
【数６９】
2013539704.tif 000071
【０２４９】
が多少なりとも大きくなると仮定することによって定量的に変換することができる。したがって、同じユークリッドノルム
【０２５０】
【数７０】
2013539704.tif 000072
【０２５１】
を有するすべての連続確率分布の中から境界点ｈ（ｔ１）およびｈ（ｔＮ）が別々に取り扱われる事前確率分布ｐ（ｈ’｜ｈ（ｔ１），ｈ（ｔＮ），Ｉ）が探し求められる。ハイパー象限［０，＋∞］Ｎ−２を台とするそれらのすべての分布の中から最も不確か、したがって最も正当であるがゆえに最大の微分Ｓｈａｎｎｏｎエントロピー（Ｌｅｂｅｓｇｕｅ基準測度のもとで）を有する分布を選択することからなる最大エントロピーの原理の適用により、一定のベクトルの数学的期待値Ｍ＝（μ１，…，μ１）Ｔおよび共分散行列ε１２Ｗ１−１／σ２を有する［０，＋∞］Ｎ−２についての条件付き多変量打ち切りガウス分布（あるいは、相補累積分布関数の値のベクトルについて［０，１］Ｎ−２で打ち切られた多変量ガウス分布）が得られる。
【０２５２】
【数７１】
2013539704.tif 000073
【０２５３】
ここでＣ１（μ１，ε１，σＳ）は、正規化定数
【０２５４】
【数７２】
2013539704.tif 000074
【０２５５】
である。
【０２５６】
したがって、２つの新たなハイパーパラメータ、すなわち我々の灌流モデルの正則化パラメータの役割を果たすスカラー数学的期待値μ１および逆分数分散ε１が、我々の大域的灌流モデルに現れる。ε１は、実験データｓおよびｓａによってもたらされる確実な定量的情報に対する事前の不確実な定性的情報Ｉｈ１の不確実さを定量化する。
【０２５７】
さらに、定義により、下記のとおりである。
【０２５８】
【数７３】
2013539704.tif 000075
【０２５９】
しかしながら、遅延τがサンプリング時点ｔｉ，ｉ＝１，Ｎにほとんど等しくないという事実を考慮するために、
【０２６０】
【数７４】
2013539704.tif 000076
【０２６１】
とすることができ、ここで
【０２６２】
【数７５】
2013539704.tif 000077
【０２６３】
かつＣ１，１（μ１，１，ε１，１，σＳ）は正規化定数
【０２６４】
【数７６】
2013539704.tif 000078
【０２６５】
であり、
【０２６６】
【数７７】
2013539704.tif 000079
【０２６７】
だけが示される。
【０２６８】
最後に、情報Ｉｈ１についてハイパー象限［０，＋∞］Ｎを台とするｈの事前確率分布を、以下のように表わすことができる。
【０２６９】
【数７８】
2013539704.tif 000080
【０２７０】
一般に、例えばＥｈ＝（μ１，ε１）またはＥｈ＝（μ１，ε１，μ１，１，ε１，１）など、ｈの事前分布のハイパーパラメータのベクトルがＥｈ（または、Ｅｂ）と称される。
【０２７１】
Ｅｈの事前確率分布を割り当てることが残る。例えば、無情報非正則Ｂａｙｅｓ−Ｌａｐｌａｃｅ−Ｌｈｏｓｔｅ−Ｊｅｆｆｒｅｙｓ分布ｐ（μ１，ε１｜Ｉｈ１，Ｍ）∝ε１−１またはｐ（μ１，ε１，μ１，１，ε１，１｜Ｉｈ１，Ｍ）∝ε１−１ε１，１−１を割り当てることができる。
【０２７２】
同じやり方で、
【０２７３】
【数７９】
2013539704.tif 000081
【０２７４】
という完全不確実な定量的情報Ｉｈ２を導入でき、より一般的には、
【０２７５】
【数８０】
2013539704.tif 000082
【０２７６】
という完全不確実情報Ｉｈｋ，ｋ＝１，２，…を導入することができる。
【０２７７】
最後に、最大エントロピーの原理を適用することによってｈの新たな事前多変量打ち切りガウス確率分布が得られる。
【０２７８】
上述のような不確実情報Ｉｈ１，Ｉｈ２，…が唯一の完全不確実情報であり、我々の問題に正当に導入することができる唯一の制約であり、それらが検証できる任意の仮説ではなく、あるいは実験によるものでなく、むしろ対照的に、それらは基本的な生理学的特性の表現にすぎず、それなしでは血行動態パラメータ、インパルス応答、または相補累積分布関数の推定の問題が実際に絶対的に無意味になりかねず、それらが我々の問題にとって論理的に必要かつ充分であることに、注目すべきである。他のあらゆる情報は、対照的に、潜在的に実験によって検証されることができる単純仮説であると考えられる。
【０２７９】
本発明によれば、それでもやはり、単純作業仮説である他の完全不確実情報を導入することができる。例えば、インパルス応答が所与の関数形式にほぼ従うにすぎない旨を、確実情報によってこの形式に正確に従うように制約するのではなく、示すことができる。この種の半定量的な完全不確実情報の追加は、提案された関数形式でインパルス応答をどの程度まで表わすことができるかを判断できるようにする。したがって、例えば＜＜ｈ（ｔ）がパラメトリックまたはセミパラメトリックモデルＭｈ：ｆ（ｔ，Θｈ）にほぼ従う＞＞という上述のような完全不確実情報Ｉｈ３を導入したいと仮定する（ここで、Θｈはこのモデルのパラメータのベクトルである）。すでに述べたように、そのようなパラメトリックモデルは、例えば下記のように与えられる２パラメータΓモデルなど、すでに提案されている。
【０２８０】
【数８１】
2013539704.tif 000083
【０２８１】
この場合に、パラメータＭＴＴを、上述のようにインパルス応答の第１のモーメント（または、相補累積分布関数の積分）を数値的に推定する必要なく、直接推定できることに注目すべきである。
【０２８２】
ｈが所与の関数形式ｆ（ｔ，Θｈ）にほぼ従う旨を示すことは、剰余のベクトルのユークリッドノルム
【０２８３】
【数８２】
2013539704.tif 000084
【０２８４】
が多少なりとも大きい旨を定量的に示すことを意味する。やはり最大エントロピーの原理を適用し、同じやり方で、ｈの事前確率分布がハイパー象限［０，＋∞］Ｎで打ち切られた多変量分布であることが分かる。
【０２８５】
【数８３】
2013539704.tif 000085
【０２８６】
ここでＣ３（ε３，σＳ，Θｈ）は、正規化定数
【０２８７】
【数８４】
2013539704.tif 000086
【０２８８】
である。
【０２８９】
同じやり方で、＜＜ｈ（ｔ）が所与の値ｈ−＝［ｈ−（ｔ１），…，ｈ−（ｔＮ）］Ｔのベクトルにほぼ従う＞＞などの不確実情報Ｉｈ４を導入することができる。
【０２９０】
最大エントロピーの原理を適用し、事前確率分布
【０２９１】
【数８５】
2013539704.tif 000087
【０２９２】
を得る。
【０２９３】
ここでＣ４（ε４，σ，ｈ−）は、正規化定数
【０２９４】
【数８６】
2013539704.tif 000088
【０２９５】
である。
【０２９６】
さらに、本発明は、インパルス応答ｈ（ｔ）（または、相補累積分布関数）についてのいくつかの情報およびそれらの対応する事前確率分布を組み合わせることを可能にする。したがって、ｐ（ｈ｜Ｅｈ１，Ｉｈ１，Ｍ），…，ｐ（ｈ｜Ｅｈｎ，Ｉｈｎ，Ｍ）がそれぞれの正則化パラメータＥｈ１，…，Ｅｈｎ（例えば、Ｅｈ１＝（ε１，μ１）、Ｅｈ３＝（ε３，Θｈ）、など）で情報Ｉｈ１，…，Ｉｈｎを変換するｎ個の事前確率分布である場合、それらｎ個の情報を考慮するｈの事前分布を、Ｅｈ＝（Ｅｈ１，…，Ｅｈｎ）およびＩｈ＝Ｉｈ１∧…∧Ｉｈｎとすることによって
【０２９７】
【数８７】
2013539704.tif 000089
【０２９８】
と表わすことができる。
【０２９９】
情報Ｉａを局所的動脈入力関数へとエンコードするために、本発明による方法は、事前確率分布
【０３００】
【数８８】
2013539704.tif 000090
【０３０１】
を割り当てることからなるステップ５２ｃをさらに含む。そのような分布は、動脈入力関数についての確実および／または不確実情報を導入して組み合わせることによって、インパルス応答ｈに関する分布と同じやり方で割り当てられる。例えば、動脈入力関数がほぼ滑らかであり、正であり、単峰であり、双峰であり、原点においてゼロであり、ゼロに漸近であり、あるいは所与の範囲を有する旨を指定することができ、もしくは所与のパラメータまたはセミパラメトリックモデルＣａ（ｔ，Θａ）にほぼ従うことだけを示すことができ、ここでΘａは、例えば１１パラメータの＜＜３ガンマ＞＞モデルなどの、パラメータのベクトルである。
【０３０２】
【数８９】
2013539704.tif 000091
【０３０３】
すでに述べたように、ガウス事前確率分布
【０３０４】
【数９０】
2013539704.tif 000092
【０３０５】
を得る。
【０３０６】
ここでＣ（εａ１，σＳａ，Θａ）は、正規化定数
【０３０７】
【数９１】
2013539704.tif 000093
【０３０８】
である。
【０３０９】
同様に、ベクトルａが局所的動脈入力関数などの所与の値のベクトルａ−＝［ａ−（ｔ１），…，ａ−（ｔＮ）］Ｔにほぼ近いことだけを示す（そうでない場合、動脈入力関数が与えられる事例に帰着する）ことができ、事前確率分布
【０３１０】
【数９２】
2013539704.tif 000094
【０３１１】
を得る、ここでＣ（εａ，σａ，ａ−）は、正規化定数
【０３１２】
【数９３】
2013539704.tif 000095
【０３１３】
である。
【０３１４】
本発明による方法は、典型的には以下のように因数分解される事前確率分布ｐ（ＥＳ，ＥＳａ，ΘＳ，ΘＳａ｜ＩＳ，ＩＳａ，Ｍ）を割り当てることからなる設定ステップ５２ｄをさらに含む。
【０３１５】
【数９４】
2013539704.tif 000096
【０３１６】
例えば、無情報事前確率分布が得られる。
【０３１７】
【数９５】
2013539704.tif 000097
【０３１８】
ステップ５２ｂ〜５２ｄにおいて導入された種々のハイパーパラメータを考慮し、結合事前確率分布５３を、ｐ（ａ，Ｅａ，ｈ，Ｅｈ，Θ，ΘＳ，ΘＳａ，ＥＳ，ＥＳａ｜ＩΘ，Ｉｈ，Ｉａ，ＩＳ，ＩＳａ，Ｍ）と書き直すことができる。
【０３１９】
以下の表現を簡単にするために、処理ユニットの設定パラメータとして入力された情報一式をＩ＝（ＩΘ，Ｉｈ，Ｉａ，ＩＳ，ＩＳａ，Ｍ）とする。
【０３２０】
上述のとおりの直接確率分布５４および結合事前確率分布５３に鑑み、ベイズ規則を適用することによってすべてのパラメータの結合事後確率分布５５が得られる。
【０３２１】
【数９６】
2013539704.tif 000098
【０３２２】
初期設定が行なわれると、今や本発明は、ベクトルΞ＝（ａ，Ｅａ，ｈ，Ｅｈ，Θ，ΘＳ，ΘＳａ，ＥＳ，ＥＳａ）のすべての要素中からθと称するべき関心対象量の推定を可能にする。例えば、θ＝ＢＦまたはθ＝τ、あるいはθ＝ｈ（ｔｉ）などである。
【０３２３】
したがって、本発明による方法は、θの事後周辺分布を評価することからなるステップ５６、すなわち
【０３２４】
【数９７】
2013539704.tif 000099
【０３２５】
を含む。
【０３２６】
この事後周辺分布から、θの推定θ＾を入手５７することができる。例えば、この分布の数学的期待値
【０３２７】
【数９８】
2013539704.tif 000100
【０３２８】
をとることによって、二次損失関数
【０３２９】
【数９９】
2013539704.tif 000101
【０３３０】
のもとでのベイズ推定量が得られる。同じやり方で、
【０３３１】
【数１００】
2013539704.tif 000102
【０３３２】
によって最大事後推定量θ＾Ｐ、すなわちＭＡＰを得ることができる。
【０３３３】
例として、各々のサンプリング時点ｔｉにおけるインパルス応答の値ｈ（ｔｉ），ｉ＝１，Ｎの事後周辺確率分布を、他のすべての時点を周辺化することによって得ることができ、
【０３３４】
【数１０１】
2013539704.tif 000103
【０３３５】
その後にｈ＾Ｑ（ｔｉ）またはｈ＾Ｐ（ｔｉ）などのそれらの値の推定を得ることができる。
【０３３６】
同じやり方で、必ずしもサンプリング時点ｔｉに等しくない任意の時点ｘにおけるインパルス応答ｈ（ｘ）の値の推定ｈ＾（ｘ）（または、Ｒ＾（ｘ））を入手５７することも可能である。この目的のために、上述のような相補累積分布関数Ｒ（ｔ）の値の表現にｈ（ｘ）を導入し、その周辺確率分布を計算すれば充分である。そのような追加の時点ｘ１，…，ｘＬを推定の問題に導入することは、考慮される時点の数が多いほど
【０３３７】
【数１０２】
2013539704.tif 000104
【０３３８】
などの積分の数値近似ならびにｄ２ｈ／ｄｔ２またはｄｈ／ｄｔの数値近似がより良好になるため、さらに望ましい。その後に、得られる推定もより良好である。
【０３３９】
したがって、
【０３４０】
【数１０３】
2013539704.tif 000105
【０３４１】
であるため、例えばここでは右矩形法などの数値積分法を適用することによって、＜＜平均＞＞推定
【０３４２】
【数１０４】
2013539704.tif 000106
【０３４３】
または最確推定
【０３４４】
【数１０５】
2013539704.tif 000107
【０３４５】
などといったこのパラメータの推定を得ることができる。
【０３４６】
さらに、本発明による方法は、パラメータθの推定の精度の推定を得、さらにはそれらの推定についての信頼区間を得るステップ５８を含むことができる。本発明によれば、そのような方法が、前記信頼区間についての賭け率を得るステップ５９をさらに備えることができる。例えば、推定θ＾Ｑの精度を、θの事後周辺確率分布の共分散行列によって定量化することができる。
【０３４７】
【数１０６】
2013539704.tif 000108
【０３４８】
したがって、例えばθに対する＜＜１シグマ＞＞における信頼（ハイパー）区間
【０３４９】
【数１０７】
2013539704.tif 000109
【０３５０】
ならびにθがこの区間に属する確率
【０３５１】
【数１０８】
2013539704.tif 000110
【０３５２】
または同等のオッズｐＪ／（１−ｐＪ）が得られる。
【０３５３】
また、パラメータＢＶ＝ＢＦ・ＭＴＴ、相補累積分布関数の値のベクトルｂ、静脈出力関数の値のベクトルｖ＝Ａｈ、または理論上の濃度−時間曲線の値のベクトルｃ＝ＢＦ・Ａｂについての推定、信頼区間、およびそれら信頼区間への賭け率を得ることも、いくつかのランダム変数の結合確率分布に鑑みてそれらの任意の関数の確率分布関数を計算できるがゆえに可能である。例えば、数学的期待値の線形性によって、相補累積分布関数Ｒ（ｔ）の値の推定が、例えば右矩形法に従うことによってインパルス応答
【０３５４】
【数１０９】
2013539704.tif 000111
【０３５５】
および
【０３５６】
【数１１０】
2013539704.tif 000112
【０３５７】
のそれらからすぐに得られる。同じやり方で、ｃおよびΘＳについての結合確率分布ｐ（ｃ，ΘＳ｜ｓ，ｓａ，Ｉ）から、理論上の灌流信号Ｓｔｈ（ｔ）の値の確率分布ｓｔｈ＝［Ｓｔｈ（ｔ１），…，Ｓｔｈ（ｔＮ）］Ｔを、ｓｔｈ＝Ψ（ｃ，ΘＳ）であるがゆえに導出することもできる。この分布から、上述した方法と同様の方法に従って、推定ｓｔｈ＾を入手５７することができ、信頼ハイパー区間５８およびそれらハイパー区間についての賭け率５９も入手することができる。
【０３５８】
本発明によれば、信頼区間またはそれらの信頼区間の賭け率が、処理ユニットの設定のセッティング６２を可能にすることができる。したがって、設定パラメータを変更し、より高品質の推定をもたらすことができる。
【０３５９】
さらに、本発明による方法は、理論上および実験による灌流強調信号の間の剰余ｒ（ｔｉ）＝Ｓ（ｔｉ）−ｓｔｈ＾（ｔｉ），ｉ＝１，Ｎを計算するステップ６０を含むことができる。したがって、本発明は、それらのベクトルの間の種々の統計量または距離Ｄ（ｓ，ｓｔｈ＾）の計算を可能にし、その最も古典的な１つは、誤差の平方和（ＳＳＥ）
【０３６０】
【数１１１】
2013539704.tif 000113
【０３６１】
である。それらの種々の統計量は、灌流モデルの実験データへの適合度を定量化できるようにする。このようにして、関心対象の各ボクセルにおけるモデルの＜＜誤差マップ＞＞が得られる。上述の理由（すなわち、過剰適合の問題）で、灌流モデルの実験データｓおよびｓａへの適合度の定量化を、関心対象の各ボクセルの灌流モデルを鑑みて実験データ（ｓ，ｓａ）の確率を計算することによってきわめて好都合に達成することができる。
【０３６２】
【数１１２】
2013539704.tif 000114
【０３６３】
この場合には、誤差マップが１−ｐ（ｓ，ｓａ｜Ｉ）にもとづくと考えられる。
【０３６４】
さらに本発明によれば、剰余ｒ（ｔｉ）が実際に独立であるかどうか、等しく分布しているかどうか、ガウスであるかどうか、などをチェックするために、Ｑ−Ｑプロットまたはポアンカレの帰還写像などの種々の統計テストまたは種々のグラフィカル診断技術の適用６１を行なうこともできる。したがって、本発明によれば、反復プロセスおよび試行錯誤によって、灌流現象のモデル化、理解、および処理を進歩させるとともに、血行動態パラメータ、インパルス応答、相補累積分布関数、動脈入力関数、または静脈出力関数の良好な推定を精密に得るために、設定プロセス５０（詳細には、理論上の灌流モデル）の修正および洗練６２を行なうことができる。
【０３６５】
次に、図８に関連して、本発明による方法を、理論上の動脈入力関数が遅延τまで絶対的な確かさおよび無限の精度で与えられると仮定される第２の応用の例について説明する。
【０３６６】
したがって、実験による灌流モデルＭを
【０３６７】
【数１１３】
2013539704.tif 000115
【０３６８】
と表わすことができ、ここですべての量は、今や既知であると仮定される理論上の動脈入力関数の値のベクトルであるａ＝［Ｃａ（ｔ１），…，Ｃａ（ｔＮ）］Ｔを除き、すでに定めたとおりである。
【０３６９】
したがって、結合直接確率分布が、ｐ（ｓ｜ａ，ｈ，Θ，ΘＳ，ＥＳ，Ｉ）と表わされ、結合事前確率分布が、ｐ（ｈ，Ｅｈ，Θ，ΘＳ，ＥＳ｜Ｉ）と表わされ、Ｉ＝（ＩΘ，Ｉｈ，ＩＳ，Ｍ）である。これらの分布が、上述のように割り当てられる。ベイズ規則は、下記のようになる。
【０３７０】
【数１１４】
2013539704.tif 000116
【０３７１】
続いて、すでに説明した同じやり方で、任意のパラメータ
【０３７２】
【数１１５】
2013539704.tif 000117
【０３７３】
の推定、信頼区間、およびそれらの信頼区間の賭け率が得られる。
【０３７４】
次に、図８に関連して、本発明による方法を、理論上の動脈入力関数が与えられず、動脈入力信号も測定されない第３の応用の例について説明する。
【０３７５】
したがって、実験による灌流モデルＭを
【０３７６】
【数１１６】
2013539704.tif 000118
【０３７７】
と表わすことができ、ここですべての量は、すでに定めたとおりである。
【０３７８】
結合直接確率分布が、やはりｐ（ｓ｜ａ，ｈ，Θ，ΘＳ，ＥＳ，Ｉ）と表わされ、結合事前確率分布は、今やｐ（ａ，Ｅａ，ｈ，Ｅｈ，Θ，ΘＳ，ＥＳ｜Ｉ）と表わされ、Ｉ＝（ＩΘ，Ｉｈ，Ｉａ，ＩＳ，Ｍ）である。ベイズ規則は、下記のようになる。
【０３７９】
【数１１７】
2013539704.tif 000119
【０３８０】
続いて、すでに説明した同じやり方で、任意のパラメータ
【０３８１】
【数１１８】
2013539704.tif 000120
【０３８２】
の推定、信頼区間、およびそれらの信頼区間の賭け率が得られる。
【０３８３】
該当の応用の例がどれであっても、公知の方法と対照的に、本発明による推定方法が、関心対象量について先験的に有する定性的、定量的、または半定量的な情報を、実験による測定によってもたらされるそれらの量についての事後情報を一意に決定するために、ベイズ確率理論へと変換することだけからなる点で、正確な方法であることに、注目すべきである。本発明は、問題を解決するために論理的に必要かつ充分な種々の不確実な定性的制約をエンコードしている最も不確かな確率分布だけを導入するため、検証されることができ、あるいは（実際には可能であっても）検証されることができない（詳細には、インパルス応答または相補累積分布関数についての）恣意的な仮説が、不要である。
【０３８４】
したがって、動脈入力関数が与えられず、せいぜい測定されるにすぎない場合でも、それらの方法は、提案された動脈入力信号が関心対象の各ボクセルの＜＜真の＞＞局所的動脈入力関数に効果的に対応することができるか否かをテストすることを可能にし、実際に、動脈入力信号が適しておらず、真の局所的動脈入力関数に対応しない場合に、問題の解であると同時に種々の事前制約を満足する（または、少なくともその確率が先験的に無視できない）パラメータ（血行動態パラメータ、相補累積分布関数、など）の組が、存在しない可能性がある。この場合、確率論は、提案された動脈入力信号を雑音と解釈し、標準偏差σａ／εａが、より適した動脈入力信号から典型的に得られる標準偏差と比べてはるかに大きくなる。
【０３８５】
したがって、本発明は、技術水準による大域的または局所的動脈入力関数の種々の選択および推定方法をテストするというきわめて興味深い応用を有する。それらの大域的または局所的動脈入力関数から得られた推定（詳細には、それらの正則化パラメータＥａの推定）が、或るボクセルから他のボクセルへと過剰に頻繁に異常であることが明らかになる場合、それらの関数の新たなより適切な（局所的）選択法を導入するか、あるいは動脈入力関数が与えられること、または動脈入力信号が事前に測定されることを必要としない方法に頼るか、のいずれかが必要であると結論付けられ、これがまさに上述した第３の方法の目的である。
【０３８６】
応用の例として、図１または２において説明したような灌流強調解析イメージングシステムによる本発明の主たる実行のステップに言及することができる。
【０３８７】
・処理ユニット４（または、前処理ユニット７）によって、患者の記録が開かれ、あるいは一連の画像が取り込まれ、関心対象の一連の画像が選択され、詳細には図５ａによって示されるように各ボクセルについての灌流信号Ｓ（ｔ）を得るための時間につれての灌流強調画像Ｉ１〜Ｉｎが選択される。
【０３８８】
・関心対象のスライスまたは領域のユーザ６による特定を可能にするためのヒューマン・マシン・インターフェイス５による画像のプレビュー。
【０３８９】
・本発明による推定方法の実行を可能にするための設定パラメータ（入力情報）からの処理ユニット４の設定。
【０３９０】
・推定すべき１つ以上の関心対象量の選択。
【０３９１】
・人間の脳などの器官についてのＢＦ＾、τ＾、またはＭＴＴ＾などの血行動態パラメータ１４の処理ユニット４による推定。
【０３９２】
・動脈入力関数が与えられないときのＥｈ、ＥＳ、ΘＳまたはＥａ、ＥＳａまたはΘＳａなどの他のパラメータの随意による推定。
【０３９３】
・随意による、インパルス応答の値のベクトルの推定ｈ＾、相補累積分布関数の推定ｂ＾、動脈入力関数の値のベクトルの推定ａ＾、静脈出力関数の値のベクトルの推定ｖ＾、理論上の濃度の値のベクトルの推定ｃ＾、理論上の信号の推定ｓｔｈ＾、または剰余の推定ｒ＾。
【０３９４】
・の前記推定による関心対象量１４が、内容を医師へと示すために、例えば各ピクセルの輝度または色が計算された値に応じて決まるマップとして最終的に表示されるよう、ヒューマン・マシン・インターフェイス５へと引き渡される。
【０３９５】
・ユーザによって選択された関心対象のいくつかのボクセルについて、インパルス応答、相補累積分布関数、動脈入力関数、静脈出力関数、理論上の濃度、理論上の信号、または剰余が、随意により表示される。
【０３９６】
・血行動態パラメータなどの関心対象の１つ以上のパラメータについての信頼マップまたは賭け率マップが随意により表示される。
【０３９７】
・ユーザによって選択された関心対象のいくつかのボクセルに関して、いくつかのインパルス応答、相補累積分布関数、動脈入力関数、などについての信頼区間またはそれら区間への賭け率が随意により表示される。
【０３９８】
・実験データと大域的非パラメータ灌流モデルとの間の１つ以上の距離の誤差マップが随意により表示され、詳細には大域的灌流モデルに鑑みた実験データの確率が表示される。
【０３９９】
・１つ以上の血行動態パラメータ、インパルス応答（または、相補累積分布関数）、または局所的動脈入力関数の異常な分布を特徴とする関心対象の前記病気の領域の補助された選択。
【０４００】
・おそらくは病変の領域に関係しており、医師による治療行為（例えば、血栓を分解するための経静脈血栓溶解）の決定の対象となりうる異常な灌流の組織領域部分の処理ユニットによる推定。
【０４０１】
・病変の容積と異常な灌流の領域との間の比など、医師による診断および治療の決定（例えば、血栓を分解するための経静脈血栓溶解）の精緻化を可能にするための何らかの量の処理ユニットによる推定。
【０４０２】
すでに述べたように、処理ユニット４の設定プロセス５０を、ユニットそのものによって実行することができる（設定プロセス５０の実行）。あるいは、前記設定は、推定を望む関心対象量に応じた結合事後確率分布のライブラリの保存および選択するこを含むことができる。このライブラリの構築を、処理ユニット４と協働することができる専用のユニットによって達成することができる。
【０４０３】
あるいは、前記設定の精緻化のために、いくつかの関心対象量についての信頼区間およびそれら信頼区間の賭け率の推定に続いて、繰り返しを行なうことができる。実験データと大域的非パラメータ灌流モデルとの間の距離の提供、詳細には大域的灌流モデルに鑑みた実験データの確率の表示は、設定の更新をさらに含むことができる。
【０４０４】
したがって、本発明は、パラメータの推定を、各ボクセルの輝度または色が計算された値に例えば線形なやり方で依存する＜＜パラメータマップ＞＞の形態で表示することを目的とする。同様に、本発明は、おそらくは、それらの推定の標準偏差を＜＜信頼度マップ＞＞の形態で表示するとともに、対応する信頼区間についての賭け率を＜＜賭け率マップ＞＞の形態で表示することを目的とする。インパルス応答、相補累積分布関数、動脈入力関数、静脈出力関数、濃度−時間曲線、灌流信号、または剰余の値のベクトルの推定に関して、本発明は、それらをユーザによって選択された各ボクセルについて時系列の形態で表示することを目的とする。最後に、本発明は、実験による信号と非パラメータ灌流モデルとの間の距離またはこのモデルに鑑みた実験データの確率を＜＜誤差マップ＞＞の形態で表示することを目的とする。
【０４０５】
図９〜１２が、本発明に従って推定された血行動態パラメータ１４あるいはそれらに関する標準偏差または確率などのいくつかの関心対象量について、マップの形態での表示の態様を説明することを可能にしている。
【０４０６】
したがって、核磁気共鳴イメージングによって分析された人間の脳について、図９が、脳血液容量ＣＢＶの推定を眺めることを可能にする。そのようなマップ（４５８×４５８ピクセル）が、おそらくは虚血性の領域８０の強調を可能にする。実際、適切なインターフェイス５によって、右後大脳動脈の領域において反対側の半球と比べてパラメータＣＢＶが明らかに増加していることを見て取ることができる。図９に示されるマップを読み取ることによって、虚血の後の血管拡張が明らかになる。
【０４０７】
図１０が、虚血性脳梗塞の場合における脳血流の推定に関するマップ（４５８×４５８ピクセル）を示すことを可能にする。マップを分析することによって、虚血の後の右後大脳動脈の領域における反対側の半球と比べたパラメータＣＢＦ（脳血流量）の減少を見て取ることができる。そのようなマップが、おそらくは虚血性の領域８０の強調を可能にする。
【０４０８】
図１１が、平均通過時間ＭＴＴの推定に関するマップ（４５８×４５８ピクセル）を示すことを可能にする。マップを分析することによって、虚血の後の右後大脳動脈の領域８１における反対側の半球と比べたＭＴＴの明らかな増加を見て取ることができる。
【０４０９】
図１２が、脳血流量が信頼区間［ＣＢＦ＾−σ＾ＣＢＦ，ＣＢＦ＾＋σ＾ＣＢＦ］に属する確率の推定に関するマップ（４５８×４５８ピクセル）を示すことを可能にする。マップを分析することによって、異常なモデルの適合を除いて、確率が０．６８のあたりに中心を有することを見て取ることができる。これは、ＣＢＦの事後確率分布がガウスの法則に従う場合に期待することができる値である。
【０４１０】
上述のマップのおかげで、本発明は、技術水準の公知の技術を用いた場合には入手することができなかったさまざまな有用な情報をユーザに提供することを可能にする。この提供は、図１または２に示したとおりの処理ユニット４を、推定したパラメータ１４を例えば図９〜１３によって示されるとおりのマップの形態でユーザ６へと示すことができるヒューマン・マシン・インターフェイス５に適したフォーマットで、外界との通信のための手段がそのパラメータ１４をもたらすことができるように、構成することによって可能にされる。
【０４１１】
本発明によれば、もたらされる情報がより多数かつより公正である。医師にとって利用可能な情報が、診断における判断および治療上の意志決定における医師の自信を高めると考えられる。
【０４１２】
本発明によるシステムの性能を向上させるために、処理ユニットが、血行動態パラメータ、相補累積分布関数、または動脈入力関数の推定が必要とされる画像のボクセルについての演算を並列化するための手段を備えることができる。これは、グラフィカル・プロセッサ・ユニット（ＧＰＵ）などのハードウェア技術、コンピュータクラスタ、または並列モンテカルロ法などのソフトウェア技術、などを使用することによって達成することができる。あるいは、本発明による処理ユニットが、離れて位置する演算手段に頼ることができる。このやり方で、演算の時間をさらに大幅に減らすことができる。
【図１】
2013539704.tif 000121
【図２】
2013539704.tif 000122
【図３】
2013539704.tif 000123
【図４】
2013539704.tif 000124
【図５ａ−５ｂ】
2013539704.tif 000125
【図６】
2013539704.tif 000126
【図７】
2013539704.tif 000127
【図８】
2013539704.tif 000128
【図９−１０】
2013539704.tif 000129
【図１１】
2013539704.tif 000130
【図１２】
2013539704.tif 000131
【国際調査報告】
2013539704.tif 000132
2013539704.tif 000133
2013539704.tif 000134
2013539704.tif 000135
2013539704.tif 000136
2013539704.tif 000137
2013539704.tif 000138
2013539704.tif 000139