JP5918532B2 - 穿通枝血管の灌流を評価する装置及び穿通枝血管の灌流を評価するための装置の作動制御方法 - Google Patents

Description

本発明は、形成外科および再建手術のための穿通枝皮弁（ｐｅｒｆｏｒａｔｏｒ ｆｌａｐｓ）の位置特定（ｌｏｃａｔｉｎｇ）および分析に関する。

関連出願の相互参照
本出願は、２００８年４月１４日出願の米国特許仮出願第６１／０４４，７７９号の優先権を主張するものであり、その全内容を参照により本明細書に組み込む。

形成外科および再建手術には、皮膚および皮下組織から構成される皮弁の位置特定および臨床評価がしばしば伴い、こうした皮弁は、独立した（ｉｓｏｌａｔｅｄ）穿通枝血管によって供給され、したがって身体の別の部位に移植するのに潜在的に適している。穿通枝は、それらの源の血管から、深部筋組織中、または深部筋組織間を通って皮膚表面へと通じている。良好に血管新生した皮弁は、優れた移植片候補である。

例えば、１９８０年代初頭から、腹部の恵皮部皮弁が、自家組織による乳房再建の標準となってきている。腹部内では、遊離脂肪の選択肢として、完全な横行腹直筋皮（ＴＲＡＭ）皮弁から、深下腹壁動脈（ＤＩＥＡ）穿通枝皮弁などの独立した穿通枝皮弁まで幅がある。穿通枝皮弁によって、患者自身の皮膚および脂肪を、恵皮部の合併症を最小限に抑えながら、信頼性の高い手法で他の組織再建部位にも移植することが可能となってきている。ランダムパターン（ｒａｎｄｏｍ ｐａｔｔｅｒｎ）の血液供給に依存していた皮弁は、間もなく有茎の軸型パターン（ａｘｉａｌ ｐａｔｔｅｒｎ）皮弁に取って代わられ、この軸型パターン皮弁では、大量の組織を信頼性高く移植することが可能である。遊離組織移植の出現によって、恵皮部と、受皮部とが適切に整合する可能性が大幅に広がることになった。穿通枝皮弁の使用が増すにつれて、ＤＩＥＡおよびその穿通枝の個々の特定の解剖学的特徴（ａｎａｔｏｍｉｃａｌ ｆｅａｔｕｒｅ）について術前に精通する必要が高まってきており、特に、腹壁への脈管供給の解剖学的構造に大幅な差があることを考えると、その必要がある。

穿通枝（ｐｅｒｆｏｒａｔｏｒ）の位置特定および評価は、労を要し、時間のかかる工程である。しばしば術前にコンピュータ断層撮影血管造影（ＣＴＡ）撮像を実施して、位置特定が行われる。かかる手法には、多額の出費が伴い、また、外科医は、事前に取得した３Ｄ様相の画像を、目下手術台の上に横たわっている患者のその時点の２Ｄ視野と、頭の中で相関させなければならないため、さらに複雑となる。したがって、より好ましい撮像様相の研究が続けられており、近年では、蛍光信号に基づき、皮膚越しに血液循環を査定するインドシアニングリーン（ＩＣＧ）蛍光撮像の使用に関心が寄せられている。近赤外スペクトル領域の放射によって励起すると、発光ピークが約８３０ｎｍのＩＣＧ蛍光が生じる。波長約８００ｎｍの励起光は、例えば、ダイオードレーザ、発光ダイオード（ＬＥＤ）、またはアークランプ、ハロゲンランプなどの他の従来の照明光源を適切な帯域フィルタと共に用いることによって生成することができる。皮膚は、この波長を透過する。

