"組織の少なくとも1つの構造的変化に関連するデータを決定するためのシステムであって、所定の技術を実行する際に、a第1の応力レベルにおける前記組織に関する第1の情報を含む第1の干渉計測信号と、第2の応力レベルにおける前記組織に関する第2の情報を含む第2の干渉計測信号とを受信し、b比較情報を生成するために、前記第1の情報と前記第2の情報とを比較し、c前記少なくとも1つの構造的変化に関連する前記データを、前記比較情報と、i前記組織又は組織内の造影剤の少なくとも一方についての1つ以上の既知の特性又はii干渉計システムの特性のうちの少なくとも1つに関連する詳細情報と、の関数として決定するように構成された処理装置を備えるシステム。前記少なくとも1つの構造的変化が前記組織の歪みである請求項1に記載のシステム。第2の応力が第1の応力とは異なる請求項1に記載のシステム。前記詳細情報は、前記組織の速度分布、前記組織の機械的特性、組織型、又は前記組織の構造のうち、少なくとも1つを含む請求項3に記載のシステム。前記機械的特性は、圧縮的特性又は弾性的特性のうちの少なくとも1つである請求項4に記載のシステム。前記処理装置は、前記組織から得られるドップラ信号に基づいて前記組織の速度分布を決定するように更に構成される請求項4に記載のシステム。組織の少なくとも1つの構造的変化に関連するデータを決定するための方法であって、第1の応力レベルにおける前記組織に関する第1の情報を含む第1の干渉計測信号と、第2の応力レベルにおける前記組織に関する第2の情報を含む第2の干渉計測信号とを受信するステップと、比較情報を生成するために、前記第1の情報と前記第2の情報とを比較するステップと、前記少なくとも1つの構造的変化に関連する前記データを、前記比較情報と、i前記組織又は組織内の造影剤の少なくとも一方についての1つ以上の既知の特性又はii干渉計システムの特性のうちの少なくとも1つに関連する詳細情報と、の関数として決定するステップと、を有する方法。前記少なくとも1つの構造的変化が前記組織の歪みである請求項7に記載の方法。第2の応力が第1の応力とは異なる請求項7に記載の方法。前記詳細情報は、前記組織の速度分布、前記組織の圧縮的特性又は弾性的特性のうちの少なくとも1つ、組織型、又は前記組織の構造のうち、1つ以上を含む請求項7に記載の方法。前記組織から得られるドップラ信号に基づいて前記組織の速度分布を決定するステップを更に有する請求項8に記載の方法。組織の少なくとも1つの構造的変化に関連するデータを決定するためのソフトウェアであって、処理装置によって実行される際に、第1の応力レベルにおける前記組織に関する第1の情報を含む第1の干渉計測信号と、第2の応力レベルにおける前記組織に関する第2の情報を含む第2の干渉計測信号とを受信するように前記処理装置を設定する第1の命令群と、前記処理装置によって実行される際に、比較情報を生成するために、前記第1の情報と前記第2の情報とを比較するように前記処理装置を設定する第2の命令群と、前記処理装置によって実行される際に、前記少なくとも1つの構造的変化に関連する前記データを、前記比較情報と、i前記組織又は組織内の造影剤の少なくとも一方についての1つ以上の既知の特性又はii干渉計システムの特性のうちの少なくとも1つに関連する詳細情報と、の関数として決定するように、前記処理装置を設定する第3の命令群と、を有するソフトウェア。前記少なくとも1つの構造的変化が前記組織の歪みである請求項12に記載のソフトウェア。第2の応力が第1の応力とは異なる請求項13に記載のソフトウェア。前記詳細情報は、前記組織の速度分布、前記組織の圧縮的特性、組織型、又は前記組織の構造のうち、少なくとも1つを含む請求項12に記載のソフトウェア。前記処理装置によって実行される際に、前記組織から得られるドップラ信号に基づいて前記組織の速度分布を決定するように、前記処理装置を設定する第1の命令群を更に有する請求項15に記載のソフトウェア。組織の少なくとも1つの弾性率に関連するデータを決定するためのシステムであって、所定の技術を実行する際に、a前記組織に関する情報を含む少なくとも1つの干渉計測信号を受信し、b前記組織の前記弾性率を、前記受信した少なくとも1つの干渉計測信号の関数として決定するように構成された処理装置を備えるシステム。前記情報は、前記組織の構造及び前記組織の組成のうち、少なくとも1つを含む請求項17に記載のシステム。前記少なくとも1つの干渉計測信号は、第1の応力レベルにおける前記組織に関する第1の情報を含む第1の干渉計測信号と、第2の応力レベルにおける前記組織に関する第2の情報を含む第2の干渉計測信号とを含み、第2の応力が第1の応力とは異なる請求項17に記載のシステム。前記処理装置は、比較情報を生成するために、前記第1の情報と前記第2の情報とを比較するように更に構成され、前記弾性率が比較情報の関数として決定される請求項19に記載のシステム。前記処理装置は、前記第1の情報又は前記第2の情報の1つ以上の関数として、数値モデルを生成するように更に構成される請求項19に記載のシステム。前記数値モデルが動的な数値モデルである請求項21に記載のシステム。前記動的な数値モデルは、前記第1の情報又は前記第2の情報の1つ以上の関数として修正可能な、制限、モデルの複雑さ又はモデルの次数のうちの少なくとも1つを含む請求項22に記載のシステム。前記モデルの複雑さ又はモデルの次数は、前記第1の情報又は前記第2の情報の少なくとも1つの関数として修正可能である請求項23に記載のシステム。前記動的な数値モデルは、詳細情報を生成するために実行され、前記詳細情報は、前記制限、前記モデルの複雑さ又は前記のモデル次数のうちの1つ以上を修正するために、前記動的な数値モデルに対して提供される請求項23に記載のシステム。前記モデルの複雑さは、複数のモデル要素を有し、該モデル要素のうちの少なくとも第1の要素は、前記モデル要素のうちの少なくとも第2の要素と関連し、該関連は前記モデル要素の前記第1の要素及び前記第2の要素の重みに基づいている請求項23に記載のシステム。前記処理装置は、前記数値モデルを使用して前記組織に関する詳細情報を生成するように更に構成され、前記詳細情報は、前記組織に加える応力に対する前記組織の応答に関連している請求項21に記載のシステム。前記処理装置は、比較情報と前記詳細情報の関数として、更に別のデータを生成するように更に構成される請求項27に記載のシステム。前記処理装置は、前記更に別のデータの関数として、前記数値モデルを修正するように更に構成される請求項28に記載のシステム。前記弾性率が前記数値モデルに基づいて決定される請求項29に記載のシステム。前記処理装置は、前記数値モデルに基づいて前記組織の歪み情報を決定するように更に構成される請求項29に記載のシステム。前記比較情報は更に、i前記組織又は組織内の造影剤の少なくとも一方についての1つ以上の既知の特性又はii干渉計システムの特性のうち、1つ以上に関連する詳細情報に依存する請求項20に記載のシステム。前記詳細情報は、前記組織の速度分布、前記組織の圧縮的特性又は前記組織の弾性的特性の少なくとも1つ、組織型、及び前記組織の構造のうち、1つ以上を含む請求項32に記載のシステム。前記処理装置は、前記組織から得られるドップラ信号に基づいて前記組織の前記速度分布を決定するように更に構成される請求項33に記載のシステム。組織の少なくとも1つの弾性率に関連するデータを決定するための方法であって、a前記組織に関する情報を含む少なくとも1つの干渉計測信号を受信するステップと、b前記組織の弾性率を、前記受信した少なくとも1つの干渉計測信号の関数として決定するステップと、を有する方法。前記情報は、前記組織の構造及び前記組織の組成のうち、少なくとも1つを含む請求項35に記載の方法。前記少なくとも1つの干渉計測信号は、第1の応力レベルにおける前記組織に関する第1の情報を含む第1の干渉計測信号と、第2の応力レベルにおける前記組織に関する第2の情報を含む第2の干渉計測信号と、を含み、第2の応力が第1の応力とは異なる請求項35に記載の方法。比較情報を生成するために、前記第1の情報と前記第2の情報とを比較するステップを更に有し、前記弾性率が比較情報の関数として決定される請求項37に記載の方法。前記第1の情報又は前記第2の情報の1つ以上の関数として、数値モデルを生成するステップを更に有する請求項37に記載の方法。前記数値モデルを使用して前記組織に関する詳細情報を生成するステップを更に有し、前記詳細情報は、前記組織に加える応力に対する前記組織の応答に関連している請求項39に記載の方法。比較情報と前記詳細情報の関数として、更に別のデータを生成するステップを更に有する請求項40に記載の方法。前記更に別のデータの関数として、前記数値モデルを修正するステップを更に有する請求項41に記載の方法。前記弾性率が前記数値モデルに基づいて決定される請求項42に記載の方法。前記数値モデルに基づいて前記組織の歪み情報を決定するステップを更に有する請求項42に記載の方法。前記比較情報は更に、i前記組織又は前記組織内の造影剤の少なくとも一方についての1つ以上の既知の特性及びii光コヒーレンス断層撮影システムの特性のうち、1つ以上の詳細情報に依存する請求項38に記載の方法。前記詳細情報は、前記組織の速度分布、前記組織の圧縮的特性又は前記組織の弾性的特性の少なくとも1つ、組織型、及び前記組織の構造のうち、1つ以上を含む請求項45に記載の方法。前記組織から得られるドップラ信号に基づいて前記組織の速度分布を決定するステップを更に有する請求項46に記載の方法。組織の少なくとも1つの弾性率に関連するデータを決定するためのソフトウェアであって、処理装置によって実行される際に、前記組織に関する情報を含む少なくとも1つの干渉計測信号を受信するように前記処理装置を設定する第1の命令群と、前記処理装置によって実行される際に、前記組織の弾性率を、前記受信した少なくとも1つの干渉計測信号の関数として決定するように前記処理装置を設定する第2の命令群と、を有するソフトウェア。しかしながら、試料の構造、組成及び生体力学的特性を同時に決定する方法が利用可能でない。このような能力は、生物医学において広い用途をもつが、特に、脆弱なプラークを見つけだして、急性心筋梗塞との関連性を理解する上で有効となろう。更に、弾性率計測法及び弾性率画像化法は一般に、多数の独立した有限要素、つまり対象領域で空間的に分布している画像ピクセルに亘る未知の歪み又は弾性率パラメータの推定を使用する。歪み又は弾性率の画像化に使用する空間分解能が高くなるほど、推定すべき独立の未知数の数が増える。パラメータ空間が大きくなるにつれて、所望の目的関数を満たすパラメータの推定についての探索が、困難な劣決定問題underdeterminedproblemとなってしまう虞がある。通常は、未知数の数が、基礎となる画像化データから一意に決定することができる数を大幅に超えると、目的関数を満たす多くの解が得られることになる。また、高次元100次元を超えるオーダーのパラメータ空間を探索するためには、一般に高い計算コストと長い計算時間がかかる。生物組織の弾性率計測法及び弾性率画像化のための従来の方法は、有限要素毎に又は対象のピクセルでの歪み又は弾性率を独立の未知数として処理する。その際、通常は、上記非特許文献6、非特許文献7に記載されるように、目的関数の最適化のためのLevenburg-Marquardt又は同様のアルゴリズムを使用する。極小値に対するロバスト性を実現するために、多重解像度法が用いられ、この方法では、少ない未知数を用いて低解像度のグリッド格子で推定値が得られ、次に、最大解像度の領域におけるパラメータ最適化の初期化のために、これらの低解像度の推定値が高解像度のグリッドにマッピングされる。これらの従来の方法は時間がかかり、対象領域が大きいと数分の処理を要する場合がある。従来の方法とは対照的に、本発明によるシステム、プロセス、及びソフトウェア構成の例示的な実施形態は、少なくとも1つの試料における歪み及び弾性率の空間分布を、高い空間分解能及び高感度で決定することができるとともに、構造及び組成の高解像度画像を同時に提供できる。本発明のシステム、プロセス、及びソフトウェア構成は広く適用可能であり、その能力は特に生物組織及び血管組織に関連している。本発明の1つの例示的な実施形態では、血管の構造及び組織の組成を決定するためにOCTが使用される。そして、この情報が、血管を表現する数値モデルを構築するために使用され、その後、任意の応力負荷に対する血管の機械的応答を予測するために、弾性率の推定を用いた有限要素モデリングが使用される。この例示的な計算とは別に、本発明による例示的なOCT弾性率計測法は、血管の機械的応答を測定するために実行される。モデリングと画像化という2つの経路は、a想定された弾性率分布に基づく予測と、b測定を、それぞれに表している。この2つの結果の差については、弾性率分布の初期推定を修正することによって最小化される誤差関数とみなすことができる。本発明によるこの例示的な方法を反復することによって、弾性率の分布及び大きさを決定することができる。この情報は、弾性率の断面画像又は三次元画像として表示される。更に、誤差関数を最小化することによって、改善された、歪みの弾性率計測法による画像を生成することができる。その結果、本発明によるシステム、プロセス、及びソフトウェア構成の例示的な実施形態では、現在の診断技術の制約、すなわち、構造及び/又は歪みは測定されるが、組織の生体力学的特性が不明であるといった制約を解消することができる。更に、本発明は、弾性率計測法のための従来法の解像度及び感度を改善する。要約すると、本発明によるシステム、プロセス、及びソフトウェア構成の例示的な実施形態によって、医学的用途や非医学的な用途のために、試料の構造、組成、歪み及び弾性率を同時に決定できるようになる。本発明の1つの例示的な実施形態では、組織の少なくとも1つの構造的変化に関連するデータを決定するためのシステム、プロセス及びソフトウェア構成が提供される。より詳細には、第1の応力レベルでの組織に関する第1の情報を含む第1の光コヒーレンス断層撮影法OCT信号と、第2の応力レベルでの組織に関する第2の情報を含む第2のOCT信号と、が受信される。第1の情報と第2の情報は、比較情報を生成するために比較される。少なくとも1つの構造的変化に関連するデータが、比較情報と、i組織についての1つ以上の既知特性及び/又はiiOCTシステムの特性に関連する詳細情報との関数として決定される。例えば、前記構造的変化は、組織の歪みである。また、第2の応力は第1の応力と異なる。前記詳細情報には、組織の速度分布、組織の機械的特性例えば、圧縮可能性compressability及び/又は弾性的性質、組織型、組織中の造影剤の光学特性、及び/又は組織の構造がある。更に、組織の速度分布は、組織から得られるドップラ信号に基づいて決定される。本発明の別の例示的な実施形態によれば、組織の少なくとも1つの弾性率に関連するデータを決定するための方法、システム、及びソフトウェア構成が提供される。例えば、組織に関する情報を含む、少なくとも1つの光コヒーレンス断層撮影法OCT信号が受信される。そして、組織の弾性率が、受信された1つ以上のOCT信号の関数として決定される。例えば、この情報には、組織の構造及び/又は組織の組成を含むことができる。OCT信号は、第1の応力レベルでの組織に関する第1の情報を含む第1のOCT信号と、第2の応力レベルでの組織に関する第2の情報を含む第2のOCT信号とを含むことができ、その際、第2の応力が第1の応力とは異なってもよい。第1の情報と第2の情報については、比較情報を生成するために比較することができ、弾性率は比較情報の関数として決定される。また、数値モデルが、第1の情報又は第2の情報のうちの1つ以上の関数として生成される。組織に関する詳細情報を、数値モデルを使用して生成することができ、この詳細情報は、組織に加える応力に対する組織の応答に関連している。前記数値モデルは動的な数値モデルとすることが可能であり、この数値モデルは、i制限、iiモデルの複雑さ、及び/又はiiiモデルの次数を含むことができ、これらについては、第1の情報及び/又は第2の情報の関数として修正可能である。モデルの複雑さ及び/又はモデルの次数は、第1の情報及び/又は第2の情報の関数として修正可能である。動的な数値モデルは詳細情報を生成するために実行され、この詳細情報は、制限、モデルの複雑さ及び/又はモデルの次数を修正するために動的な数値モデルに与えられる。モデルの複雑さには、複数のモデル要素が含まれ、該モデル要素のうち、少なくとも第1の要素は、該モデル要素における第1の要素及び/又は第2の要素の重みに基づいて、複数の要素と関連をもつ。また、前記比較情報及び詳細情報の関数として、更に別のデータが生成される。数値モデルは、この更なるデータの関数として修正される。そして、弾性率が数値モデルに基づいて決定される。組織の歪み情報は数値モデルに基づいて得られる。前記比較情報は更に、詳細情報、つまり、i組織の少なくとも1つの既知の特性、及び/又は、iiOCTシステムの特性についての情報に依存する。この詳細情報には、組織の速度分布、組織の圧縮性/弾性的性質、組織型、組織内の造影剤の光学特性、及び/又は組織の構造が含まれる。組織の速度分布は、組織から得られるドップラ信号に基づいて決定される。更に、生体力学的画像化のための従来の方法及び技術とは対照的に、本発明の別の例示的な実施形態では、同じ型の組織が同様の、そしておそらくはほぼ同一の機械的な特性を有することや、対象領域を、区別できる組織構成要素例えば、繊維性、脂質性、石灰化などに区分化するために、組織の高コントラストがOCT法に利用可能であることを考慮している。本発明によるプロセスの、このような例示的な実施形態は、組織間に存在する境界を保持したままで、パラメータの探索空間を低減させることができる。例えば、低解像度グリッドでの推定とは異なり、この技術によって、各要素中の組織型の部分的なボリューミングが最小化され、これにより歪み又は弾性率の急峻な空間的勾配が保たれる。また、生体力学モデルがデータに十分良好に適合していない要素の適応的メッシュの精緻化は、弾性率計測法及び弾性率画像化法に対して有益に寄与することとなる。本発明の上記の目的及び他の目的、特徴及び利点は、添付の特許請求の範囲と併せて、以下に記載する本発明の詳細な説明を読む際に、更に明白となるであろう。本発明の更に別の目的、特徴及び利点については、本発明の例示的な実施形態を示す図面を参照して、以下の発明の詳細な説明から明らとなるであろう。本発明の特定の例示的な実施形態では、試料の構造、組成、歪み及び/又は弾性率を決定するために、OCT法を有限要素モデリング法と組み合わせたハイブリッド技術を利用することができる。図1は、正常な血管壁4内に埋設された脂質プール3を有する、罹患動脈の断面を例示的に示す。内腔5内の血圧変動によって、血管及びプラーク形状の変形が生じる。例えば、図1では、腔内の圧力が上昇すると、点線の外形線1が実線の外形線2の位置へと変形する。ここに記載するOCT弾性率計測法の例示的な実施形態では、血管及びプラークの境界の変位を追跡することはもとより、組織自体内で発生する生体力学的な歪みをも推定することができる。本発明によるOCT弾性率計測法の1つの例示的な実施形態によれば、試料内の散乱体の速度例えば大きさ及び方向を決定することができる。また、OCTの構造的な画像と組織の弾性率推定に基づいて、ある例示的な有限要素モデリング法を用いることで、対応する速度分布を予測することが可能になる。これら2つの例示的な速度分布間の差は、弾性率の初期推定について最適化を繰り返すことによって最小化される誤差関数とみなせる。そして、結果的に得られた弾性率分布が、画像形態で視覚化される。更に、歪みを画像形態で視覚的に表現するために、OCT弾性率計測法のデータを使用することができる。図2は、本発明に従い、弾性率分布を決定するために、OCT弾性率計測法と有限要素の例示的なモデリング法とを組み合わせて、別々のステップとして図示的に示す一般的なフローチャートである。処理の最初には、例えば、動脈が心臓周期に亘って動的な変形を受ける際に、OCTの画像化及び取得がステップ50で行われる。同時に、ステップ55において、腔内の圧力がデジタル化されるか、対応するOCTフレームで記録されるか、この両方が行われる。ステップ65において、1つの圧力レベルから取得されたOCT画像が、罹患血管のための幾何形状モデルの基礎となり、これは、ステップ70において、有限要素モデリングFEMによる数値シミュレーションのためにメッシュ化される。ステップ60において組織の速度場を得るために、時間の関数としてのOCT画像データの変化が、運動推定のための例示的な方法を用いて追跡される。次に、ステップ90において、測定された速度場から、これに対応する歪みの固有値及び固有ベクトルが計算され又は決定される。ステップ95において、得られた組織の歪みの画像が、画像表示される。また、ステップ60で生成された、推定済みの組織の速度は、モデルに基づいた弾性率70の決定に対して、その基礎となる。ステップ75では、測定された圧力負荷ステップ55及び既知の境界条件を、弾性率の値についての既定の分布でもって初期化した有限要素メッシュに適用することによって、数値モデルが実行され、予測の速度分布が得られる。ステップ80において、2乗誤差測定squared-error-measure法を使用することによって、ステップ75で得られたモデル予測による速度が、測定された速度60と比較される。ステップ80での比較によって生じる誤差に基づいて、弾性率の値が更新されて、ステップ65にてモデルが再構築され、次にステップ70において再シミュレーションが行われて、新しい予測速度の組が得られる。このプロセスは、弾性率の推定値が指定された許容レベルに収束するまで続けられる。収束の後に、ステップ85では、最終的に得られた例示的な弾性率分布が画像表示される。図3は、本発明による例示的な方法、つまり、例示的なOCT法を使用して弾性率分布を決定するために速度分布を使用する方法についてフローチャートを示す。より詳細には、ステップ100のOCT画像取得とステップ105の腔内圧力の記録が同時に実行される。本発明のこの例示的な実施形態では、ステップ110において、多重解像度変分法を使用して、動的なOCTデータセットが処理され、その結果は、2つの画像化時点間の組織変位のロバストな推定となる。組織歪みの固有値と固有ベクトルは、ステップ140において速度推定から決定され、そして、ステップ145において視覚的に画像表示される。例示的な弾性率推定では、ステップ115において、血管及びプラークの表面を抽出するために、参照時点でのOCTデータが区分化される。また、これらの表面については、動脈特有の形状を画定するために、三次元で再構築し得る。そして、この血管形状がメッシュ化され、境界条件が適用されて、ステップ120においてメッシュ要素が初期の弾性率の値を割り当てられるが、これは、ステップ125の有限要素モデリング法で更に使用するためである。この例示的な処理/方法から、2乗誤差測定値を決定するために、測定された組織速度と共に使用される組織の速度予測の組がステップ130で得られ、この2乗誤差測定は、数値モデルで使用される弾性率値及び境界条件の更新を推し進める。弾性率の更新と数値シミュレーションの例示的な方法は、2乗誤差測定値を最小化して、指定の許容レベルに収束し得る弾性率の推定値を生成するために、反復的に続行される。最終的な弾性率分布は、ステップ135において視覚的に画像表示される。本発明の別の例示的な実施形態では、加える応力を変えて測定しながら、試料の構造及び組成を決定するためにOCT法を実行することができる。OCT画像取得レートについては、個々の画像内における運動への顕著な人為的影響を回避できるように十分大きい。加えた応力の1つの値について決定される構造及び組成は、組織を表わす数値モデルを生成するために使用され、そして、数値モデリングは、測定された応力変化及び弾性率分布の初期推定を組み込むことで、第2の応力に対する構造を予測するために使用される。対応する第2の応力で取得されたOCT画像は予測された構造と比較され、予測された構造と測定された構造との差が、弾性率分布の初期推定を反復的に最適化することによって最小化される。この例示的な実施形態では、数値モデリング、例えば、最適化された弾性率分布に基づくモデリングが、弾性率分布と歪みを統合した最終的な決定のために使用される。この結果が、画像形態で視覚的に表示される。本発明の更に別の例示的な実施形態では、試料の構造及び組成を決定するためにOCT法が実行され、この情報からは、組織を表現する数値モデルが生成される。数値モデリングは、試料の弾性率分布の推定に基づき、速度分布を予測するために使用され、該速度分布は、加えた応力に対する応答として試料内に生じるものである。そして、試料中を移動する散乱体からのOCTビームの反射によって生じるドップラ周波数偏移は、画像強度データに加えて、試料中の深さ分解速度分布を決定するために使用される。モデルによる速度予測値と、OCTドップラデータ及び画像強度データからの速度測定値との差は、初期弾性率分布を反復的に最適化することによって最小化される。そして、得られた弾性率分布が、画像形態で可視化される。これとは独立に、試料中の歪みを視覚的に表現するために、ドップラOCTデータを使用することができる。図4は、本発明による別の例示的な方法のフローチャートを示し、この方法では弾性率分布及び歪み分布を決定するために構造及び構造的変形を使用する。ステップ200において、例示的なOCT画像化法が実行され、同時にステップ205で腔内圧力の記録が実行される。得られたデータは、ステップ220における画像及びドップラ情報を有する探索データセットと、ステップ210における数値モデル構築のために処理される参照データセットと、に分けられる。そして、ステップ215では、参照データセットと探索データセットとの間で、有限要素法シミュレーションとモデルの剛体運動推定との統合のために、参照形状と圧力負荷が使用される。ステップ225では、推定された剛体モデルの変換が、参照強度データの再サンプリングを行うために、モデル予測によるメッシュ変形例えば、参照強度データを、参照する探索データセットフレームへと効果的に歪ませることと組み合わされる。ステップ230では、歪められた参照データと測定された探索データが、OCTで指定される目的関数内で結合される。未知の弾性率値は、目的関数を最大化するために、モデル構築ステップ210にて反復的に更新される。弾性率の推定値が収束すると、ステップ235、240において、数値モデルからの対応する弾性率及び歪みの分布が出力され、視覚的に画像表示される。光コヒーレンス弾性率計測法光コヒーレンス弾性率計測法は、好ましくは、超音波弾性率計測法の基礎をなす原理と同じ原理に基づく。例えば、組織が機械的負荷の下で画像化されると、その変位は、巨視的な構造例えば組織界面に対応する画像特徴に現れる。また、運動はコヒーレント画像のスペックルに現れるが、これは微視的な組織散乱体の空間分布が負荷の下で変化するためである。巨視的構造及び微視的スペックルからの運動推定では、連続した画像同士の間で画像特徴が良好に保たれると仮定している。このため、所望の速度推定は、異なる負荷条件下で取得した、参照画像内のブロックと、探索画像内のブロックとの間の類似性測度を最大化し得る。OCTで得られる干渉画像は、ビーム減衰をモデル化する指数関数的な減衰項と、空間的な畳み込みとの積によって近似することができる。座標xは側方の走査方向に対応し、そして座標yは軸上の走査方向に対応しており、前者はサンプルビーム軸に直交する方向を向き、後者はサンプルビーム軸に対して平行な方向を向いている。上付きバー記号を付したパラメータμsは、試料の散乱による平均減衰量であり、σbx,yは、様々な後方散乱断面積を有する、点の分布として試料内の後方散乱をモデル化しており、hx,yは、OCTシステムの点広がり関数PSFを表している。OCTのPSFは、ソース検出光学系の自己相関関数Γyと、瞳関数pxとの間での、分離可能な空間的不変積として近似できる。ガウス型ビームの場合、以下のようになる。ここで、ESは光源の電界ベクトルの振幅、0は光源の自由空間中心波長、Δは光源のFWHM半値全幅でのスペクトルバンド幅、fは対物レンズの焦点距離、Dは対物レンズの入射瞳におけるサンプルビームの1/e2強度直径である。この画像形成モデルに基づいて、参照画像における位置x0,y0から探索画像内の新位置x0u,y0vへと変位を受ける1つの点状散乱体は、これに対応する以下の参照及び探索の干渉像を有することになる。従来の速度測定を用いて、組織の運動は、式5及び6のエンベロープ又は複素振幅のサブブロック間での相関係数を最大化することによって概ね推定される。各画像は、事前に決められたサイズのブロックへと更に分割できる。各参照ブロックについて、全ての探索画像ブロックとの相互相関が計算されて、相関係数が、参照位置と探索位置との間の相対変位の関数として得られる。例えば、本発明の1つの例示的な実施形態によれば、探索画像内で最もマッチングするブロックは、正規化された相互相関関数を最大にし、このブロックと参照ブロックとの間の相対的なオフセットが所望の速度推定を与える。この手順は、参照位置x,y、平均強度μR及びμSを有するM×N個のサブブロックに関して、式7及び8で数学的に表現される。重複するサブブロックは、参照画像内のより細密な格子で速度を推定するために使用される。平行移動するインパルス散乱体の場合に、式5~8による速度推定は、のときに、インパルス応答のシフトを追跡することができる。非剛性の変形を受ける散乱体の集団を含む現実の組織の場合には、各散乱体からのコヒーレントインパルス応答が後方散乱信号の干渉パターンを生成し、この信号は、時間的に連続する画像間で単純に平行移動しそうにない。速度推定は、干渉ノイズ、つまり式8で最大に達成可能とされる相関係数を低減させるノイズに敏感となる。上記のような非相関効果は、相関ウィンドウ窓のサイズが対象となる変形構造と比べて広いときは常に生じ、また負荷によって生起される歪みが大きい場合に生じることがある。いずれの場合にも、機械的負荷の影響は、単純なスペックルの平行移動ではモデル化できず、その理由は、空間歪みが根本的な散乱体分布で生じるためである。特定の現実的な状況では、画像ノイズと非相関性は、相関表面x,yu,v内の真の変位における相関値を低減させるだけではなく、真の変位における相関値を超える値をもった複数の極大値又は偽ピークに加えて、相関ピークの位置をシフトさせるジッタをもたらす。Aライン間の一次元相関は、理想的には真の変位において単一の明確なピークとなるべきである。しかし、スペックル非相関性とノイズのために、一般に複数のピークが相関関数に存在し、そのうちの最も高いピークは真の変位よりも非常に小さい速度に存在する。また、二次元の運動推定の場合には、理想的な相関関数が、単一の明確なピークを示さねばならないが、複数の極大値が存在することがある。相関ウィンドウ全体に及ぶ境界等の特徴を有する画像からの二次元推定の場合には、境界の接線方向の速度成分を決定することが難しくなる。このような場合は、相関関数は明確なピークを含まず、相関値は、境界の接線方向を向いて幅広い範囲の変位に亘って高くなる。一次元及び二次元のいずれの場合にも、得られる速度推定は、血管OCEでの使用にとって極度にノイズが多い歪み推定をもたらす虞がある。ロバストなコヒーレンス弾性率計測法速度推定を改善するための方策には、ノイズ抑圧のための画像シーケンス列のブラーリング、広い相関ウィンドウの採用、及び相関最大化による推定後の速度場のスムージングがある。これらの方策は、ある観察に基づくと、速度及び歪みの推定をある程度改善するが、弾性率計測法に対するOCTの空間分解能の利点を損なうこともある。例えば、画像シーケンスのブラーリングは、ノイズだけでなく、運動追跡に有用とされる精緻な画像特徴をも除去する場合がある。広い相関ウィンドウは、速度場の微細な変化を追跡する能力を低下させる虞があり、また、式7及び8に仮定した並進スペックルモデルの破綻につながる可能性がある。メジアンフィルタ又は他のスムージングカーネルを用いた速度又は歪みのフィルタリングは、測定が行われた後の測定値に作用する。このため、このような方法は、速度及び歪みの推定を改善するために、基礎をなす相関関数内に存在する情報を使用できない可能性がある。より好ましい推定法は、相関を最大化する手順自体の間に、データ駆動型の速度フィルタリングを可能にすることである。そのような例示的な方法は、下記に記載する変分法である。速度推定問題は、変分エネルギーの最小化として提示され、これは、拍動性動脈壁内の速度場に関する従来の知識を用いて、推定にロバスト性を与える一方で、相関関数に存在する速度情報を有効に利用することを目的とする。この手法では、最大解像度データからの利用可能な全情報を維持するために、画像のスムージングを使わない。全体的な変分エネルギー関数は、以下のようになる。このエネルギーは、未知の速度場V=[u,v]に依存し、データ忠実度EDV、歪み場の平滑度ESV、及び動脈壁の非圧縮性EIVを調整する3つの項の重み付けされた組み合わせである。これら各項の関数形は、以下の通りである。ここで、式10の表式x,yVは、式8に示した相関係数関数である。歪み平滑度項及び組織非圧縮性項が存在しない場合において、データ忠実度項の最小化は、相関関数の最大化と同じであり、これは従来の速度測定からの推定と同一の速度推定をもたらす。歪み平滑度項及び組織非圧縮性項は、動脈の組織生体力学に関する従来の知識から逸脱しないように、速度推定を制限する。隣接する参照位置からの相関関数の情報は、非相関性、偽ピーク、及び高い相関係数値の明確に定義されない領域について、ロバスト性を与えるために、効果的に組み合わされる。歪み平滑度項は、動脈の速度場の二次導関数を、壁に亘って滑らかに変化させるのに対し、非圧縮性モデルは、u及びvの速度場の挙動と結び付くことで、壁内部の点が非圧縮性から大きく逸脱しないようにする。所望の速度場推定は、全体的な変分エネルギーを最小化する。エネルギー最小化問題に対する数値解を得るためには、式13の連続的表現を以下のように離散化する。ここで、列x方向及び行y方向の離散的な速度成分は、順序付けられた以下の縦ベクトルとしてそれぞれに表現される。ここで、kは、対象のk番目の参照位置についての、順序付けの添え字であり、[ik,jk]は、参照画像行列IR[i,j]内のこの位置での行座標及び列座標である。離散化されたデータ忠実度項は、以下の式で表される。ここで、離散的な領域における相関係数関数は、参照画像行列をIR[i,j]、探索画像行列をIS[i,j]とした場合に、以下の通りである。参照画像行列IR[i,j]及び探索画像行列IS[i,j]については、直線的な格子上で定義された、一定間隔を置いた点集合からサンプリングされる。実際には、全画像について秒オーダーの迅速な計算を実現するために、二次元FFTを使用する高速な正規化された相互相関近似を使用して、式17の分子と、式17の分母についての事前に計算される移動和を計算する。離散化された歪み平滑度項及び非圧縮性項は、それぞれ下式のようになる。ここで、D2r及びD2cは、列ベクトルu及びvにおける近隣位置からの速度に作用する第二次の行差分行列及び列差分行列である。行列D1r及びD1cは、これらに対応する第一次の行差分演算子及び列差分演算子である。M×Nの矩形状領域で定義された二次元の速度場から生成される、順序付けられた速度ベクトルの場合に、第一次の行差演算子D1r及び第一次の列差演算子D1cは、以下のように定義される。ここで、D1M-1NはM-1×Nの第一次の差分行列であり、IMはM×Mの単位行列である。これに対応する第二次の行差分演算子及び列差分演算子は、それぞれ以下のようになる。ここで、D2M-2NはM-2×Nの第二次の差分行列であり、IMはM×Mの単位行列である。式14を最小化するために、第一次の変分を導出して、オイラー方程式を得る。ここで、データ忠実度項の第一次変分は以下のように定義される。式24のオイラー方程式は、以下の発展方程式を作ることで反復的に解くことができる。ここで、Aは以下の通りである。tは、反復毎にとった時間ステップである。時間tにおける更新された速度推定を解くためにこれを整理すると、以下の行列及びベクトルの方程式を得る。定常状態では時間導関数が消え、これにより得られる速度推定は式24を満たす。実際には、速度場の初期推定から開始して、式27の各反復においてLU分解法を使って、更新された速度推定を解くことが可能である。非整数の速度推定では、双三次補間を用いて、データ忠実度項で必要な勾配を計算する。この例示的な方法は、速度場の大きさについての最大変化が0.01未満となるまで続行される。大域的最小値への収束のための多重解像度推定式27の解は、変分エネルギー関数の極小値へと収束することとなる。よって、良好な大域的収束性を保証するために、未知の速度場は大域的最小値の近くで初期化する必要がある。これを実現するために、図5のブロック図として示した、本発明による例示的な多重解像度技術が使用可能である。例えば、入力される最大解像度の参照画像及び探索画像は先ず、ステップ255において10分の1にダウンサンプリングされて、低解像度シーケンスが得られ、ステップ260では、それから、式7、8の相関最大化によって、速度場の低解像度初期推定が得られる。この推定は、ステップ265において、低解像度領域に適用される変分法の初期化に使用される。そして、ステップ270において、速度についてのロバストな低解像度推定が、高解像度領域にマップ化され、高解像度の探索領域を定義するために使用されるが、これは、ステップ270において、対象の各参照位置における最大解像度相関関数を計算するためである。また、変分法から得た低解像度推定は、最大解像度相関関数による速度場の反復的な推定にとって、良好な初期推測として役立つ。そして、ステップ275で得られた最大解像度の速度推定は、ステップ280での表示や後続の歪み計算に使用される。弾性率の決定及び改良された歪みマッピングOCTに基づく組織弾性画像化については、組織弾性率及び歪み分布を併せて推定するために、図4に示すような一元化された計算フレームワークを利用するが、これは、以下の基本手順から構成される。1.参照OCTデータセットによる数値的再構築ステップ210,215。a.OCT区分化及び組織分類。b.弾性率の値及び境界条件の初期化。c.測定の圧力負荷の適用。この負荷は、研究対象の系に固有例えば、心臓周期中の血管内圧の正常な変動であってもよく、また外部的に制御されてもよい。d.血管変形の予測を得るための、例えば有限要素モデリングFEMなどの数値シミュレーション。2.参照OCT画像データと探索OCT画像データとの間での、モデルの剛体並進及び回転の推定ステップ215。3.剛体変換と、参照OCTデータセットを歪ませるための数値予測された変形場との組み合わせステップ225。4.OCT固有のデータ忠実度項の計算ステップ230。5.弾性率の更新と、OCTデータ忠実度を最大にするための剛体変換、及び更新された数値モデルのシミュレーションステップ210。6.弾性率ベクトル推定値の収束後における、最終的な弾性率及び歪みの表示ステップ235,240。例えば、取得された参照OCTデータは最初、血管壁の形状を画定するために区分化されるステップ210。勾配に基づく能動輪郭を用いて内腔境界が抽出され、これは、血管のOCT画像化において、大きな強度勾配の値を示す。好ましい実施形態では、レベル設定に基づく能動輪郭が使用される。このレベル設定法では、動脈表面が、高次元の埋め込み空間のゼロレベルセットとしてモデル化される。閉じた平面曲線Cp:[0,1]→R2の場合には、埋め込み表面embeddingsurface、u:→R2→Rは、曲線に対する符号付きの距離関数によって表される。uの各点に割り当てられる値は、曲線C上での最も近い点までの符号付き距離である曲線の内側で負、外側で正である。画像内の物体境界に最も適合する曲線を特定するために、以下の曲線の発展方程式が離散化されて、勾配降下法によって反復的に解かれる。gは画像勾配値の関数であり、κはユークリッド曲率であり、は、単位法線ベクトルである。曲率項により、大きな画像勾配が存在する場合以外は、曲線が滑らかとなり、また、画像から導出される曲線伝播項は、物体境界において曲線を大きな勾配に向けて引きつける。uはCの間接的表現であるため、式28を解くことは以下の式を解き、曲線Cを抽出するためにゼロレベルセットを選択することと同じである。式29のレベルセット発展方程式は位相幾何学的に柔軟であり、血管分岐の複雑な変化を自動的に受け入れる。エントロピーを満たす上流有限差分法は、数値の安定性及び速さをもって2を解くために、上述した狭帯域の更新法と共に使用される。OCTデータにおける血管壁の最も外側の境界は、強度閾値処理及び計算幾何学に基づいて定義することができる。画像化データには、OCTシステムの測定されたノイズフロアを超える信号強度を有する全ての点を探すために、最初に閾値が設定される。そして当業者には理解されるように、これらの点から幾何学的な凸包が形成されて、得られた表面が外側の血管境界を画定するために使用される。次に、内側及び外側の血管表面に含まれるOCTボクセルが、管理された最大事後MaximumAPosteriori:MAP分類、及び実験的に導出されたクラス条件付きclasso-conditionalの強度確率密度関数に基づいて、脂質に富む組織、繊維性組織又は、石灰化組織の領域に割り当てられる。本発明の例示的な一実施形態では、OCT及び組織学からのトレーニングデータが記録され、脂質に富むL組織、繊維性F組織及び石灰化C組織の領域が、OCTデータにおいて、これに対応する組織学データに基づいて位置される。各組織クラスに対して、OCT画像強度値が抽出されて、クラス条件付きの強度確率分布を近似する頻度ヒストグラムが生成される。この手順により、OCT強度Iに対して、任意のボクセルPI|F、PI|L及びPI|Cにおける3つのクラス条件付き確率分布が得られる。ベイズ式を使用すると、ボクセルが組織クラスF、L又はCに属する確率は、その強度値をgとして、下式で与えられる。先の組織クラスの確率PF、PL、PCは、各病変型の頻度に関する生体内の観察から導出することができ、これらの観察は各ボクセルについてもあてはまると仮定される。クラス条件付きの確率値から、任意の組織クラスを、あるボクセルに割り当てることに関連する分類誤差の確率が計算される。分類誤差の確率を最小化する組織割当は、ボクセル組織クラスに対して選択される。区分化と分類を行った後で、血管形状、及び、同様の弾性を有するプラーク内ボクセル間の境界面を規定する表面を得ることができる。そして技術的に理解されるように、これらの表面は、市販されるか又はカスタムコード化された有限要素分析プログラムにおいて、有限要素モデルを作成するために使用される。このプログラム又は当業界で公知の代替のメッシュ生成手順を用いて、血管及びプラーク内成分の有限要素メッシュが生成される。境界条件については、動脈入口における内腔輪郭の固定点と、径方向に自由に移動できる内腔重心の反対側の点を定義することによって適用される。測定された圧力負荷は、血管モデルの内腔表面に適用されて、対応する画像領域におけるOCT組織分類の結果と、血管の生体力学に関する文献から得た脂質組織、繊維性組織及び石灰化組織の平均弾性率値に基づいて、各メッシュ要素に対して初期弾性率推定値が割り当てられる。本発明による技術の1つの例示的な実施において、各メッシュは、入力OCTデータが二次元画像であるか、三次元ボリュームであるかに応じて、二次元の四辺形又は三次元のアイソパラメトリックな線形弾性有限要素のいずれかより構成される。変位uに基づく要素又は変位/圧力u/pに基づく要素は、ポアソン比0.5の手法でのロッキングの回避に使用できる。しかし、OCTによる組織弾性画像化のための一元的なフレームワークについては、このような特定の構造要素又は構成的材料モデルの使用に限定されない。例えば、非線形Mooney-Rivlin歪みエネルギー関数を、本実施形態において使用する線形材料モデルに置き換えことができる。モデル構築後に、有限要素シミュレーションを行って、探索OCTデータセットにおける変形メッシュのノード位置が予測される。剛体の並進及び回転は、参照画像と探索画像との間で血管内に生じ、上記したFEMモデルは、適用された境界条件のために、これを反映しない。よって、当業界で公知のように、多次元画像相関法のための標準的な方法を用いて、参照データセットと探索データセットとの間で、有限要素モデルにおける剛体の並進及び回転が計算される。そして、推定された剛体の並進がFEMで予測されるメッシュ変形に適用されて、参照データセットと探索データセットとの間で有限要素モデルの全体的な挙動が得られるステップ215を参照。FEMのノード変位は、剛体運動を受けるときにメッシュ生成し、変形は、参照データの空間座標を、探索データセットの予測位置に対してマッピングにより位置付ける歪み場を規定する。この歪み場は、多次元補間を用いてOCT参照画像化データを空間的に再サンプリングする225ために使用される。好ましい実施形態では、線形補間カーネルがその単純性及び計算速度を理由に使用される。しかし、原理上は、これに代えて、スプライン補間又は三次補間等、任意の多次元補間法が使用可能である。結果としての歪められた参照データは、探索データ及びOCTドップラ測定値と共に、以下の項の線形結合によって定義されるOCT特有の目的関数を計算するために使用されるステップ230を参照。1.歪められた参照画像と探索画像中の測定されたOCTデータの強度間での二乗差についての負の和。2.歪められた参照画像と探索画像中の測定されたOCTデータの強度間での正規化された相互情報量。3.歪められた参照画像と探索画像中の測定されたOCTデータの強度間での相関係数の要素に亘る総和。4.加えられた圧力/変位負荷による、組織又は境界の造影剤の光学特性についてモデル予測の変化と測定された変化との間の二乗誤差のノードに亘る総和。5.画像化ビーム方向における、光学的ドップラ速度測定値と、FEM予測の変位との間の内積のノードに亘る総和。この例示的な目的関数については、最初に弾性率及び剛体変換パラメータを更新し、次に、参照フレームと探索フレームとの間の血管変形の改善された予測を得るために、各反復においてFEMモデルを繰り返し実行することによって、反復的に最大化される。当業者には理解されるように、全体的な目的関数を最大化にするためには、多次元の制限付き共役勾配法を使用することができる。この例示的な方法による制限は、各要素のとりうる弾性率値の範囲を制限するために課され、これは、生体力学的挙動生体力学試験から先験的に測定される、つまり、OCT画像強度に基づいて要素に割り当てられる組織クラスでの典型的な挙動を基礎とする。制限された関数を最大化するためのこの手順は、代用可能な多くの可能な技術のうちの1つである。本発明の1つの例示的な実装では、最大弾性率値の絶対変化が0.0001を下回る場合に、推定手順が完了したとみなせる。収束時には、最終的な弾性率分布と、これに対応する有限要素法シミュレーションからの、対応する歪み及び応力の場について、入力データの次元に応じて、色分けによりマップ化された画像又はボリュームレンダリングとして表示することができる235,240。自動レンジ調節本発明はまた、自動レンジ調節技術を使用することができ、これには、2002年4月30日に出願された同時係属中の特許文献4に説明する処理法が含まれ、この出願の全ての開示が参照によりここに援用される。1つの例示的な実施形態では、自動レンジ調節メカニズムにより、本発明の技術を血管画像化に適用できるようになり、これにより、画像カテーテルを血管内腔の中心に位置させる必要がなくなる。この例示的な実施形態では、自動レンジ調節メカニズムのフィードバック信号を本発明の画像化メカニズムに組み込むことが望ましく、これは、例えば、血管構造の精度の高い決定を容易にし、かつ維持するためである。実施例以下の説明は、本発明による方法の例示的な実施形態について、実験的な検査に関する詳説である。具体的には、例示的な多重解像度変分法が、円形の内包物を含む組織ブロックの軸方向における圧縮中に、シミュレーションによるOCT画像化において実行された。図6は、組織ブロック及び円形の内包物のために使用した、有限要素モデルの形状300、及び、これに対応する有限要素メッシュ305を示す。連続する干渉画像は、式1~6で説明したように、指数関数的な減衰項と、コヒーレントOCT点広がり関数と試料内を移動する点状散乱体から生じる後方散乱分布との間の畳み込みと、の積を計算することによって生成された。組織ブロック中の離散的な点における後方散乱値については、ブロック中の散乱体に対する自由度10と、円形の内包物中の散乱体に対する自由度2でもって、独立した一様な乱数変数としてシミュレーションが行われた。これらの値は、内包物に対するブロック内の平均後方散乱値を大きくするために経験的に選ばれた。得られたコントラスト差は、OCT画像内の脂質と正常な動脈壁との間で観察されたコントラスト差を模倣する。圧縮中における組織の散乱体の運動については、組織形状の有限要素モデリングから、変位場を使用してシミュレーションが行われた。円形の内包物315が組み込まれた二次元の矩形状をした幾何形状310が定義された。全シミュレーションにおいて内包物の直径は500μmマイクロメートルであった。不動点境界条件がブロック320の下端中央に課され、緩い境界条件が下端325における他の全ての点に与えられた。軸方向下流の変位負荷が上面330に加えられ、5回の時間ステップでブロックが0.15mmだけ圧縮されて、全体的な歪みが4.3となった。各時間ステップでは、負荷により、ブロック内に約0.9の歪みが生起された。ブロックと内包物はともに、ほぼ非圧縮の線形弾性材料ν=0.495としてモデル化された。全てのシミュレーションについて、ブロックには単位弾性率が割り当てられる一方で、内包物の弾性率は、線維組織に埋め込まれた脂質に富む病巣又は石灰化の病巣を表現するように変更された。線維プラークに対する脂質の係数比は約0.0001で、線維プラークに対するカルシウムの係数比は約5である。有限要素モデリングは、ADINA8.0マサチューセッツ州ウォータータウンを使用して、9ノードの四辺形での二次元平面歪要素からなるメッシュ335,340を用いて実施された。メッシュ密度は、各要素のエッジの長さによって定義され、内包物340の内外では0.025mm、周辺ブロック335では0.1mmであった。各シミュレーションモデルは、約3200の要素及び13000のノードから構成された。各時間ステップで計算される変位場は、連続するOCTフレーム間での、組織の散乱体の速度ux,y及びvx,yを表現するために使用された。後方散乱場σbx,y,t1は、時刻t1での移動する組織散乱体に関連しており、下式を使用して計算された。σbx,y,t1=σbx-ux,y,y-vx,y,t28実際には、第1フレームにおける組織散乱体の場が、最初にアップサンプリングされて、次に式28の線形補間を用いて不均一に再サンプリングされることで、後続フレーム中の組織散乱体の場が得られる。式2~4におけるコヒーレントPSF近似を用いることにより、本OCTシステムの軸方向のフリンジパターンの測定に基づいて、PSFのシミュレーションが可能である。図7のグラフ350は、シミュレーションを行ったOCT点広がり関数PSF365に対する、フリンジ分解OCTPSF360の測定を示す。シミュレーションと測定は、ある一致した特徴を有する。このコヒーレントPSFについて、式28からの組織散乱体の場との畳み込みを行った後で、軸方向の指数関数的な減衰との乗算が行われ、ダウンサンプリングが行われることで、ピクセルサイズが1μm×25μmとされる一連のシミュレーション画像Ix,y,tが得られる。このようなシミュレーションに存在するスペックル又は乗算ノイズの度合いは、以下のノイズモデルを使用することより、更に制御された。Inx,y,t=Ix,y,tnIx,y,t29ここで、nは、平均値がゼロで分散がσn2とされた一様分布の乱数変数である。図7において、355は、図6の300、及びOCTPSF365におけるFEM形状から生成されたシミュレーションOCT画像の復調結果を示す。内包物370が、組織ブロック375内に見え、正常な動脈壁中の脂質プールに近い外観を有する。上記の方法は、迅速な処理のためにCプログラミング言語を使用して実装することができ、可視化及び分析はMATLAB米国マサチューセッツ州ナティック所在MathWorksで実行することができる。数秒程度という高速で全画像の計算を実現するために、式8の相関係数関数に対して高速な正規化された相互相関近似を使用したが、これには、この式の分子、及び分母について事前に計算される移動和を計算するために二次元FFTを使用する。相関の最大化に基づく従来の速度推定技術の場合、各相関関数は、ピーク周辺で双三次補間を用いて50倍のアップサンプリングが行われたが、これは軸方向で0.02μm、横方向で0.5μmのサブピクセルシフトを検出するためである。従来の方法では、相関関数を計算するために、81×7ピクセル81×175μmの参照ブロックサイズと361×21ピクセル361×525μmの探索領域が使用された。これらのパラメータは、速度推定の誤差を最小化する要求に対して、速度の空間変化の検出感度に対する要求を均衡させるために経験的に決定された。従来の運動追跡からの速度及び歪み推定についてのメジアンフィルタリングは、本発明による例示的な変分法から得た結果と比較する前に、5×5のカーネルで実行された。変分法の例示的な実施形態における低解像度ステップでは、各相関関数を計算するために、15×11ピクセル75×1375μmの参照ブロックサイズと、61×41ピクセル305×5125μmの探索領域が使用された。最大解像度では、25×7ピクセル25×175μmの参照ブロックサイズと、101×21ピクセル101×525μmの探索領域が使用された。変分法の重み付けパラメータは、シミュレーションされたOCTシーケンスのサブセットに基づいて、経験的にa=1、b=20、c=0.1に決定された。そして、これらの値は後述するOCT実験の速度推定に使用された。速度推定の後、変形行列Fが、以下の表式に対する第一次の有限差分近似を使用して速度場の全ての点で計算された。歪みを微小とする近似を仮定した場合に、変形行列は、以下の関係によって歪み行列E及び単位行列Iに関連付けられる。軸方向圧縮のシミュレーションOCT画像化では、歪み行列の軸方向における歪み成分yyの結果を示し、また横方向の伸張における画像化実験では、横方向の歪み成分xxの結果を示す。速度及び歪み場の推定の誤差は、正規化された二乗平均平方根の誤差測定に基づいて評価された。ここでNは推定値の合計数、vkは軸方向又は横方向の推定速度、vk,realは軸方向又は横方向の実際の速度、kは速度場中のk番目の対象点での軸方向又は横方向の推定歪み、k,realは、速度場中のk番目の対象点での軸方向又は横方向の実際の歪みである。シミュレーション実験500μmの内包物の検出能力が、スペックル非相関性の関数として検査され、ロバストな推定から得た値に対して、従来の運動推定から得た軸方向の速度及び歪みの推定におけるRMS二乗平均平方根誤差が比較された。vk,real及びk,realの値は、有限要素モデリングから直接得た。速度は、下方及び右方への変位にそれぞれ対応した、正の軸方向速度及び正の横方向速度についてピクセル単位で報告された。図8、図9はそれぞれ、弾力的な内包物と固い内包物の運動シーケンスにそれぞれ対応するFEMによる軸方向速度場400,450を示す。また、これらは、従来の運動追跡405,455と、変分法フレームワーク410,460内でのロバストな推定から得た速度測定結果の軸方向成分を示す。これらの例では、内包物の位置は、符号415,420,425,465,470,及び475に対応している。この2つの方法から推定された速度場は、FEMから得た真の速度と定性的に同じである。変分法による推定値は、従来の追跡よりも著しく滑らかとなる。更に、軸方向速度のRMS誤差は、ロバストな推定よりも従来の速度測定による方が大きい。固い内包物の場合、従来の速度測定による軸方向速度場のRMS誤差が1.60であったのに対し、ロバストな推定のRMS誤差は1.04であった。同様の結果は、弾力的な内包物についても得られた従来法の1.83に対し、変分法では1.40。図10、図11は、弾力的な内包物と堅い内包物のそれぞれのシミュレーションに対応する軸方向の歪み場を示す。フレーム500,550は、有限要素モデリングによる真の軸歪みを示し、フレーム505,555は、従来の運動追跡による軸方向歪みの推定を示し、フレーム510,560は、多重解像度変分法による軸方向歪みの推定を示す。これらの画像における内包物の位置は、符号515,520,525,565,570,及び575に対応している。歪み場測定に関する速度ノイズの影響は、これらの例から明らかである。従来のロバストな速度測定による推定は、FEMで導出される変位と定性的に同様であるが、歪み推定は、真の歪に対して同様ではない。軸方向歪みは軸方向速度場の導関数であり、このため、速度推定に存在するノイズは、導関数演算子のハイパス高域通過特性のため、歪み画像内で強調される。メジアンフィルタリングを用いた場合であっても、従来の運動追跡から導出される歪みについて、内包物内の歪みの差の程度と大きさを充分に評価することは困難である。これとは対照的に、ロバストな歪み推定を用いれば、内包物の位置及びサイズを決定することが可能である。更に、変分法によって、弾力的な内包物と固い内包物とを視覚的に区別することが可能である。定量的には、従来法と変分法との間のRMS歪み誤差の差は顕著である。例えば、固い内包物の場合、従来の手法ではRMS歪み誤差が109.1であるが、変分法ではRMS歪み誤差が27.5である。従来の追跡から得られる歪みは、解釈が困難であるのに対し、ロバストな推定では、病変型の間での比較について、歪み測定をより簡単にかつ高精度をもって解釈可能となる。画像化実験中心波長1310nm及び帯域幅70nmの広帯域光源の光が、干渉計内で参照場とサンプル場に分かれる。サンプル場は、参照アームの光路長に対応する深さで組織を調べるために、走査光学装置によって合焦される。サンプルアームからの戻り後方散乱光が参照場と混合されて干渉信号が生成され、該信号は、25μm横方向×1μm軸方向のピクセルを生成するためにデジタル化される。干渉信号の振幅は、参照アームに規定される走査深さにおける組織構造及び光学特性に関する情報を含んでいる。XY断面における組織構造については、参照アームの光路長を変更することにより軸Y方向で探索され、検体を横切るサンプルビームを掃引することによって横X方向で探索される。この例示的な方法では、横方向500ピクセル×軸方向2500ピクセルからなる画像フレームが、250ミリ秒で取得される。検視解剖で得たヒトの大動脈の正常断片は、リン酸緩衝生理食塩水中で37°Cに暖められた。画像化は、断片を取り出して24時間以内に実施された。円筒状の大動脈の断片が、縦方向に切断されて開かれ、管腔の表面を検体光束に露出させる矩形の組織検体が得られた。横方向の機械的負荷が、無処置の大動脈断片の周縁部における伸張に近くなるように、縦方向に切断された端が検体ホルダに固定された。検体は検体ホルダ内に水平に取り付けられ、これにより、その一端がしっかりと固定され、他端がマイクロマニピュレータ、つまり、水平軸に沿って一次元の平行移動ができる極微操作装置に固定された。試料内の画像化位置は、検体光束と同期する照準ビームレーザダイオード、635nmを視覚化することによって監視された。走査位置は検体中心に設定され、サンプルアーム内の走査方向は、散乱体の変位が画像化面内に制限されるように位置合わせされた。この画像化実験において、速度及び歪みの推定が、静止又は横方向の伸張のいずれかの状態下にある正常な大動脈の検体で調べられた。静止状態の場合、グランドトゥルースの速度場及び歪み場はともにゼロとすべきである。このため、速度及び歪みの推定の、ゼロからの偏差が1つの性能評価基準として決定された。横方向の伸張の場合、予測された歪み分布は均一であり、推定された歪みの標準偏差を、試料内の歪みの均一性についての評価基準として報告する。静止状態では、従来法で推定された平均横方向速度が、フレーム当たり0.17μmであり、これは、名目上用いられる並進量、フレーム当たり0ピクセルと比較される。場についての横方向速度推定の標準偏差は、フレーム当たり3.94μmであった。軸方向では、平均偏差及び標準偏差がそれぞれ0.66μm及び4.28μmであった。変分法による速度推定は、横方向でフレーム当たり0.0017μm軸方向ではフレーム当たり0.54μmの平均速度を示し、これは、従来法の結果よりも、実際の並進量に近かった。また、横方向でフレーム当たり0.16μm軸方向ではフレーム当たり0.58μmの速度標準偏差は従来の方法の結果よりも小さかった。場の歪み推定の平均偏差及び標準偏差は、従来法の場合に横方向ではそれぞれ-0.043と5.5軸方向では-0.24と11.26であったが、これは場の名目上用いられる歪み0と比較される。これに対して、ロバストな推定では、一般に、より良好な歪み推定性能を示す。歪み推定の平均偏差及び標準偏差は、横方向ではそれぞれ-0.0038と0.12軸方向では-0.024と0.12であった。わずかな横方向の伸張による横方向の移動を受けた大動脈断片についての例示的なOCT画像を、図12に600で示す。符号605で示した領域がマスクされており、これは大動脈の組織610の表面より上の空気に対応している。図13は、速度場測定の横方向成分を示す。フレーム650は、従来の速度推定手法によって得られ、またフレーム655は、本発明による変分法の例示的な実施形態を使用している。組織より上の空気に対応する領域660,665がマスクされ、領域670,675は大動脈組織に対応している。図14は、推定された歪み場の横方向成分を示す。フレーム700は従来法によって得られ、またフレーム705は、本発明による変分法の例示的な実施形態から得られた。組織より上の空気に対応する領域710,715がマスクされ、領域720,725は大動脈組織に対応する。ロバストな推定から得た速度は、従来の推定から得た推定値よりも、横方向及び軸方向において滑らかに変化することが認められた。従来法では、歪みの平均偏差及び歪みの標準偏差は、横方向で0.22、47.26であり、軸方向で-3.46と85.08であった。これに対して、ロバストな推定による歪み場は、平均偏差と標準偏差が、横方向では-0.026と1.32であり、軸方向では-0.02と0.14であることを示した。これらの結果は、ロバストな歪み計測から得た推定値が、従来の歪み計測から得た推定値よりも、予測される均一な歪み場に対して、より合致していることを示している。本発明の更に別の例示的な実施形態によれば、例示的なフレームワークが、効率的なパラメータ削減策に使用可能な生体力学特性のOCTベースの推定に提供される。このような方策は、基礎データに存在し得る機械的特性の急峻な勾配を維持しつつ、迅速な手術中での推定を支援することができる。通常、パラメータを減らすことと空間分解能とはトレードオフ二律背反の関係にある。組織領域同士を正しく判別するためにOCTの軟組織のコントラストを使用するとともに、同種の組織が類似した力学的挙動を示すという知見によって、このトレードオフを低減させるか又は最小にし、計算性能を改善することが好ましい。これを実現するための本発明による方法の1つの例示的な実施形態が、図15に示されており、以下の通りである。a.取得したOCT参照フレームを主要な組織領域に区分化するステップ805。区分化は、同様の平均強度、強度分散、Aラインの深さの関数としての勾配又は他の候補となる強度計量を用いたクラスタリング領域に基づいて自動化することができる。異なる領域の数は、推定すべき未知パラメータの数例えば、約2~6個の未知パラメータであり得る。b.ステップ810において、特定された組織の形状及び境界が、粗い有限要素モデルメッシュを生成するために使用されるステップ810。組織境界は、メッシュ生成ステップへの入力として使用されるため、要素が1つ以上の組織型に跨ってはならない。各要素は1つの組織型に属し、該要素はその組織についての対応する未知の弾性率に割り当てられる。c.メッシュの変位境界条件を推定するステップ815。d.ステップ815で得た境界条件を使用することによって、所望の目的関数を最小化する値を特定するために、上記非特許文献6乃至8に記載されるように、低減された弾性率パラメータの組に亘って反復される最適化法で前進型有限要素法シミュレーションを拡張するステップ820。e.各要素において目的関数の値を決定するステップ825。目的関数の値が所定の許容閾値を超える要素は、有限要素モデルがデータに十分に適合していない位置に対応する。このような要素は、更なる精緻化のためのリストを形成することができる。f.精緻化リストの要素を更に分割して、元の要素が占めていた空間に及ぶ、より小さな新要素を形成するステップ830。このような精緻化された要素はそれぞれ、元の要素組織クラスについて推定された弾性率に初期化される、独立した弾性率の値をもつことができる。g.元の組織クラスの割当に特有とされる限定範囲に亘って値が変化するように、精緻化された要素に対して制限条件付き最適化法を用い、ステップ820~825を繰り返すステップ835。h.収束時に、最終的な歪み及び弾性率の分布を表示するステップ840。メッシュの精緻化のための閾値を適切に調整することによって、未知パラメータの数を、従来法の最大次元数、つまり全要素が独立の未知数として処理される次元数から大幅に減らすことができる。上記説明は、本発明の原理を示したに過ぎない。記載した実施形態を様々に変更及び変形することは、ここに記載した教示を鑑みれば、当業者にとって明らかであろう。例えば、ここに記載した本発明は、2003年10月27日出願の特許文献3及び2003年1月24日出願の特許文献5に記載した例示的な方法、システム及び装置において有用である。これら出願の開示は、その全体をここに参照することで援用される。よって、当業者には分かるように、ここに明確に図示又は記載されない、数多くのシステム、構成及び方法を工夫し、本発明の原理を具体化できるが、これらは本発明の趣旨及び範囲に含まれる。また、上記に示した全ての刊行物、特許及び特許出願は、その全体が参照によりここに援用される。"