(57)【要約】【課題】  信号にアノテーション付けするための方法を提供する。【解決手段】  ヒト被験体の心筋に近接する一対の電極からバイポーラ信号を受信すること、及び一対の電極のうちの選択された一方からユニポーラ信号を受信すること、を含む、方法。方法は、ユニポーラ信号の局所ユニポーラ最小導関数、及びユニポーラ最小導関数の発生時刻を計算することを更に含む。方法は、バイポーラ信号のバイポーラ導関数を計算すること、バイポーラ導関数と局所ユニポーラ最小導関数との比を求めること、及び比が予め設定された閾値の比の値を超える場合は、心筋の活性化時間として発生時刻を特定すること、を更に含む。【選択図】  図１
2018023787.tif 000002
【特許請求の範囲】
【請求項１】
ヒト被験体の心筋に近接する一対の電極からバイポーラ信号を受信すること、
前記一対の電極のうちの選択された一方からユニポーラ信号を受信すること、
前記ユニポーラ信号の局所ユニポーラ最小導関数、及び前記ユニポーラ最小導関数の発生時刻を計算すること、
前記バイポーラ信号のバイポーラ導関数を計算すること、
前記バイポーラ導関数と前記局所ユニポーラ最小導関数との比を求めること、並びに
前記比が予め設定された閾値の比の値を超える場合は、前記心筋の活性化時間として前記発生時刻を特定すること、を含む、方法。
【請求項２】
前記バイポーラ導関数が予め設定されたバイポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、請求項１に記載の方法。
【請求項３】
前記局所ユニポーラ最小導関数が予め設定されたユニポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、請求項１に記載の方法。
【請求項４】
前記発生時刻に対応する信頼水準が予め設定された信頼水準を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、請求項１に記載の方法。
【請求項５】
前記バイポーラ信号の振幅が予め設定されたバイポーラ信号の閾値を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、請求項４に記載の方法。
【請求項６】
ヒト被験体の心筋に近接して配置されるように構成された一対の電極と、
プロセッサであって、
前記一対の電極からバイポーラ信号を受信し、
前記一対の電極のうちの選択された一方からユニポーラ信号を受信し、
前記ユニポーラ信号の局所ユニポーラ最小導関数、及び前記ユニポーラ最小導関数の発生時刻を計算し、
前記バイポーラ信号のバイポーラ導関数を計算し、
前記バイポーラ導関数と前記局所ユニポーラ最小導関数との比を求め、かつ
前記比が予め設定された閾値の比の値を超える場合は、前記心筋の活性化時間として前記発生時刻を特定するように構成されたプロセッサと、を含む、装置。
【請求項７】
前記バイポーラ導関数が予め設定されたバイポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、請求項６に記載の装置。
【請求項８】
前記局所ユニポーラ最小導関数が予め設定されたユニポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、請求項６に記載の装置。
【請求項９】
前記発生時刻に対応する信頼水準が予め設定された信頼水準を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、請求項６に記載の装置。
【請求項１０】
前記バイポーラ信号の振幅が予め設定されたバイポーラ信号の閾値を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、請求項９に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
（関連出願の相互参照）
本出願は、参照により本明細書に組み込まれる２０１６年８月１１日出願の米国特許仮出願第６２／３７３，４５８号の利益を主張するものである。本出願は、本出願と同日付で出願された「Ｃｌａｓｓｉｆｙｉｎｇ  ＥＣＧ  Ｓｉｇｎａｌｓ」と題された出願に関連する。
【０００２】
（発明の分野）
本発明は、広義には心電計（ＥＣＧ）信号に関し、より詳細には信号にアノテーション付けするための方法に関する。
【背景技術】
【０００３】
心臓の電気信号のマッピング及びイメージングは、典型的には、カテーテルのＥＣＧ信号によって示される局所活性化時間（ＬＡＴ）と信号の空間位置との組み合わせに基づく。こうした方法は、Ｂｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒ（Ｄｉａｍｏｎｄ  Ｂａｒ，Ｃａ）製のＣＡＲＴＯ（登録商標）３システムで用いられている。
【０００４】
既存の活動測定アルゴリズムは、本質的に遠距離場の干渉の影響を受けないバイポーラ信号の処理に基づく。歴史的には、ＬＡＴは基準アノテーションと、マッピングカテーテルチャネル内のバイポーラ信号における最初の活動と、の間の時間差によって測定される。しかしながら、最初の活性化は、２つのユニポーラ信号からの接近場（approaching fields）構成であり、活動波の開始とは必ずしも一致しないため、バイポーラ信号における最初の活性化を測定すると、ＬＡＴの測定値は曖昧になる。
【０００５】
本特許出願に参照により組み込まれる文献は、いずれかの用語がこれらの組み込まれた文献において、本明細書に明確に又は暗示的になされる定義と矛盾する形で定義されている場合には、本明細書中の定義のみが考慮されるべきである点を除いて、本出願の一部とみなされるものとする。
【発明の概要】
【課題を解決するための手段】
【０００６】
本発明の実施形態は、
ヒト被験体の心筋に近接する一対の電極からバイポーラ信号を受信すること、
一対の電極のうちの選択された一方からユニポーラ信号を受信すること、
ユニポーラ信号の局所ユニポーラ最小導関数、及びユニポーラ最小導関数の発生時刻を計算すること、
バイポーラ信号のバイポーラ導関数を計算すること、
バイポーラ導関数と局所ユニポーラ最小導関数との比を求めること、並びに、
比が予め設定された閾値の比の値を超える場合は、心筋の活性化時間として発生時刻を特定すること、を含む、方法を提供する。
【０００７】
開示された実施形態は、バイポーラ導関数が予め設定されたバイポーラ導関数の閾値を下回る場合は、心筋の活性化時間として発生時刻を割り当てることを含む。
【０００８】
更なる開示された実施形態は、局所ユニポーラ最小導関数が予め設定されたユニポーラ導関数の閾値を下回る場合は、心筋の活性化時間として発生時刻を割り当てることを含む。
【０００９】
なお更に開示される実施形態は、発生時刻に対応する信頼水準が予め設定された信頼水準を上回る場合は、心筋の活性化時間として発生時刻を割り当てることのみを含む。典型的には、方法は、バイポーラ信号の振幅が予め設定されたバイポーラ信号の閾値を上回る場合は、心筋の活性化時間としての発生時刻を割り当てることのみを含み得る。
【００１０】
本発明の一実施形態によれば、更に、
ヒト被験体の心筋に近接して配置されるように構成された一対の電極と、
一対の電極からバイポーラ信号を受信し、
一対の電極のうちの選択された一方からユニポーラ信号を受信し、
ユニポーラ信号の局所ユニポーラ最小導関数、及びユニポーラ最小導関数の発生時刻を計算し、
バイポーラ信号のバイポーラ導関数を計算し、
バイポーラ導関数と局所ユニポーラ最小導関数との比を求め、かつ
比が予め設定された閾値の比の値を超える場合は、心筋の活性化時間として発生時刻を特定するように構成されたプロセッサと、を含む、装置が提供される。
【００１１】
以下の本開示の実施形態の詳細な説明を図面と併せ読むことで本開示のより完全な理解が得られるであろう。
【図面の簡単な説明】
【００１２】
【図１】本発明の一実施形態によるアルゴリズムの概略的なブロック図である。
【図２】本発明の一実施形態による、バイポーラ信号、及びユニポーラ正極信号によって測定される活動の例である。
【図３】本発明の一実施形態による、ベースラインのふらつきの除去を図示するグラフである。
【図４】本発明の一実施形態による、ベースラインのふらつき除去システムのブロック図である。
【図５】本発明の一実施形態による、２つのガウスフィルタのグラフである。
【図６】本発明の一実施形態によるアノテーション検出器ブロックの概略的なブロック図である。
【図７】本発明の一実施形態による、ユニポーラ及びバイポーラ信号、並びにこれらの導関数のグラフである。
【図８】本発明の一実施形態による、アノテーションアルゴリズムの第１拒絶フェーズを図示するグラフを示す。
【図９】本発明の一実施形態による、局所及び遠距離場の候補アノテーションを図示するグラフを示す。
【図１０】本発明の一実施形態による、候補アノテーションのマージ、及びどのように拒絶判定基準が用いられるかを図示するグラフである。
【図１１】本発明の一実施形態による、ユニポーラ導関数のファジー関数のグラフである。
【図１２】本発明の一実施形態による、ユニポーラ信号のセグメンテーションを図示するグラフである。
【図１３】本発明の一実施形態による、ユニポーラ持続時間のファジー関数のグラフである。
【図１４】本発明の一実施形態による、ユニポーラ振幅のファジー関数のグラフである。
【図１５】本発明の一実施形態による、ユニポーラ持続時間と振幅との比のファジー関数のグラフである。
【図１６】本発明の一実施形態による、バイポーラ振幅のファジー関数のグラフである。
【図１７】本発明の一実施形態による、装置を用いた侵襲性医療処置の概略図である。
【発明を実施するための形態】
【００１３】
概説
バイポーラ信号とは異なり、ユニポーラ信号における活動開始は、信号の急な下方偏向によって明確にマークされるが、バイポーラ信号とは異なり、ユニポーラは遠距離場活性化に対して敏感である。
【００１４】
本発明の実施形態は、２種類のＥＣＧ信号（バイポーラ信号、及びそれに付随するユニポーラ信号のうちの１つ）の特性を組み合わせるように作用する波面アノテーションアルゴリズムを用いて、正確な信号アノテーションを生成する。本発明者らは、本アルゴリズムが、遠距離場の干渉の影響を受けない正確なアノテーションを提供することを実証した。
【００１５】
波面アノテーションアルゴリズムは、比較的短時間に多数のＬＡＴポイントの捕捉及びアノテーション付けを可能にする自動的で信頼性の高いアノテーションポイントの検出を実現する。
【００１６】
本発明の実施形態は、ヒト被験体の心筋に近接する一対の電極からバイポーラ信号を受信することを含む。一対の電極のうちの選択された一方からユニポーラ信号もまた受信され、ユニポーラ信号の局所ユニポーラ最小導関数が計算される。また、ユニポーラ最小導関数の発生時刻も計算される。バイポーラ信号のバイポーラ導関数が計算され、バイポーラ導関数と局所ユニポーラ最小導関数との比が求められる。
【００１７】
比が予め設定された閾値の比の値を超える場合は、心筋の活性化時間として（ユニポーラ最小導関数の）発生時刻が同定される。
【００１８】
図１は、本発明の一実施形態による波面アノテーションアルゴリズムの概略的なブロック図である。アルゴリズムの入力値は、典型的には、５００Ｈｚをカットオフした低域通過フィルタ及び電力遮断フィルタ（power rejection filter）を辿る、アルゴリズムを操作するプロセッサ２０に提供される単一のバイポーラ信号、及びそのユニポーラ信号のうちの１つからなる。プロセッサ２０の操作の更なる詳細は、下記の図１７を参照して提供される。ユニポーラ信号の極性は、既知であると考えられる（すなわち、これは正極又は陰極のいずれかから導かれる）。プロセッサは、スタンドアロン型プロセッサ、及び／又は典型的にはコンピュータを操作する汎用プロセッサであってもよい。アルゴリズムは、本明細書で要約される多数の段階からなる。
【００１９】
前処理段階２２は、ベースラインのふらつきの除去、低域通過フィルタリング、及び微分の任意の順序を含む。ベースラインのふらつきの除去は、アーチファクトであり、機械的なカテーテルの運動又は呼吸などの様々な理由に起因する付加的な低周波信号の除去を含む。この低周波信号は、推定された信号の導関数を変更させ得るため、典型的には除去される。
【００２０】
特徴抽出段階２４は、処理後信号を用いて、全ての候補アノテーションに関して特徴を抽出する。
【００２１】
第１のアノテーション検出器段階２６は、特徴のサブセットに基づいて候補アノテーションの排除を行う。
【００２２】
次に、ペア排除段階２８で、必要な特徴閾値を通過するが、別の非常に近接した活性化と比較すると重要ではない候補アノテーションは捨ててもよい。
【００２３】
最後に、第２のアノテーション検出器段階３０で、各候補アノテーションにその特徴値に基づいたスコアが付与される。スコアの閾値を上回る候補アノテーションのみが、有効なアノテーションであると考えられ、これらのタイミング及び特徴は、候補アノテーションのマップ生成といったプロセッサの更なる操作においてプロセッサによって用いられる。
【００２４】
アルゴリズムの要素については下記でより詳細に説明する。
【００２５】
アルゴリズムの核心部は、３つの基本的な観察に依存する。
・ユニポーラ活動は信号の振幅における急な下方偏向によってマークされる。これらの偏向は、活動速度信号における局所最小値として容易に同定が可能である（即ちユニポーラ信号の導関数）。しかしながら、全ての局所最小値が活動の指標となるわけではなく、いくつかは、ノイズ又は遠距離場活動の結果である。したがって、アルゴリズムの目的は、局所の実際の活性化と関連付けられる速度信号における局所最小値と、そうではないものと、を区別することにある。これは、次の観察を用いることで可能となる。
・遠距離場活性化は、ユニポーラのペアの電位にほぼ等しく影響を与える。したがって、これらのバイポーラ対応物は、遠距離場活性化の間に残留活動しか有さない。これは、電極のうちの一方が活性化部位（activation）に近接し、もう一方が比較的遠い（数ｍｍで十分である）局所活動測定の間には当てはまらない。この場合、バイポーラ信号は、ユニポーラ信号の勾配変化と同時的な勾配変化を更に示す（図２及び下記の説明を参照されたい）。この現象を用いて、近距離場活性化から生じるユニポーラの急な偏向と、遠距離場活性化から生じるものと、を区別することが可能である。
・候補アノテーションの複数の特徴を組み合わせることにより、より強固なアノテーションの検出を実施するための方法が提供される。例えば、局所活性化の結果としてのバイポーラ振幅の変化を、特徴として用いることができる。しかしながら、時には遠距離場活性化はバイポーラコンポーネントを有し得る（例えば、心室弁付近の下部心房（inferior atria）をマッピングする時）。したがって、バイポーラ振幅のみに基づく意思決定は失敗する場合がある。しかしながら、信号の追加的な特徴が用いられると、意思決定はより強固なものになり得る。１つのこうした特徴は、活性化部位周囲のユニポーラ振幅であり得る（図２）。
【００２６】
図２は、本発明の一実施形態による、バイポーラ信号、及びユニポーラ正極信号によって測定される活動の例である。グラフ４０はバイポーラ信号を示し、グラフ４４はユニポーラ信号を示す。区域Ａにおける左側の急な下方偏向は、ユニポーラ及びバイポーラ信号において同時発生的な近距離場活動である。遠距離場の心室活性化中に区域Ｂに示されるように、ユニポーラ信号は変化するが、バイポーラ活動は無視できる程度である。本発明の実施形態は、上記で例示されたものと類似した信号の複数の特徴を用いて、局所活性化と遠距離場活性化とを分離するのを補助する。例えば、区域Ａでは、ユニポーラ振幅及びその速度はバイポーラ信号と類似しているが、一方で、区域Ｂでは、バイポーラ信号と比較してユニポーラ信号振幅ははるかに大きく、その速度ははるかに速い。
【００２７】
以下の説明は、図１に例示されたアルゴリズムの要素を説明する。
【００２８】
前処理段階２２及び特徴抽出段階２４（図１）
前処理段階及び特徴抽出段階の目的は、ユニポーラ及びバイポーラ信号における干渉を除去及び減衰させ、同時に後続の段階で用いられる信号のこれらの特徴を維持して、かつ強調することにある。分かりやすさのために、本明細書で説明される作用は段階２２及び２４で発生すると仮定してあるが、これらの作用のうちの少なくともいくつかは、アルゴリズムのその他の段階で発生する場合があることを理解されたい。我々が保持したいと欲する特徴は、活性化のモルホロジーであり、その理由は、これが勾配の変化を反映するためである。典型的に捨てられる特徴は、勾配の測定を誤らせ得る付加的な信号として、かつ高周波ノイズとしても作用するベースラインのふらつきである。段階２２及び２４は、４つの亜段階に分割される。
【００２９】
１．ユニポーラ前処理亜段階
ユニポーラ前処理段階は、以下の工程を順番に適用することからなる。
１．１ＫＨｚでのベースラインの推定及びサブトラクション（中央値フィルタ及び低域通過フィルタ（ＬＰＦ）を使用した）
２．８ＫＨｚへのアップサンプリング（因数８によるサンプリング及び保持、又はその他のアップサンプリング技法）
３．第１の平滑化フィルタ−ＬＰＦ  ＦＩＲ（−６ｄｂ＠１５５Ｈｚ、１４５タップ）。フィルタは、５００Ｈｚの６５タップ等リプルフィルタとガウス１０ミリ秒窓とのコンボリューションである。
４．第２の平滑化フィルタ。用いられるフィルタは、典型的には５００Ｈｚの低域通過フィルタシステムである、アンチエイリアスＬＰＦである。
５．導関数
【００３０】
工程５の導関数は、第１のアノテーション検出器段階２６（図１）におけるユニポーラアノテーション検出器（フェーズＩ）への入力値として用いられる。工程４の追加的にフィルタリングされた信号の出力値は、アルゴリズムの特徴抽出段階２４のために用いられる。
【００３１】
２．バイポーラ前処理亜段階
バイポーラ前処理段階は、以下の工程を順番に適用することからなる。
１．１ＫＨｚでのベースライン推定及びサブトラクション（中央値フィルタ＋ＬＰＦ）
２．８ＫＨｚへのアップサンプリング（因数８によるサンプリング及び保持）
３．平滑化−ＬＰＦ  ＦＩＲ（−６ｄｂ＠３１０Ｈｚ、１１３タップ）フィルタは、５００Ｈｚの６５タップ等リプルフィルタとガウス６ミリ秒窓とのコンボリューションである。
４．導関数
【００３２】
バイポーラ前処理段階（バイポーラ導関数）の最終出力値は、上記のユニポーラアノテーション検出器（フェーズＩ）への入力値として用いられる（図１）。
【００３３】
３．ベースラインのふらつき推定亜段階
心臓内（ＩＣ）信号は、カテーテルの運動、被験者の運動、及び組織との界面を変化させる呼吸から生じる付加的なベースラインのふらつき信号を含む場合がある（図３及び下記の説明を参照されたい）。これらの運動アーチファクトが含む成分は、ほとんどが低周波成分である。しかし、近距離場活動信号もまた、これらのスペクトル帯に有意なエネルギーを含む場合がある。したがって、高域通過ＩＩＲ又はＦＩＲフィルタによる従来の除去方法には問題があり、ＩＣ信号に歪み及び形態変化を引き起こす場合がある。結果として、我々が使用する選択されたアプローチは、ベースラインのふらつきの推定（図３）、及び信号からそれを減算することに基づく。
【００３４】
図３は、本発明の一実施形態による、ベースラインのふらつきの除去を図示するグラフである。ユニポーラ信号５０は、当初はベースラインのふらつきに寄与する低周波アーチファクトによって汚染されている。ベースラインを推定する目的は、続いて信号から減算されるベースラインを計算することである。図では、計算されたベースライン５４がユニポーラ信号にオーバーレイされている。ベースラインのふらつきはユニポーラ導関数の推定にノイズを加えるおそれがあり、したがってアノテーション検出に影響を与え得るため、ベースラインのふらつきを拒絶することが重要である。
【００３５】
ベースラインのふらつき推定、及びオリジナルからそれを減算することは、図４に図示される通りに一連の２つのフィルタを用いて近距離場活動を除去することにより達成される。
【００３６】
図４は、本発明の一実施形態による、ベースラインのふらつき除去システムのブロック図である。典型的には６０ｍｓの窓を有する中央値フィルタ６０は、原信号からの活動を除去するように設計され、一方で、一実施形態では、約１０Ｈｚの典型的なカットオフを有する８９タップＦＩＲハニングフィルタであるＬＰＦ  ６４は、中央値フィルタから生じた縁部を平滑化するように設計されている。最終的に、否定６８、続いて加算７２の処理によってベースライン推定を原信号から減算して、ベースラインのふらつきを含まない信号を得る。
【００３７】
４．平滑化導関数亜段階
図５は、本発明の一実施形態による、２つのガウスフィルタのグラフである。信号における鋭い偏向ポイントの検出は、信号の速度に基づくものであり、したがって導関数アプローチを用いる。しかしながら、導関数は高域通過フィルタとして作用するため、高周波ノイズを強化する。したがって、導関数推定のノイズを減少させるために我々は平滑化関数を用いる。我々が用いる平滑化関数は、図５に図示されるユニポーラガウス関数８０及びバイポーラガウス関数８４を含む、正規化されたゼロ平均のガウス関数である。これらのユニポーラ及びバイポーラガウスフィルタは、それぞれ±２ｍｓ及び±１ｍｓの時間窓のエネルギーの９０％を有する。したがって、これらの値よりも大きい活性化又は接近する遠距離場は、事実上無視され、導関数値には影響しない。
【００３８】
アノテーション検出器−Ｉ段階２６（図１）
ここで、図６〜９を参照する。図６はアノテーション検出器−Ｉ段階２６の概略的ブロック図であり、図７は、ユニポーラ及びバイポーラ信号、並びにその導関数のグラフであり、図８は、アノテーションアルゴリズムの第１の拒絶フェーズを図示したグラフを有し、図９は、本発明の一実施形態による、局所及び遠距離場の候補アノテーションを図示するグラフを有する。
【００３９】
図６を参照すると、下記の表Ｉは、検出器で用いられるパラメータ、及びブロック図での対応する頭字語を提供する。
【００４０】
【表１】
2018023787.tif 000003
【００４１】
図７は、ユニポーラアノテーション周辺のゼロのバイポーラ勾配の例を示す。グラフは、ユニポーラ遠位信号１００、その導関数１０２、その局所活性化（Ａ）、並びにバイポーラ信号１０４、及びその導関数１０６を示す。ユニポーラ偏向ポイント（Ａ）では、バイポーラ導関数はほぼゼロであり、バイポーラ振幅の大きな変化を示唆しないことに留意されたい。
【００４２】
図８は、図１のアノテーションアルゴリズムの第１の拒絶段階を図示したグラフを有する。上のグラフ１１０はユニポーラ信号を示し、下のグラフ１１４はその平滑化導関数を示す。黒色ドット１１８は閾値未満の導関数信号の最小値であり、可能なアノテーションポイントとして更に検討され、一方で灰色ドット１２２は拒絶される閾値を超える最小値を示す。
【００４３】
図９は、本明細書で説明するバイポーラ導関数とユニポーラ導関数との比の特徴を用いた局所（Ａ）候補アノテーションと遠距離場（Ｂ）候補アノテーションとの分離を図示する。図は、ユニポーラ１３０及びバイポーラ１３２信号、並びにユニポーラ１３６及びバイポーラ１３８導関数を示す。局所活性化では、ユニポーラ導関数の変化は、２ｍｓの活動窓（activity window）１４０によって図示されるように、バイポーラ導関数の変化によって達成される。しかしながら、窓１４４によって図示されるように、これは遠距離場によって導かれた偏向（Ｂ）には当てはまらず、したがって、遠距離場の場合におけるバイポーラ及びユニポーラ勾配の変化の間の比は、必要な比の閾値を下回る。
【００４４】
図６に戻ると、アノテーション検出器−Ｉブロックの入力値は、試験されている適切なユニポーラ信号導関数、その極性、及びその平滑化バイポーラ導関数である。ブロックの出力値は、アノテーション指標及びこれらの勾配値である（アノテーション指標におけるユニポーラ導関数値）。勾配値はアノテーションのスコアとしての役割を果たす。
【００４５】
本発明の一実施形態では、閾値（典型的には−０．０１ｍｖ／ｍｓ）を下回る最小値ポイントを発見することによりブロック９０及び９２にユニポーラ信号の右下がり勾配における偏向ポイントが検出されるが、これは図８を更に参照されたい。活動は、以下の他の２つに加えて、この条件を典型的に満たす。
１．ユニポーラ偏向ポイント（典型的には±２ｍｓ）周囲の時間窓のバイポーラ平滑化導関数信号（Ｓ−ＢＩＰ）の値は、閾値ＴＨ−ＢＩＰを負の様式で超えるべきである。
したがってＳ−ＢＩＰ＜ＴＨ−ＢＩＰである。一実施形態では、ＴＨ−ＢＩＰは典型的には約０．００８ｍｖ／ｍｓである。
２．このバイポーラ平滑化導関数値とユニポーラ平滑化導関数勾配値との間の比はＴｈ−Ｒａｔｉｏよりも高く、典型的には約０．２であるべきである。
【００４６】
＃１及び＃２は、ブロック９４及び９６、並びに判定９８で評価される。
【００４７】
図６を参照すると、バイポーラ導関数値（Ｓ−ｂｉｐ）は、正極と負極とでは異なる方法で計算される。開示された実施形態では、正極については、それは２ミリ秒の時間窓内の最小値であり、負極については、それはその時間窓内の最大値の負値である。アノテーションポイントで導関数ではなく時間窓を用いる理由は、特定の病理及び／又は向き（カテーテルと波の伝播の方向の）において、ユニポーラ活性化間の活動の時間遅延が相殺され得るために、所与のポイントでのバイポーラ信号が小さい又は更にはゼロであり得るからである（図７）。正極における先端の活動はバイポーラ信号において右下がり勾配として認められ、一方、負極での活動はバイポーラ信号において右上がり勾配として認められるため、その値は、正極と負極では異なる方法で計算される。
【００４８】
ユニポーラ導関数とバイポーラの導関数との比は、分類の判定基準としても用いられ得るが、それは、この判定基準が近距離場活動と遠距離場活動とを区別し得るためである。近距離場活動では、右下がり勾配の活動の少なくともいくらかは、典型的にはバイポーラ信号で表され、一方で、遠距離場の場合は、バイポーラ信号は残留活動のみを有し得る。
【００４９】
ペア排除段階２８（図１）
アルゴリズムのペア排除段階は、単一の活動から生じる２つのアノテーションをマージする役割を担う。このスプリット現象は、何らかの理由で、もう一方の電極で記録された活動、又は一方の電極に他方よりも強い影響を与える遠距離場活動のいずれかによって、近距離場活動の右下がり勾配が瞬間的な右上がり勾配を含む場合に発生し得る。瞬間的な右上がり勾配は、信号の導関数に２つの最小値を生じさせ、これらが十分に強い場合は、これは２つのアノテーションを生じさせる。これらのケースを除外するために、我々は右上がり勾配に起因する信号の変化を評価する。
【００５０】
過度に離れていない（典型的には５０ｍｓ未満）同じユニポーラ信号における全てのアノテーションペアを、スプリットに関して分析する。ユニポーラ導関数信号の２つの候補アノテーション間のセグメントを右上がり勾配に関して分析する。右上がり勾配の振幅が有意であると考えられる時は、２つのアノテーションは維持される。そうでない場合は、より小さな右下がり勾配を有するアノテーションは捨てられる。
【００５１】
図１０は、本発明の一実施形態による、候補アノテーションのマージ、及びどのように拒絶判定基準が用いられるかを図示するグラフ１５０である。グラフは、ユニポーラ導関数信号及び２つの可能なアノテーション（Ａ［ｉ］及びＡ［ｉ＋１］）と記された円）を示す。
【００５２】
ペア排除ブロック３８の目的は、最小の導関数振幅と、２つの可能なアノテーション間のピークＰと、の間の右上がり勾配振幅の変化（垂直の両矢印でマークされている）が有意であるか否かを判定することである。変化が有意であると考えられる場合は両方のアノテーションが維持され、そうでない場合は、より弱い活性化（Ａ［ｉ］）が捨てられる。
【００５３】
したがって、捨てられるアノテーションＡ［ｉ］に関しては、５０ｍｓの時間窓Ａ［ｉ＋１］内でより強い勾配を有する、任意の隣接する候補アノテーション間での、ピーク振幅（Ｐ）に対する相対的な変化が考えられる。ピークが有意に高い場合は、このポイントは拒絶されない。数学的条件では、一実施形態では、（Ｐ−Ａ［ｉ］）／（０．０２−Ａ［ｉ］）の値が０．５未満である場合、アノテーションＡ［ｉ］が捨てられる。即ち、５０ｍｓの時間窓内の１つ又は２つ以上のアノテーションが上記のルールに従う場合は、アノテーションＡ［ｉ］は拒絶される。
【００５４】
アノテーション検出器ＩＩ段階３０（図１）
前フェーズを通過した候補アノテーションを、追加的な特徴及び測定基準を用いてこのブロックで再評価する。このブロックを通過して、更にユーザーのバイポーラ電圧で制御された閾値も通過するアノテーションのみが、有効なアノテーションと考えられる。各アノテーションに関して、複数の特徴を計算する。各特徴値に対して、特徴の信頼値に対応する０〜１の範囲のファジースコアを付与する。最後に、全てのスコアを組み合わせてこれらの値をグローバルスコア閾値に対して試験する。グローバルスコア閾値を通過する、即ち高い信頼値を有するアノテーションは有効なアノテーションとみなされ、通過しない、即ち低い信頼値を有するアノテーションは拒絶される。
【００５５】
本明細書で説明するファジー関数は、本発明の一実施形態で用いられる関数の例である。しかしながら、その他のファジー関数又はその他の確率論／関数も、当業者には明らかであり、そのような関数は全て本発明の範囲内に含まれると想定される。更に、特定の要件では、複数のファジースコアを用いてもよい（例えば、強い又は小さいバイポーラ信号をハイライトするファジー関数など）。
【００５６】
全てのファジー関数は、０〜１の間で有界である。
【００５７】
ブロックが用いる特徴は、以下の通りである。
１．ユニポーラ導関数値
２．ユニポーラ勾配の持続時間ｓ２
３．時間窓ｓ２でのユニポーラ勾配の振幅
４．上記の持続時間と振幅との間の比
５．時間窓−ｓ２におけるバイポーラ信号の振幅
【００５８】
５つの特徴について下記で説明する。
【００５９】
１．ユニポーラ導関数
図１１は、本発明の一実施形態による、ユニポーラ導関数のファジー関数のグラフ１６０である。グラフは、導関数に割り当てられるスコアｆ（ｓ１）を提供し、本明細書では導関数値はｓ１と称される。グラフに示すように、−０．０７未満の導関数値には１のスコアが与えられ、−０．０７を超える値は、勾配の−０．０１８で０．５のスコアが到達されるように直線状に減少する。−０．０１未満の導関数値には０のスコアが与えられる。
【００６０】
ユニポーラ導関数ｓ１は両方の検出器段階で用いられるが、０．０１ｍｖ／ｍｓの二分閾値（dichotomy threshold）を有する第１の段階とは異なり、本明細書ではその値はスコアｆ（ｓ１）を提供するために用いられる。スコアが高いほど、この特徴のみに従ってこれが有効なアノテーションである可能性がより高くなる。
【００６１】
２．ユニポーラ活動のセグメンテーション及び持続時間
図１２は、本発明の一実施形態による、ユニポーラ信号のセグメンテーションを図示するグラフを示す。セグメンテーションについて下記でより詳細に説明する。ユニポーラ信号１７０及びその導関数１７４は、候補アノテーションの時間指標１７６（黒点）の周囲に図示されている。閾値を示す水平点線１８０は、両方向の探索セグメント（search segment）（典型的には約±２５ｍｓ）をマークする。一実施形態では、セグメント値は、アノテーションポイントにおける絶対最大ユニポーラ導関数値の２０％として定義される。セグメントＡ、Ｂは、信号導関数が閾値未満である検索窓内の時間間隔をマークする。この実施例の最終セグメントは、セグメントＡか、あるいは特定の条件（本明細書の下記で説明する）が満たされる場合は、開始Ａで始まりＢで終結する結合セグメントか、のいずれであってもよい。
【００６２】
我々がユニポーラ信号から導く特徴は、候補アノテーション周囲の右下がり勾配セグメントの持続時間ｓ２である。この目的は、その初期下降から右上がり勾配を開始するまでの間にユニポーラの右下がり勾配を検出することである。その動機は、そのセグメントにおける持続時間の特性、振幅、及びこれらの関係などの信号の特徴を検査して、これらを分類器の基準として用いることである。本発明者らは、この課題のためにいくつかの方法を検討し、それら全てが明白な単一勾配のケースでは効果があったが、勾配セグメント内で勾配の傾向の変化、及び局所的ピークを有する複雑なケースでも効果があるために、本明細書で説明される方法が選択された。
【００６３】
図１２を参照すると、セグメンテーションは以下の工程を介してユニポーラ導関数を分析することに基づく。
１．候補アノテーションの時間指標１７６上に集中したユニポーラ信号導関数１７４の５０ｍｓセグメントを、それに対してセグメントが画定され得る閾値線１８０の最大スパンとみなす。図１２では、スパンは線１８０の端点１８２と１８６との間にある。我々は、ユニポーラ信号の右下がり勾配セグメントは、この５０ｍｓの時間窓内に制限されているものとみなす。セグメントがこれよりも大きい場合は、これは、この５０ｍｓの区画に強制的に限定される（force-bounded）。
２．導関数セグメントの振幅を、一定の閾値に対して比較する。本明細書で説明する実施形態における閾値は、候補アノテーション時間においてユニポーラ導関数値の２０％であると推定される。その値未満である線１８０上のセグメントは、図１２で２つのセグメントＡ及びＢとしてマークされている。
３．次の工程は、セグメントの境界を計算して、これらのセグメントにおける信号の絶対値の合計に対応する、各サブセグメントの個別の導関数未満領域を合計することである。
４．セグメントのマージ。間を占めるセグメント、及びこれらの面積に基づいて、最終セグメントが主セグメント（Ａ）又は追加的セグメント（Ｂ）のどちらを含むべきか決定する。一実施形態では、セグメントを結合するために、隣接する端点は互いから１ｍｓ以下でなければならず、追加的セグメント（Ｂ）の面積は、Ｂの信号δがＡの信号δの３０％未満となるように、主セグメントの３０％未満であるべきである。
【００６４】
本明細書ではｓ２と称される、上記の工程から求められた持続時間に、続いて、図１３に関連して下記で説明されるファジー関数を用いたスコアｆ（ｓ２）を割り当てる。
【００６５】
図１３は、本発明の一実施形態による、ユニポーラ持続時間のファジー関数のグラフ１９０である。２ｍｓ未満の非常に短い勾配は、実際の活性化から生じたものではないことが疑われ、非常に長い活性化は、おそらく遠距離場のイベントである。更に、局所的な有効活性化のユニポーラ持続時間は、短すぎてもいけず、長すぎてもいけない。上記の観察は、スコアｆ（ｓ２）を提供する図１３のファジー関数に要約される。関数ポイント１９２、１９４は、｛２，０．５｝、｛１９，０．５｝であり、勾配はそれぞれ０．５及び−０．５である。
【００６６】
３．ユニポーラ振幅
図１４は、本発明の一実施形態による、ユニポーラ振幅のファジー関数のグラフ２００である。ユニポーラ振幅は、検出された活動セグメント（ピーク対ピーク）の持続時間ｓ２におけるユニポーラ信号の振幅である（本明細書ではｓ３と称される）。一実施形態では、ファジー関数勾配は、ポイント２０２、２０４｛０．１，０．５｝、｛０．４２，１｝で交差する。ファジー関数から導かれたスコアｆ（ｓ３）は、信号の振幅よりも高い。即ち、高いスコア及び高い振幅の場合、遠距離場信号が大きな振幅を有していない限り、信号は局所活性化に由来している可能性が高い。
【００６７】
４．ユニポーラ持続時間と振幅との比
図１５は、本発明の一実施形態による、ユニポーラ持続時間と振幅との比のファジー関数のグラフ２１０である。活動が長くなるほど振幅は小さくなり、これは誤ったアノテーションである可能性が高いため、ユニポーラ持続時間と振幅との比は、高い比の値を除外する。一実施形態では、ファジー関数線の式は、
ｆ（ｓ４）＝−０．０１８４・ｓ４＋１．２８３                    （１）であり、
式中、ｓ４は、持続時間と振幅との比であり、
ｆ（ｓ４）は比に割り当てられるスコアである。
【００６８】
５．バイポーラ振幅
図１６は、本発明の一実施形態による、バイポーラ振幅のファジー関数のグラフ２２０である。ユニポーラ活動セグメント内のバイポーラ振幅（ピーク対ピーク）であるｓ２もまた、候補アノテーションの可能性をスコア付けするために用いられる。値が高いほど、これが真の活性化である可能性が高まる。
【００６９】
ファジー関数の式は、
ｆ（ｓ５）＝２５・ｓ５，０≦ｓ５≦０．０４；ｆ（ｓ５）＝１，ｓ５＞０．０４    （２）であり、
式中、ｓ５はバイポーラ振幅であり、
ｆ（ｓ５）は振幅に割り当てられるスコアである。
【００７０】
振幅は、ガウス低域通過フィルタ及びアンチエイリアシングフィルタの後で、ベースラインによって拒絶されたバイポーラによって平滑化された信号に基づいて計算される。
【００７１】
６．最終スコア
上述したように、各特徴はスコアを付与され、スコアはグローバルスコアの生成時に共に用いられる。この意図は、特徴は、アノテーションの取り入れ又は除外において相互に支援し得るというものである。一実施形態では、我々が用いたスコア法は、以下に定義される通りである。
【００７２】
【数１】
2018023787.tif 000004
式中、ＧＳはグローバルスコアである。
【００７３】
有効であるとみなされるアノテーションの場合、ＧＳの値は、例えば０．８の特定の閾値を通過するべきである。
【００７４】
上記で例示したものとは異なるが、同等の結果を有するグローバルスコアを本発明の実施形態で用いてもよいことが、当業者には明らかになるであろう。こうしたグローバルスコアは、個々のスコアの加重平均、及び／又は個々のスコアのドット積のほぼ全ての組み合わせを含むことができる。こうしたグローバルスコアはまた、ファジー特徴のサブセットに基づくスコア構成を含み得る。本発明の範囲は、こうしたグローバルスコアの全てを含む。
【００７５】
バイポーラ振幅のフィルタリング
いくつかの実施形態では、アルゴリズムの最終段階は、低いバイポーラ振幅を有する場合に検出されたアノテーションを排除する能力をユーザーに提供するように設計されている。必要な振幅の閾値は、ユーザーによって制御される。バイポーラ振幅のフィルタリングは、処理後段階を超えた各アノテーションのバイポーラ振幅と、閾値とを比較する。閾値を超えるバイポーラ振幅を有するアノテーションのみがシステムに渡される。（ユーザーがこの段階をスキップすることを望む場合は、ユーザーは閾値を０に設定して、この段階のルールを排除してもよい）。
【００７６】
各アノテーションのバイポーラ振幅は、ピーク対ピーク振幅、除去されたベースライン、アノテーション時間を中心とした１４ｍｓ窓の１ＫＨｚバイポーラ信号（最大ユニポーラ速度ポイント）を測定することにより定義される。一実施形態では、バイポーラ振幅閾値のシステムデフォルト値は、３０マイクロボルトに設定されている。
【００７７】
このバイポーラ振幅は固定間隔で測定されるため、このバイポーラ振幅は、ファジー制御されたバイポーラ振幅（上述）とは異なる。ファジー分類器は、ユニポーラ活性化の動的セグメントを用いるため、いくつかの実施形態では、動的セグメントは分類器としてより有意性があり得る。更に、この分類器は、ユーザー制御される二値（dichotomic）閾値として用いられる。
【００７８】
アルゴリズム最終出力値
ファジースコア及びバイポーラのユーザーによって制御されるバイポーラ振幅を通過する全てのアノテーションは、プロセッサによって用いることのできる有効なアノテーションとみなされる。
【００７９】
一実施形態では、各アノテーションは以下の特徴を有するべきである。
１．アノテーション時間指標
２．ユニポーラ及びバイポーラ導関数値
３．ファジースコア
４．ユニポーラによって検出された（unipolar detected）右下がり勾配セグメント持続時間
５．そのセグメント内のユニポーラ振幅
６．そのセグメント内のバイポーラ振幅
７．ユーザーによって制御された値のバイポーラ振幅
更に、以下を含む追跡ファイル（trace files）が提供されてもよい。
１．各特徴に対する特定のファジースコア
２．ユニポーラセグメントの開始及び終了時間指標
【００８０】
図１７は、本発明の一実施形態による、装置４００を用いた侵襲性医療処置の概略図である。手順は、医療専門家４０２によって実施され、また一例として、本明細書の下記の説明における手順は、ヒト患者４０６の心臓４０４からＥＣＧ信号を取得することを含むと想定される。
【００８１】
信号を取得するために、専門家４０２は、患者の内腔内に予め配置されたシース４１０にプローブ４０８を挿入する。シース４１０は、プローブの遠位端４１２が、シースの遠位端４１４から抜け出た後に患者の心臓に侵入して、心臓の組織に接触し得るように配置される。
【００８２】
プローブ４０８は、患者の心臓に挿入することが可能で、かつ典型的には磁気追跡システム及び／又はインピーダンス測定システムを用いて追跡可能な任意のタイプのカテーテルを含んでもよい。例えば、プローブ４０８は、ラッソーカテーテル、シャフト様カテーテル、若しくはＢｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒ（Ｄｉａｍｏｎｄ  Ｂａｒ，ＣＡ）製ｐｅｎｔａＲａｙカテーテル、又はこれらのカテーテルに概ね類似するカテーテルを含んでもよい。Ｂｉｏｓｅｎｓｅ  Ｗｅｂｓｔｅｒは、本発明の実施形態で使用可能な磁気追跡システム及びインピーダンス測定システムも製造している。
【００８３】
プローブ４０８は、本明細書で説明されるアルゴリズムの実施中にプロセッサ２０によって使用されるＥＣＧ信号を取得するために用いられる少なくとも２つの電極４１１を含む。
【００８４】
装置４００はプロセッサ２０（図１）によって制御され、プロセッサは、典型的にはフィールド・プログラマブル・ゲート・アレイ（ＦＰＧＡ）として構成されるリアルタイムノイズリダクション回路４２０、続いてアナログ・デジタル（Ａ／Ｄ）信号変換集積回路４２４を含み得る。プロセッサは、信号をＡ／Ｄ回路４２４から別のプロセッサに渡すことができ、また、本明細書に開示されるアルゴリズムを実施するようにプログラム可能である。
【００８５】
プロセッサ２０は、装置の操作コンソール４３０に設置されている。コンソール４３０は、専門家４０２によって用いられてプロセッサと通信するための制御装置４３２を含む。手順の間、プロセッサ２０は、ＥＣＧ信号を取得し、並びに本明細書で開示されるアルゴリズムを実施するために、モジュールバンク４４０内のＥＣＧモジュール４３６と通信する。
【００８６】
ＥＣＧモジュール４３６は電極４１１からＥＣＧ信号を受信する。一実施形態では、信号は、モジュール４３６で低ノイズ前置増幅器４３８を通り、帯域通過フィルタ４４０を介して、主増幅器４４２に転送される。モジュール４３６は、本明細書で説明されたアルゴリズムをプロセッサによって実装するために、ＥＣＧ信号のデジタル化された値をプロセッサ２０に転送するアナログ・デジタル変換器（ＡＤＣ）４４４を更に含む。典型的には、プロセッサ２０は前置増幅器４３８、フィルタ４４０、増幅器４４２、及びＡＤＣ  ４４４の操作を制御する。
【００８７】
かくして、ＥＣＧモジュール４３６によって、電極４１１が受信した、本明細書で参照されたＥＣＧ信号を含むＥＰ（電気生理学的）信号を、プロセッサ２０が取得して分析することが可能となる。信号は、典型的には、表示画面４５０にリアルタイムで更新される電圧／時間グラフとして専門家４０２に示される。
【００８８】
プロセッサ２０及びモジュールバンク４４０のソフトウェアは、例えばネットワークを介した電子的形態でプロセッサにダウンロード可能である。あるいは又は更には、ソフトウェアは、例えば、光学的、磁気的、又は電子的記憶媒体のような非一時的有形媒体上で提供され得る。
【００８９】
装置１２を操作するために、モジュールバンク５０は、典型的には、プロセッサがプローブ４０８の遠位端を追跡することを可能にする１つ又は２つ以上の追跡モジュールなどの、上記のＥＣＧモジュール以外のモジュールを含む。分かりやすさのために、こうしたその他のモジュールは図１に図示されていない。全てのモジュールはハードウェア及びソフトウェア要素を含み得る。
【００９０】
電極４１１によって取得されたＥＣＧ信号を示す表示画面４５０に加えて、更に、本明細書で説明するアルゴリズムの結果４５２を表示画面上でアルゴリズムのユーザーに示してもよい。
【００９１】
上に述べた実施形態は例として挙げたものであり、本発明は上記に具体的に示し、説明したものに限定されない点は理解されよう。むしろ、本発明の範囲は、上記されている種々の特徴の組み合わせ及び部分的組み合わせと、前述の説明を読むことに基づいて当業者が想起するであろう、先行技術に開示されていない変形例及び修飾と、の両方を含む。
【００９２】
〔実施の態様〕
（１）  ヒト被験体の心筋に近接する一対の電極からバイポーラ信号を受信すること、
前記一対の電極のうちの選択された一方からユニポーラ信号を受信すること、
前記ユニポーラ信号の局所ユニポーラ最小導関数、及び前記ユニポーラ最小導関数の発生時刻を計算すること、
前記バイポーラ信号のバイポーラ導関数を計算すること、
前記バイポーラ導関数と前記局所ユニポーラ最小導関数との比を求めること、並びに
前記比が予め設定された閾値の比の値を超える場合は、前記心筋の活性化時間として前記発生時刻を特定すること、を含む、方法。
（２）  前記バイポーラ導関数が予め設定されたバイポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、実施態様１に記載の方法。
（３）  前記局所ユニポーラ最小導関数が予め設定されたユニポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、実施態様１に記載の方法。
（４）  前記発生時刻に対応する信頼水準が予め設定された信頼水準を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、実施態様１に記載の方法。
（５）  前記バイポーラ信号の振幅が予め設定されたバイポーラ信号の閾値を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、実施態様４に記載の方法。
【００９３】
（６）  ヒト被験体の心筋に近接して配置されるように構成された一対の電極と、
プロセッサであって、
前記一対の電極からバイポーラ信号を受信し、
前記一対の電極のうちの選択された一方からユニポーラ信号を受信し、
前記ユニポーラ信号の局所ユニポーラ最小導関数、及び前記ユニポーラ最小導関数の発生時刻を計算し、
前記バイポーラ信号のバイポーラ導関数を計算し、
前記バイポーラ導関数と前記局所ユニポーラ最小導関数との比を求め、かつ
前記比が予め設定された閾値の比の値を超える場合は、前記心筋の活性化時間として前記発生時刻を特定するように構成されたプロセッサと、を含む、装置。
（７）  前記バイポーラ導関数が予め設定されたバイポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、実施態様６に記載の装置。
（８）  前記局所ユニポーラ最小導関数が予め設定されたユニポーラ導関数の閾値を下回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることを更に含む、実施態様６に記載の装置。
（９）  前記発生時刻に対応する信頼水準が予め設定された信頼水準を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、実施態様６に記載の装置。
（１０）  前記バイポーラ信号の振幅が予め設定されたバイポーラ信号の閾値を上回る場合は、前記心筋の前記活性化時間として前記発生時刻を割り当てることのみを含む、実施態様９に記載の装置。
【図１】
2018023787.tif 000005
【図２】
2018023787.tif 000006
【図３】
2018023787.tif 000007
【図４】
2018023787.tif 000008
【図５】
2018023787.tif 000009
【図６】
2018023787.tif 000010
【図７】
2018023787.tif 000011
【図８】
2018023787.tif 000012
【図９】
2018023787.tif 000013
【図１０】
2018023787.tif 000014
【図１１】
2018023787.tif 000015
【図１２】
2018023787.tif 000016
【図１３】
2018023787.tif 000017
【図１４】
2018023787.tif 000018
【図１５】
2018023787.tif 000019
【図１６】
2018023787.tif 000020
【図１７】
2018023787.tif 000021
【外国語明細書】
2018023787000001.pdf