ＩＣＧは、血液タンパク質と強く結合し、以前から心拍出量測定、肝機能評価、および眼血管造影に使用されてきており、有害反応はほとんど見られない。ＩＣＧ蛍光信号の評価を用いて、穿通枝の位置を特定することができる。穿通枝付近の皮膚表面には、一般に周囲組織よりも多量の血液が、より速い血流で蓄積されており、ＩＣＧを注入すると、穿通枝は、周囲組織よりも明るく、かつ急速に蛍光を発する傾向がある。この迅速で高強度の蛍光によって、穿通枝の位置を視覚的に特定することが可能となる。しかし、外科医は、優れた臨床判断をサポートするために、穿通枝の位置特定だけでなく、その評価および比較にも関心を寄せている。外科医は、いくつかの穿通枝のうちのいずれが最良の移植片候補であるか決定しなければならない。そのためには、蛍光が迅速に蓄積し、消散する様子を単に視覚的に観察するだけでは十分でない。例えば、連続注入によって残留ＩＣＧが組織内に蓄積し、注入ごとに背景輝度が徐々に上がるため、最良の穿通枝候補を視覚的に識別することがさらに容易でなくなるという傾向がある。さらに、ＩＣＧは、数分にわたって非常にゆっくり移動する場合があり、そのため、かかるオンザフライ分析は非常に困難となり、主観的となってしまう。外科医は、以下の事項によって査定を行うことになる。
１）ＩＣＧ結合血液が組織中に存在する量。
２）ＩＣＧ結合血液が組織中に存在する時間。
３）ＩＣＧ結合血液が組織中を移動する速度。
４）ボーラス注入後、解剖学的領域が発光する順序。

上記事項は、主観に基づいて答えることは困難である。したがって、穿通枝の位置を特定し、評価する客観的な標準を適用するための、より高度な画像処理および表示方法が求められている。

上述したように、本発明は、穿通枝皮弁の位置を特定し、分析する方法および装置を提供することにある。

本発明の一態様によれば、形成外科および／または再建手術のための、ＩＣＧ蛍光血管造影撮像を用いて穿通枝血管（ｐｅｒｆｏｒａｔｏｒ ｖｅｓｓｅｌ）を術前に識別する方法が開示される。本方法は、穿通枝位置をハイライトし、様々なメトリクス（ｍｅｔｒｉｃｓ）を計算することによって、穿通枝候補の視覚的識別を可能にする時間分解画像処理を含む。外科医は、以下の処理行為の少なくとも１つに従って、時系列分析を行い、メトリクスを出力するアルゴリズムの結果を選択し、比較することが可能となる。

ピクセルバイピクセル（ｐｉｘｅｌ−ｂｙ−ｐｉｘｅｌ）ベースで時間積分蛍光を求めること。
時間積分蛍光を経過時間で除算して、平均蛍光を計算すること。
蛍光の増大／流出率を求めること。
ピーク蛍光に達する経過時間を求めること。

様々な画像処理ステップによって、全取得時間にわたって、または選択した一時的な下位範囲で観察される入力シーケンスの各ピクセルについて、画像ピクセルを独立に処理し、特有の数値メトリクスを計算する。したがって、各画像出力は、入力画像シーケンスのフレームと同じ寸法、すなわち同じピクセル数および配列を有する数値アレイとなる。したがって、処理した画像を、撮像領域にわたって計算したピクセル値の３次元表示、例えば等高線図として、または色分けした２次元画像もしくは起伏図として表示することができる。かかる画像表示によって、画像特徴の迅速な理解が促進され、また、画像上の領域間の比較、この場合は、皮膚下の穿通枝の位置の比較が容易になる。

本発明の上記およびその他の特徴および利点は、以下の本発明の詳細な説明からよりよく理解されよう。また、以下に示す図は、本発明のある例示的な実施形態を示し、同じ参照番号は同じ要素を示す。これら図示の実施形態は、本発明の例示として理解すべきであり、いかなる形でも本発明を限定するものではない。

ＩＣＧ蛍光を観察するカメラシステムの概略図である。 ピクセル値を時間で積分した、皮膚のある領域のＩＣＧ蛍光画像である。 ピクセル値を経過時間で逆に重み付けして時間で積分した、皮膚のある領域のＩＣＧ蛍光画像である。 ピクセル値を蛍光の増大率によって求めた、皮膚のある領域のＩＣＧ蛍光画像である。 ピクセル値を最大蛍光に達する経過時間によって求めた、皮膚のある領域のＩＣＧ蛍光画像である。 ピクセル値をピーク蛍光によって求めた、皮膚のある領域のＩＣＧ蛍光画像である。 可変コントラスト伝達関数（ｖａｒｉａｂｌｅ ｃｏｎｔｒａｓｔ ｔｒａｎｓｆｅｒ ｆｕｎｃｔｉｏｎ）で処理した蛍光画像のオーバーレイを示す図である。 別の可変コントラスト伝達関数で処理した蛍光画像のオーバーレイを示す図である。 別の可変コントラスト伝達関数で処理した蛍光画像のオーバーレイを示す図である。

