(57)【特許請求の範囲】
【請求項１】
光学コヒーレントトモグラフィ（ＯＣＴ）カテーテルを自動較正するためのシステムであって、
カテーテル内の画像を撮影するように構成される撮像部と、
前記カテーテルに特有の前記カテーテルの複数の構成要素からの複数の内部カテーテル反射線の位置を識別するように構成されるテンプレートと
を含み、
前記画像と前記テンプレートとを照合することにより、前記画像における前記内部カテーテル反射線の位置を識別し、前記画像と前記テンプレートとのテンプレート一致が識別される場合、Ｚ−オフセットを調節することによって、前記内部カテーテル反射線をシフトさせ、
前記テンプレート一致が識別されない場合、前記Ｚ−オフセットを調節し、前記テンプレートを用いて評価されるべき新しい位置において、前記カテーテルの第２の画像を撮影する、システム。
【請求項２】
前記テンプレートは、特定の前記カテーテルとともに記憶される、または動作可能にアクセス可能である、請求項１に記載のシステム。
【請求項３】
前記テンプレートは、前記カテーテルの前記画像から直接導出される、請求項１に記載のシステム。
【請求項４】
前記テンプレートを前記画像のあるバージョンと照合し、ピークを識別するように前記システムが構成される、請求項１に記載のシステム。
【請求項５】
前記画像のバージョンは、バイナリである、請求項４に記載のシステム。
【請求項６】
最も強い前記テンプレート一致に対応する前記ピークを選択するように前記システムが構成される、請求項４に記載のシステム。
【請求項７】
前記ピークは、一定の値を上回る、請求項６に記載のシステム。
【請求項８】
前記Ｚ−オフセットを調節することは、参照アーム経路長を調節するための可変遅延線の機械的シフトを含む、請求項６に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願の引用）
本願は、米国特許出願第１３／４０２，７９８号（２０１２年２月２２日出願）の利益およびそれを基礎とする優先権を主張する。該出願の全体は、参照により本明細書に引用される。
【０００２】
（発明の分野）
本発明は、概して、較正システムおよび使用方法に関し、より具体的には、光学撮像システムのための較正システムに関する。
【背景技術】
【０００３】
正確な光学コヒーレントトモグラフィ（ＯＣＴ）測定または次元分析は、表示されるトモグラフィ画像が、物理的空間を正確に表すこと（すなわち、画像画素から物理的ｍｍへの変換）を要求する。本要件は、さまざまな物体屈折率（カテーテル光学、シース、管腔、組織）、および干渉計の試料経路と参照経路とにおける不一致による関連画像特徴の不定の場所（ｚにおける）等の要因によって複雑となる。
【０００４】
最近の方法は、ユーザが、カテーテルシースの外径が画面上の固定目盛りと整列するまで、Ｚ−オフセット位置（参照アーム経路長）を調節することによって、手動で画像を較正することを要求する。本方法は、時間がかかり、オペレータエラーを招く可能性がある。加えて、カテーテルが、元の較正された位置からシフトされると、較正は、時変機械歪み（例えば、ＰＩＭケーブルの引き戻し動作または操作）または温度変化（室温対体温）により、破棄され得る。
【発明の概要】
【課題を解決するための手段】
【０００５】
本発明は、前述の識別される問題に対処し、自動較正に関する。カテーテルは、本発明に従って、カテーテル内の画像を撮影し、テンプレートを利用して、カテーテル反射線の位置を識別することによって較正されることができる。カテーテルを較正するための別の方法は、カテーテル引き戻しを撮像し、画像データを整列させ、画像取得中の反射位置に対するいかなるシフトも補正されることを確実にし、各画像フレームを通して、反射位置を追跡する。さらにカテーテルを較正する別の方法は、カテーテルの画像を撮影し、カテーテルの反射線を追跡し、画像内の反射線を検出せず、グラフ検索ステップまたはテンプレート照合ステップを用いて、喪失トラック画像を再取得する。
【０００６】
本発明の前述および他の特徴、利点、ならびに目的は、以下の本開示を参照することによって、より明白となるであろう。付随の図と併せて読まれる、例示的実施形態の以下の説明は、限定ではなく、単に、例証である。
【図面の簡単な説明】
【０００７】
以下の簡単な説明は、図のより完全な理解のために提供されるが、本発明による実施形態は、必ずしも、示される精密な配列および構成に限定されないことを理解されたい。
【図１】図１は、一実施形態による、較正のために使用される、カテーテル反射を実証する試料フレームを示す、グラフである。
【図２Ａ】図２Ａは、一実施形態による、ビン３４０においてピークを識別する、テンプレート照合アルゴリズムのグラフであり、図２Ｂは、一実施形態による、ビン３４０においてテンプレート一致を示す、領域Ａが拡大されたグラフである。
【図２Ｂ】図２Ａは、一実施形態による、ビン３４０においてピークを識別する、テンプレート照合アルゴリズムのグラフであって、図２Ｂは、一実施形態による、ビン３４０においてテンプレート一致を示す、領域Ａが拡大されたグラフである。
【図３】図３は、一実施形態による、１つのフレームを通した反射運動（蛇行性引き戻しによる運動）の拡大画像である。
【図４】図４は、一実施形態による、グラフ検索アルゴリズムである。
【図５】図５は、一実施形態による、カテーテル引き戻しシーケンスの画像フレームおよび画像フレームを通して維持されるトラックである。
【図６】図６は、一実施形態による、カテーテル引き戻しシーケンスの画像フレーム４０であり、カテーテル反射線は、存在せず、トラックが、喪失されている。
【図７】図７は、一実施形態による、カテーテル引き戻しシーケンスの画像フレーム５０であり、カテーテル反射が再現し、トラックが、自動較正のために再取得される。
【図８Ａ】図８Ａは、一実施形態による、テンプレート照合ステップ１００の第１のモードを表す、流れ図であり、図８Ｂは、代替実施形態による、テンプレート照合ステップ１００の第１のモードを表す、流れ図である。
【図８Ｂ】図８Ａは、一実施形態による、テンプレート照合ステップ１００の第１のモードを表す、流れ図であって、図８Ｂは、代替実施形態による、テンプレート照合ステップ１００の第１のモードを表す、流れ図である。
【図９】図９は、一実施形態による、再生追跡ステップ２００の第２のモードを表す、流れ図である。
【図１０】図１０は、一実施形態による、喪失トラック再取得ステップ３００の第３のモードを表す、流れ図である。
【図１１】図１１は、一実施形態による、再生方法４００における自動較正を表す、流れ図である。
【図１２】図１２は、一実施形態による、自動較正初期ロック方法５００を表す、流れ図である。
【図１３Ａ】図１３Ａ−１３Ｄは、平均振幅または勾配を使用して、強い反射を検索し、検出されるまで、Ｚ−オフセットを調節する、初期ロックモードにおける第１のステップのグラフである。
【図１３Ｂ】図１３Ａ−１３Ｄは、平均振幅または勾配を使用して、強い反射を検索し、検出されるまで、Ｚ−オフセットを調節する、初期ロックモードにおける第１のステップのグラフである。
【図１３Ｃ】図１３Ａ−１３Ｄは、平均振幅または勾配を使用して、強い反射を検索し、検出されるまで、Ｚ−オフセットを調節する、初期ロックモードにおける第１のステップのグラフである。
【図１３Ｄ】図１３Ａ−１３Ｄは、平均振幅または勾配を使用して、強い反射を検索し、検出されるまで、Ｚ−オフセットを調節する、初期ロックモードにおける第１のステップのグラフである。
【図１４】図１４Ａ−１４Ｂは、長方形画像内のＡ走査を横断して反射を整列させる、初期ロックモードにおける第２のステップのグラフである。
【図１５】図１５Ａ−１５Ｂは、勾配および振幅を使用して、カテーテル反射を識別する第３のステップのグラフである。
【図１６】図１６Ａ−１６Ｂは、反射のテンプレートを生成し、後の使用のために記憶する、初期ロックモードにおける第４のステップのグラフである。
【図１７】図１７は、初期Ｚ−オフセット較正の流れ図である。
【図１８】図１８Ａ−１８Ｄは、テンプレートをシフトさせ、場所を検索し、ミラー信号を伴う「完全テンプレート」を作成する、ライブ追跡モードにおける第１のステップのグラフである。
【図１９】図１９は、テンプレートおよびテンプレートサイズに限定される全Ａ走査の相関係数を算出するためのライブ追跡モードにおける第２のステップのためのグラフおよび式である。
【図２０】図２０は、各Ａ走査に対する最大相関を見つけ、上位「ｎ」個のＡ走査の中央値を求めるためのライブ追跡モードにおける第３のステップのグラフである。
【図２１】図２１Ａ−２１Ｂは、デジタルシフトを画像および走査変換に適用することによる、第４のステップにおける較正前後のＯＣＴ画像である。
【図２２】図２２は、一実施形態による、ライブモード追跡の流れ図である。
【図２３】図２３Ａ−２３Ｂは、大検索領域を伴う全Ａ走査を横断して、最大相関を使用して反射追跡の開始位置を決定するための再生追跡モードにおける第１のステップのグラフである。
【図２４】図２４は、同一の相関技法を使用するが、小検索領域を伴う、Ａ走査毎に、後方反射および前方反射を追跡する、開始位置からの再生追跡モードの第２のステップのグラフである。
【図２５】図２５Ａ−２５Ｂは、継続して、データセット内の全フレームを通して、Ａ走査毎に追跡し、後の整列およびビデオモニタまたは他の物理的表示デバイス上で閲覧者への表示のために、反射位置を記憶する、再生追跡モードにおける第３のステップのグラフである。
【図２６】図２６は、再生モード追跡アルゴリズムの流れ図である。
【図２７】図２７は、自動較正の一実施形態の流れ図である。
【発明を実施するための形態】
【０００８】
一般に、本発明による自動較正システムおよび方法は、内部カテーテル反射を検出し、内部カテーテル反射をシフトさせることにより、画像を較正する、再現可能な方法を提供する。一実施形態では、内部カテーテル反射は、光ファイバケーブル端、ミラー、レンズ、シース、流体、生物学的脈管、または反射を生じさせる任意の他の物体等による反射を含む。内部カテーテル反射は、機械的および／またはデジタル的にシフトされることができる。概して、自動較正は、第１のモード、第２のモード、および第３のモードを備える。較正システムおよび方法は、初期較正後、連続的なフレーム毎ベースで較正を更新および維持する。
【０００９】
光学コヒーレントトモグラフィ（ＯＣＴ）システムは、時としてスペクトルドメインＯＣＴ（「ＳＤ−ＯＣＴ」）としても知られるフーリエドメインＯＣＴ（「ＦＤ−ＯＣＴ」）、または時間ドメインＯＣＴ走査（「ＴＤ−ＯＣＴ」）を含み得、そこでは、干渉計の参照アームにおける光の光学経路長が、撮像深度範囲に対応する距離にわたって高速で走査される。ＯＣＴシステムは、高偏光感度または高位相感度であり、適宜、調節され得る。代替として、撮像システムは、限定されないが、分光法、（蛍光、吸収、散乱、およびラマン分光を含む）を含む、任意の他の光学撮像ベースのシステムであり得る。
【００１０】
（ＯＣＴ深度較正および自動化された範囲調節）
円形および円筒形ＯＣＴ走査デバイス、すなわち、回転カテーテル走査デバイスが、本質的極座標系（例えば、長さおよび幅ではなく、半径および角度）で物理的空間をサンプリングする。円形および円筒形ＯＣＴ走査デバイスは、円筒形状断面、例えば、気道および血管管腔を伴う生理学的構造を撮像するために適用される。画像のデジタル表現（すなわち、数値を表す画素のアレイ）は、本質的に長方形である。意図的または人工的に生成されたＯＣＴ画像特徴を検出および使用するための方法は、自動的に、極座標（「レーダ状」）ＯＣＴ画像内の深度範囲を調節することを含む。
【００１１】
極座標ＯＣＴ画像は、ビデオモニタまたは他の物理的表示デバイス上において閲覧者に表示する前に、その長方形表現から変換される。加えて、定量値（例えば、管腔直径、管腔面積、円周等）が、極座標画像上で測定されるべき場合、長方形から極座標への変換は、全次元（半径方向および角度）において画素間の相対的距離を保持する。概して、ＯＣＴ深度走査（直交座標のｙ軸）は、直接、半径にマッピングし、ＯＣＴ円周方向走査（直交座標のｘ軸）は、２＊Ｐｉラジアン（または、３６０°）極座標角度のうちのある増分にマッピングする。
【００１２】
例えば、ｙ＝０（長方形画像の最上行）は、半径＝０（極座標画像の中心）にマッピングし、ｙ＝ｙｍａｘ（長方形画像の最下行）は、半径＝ｙｍａｘ（極座標画像の周縁）にマッピングする。同様に、ｘ＝０（長方形画像の左列）は、角度＝０°にマッピングし、ｘ＝ｘｍａｘ／２は、約１８０°にマッピングし、ｘ＝ｘｍａｘは、約３５９°の角度にマッピングする。
【００１３】
極座標画像内の正確な定量的次元測定の場合、半径＝０への画素のマッピングは、撮像プローブの回転軸の中心における実際の物理的空間を表し、そうでなければ、極座標画像は、半径方向に人工的に歪められる（拡大または縮小される）であろう。しかしながら、任意のＯＣＴ画像では、ｙ＝０における画素は、必ずしも、この要件を充足せず、極座標表現へのマッピング前に、これが充足されるまで、ｙ−次元においてシフトされなければならない。試料の経路長対干渉計の参照アームの変位差（制御または非制御のいずれか）は、ｙ−次元において画素をシフトさせるであろう。
【００１４】
制御不能変位は、円筒形または螺旋走査光ファイバＯＣＴカテーテルを使用するときに生じ得る。例えば、カテーテルが、縦方向に押し引きされると、光ファイバケーブルは、圧縮または伸長され得、したがって、経路長変位を被る。
【００１５】
本方法は、概して、自動的に、静止しているが制御不能変位によるものではない画像特徴を検索することによって非制御変位効果を認識することと、極座標表現が、正確な次元測定のために使用され得るように、連続ＯＣＴ画像データを較正することとを含む。一実施形態では、本方法はさらに、ビデオモニタまたは他の表示デバイス上での表示に先立って、画像内の画像特徴を除去することを含む。
【００１６】
本方法によって使用される画像特徴は、カテーテル自体内（撮像される被写体または周囲環境内ではなく）で生成され、カテーテルの３６０°回転全体を通して、幾分、深度が安定し、強度が一貫して見える。これらの画像特徴は、光学構成要素間の界面における後方反射（「ゴースト線」または「エコーアーチファクト」として知られ、これらは、回転部分の光学軸に沿って生じ、したがって、カテーテルの１回転の過程にわたって、経路長変位差が生じないとき、極座標画像内で均一な円に見える）、あるいは、静止した（非回転）カテーテルシース（断面外形が円形であって、また、回転部分と機械的に同心である場合）の境界から、またはその中からの反射が含まれるが、それらに限定されない。
【００１７】
本明細書に開示される実施形態は、適切な較正を達成するための必要とされるシフトを識別するために複数の後方反射を利用する、自動較正のための３つの方法を含む。３つの方法の各々の間には、重複が存在し得るが、３つの方法の各々は、説明目的のためだけに、別個のセクションにおいて記述されるが、３つの方法の各々は、代替構成、方法、パラメータ等において組み合わせられ得る。第１の方法は、平均的かつ一般的自動較正実装である、Ｚ−オフセットの自動較正を含む。第２の方法は、テンプレート照合およびグラフ検索方法を含む、Ｚ−オフセットの自動較正を含む。第３の方法は、完全テンプレート相関を含む、Ｚ−オフセットの自動較正である。
【００１８】
（方法１：Ｚ−オフセットの自動較正および平均化）  一実施形態では、自動認識および較正方法におけるステップは、（１）ｘ−（すなわち、角度）次元に沿ってＯＣＴ画像フレームを平均化することにより、ｙ−次元（すなわち、半径）において回転安定する特徴を、被写体または周囲環境によって生成される他の画像特徴に対して選択的に強めることを含む。平均化ステップの有効性は、周囲画素に対して高強度を有する画像特徴を選択することによって、かつ被写体／環境特徴（雑音）が、強い円周方向対称性を有していない場合、改善される。一実施形態では、本方法はさらに、（２）ピーク検索、相関、閾値化、または他のパターン認識アルゴリズムを使用して、画像特徴を見つけることを含む。画像特徴を見つけることの有効性は、非制御経路長変位が生じ得る範囲が、先験的に既知である場合、したがって、必要とされる検索空間を限定する場合、改善される。一実施形態では、本方法はさらに、（３）ステップ２で見つかった画像特徴のｙ−値を、回転軸に対するその画像特徴の実際の物理的場所を表す、事前に較正されたｙ−値と比較すること、または、スペクトルドメインＯＣＴを使用する場合、その特徴の既知の「共役画像」または「エイリアス画像」の場所と比較することを含む。一実施形態では、本方法はさらに、（４）検索された画像特徴と事前に較正された画像特徴との間の差異だけ、ｙ−次元におけるＯＣＴ画像画素をシフトさせることによって較正することを含む。複数の特徴が、アルゴリズムの有効性を改善するために使用されることができる。ｙ−次元において長方形画像をシフトさせた後、極座標画像座標へのマッピングが、行なわれ得る。較正された極座標画像の中心に対して測定された半径は、物理的空間内の回転軸に対して測定された実際の半径を表す。カテーテルによるいくつかの画像特徴は、ビデオモニタまたは他の物理的表示デバイス上での被写体／環境特徴の効果的かつ混乱のない表示のために望ましくない。例えば、カテーテル画像特徴は、被写体／環境特徴に重複し得る。
【００１９】
一実施形態では、画像特徴を除去する（または、目立たなくする）ステップは、半径方向ｙ−方向において、かつ全列／角度にいて画像特徴範囲を切り取ることと、全列／角度に対して切り取られた領域のすぐ内側および外側（上方および下方）の画素の平均値を計算することと、この平均化された行／円周を切り取られた場所内に挿入することとを含む。切り取り動作はまた、被写体／環境特徴を除去し得、半径方向次元内の画像を歪ませる。この歪みは、測定ツールが、次いで、画素が切り取られた場所およびクロップされていない場所を考慮しなければならないので（または、クロップされていない画像上で測定を行わければならない）、そのような画像上の正確な定量値の測定をより複雑にする。
【００２０】
前述の較正実施形態では、較正方法は、反射を識別するためにフレームにわたって平均化し、次いで、特徴場所に基づいて画像をデジタル的に調節し、画像内の全Ａ走査に一定シフトを適用する。Ｚ−オフセットの自動較正のための代替方法は、画像内に現れる内部カテーテル特徴を使用して、必要とされるシフトを識別するものであり、この方法は、画像特徴を見つけるために、フレームにわたって画像強度を平均化しないが、テンプレートの形態における反射のパターンを使用して、初期ロックアルゴリズムにおいて、反射の位置を識別する。「テンプレート」とは、概して、カテーテル反射パターンを指す。この方法は、線毎のシフトを適用し、再生モードアルゴリズムにおける測定のために、全てのＡ走査が適切に整列されることを確実にする。Ｚ−オフセットの自動較正のためのバイナリテンプレート照合およびグラフ検索のこの方法は、以下により詳細に説明される。
【００２１】
（方法２：Ｚ−オフセットの自動較正およびテンプレート照合およびグラフ検索）  一実施形態では、Ｚ−オフセットの自動較正は、初期ロックステップを含む、カテーテル反射を較正する第１のモードを含む。初期ロックは、テンプレートを利用して、図１に示されるように、特定のカテーテルに特有のカテーテル反射線の位置を識別することを含む。各カテーテルのためのテンプレートは、特定のカテーテルとともに記憶される。一実施形態では、テンプレートは、無線周波数識別（ＲＦＩＤ）チップ上に記憶され、代替として、テンプレートは、コンピュータチップ等上に記憶され得る。代替として、ＲＦＩＤまたはコンピュータチップは、除去可能であり、次いで、ポータブルテンプレートまたは医療記録として使用され得る。代替として、カテーテルが、再使用のために承認される場合、患者インターフェースモジュールまたはＰＩＭは、特定のカテーテルに対するテンプレートに関する具体的情報をダウンロードし得る。このテンプレート情報は、記憶され、かつカテーテルモニタ上で追跡され、使用回数または使用時間を同様にカテーテル上に記憶された所定の量に限定し得る。一実施形態では、ＲＦＩＤチップは、情報を記憶し、インターフェースデバイスと通信するために、カテーテルの近位端のコネクタ上に搭載されたＭａｘｗｅｌｌ  ＭＥ１またはＭＥ２  ＲＦＩＤチップであり得る。代替実施形態では、カテーテルは、コネクタの第１のＲＦＩＤチップから１８０度に搭載された第２のＲＦＩＤチップ（図示せず）を有し得、したがって、カテーテルは、２つ以上の円周方向配向において、インターフェースデバイスに接続されることができる。ＲＦＩＤチップは、１２８バイト、代替として、１Ｋバイト、代替として、２Ｋバイト、代替として、４Ｋバイトのメモリを有し、例えば、カテーテル製造番号、名称、製造業者、またはモデル、較正係数、撮像要素感度、タイムゲインコントロール、ポストアンプゲイン、許容使用回数、地理的許容使用場所、ブートモード、パルス幅、またはカテーテルの有効期限を含む、カテーテル特有情報を記憶し得る。
【００２２】
テンプレートは、勾配画像のバイナリバージョンとコンボリューションされ、ピークは、図２Ａに示されるように、テンプレート照合ステップにおいて識別される。テンプレート照合ステップは、最も強いテンプレート一致に対応するピークを選択することを含む。一実施形態では、ピークは、最も強いテンプレート一致が識別されることを確実にするために、一定の値を上回る。テンプレート一致が識別されない場合、Ｚ−オフセット位置は、調節され、新しい位置における画像が、テンプレート照合ステップまたはアルゴリズムを用いて評価される。カテーテル反射線の位置が識別されると、カテーテル反射線は、Ｚ−オフセットを調節することによってシフトされ、反射をその所望の場所に移動させ、画像／カテーテルは、ここで、図２Ｂに示されるように、較正される。一実施形態では、Ｚ−オフセットの調節は、可変遅延線（ＶＤＬ）の機械的シフトを含む。テンプレート照合アルゴリズムが、強いテンプレート一致を識別できない場合、ユーザは、警告され、自動較正を再試行するか、または手動で較正するかの選択肢が与えられる。
【００２３】
テンプレート照合およびグラフ検索のためのＺ−オフセットの自動較正の方法２の一実施形態では、カテーテル反射を較正する第２のモードは、再生追跡を含む。再生追跡は、概して、取得の間の反射位置に対するいかなるシフトも、適切な分析および／または測定を可能にするように補正されることを確実にするために、カテーテル引き戻しを記録後、画像データを整列させることを含む。記録されたデータセットを通した反射の追跡は、引き戻しの間の動きにより、若干困難であり、反射の位置は、単一フレームにわたって、有意に変動し得る。図３は、蛇行性引き戻しの間のカテーテル反射の有意なシフトの実施例を実証する。グラフ検索ステップまたはアルゴリズムは、各フレームを通して反射を追跡するために利用される。グラフ検索初期が、前述のテンプレート照合ステップまたはアルゴリズムによって識別される。初期ロックが取得されると、反射線は、その振幅および相対的位置に基づいて、画像を通して追跡される。一実施形態では、線は、一定距離を維持することによって、その振幅および相対的位置に基づいて、画像を通して追跡される。グラフ検索ステップは、第１のＡ走査から開始し、次いで、次のステップのための方向を決定するために、続くＡ走査における近傍画素を見る。本プロセスは、次いで、各Ａ走査に対して繰り返される。アルゴリズムはまた、ステップ方向を決定する前に、次のＡ走査を評価し、次の「ｎ」回のＡ走査を先読みすることを含む「先読み」を可能にする。図３における反射をトレースする黒色線は、グラフ検索アルゴリズムの結果を実証する。図４は、グラフ検索アルゴリズムの若干より詳細な説明を提供する。線が、フレーム全体を通してトレースされると、それらはデジタル的に整列され、フレームが、次いで、表示および測定のために適切に較正される。
【００２４】
一実施形態では、カテーテル反射を較正するための第３のモードは、喪失トラック再取得ステップを含む。喪失トラック再取得ステップは、反射線が前の画像フレーム内で検出されない場合、トラックを再取得することを含む。図５−７に示されるように、反射の追跡は、蛇行性引き戻しの間、喪失され、次いで、反射が再現すると、再取得される。喪失トラックを再取得するために、グラフ検索ステップ２２０およびテンプレート照合アルゴリズムの両方が、使用され得る。グラフ検索ステップは、許容可能な解の領域を拡大し、再現反射に対して、より多くの数のビンを検索する。「保護帯域」が、可能な検索領域を限定するために識別され、脈管壁からの明るいはね返りにロックオンしないように防止する。テンプレート照合ステップも、初期ロックステップにおいて説明されるように行なわれ得る。トラックが再取得されると、アルゴリズムは、再生追跡モードに移行して戻り、各後続フレームを通して、反射を追跡し続ける。
【００２５】
図８Ａを参照すると、第１のモードおよび初期ロックステップ１００の例証方法が、示される。プロセスは、ステップ１１０から開始し、カテーテル内のＯＣＴ画像を撮影する。次いで、ステップ１２０は、テンプレートを利用して、図２Ｂに示されるように、モジュールまたは他のソフトウェアデバイスから、特定のカテーテルに特有のカテーテル反射線の位置を識別する。各カテーテルのためのテンプレートは、特定のカテーテルとともに記憶される、または動作可能にアクセス可能である。一実施形態では、テンプレートは、無線周波数識別（ＲＦＩＤ）チップまたはトランスポンダまたはタグ上に記憶され、代替として、テンプレートは、コンピュータチップ、キャッシュ、フラッシュドライブ、および任意の他の記憶媒体上に記憶され得る。テンプレート照合ステップまたはアルゴリズム１３０は、テンプレートを勾配ＯＣＴ画像のバイナリバージョンとコンボリューションし、ピークが、図２Ａに示されるように、識別される。テンプレート照合ステップは、最も強いテンプレート一致に対応するピークを選択するためのステップ１４０を含む。一実施形態では、ピークは、最も強いテンプレート一致が識別されることを確実にするために、一定の値を上回る。テンプレート一致が、判定１５０において識別されない場合、方法は、ステップ１７０に進み、Ｚ−オフセット位置が、調節され、新しい位置における画像が、テンプレート照合ステップまたはアルゴリズム１３０を用いて、評価される。テンプレート照合アルゴリズムが、強い一致を識別できない場合、ユーザは、警告され、自動較正を再試行するか、または手動で較正するかの選択肢が与えられる。テンプレート一致が、判定１５０において識別される場合、カテーテル反射線の位置が、識別され、カテーテル反射線は、Ｚ−オフセットを調節することによってシフトされ、反射をその所望の場所に移動させ、画像／カテーテルは、ここで、図２Ｂに示されるように、較正される。一実施形態では、Ｚ−オフセットの調節は、可変遅延線（ＶＤＬ）の機械的シフトを含む。
【００２６】
図８Ｂを参照すると、第１のモードの代替方法、初期ロックステップ１００ｂが、示される。プロセス１００ｂは、プロセス１００ａと同様に、ステップ１１０から開始し、カテーテル内のＯＣＴ画像を撮影する。次いで、ステップ１２０は、図２Ｂに示されるように、モジュールまたは他のソフトウェアデバイスからのテンプレートを利用し、特定のカテーテルに特有のカテーテル反射線の位置を識別する。各カテーテルのためのテンプレートは、特定のカテーテルとともに記憶されているか、または動作可能にアクセス可能である。一実施形態では、テンプレートは、無線周波数識別（ＲＦＩＤ）チップまたはトランスポンダまたはタグ上に記憶され、代替として、テンプレートは、コンピュータチップ、キャッシュ、フラッシュドライブ、および任意の他の記憶媒体上に記憶され得る。テンプレート照合ステップ１３０は、テンプレートを勾配ＯＣＴ画像のバイナリバージョンとコンボリューションすることにより、ピークが、図２Ａに示されるように、識別される。テンプレート照合ステップは、最も強いテンプレート一致に対応するピークを選択するためのステップ１４０を含む。一実施形態では、ピークは、最も強いテンプレート一致が識別されることを確実にするために、一定の値を上回る。テンプレート一致が、判定１５０において識別されない場合、方法は、ステップ１７０に進み、Ｚ−オフセット位置が、調節され、新しい位置における画像が、テンプレート照合ステップ１３０を用いて評価される。テンプレート照合アルゴリズムが、強いテンプレート一致を識別できない場合、ユーザは、警告され、自動較正を再試行するか、または手動で較正するかの選択肢が与えられる。テンプレート一致が、判定１５０において識別される場合、カテーテル反射線の位置は、識別され、カテーテル反射線は、Ｚ−オフセットを調節することによってシフトされ、反射をその所望の場所に移動させ、画像／カテーテルは、ここで、図２Ｂに示されるように、較正される。
【００２７】
図９を参照すると、第２のモードおよび再生追跡２００の例証方法が、示される。再生追跡方法２００は、ステップ２１０から開始し、カテーテル引き戻しまたは押し進めを記録する。再生モード追跡の目的は、画像データをデジタル的に整列させ、画像取得中の反射位置に対するいかなるシフトも補正されることを確実にし、適切な分析および／または測定を可能にすることである。記録されたデータセットを通した反射の追跡は、引き戻しの間の動きにより、若干、困難であり、反射の位置は、単一フレームにわたって、有意に変動し得る。図３は、蛇行性引き戻しの間のカテーテル反射の有意なシフトの実施例を実証する。再生追跡方法２００は、グラフ検索アルゴリズムを利用して、各画像フレームを通して反射を追跡する。ステップ２２０におけるグラフ検索アルゴリズムの開始に先立って、反射の初期位置が、前述のテンプレート照合アルゴリズム１３０を使用して、識別される。初期ロック２３０が取得されると、ステップ２４０は、その振幅および相対的位置に基づいて、画像を通して反射線を追跡する。一実施形態では、線は、一定距離を維持することによって、線の振幅および相対的位置に基づいて、画像を通して追跡される。ステップ２５０は、第１のＡ走査から開始し、次いで、続くＡ走査内の近傍画素を見、次のステップのための方向を決定し、各Ａ走査に対して繰り返される。ステップ２６０は、ステップ方向を決定する前に、次のＡ走査を評価し、次の「ｎ」回のＡ走査を先読みすることを含む「先読み」を可能にする。ステップ２７０は、ステップ２５０および２６０からの情報を使用して、現在のＡ走査における反射の方向を決定する。２８０では、アルゴリズムは、次のＡ走査にインクリメントし、再生における全Ａ走査が評価されるまで、同一の処理を繰り返す。図３における反射をトレースする黒色線は、グラフ検索アルゴリズムの結果を実証する。図４は、グラフ検索アルゴリズムの若干より詳細な説明を提供する。ステップ２９０は、フレーム全体を通して反射線をトレースし、それらは、デジタル的に整列され、フレームが、表示および測定のために適切に較正される。
【００２８】
図１０を参照すると、第３のモードおよび喪失トラック再取得方法３００の例証方法が、示される。喪失トラック再取得方法３００は、反射線が前の画像フレーム内で検出されない場合、ステップ３１０から開始する。一実施形態では、反射線は、カテーテルの蛇行性引き戻しの間に喪失され得る。反射が、ステップ３２０において再現すると、喪失トラックは、前述のように、グラフ検索ステップ２２０およびテンプレート照合アルゴリズム１２０を用いて、再取得され得る。グラフ検索ステップ２２０は、ステップ３３０において、許容可能な解の領域を拡大し、再現反射に対して、より多くの数のビンを検索する。ステップ３４０では、「保護帯域」が、識別され、可能な検索領域を限定し、次いで、脈管壁からの明るいはね返りへのロックオンが、防止される。テンプレート照合ステップ１２０も、前述の初期ロックステップに説明されるように行なわれ得る。ステップ３５０では、トラックが再取得されると、アルゴリズムは、再生追跡モードに移行して戻り、各後続フレームを通して、反射を追跡し続ける。
【００２９】
図１１を参照すると、再生方法４００において、自動較正の代替実施形態が示される。方法４００は、概して、引き戻しまたは押し進めデータ４０２を取得すること、および閾値画像４０４を得ることを含む。Ｂ走査データ４１２、雑音推定値４１４、または現在のＺ−オフセット４１６等のいくつかの入力４１０が、閾値画像に結合され得る。次に、ステップ４２０は、グラフ検索アルゴリズム４２０であり、グラフ検索アルゴリズム４２０は、許容可能シフトに対して、テンプレートと画素振幅との差異を算出すること４２２と、ｎ個の先読みに対する差異を算出すること４２４と、最小振幅差異を伴うシフト方向を見つけること４２６と、テンプレート場所を更新し、反射振幅を記憶すること４２８とを含む。次に、判定４３０は、全Ａ走査が処理されたかどうか決定する。そうである場合、方法は、ステップ４４０に進み、将来の整列および表示のために全Ａ走査シフトを記憶し、次いで、自動較正のための再生方法は、完了する４８０。さらなる出力４９０は、再生方法４００のためのＡ走査シフト４９２のデータを含み得る。全Ａ走査が、判定４３０において処理されていない場合、ステップ４５０は、反射が予期される振幅値において検出されるかどうか決定する。判定４５２は、ロックが喪失されているかどうか決定する。ロックが喪失されていない場合、ステップ４５４は、次のＡ走査およびグラフ検索アルゴリズム４２０へのステップ、および許容可能シフトに対してテンプレートと画素振幅値の差異を算出するステップ４２２に進む。ロックが喪失されている場合、喪失トラック再取得ステップ４６０が、実装される。喪失トラック再取得ステップ４６０は、ステップ４６２から開始し、保護帯域領域を算出し、次いで、前述のように、ステップ４６４において、テンプレート照合アルゴリズムのグラフ検索ステップをより大きな領域に適用する。次に、判定４６６は、ロックが再取得されたかどうか決定する。ロックが再取得された場合、グラフ検索アルゴリズム４２０が、始動され、ステップ４２２は、許容可能シフトに対して、テンプレートと画素振幅値の差異を算出する。ロックが再取得されない場合、判定４６８は、ロック喪失時間を超えたかどうか決定する。ロック喪失時間を超えた場合、ステップ４７２は、ユーザに、較正が完了されることができないことを警告し、自動較正の再生方法は、ステップ４８０において、ユーザによって再始動または調節されるために完了する。ロック喪失時間を超えていない場合、ステップ４７０は、次の画像に進み、ロックの取得を再試行し、トラック再取得ステップ４６０のステップ４６２に進み、再び、保護帯域領域を算出する。
【００３０】
図１２を参照すると、自動較正初期ロック方法の代替実施形態５００が、示される。方法５００は、概して、オンモード画像を選択することを含み、決定５０４に進み、第１の画像がオンモードであるかどうか決定する。画像の第１のインスタンスがカテーテルに対してオンモードではない場合、ステップ５０６は、Ｚ−オフセット位置をシフトさせずに、ロックの再取得に進む。画像の第１のインスタンスがオンモードである場合、前述のように、テンプレート照合ステップ５１０に進む。テンプレート照合ステップ５０１は、ステップ５１２から開始し、画像をバイナリＢ走査に変換し、ステップ５１４に進み、ＸおよびＹ勾配を算出し、ステップ５１６に進み、勾配をテンプレートとコンボリューション（順方向およびＣＣ）し、ステップ５１８において、ピーク５１８を見つける。テンプレート照合ステップ５１０は、終了し、判定５２０に進み、ピーク閾値が得られたかどうか決定する。ピーク閾値が得られている場合、ステップ５２２は、領域内の信号対雑音比閾値を伴うピークを見つける。ピーク閾値が得られていない場合、判定５２４は、全Ｚ−オフセット位置が評価されたかどうか決定する。ステップ５２２は、決定５３０に進み、少なくとも２つのピークが見つかったかどうか決定する。少なくとも２つのピークが見つかっていない場合、判定５２４が、始動される。少なくとも２つのピークが見つかっている場合、ステップ５３２は、Ｚ−オフセットを＋ｎビンにシフトさせる。次いで、ステップ５３４は、Ｚ−オフセット位置が到達されることを待ち、判定５４０に進む。決定５４０は、信号がビン＋ｎに移動されているかどうか決定する。信号がビン＋ｎに移動されている場合、ステップ５４２は、ピークを所望の場所にシフトさせる。信号がビン＋ｎに移動されていない場合、決定５５０は、信号がビン−ｎに移動されているかどうか決定する。信号がビン−ｎに移動されている場合、ステップ５４２に進む。信号がビン−ｎに移動されていない場合、決定ステップ５２６に進み、Ｚ−オフセットをＡ走査サイズの１／４シフトさせる。判定５２４も、全Ｚ−オフセット位置が評価されていない場合、ステップ５２６に進む。全Ｚ−オフセット位置が評価された場合、ステップ５２８は、ユーザに、自動較正に問題が生じ、方法５００を再試行するように警告する。次いで、ステップ５２８は、決定５６０に進み、ユーザが、再び、自動較正方法を要求しているかどうか決定する。ユーザが、自動較正方法を再び要求している場合、ステップ５６２は、Ｚ−オフセットを開始位置にリセットする。ステップ５６２は、次いで、ステップ５６４に進み、Ｚ−オフセット位置が到達されることを待ち、次いで、テンプレート照合ステップ５１０、および画像をバイナリＢ走査に変換するステップ５１２に進む。
【００３１】
ステップ５４２は、ステップ５４４に進み、ユーザに、手動較正の再試行または自動較正の容認のために、自動較正完了を促す。ステップ５４４は、次いで、決定５７０に進み、ユーザが較正をＯＫしているかどうか決定する。ユーザが、較正を容認する場合、初期ロック方法は、ステップ５７２において完了する。初期ロックのための出力５８０は、現在のＺ−オフセット５３２、反射場所５８４、または反射振幅５８６を含む。ユーザが、較正を容認しない場合、ステップ５７４は、グラフィカルユーザインターフェース（ＧＵＩ）を通して、手動較正モードに移行する。次いで、決定５９０は、ユーザに較正を完了させる。ユーザが較正を完了すると、ステップ５９２は、較正がさらに要求されない限りＺ−オフセットを維持する。ユーザが較正を完了しない場合、ステップ５７４は、ユーザによるさらなる試行のために、手動較正モードに移行する。
【００３２】
（Ｚ−オフセットの自動較正方法３：自動テンプレート生成および完全相関）  代替として、テンプレートは、前の方法２におけるように、メモリチップまたはＲＦＩＤ上に記憶されなくてもよく、テンプレートは、方法３に関して以下に説明されるように、初期ロックモードプロセス中、自動的に生成され得る。前の方法は、テンプレート照合のためにバイナリテンプレートを利用したが、方法３は、振幅情報および複素共役信号情報を伴うテンプレートを生成する。振幅情報の利用およびミラー信号の生成は、正確な反射線にロックオンする可能性を増加させる。加えて、方法３ステップまたはアルゴリズムは、初期ロックおよび再生モード中、自動較正を行い、方法３はまた、ライブモード中（カテーテルが、撮像しているが、データを記録していないとき）、較正を維持する。
【００３３】
一実施形態では、第１のモードは、初期ロックであり、前の方法で説明されるように、内部カテーテル反射（ファイバ、ミラー、およびレンズ）を利用して、較正位置に到達するために、必要とされるＶＤＬシフトを識別する。初期ロックＺ−オフセット較正は、反射パターンまたはテンプレートが決定されるステップである。初期Ｚ−オフセット較正において識別されたテンプレートは、各後続較正モードにおいて利用され、反射のシフトを追跡し、要求または実装に応じて、アナログおよびデジタルシフトを適用する。テンプレート領域は、勾配および振幅情報を使用して識別される。テンプレートが識別され、後の使用のために記憶されると、ＶＤＬシフトが適用され、カテーテルは、ライブおよび再生モードの間の較正の準備ができる。この位置における容認可能誤差は、ＯＣＴ画像内のカテーテル反射線の参照パターンにロックオンするための次のモード（ライブ撮像の間のＺ−オフセットの維持）の能力によって決定されるであろう。
【００３４】
図１３−１６は、初期Ｚ−オフセット較正のための代替ステップの概略を提供する。図１３Ａ−１３Ｄは、平均振幅または勾配を使用して、強い反射を検索し、検出されるまでＺ−オフセットを調節する、初期ロックモードにおける第１のステップのグラフである。図１３Ａは、ファイバ、ミラー、レンズ、シースに認められる反射、およびテンプレートに対する反射が存在することを示す。アルゴリズムは、画像のＸおよびＹ勾配を利用して、反射が存在するかどうか決定する。図１３Ｂは、Ｘ勾配の変化であり、図１３Ｃは、Ｙ勾配の変化である。反射が見つからない場合、ＶＤＬは、次の可能なＺ−オフセット場所にシフトされる。強い反射が識別されると、若干の正ＶＤＬシフトが、反射の配向を検証するために適用される。反射が、画像の中心に向かってシフトする場合、反射は、較正のために、適切に配向されている。反射が、外向きにシフトする場合、画像は、複素共役信号であって、大幅な負ＶＤＬシフトが、画像をアンラップ（ｕｎｗｒａｐ）するために適用される。画像が、適切に配向されたと決定された後、反射は、ウィンドウの中心にシフトされ、完全テンプレート領域が可視であることを確実にする。シフトステップの各々の間の任意の時点において、反射線が、検出されなくなる場合、アルゴリズムは、新しいＺ−オフセットを適用し、強いリターンを探すことをやり直す。図１７は、アルゴリズムおよびユーザ相互作用のより詳細なバージョンを提供する。
【００３５】
反射が画像ウィンドウ内の中心に置かれると、テンプレートが、算出される。テンプレートは、ファイバ反射から開始し、レンズ反射で終了する、平均振幅値に対する画素番号のアレイである。テンプレートを識別することにおける第１のステップは、図１４Ａに示されるように、単純グラフ検索アルゴリズムを使用して、第１の強い反射に基づいて、画像を整列させることである。長方形画像におけるＡ走査を横断する反射の整列は、初期ロックＺ−オフセット較正モードのための第２のステップである。画像整列は、図１４Ｂに示されるように、反射が、長方形画像内で直線であり、その勾配および振幅によって容易に識別されるであろう可能性を増加させるように行なわれる。画像整列が完了すると、内部カテーテル反射が、４つの特性、すなわち、（１）強いＸおよび弱いＹ勾配（Ｙ軸に沿ったＡ走査を前提）；（２）少なくとも１５ｄＢを上回る一貫した高信号対雑音比（ＳＮＲ）、（３）最初の反射（ファイバ）から最後の反射（レンズ）までの最大距離、および（４）最小でも少なくとも２つの反射を使用して識別される。
【００３６】
図１５Ａ−１５Ｂに示されるように、初期ロックモードにおける第３のステップは、ＸおよびＹ勾配を使用して、カテーテル反射を識別することである。初期ロックのステップ４は、反射のテンプレートを生成し、後の使用のために、テンプレートを記憶する。テンプレート領域は、次いで、図１６Ａに示されるように、最初の反射から開始し、最後の反射で終了するＡ走査（すなわち、角度）（両側に５画素のマージンを伴う）を横断した平均として定義される。弱いレンズ反射による小テンプレートを防止するために、最小テンプレートサイズは、１００ビンである。代替として、最小テンプレートサイズは、少なくとも約２０〜１０００ビンである。したがって、レンズ反射が検出されない場合、ファイバ反射後の１００ビンが、テンプレート領域として選択される。この最小距離閾値は、ファイバからシースの内径（ＩＤ）までの最小距離に基づいて決定され得る。テンプレートの精度は、ファイバ反射を特定することに大きく依存する。ファイバ反射が見つからない場合、テンプレートは、シースまたは脈管領域を含み、不適切な較正をもたらし得る。
【００３７】
図１６Ｂに示されるように、テンプレートが見つかると、後続自動較正アルゴリズムにおいて後に使用するために記憶される。ＶＤＬは、次いで、所定の場所におけるファイバ反射の位置にシフトされ、システムは、ライブモードに移行する。参照パターンが、この初期ＶＤＬ位置において見つからない場合（可能性が低いが、依然として、非ゼロ確率である）、値が正しくないと仮定され、バックアップ初期推定検索が、システムによって行なわれる。バックアップ検索手技は、パターン認識アルゴリズムが、参照パターンにロックオンを試みながら、所定の方式において、ＶＤＬ位置を通して掃引するものである。このバックアップ検索手技に入ることは、追加の時間が、正しい初期較正を見つけるために使用されるので、望ましくない。バックアップ検索が失敗する（すなわち、参照パターンにロックオンできない）場合、システムは、カテーテルまたはＰＩＭへのその接続に不具合があると仮定し、ユーザに、代替カテーテルを使用するように通知する。
【００３８】
図１７は、初期ロックＺ−オフセット較正方法６００の流れ図である。第１のステップは、６９０であり、ユーザは、オンモードの画像を選択する。これが、決定６９２において、オンモードにある第１の画像である場合、アルゴリズムは、６１２に進み、強い反射を検索する。Ｂ走査、雑音推定値、および現在Ｚ−オフセット等の種々の入力６１０が、ステップ６１２と結合され得る。第１の画像がオンモードではない場合、アルゴリズムは、完了ボックス６９４に移動し、単に、ＶＤＬを最後のテンプレートが取得された位置にシフトさせることによって、テンプレートを再取得する。ステップ６１２後、判定６２０は、反射線が見つかったかどうか決定する。反射が見つかった場合、ステップ６２２は、＋Ｚ−シフトのために、反射を左にｎ個のビンだけシフトさせる。一実施形態では、反射のシフトは、約２５〜１００ビンであり得る。反射が見つからなかった場合、ステップ６２４は、新しいＺ−オフセットにシフトする。ステップ６２２後、ステップ６２６は、前述のように、平均振幅または勾配を使用して、ある閾値を満たす強い反射を見つける。次いで、決定６３０は、反射線移動が左に行なわれたかどうか決定する。反射線が移動した場合、ステップ６３２は、特定のビン番号にシフトする。反射線が移動されなかった場合、ステップ６０２は、ＶＤＬを開始場所の複素共役（ＣＣ）にシフトさせ（−Ｚ−シフト）、強い反射６１２を見つけるためにリセットする。ステップ６３２後、ステップ６３４は、強い反射を見つけ、判定６４０に進み、線が特定のビンにおいて検出されたか決定する。線が特定のビンにおいて検出された場合、ステップ６４２は、テンプレートを算出する。線が検出されなかった場合、判定６５０は、全Ｚ−オフセットが試行されたかどうか決定する。全Ｚ−オフセットが試行されていない場合、ステップ６２４は、新しいＺ−オフセット位置にシフトし、ステップ６１２は、強い反射を見つける。全Ｚ−オフセットが試行された場合、ステップ６５２は、ユーザに、プログラムを較正できないことを警告する。ステップ６５２後、判定６５４は、ユーザに、再試行するように警告する。ユーザが、再試行を選択する場合、ステップ６２４は、新しいＺ−オフセットにシフトし、ステップ６１２は、強い反射を見つける。ユーザが、再試行を選択しない場合、ステップ６６０は、手動較正を可能にし、ユーザに、手動較正モードを警告する。ステップ６６０後、結果６６２は、初期較正およびライブ較正モードの完全移行をもたらす。出力６７０は、最終Ｚ−オフセットおよびテンプレートに対して提供され得る。
【００３９】
ステップ６４２後、判定６４４は、少なくとも２つの反射が見つかるかどうか決定する。少なくとも２つの反射が見つかる場合、ステップ６４６は、較正された場所にシフトする。少なくとも２つの反射が見つからない場合、決定６５０は、全Ｚ−オフセットが試行されたかどうか決定を試みる。一実施形態では、タイマ６８０が、決定６８０と結合され、自動較正時間を超えたかどうか決定し得る。自動較正時間を超えた場合、ステップ６８４は、ユーザに、プログラムが較正できないことを警告する。ステップ６８４後、決定６８６は、ユーザが再試行することを可能にする。ユーザが再試行を選択する場合、ステップ６２４は、新しいＺ−オフセットにシフトし、強い反射を見つける。ユーザが再試行を選択しない場合、ステップ６６０は、手動較正を可能にする。
【００４０】
第２のモードの代替実施形態では、ライブモード追跡ステップが、採用され得る。ライブモード自動較正中、初期ロックステップ中に算出されたテンプレートが、ビデオモニタまたは他の表示デバイス上の画面上に表示される全フレームに対して初期ロック較正位置を維持するために利用される。一実施形態では、全フレームに対する初期ロック較正位置は、少なくとも３０フレーム／秒（ｆｐｓ）、代替として、約１０〜５０ｆｐｓのレートであり得る。カテーテルシステムは、カテーテルが体内に挿入されたときの温度変化またはカテーテルが縦方向に押し引きされたときのファイバにかかる機械的歪みによって生じる光学経路長のシフトにより、非較正になり得る。ライブモードアルゴリズムは、テンプレートを使用して、カテーテル反射の位置を検出し、デジタルおよびアナログ較正設定を更新し、適切な較正設定を維持する。わずかなシフトのみ、較正位置を維持するために必要である場合、デジタルシフトが、表示に先立って、画像に適用される。しかしながら、システムが、有意に較正されないままである、または大きなシフトが、較正を維持するために必要である場合、Ｚ−オフセット更新が、適用される（ＶＤＬシフト）。
【００４１】
反射は、テンプレートとＡ走査との間の最大相関を見つけること（すなわち、テンプレート照合）によって、ライブモード追跡の間に識別される。反射を識別するための検索領域は、フレーム毎の最大予測シフトに基づいて限定される。テンプレート照合アルゴリズムは、ラップされた複素共役信号を考慮するために、検索位置に基づいて、テンプレートを修正するので、大部分の標準的テンプレート照合実装と若干異なる。相関の算出に先立って、図１８Ａ−１８Ｄに示されるように、ミラー複素共役信号を含む、「完全テンプレート」が、生成される。完全テンプレートを算出するために、最初に、オリジナルテンプレートが、図１８Ａに示されるように、検索位置にシフトされ、次に、信号が、図１８Ｂに示されるように、同一の信号のミラーバージョンと加算される。完全テンプレートおよび各Ａ走査の相関係数が、次いで、図１９に示されるように、算出される。本プロセスは、次いで、検索領域内の可能なシフト位置毎に繰り返される。全相関が算出されると、各Ａ走査に対する最大相関の位置が、図２０に示されるように、見つけられる。最終反射位置は、上位「ｎ」個の相関の中央値位置として割り当てられる。較正シフトが識別されると、画像は、デジタル的に調節され、図２１Ａ−２１Ｂに示されるように、その較正された位置に反射を戻す。シフトが、「ｎ」個のフレームに対して事前に決定された閾値を超える場合、Ｚ−オフセット（ＶＤＬ）への更新が、適用される。本アルゴリズムは、ライブモードにおいて、画面上に表示される全ての画像に繰り返され、前のフレームの相関一致およびＺ−オフセットを利用して、次のフレームのための検索領域を決定する。
【００４２】
図２２は、ライブモード追跡プロセス７００のためのアルゴリズムおよびユーザ相互作用の流れ図を提供する。ライブモード較正プロセスでは、較正は、「画像オフ」が選択されるか、またはカテーテル縦方向引き戻しが、始動されるまで、継続する。引き戻し開始直前の反射の位置が、以下に説明されるように、再生モード自動較正設定において使用するために記憶される。前述のように、Ｂ走査、現在のＺ−オフセット等の種々の入力７１０が、ライブモード追跡プロセスと結合され得る。ステップ７１２は、初期ロックの自動較正が完了したかどうか決定する。次いで、ステップ７１４は、全許容可能シフト位置に対して完全テンプレートを算出する。次いで、ステップ７１６は、Ａ走査およびテンプレートのサブセットに対する相関を算出する。次いで、ステップ７１８は、走査毎の最大相関を見つける。次いで、ステップ７２０は、上位ｎ個の相関の中央値シフトを見つける。次いで、判定７３０は、相関が特定の閾値を上回るかどうか決定する。相関が特定の閾値を上回る場合、ステップ７３２は、移動平均に基づいて、相関閾値を算出する。相関が特定の閾値を上回らない場合、ステップ７３４において、インクリメンタルロック喪失カウンタが進められる。ステップ７３４後、判定７４０は、ロック喪失カウンタ閾値を超えたかどうか決定する。ロック喪失カウンタ閾値を超える場合、判定７４２は、ユーザが自動較正を「オン」として選択したかどうかチェックし、ユーザにエラーを警告する必要があるかどうか決定する。ロック喪失カウンタ閾値を超えていない場合、ステップ７４４は、次の画像にインクリメントし、その後、ステップ７１４が続き、ライブモード追跡のための全許容可能シフト位置に対して完全テンプレートを算出する。判定７４２では、自動較正が、ユーザによって「オン」に選択される場合、判定７５０は、ロック喪失カウンタが閾値を１上回るかどうか決定する。自動較正が、ユーザによって「オン」に選択されていない場合、ステップ７４４は、次の画像にインクリメントするステップに進む。ロック喪失カウンタ閾値を１上回らない場合、ステップ７４４は、次の画像にインクリメントする。ロック喪失カウンタ閾値が閾値を１上回る場合、ステップ７５２は、ユーザに、自動較正ロックが喪失されたことを警告し、それによって、プログラムは、コンピュータのウィンドウまたは状態バーにフェードインすることができる。
【００４３】
ステップ７３２後、ステップ７３６は、ロック喪失を０にリセットする。次いで、判定７６０は、自動較正がユーザによって「オン」に選択されたかどうか決定する。自動較正がユーザによって「オン」に選択されている場合、ステップ７６２は、表示のために、デジタルシフトを現在の画像に適用する。自動較正がユーザによって「オン」に選択されていない場合、ステップ７６４は、次の画像のために、反射位置を更新する。ステップ７６４後、ステップ７４４は、プログラムを次の画像にインクリメントする。ステップ７６２後、判定７７０は、デジタルシフト閾値が満たされたかどうか決定する。デジタルシフト閾値が満たされた場合、ステップ７７２は、インクリメンタルデジタルシフトカウンタを進める。デジタルシフト閾値が満たされていない場合、ステップ７７４は、デジタルシフトカウンタを０にリセットし、次いで、ステップ７８４に進み、次の画像のために反射位置を更新する。ステップ７７２後、判定７８０は、デジタルシフトカウンタ閾値を超えたかどうか決定する。デジタルシフトカウンタ閾値を超えた場合、ステップ７８２は、ＶＤＬシフトを適用する。デジタルシフトカウンタ閾値を超えていない場合、ステップ７４４は、ライブモード追跡プロセスのために、次の画像にインクリメントする。ライブモード較正プロセスでは、較正は、「画像オフ」が選択されるか、またはカテーテル縦方向引き戻しが始動されるまで、継続する。カテーテル引き戻し開始直前の反射位置は、以下に説明されるように、再生モード自動較正設定において使用するために記憶される。
【００４４】
第３のモードの代替実施形態では、再生モード追跡が、ユーザが画像データセットを記録した後に生じる。再生モード追跡は、データセット内の全てのＡ走査に自動較正を行なう。ライブモード追跡と同様に、再生モードは、限定されたシフト場所におけるテンプレートおよび画像Ａ走査の相関を利用して、反射の位置を決定する。反射の初期位置の識別は、データセットの第１のフレームが、他のフレームと異なって処理されるようなものである。データセットの第１のフレームでは、全許容可能シフトおよび全Ａ走査に対する相関が、図２３Ａ−２３Ｂに示されるように、最大相関を見つけるために算出される。最大相関点から、アルゴリズムは、次いで、各Ａ後方および前方走査を通してトレースし、図２４に示されるように、各可能なシフト毎の相関を算出する。第１の検索のための許容可能シフト領域は、ライブモードから再生モードへの移行の間に生じ得る、突然のジャンプを可能にするために広範である。開始位置が決定されると、Ａ走査毎の検索は、Ａ走査間の時間および可能な移動は、フレーム毎の移動に対して小さいことを前提として、小領域に限定される。例えば、第１のフレームに対して、検索領域は、最後の場所から−５０〜＋５０画素に設定され得、最大相関が見つかる場合、検索領域は、Ａ走査毎に−１〜＋１画素に限定される。代替として、検索領域は、少なくとも約−５００〜＋５００画素、代替として、少なくとも約−４００〜＋４００、代替として、約−３００〜＋３００、等に設定され得る。検索領域は、約−１０〜＋１０、代替として、約−５〜＋５、代替として、約−０．１〜＋０．１に限定され得る。第１のフレームが、図２５Ａに示されるように、完全にトレースされると、アルゴリズムは、次のフレームに移動し、第１のＡ走査から開始し、最後のＡ走査内の反射位置に基づいて、検索領域を限定する。図２５Ｂは、反射位置が、ビデオモニタまたは他の表示デバイスの画面上での画像の表示に先立って、整列のために、テンプレート一致位置を記憶することを通して、後の整列および表示のために記憶されることを示す。これは、全フレーム内の各Ａ走査に対して繰り返される。
【００４５】
再生モード較正プロセス８００の詳細な流れ図は、図２６に提供される。再生モード較正は、画像内の全Ａ走査を検索し、ピーク相関を識別することによって始動する。ピークから、相関トラッカが、前方および後方に追跡し、第１のフレーム内のＡ走査毎に反射位置を識別する。この検索は、第１のフレームに適用され、強い初期ロックを保証する。第１のフレーム後の続くフレームの各々は、限定された検索領域を用いて、Ａ走査毎に追跡される。手動モード８４０は、ロック喪失カウンタ閾値を超え、ユーザが手動モードを選択した場合、移行される。代替として、手動モード８４０は、押し引きデータが記録された場合、選択され得る８４２。手動モードが選択された場合、再生モード自動較正は、ステップ８４４において適用されないであろう。手動モードが選択されていない場合、ステップ８５０は、カテーテルの押し引きに基づいて、第１のフレームのための許容可能シフトを識別する。次いで、ステップ８５２は、全許容可能シフト位置に対する完全テンプレートを算出する。次いで、ステップ８５４は、全Ａ走査およびテンプレート位置に対して相関を算出する。次いで、ステップ８５６は、Ａ走査毎に最大相関を見つける。次いで、ステップ８５８は、最大相関および対応するＡ走査を見つける。次いで、ステップ８６０は、相関追跡アルゴリズム８１０を画像内の各Ａ走査に適用する。
【００４６】
相関追跡アルゴリズム８１０は、現在のＡ走査内の許容可能テンプレートシフトに対する相関を算出するステップ８１２から開始する。次いで、ステップ８１４は、そのＡ走査に対する最大相関を見つける。次いで、ステップ８１６は、移動平均の相関信頼閾値を算出する。次いで、判定８２０は、相関が特定の信頼閾値を上回るかどうか決定する。相関が特定の信頼閾値を上回る場合、ステップ８２２は、ロック喪失カウンタを０にリセットする。相関が特定の信頼閾値を上回らない場合、ステップ８２４は、ロック喪失カウンタをインクリメントするステップに進む。ステップ８２４後、判定８３０は、ロック喪失カウンタ閾値を超えるかどうか決定する。ロック喪失カウンタ閾値を上回らない場合、ステップ８３２は、前のＡ走査のトラック位置を使用する。ロック喪失カウンタ閾値を上回る場合、ステップ８３４は、ユーザに警告し、手動モードに移行する。ステップ８２２および８３２は両方とも、ステップ８３６に進み、トラック位置を更新し、次のＡ走査に進む。ステップ８３６後、ステップ８１２は、現在のＡ走査内の許容可能シフトに対する相関を算出する。
【００４７】
相関追跡アルゴリズムを画像内の各Ａ走査に適用するステップ８６４後、判定８７０は、最後のフレームであるかどうか決定する。最後のフレームである場合、ステップ８７２は、表示のために、較正位置を記憶する。さらにフレームが存在する場合、判定８８０は、手動モードへの移行が必要とされるか、または命令されるかを決定する。手動モードへの移行が選択される場合、ステップ８６２は、次のフレームにインクリメントする。手動モードへの移行が選択されない場合、ステップ８８２は、ユーザに警告し、手動モードに移行する。ステップ８７２後、ステップ８７４は、再生モード較正が完了したかどうか決定する。Ｂ走査、現在のＺ−オフセット、およびカテーテルのための押し引きインジケータ等の種々の入力８９０が、再生モードプロセスと結合され得る。
【００４８】
概して、自動較正９００の一実施形態が、図２７に示される。前述の較正方法のいずれも、初期較正後、フレーム毎に継続的に較正を更新し、それを維持するために使用され得る。ステップ９０２は、前述のように、初期自動較正または手動較正を行なう。ステップ９０４は、較正位置を示す、少なくとも１つのパラメータを監視する。判定９０６は、較正が既存のフレームまたは後続フレームにおいて更新される必要があるかどうか決定する。較正がフレームに対して更新される必要がない場合、自動較正は、ステップ９０４において、較正位置を示すパラメータの監視を継続する。較正がフレームに対して更新される必要がある場合、ステップ９０８は、自動的に、較正を更新する（画像をデジタル的にシフトさせる、ｚ−オフセットシフトを適用する、または前述の方法のいずれか等）。フレームは、Ａ走査またはフレームのセットであり得る。
【００４９】
流れ図例証の各ブロックおよび流れ図例証におけるブロックの組み合わせは、本明細書に開示されるモジュール、システム、および方法の任意の部分と同様に、コンピュータプログラム命令によって実装されることができることを理解されるであろう。これらのプログラム命令は、プロセッサに提供され、プロセッサ上で実行する命令が、流れ図ブロックまたは複数のブロックにおいて規定され、あるいは本明細書に開示される組織分類器、撮像機、制御モジュール、システム、および方法に関して説明される作用を実装するための手段を作成するような機械をもたらし得る。コンピュータプログラム命令は、プロセッサによって実行され、一連の動作ステップをプロセッサによって行なわせ、コンピュータ実装プロセスをもたらし得る。コンピュータプログラム命令はまた、動作ステップのうちの少なくともいくつかを並行して行なわせ得る。さらに、ステップのうちのいくつかはまた、マルチプロセッサコンピュータシステム内に生じ得るように、２つ以上のプロセッサにわたって行なわれ得る。加えて、１つ以上のプロセスはまた、本発明の範囲または精神から逸脱することなく、他のプロセスと並行して、または図示されるものと異なるシステム内でさえ、行なわれ得る。
【００５０】
コンピュータプログラム命令は、限定されないが、ＲＡＭ、ＲＯＭ、ＥＥＰＲＯＭ、フラッシュメモリまたは他のメモリ技術、ＣＤ−ＲＯＭ、デジタル多用途ディスク（ＤＶＤ）または他の光学記憶、磁気カセット、磁気テープ、磁気ディスク記憶または他の磁気記憶デバイス、あるいは所望の情報を記憶するために使用されることができ、かつコンピューティングデバイスによってアクセスされることができる、任意の他の媒体を含む、任意の好適なコンピュータ読み取り可能な媒体上に記憶されることができる。
【００５１】
カテーテル引き戻しは、カテーテルを撮像されている領域の近位端から遠位端に引っ張ることによって行なわれ得ることを理解されるであろう。また、前述の血管内撮像技法はまた、患者血管系内に挿入可能なカテーテルを使用する他のタイプの撮像技法とともに使用されることができることを理解されるであろう。例えば、血管内撮像技法は、患者組織の１つ以上の測定可能特性（例えば、血管内磁気共鳴画像診断、分光法、温度マッピング等）を査定するために構成および配列される任意の撮像技法とともに使用されることができる。
【００５２】
本明細書に説明されるシステムおよび方法は、多くの異なる形態で具現化され得、本明細書に記載される実施形態に限定されるものと解釈されるべきではない。故に、開示されるシステムおよび方法は、全体的にハードウェア実施形態、全体的にソフトウェア実施形態、またはソフトウェアおよびハードウェア側面を組み合わせる実施形態の形態をとってもよい。本明細書に説明されるシステムおよび使用方法は、各デバイスがプロセスまたは方法の少なくとも一部を行なう、プロセッサまたはコンピューティングデバイスの任意の組み合わせを含む、コンピュータ等の任意のタイプのコンピューティングデバイスを使用して行なわれることができる。
【００５３】
好適なコンピューティングデバイスは、典型的には、マスメモリを含み、典型的には、デバイス間の通信を含む。マスメモリは、コンピュータ読み取り可能な媒体、すなわち、コンピュータ記憶媒体のタイプを表す。コンピュータ記憶媒体として、コンピュータ読み取り可能な命令、データ構造、プログラムモジュール、または他のデータ等の情報の記憶のための任意の方法または技術において実装される、揮発性、不揮発性、除去可能、および非除去可能媒体が挙げられ得る。コンピュータ記憶媒体の実施例として、ＲＡＭ、ＲＯＭ、ＥＥＰＲＯＭ、フラッシュメモリ、または他のメモリ技術、ＣＤ−ＲＯＭ、デジタル多用途ディスク（ＤＶＤ）または他の光学記憶、磁気カセット、磁気テープ、磁気ディスク記憶または他の磁気記憶デバイス、無線周波数識別タグまたはチップ、あるいは所望の情報を記憶するために使用されることができ、かつコンピューティングデバイスによってアクセスされることができる、任意の他の媒体が挙げられる。デバイスまたはシステムの構成間の通信は、有線および無線（例えば、ＲＦ、光学、または赤外線）通信の両方を含むことができる。
【００５４】
本発明は、種々の実施形態に関連して説明されたが、本発明は、さらなる修正が可能であることを理解されるであろう。本願は、本発明の一般的原理に従い、本発明が属する分野において既知および慣行されている範囲内での本開示からの逸脱を含む、本発明のいかなる変形例、使用、または適応も網羅することを意図する。
本発明は一側面において以下の発明を包含する。
（発明１）
カテーテルを較正する方法であって、
カテーテル内の画像を撮影することと、
テンプレートを利用して、前記カテーテル反射線の位置を識別することと
を含む、方法。
（発明２）
前記テンプレートは、前記特定のカテーテルとともに記憶される、または動作可能にアクセス可能である、発明１に記載の方法。
（発明３）
前記テンプレートは、前記カテーテルの画像から直接導出される、発明１に記載の方法。
（発明４）
前記テンプレートを前記画像のあるバージョンと照合し、ピークを識別することをさらに含む、発明１に記載の方法。
（発明５）
前記画像のバージョンは、バイナリである、発明４に記載の方法。
（発明６）
最も強いテンプレート一致に対応する前記ピークを選択することをさらに含む、発明４に記載の方法。
（発明７）
前記ピークは、一定の値を上回る、発明６に記載の方法。
（発明８）
テンプレート一致が識別されない場合、前記Ｚ−オフセット位置を調節し、前記テンプレート照合ステップを用いて評価されるべき新しい位置において、前記カテーテルの第２の画像を撮影することをさらに含む、発明１に記載の方法。
（発明９）
前記カテーテル反射線の位置を識別することと、テンプレート一致が識別された場合、前記Ｚ−オフセットを調節することによって、前記カテーテル反射線をシフトさせることとをさらに含む、発明４に記載の方法。
（発明１０）
前記Ｚ−オフセットを調節することは、可変遅延線の機械的シフトを含む、発明６に記載の方法。
（発明１１）
カテーテルを較正する方法であって、
カテーテル引き戻しを撮像することと、
画像取得中の反射位置に対するいかなるシフトも補正されることを確実にするために、画像データを整列させることと、
各画像フレームを通して、前記反射位置を追跡することと
を含む、方法。
（発明１２）
前記反射位置を追跡するステップは、グラフ検索アルゴリズムを利用することを含む、発明１１に記載の方法。
（発明１３）
前記反射位置を追跡するステップは、テンプレート照合ステップを通して、前記反射の初期位置を識別する、発明１１に記載の方法。
（発明１４）
前記反射位置を追跡するステップは、初期ロックが取得されると、振幅および相対的位置に基づいて、前記画像を通して前記反射線を追跡することをさらに含む、発明１３に記載の方法。
（発明１５）
前記線は、一定距離を維持することによって、前記線の振幅および相対的位置に基づいて、前記画像を通して追跡される、発明１４に記載の方法。
（発明１６）
第１のＡ走査時、次のステップの方向を決定するために、近傍画素が続くＡ走査において検討され、各Ａ走査に対して繰り返される、発明１５に記載の方法。
（発明１７）
前記ステップ方向を決定する前に、前記次のＡ走査を評価することと、次の「ｎ」個のＡ走査を先読みすることとを含む先読みステップをさらに含む、発明１６に記載の方法。
（発明１８）
フレーム全体を通して前記反射線をトレースすることと、前記反射線をデジタル的に整列させることとをさらに含む、発明１７に記載の方法。
（発明１９）
カテーテルを較正する方法であって、
カテーテルの画像を撮影することと、
前記カテーテルの反射線を追跡することと、
画像において前記反射線を検出しないことと、
グラフ検索ステップまたはテンプレート照合ステップを用いて、喪失トラック画像を再取得することと
を含む、方法。
（発明２０）
前記グラフ検索ステップは、許容可能な解の領域を拡大し、再現反射に対して、より多くの数のビンを検索する、発明１９に記載の方法。
（発明２１）
可能な検索領域を限定する保護帯域を識別し、脈管壁からの明るいはね返りへのロッオンを防止することをさらに含む、発明２０に記載の方法。
（発明２２）
前記テンプレート照合ステップは、前記喪失トラック画像を再取得し、前記カテーテル反射を追跡することを各後続フレームを通して継続するために移行して戻る、発明１９に記載の方法。
（発明２３）
深度較正および自動化された範囲調節の方法であって、
静止しているが、制御不能変位によるものではない画像特徴を検索することによって、非制御変位効果を認識することと、
極座標表現が正確な次元測定のために使用され得るように、連続した画像データを較正することと
を含む、方法。
（発明２４）
表示に先立って、前記画像内の画像特徴を除去することをさらに含む、発明２３に記載の方法。
（発明２５）
初期較正モードのための方法であって、
内部カテーテル反射を利用して、較正位置に到達するために必要とされる可変遅延線シフトを識別することと、
前記カテーテル反射の初期テンプレートを決定することと、
各後続較正モードにおいて、前記初期テンプレートを利用し、前記反射のシフトを追跡することと、
前記適切なシフトを適用することと、
を含む、方法。
（発明２６）
平均振幅または勾配を使用して、強い反射を検索することと、前記強い反射が検出されるまで、Ｚ−オフセットを調節することとをさらに含む、発明２６に記載の方法。
（発明２７）
画像を横断して、前記強い反射を整列させることをさらに含む、発明２７に記載の方法。
（発明２８）
前記勾配および振幅の変化を使用して、前記カテーテル反射を識別することと、閾値勾配および振幅情報を識別し、前記テンプレート領域を決定することとをさらに含む、発明２７に記載の方法。
（発明２９）
反射のテンプレートを生成し、他の自動較正モードにおける後の使用のために記憶することをさらに含む、発明２８に記載の方法。
（発明３０）
全画像フレームに対して、初期ロック較正位置を維持することをさらに含む、発明２９に記載の方法。
（発明３１）
前記反射のテンプレートを使用して、前記カテーテル反射の位置を検出し、較正設定を更新し、適切な較正設定を維持することをさらに含む、発明３０に記載の方法。
（発明３２）
前記反射のテンプレートをシフトさせ、場所を検索し、ミラー信号を伴う完全テンプレートを作成することをさらに含む、発明３１に記載の方法。
（発明３３）
前記テンプレートおよび全Ａ走査に対する相関係数を算出することをさらに含む、発明３２に記載の方法。
（発明３４）
各Ａ走査に対する最大相関を見つけ、上位Ａ走査の中央値を求めることをさらに含む、発明３３に記載の方法。
（発明３５）
デジタルシフトを画像および走査変換に適用することをさらに含む、発明３４に記載の方法。
（発明３６）
再生モード追跡較正のための方法であって、
第１のフレームに対するカテーテルの反射の初期位置を識別することと、
前記第１のフレームに対する全許容可能シフトおよび全Ａ走査に対する相関を算出し、最大相関を識別することと、
後続のＡ走査に対して前記最大相関を後方および前方に利用して、各可能なシフトに対する相関を算出することと、
データセット内の全フレームを通して、前記Ａ走査毎に追跡し、後の整列および表示のために、前記反射位置を記憶することと、
を含む、方法。
（発明３７）
自動較正のための方法であって、
較正位置に対してカテーテルの初期較正を行なうことと、
前記較正位置を示す少なくとも１つのパラメータを監視することと、
前記較正位置が画像フレームでの更新を必要とするかどうか決定することと、
前記較正位置を更新することと
を含む、方法。
（発明３８）
前記画像フレームは、Ａ走査である、発明３７に記載の方法。
（発明３９）
前記画像フレームは、一式の画像フレームを含む、発明３７に記載の方法。
【図１】
0006475496.tif 000002
【図２Ａ】
0006475496.tif 000003
【図２Ｂ】
0006475496.tif 000004
【図３】
0006475496.tif 000005
【図４】
0006475496.tif 000006
【図５】
0006475496.tif 000007
【図６】
0006475496.tif 000008
【図７】
0006475496.tif 000009
【図８Ａ】
0006475496.tif 000010
【図８Ｂ】
0006475496.tif 000011
【図９】
0006475496.tif 000012
【図１０】
0006475496.tif 000013
【図１１】
0006475496.tif 000014
【図１２】
0006475496.tif 000015
【図１３Ａ】
0006475496.tif 000016
【図１３Ｂ】
0006475496.tif 000017
【図１３Ｃ】
0006475496.tif 000018
【図１３Ｄ】
0006475496.tif 000019
【図１４】
0006475496.tif 000020
【図１５】
0006475496.tif 000021
【図１６】
0006475496.tif 000022
【図１７】
0006475496.tif 000023
【図１８】
0006475496.tif 000024
【図１９】
0006475496.tif 000025
【図２０】
0006475496.tif 000026
【図２１】
0006475496.tif 000027
【図２２】
0006475496.tif 000028
【図２３】
0006475496.tif 000029
【図２４】
0006475496.tif 000030
【図２５】
0006475496.tif 000031
【図２６】
0006475496.tif 000032
【図２７】
0006475496.tif 000033