本発明は、いかなる切開も行う前の、非侵襲的方法による穿通枝皮弁の穿通枝血管の位置の術前判定を対象とする。

図１は、ＩＣＧ蛍光撮像による、術中、特に術前の、皮膚越しの非侵襲的な組織灌流判定に用いる装置を概略的に示す。赤外光源、例えば、ピーク発光が約７８０〜８００ｎｍの、ＩＣＧの蛍光を励起する１つまたは複数のダイオードレーザ、またはＬＥＤが、ハウジング１内に配置されている。蛍光信号は、適切な近赤外感度を有するＣＣＤカメラ２によって検出され、かかるカメラは、いくつかの製造業者（Ｈｉｔａｃｈｉ、Ｈａｍａｍａｔｓｕなど）から市販されている。ＣＣＤカメラ２には、ビューファインダ８を設けることができるが、画像はまた、術中に外部モニタ上で見ることもでき、このモニタは、電子画像処理および評価システム１１の一部とすることができる。

発散ビームまたは走査ビームでよい光ビーム３が、ハウジング１から発して、関心領域４、すなわち適切な穿通枝血管を備えた皮弁が位置すると思われる領域を照明する。関心領域は、約１０ｃｍ×１０ｃｍでよいが、外科的要件、および利用可能な照度、およびカメラ感度に基づいて変えることができる。

励起光がカメラセンサに達するのを阻止し、一方で蛍光光は通過可能となるように、典型的にはフィルタ６がカメラレンズ７の前に配置される。フィルタ６は、約８１５ｎｍよりも大きい波長しか透過しないＮＩＲ長波パスフィルタ（カットフィルタ）、または好ましくは８３０から８４５ｎｍのピーク波長を透過し、約１０ｎｍから２５ｎｍ、すなわち励起波長帯以外の半値全幅（ＦＷＨＭ）透過窓を有する帯域フィルタでよい。カメラ２はまた、関心領域のカラー画像を取得するように設計することができ、したがって蛍光画像とカラー画像との実時間相関が可能である。

本発明の例では、図１に示す装置を用いて、手術前に穿通枝血管を識別／位置特定するものであり、これは、外科医が、再建術中に使用すべき最良の皮弁または皮弁帯を選択する一助となる。

他の術後用途では、この装置は、以下の用途で使用することができる。
吻合開通、ならびに動脈流および静脈流の確認。これは、潜在的に皮弁不全をなくす結果を向上させ得るものであり、皮弁不全は、動脈流が低下し、灌流（ｐｅｒｆｕｓｉｏｎ）が不十分となり、静脈還流もやはり低下して、鬱血が生じる結果として生じる可能性がある。

完全な組織灌流の視覚化および確認。皮弁生着には、皮弁全体と元の組織との微小血管灌流が重要となるためである。

本発明では、画像処理および表示技術によって穿通枝の位置が視覚化され、それによって、穿通枝候補の客観的な視覚的識別が容易になる。ＩＣＧを注入し、ＩＣＧ蛍光が灌流し流出する全サイクルを、画像形成装置によって取得する。画像取得後、全シーケンス、または画像のいくつかの一時的な下位範囲を画像処理アルゴリズムによって処理し、このアルゴリズムは外科医が選択することができる。

蛍光測定値の処理結果は、例えば、疑似カラー画像または等高線図として視覚化することができ、それによって適用したアルゴリズムに従った迅速な視覚的評価が可能となる。例えば、各ピクセルの蛍光強度は、青（「低温」スポットまたは低蛍光強度もしくは低蛍光率）から赤（「高温」スポットまたは高蛍光強度もしくは高蛍光率）までの様々なスペクトルカラーとして描画することができる。他のスペクトル関連事項も容易に適応される。出力は、元の解剖学的画像上に半透明で重ね合わせて示すことができる。こうすることによって、「高温」スポットと、下にある解剖学的構造とを視覚的に相関させることが可能となる。「高温」スポットの意味は、積分強度、重み付けありまたは重み付けなし、増大率または流出率など、使用するアルゴリズムによって変わる。

使用者には、「高温」から「低温」へのカラーマッピングについて対話式制御が与えられ、実時間でこのカラーマッピングを変化させて、各アルゴリズムの出力メトリクスのダイナミックレンジのより細かい、またはより粗い下位範囲を調べることができる。カラー窓を広げると、最も高温の領域がまずハイライトされ、次いでより低温の領域がハイライトされる。この種の調整は、取得したピクセルと、表示されている画像のピクセルとの明度またはコントラストのマッピングを変えることによって行うことができる。かかるマッピング機能は、標準の画像形成プログラムに含まれ得る。現時点で使用しているメトリクスに基づいたこうした窓掛け処理は、穿通枝間での識別の一助となり、知覚を高め、注入したＩＣＧの動きを外科医がよりよく理解できることになる。

本発明はまた、患者の皮膚上の異なる２つの位置からの２つのシーケンスの同時表示および評価をサポートする。これによって、画像形成システムの視野よりも広い間隔を置いて離れている皮弁候補間の比較が可能となる。

図２は、適切な穿通枝血管が識別されると思われる、患者の皮膚のある領域の画像を示している。各ピクセルは、画像シーケンスの露光時間にわたる蛍光強度の時間積分を表している。このモードは、典型的には、画像処理の「積分モード（ｉｎｔｅｇｒａｔｉｏｎ ｍｏｄｅ）」と呼ばれ、多くの画像処理プロセッサで、このモードが標準機能として搭載されている。実際には、画像シーケンスの各フレーム中に取得されたピクセル強度（ＣＣＤ内で収集された電荷）を、例えば画像処理プロセッサ内で、ピクセルバイピクセルベースで加算し、フレームの数で除算し、その後、その総和を例えば１から２５５（８ビット）の固定したダイナミックレンジに正規化することができる。その概念は、画像中のより明るいピクセルが、予め設定された期間にわたってＩＣＧ搬送血液がより大量に流れ込んでいる皮膚領域を表す、ということである。図２では、穿通枝血管２４が最も高い積分蛍光強度を示し、２６で示す別の穿通枝血管が、それよりも弱い蛍光強度を示している。

この画像の透過性は、画面右上の医師によるマーカ２２が、ＩＣＧ蛍光画像の透過性カラーオーバーレイを透かして見えるように設定されている点に留意されたい。

図３は、ピクセル値を経過時間で逆に重み付けして時間で積分した同じ皮膚領域のＩＣＧ蛍光画像を示している。この画像処理アルゴリズムは、先に説明した積分と類似しているが、各ピクセルの測定強度を直接加算するのではなく、測定強度値を、ＩＣＧ蛍光の観測開始後の経過時間でまず除算し、その後加算する。このようにすると、より初期の蛍光信号には、後から取得した蛍光信号よりも高い重要性が与えられる。「最も高温」のピクセルは、画像フレームのシーケンス中で他のピクセルよりも早く蛍光を発するピクセルであり、他のピクセルにはＩＣＧボーラスが後に到達する。図２と同じ穿通枝血管３４が識別され、ここでは別の血管３６は、かろうじて識別可能である。

図４は、やはり同じ皮膚領域のＩＣＧ蛍光画像を示し、この画像のピクセル値は、蛍光強度の増大率によって求めている。この画像処理アルゴリズムでは、画像の各ピクセルについて、ピクセル強度対経過時間の傾きが計算される。例えば、各ピクセルには、最も低い強度値（ベースライン）と、最も高い強度値（または、別の比較的高い強度値）とを割り当てることができる。各ピクセルについて、ピクセル強度がベースラインと交差する時間と、ピクセル強度が高強度値と交差する時間とが記録される。この情報から、画像処理アルゴリズムによって、画像中の各ピクセルについて増大率を計算すると、「より高温の」ピクセルはより大きい傾きを有し、すなわち、それらのピクセルは、「より低温の」ピクセルよりも速く高強度値に達することになる。したがって、画像処理アルゴリズムのこの実施形態は、ＩＣＧボーラスが穿通枝血管に達する速度をハイライトする。図４では透過性がオフになっており、したがって外科医の器具は本画像では見られない。

ここでは参照符号４４で示す、先に識別された穿通枝血管は、遙かに良好に画定されており、血管４６（先に２６および３６として示したもの）および別の血管４８も同様に良好に画定されている。

図５は、ピクセル値を最大蛍光に達する経過時間によって求めた、同じ皮膚領域のＩＣＧ蛍光画像を示している。時間による変化率を示す図４とは異なり、図５の画像処理アルゴリズムは、ピクセルがそれらの最大強度に達する時間を示し、「より高温の」ピクセルが、より低温のピクセルよりも素早くそれぞれのピーク蛍光強度に達している。したがって、このアルゴリズムは、穿通枝がそれらのピーク強度に達する順序で画像領域をハイライトしている。この画像では、先に識別された穿通枝血管２４、３４、４４がやはり明白に識別可能であり、血管５６および５８は、図４の血管４６および４８に対応している。

図６は、ピクセル値を各ピクセルのピーク蛍光値によって求めた、同じ皮膚領域のＩＣＧ蛍光画像を示している。より高い（「高温」）蛍光強度値６４は、より高いＩＣＧ濃度を示し得るか、または穿通枝血管が皮膚表面により近く位置するために示されることもあり、皮膚表面では、励起光の吸収／蛍光応答が低減する。図４および５では明白に見えていた血管６６、６８は、背景からかろうじて識別可能である。

例えば図２および６に示す画像は、上述の様々なアルゴリズムを用いて処理したピクセル値を表示ピクセル強度に１：１でマッピングする線形コントラスト伝達関数を用いて描画しているが、画像はまた、画像の視覚差を強調するように可変コントラスト伝達関数を用いて（等高線図または疑似カラーオーバーレイとして）描画することもできる。さらに、オーバーレイ画像には、以下でＡＣＲ（積算または時間積分強度比率（ａｃｃｕｍｕｌａｔｅｄ ｏｒ ｔｉｍｅ−ｉｎｔｅｇｒａｔｅｄ ｉｎｔｅｎｓｉｔｙ ｒａｔｉｏ））ラベルと称するラベルを付すことができ、こうしたラベルによって、解剖学的構造の２つ以上の領域間での量的比較が容易になる。

可変コントラスト伝達関数のダイナミックレンジおよび傾きを改変すると、画像の絶対ピクセル値が変わるため、ＡＣＲラベルによって、選択された任意のオーバーレイ技術（例えば積算／時間平均強度など）で測定した異なる画像領域における相対的な灌流を使用者が比較することが可能となる。

以下の手法を用いて、ＡＣＲラベル値を計算する。分かりやすいように、積算強度をオーバーレイ技術として選択するが、利用可能なオーバーレイ技術のいずれを用いても同じ手法を取ることができる。
１）画像シーケンス中の全画像の全ピクセルについて、積算強度を時間窓にわたって計算する。
２）積算強度を、選択したラベルの領域（例えば、５×５ピクセル正方行列）にわたって平均する。
３）平均強度を、全画像の積算強度の最大値に正規化する。
４）正規化した平均を、伝達関数の最大値が１００％を表すようにスケーリングする。

この手法に従うことによって、伝達関数の傾きが改変されても、異なる２つのＡＣＲラベルの相対比率は変わらないままとなる。図７および図８は、異なる２つのコントラスト関数について、上述のアルゴリズムのうちの１つを用いて処理した画像シーケンスからの蛍光画像（グレースケール画像の上部）と、シーケンスからの積算強度をカラー（低値の青から高値の赤）で描画した疑似カラーオーバーレイ画像とを示す。図７のピクセル値は、第１のコントラスト伝達関数を用い、２つの領域についてそれぞれ強度５２％および７２％として処理しており、ラベル付けした２つの領域間の比率は５２／７２＝．７２に対応する。図８の第２のオーバーレイ画像は、異なるコントラスト伝達関数を用いて処理した同じピクセル値を示し、２つの領域の強度は、ここではそれぞれ９９％および７１％とラベル付けしてある。しかし、それらの相対比率は、７１／９９＝０．７２で本質的に変わらないままとなっている。

使用者は、制御領域が１００％でラベル付けされるように伝達関数を改変することができ、次いで、他の全ての領域をこの制御領域と比較することができる。

図９は、オーバーレイが透過性であることを示し、ここでは、積算強度ピクセルは、伝達関数ランプ（ｔｒａｎｓｆｅｒ ｆｕｎｃｔｉｏｎ ｒａｍｐ）の底部がピクセル値横軸と交差する点よりも低い値を有する。さらに、この例では、画像領域の１２％（下部ウィンドウ右下の被写域数字（ｃｏｖｅｒａｇｅ ｎｕｍｂｅｒ））が、最大積算強度の５２％よりも高い積算強度を有することを示している。この図は、いくつかの領域がそれらの５２％等高線で囲まれている様子を示している。

上述の実施形態は、近赤外スペクトル領域での励起後にＩＣＧによって経皮的に発せられた蛍光信号を検出する。しかし、組織が光を透過するスペクトル領域で励起し、蛍光を発することができるその他の染料も使用することができることが、当業者には理解されよう。

動脈血流、すなわち穿通枝血管（複数可）への血液供給の例に関して本発明を説明してきたが、本方法はまた、ピーク強度からベースラインに下がる変化率を定量化し、表示することによって、静脈鬱血による移植片不全を検出することもできる。ここでは、灌流領域の静脈還流がハイライトされることになる。

本発明には、様々な改変形態および代替形態が許容されるが、その特定の例を図面に示し、本明細書で詳細に説明してきた。しかし、本発明は開示の特定の形態または方法に限られるものではなく、本発明は、添付の特許請求の範囲の趣旨および範囲に含まれる全ての改変形態、均等物、代替形態をも包含するものであることを理解されたい。

Claims (6)

1. 穿通枝血管の灌流を評価するための装置の作動制御方法であって、
   検出手段が、ＩＣＧボーラスが血流に加えられた組織からの蛍光応答を検出するステップと、
   取得する手段が、所定の時間にわたって前記蛍光応答の時間的な画像シーケンスを取得するステップと、
   選択する手段が、前記画像シーケンスを処理するためのアルゴリズムを選択するステップと、
   画像を処理する手段が、前記選択したアルゴリズムを用いて、前記画像のピクセル値の時間積分強度または強度の時間微分を算出するように、前記画像シーケンスを処理するステップと、
   表示する手段が、前記時間積分強度または強度の時間微分を、カラー画像として表示するステップと、
   適用する手段が、コントラスト伝達関数を前記時間積分強度または前記強度の時間微分に適用するステップとを含み、
   前記コントラスト伝達関数は、前記時間積分強度または前記強度の時間微分を、異なる灌流特性を異なる色で表すオーバーレイ画像に変換する線形関数または非線形関数を表すことを特徴とする方法。
2. 前記時間的な画像シーケンスは、各ピクセルの前記時間積分強度値を、ＩＣＧボーラス注入後の経過時間で除算することによって処理されることを特徴とする請求項１に記載の方法。
3. 前記コントラスト伝達関数は、異なる傾きの領域を有する非線形関数であり、前記傾きと、前記異なる傾き間の移行とを前記血管の前記灌流を評価する手順の間、臨床医によって選択することができることを特徴とする請求項１乃至２のいずれか一項に記載の方法。
4. 前記時間積分強度または前記強度の時間微分の数値を、前記オーバーレイ画像で表示するステップをさらに含むことを特徴とする請求項１乃至３のいずれか一項に記載の方法。
5. １つ以上の選択された画像領域について、時間積分強度比率を計算するステップをさらに含み、
   前記時間積分強度比率は、
   全ピクセルについて前記時間積分強度を計算するステップと、
   前記時間積分強度を、前記選択された画像領域にわたって平均するステップと、
   前記平均された時間積分強度を、全画像の前記時間積分強度の最大値に正規化するステップと、
   前記正規化された強度を、前記コントラスト伝達関数の前記最大値を用いてスケーリングするステップと、
   を用いて計算されることを特徴とする請求項１乃至４のいずれか一項に記載の方法。
6. 穿通枝血管の灌流を評価する装置であって、
   ＩＣＧボーラスが血流に加えられた組織からの蛍光応答を検出する手段と、
   所定の時間にわたって前記蛍光応答の時間的な画像シーケンスを取得する手段と、
   前記画像シーケンスを処理するためのアルゴリズムを選択する手段と、
   前記選択したアルゴリズムを用いて、前記画像のピクセル値の時間積分強度または強度の時間微分を算出するように、前記画像シーケンスを処理する手段と、
   前記時間積分強度または前記強度の時間微分を、カラー画像として表示する手段と、
   コントラスト伝達関数を前記時間積分強度または前記強度の時間微分に適用する手段とを備え、
   前記コントラスト伝達関数は、前記時間積分強度または前記強度の時間微分を、異なる灌流特性を異なる色で表すオーバーレイ画像に変換する線形関数または非線形関数を表すことを特徴とする装置。

Images