(57)【要約】要約書なし【特許請求の範囲】
【請求項１】
血管内の管腔を検出する方法であって、
血管内へ音響エネルギーを伝送することと、
前記血管内の前記音響エネルギーの伝搬の関数として前記音響エネルギーの反射の信号レベルを評価することと、
前記反射内の低下信号レベルの少なくとも１つの領域を識別することと
を含む方法。
【請求項２】
請求項１記載の方法において、前記評価する行為は前記反射を復調することを更に含む方法。
【請求項３】
請求項２記載の方法において、前記復調する行為は
（ｉ）  前記反射の第１周波数成分を導出することと、
（ｉｉ）  前記反射の第２周波数成分を導出することと、
（ｉｉｉ）  前記第１周波数成分と前記第２周波数成分との部分集合のみを保持するために前記第１周波数成分と前記第２周波数成分とをフィルタすることとを含む方法。
【請求項４】
請求項３記載の方法において、第１周波数を前記導出する行為は前記伝送された音響エネルギーの中心周波数に関係した正弦関数を前記反射に乗算することを含み、かつ第２周波数を前記導出する行為は前記伝送された音響エネルギーの中心周波数に関係した余弦関数を前記反射に乗算することを含む方法。
【請求項５】
請求項４記載の方法において、前記第１周波数成分と前記第２周波数成分とを前記フィルタする行為は差（ベースバンド）周波数成分のみを保持するために前記成分を低域通過フィルタすることを含む方法。
【請求項６】
請求項３記載の方法であって、
平方パラメータを生じるために前記フィルタされた成分の部分集合を平方することと、
平方デシメート化パラメータを生じるために前記平方パラメータを第１サンプリング・レートへデシメートすることと
を更に含む方法。
【請求項７】
請求項６記載の方法であって、前記平方デシメート化パラメータに対する少なくとも１つの深さ依存関数を更に含み、前記少なくとも１つの深さ依存関数は前記血管内の前記音響エネルギーの伝搬損失を少なくとも一部分保証する方法。
【請求項８】
請求項２記載の方法であって、
複数の伝搬間隔中のパワーを決定することと、
前記複数の伝搬間隔の少なくとも１つ内の少なくとも１つの最小を識別することと
を更に含む方法。
【請求項９】
請求項８記載の方法であって、前記反射内の異なった時点に決定された前記少なくとも１つの最小の平均値を決定することを更に含む方法。
【請求項１０】
血管内の少なくとも１つの壁を検出する方法であって、
血管内へ音響エネルギーを伝送することと、
前記血管内の管腔と関連した少なくとも１つの領域を検出することと、
前記管腔に対する血管の前記少なくとも１つの壁の位置を検出することと
を含み、
位置を検出する行為は前記伝送する行為から導出されたＡモード・データを分析することを含む
方法。
【請求項１１】
請求項１０記載の方法において、位置を検出する行為は、
（ｉ）  前記管腔の前記領域と関連したＡモード信号レベルを決定することと、
（ｉｉ）  Ａモード信号レベルが前記管腔の前記領域と関連した前記レベルとより高い信号レベルとの間で変動する少なくとも１つの領域を検出することと、
（ｉｉｉ）  前記少なくとも１つの壁を変動の前記少なくとも１つの領域と関連させることと
を更に含む方法。
【請求項１２】
請求項１１記載の方法において、前記管腔と関連した少なくとも１つの領域を検出する行為は前記Ａモード信号の少なくとも部分を復調することを含む方法。
【請求項１３】
請求項１２記載の方法において、前記復調する行為は
前記Ａモード信号の第１周波数成分を導出することと、
前記Ａモード信号の第２周波数成分を導出することと、
前記Ａモード信号を表すベースバンド成分を得るために前記第１周波数成分と前記第２周波数成分とをフィルタすることと
を含む方法。
【請求項１４】
血管の直径を決定する方法であって、
血管内へ音響エネルギーを伝送することと、
前記血管内の管腔と関連した領域を検出することと、
前記管腔に対する血管の第１壁の位置を検出することと、
前記管腔に対する血管の第２壁の位置を検出することと、
前記第１壁の前記位置と前記第２壁の前記位置とに少なくとも一部分基づいて前記血管の少なくとも部分の直径を決定することと
を含む方法。
【請求項１５】
請求項１４記載の方法において、位置を前記決定する行為は、各々、
（ｉ）  前記管腔の前記領域と関連したＡモード信号レベルを決定することと、
（ｉｉ）  Ａモード信号レベルが前記管腔の前記領域と関連した前記レベルとより高い信号レベルとの間で変動する少なくとも１つの領域を検出することと、
（ｉｉｉ）  前記第１壁又は前記第２壁の位置を変動の前記少なくとも１つの領域とに関連させることと
を更に含む方法。
【請求項１６】
請求項１５記載の方法において、前記第１壁は前記血管の前壁を含み、かつＡモード信号レベルが変動する少なくとも１つの領域を検出する行為は前記管腔の深さよりも浅い深さでこのような少なくとも１つの領域を検出することを含む方法。
【請求項１７】
生被検体の循環系内の血行力学性質を査定する方法であって、
前記生被検体の血管から第１パラメータを測定することと、
前記血管を少なくとも部分的に圧縮することと、
前記圧縮する行為の間中前記血管から第２パラメータを測定することであって、前記第２パラメータを測定する行為は
（ｉ）  前記血管内へ音響波を伝送することと、
（ｉｉ）  前記音響波の少なくとも１つの反響を受信することと、
（ｉｉｉ）  前記血管の直径の推定値を導出するために前記少なくとも１つの反響を分析することと
を含む前記第２パラメータを前記測定する行為と、
前記第２パラメータに少なくとも一部分基づいて校正関数を導出することと、
前記校正関数を使用して第１パラメータを校正することと
を含む方法。
【請求項１８】
請求項１７記載の方法において、校正関数を導出する行為は
前記第２パラメータ内の少なくとも１つのアーチファクトを識別することと、
前記少なくとも１つのアーチファクトに基づいて校正関数を発生することと
を含む方法。
【請求項１９】
生被検体の血管内の血行力学性質を測定する装置であって、
前記血管と関連した少なくとも第１血行力学パラメータを測定するために適合した第１トランスジューサと、
前記血管内へ超音波エネルギーを伝送しかつ前記血管から反響を受信するために適合した第２トランスジューサと、
前記第２トランスジューサに動作上接続されかつ前記反響に基づいて前記血管内の管腔を検出し、かつ前記管腔に基づいて前記血管の直径を推定するように構成された信号プロセッサであって、前記血管の直径の前記推定値は前記第１血行力学パラメータの前記測定を校正するために少なくとも一部分使用される前記信号プロセッサと
を含む装置。
【請求項２０】
請求項１９記載の装置であって、前記第２血行力学パラメータを測定する間に前記血管を圧縮するために適合した圧平デバイスを更に含む装置。
【請求項２１】
請求項２０記載の装置において、前記血管は人間の橈骨動脈を含む装置。
【請求項２２】
請求項２１記載の装置において、前記第１トランスジューサは前記橈骨動脈に近接して配置された圧力トランスジューサを含み、かつ前記第２トランスジューサは前記橈骨動脈にまた近接してまた配置された超音波トランスジューサを含む装置。
【請求項２３】
請求項２２記載の装置において、前記プロセッサは前記反響から導出されたＡモード信号を分析するために適合している装置。
【請求項２４】
生被検体に治療を施す方法であって、
圧力データを測定するために有用である前記被検体の血管を選択することと、
非侵襲的に前記選択された血管から圧力データを測定することと、
前記血管内の管腔の相対的位置と前記血管の少なくとも１つの壁の相対的位置とを検出することによって前記血管から血行力学パラメータを測定することと、
前記血行力学パラメータの前記測定から得られたデータから校正関数を発生することと、
前記血管内の血圧の表示を生じるために前記血圧データに前記校正関数を適用することと、
前記表示に基づいて前記被検体に治療を施すことと
を含む方法。
【請求項２５】
請求項２４記載の方法において、前記血管は人間の橈骨動脈であり、かつ圧力データを測定する行為は前記橈骨動脈に近接した皮膚に当てられたセンサを使用して圧力波形を得ることを含む方法。
【請求項２６】
請求項２５記載の方法において、血行力学パラメータを測定する行為は
音響波を発生することと、
前記橈骨動脈内へ前記音響波を伝送することと、
前記音響波の反響を受信することと、
前記管腔の位置を検出するために前記反響を分析することと、
前記管腔を前記検出する行為に基づいて少なくとも１つの動脈壁の位置を検出することと、
前記管腔の位置と前記少なくとも１つの動脈壁の位置とに基づいて前記動脈の直径を決定することと
を含む方法。
【請求項２７】
生被検体の循環系内の第１パラメータの実値を決定する方法であって、
非侵襲的に循環系と関連した第１パラメータの値を測定することであって、第１パラメータの測定された値は前記第１パラメータの実値とは異なっている第１パラメータを前記測定することと、
循環系の血行力学性質を変更するために前記循環系にストレスを印加することであって、ストレスの前記印加は前記循環系の少なくとも１つの血管内の管腔の位置を検出することによって少なくとも一部分決定されるストレスを前記印加することと、
ストレスが印加される間に循環系と関連した第２パラメータを測定することと、
循環系と関連した第３パラメータを測定することと、
測定された第２パラメータに少なくとも一部分基づいて校正関数を発生することと、
校正関数を使用して測定された第１パラメータを補正することと
を含む方法。
【請求項２８】
請求項２７記載の方法において、
（ｉ）  前記第１パラメータを測定する行為は前記循環系の少なくとも部分内の圧力を測定することを含み、
（ｉｉ）  前記第２パラメータを測定する行為は前記循環系の前記少なくとも部分の直径を測定することを含み、かつ
（ｉｉｉ）  前記第３パラメータを測定する行為は前記管腔内を流れる血液の速度又は運動エネルギーを測定することを含む
方法。
【請求項２９】
血行力学パラメータを拡大縮小する方法であって、
基準波形を用意することと、
前記血行力学パラメータを含む未補正の波形を発生することと、
前記基準波形に少なくとも一部分基づいて前記未補正の波形に対する少なくとも１つの補正を決定することと、
前記未補正の波形を拡大縮小するために前記未補正の波形に前記補正を適用することと
を含む方法。
【請求項３０】
請求項２９記載の方法であって、前記未補正の波形と前記拡大縮小された波形とに少なくとも一部分基づいて伝達関数を計算することを更に含む方法。
【発明の詳細な説明】
【０００１】
本願は、２０００年３月２３日に出願された「生被検体の循環系内の血行力学性質を査定する方法と装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐａｒａｔｕｓ  Ｆｏｒ  Ａｓｓｅｓｓｉｎｇ  Ｈｅｍｏｄｙｎａｍｉｃ  Ｐｒｏｐｅｒｔｉｅｓ  ｗｉｔｈｉｎ  ｔｈｅ  Ｃｉｒｃｕｌａｔｏｒｙ  Ｓｙｓｔｅｍ  ｏｆ  ａ  Ｌｉｖｉｎｇ  Ｓｕｂｊｅｃｔ）」と題するかつ本願の讓受人に讓受され、かつ言及することによってその内容をその全体にわたって本明細書に組み入れてある米国特許出願第０９／５３４，９００号の一部継続出願である。
【０００２】
（関連出願）
本願は、２００１年３月２２日に出願され、本願の讓受人に讓受された「血管位置を含む血行力学パラメータの非侵襲的査定方法と装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐａｒａｔｕｓ  Ｆｏｒ  Ｔｈｅ  Ｎｏｎｉｎｖａｓｉｖｅ  Ａｓｓｅｓｓｍｅｎｔ  Ｏｆ  Ｈｅｍｏｄｙｎａｍｉｃ  Ｐａｒａｍｅｔｅｒｓ  Ｉｎｃｌｕｄｉｎｇ  Ｂｌｏｏｄ  Ｖｅｓｓｅｌ  Ｌｏｃａｔｉｏｎ）」と題する米国特許出願第０９／８１５，９８２号に関する。
【０００３】
（発明の背景）
（１．発明の分野）
本願は、一般に、生被検体の循環系と関連したパラメータを監視する方法と装置、特に、動脈血圧の被侵襲的監視に関する。
【０００４】
（２．関連技術の説明）
動脈血圧測定
いくつかの周知の技術が被検体の動脈血圧波形を非侵襲的に監視するためにこれまでに使用されてきている。すなわち、聴診法（ａｕｓｃｕｌｔａｔｉｏｎ）、振動測定法（ｏｓｃｉｌｌｏｍｅｔｒｙ）、圧力測定法（ｔｏｎｏｍｅｔｒｙ）これである。聴診法技術と振動測定法技術の両方は、被検体の上腕動脈を閉塞（ｏｃｃｌｕｄｅ）する標準膨張性アーム・カフ（ａｒｍ  ｃｕｆｆ）を使用する。聴診法技術は、カフをゆっくりと収縮させるに従って起こる或る種のコロトコフ（Ｋｏｒｏｔｋｏｆｆ）音を監視することによって被検体の心収縮期圧（ｓｙｓｔｏｌｉｃ  ｐｒｅｓｓｕｒｅ）と心拡張期圧（ｄｉａｓｔｏｌｉｃ  ｐｒｅｓｓｕｒｅ）を決定する。他方、振動測定法技術は、カフを収縮させるに従ってカフ内に起こる実圧変化を測定することによって、これらの圧力ばかりでなく被検体の平均圧を測定する。両技術は圧力を決定する、しかし一方でそれは被検体の実血圧波形を再現することができない。それゆえ、真の連続的、拍動間血圧（ｂｅａｔ−ｔｏ−ｂｅａｔ  ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）監視は、これらの技術を使用して達成することはできない。
【０００５】
上に短く説明した種類の閉塞性カフ器械は、一般に、被検体の血圧の長期傾向を検知（ｓｅｎｓｉｎｇ）するにはいくらか有効であってきた。しかしながら、このような器械は短期血圧変動を検知するには一般に有効でなかった。短期血圧振動を検知することは、外科を含む多くの医療応用に決定的に重要である。
【０００６】
動脈圧力測定法の技術は、医療技術上また周知である。動脈圧力測定法の理論に従えば、橈骨動脈（ｒａｄｉａｌ  ａｒｔｅｒｙ）のような、充分な骨性支持器を伴う浅動脈（ｓｕｐｅｒｆｉｃｉａｌ  ａｒｔｅｒｙ）内の圧力は、経壁圧（ｔｒａｎｓｍｕｒａｌ  ｐｒｅｓｓｕｒｅ）がゼロに等しいとき、圧平（ａｐｐｌａｎａｔｉｏｎ）掃引の間中に正確に記録される。用語、「平圧」は、動脈に印加される圧力を変動させるプロセスを称する。圧平掃引は、動脈への圧力を過圧縮から不足圧縮へ又はこの逆に変動させる間の時間間隔を称する。漸減（ｄｅｃｒｅａｓｉｎｇ）圧平掃引の開始では、動脈は「ドッグ・ボーン（ｄｏｇ  ｂｏｎｅ）」形状に過圧縮されるので、圧力パルスは記録されない。掃引の終で、動脈は不足圧縮されるので、最小振幅圧力パルスが記録される。掃引内では、動脈壁張力が圧力計表面に平行である間中に圧平が起こると想定する。ここで、動脈圧はその表面に垂直でありかつ圧力測定計センサによって検出された唯一のストレスである。この圧力では、得られた最大ピーク間振幅（「最大拍動（ｍａｘｉｍｕｍ  ｐｕｌｓａｔｉｌｅ）」圧はゼロ経壁圧に相当すると想定する。この理論は、図１にグラフで例示してある。図１で、骨又は他の堅牢な部材が動脈の下にあると想定する。
【０００７】
圧力測定法技術を実施する１先行技術のデバイスは、末梢動脈（ｐｅｒｉｐｈｅｒａｌ  ａｒｔｅｒｙ）、例えば、橈骨動脈の上に横たわる組織に当てるミニアチャ圧力トランスジューサの堅牢なアレーを含む。それらのトランスジューサは、各々がその下に横たわる被検体の組織内の機械力を直接検知し、かつ各々がその下に横たわる動脈の小部分のみをカバーするように寸法を与えられている。アレーは、組織に対して圧迫されて、その下に横たわる動脈を圧平し、それによって動脈内の拍動間圧力変動を組織を通してトランスジューサの少なくとも或るものに結合させる。被検体上のアレーの位置にかかわらず、少なくとも１つのトランスジューサが動脈をおおって常にあることを保証するためにいろいろと異なったトランスジューサのアレーが使用される。しかしながら、圧力計のこの型式は、いくつかの欠点を抱えている。第１に、離散トランスジューサのアレーは、一般に、検知される動脈の上に横たわる被検体の組織の連続輪郭と自動的には両立しない。これは、歴史的に、結果のトランスジューサ信号に不正確さを生じることに立ち至る。更に、或る場合には、この非両立性は、組織外傷及び神経破壊を引き起こすおそれがあり、かつ末梢組織への血流量（ｂｌｏｏｄ  ｆｌｏｗ）を制限するおそれがある。
【０００８】
他の先行技術は、動脈の上に横に単一圧力測定法センサをより正確に置き、それによってセンサを動脈内の圧力変動により完全に結合することを求めている。しかしながら、このようなシステムは、センサが幾何学的には「心出しされている（ｃｅｎｔｅｒｅｄ）がしかし信号結合のためには最適には位置決めされない位置にセンサを置くことがあり、かつ測定の間中の被検体の動きに因って、典型的に、比較的頻繁に再校正又は再位置決めを必要とする。
【０００９】
圧力測定法システムは、監視される被検体上の圧力トランスジューサの配向にまた普通極めて敏感である。特に、このようなシステムは、トランスジューサと動脈との間の角関係が「最適」入射角から変動するとき正確性の劣化を示す。これは重要な考慮すべき問題である。というのは、２つの測定がデバイスを動脈に対しておよそ精確に同じ角に置く又は維持するためしはないからである。横センサ位置決めへの前述のアプローチの多くは、手頸表面の湾曲のような、被検体の解剖学的特徴を償うのに適合していない位置決め機構が多くの場合原因で、横位置にかかわらず動脈と一定角関係を維持することができないという同じような難点がある。
【００１０】
動脈圧力測定システムの他の重大な欠点は、一般に、動脈壁圧縮のレベルを連続的に監視しかつゼロ経壁圧という最適レベルへ調節することができないというそれらの不能力性である。一般に、動脈壁圧縮の最適化は、周期的再校正によってのみ達成されている。これは被検体監視機能の中断を必要としており、中断は決定的に重要な期間中に起こるおそれがある。この不能力性は、臨床環境での圧力計の受入れを厳しく制限する。
【００１１】
圧力測定法アプローチの更に制限は、血管の内部から血管をおおう皮膚の表面上の測定点への不完全圧力パルス転送に関する。特に、動脈圧縮の最適レベルが達成されるときでさえも、血管壁を通して及び組織を通しての動脈血圧の皮膚表面への不完全かつ複雑結合があるので、血管内で起こる圧力変動の大きさが皮膚上に置かれた圧力測定法のセンサ（圧力トランスジューサ）によって測定されたものと異なっているようなことになる。それゆえ、皮膚で測定されたいずれの圧力信号又は圧力波形も動脈内の真圧と必然的に異なる。手頸の動脈壁、組織、筋系（ｍｕｓｃｕｌａｔｕｒｅ）、腱、骨、皮膚の肉体的応答をモデルすることは、一筋縄で行くものではなくて、各個人毎に本来的に不確定性と異常性を含む。これらの不確定性と異常性は、予期できない誤りを圧力測定法センサを経由して行われた血圧のいずれの測定にも導入する。
【００１２】
圧力測定法による圧力測定を校正する１先行技術の方法は、振動測定法デバイス（すなわち、圧力カフ又はこれに類似のもの）を利用して「実」圧情報を周期的に得、この情報は次いで圧力測定法測定を校正するために使用される。このアプローチは、デバイス校正を維持するために、進行中校正事象を、特にカフの膨張／収縮を遂行する必要があるという難点がある。このような校正事象は、注意を散らせ、厄介であり、かつ比較的長い周期性で以てしか実際に遂行することができない。更になお、この技術は、循環系内に起こる実血行力学変化の測定に基づいて校正するのではなく、むしろ実変化の表現であることもないこともある外部測定に基づいている。介在した組織等に因る血管からセンサへの不完全パルス転送に対する補正用機構は、いずれにしても、備わっていない。
【００１３】
他の先行技術による校正技術は、測定された血行力学パラメータに対するオフセット又は補正を発生するために、血管内を流れる血液に擾乱を伝送又は誘発し、かつその後、測定されたパラメータ（例えば、血圧波形）中のその信号の成分を検知することを企図している。例えば、ヴァイタル・インサイト社（Ｖｉｔａｌ  Ｉｎｓｉｔｅ，Ｉｎｃ．）に讓受された「血圧を決定するために誘発擾乱を測定する装置と方法（Ａｐｐａｒａｔｕｒｓ  ａｎｄ  Ｍｅｔｈｏｄ  ｆｏｒ  Ｍｅａｓｕｒｉｎｇ  ａｎ  Ｉｎｄｕｃｅｄ  Ｐｅｒｔｕｒｂａｔｉｏｎ  ｔｏＤｅｔｅｒｍｉｎｅ  Ｂｌｏｏｄ  Ｐｒｅｓｓｕｒｅ）」と題する米国特許第５，５９０，６４９号（’６４９特許）を見られたい。’６４９特許のアプローチの下で、血圧の変化から表向き生じる種々の血行力学パラメータの変化がモデルされかつデバイス内に記憶され、かつ圧力測定法センサから得られたデータと比較される。しかしながら、このアプローチは、校正オフセットが評価される被検体の血行力学パラメータの直接測定によっては決定されないで、血圧と擾乱波速度との間の関係をモデルすることによって決定される。すなわち、速度と位相が血圧の変化と或る種の関係を持つようにモデルされる。したがって、理論上は、擾乱波の速度／位相の観察された変化が評価される被検体内の実血圧の推定を発生するために使用することができる。このシステムの限界は、多くの複雑な、非線形、相互依存パラメータを正確にモデルするばかりでなく、これらの多くのパラメータの時間分散を予測する能力によって明らかに支配される。
【００１４】
循環系の血行力学と疾患
血行力学の科学、又は体内を流れる液体（血）流量の分析は、循環系の疾患又は同系内の欠陥を検出し及び／又は診断するために有効に現在使用される。例えば、弁膜疾患、心臓の構造的欠陥、静脈疾患、低下心機能、及び動脈疾患は、いかに血液が循環系の種々の部分を通して流れるか調べることによって査定される。特に注目されるのは、狭窄症（ｓｔｅｎｏｓｉｓ）（すなわち、動脈斑（ａｒｔｒｉａｌ  ｐｌａｑｕｅ）等に因る有効断面積の遮断又は縮小）のような動脈疾患の分析である。生被検体の血管内の狭窄症の程度が変動するに連れて、循環系のパラメータに及び被検体の全面的健康に或る種の変化が起こる。図２に例示したように、仮想血管内狭窄症の変動する程度は、その血管を全体的に比例する程度に閉塞することになる。すなわち、狭窄症がないならば閉塞と付随症状を生じない一方、完全な狭窄症は完全な閉塞を生じ、血管を通して血液は流れず、かつ非常に差し迫った症状が被検体に現れる。それらの間のどこかに入る狭窄症のレベルでは、反応（ｒｅｓｐｏｎｓｅ）はいくぶんもっと複雑である。例えば、被検体は、所与の血管の有効断面積を非常に顕著に縮小させる狭窄症を患うが、なお運動（ｅｘｅｒｃｉｓｅ）の普通のレベルの下ではなんらかの症状があっても非常に僅かしか現れないことがある。しかしながら、同じ被検体が運動を増すと共に劇的症状を示すことがあり得る。患者が更にそれに励むと、波及を受けた組織はより高い代謝要求を生じてこれが灌流（ｐｅｒｆｕｓｉｏｎ）の増大を必要とする。普通、血管拡張（ｖａｓｏｄｉｌａｔｉｏｎ）と側副枝化（ｃｏｌｌａｔｅｒａｌｉｚｅｄ）血流量は、より高い容積要求を満たすために容積流量を増す補償機構を用意する。しかしながら、血管は顕著に狭窄しているから、補償機構は普通の無運動要求を満たすために既に利用されている。結果として、身体は容積要求を増すことができない。というのは、血管の狭窄（縮小）面積の抵抗を克服することに関連したエネルギー損失を最小限にする術（ｗａｙ）を持たないからである。容積流量が増えないならば、増大した代謝要求は満たされず、抹消組織は虚血症（ｉｓｃｈｅｍｉｃ）になる。
【００１５】
内圧（Ｐ）、及び血質量流量（Ｑ）又は血液速度（ｖ）を有する流体系として狭窄症動脈をモデルすることによって、次のように動脈内の血液の流量を説明するために周知のベルヌーイの式の修正バージョンを応用することができる。
【００１６】
【数１】
2004521660.tif 000002
【００１７】
それゆえ、前掲の関係は、他が既知のとき、１つの血行力学パラメータを査定するために使用されることがある。例えば、動脈内の狭窄症を横断しての圧力勾配（ΔＰ）は、狭窄症を通して流れる血液の速度についてデータを得、次いでこの速度データを式（１）内に使用することによって推定される。速度データは、スペクトル・ドップラー超音波（ｓｐｅｃｔｒａｌ  Ｄｏｐｐｌｅｒ  ｕｌｔｒａｓｏｕｎｄ）のような、多くの周知の技術によっても得ることができる。
【００１８】
しかしながら、動脈内に存在する狭窄症及び他のこのような疾患の酷しさを査定するに当たってのそれらの有用性にかかわらず、先行技術の血行力学評価技術は、いずれかの所与の時刻に動脈内の絶対血圧を査定することは有効にはできない。理論上、循環系の反応の正確モデルは、（真の動脈圧のような）系内のパラメータの値を他のパラメータの既知の値又は測定された値に基づいて推定するために使用することができる。しかしながら、承知のように、生器官の、特に人間の循環系は、極端に複雑で、誇張なしに数千の相互接続血管を備える。この系は、なかでも、夥しい毛細管、静脈、動脈であって各々がそれら自体の独特の肉体的性質を有するものを含む。更になお、血管の前述のカテゴリの各々内で、個々の要素（ｃｏｎｓｔｉｔｕｅｎｔ）は、著しく異なった性質を有しかつ循環系内で反応する。例えば、人間の身体内の２つの動脈は、（ｉ）それらの長さに沿って異なる点で異なる直径を有し、（ｉｉ）他よりも静脈と毛細管に多く又は少なく供給し、（ｉｉｉ）多く又は少なく弾性を有し、及び（ｉｖ）それらと関連した狭窄症を多く又は少なく有する。
【００１９】
血管の各々の性質と反応はまた、身体内への麻酔薬の導入のような、種々の内部刺激及び／又は外部刺激によって異なって影響されることがある。呼吸のような身体内の普通の自律反応でさえも循環系内の圧力に影響し、したがって、考慮する必要があるといえる。
【００２０】
これらの制限を考慮すると、血圧推定に使用するために人間の循環系をその流体力学の面から正確にモデルすることは、不可能ではないにしても、極端に困難になる。たとえ仮説循環系を正確にモデルすることができても、このようなモデルの応用は、各被検体の特定肉体的性質と反応に因り被検体から被検体への顕著な変動性に影響されやすい。それゆえ、このようなアプローチは、最良であっても循環系の挙動の大まかな近似を形成する希望をすることしかできず、したがって、これまでの所、生被検体内の血圧を正確に決定するには有効でないことが判明している。
【００２１】
前述に基づいて、必要であるのは、生被検体内の、血圧を含む、血行力学パラメータを査定する改善された方法と装置である。このような方法と装置は、理想的には非侵襲性であることになり、連続的に又はほとんど連続的に自己校正を行うことになり、かつ被検体が意識がある及び麻酔を掛けられているの両方のときのような、種々の異なった被検体生理的環境の下で有効である及び信頼性結果を生じるの両方であることになろう。最後に、このような改善された方法と装置は、査定される各特定被検体から測定されたパラメータに主として基づいていることになり、それによって各個に特有の校正を考慮することになろう。
【００２２】
（発明の概要）
本発明は、生被検体内の、血圧を含む、血行力学パラメータを査定する改善された方法と装置によって前述の必要を満足させる。
【００２３】
本発明の第１態様では、循環系内血圧を含む血行力学性質を査定する方法が開示される。本方法は、一般に、被検体の血管から第１パラメータを測定するステップ、その血管から第２パラメータを測定するステップ、第２パラメータに基づいて校正関数を導出するステップ、及び校正関数を使用して第１パラメータを補正するステップを含む。いったん校正されると、第２パラメータは連続的に又は周期的に監視され、そのパラメータの変化が注目の血行力学性質の変化を表示するために使用される。第１模範的実施の形態では、第１パラメータは圧力波形を含み、及び第２パラメータは血管内血液の全体流量運動エネルギーを含む。圧力波形の測定の間中、血管と循環系内血行力学性質に変化を誘発するように血管は圧平（圧縮）され、次いで、このような圧平の間中の運動エネルギーが測定され、かつ圧力波形中の１つ以上のアーチファクトを識別するために使用される。次いで、補正関数がこれらのアーチファクトに基づいて発生され、かつ測定された波形に適用されて、血管内の実圧を表す補正された又は校正された波形を発生する。第２模範的実施の形態では、校正関数を導出するために血管内の血流量の最高速度が音響信号を使用して決定される。
【００２４】
本発明の第２態様では、１つ以上の測定されたパラメータを使用して、生被検体の血管から得られた圧力信号を校正する改善された方法が開示される。一般に、本方法は、血管から圧力波形を測定すること、血管から少なくとも周期的に第２パラメータを測定すること、第２パラメータに基づいて校正関数を導出すること、及び導出された校正関数を使用して第１パラメータを補正することを含む。１模範的実施の形態では、本方法は、被検体の血管から圧力波形を測定すること、少なくとも１回同じ血管から第２パラメータを測定すること、第２パラメータに基づいて圧力波形中の少なくとも１つのアーチファクトを識別すること、測定された第２パラメータと少なくとも１つのアーチファクトと関連した少なくとも１つの性質とに基づいて校正関数を導出すること、血管内の血圧の校正された表示を発生するために少なくとも１回校正関数を圧力波形に適用すること、及び時間に従う血圧の変動を識別するために第２パラメータを連続的に監視することを含む。
【００２５】
本発明の第３態様では、生被検体の循環系の血行力学的反応を特徴表示する改善された方法が開示される。本方法は、一般に、或る種の条件の下で血管と関連した第１パラメータと第２パラメータとの間の第１機能関係を導出するステップ、これらの或る種の条件の下で非侵襲的に第１パラメータと第２パラメータを測定するステップ、測定されたパラメータの少なくとも１つ中の少なくとも１つのアーチファクトを識別するステップ、及び少なくとも第１機能関係と少なくとも１つのアーチファクトに基づいて第１パラメータの測定を拡大縮小するステップを含む。
【００２６】
本発明の第４態様では、誤り成分を有する血行力学パラメータ測定を校正する改善された方法が開示される。一般に、本発明は、血管と関連した血行力学パラメータを測定すること、第１パラメータと関連した誤り源を識別すること、誤り源に基づいて校正関数を発生すること、及び校正関数を使用して測定された血行力学パラメータを補正することを含む。１模範的実施の形態では、本方法は、血管から圧力波形を測定すること、呼吸作用に因る時間を通しての血管内の血液の運動エネルギー（又は最高速度）と関連した周期的変動を識別すること、時間を通しての運動エネルギーの変動との同期に基づいて校正関数を発生すること、及び周期的変動に対して波形を補正するために校正関数をその波形に適用することを含む。この呼吸作用は、圧力信号から、及び潜在的には他の信号からも充分に検出可能である。
【００２７】
本発明の第５態様では、生被検体の血管内の血行力学性質を測定する改善された装置が開示される。本装置は、一般に、血管と関連した第１血行力学パラメータを測定する第１トランスジューサ、血管と関連した第２血行力学パラメータを測定する第２トランスジューサ、及び第１トランスジューサと第２トランスジューサに動作上接続されて、第２トランスジューサによって生じた信号に基づいて校正関数を発生し、かつ補正関数を第１トランスジューサによって生じた信号に適用する信号プロセッサを含む。１模範的実施の形態では、血管は人間の橈骨動脈を含み、かつ本装置は、それの付近に非侵襲的に配置された圧力トランスジューサ、それの同じく付近にまた配置された音響トランスジューサ、血管を圧平するために使用される圧平デバイス、及び血管の圧平の間中圧力トランスジューサと音響トランスジューサからの信号を処理するプロセッサを含む。音響トランスジューサは、音響発射（ａｃｏｕｓｔｉｃ  ｅｍｉｓｓｉｏｎ）を血管内へ伝送しかつそれからの反響を受信する。次いで、圧平の間中の血液の速度及び／又は運動エネルギーに関する情報が、測定された圧力波形を校正するためにこの波形に印加される。
【００２８】
本発明の第６態様では、血行力学の査定、モデリング、及び校正についての前述の方法を実施する改善されたコンピュータ・プログラムが開示される。１模範的実施の形態では、このコンピュータ・プログラムは、Ｃ＋＋ソース・コード・リスティングのオブジェクト・コード表現を含み、このオブジェクト・コード表現はマイクロコンピュータ・システムの記憶デバイス上に配置され、かつマイクロコンピュータ・システムのマイクロプロセッサ上でランするように適合している。コンピュータ・プログラスは、マイクロコンピュータのディスプレイと入力デバイスに動作上結合されたグラフィカル・ユーザ・インタフェース（ＧＵＩ）を更に含む。マイクロコンピュータに供給された測定されたパラメータ・データに基づいて本明細書に説明された血行力学査定、モデリング、及び校正方法論を実施する１つ位上のサブルーチン又はアルゴリズムがプログラムに含まれる。第２模範的実施の形態では、コンピュータ・プログラムは、前述の血行力学測定装置と関連したディジタル信号プロセッサ（ＤＳＰ）の（埋込みプログラム・メモリのような）記憶デバイス内に配置された命令集合を含む。
【００２９】
本発明の第７態様では、前述の方法に従って得られたパラメータ・データを分析するかつ前述のコンピュータ・プログラムを利用する改善された装置が開示される。１模範的実施の形態では、本装置は、プロセッサ、不揮発性記憶デバイス、ランダム・アクセス・モリ、入力デバイス、ディスプレイ・デバイス、及び１つ以上のセンシング・デバイスに動作上結合された直並列データ・ポートを有するマイクロコンピュータを含む。分析下にある被検体から得られたデータは、直列又は並列データ・ポートを経由してマイクロコンピュータに入力される。記憶デバイス上に記憶されたコンピュータ・プログラムのオブジェクト・コード表現は、マイクロコンピュータのランダム・アクセス・メモリ内へロードされ、かつユーザによって入力デバイスを経由して入力されたコマンドと関連して入力データを分析する要求に従ってプロセッサ上で実行される。
【００３０】
本発明の第８態様では、前述の方法を使用して被検体に治療を施す改善された方法が開示される。本方法は、一般に、圧力データを測定するために有用な被検体の血管を選択するステップ、非侵襲的に被検体の圧力データを測定するステップ、校正関数を発生するステップ、血管内の血圧の校正された表示を生じるために校正関数を測定された圧力データに適用するステップ、及び校正された推定値に基づいて被検体に治療を施すステップを含む。１模範的実施の形態では、血管は人間の橈骨動脈を含み、かつ本方法は圧力トランスジューサを経由して橈骨動脈から圧力波形を測定すること、少なくとも１つの血行力学パラメータを測定するために音響波を使用すること、測定された血行力学パラメータに少なくとも一部分基づいて校正関数を導出すること、被検体内の１つ以上の医療状態を診断するために有用な血圧の校正された表示を導出するために校正関数を使用して圧力波形を校正すること、及び校正された表示に一部分基づいて被検体に対する治療の行程を用意することを含む。
【００３１】
本発明の第９態様では、血管内に配置された管腔（ｌｕｍｅｎ）の改善された検出方法が開示される。本方法は、一般に、血管内へ音響エネルギーを伝送すること、血管内でのエネルギーの伝搬の関数として音響エネルギーの反射の信号レベルを評価すること、及び反射内の低下した信号レベルの少なくとも１つの領域を識別することを含む。１模範的実施の形態では、パワー計量（ｐｏｗｅｒ  ｍｅｔｒｉｃ）がＡモード包絡線から導出される。深さの関数として後方散乱パワーの積分がアーチファクト（例えば、「プラトー（ｐｌａｔｅａｕｓ）」）を識別するために利用され、アーチファクトは管腔の深さを検出するために、適当な制約付きで、使用される。第２実施の形態では、後方散乱エネルギー（例えば、Ａモード信号）が、それらの中の低下信号レベルの領域を識別するために、事前選択された深さ間隔にわたって遂行されるパワー計算を使用して分析される。
【００３２】
本発明の第１０態様では、血管の少なくとも１つの壁を検出する方法が開示される。本方法は、一般に、血管内へ音響エネルギーを伝送すること、血管内の管腔と関連した少なくとも１つの領域を検出すること、及び管腔に対する血管の少なくとも１つの壁の位置を検出することを含み、位置を検出することは伝送された音響エネルギーから導出されたＡモード・データを分析することを含む。１実施の形態では、壁は、所与の深さでの積分されたパワーを平均管腔パワーと比較することによって検出される。他の実施の形態では、包絡線平方（ｅｎｖｅｌｏｐｅ  ｓｑｕａｒｅｄ）信号のレベル（振幅）が管腔のそれと比較される。
【００３３】
本発明の第１１態様では、血管の直径を決定する方法が開示される。本方法は、一般に、血管内へ音響エネルギーを伝送すること、血管内の管腔と関連した領域を検出すること、管腔に対する血管の第１壁の位置を検出すること、管腔に対する血管の第２壁の位置を検出すること、及び第１壁と第２壁の位置に少なくとも一部分基づいて血管の少なくとも部分の直径を決定することを含む。
【００３４】
本発明の第１２態様では、管腔検出を使用して生被検体と関連した１つ以上の血行力学パラメータを査定する改善された方法が開示される。本方法は、一般に、血管内へ音響エネルギーを伝送すること、血管と関連した管腔領域を検出すること、管腔の検出に基づいて血管と関連した少なくとも１つの壁を検出すること、及び注目の血行力学パラメータを推定するために印加されたストレスの関数として少なくとも１つの壁の位置の変化を分析することを含む。１模範的実施の形態では、本方法は、（管腔を経由して）前動脈壁と後動脈壁を検出するためにＡモード信号を分析し、かつ動脈心拡張期血圧を決定するために圧力波形上に圧平圧力の関数として動脈直径の変動をマップすることを含む。
【００３５】
本発明の第１３態様では、生被検体の組織内の血管位置に関する改善された方法が開示される。本方法は、一般に、音響エネルギーを発生すること、血管の付近の被検体の組織内へ音響エネルギーを伝送すること、血管による後方散乱エネルギーの少なくとも部分を受信すること、後方散乱エネルギーに基づいて血管と関連した少なくとも１つのアーチファクトを検出すること、及び検出されたアーチファクトに少なくとも一部分基づいて血管の位置を決定することを含む。１模範的実施の形態では、血管は人間の橈骨動脈を含み、及び音響エネルギーは被検体の手頸／前腕の内部領域を経由して動脈内へ伝送される超音波エネルギーを含む。動脈壁と管腔から後方散乱されたエネルギーは、超音波トランスジューサによって受信され、かつ動脈に対するそのトランスジューサの位置及び／又は最適位置を識別するために積分されたパワー計量を使用して分析される。
【００３６】
他の実施の形態では、積分されたパワー及び／又は信号レベルは、管腔の位置を検出するために組織内深さの関数として分析される。次いで、動脈壁の位置が管腔に対して決定され、それによって、動脈の直径の測定をできるようにする。動脈に対するトランスジューサの位置は壁分離（直径）を最大限にするように調節され、それによって動脈の直上にトランスジューサ（及びいずれかの関連圧平デバイス）を有効に維持する。血管に対していずれか他の所望配向で、例えば、固定オフセットを伴って、血管の直径の関数として変動するオフセットを伴って等、トランスジューサを維持するために、又は信号品質のような他のパラメータ、嚢胞性構成要素（ｃｙｓｔｉｃ  ｃｏｍｐｏｎｅｎｔ）又はクラッタ（ｃｌｕｔｔｅｒ）、等々の存在に基づいて、トランスジューサの位置を制御するために、アルゴリズムがまた調節されることがある。
【００３７】
（発明の詳細な説明）
図面を参照するが、これらの図面を通じて同様の番号は同様の部品を指す。
【００３８】
注意するのは、本発明は人間被検体の橈骨動脈（すなわち、手頸）を経由して循環系の血行力学パラメータを査定する方法と装置の面で本明細書に説明されるが、本発明は人体の他の位置でこのようなパラメータを監視するばかりでなく、他の温血種上のこれらのパラメータを監視するためにもまた具体化される又は適合しているといえることである。全てのこのような適合及び代替実施の形態は、本明細書に添付された特許請求の範囲に属すると考えられる。
【００３９】
概  観
１基本的態様では、本発明は、生被検体の循環系上に人工的に「ストレス」を誘発することによってその被検体内の血行力学パラメータを査定する方法を含む。これらのストレスに対する循環系の反応は、知られているか又は決定可能であり、かつ観察されたデータで以てアーチファクト又はマーカを識別するのに有効である。これらのマーカは、その後、前述の血行力学パラメータの測定を校正するために使用される。
【００４０】
例えば、下により詳細に説明するように、本発明は、人間の橈骨動脈をおおって配置された圧力測定法センサ又は表面圧力センサから得られた血圧波形を校正するのに有効であり、未校正の圧力波形は、橈骨動脈自体内に実際に経験される波形から潜在的に実質的に変動する。１実施の形態では、動脈に加えられた「ストレス」は圧平（すなわち、圧縮）であり、かつ圧平のエリアを通して流れる血液の速度は速度プロファイル中のマーカを識別するために監視される。これらのマーカは、なかでも、動脈壁を横断する近ゼロ経壁圧の状態に相当する。このようにして、真動脈圧の正確測定値が非侵襲的に得られることがある。しかしながら、納得するように、本明細書に説明する本発明は、血管内の２つの位置間の圧力差、静脈壁又は動脈壁のコンプライアンス、心室収縮の強度の変動、等々のような他の血行力学性質を査定するのにまたたやすく使用され、したがって、動脈血圧の測定に限定されない。
【００４１】
血行力学性質を査定する方法
図３を参照して、本発明に従って循環系内血圧を含む血行力学性質を査定する方法を説明する。図３に示したように、本方法３００の第１ステップ３０２は、被検体の血管から第１パラメータを測定することを含む。本文脈の中では、測定されるパラメータは、本明細書に図３ａに関して下に更に詳細に説明するように、被検体の橈骨動脈の付近に配置された圧力センサ又は圧力トランスジューサから導出された血圧波形であることになる。しかしながら、納得するように、先に注意したように他の血行力学パラメータもまた測定される。第１パラメータの測定に黙示されているのは１つ以上の誤り源の存在である。すなわち、パラメータの測定された値は、被検体の循環系に存在するパラメータの実値の完全なる表示ではない、つまり、それと異なっている。動脈血圧の場合には、実値は、侵襲的動脈カテーテル法（ｉｎｖａｓｉｖｅ  ａｔｅｒｉａｌ  ｃａｔｈｅｔｅｒｉｚａｔｉｏｎ）のＡ線（Ａ−ｌｉｎｅ）技術又は「ゴールド・スタンダード（ｇｏｌｄ  ｓｔａｎｄａｒｄ）技術によって測定されることがあるように、動脈自体内に存在する値である。このような誤り又は差の理由は、図３ａを参照して下に更に詳細に論じる。
【００４２】
次に、図３のステップ３０４で、ストレスが血管上に誘発され、このことが血管の血行力学性質を（少なくとも局所的に）変更し、それによって血管又は循環系と関連した他のパラメータに全体的に変化を誘発する。図３ｂに関して論じるように、このストレスは、１実施の形態では、時間の関数として血管を圧平する又は不定に圧縮し、それによって、なかでも、圧平の領域内の血液の質量流量（Ｑ）、速度（ｖ）又は速度勾配、及び運動エネルギー（ＫＥ）に変化を誘発する。しかしながら、注意するのは、例えば、動脈断面積に影響する（カフ状デバイスで以て起こるであろうような）環状面閉塞（ｃｅｒｃｕｍｆｅｒｅｎｔｉａｌ  ｏｃｕｌｓｉｏｎ）、又は動脈のコンプライアンスに影響する被検体内への化学物質の局所化導入のような、類似の又は他の血行力学性質に影響するために、先に説明した圧平ストレス以外のストレッサ（ｓｔｒｅｓｓｏｒｓ）を被検体に印加することもあることである。多くのこのようなストレッサ／血行力学パラメータ組合わせは、本発明と両立して使用されるといえる。
【００４３】
ステップ３０６で、下のステップ３０８での校正関数の導出を容易にするために血管と関連した第２パラメータが測定される。本明細書に図３ｃに関してより詳細に論じるように、第２パラメータは、１実施の形態では、全血流量運動エネルギーを含む。というのは、このパラメータはステップ３０４におけるストレッサの印加の関数として或る種の容易に識別される「アーチファクト」を示すからである。本明細書に使用されるように、用語「アーチファクト（ａｒｔｉｆａｃｔ）」とマーカは、同義に使用され、かつデータ集合内に存在するいずれもの識別可能特徴又は関係を称する。しかしながら、同じ又は他のアーチファクトを示す他のパラメータも、例えば、最高血液速度、最大血管断面積、及び最大血質量流量を含む校正関数を導出するために使用されることがある。
【００４４】
図３のステップ３０８で、次に、校正計量又は校正関数がステップ３０６で導出されたパラメータ情報に基づいて導出される。特に、１つ以上のアーチファクト又はマーカがパラメータ・データ中で識別され、これらのアーチファクトは上のステップ３０２の第１パラメータの実値と測定された値との間に或る種の関係が存在するときを表示する。本明細書で図３ｄを参照して論じるように、校正関数を導出するプロセスの１実施の形態は、印加されたストレッサ（圧平）の領域内の全血流量運動エネルギーを測定し、かつ（血管の断面積のパーセンテージ減少に相関した）圧平の関数として心収縮期速度及び／又は心拡張期速度のプロファイル中の変化を識別する。
【００４５】
図３の方法のステップ３１０では、補正された又は校正された測定を発生するために、ステップ３０８で導出された校正関数がステップ３０２の第１パラメータの測定に適用される。第１パラメータが時間の関数として連続的に（又は周期的に）測定されるならば、ステップ３０８の相関関数は適当であるに従って連続的に又は周期的に適用され、それによって進行様式又は連続様式で第１パラメータの校正された測定を発生することに注意されたい。しかしながら、種々の異なったファクタに因り、実「Ａ線」動脈血圧と必要とされた校正関数の拡大縮小又は大きさの両方は時間の関数として変動することがあり、それゆえ、先に計算された校正関数に基づいたいずれの「校正された」測定も誤りがあることになる。１代替では、ストレッサを再印加（例えば、他の圧平掃引を遂行）し、更新された補正関数を発生し、かつこれを第１パラメータの測定された値に適用することによって、ユーザは簡単に周期的に再校正することもできる。
【００４６】
しかしながら、本明細書で下に更に詳細に説明するように、本発明は、都合良く、第１時刻ｔ１に校正関数を発生し、次いで、測定されたパラメータの変動を表示するために第２血行力学パラメータ（例えば、最高速度、最大運動エネルギー、最大面積、又は最大流量）を連続的に監視する能力を備える。これは、更に下で説明するように、真動脈圧と外圧との間に所定関係を確立するように動脈に印加される外圧を制御することによって方法３００のステップ３１２で完遂される。
【００４７】
ステップ３１２で、動脈内で動脈壁を横断する圧力（すなわち、「経壁圧」）を所望値に又は近くに維持するように、動脈に印加される圧力が第１パラメータの選択された値へ制御される。このプロセスは、本明細書では、特定値へ「サーボする（ｓｅｒｖｏｉｎｇ）」と称す。本明細書に図５ａ〜５ｄを参照して詳細に論じるように、このサーボすることは圧平のエリアに特定血流量運動エネルギーを発生し、次いで、この運動エネルギーの変化が真動脈圧の変化を識別するために使用される。この「連続校正」は、本発明の望ましい属性である。というのは、被検体の血管を用いた血行力学の連続的正確測定は、特に血圧が被検体の治療における瞬間瞬間の決断の基礎として使用されるような外科又は他の生命治療展開の文脈の中で決定的に重要である。
【００４８】
図３ａを参照して、生被検体内の１つ以上の血行力学パラメータを測定する（図３のステップ３０２）方法の１実施の形態を説明する。図３ａの方法３２０の第１ステップ３２２は、測定用１つ以上の血行力学パラメータを選択することを含む。測定されるパラメータの選択は、被検体の血圧又は動脈に関する狭窄症の酷しさのような、査定される状態のだけでなく、また下のステップ３２４で選択される監視位置の関数である（すなわち、離散点に肉体的に位置している動脈内の局所化狭窄症の場合のように、或る種のパラメータは或る定まった位置でのみ測定されることがある）。
【００４９】
次に、ステップ３２４で、被検体の体内の血管が監視用に選択される。その立入り易さ（ａｃｃｅｓｓｉｂｉｌｉｔｙ）及び皮膚の表面に比較的近接していることに因り、人間の橈骨動脈は循環系内の血行力学パラメータを監視するための優れた位置である。もっとも、承知のように人間（又は他の種）上の他の位置もこの目的に使用されることがある。上に注意したように、監視する位置はまた、査定される状態の又は遂行される監視の型式に関係する又はこの型式によって決定されることがある。もちろん、多数の監視位置が、逐次又は並列であろうと、本発明の方法で以て採用される。
【００５０】
人間の橈骨動脈に関しては、更に注意するのは、橈骨動脈が糖尿病に因る狭窄症及び石灰化（ｃａｌｃｉｆｉｃａｔｉｏｎ）を含む動脈疾患によって最小限にしか影響されないと逸話的証拠が示唆していることである。この観測された挙動に対する理由は、この論議の範囲を超えている。しかしながら、この挙動は、図３ｂに関して本明細書に行われる圧平ストレスについての論議にいくぶん重要である。というのは、内側石灰化（ｍｅｄｉａｌ  ｃａｌｃｉｆｉｃａｔｉｏｎ）のような事前存在動脈疾患があることは動脈血圧を正確に測定する能力に強い影響を与えるおそれがある。印加ストレスとして制御された圧平を利用する血圧測定を遂行するとき橈骨動脈を選択することによって、ユーザは、事前存在する狭窄症や石灰化に関係する多くの潜在的誤り源から有効に隔離される。
【００５１】
次に、図３ａのステップ３０２で、選択されたパラメータを測定又は検知する能力のある１つ以上のパラメータ・センサが選択された血管の付近に配置される。人間の橈骨動脈上の血圧を測定する場合、圧力センサ（トランスジューサ）が、橈骨動脈の真上にあるように、手頸の内側表面上で皮膚と接触して物理的に配置される。トランスジューサは、それの表面に印加された圧力に既知の関係で圧力信号を生じる能力のある良く理解された圧電型、又はその他の型式の１つであってよい。最適信号結合と検知が達成されるというようにトランスジューサを位置決めする方法と装置は、血圧測定技術上また周知であり、したがって、本明細書で更に説明しないことにする。手頸の皮膚と接触している間トランスジューサは、下に横たわる組織／動脈の低又はゼロ圧縮の状態に初期的に維持され、理由は本明細書に更に充分に説明する。
【００５２】
ステップ３２８で、信号が時間の関数としてトランスジューサから測定される。信号は、監視の所望頻度に依存して、離散的に（例えば、所定間隔で）又は連続的に測定されることがある。先に説明した模範的圧力トランスジューサの場合は、連続測定によっての出力信号は、時変（ｔｉｍｅ  ｖａｒｉａｎｔ）波形を含むことになる。動脈血圧の場合、波形は、一般に、実「ゴールド・スタンダード」動脈圧を追跡することになるが、心周期（ｃａｒｄｉａｃ  ｃｙｃｌｅ）の種々の位相に従う圧力変化に伴って変動する誤り又はオフセットをなお含むことになる。この時変的、非線形誤り、つまり測定された圧力波形と実血圧波形との間の「可変誤り」は測定プロセスに追加の複雑性を出現させるが、これは、下に更に詳細に説明するように、本発明がそれを克服するのに特に適合している１つである。
【００５３】
図３ｂを参照して、被検体の循環系に１つ以上のストレスを誘発する（図３のステップ３０４）方法の１実施の形態を詳細に説明する。方法３３０の第１ステップ３３２で、印加されるストレスが選択される。本明細書に使用されるように、用語「ストレス」（又は「ストレッサ」）は、被検体の循環系内に人工的に誘発されるいずれかの肉体的又は生理的変化を称する。本実施の形態では、印加されるストレスは、圧平、つまり時間の関数としての選択された（橈骨）動脈の肉体的圧縮を含む。本明細書に使用された、圧平「掃引」は、一般に、動脈の上に横たわる皮膚の表面に全体的に垂直な方向に動脈へ増大する又は減少する圧力を定常印加することを称する。圧平の概念は、一方の手頸を他方の手の親指と人指し指との間に、親指をその手頸の内側部分の頂上に置いて、その動脈が閉塞されるまで、橈骨動脈上の圧力をゆっくりと増すことによって簡単に例示される。しかしながら、納得するように、一般的な主張（ｐｒｏｐｏｓｉｔｉｏｎ）としては、本明細書に使用される圧平は、（ｉ）時間を通して圧縮を増加又は減少させる連続的線形速度で、（ｉｉ）時間を通して圧縮を連続的非線形（例えば、対数関数的に）増加又は減少させる、（ｉｉｉ）不連続的又は区分的連続的線形又は非線形圧縮、（ｉｖ）圧縮と弛緩を交互させる、（ｖｉ）正弦波関数又は三角関数、又は（ｖｉ）（「ランダム・ウォーク」のような）ランダム運動のような、種々の異なった形を取ることがある。全てのこのような形は、用語「圧平」に包含されると考えられる。
【００５４】
図４ａと４ｂを参照して、平圧の血行力学的作用を詳細に説明する。直ちに納得するように、人間の橈骨動脈のような動脈４００の圧平を増す結果、その動脈の有効断面を減少させる。先に説明した動脈狭窄症に類似して、圧平４０２は、動脈内の流体面積（ｆｌｏｗ  ａｒｅａ）を減少させ、それによって一定容積流を維持するために絞り（ｒｅｓｔｒｉｃｔｉｏｎ）を通して血液速度（ｖ）が上昇する結果となる。この関係は、流体力学技術で良く理解されている。減少する流体面積を横断しての速度のプロファイルは、図４ａの速度勾配４０４によって例示されたように、また充分に変更される。それゆえ、圧平の結果、絞り流体面積を横断して、より高い最高速度、より急な速度勾配、かつより大きなエネルギー又はより急な圧力勾配を生じる。
【００５５】
図４ｂは、動脈の流体面積のパーセント減少の関数としての動脈内のピーク又は最高流速を例示する。全血流量運動エネルギーが同様に面積に関係するのは速度とのその関係に一部分因っているが、もっとも導出するのは下により詳細に説明するように或る程度もっと困難である。現在の議論において簡単化のために、流体面積のパーセント減少は組織（皮膚）表面に印加された圧平圧力に正比例すると想定するが、もっとも実際にはこの関係は下に更に説明するように実質的にもっと複雑であることに注意されたい。更に、図４ｂは、「定常状態」動作を一般的に例示しており、かつ、下で更に詳細にまた論じる、例えば、正規心周期に因る圧力の変動の影響は調べない。
【００５６】
図４ｂに例示したように、低圧平圧力の領域４１０では、流体面積のパーセンテージ減少は小さく、かつ流速と勾配への影響は最小限である。容積血流量（Ｑ）は、影響されない。圧平圧力が増大する（領域４１２）に連れて、流体面積は更に減少させられ、かつ容積流量が維持される間に血液速度、速度勾配、圧力勾配、運動エネルギーは相当して増大し始める。圧平圧力が更に増大するに連れて、流体面積は実質的に減少させられ、速度、速度勾配、血厚勾配、及び運動エネルギーは実質的に上昇、増大する一方、正規代謝要求の下での容積流量を依然維持する。狭窄症動脈では、この領域４１４は、「亜危険（ｓｕｂ−ｃｒｉｔｉｃａｌ）」狭窄症、すなわち、被検体の過剰容量（ｅｘｃｅｓｓ  ｖｏｌｕｍｅｔｒｉｃ  ｃａｐａｃｉｔｙ）が顕著に減少させられ、一般に、付随症状が少ない又はない狭窄症のレベルに相当する。「亜危険」という名称は、正規代謝要求の下で適当な組織灌流があることで以て無症候性である事実を称する。
【００５７】
流体面積が更に減少すると、「危険」領域４１６として知られている所を通して遷移を生じる。危険領域では、流体面積が非常に減少させられるので高まった流体抵抗を克服するには不充分なエネルギーしかなく、かつ容積流量はもはや維持されない。これらの領域４１４、４１６の間で、速度「ピーク」４２０が形成される。逸話的証拠は、このピーク４２０が動脈直径の大略５０％減少点（流体面積の大略７５％減少に相当する）で起こることを示唆する。結果として、狭窄症面積に対して抹消の血液速度、容積流量、及び流量運動エネルギーは、流体面積が更に減少するのと共に急激に低下する。動脈が全く閉塞されかつ流体面積がゼロに接近する（領域４１８）に連れて、容積流量（Ｑ）は、血液速度と血流運動エネルギーと同じように、ゼロに接近する。
【００５８】
図４ｂを調べることは、１印加ストレスに対する循環系の１反応を特徴表示する面で重要な情報を生じる。特に、圧平の関数としての速度の挙動、かつ最も著しくは速度プロファイル中の最高速度の上昇は、動脈内圧が外部圧平によって動脈の壁に印加された圧力と事実上等しい点の識別をできるようにする。この状態は、本明細書では、「ゼロ経壁圧」の状態と称する。図４ａと４ｂに例示した圧平掃引の間中、動脈壁の外側に印加された外部圧力が動脈内圧によってちょうど相殺される点に到達する。この点に到達するまでは、流体面積に顕著な減少（かつ、先に説明したように、その結果に付随する速度や、速度勾配や、容積流量（Ｑ）や、運動エネルギーの変化）は、起こらない。しかしながら、印加圧力が動脈内圧を超えると、動脈の直径と断面積は減少され始め、最高流速と速度勾配が上昇、急峻化し始める（図４ｂの領域４１４）。最高速度（と運動エネルギー）の上昇は、本実施の形態では経壁圧が概略的に平衡する点の「マーカ」として使用される。
【００５９】
しかしながら、先に論じたように、循環系は静止系ではなく、むしろ動的でありかつ、正規心周期ばかりでなく、これと（下に論じる）呼吸のような他のファクタの両方に因って、顕著な動脈内圧力揺動（ｉｎｔｒａ−ａｔｅｒｉａｌ  ｐｒｅｓｓｕｒｅ  ｆｌｕｃｔｕａｔｉｏｎ）を受ける。それゆえ、血行力学性質を、特に動脈内圧を測定するとき、このような圧力揺動をまた、考慮しなけれならない。
【００６０】
図５ａ〜５ｄを参照して、前述の動的圧力揺動の下での循環系の反応を説明する。図５ａは、人間の正規心洞調律（ｎｏｒｍａｌ  ｓｉｎｕｓ  ｃａｒｄｉａｃ  ｒｈｙｔｈｍ）５００を例示する。医療技術上周知であるように、この洞調律５００中に心臓内の種々の心室機能に相当する心収縮期５０２と心拡張期５０４の両方がある。これらは、洞調律５００中の、及び動脈内圧の最大値と最小値を事実上表す。
【００６１】
図５ｂは、図５ａについて前述した洞調律５００の、及び図５ｃによる動脈に印加された外部圧平圧力の関数としての動脈壁の変位を例示する。２つの対向する動脈壁５１０、５１２を明確のために図５ｂに例示する。もっとも、これらは圧力反応の面では事実上互いに鏡像である。技術上周知であるように、動脈壁は、典型的に（健康な人間では）顕著な弾性（ｅｌａｓｔｉｃｉｔｙ）と反発性（ｒｅｓｉｌｉｅｎｃｙ）を有する実質的にコンプライアンスな管（ｖｅｓｓｅｌ）である。それゆえ、管内の圧力を増すと、動脈の対向する壁５１０、５１２は、外向きに反れる傾向があって、膨張している風船のように、大いに動脈の直径を増大する。同様に、動脈内圧を減少すると、動脈壁の反発性が直径を減少させる。人間の動脈が、正規心周期の間中或る程度周期的に拡張及び収縮することは周知である。
【００６２】
図５ｂに示したように、動脈内血圧の変動は、動脈の壁を心収縮期圧５０２に相当する最大直径５１６へと外側へ向けて反らし、及び動脈壁を心拡張期圧５０４に相当する最小直径５１８へと凹ませる（ｃｏｌｌａｐｓｅ）。動脈の外側に印加された圧平圧力（図５ｃ）が増大するに連れて、ゼロ経壁圧の先に説明した状態に心収縮期圧と心拡張期圧の両方によって逐次到達する。特に圧平圧力を増大するに伴って、心拡張期（より低圧）状態５２０でのゼロ経壁圧がまず達成され、これに心収縮期（より高圧）状態５２２でのゼロ経壁圧が続く。圧平圧力をゼロ経壁圧状態を超えて増大する際の、心拡張期状態５２０をまず考察すると、動脈の有効直径（と流体面積）が漸次減少し始め、その結果、先に説明したように、血流量勾配が急峻になり、ピーク血液速度が上昇し、及び運動エネルギーが増大する。圧平圧力を心拡張期圧の上へ充分に増大するに連れて、動脈は、心周期の心拡張期部分の間中点５２７でより完全に閉じる。同様に、圧平圧力を増すと、動脈の直径は、心収縮期状態５２２でまた減少し始めて、動脈が心周期の心拡張期部分５２７と心収縮期部分５２９の両方の下で完全に凹むまで、同様の結果を生じる。前述の挙動に基づいて、心拡張期状態５２０と心収縮期状態５２２の両方に対して、最高血液速度の変動とパーセント流体面積の減少に関して、２つの曲線を作図することができる（図５ｄ）。心収縮期状態５２２の速度「ピーク」５２４は心拡張期状態５２０に対する相当するピーク５２６よりも圧平の高いレベルで起こり、それは大きな外圧は後者に反して前者では動脈を凹ませるために印加されなければならないからであることに注意されたい。更に注意するのは図５ａの心収縮期最大と最小との間及び心拡張期最大と最小との間にある圧力で、図５ｄのそれらと類似の曲線の系列を作図することができ、曲線のこのような系列は全心周期の間中に動脈の挙動及び関連した血行力学パラメータを特徴表示するのに有効であることである。
【００６３】
図４ｂの速度曲線で以てのように、図５ｄの曲線は、心周期の心拡張期部分と心収縮期部分の両方（又はそれらの間の他のいずれもの部分）の間中、ゼロ経壁圧が達成される点を圧平掃引の間中マークするために有用である。この情報の効用と応用は、本明細書に図３ｄ〜３ｃを参照して詳細に説明する。
【００６４】
圧縮ストレス又は圧平ストレスの前述の模範的印加は、なかでも、選択された動脈内血圧の測定に有用であるが、納得するように、ストレスの他の型式を循環系内に反応を誘発するために印加することもある。これらのストレスと関連したアーチファクト又は「マーカ」は、圧平ストレスに対するものと全体的に類似した様式で、すなわち、マーカ又は既知の関係の存在を一般には循環系又は特には血管内の或る種の血行力学状態と相関させることによって、利用することもできる。それゆえ、図３ｂの方法は、圧縮ストレスの使用に決して限定されない。
【００６５】
図３ｂに再び帰って、ストレスを選択された血管に印加する方法３３０の第２ステップ３３４は、このようなストレスを印加することができる機構を用意することを含む。上に説明した圧平の文脈の中で、動脈を圧平する手段として（圧力波形を測定するために使用された）前述の圧力トランスジューサを使用することに特に効用がある。というのは、この配置（ａｒｒａｎｇｅｍｅｎｔ）は測定を精確に圧平の点で行うことをできるようにするからである。この型式の圧平機構は、本明細書に図１２に関して説明する。しかしながら、承知のように、分離圧力トランスジューサと圧平機構、又は他の構成（ｃｏｎｆｉｇｕｒａｔｉｏｎ）も本発明と一緒に使用されることがある。
【００６６】
図３ｂのステップ３３６で、用意された機構は、選択された動脈にストレスを印加するために利用される。圧平の特殊な場合に、圧力の連続的に高まるレベルで圧力トランスジューサが手頸の皮膚に対してアサート（ａｓｓｅｒｔ）されるというように、先に説明した圧平「掃引」が印加され、それによって下に横たわる動脈を圧縮する。図３ａの方法３２０で以てのように、動脈をおおう圧平デバイスの最適配置と最適配向は、いずれか種々の充分に理解された先行技術を使用して決定することができ、これらの技術は、讓受人の、１９９９年、６月２９日に出願された「動脈血圧の非侵襲的決定方法と装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐａｒａｔｕｓ  Ｆｏｒ  Ｔｈｅ  Ｎｏｎｉｎｖａｓｉｖｅ  Ｄｅｔｅｒｍｉｎａｔｉｏｎ  Ｏｆ  Ａｒｔｅｒｉａｌ  Ｂｌｏｏｄ  Ｐｒｅｓｓｕｒｅ）」と題する同時係属米国特許出願第０９／３４２，５４９号、及び２００１年３月２２日に出願された「血管位置を含む血行力学パラメータの非侵襲的査定方法と装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐａｒａｔｕｓ  ｆｏｒ  ｔｈｅ  Ｎｏｎｉｎｖａｓｉｖｅ  Ａｓｓｅｓｓｍｅｎｔ  ｏｆ  Ｈｅｍｏｄｙｎａｍｉｃ  Ｐａｒａｍｅｔｅｒｓ  Ｉｎｃｌｕｄｉｎｇ  Ｂｌｏｏｄ  Ｖｅｓｓｅｌ  Ｌｏｃａｔｉｏｎ）」と題する米国特許出願第０９／８１５，９８２号に説明した技術を含み、これらの両方は言及することによってそれらの内容を全面的に本明細書に組み入れてある。しかしながら、注意するのは、前述の方法３２０は、（例えば、治療個別管理（ｉｎｄｉｖｉｄｕａｌ  ａｄｍｉｎｉｓｔｅｒｉｎｇ  ｔｒｅａｔｍｅｎｔ）に従う手動による）非最適トランスジューサ配置で以てさえも、このような場合の信号結合が妥当である限り、多分利用されることである。
【００６７】
図３ｃを参照して、校正関数の導出を容易にするために血管と関連した第２血行力学パラメータを測定する（図３のステップ３０６）方法の１実施の形態を説明する。方法３４０の第１ステップ３４２で、測定される第２血行力学パラメータが選択される。このパラメータの選択は、或る点では、測定される第１血行力学パラメータの選択（図３ａ）ばかりでなく、循環系へのストレスの印加（図３ｂ）に結合されている。先に説明したように血圧測定及び圧縮ストレスの使用（圧平）の文脈の中で、いくつかの「第２」血行力学パラメータは、思うに、血液速度、全血流量運動エネルギー、及び血容積流量ばかりでなく、それのらのいずれかの変形又は組合わせを、これらに限定することなく、含む校正関数を発生するために使用されることがある。全血流量運動エネルギーは、測定するために特に有用な１つのパラメータである。というのは、それがゼロ経壁圧状態又はその他の有用な関係の１つ以上のたやすく観察可能なマーカを含むからである。全血流量運動エネルギーは、或る種の他のパラメータよりもまた誤りを生じることが少ない。というのは、それが全血管を横断して得られた速度情報ばかりでなく、振幅情報を利用するからである。更に、運動エネルギーのピーキング（ｐｅａｋｉｎｇ）は、最最高速度のような他のパラメータのピーキングよりもめざましい。
【００６８】
次に、ステップ３４４で、選択された「第２」パラメータが適当なセンサ又は測定技術を使用して測定される。運動エネルギーの測定又は血液速度の測定の場合、一般にこれらのパラメータを非侵襲的に測定するいくつかの周知の技術が存在する。特に注意するのは、血液速度を測定するために音響エネルギー（例えば、超音波）を使用することである。特に、音響測定技術は、一般に、速度測定上周知のドップラー原理を採用し、血管内を流れる血液によって反射された反響と関連した周波数偏移が血液速度の測定を行うために分析されるようになっている。音響による血液速度測定の多くのいろいろな変形が存在し、これらには連続波（ＣＷ）及び音響パルス（パルス・ドップラー）の使用がある。このような技術は、周知かつ充分に理解されており、したがって、更にここでは説明しないことにする。
【００６９】
同様に、音響測定技術は、被検体血管内を流れる血液の運動エネルギーの測定を導出するために使用されることがある。注意するのは、圧平の間中血管内に作り出される複雑な血液速度勾配（図４ａ）の結果、血管内の血液の運動エネルギーの計算は、全体として、上に説明した最高血液速度の平方に単に比例するのではない。むしろ、その運動エネルギーの推定は、この勾配の複雑性を捕らえる加算技術又は積分技術の適用を必要とすることである。血液運動エネルギーを計算するこのような加算技術／積分技術は、技術上周知であり、したがって、本明細書に更に説明しない。
【００７０】
他の実施の形態では、所望マーカが示されるに当たっての圧平（外部）圧力は、先に本明細書にその内容を組み入れてある讓受人の同時係属米国特許出願第０９／３４２，５４９号に説明された時間−周波数方法論を使用して決定されることがある。この時間−周波数方法論を使用して、経壁圧がゼロに等しいに当たっての圧平圧力は、動脈内で反射された音響エネルギーの時間−周波数表示を作図することによって決定することができる。時間−周波数分布が最大限であるとき、ゼロ経壁圧状態が達成される。それゆえ、時間−周波数分布は、本発明の目的のためになお他のマーカとして働く。
【００７１】
なお他の実施の形態では、第２血行力学パラメータを監視するためにいわゆる音響「Ａモード」が使用されることがある。このアプローチでは、音響波が発生されかつ血管内へ伝送され、伝送からの反射つまり反響が受信されかつ伝送の時刻と受信の時刻との間の関係を決定するために分析される。このような分析を通して、時間上の異なった点、及び心周期内の異なった点での動脈の相対直径を決定することができる。周知の時間領域反射率計（ＴＤＲ）に類似して、Ａモード技術は、１実施の形態では、異なった音響性質の材料間の種々の境界（例えば、「近」動脈壁／組織境界、「近」動脈壁／血流（ｂｌｏｏｄ  ｓｔｒｅａｍ）境界、血流／「遠」動脈壁境界、等）を横切る音響波の遷移によって発生された反射を利用する。特に、これらの反射の相対タイミングが種々の境界間の距離を決定するために分析される。異なった媒体を通しての音響波の伝搬速度を知って、反射境界（すなわち、組織厚さ、動脈直径等）を決定することができる。図５ｂによって、（心周期の心収縮期部分と心拡張期部分の両方の下で）動脈壁の反れが圧平圧力の関数として変動することを思い出すならば、動脈直径（及びこれに関係した面積）の変化は、上昇する最高速度を使用することと類似して、ゼロ経壁圧状態又は重要な他の状態を識別するためにこれらの状態の「マーカ」として使用することもできる。
【００７２】
Ａモード音響を使用して血行力学パラメータを査定する他の技術もまた採用することができる。このような代替技術は、本明細書に図１７〜３９に関して更に詳細に論じる。
【００７３】
更に納得するように、音響信号技術上周知である「Ｍモード」（運動モード）と「Ｂモード」（輝度モード）を含み、他の音響モード性（ａｃｏｕｓｔｉｃｍｏｄａｌｉｔｙ）も本明細書に説明した本発明と一緒に採用することができる。
【００７４】
第２血行力学パラメータ及びそれと関連したマーカを測定する前述の実施の形態の各々での音響波の使用にかかわらず、納得するように、他の非音響技術がこのようなマーカを識別するために応用されることがある。例えば、干渉計法（ｉｎｔｅｒｆｅｒｏｍｅｔｒｙ）を使用するような、動脈直径／面積を正確に測定する他の方法がゼロ経壁圧状態を識別するために採用されることがある。全てのこのような技術は、本発明の範囲内に入ると考えられる。
【００７５】
図３ｄを参照すると、所望血行力学状態が達成される当たってのストレッサ大きさ（例えば、ゼロ経壁圧が達成されるに当たっての圧平圧力）が実圧つまり真動脈圧に相関させられる。印加されたストレスと血管によって実際に感じられたストレスとの間に高度の結合がある簡単な場合には、測定されたストレスを実ストレスと等しいとすることができる。特に、圧平（圧縮）ストレスが印加されるに当たっての動脈血圧測定の文脈の中で、圧平デバイスによって印加されかつ関連した圧力トランスジューサによって検知された圧力は、アーチファクト状態又は「マーカ」状態が観察されるとき実動脈圧と等しいとすることがきできる。例えば、増大する血液運動エネルギーが先に論じたようにゼロ又は近ゼロ経壁圧の状態に相関するならば、運動エネルギーのこのような増大が観察されたときに動脈壁に印加された圧力は真動脈内圧に等しいであろう。それゆえ、印加点（例えば、皮膚）と動脈壁との間に印加された圧力が非常に高かったならば、印加点に印加された圧力は動脈壁に印加された圧力に近似しているであろう、したがってまた、動脈内圧に近似しているであろう。
【００７６】
しかしながら、先に論じたように、動脈壁と圧力トランスジューサとの間に介在した組織、腱、及び皮膚は、多くの場合、トランスジューサ（又は圧平機構）によって印加された圧力と動脈壁によって実際に感じられた圧力との間に複雑な関係を作り出す。簡単に述べると、皮膚に印加された圧力のいくらかは、この介在した材料を圧縮するために使用される。それゆえ、外部的に印加された圧力の部分のみが実際には動脈壁によって感じられる。更に、注意するのは、組織は血管の下及び骨の上にまた存在し、同じくまたこの組織を圧縮するのにいくらかの損失が起こる。
【００７７】
したがって、動脈圧の測定が外部的に印加された圧力（例えば、圧力トランスジューサによって測定された圧力）に基づくとき、所与の被検体についてその組織コンプライアンス及び結合の程度に依存して、このような測定に或る量の誤りが導入されることになる。
【００７８】
この問題への１先行技術アプローチは、介在した材料の反応を（例えば、線形力定数を有するばねの系として）モデルし、かつ圧力トランスジューサによって測定された圧力をこのモデルに基づいて補正することであった。しかしながら、このアプローチは使用されたモデルほどに良いに過ぎない。異なった組織厚さ値、密度値、及びコンプライアンス値（ばかりでなく腱と骨の互いに対する及び動脈に対する位置）を持つ異なった被検体は異なって反応し、これらの差はこのようなモデルでは償われない。更になお、単一被検体についてさえ、その被検体の組織及び動脈の反応の変化は、時間を通じて起こるか又は外部的に誘発されたストレスの関数として起こることがある。例えば、麻酔薬が被検体の循環系に導入されると、所与の動脈は実質的によりコンプライアンス性のものになり、それによってその反発性の多くを失うことになる。コンプライアンスのこの変化は、実動脈圧と測定された動脈圧との間の関係を変更し、したがって、それに基づくいずれの血圧推定値の正確性をも低下させる。
【００７９】
対照的に、本発明の方法論は、この重要な制限を、先行技術におけるように被検体を包括的にモデルするのに反して、各被検体についてその介在した組織と材料の実反応を測定することによって克服する。特に、本発明は、組織と動脈の実圧縮に基づいてこれらの構成要素のコンプライアンスの機能表示を発生する。
【００８０】
図３ｄに例示した方法３５０の模範的実施の形態では、前述した「Ａモード」音響伝送が圧平デバイスと動脈内壁との間に介在した構成要素の各々の圧縮を監視するために使用される。これらの構成要素の圧縮（ステップ３５２）は、一般に、それら構成要素の個々の材料性質に従って進行するが、これの材料性質は未知かつ相互依存性でありかつしたがってモデルするには甚だ複雑である。しかしながら、これらの構成要素の実圧縮の直接観察を行うことによって、外部的に印加された圧力と動脈壁と動脈内の加圧された流体（血液）との界面に印加された力との間に存在する伝達関数は、各個にかつ圧平される指定位置について近似的に決定することができる。図６に例示したように、動脈の内壁と動脈の上の皮膚の表面との間の領域は、皮膚６０２、組織６０４、及び動脈壁６０６のようないくつかの離散領域に副分割されるといえる。皮膚の表面６０８と皮膚／組織境界６１０との間の、皮膚／組織境界と組織／動脈境界６１２との間の、皮膚／動脈境界と動脈／血液境界６１４との間の、及び動脈／血液境界６１４と血液／動脈境界６１６との間の、それぞれ、距離ｄ１、ｄ２、ｄ３、及びｄ４は、Ａモード音響伝送を使用してステップ３５４で測定され、この音響伝送は、先に説明したように、これらの境界からの反射を識別する。更に、骨６２０と腱６１８との相対位置は、転送損失に大きな影響を及ぼす。事実、腱６１８と組織６２２との間、及び骨６２０と組織６２２との間にいろいろな種類の回復ばね力が存在する。圧力転送の損失は、これらの回復力を克服することばかりでなく、組織のコンプライアンスと少なくとも部分的に関連している。それゆえ、圧平の間中（及び心周期の特別部分の間中）、動脈直径（及び流体面積）と印加された外圧との間の伝達関数は、ステップ３５６によって開発される、特に、心周期の心拡張期部分と心収縮期部分に対しては、図７ａ〜７ｂに例示したように異なった伝達関数が存在する。（実動脈内圧に対して測定された）低印加圧力では、介在した組織、下に横たわる組織／腱、動脈壁、等の比較的小さな圧縮が起こっており、それゆえ、印加圧力の更に増大は、一般に、これらの構成要素の更なる圧縮に不釣合いに寄与する。印加圧力のより高い値では、介在した構成要素は実質的に圧縮され、かついくらかでも更に増大する圧平圧力の比較的小さい分数は介在したかつ下に横たわる構成要素の圧縮に使用される。それゆえ、一般に、「伝達関数」又は転送された圧力の印加された圧力に対する比は、図７ｂに例示したように、印加された圧力の関数として増大する。動脈内の加圧された血液とトランスジューサとの間に介在した自由浮動する非圧縮性材料という論理的場合には、伝達分数は１：１であることになり、完全結合を表示する。
【００８１】
それで、前述の導出された伝達関数は、圧力検知によって測定された不完全圧力転送に介在組織等によって導入された誤りを、ステップ３５８によって注目の領域を識別することによって、補正するために利用することができる。例えば、６０ｍｍＨｇ（８ｋＰａ）の圧力が測定されるとき心周期の心拡張期部分の間中に動脈内ゼロ経壁圧状態が達成されるならば、真心拡張期圧は或るパーセンテージより高いことになり、ここにパーセンテージは圧力転送損失の程度によって決定される。その監視する位置に対する伝達分数は、圧力測定センサの表面に転送される脈管圧（ｉｎｔｒａｖａｓｃｕｌａｒ  ｐｒｅｓｓｕｒｅ）の分数又はパーセンテージである。
【００８２】
伝達関数及び／又は伝達分数は、測定後に、図１２の装置に関して下に更に説明するように、ディジタル・ランダム・メモリ内のルックアップ・テーブルのような、いずれか種々の異なったフォーマットで表示されかつ記憶されることがあることに注意されたい。更に、直ちに承知のように、上に説明した方法３５０は１つ以上の離散圧力状態（すなわち、心収縮期圧及び／又は心拡張期圧）に対する伝達分数を決定するために使用されるが、伝達分数は圧力の範囲に対してたやすく決定することができ、それによって、下に更に詳細に説明するように、圧力の関数として伝達関数を形成する。それゆえ、被検体の血圧が変動するならば、本発明は、いずれかの圧力範囲内の圧力の測定された値を補正するためにこの伝達関数を利用する。
【００８３】
同様に、納得するように、Ａモード音響技術以外の伝達関数／伝達分数を決定する方法も、単独で又はＡモード技術と一緒に利用されることがある。
【００８４】
要するに、図３ｄの方法３５０は、被検体に対する印加ストレス（例えば、圧力）と監視される位置に関係したような伝達関数／伝達分数を決定し、かつパラメータの実値を決定するために、指定された「マーカ」点で測定されたパラメータを、伝達関数／伝達分数を使用して校正することを伴う。血圧監視の場合、このプロセスは、運動エネルギー項が増大し始める（又はその代わりに、最高血液速度が上昇し始める、流体面積が減少し始める、又は或る他の所望状態が観察される）に当たっての圧平圧力を印加し、次いで、実動脈内圧を決定するために、伝達分数を使用して、測定された圧力の測定値を補正することを伴う。全心周期にわたって考えるとき、この方法３５０は、全測定された圧力波形８００を真動脈内圧８０１へ校正するためにこの波形に適用される拡大縮小関数又は「伸縮（ｓｔｒｅｔｃｈｉｎｇ）」関数を生じ、それによって図８に示したような「校正された（ｃａｌｉｂｒａｔｅｄ）」波形８０６を生じる。注意するのは、考えられている測定された圧力波形の部分（例えば、心拡張期部分、心収縮期部分、又はそれらの間）に依存して実動脈内圧又はＡ線（Ａ−ｌｉｎｅ）動脈内圧の測定された圧力に対する比が変動することである。この概念は図８にグラフで例示してあり、そこでは心周期の心収縮期部分での振幅の比Ｒ１８０２は、心拡張期部分出の振幅の比Ｒ２８０４に等しくない。
【００８５】
注意を要するのは、或る種の環境及び個々の被検体が図３ｄに関して説明したような伝達関数の決定と適用を必要とするが、本発明の一般的方法論は或る場合には伝達関数がなくても潜在的に適用されることである。例えば、被検体が皮膚と動脈壁との間に高度の結合を有する場合、皮膚表面に置かれたトランスジューサを経由して測定された圧力に関連した誤りは、全測定の小さな分数を構成するに過ぎず、したがって、或る種の監視環境では許容可能であることになる。それゆえ、伝達関数の計算と適用は、あらゆる環境の下での本発明の要件ではない。
【００８６】
図３ｅを参照して、血行力学パラメータを連続的に校正する方法を説明する。上で図３ｄを参照して論じたように、伝達関数は、介在した組織、動脈壁等の圧縮に対して測定された圧力波形を補正するために有用である。この伝達関数は、橈骨動脈のような、被検体上の所与の監視位置で遂行される圧平掃引の間中に得られる。しかしながら、被検体の動脈血圧の連続監視をできるようにするためには、測定されたパラメータの変化を正確に監視しかつ校正事象（例えば、圧平掃引）と事象との間で拡大縮小することができる機構が必要である。
【００８７】
先に論じたように、先行技術の校正アプローチは、測定された圧力波形を「連続的に校正するために（アスクレトリ（ａｓｃｕｌａｔｏｒｙ）カフ測定のような）周期的校正事象に頼った。これらのシステムを参照して使用された用語「連続的に」は、いくぶん不適当な表し方である。というのは、実際に起こるのは拡大縮小関数の（連続的ではなくむしろ）周期的更新だからである。このアプローチは、少なくとも１つの深刻な欠陥を露呈し、それは周期的校正更新と更新との間隔中の校正の欠如である。監視される被検体の活動次第で、それらの真動脈血圧は、短期間に非常に顕著に、或る場合には急速に又はプロンプトする様相で変動することがある。例えば、外科手術の間中、動脈切除のような外科医による行為（ａｃｔｉｏｎ）は、被検体の動脈血圧を含む被検体の循環系に深い影響を及ぼすことがある。同様に、麻酔薬から結果する被検体の体内の動脈の多く内のコンプライアンスの変化に大部分起因して、誘発前（ｐｒｅ−ｉｎｄｕｃｔｉｏｎ）（すなわち、麻酔前）血圧値と誘発後血圧値との間の差は劇的に異なっている。
【００８８】
先行技術のアプローチは動脈内に起こる実血行力学性質を決して監視しないから、もし真動脈血圧のこのような重大な変化が周期的校正事象の間に起こるならば、それらは、多くの場合、検出されないままである。むしろ、このような先行技術のアプローチは、血圧を典型的に圧力測定法的に監視し、これらの測定は、潜在的に、実動脈圧と非常に異なっている。先行技術のシステムは、（真血圧に接近していることもあり、ないこともある）圧力測定法による血圧の測定された変化を償うために、典型的に拡大縮小率又は校正を調節する。この方法の結果は、事実は校正されないで、広く変動する拡大縮小成分を含む、いわゆる、「校正された」血圧値を生じることを含む。校正事象と事象との間の実圧つまり真動脈血圧を追跡することのこの失敗は、外科手術の間中のような、血圧の分刻み（ｍｉｎｕｔｅ−ｔｏ−ｍｉｎｕｔｅ）測定が決定的に重要であるといえる場合に壊滅的であるおそれがある。
【００８９】
本発明の方法論は、下に更に詳細に説明するように、第１、つまり、「主（ｐｒｉｍａｒｙ）」測定された血行力学パラメータ（例えば、血圧）の変化を追跡するために先に説明した測定された「第２」血行力学パラメータを使用することによって先行技術の前述の重大な限界を克服する。
【００９０】
１実施の形態では、図３ｅの方法３７０によって例示したように、血液の運動エネルギーが前述の音響（又は他の）技術を使用して監視される一方、ゼロ経壁圧状態（又は重要であると判定された或る他の状態）が動脈内で維持される。特に、図１２に関して下に説明する実施の形態ではまた圧力トランスジューサと超音波トランスジューサを含む圧平デバイスが所望圧力状態を維持するようにステップ３７２で監視される動脈の上の皮膚に対して「サーボ（ｓｅｒｖｏ）」され又は連続的に変調される。測定された（未校正の）主パラメータ、ここでは、圧力がステップ３７４によって同時に時間の関数として監視される。特定応用次第で、ステップ３７２の変調は、心周期の心拡張期部分の間中、又はこれに代わりその周期の心収縮期部分の間中指定値に経壁圧を維持するように制御されることがある。なお他の代替実施の形態では、圧平デバイスは、（１つ以上の完全な心周期にわたって計算された）平均経壁圧を所定値に維持するように変調され又はサーボされることがある。サーボすることは、所望最高血液速度状態、又は所望最大断面積状態を維持するためにまた実施されることがある。多くの他のこのような「ターゲット」サーボ値で置換しても等しく首尾良く行くことがあり、かつこのような値ばかりでなく（この値が基づいているパラメータ関係（例えば、それを動作させることが望ましい図５ｄの最高速度対流体面積プロット領域）がユーザのニーズ及び本方法を採用する特定応用によってもっぱら決定される。
【００９１】
次に、ステップ３７６で、第２血行力学パラメータが適切な技術を使用して時間の関数として測定される。本実施の形態では、全運動エネルギー（又は最高血液速度）は、先に説明した型式の音響ドップラー・システムを使用して測定される。
【００９２】
ステップ３７８で、ステップ３７６で測定された第２パラメータが、主パラメータの変化を識別するために分析される。例えば、圧平デバイスが心周期の心拡張期部分においてゼロ経壁圧を維持するようにサーボされるとき、運動エネルギーの変化は、動脈内血圧の変化を追跡するために使用される。この分析の結果は、圧平デバイスの更に調節が必要とされるかどうか決定するためにステップ３８０によって所定許容判定基準又は制御判定基準と比較される（ステップ３８２）。例えば、ステップ３７８〜３８２で全血流量運動エネルギーに顕著な増大又は増大の速度が観察され（それによって、動脈壁によって感じられた圧平圧力が真動脈圧を超えていた）としたならば、圧平圧力は、運動エネルギーのより僅かな増大又はその増大の速度の僅かな上昇によって反映される近ゼロ経壁圧状態を維持するように、減少させることができるであろう。納得するように、制御システム技術上周知の型式のファジイ論理又はＰＩＤコントローラを含む、１つ以上の他のパラメータの測定に基づいて１つのパラメータを制御する制御方式のいずれの型式も、所望挙動を成し遂げるために使用することができる。
【００９３】
前述にもかかわず、納得するように、図３ｅの方法を使用する第Ｉ血行力学パラメータの連続校正は、例えば、最高血液速度及び／又は動脈断面（流体）面積を含む追加の又は他の第２パラメータを使用して完遂されることがある。
【００９４】
また再び注意するのは、先行技術アプローチと対照的に、上に論じた図３〜３ｅの技術は、都合良く、監視されている被検体の循環系内のパラメータのモデリングも推定も伴わない。全ての情報は、選択された位置での被検体の直接測定を介して導出され、したがって、その個人とその位置に特に適していることである。
【００９５】
循環系の血行力学反応を特徴表示する方法
図９を参照して、生被検体の循環系の血行力学反応を特徴表示する方法を開示する。図９に例示したように、方法９００の第１ステップ９０２は、印加ストレスに関係して血管と関連した第１パラメータと第２パラメータとの間の第１機能関係を導出することを含む。動脈血圧測定の文脈の中で、ステップ９０２で導出された第１機能関係は、本明細書に先に説明したように、動脈断面積（圧平圧力）と全血流量運動エネルギーとの間の関係を含む。もっとも、納得するように、多くのいろいろな異なった機能関係でそれを置換することがある。例えば、最高血液速度と流体面積との間、又は容積血流量（Ｑ）と流体面積との間の機能関係を、所望ならば、使用してよい。
【００９６】
次に、ステップ９０４で、上のステップ９０２で導出された機能関係中に現れる１つ以上のアーチファクト又はマーカが識別される。先に説明した動脈血圧測定の場合、アーチファクトは、ゼロ経壁圧の状態が心拡張期状態及び／又は心収縮期状態に対して達成された後増大する運動エネルギー又は上昇する血液速度を含む。これらのアーチアクトは、本明細書に図３ｃ〜３ｅに関して先に説明した校正関数に対する点を含む。これに代わり、Ａモード超音波技術又は他の類似の技術によって測定される際に、壁直径が心周期の心収縮期部分又は心拡張期部分で減少し始める点は、ゼロ経壁圧のマーカを構成することがある。
【００９７】
次に、ステップ９０６で、上のステップ９０２からの機能的に関係したパラメータの１つが印加されたストレスの関数として非侵襲的に測定される。上で参照した例では、この測定は、時間（かつ圧平圧力）の関数として動脈内の血液速度を測定し、かつそれから全流量運動エネルギーを導出することを含むことになる。
【００９８】
最後に、ステップ９０８で、校正「関数」（これは理論上１つほどの数少ないデータ点であることがある）が循環系と関連した選択されたパラメータの測定された反応にステップ９０４で識別されたアーチファクトに基づいて適用され、それによって、そのパラメータの反応の校正された特性表示を生じる。血圧については、選択されたパラメータは圧力測定法で測定された（すなわち、未校正の）圧力を含み、及び校正された特性表示は、なかでも、血液の運動エネルギーが増大し始める点で決定された校正された（又は「真の」）動脈血圧を含む。
【００９９】
更になお、（組織コンプライアンスに原因がある不完全信号転送に因る誤りのような）潜在的誤りの影響が、ステップ９０８の部分として、また充分に償われることがある。
【０１００】
呼吸を含む、周期的誤り源に対する校正方法
図１０を参照して、周期的誤り源に対して血行力学測定を校正する方法を開示する。方法１０００の第１ステップ１００２は、血管と関連した第１血行力学パラメータを測定することを含む。先に説明したように、このパラメータは、動脈血圧、又は微分された圧力のような他の圧力を含むことがある。動脈血圧測定の場合、このパラメータは、圧力測定法の圧力トランスジューサを使用して測定された未校正の圧力波形である。
【０１０１】
次に、ステップ１００４で、第２血行力学パラメータが、先に説明したように、被検体について測定される。この第２血行力学パラメータは、運動エネルギー、最高血液速度、最大動脈直径、最大流体面積等を含むことがある。１実施の形態では、運動エネルギーは、ドップラー超音波を使用して行われる血液速度の測定に基づいて計算される。
【０１０２】
次に、ステップ１００６で、第１パラメータと関連した周期的誤り源が第２パラメータ中で識別される。１模範的場合では、図１１に例示した、周期的誤り源は、監視されている被検体の呼吸に関係する。図１１に示したように、人間の橈骨動脈内を流れる血液の速度と運動エネルギーは、時間と共に変動する周期的成分を含む。この周期的挙動は、被検体の呼吸周期に実質的部分を起因しており、かつ典型的心周期よりも遥かに低い周波数で起こる。それゆえ、正規心周期１１０２は、周期的呼吸変動１１０４によって「振幅変調」されると考えることができる。
【０１０３】
呼吸周期的変動の起源は、横隔膜が上下するに連れて起こる変動する圧力に関係する。吸息と共に、横隔膜は当然下がり、腹腔内圧（ｉｎｔｒａ−ａｂｄｏｍｉｎａｌ  ｐｒｅｓｓｕｒｅ）を上昇させ、かつ胸腔内圧（ｉｎｔｒａ−ｔｈｏｌａｃｉｃ  ｐｒｅｓｓｕｒｅ）を低下させる。腹腔から右心房への圧力差の増大は、右心房へ帰る容積流量を増大させる。呼気と共に、横隔膜は上がり、腹腔内圧は低下し、かつ胸腔内圧は上昇する。結果は、下肢から腹腔へのより多くの静脈環流であるが、しかし右心房への環流は少なくなる。容積と圧力の巡回変化は循環系全体を通してあらゆる所に反映する。というのは、それが閉じた系であるからである。
【０１０４】
上述の巡回呼吸変化は、なかでも、心拡張期圧と心収縮期圧に対して変動流速と変動運動エネルギーを生じる。正常成人では、本願の出願人によって得られた逸話的証拠は、このような変動の大きさが酷しい場合で２０ｍｍＨｇ（２．７ｋＰａ）以上の程度であると示唆する。心収縮期圧と心拡張期圧の分数又はパーセンテージを取り挙げると、圧力のこの変動は、被検体が外科手術の間中のような特に歩行困難であるときに測定された低目の心拡張期圧に対して重大である。
【０１０５】
これらの変動は、本発明において、必要とされるとき、第２パラメータ（例えば、速度、運動エネギー、又は面積）の測定からの校正関数の導出に同期させることによって償われる。特に、方法１０００のステップ１００８で、呼吸変動の周期性が分析され、決定され、かつこの情報が校正関数の導出を周期（「搬送」）呼吸波形上の共通点に同期させるために使用される。呼吸成分とその周期性の識別は、信号処理技術上周知の多くのアルゴリズムのいずれか１つを使用して完遂される。したがって、このようなアルゴリズムは、本明細書に更に説明しないことにする。注意するのは、被検体の呼吸速度及び／又は呼吸の「深さ」は時間と共に変動するということができ、それによって動脈内の圧力変動／流量変動の周期性と大きさに影響するから、呼吸作用の周期性を連続的に（又は少なくとも頻繁に）計算するべきことである。
【０１０６】
次に、ステップ１０１０で、校正関数がステップ１００８の結果によって規定された周期性で行われた第２パラメータの測定に基づいて開発される。例えば、一連の血液速度測定は７秒目毎に行われ（各測定は呼吸波形上の同じ相対点であるが、しかし時間的に変位した点に対応する）、かつこの情報は図３〜３ｅを参照して本明細書に先に説明した運動エネルギー値及び校正関数又は「伸縮」関数を導出するために使用される。
【０１０７】
最後に、ステップ１０１２で、ステップ１０１０の伸縮関数がステップ１００２の測定された（すなわち、未校正の）波形に適用される。呼吸波形に対する類似の点で第２パラメータを測定することの効用で、全呼吸周期を横断しての呼吸の作用は償われることに注意されたい。それゆえ、導出された伸縮関数は、第２パラメータが実際に測定された時点に対応する波形の部分にのみ適用されるのとは対立的に、全未校正の圧力波形）に適用されるといえることに注意されたい。圧力転送が呼吸に伴う心収縮期圧変動の回りで、及び呼吸に伴う心拡張期圧変動の回りで比較的線形であると想定すると、他の補正は必要ないことになる。非線形性が充分に顕著ならば、呼吸周期の異なった位相で補正率（ｃｏｒｒｅｃｔｉｏｎ  ｆａｃｔｏｒ）を計算することによって、追加補正を適用することができる。これは、真動脈血圧の（周期的によるのとは対立的に）連続的に校正された表示を用意することの重要な利点を表す。
【０１０８】
納得するように、前述の論議は呼吸系作用に因る周期的誤りの面で行われているが、周期的又は非周期的な、他の型式の誤りも、図１０に例示したような本発明の方法論を使用して、償うことができる。例えば、被検体の心臓内不整脈（ａｒｒｈｙｔｈｍｉａ）の作用は、校正関数の導出の間中に識別しかつ償うことができる。被検体の心臓内不整脈は、波形内の非周期的成分を識別する目的に特に適合した信号処理アルゴリズムを使用して識別することができ、このようなアルゴリズムは信号処理技術の当業者に周知である。非周期成分誤りの多くの他の型式もまた図１０の方法と関連して識別することができる。
【０１０９】
血行力学査定装置
図１２を参照して、生被検体の血管内血行力学性質を測定する装置を説明する。例示した実施の形態では、その装置は、人間の橈骨動脈内血圧の測定に適合している。もっとも、納得するように、他の血行力学パラメータ、他の監視サイト、及び生器官の他の型式も、本発明の最も広い意味で本発明と一緒に利用されることがある。図１２の装置１２００は、基本的に、圧力測定法によって橈骨動脈から血圧を測定する圧力トランスジューサ１２０２、その動脈上の圧平（圧縮）の程度を変動させるためにトランスジューサ１２０２に結合された圧平デバイス１２０４、音響発射が血液速度（と運動エネルギー）を作り出すために使用され、これらの音響発射とそれの反射を発生する音響トランスジューサ１２０６、圧力トランスジューサ１２０２、音響トランスジューサ１２０６によって発生された信号を分析しかつそれらの信号に基づいて校正関数を発生するためにこれらのトランスジューサに動作上接続された信号プロセッサ１２０８、動脈内への伝送用音響信号を発生しかつ音響トランスジューサ１２０６から信号を受信する信号発生器／受信器１２１０、及び動脈に印加される圧平圧力の程度を制御するために圧平デバイス１２０４と信号プロセッサ１２０８に動作上接続されたコントローラ１２１１を含む。
【０１１０】
圧力トランスジューサ１２０２は、本実施の形態では、圧電トランスジューサ素子であって、その検知面１２１２に印加された圧力に関数関係にある（例えば、比例した）電気信号を発生する。同様に、音響トランスジューサ１２０６は、圧電（セラミック）デバイスを含み、このデバイスはモードに依存して音響波及び／又はパルスを発生しかつこれらを受信する能力がある。例示の実施の形態では、音響トランスジューサ１２０６は、８ＭＨｚに中心を置いた超音波周波数を発生するように同調させられる。もっとも、変動帯域幅を持つ他の中心周波数を使用してもよい。信号発生器／受信器１２１０は電気信号又は電気パルスを発生し、これらは音響トランスジューサ１２０６に供給されかつ血管内へ放射される音響エネルギーに変換される。この音響エネルギーは、動脈内を流れる血液を含む動脈内の種々の構造物ばかりでなく、動脈の付近の組織及び他の肉体上の構成要素によって反射される。これらの音響反射（反響）は、音響トランスジューサ１２０６によって受信されかつ電気信号に変換され、電気信号は、次いで、信号発生器／受信器１２１０によって（例えば、ＡＤＣを使用して）ディジル形式に変換され、分析のために信号プロセッサ１２０８に送られる。採用された音響分析の型式とモードに依存して、信号プロセッサ１２０８は、受信された信号を分析するために（埋め込まれた又は外部記憶デバイスに記憶されたのどちらかの）そのプログラムを利用する。例えば、そのシステムが最高血液速度を測定するために使用されるならば、受信された反響は、なかでも、ドップラー周波数偏移について分析される。これに代わり、動脈直径（面積）が測定されるならば、前述のＡモードに適当な分析が採用される。
【０１１１】
校正「掃引」の間中、コントローラ１２１１は、所定プロファイルに従って動脈（と介在した組織）を圧平するように圧平デバイスを制御する。この掃引の間中、音響信号は、組織の進行中の圧平に直接近接する領域内へ、好適には、動脈内へ伝送されかつこれから受信される。速度データ、運動エネルギー・データ、及び／又は動脈直径データは、受信された反響から抽出及び／又は導出され、かつ心周期の選択された部分に対する圧平圧力の関数として記録される。信号プロセッサ１２０８及び関連したアルゴリズムは、１つ以上のマーカを識別し、かつ連続監視が測定されたマーカに基づいて起こるに当たっての所望印加圧力を決定する。例えば、図４ｂに示した最高速度中のピークがマーカとして選択されたとしたならば、アルゴリズムは、このピークを識別しかつこのピークに相当する圧力データを識別することになるであろう。続く血圧監視の間中、コントローラ１２１１は、ターゲット圧力又はプログラマ／ユーザによって選択されたいずれかの値を維持するように圧平デバイス１２０４（本実施の形態では、圧力トランスジューサ１２０２）の位置をサーボすることになるであろう。測定されるパラメータ（例えば、全血液流運動エネルギー）のその後の変化は動脈内の実血圧の変化を識別するために使用され、それによって一連の校正掃引を引き続き行う必要を回避する。
【０１１２】
オプショナルに、装置１２００は、信号源とセンサとの間に介在した組織及び他の肉体的構成要素の伝達関数を測定するようにまた構成される。上に図７に関して説明したように、動脈内の血液と圧力センサとの間にエネルギーの不完全転送又は分数転送がある。この問題点（ｉｓｓｕｅ）に取り組むために、図１２の装置１２００は、圧平されたときの動脈直径とその付近の肉体的構成要素の相対圧縮に関するＡモード分析又は他の技術から得られたデータを利用する伝達関数アルゴリズム（図示してない）を含む。それゆえ、校正掃引の間中、装置１２００は、変動する圧平圧力の関数として動脈と構成要素の相対圧縮を決定する伝達関数アルゴリズムによって使用されるＡモード・データ又は他の匹敵するデータを記憶する。伝達関数（例えば、圧平圧力の関数としての動脈直径の変化）は、そのアルゴリズムによって発生され、かつルックアップ・テーブル又は数学関数のような、多くの異なった有様で記憶される。校正掃引に続いて、装置１２００が識別されたマーカから導出された所望印加圧力にサーボするに連れて、補正が測定された圧力に伝達関数に基づいて課せられる。例えば、システムが圧力トランスジューサ１２０２によって測定された６０ｍｍＨｇ（１９．８ｋＰａ）の心拡張期圧へサーボしつつあるならば、動脈内の圧力の真値は、伝達関数に従って６０ｍｍＨｇ（１９．８ｋＰａ）よりもいくらか高い値へ補正される。
【０１１３】
図１３を参照して、図１２の血行力学パラメータを査定する装置の１特別実施の形態を説明する。図１３の実施の形態では、装置１３００は内蔵ユニットを含み、このユニットは組合わせ圧力トランスジューサ１３０２と圧平デバイス１３０４、音響トランスジューサ１３０６、（ファジイ論理アルゴリズムのような）制御マイクロ・コードを備えたマイクロコントローラ１３０８、埋込みメモリ１３１２と命令集合（校正ルックアップ・テーブルを含む）を備えたディジタル信号プロセッサ（ＤＳＰ）１３１０、信号発生器と受信器ユニット１３１４、記憶デバイス１３１６、ディスプレイ・デバイス１３１８、及び電源１３２０を有する。この実施の形態では、マイクロコントローラ１３０８は、初期圧平「掃引」が遂行されるように、組合わせ圧力トランスジューサ１３０２／圧平デバイス１３０４の動作を制御するために使用される。特に、圧力トランスジューサ１３０２は、橈骨動脈をおおって、被検体の手頸の内側皮膚と連絡して置かれ、かつ図１３に例示したようにあるべき場所に固定される。橈骨動脈から未校正の血圧測定が開始され、その直ぐ後でマイクロコントローラ１３０８は圧平機構１３０４に指令してトランスジューサ１３０４を被検体の手頸に対して押し圧すると共に圧力を増させ、それによってその下に横たわる動脈を圧平させる。動脈が圧平されるに連れて、音響トランスジューサ１３０６は、動脈をおおう皮膚と連絡するようにまた押し付けられ、かつ信号発生器１３１４は一連の音響パルスを発生し、これらのパルスは皮膚を通して動脈内へ伝送される。動脈の圧平が続くに連れて、信号発生器／受信器ユニット１３１４は、動脈内の血液及び他の構成要素から音響トランスジューサ１３０６を経由して反響を受信し、かつ受信された反響に関係する出力信号を発生する。この出力信号は、処理され、かつ次いで、ＤＳＰ又は類似の処理エンジンによるその後の分析のためにディジタル化される。同様に、圧力トランスジューサ１３０２の出力信号は、ディジタル化され、かつＤＳＰに入力される。ディジタル化された信号は、次いでＤＳＰ内の埋込みプログラムを使用して分析され、このプログラムは本明細書に後に説明するコンピュータ・プログラムの機械コード表現である。ディジタル信号プロセッサの出力は、補正された圧力波形であって、これが、次いで、ユーザに表示されるためにディスプレイ・デバイス１３１８に（使用されたデバイスの型式次第でディジタル形式でかアナログ形式でかいずれにせよ）供給される。オプショナルに、ＤＳＰの出力は、記憶デバイス１３１６内の１つ以上の記憶位置に記憶される及び／又は外部デバイスへ出力される。
【０１１４】
注意するのは、本明細書に説明した装置１２００、１３００は、種々の異なった構成要素を使用して、種々の構成に組み立てられることがあり、かつ種々の異なった血行力学パラメータを測定することがあることである。模範的な制御、信号発生／処理、及び圧平の機構及び回路は、先に本明細書にその内容を組み入れてある、讓受人の「動脈血圧の非侵襲的決定方法と装置（Ｍｅｔｈｏｄ  Ａｎｄ  Ａｐｐｒａｔｕｓ  Ｆｏｒ  Ｔｈｅ  Ｎｏｎｉｎｖａｓｉｖｅ  Ｄｅｔｅｒｍｉｎａｔｉｏｎ  Ｏｆ  Ａｒｔｅｒｉａｌ  Ｂｌｏｏｄ  Ｐｒｅｓｓｕｒｅ）」と題する同時係属米国特許出願第０９／３４２，５４９号に説明されている。
【０１１５】
コンピュータ・プログラムと関連装置
血行力学査定、モデリング、及び校正の前述した方法を実施するコンピュータ・プログラムを説明する。１模範的実施の形態では、そのコンピュータ・プログラムは、図３〜３ｅ、９、及び１０の方法論を、個々に又はそれらの組合わせのどちらかで実施するＣ＋＋ソース・コード・リスティングのオブジェクト（「機械」）コード表現を含む。Ｃ＋＋言語が本実施の形態に使用されるが、承知のように、例えば、ＶｉｓｕａｌＢａｓｉｃＴＭ、Ｆｏｒｔｒａｎ、及びＣ＋を含む、他のプログラミング言語も使用されることもある。ソース・コード・リスティングのオブジェクト・コード表現は、図１４ａ〜ｂに例示したように、コンパイルされかつコンピュータ技術上周知の型式の媒体記憶デバイス上に配置される。このような媒体記憶デバイスには、光ディスクやＣＤ・ＲＯＭや磁気フロッピー（Ｒ）・ディスクや「ハード」ドライブや磁気バブル・メモリさえもが、これらに限定することなく、ある。コンピュータ・プログラムは、プログラミング技術上周知の型式のグラフィカル・ユーザ・インタフェース（ＧＵＩ）を更に含み、これは（図１５に関して下に説明するように）ホスト・コンピュータ又はプログラムがその上でランされる装置のディスプレイと入力デバイスに動作上結合される。
【０１１６】
全体構造の面で、プログラムは、１実施の形態では、ホスト・コンピュータに供給された測定されたパラメータに基づいて本明細書に説明した血行力学査定、モデリング、及び校正方法論を実施する一連のサブルーチン又はアルゴリズムで構成される。第２実施の形態では、コンピュータ・プログラムは、図１２又は１３の前述した血行力学測定装置と関連したディジタル信号プロセッサ（ＤＳＰ）又はマイクロプロセッサの、埋込み記憶デバイス、すなわち、プログラム・メモリ内に配置されたアセンブリ言語／マイクロコード化命令集合を含む。
【０１１７】
本明細書に開示したようにパラメータ・データを分析しかつ血行力学パラメータの校正された値を発生する能力のある装置の１実施の形態を説明する。コンピューティング・デバイス１５００はマザーボード１５０１を含み、このマザーボードは中央処理ユニット（ＣＰＵ）１５０２、ランダム・アクセス・メモリ（ＲＡＭ）１５０４、及び（直接メモリ・アクセス・コントローラのような）メモリ・コントローラ１５０５を有する。（ハード・ディスク・ドライブ又はＣＤ−ＲＯＭのような）記憶デバイス１５０６、（キーボード又はマウスのような）入力デバイス１５０７、及び（ＣＲＴ、プラズマ、ＴＦＴディスプレイのような）ディスプレイ・デバイス１５０８ばかりでなく、ホスト構成要素と周辺構成要素の動作を支援するために必要なバスがまた備わる。直列又は並列入出力ポート１５１１が、装置１５００へのかつこれからのデータ及び／又は制御信号の転送のためにまた含まれる。
【０１１８】
血行力学パラメータを査定するために有効な上述のコンピュータ・プログラムは、パラメータ査定の間中にＣＰＵ１５０２によって使用されるためにＲＡＭ１５０４及び／又は記憶デバイス１５０６に機械読取り可能オブジェクト・コード表現の形で記憶される。ユーザ（図示してない）は、システム・オペレーションの間中にプログラム・ディスプレイと入力デバイス１５０７を経由してコンピュータ・プログラムと関連測定機器用１つ以上の機能モードを選択することによって注目の血行力学パラメータを査定する。特に、動脈血圧測定の場合、ユーザは、被検体の選択された血管上に必要なパラメータ・センサを置き、かつ連続的又は所定の間隔でのいずれかでそれらのセンサによって出力されたデータを受け取るようにコンピュータ・プログラムを構成する。コンピュータ・プログラムは、信号が装置１５００に供給されるならば、先に説明した分析を遂行し、かつ校正された信号をディスプレイ・デバイス上に、又はシステム自体のディスプレイ上に表示されるように発生する。ルックアップ・テーブル又は類似の機構が、図１２を参照して先に説明したように、校正を容易にするためにコンピュータ・メモリ・デバイス又は記憶デバイス内に記憶される。プログラムによって発生されたこのような校正された測定は、後の検索のために記憶デバイス１５０６にオプショナルにまた記憶されるか、又は所望ならば、直列又は並列ポート１５１２を経由して、プリンタ、データ記憶ユニット、その他の周辺構成要素のような外部デバイスへ出力される。更になお、装置１５００は、ネットワーク・インタフェース・カード（ｎｅｔｗｏｒｋ  ｉｎｔｅｒｆａｃｅ  ｃａｒｄ；  ＮＩＣ）又は類似のインタフェース（図示してない）を経由して他のコンピューティング・デバイス又はデータベースにネットワークされることがあり、このことによって装置１６００によって発生されたデータは遠隔で分析又は記憶されることがある。このような遠隔デバイスへの伝送は、ローカル・エリア・ネットワーク（ＬＡＮ）や、イントラネットや、インターネットや、光ファイバ・システムや、ラジオ周波数（無線）デバイスのような種々の周知の方法を使用して完遂されるといえる。
【０１１９】
なお他の実施の形態では、装置は、（パーソナル・ディジタル・アシスタント、つまりＰＤＡのような）パーソナル・コンピューティング・デバイスを含み、このデバイスは圧力センサと音響センサから入力データを受信しかつそのデータを分析して血圧の校正された測定を生じるために適合している。また納得するように、ラップ・トップコンピュータ、計算器、及びパーソナル・オーガナイザのような他のポータブル・デバイスも本発明のコンピュータ・プログラムをランさせるように構成することもできる。更になお、入力センサ・データをこれらのデバイスに伝送する種々の異なった方法を使用することもでき、これにはネットワークされたコンピュータ、又は無線データ・リンクがある。
【０１２０】
治療を施す方法
図１６を参照して、血行力学パラメータを査定する前述の方法を使用して被検体に治療を施す方法を開示する。図１６に例示したように、方法１６００の第１ステップ１６０２は未校正の血行力学パラメータを非侵襲的に監視することを含む。血圧の場合、橈骨動脈に適用される模範的圧力トランスジューサは、本明細書に図３ａに関して説明したように使用される。
【０１２１】
次に、図１６のステップ１６０４で、ストレスが血管上に誘発され、血管はその血行力学性質を（少なくとも局所的に）変更し、それによって、全体的に血管又は循環系と関連した他のパラメータに変化を誘発する。本明細書に図３ｂに関して先に論じたように、このストレスは、１実施の形態では、血管を圧平する、つまり時間の関数として可変的に圧縮することを含み、それによってなかでも、圧平の領域で血液の容積流量（Ｑ）、速度（ｖ）、及び運動エネルギー（ＫＥ）に変化を誘発する。しかしながら、再び注意するのは、思うに他のストレッサも印加され、これらが同様の又は他の血行力学性質に影響することがあることである。
【０１２２】
次に、ステップ１６０６で、血管と関連した第２パラメータが下のステップ１６０８での校正関数の導出を容易にするために測定される。本明細書に図３ｃに関して論じたように、第２パラメータは、１実施の形態では、全血流量運動エネルギー又は最高血液速度を含む。というのは、これらのパラメータは、ステップ１６０４でのストレッサの印加の関数として或る種の容易に識別される「アーチファクト」を示すからである。しかしながら、同じ又は他のアーチファクトを示す他のパラメータも、校正関数を導出するために使用されることがある。
【０１２３】
図１６のステップ１６０８で、校正計量又は校正関数がステップ１６０６で導出されたパラメータ情報に基づいて次に導出される。特に、１つ以上のアーチファクト又はマーカがパラメータ・データ中で識別され、これらのアーチファクトは上のステップ１６０２の第１パラメータの実値と測定された値との間に或る種の関係が存在するときを表示する。本明細書に図３〜３ｅを参照して論じたように、校正関数を導出するプロセスの１実施の形態は、（ｉ）印加されたストレッサ（圧平）のエリアに近接する運動エネルギーの又は最高血液速度のプロファイルを測定し、かつ（血管の断面積のパーセンテージ減少に相関した）圧平の関数としてこれらのプロファイル中の速度を上昇させる又は運動エネルギーを増大させる領域を識別すること、及び（ｉｉ）圧力測定の領域内の組織と他の肉体的構成要素について伝達関数を測定することを含む。
【０１２４】
図１６の方法のステップ１６１０で、ステップ１６０８で導出された校正関数は、補正された又は校正された測定を発生するためにステップ１６０２の第１パラメータの測定に適用される。第１パラメータが時間の関数として連続的に（又は周期的に）測定されるならば、ステップ１６０８の補正関数が適当であるに従って連続的に又は周期的に適用されるか、又はこれに代わり第２血行力学パラメータが校正関数の変化を表示するために（先に説明した圧力サーボする間中になど）監視されることに注意されたい。
【０１２５】
最後に、ステップ１６１２で、第１パラメータの校正された測定が被検体に治療を施す基礎として使用される。例えば、血圧測定の場合、校正された心収縮期血圧値と心拡張期血圧値が発生されかつ表示されるか、そうでなければ、外科手術の間中のような実時間に、保険管理プロバイダー（ｈｅａｌｔｈ  ｃａｒｅ  ｐｒｏｖｉｄｅｒ）に提供される。これに代わり、このような校正された測定が延長期間にわたって収集されかつ被検体の循環系の状態又は反応における長期傾向について分析されることがある。
【０１２６】
代替（後方散乱）実施の形態
先に論じたように、（「Ａモード」超音波エネルギーのような）後方散乱音響エネルギーは、例えば、動脈直径と終末心拡張期血圧を含む血行力学パラメータを決定するために使用されることがある。本明細書に使用されるように、用語「Ａモード」は、通例のＡモード（すなわち、繰り返し線を経由する振幅対深さ（時間）表示）及びＭモード（すなわち、時間次元にわたって分布したグレースケール画像）とＢモード（すなわち、深さ（時間）対幅を表示する「ステアード（ｓｔｅｅｒｅｄ）」グレースケール画像）のような他の関連モード性の両方を含み、このような代替モード性は超音波技術における当業者に周知である。それゆえ、納得するように、次の議論はＡモードの面で行われるが、他の関連モード性もまた充分に採用されることがある。
【０１２７】
図１７〜３７は、本発明に従って血行力学パラメータを測定する方法と装置のなお他の実施の形態を例示する。この実施の形態は、Ａモード超音波エネルギーを使用して、その結果生じる反射されたエネルギーの波形の分析に基づいて動脈壁位置を決定し、壁位置は検知した圧力波形中の「マーカ」又はアーチファクトの識別を援助するために使用され、マーカ又はアーチファクトは血圧測定に相関している。本実施の形態のＡモード超音波エネルギーは、都合良く、本明細書に先に説明したドップラー・スペクトル・モード性を相補するが、所望ならば単独でもまた使用されてよく、それによってシステム設計者と装置の最終エンド・ユーザの両方に設計と操作の両面で大いに強化された柔軟性を与える。特に、本発明によって少なくとも３つの異なった動作方式が企図され、それは、（ｉ）Ａモード単独、（ｉｉ）ドップラー・ベースト検出単独、及び（ｉｉｉ）Ａモードとドップラー・モード性を一緒にした組合わせ、を含む。必要とするのではないが、（ｉｉｉ）におけるような相補モード性は、器械で以て得られた結果の信頼性と堅牢性を、都合良く、強化する。
【０１２８】
血管がますます圧縮される（圧平）に従って、種々のこのような血行力学的マーカが検知された圧力波形中に起こることが本願の讓受人によって観察されている。１態様では、波形の終末心拡張期部分中に或る種の「フラット・スポット」が起こり、そこでは時間の関数としての検知された圧力の変化が低い（又は、少なくとも、波形の周囲部分より低い）。これらのフラット・スポットは、血管を取り囲む組織及び血管自体内の種々の現象に相関し、これらの現象は血管の壁と関連したフープ（ｈｏｏｐ）ストレスを克服することを含む。本明細書に先に論じたように、人間の橈骨動脈のような実質的に円筒中空体と関連したフープ・ストレスは、更なる圧縮に対して逆力（ｒｅｔａｒｄｉｎｇ  ｆｏｒｃｅ）として作用する。換言すると、有限力（つまり所与の面積にわたる圧力）を血管の断面積が顕著に変形する前に血管の外面に印加しなければならない。変形のこの点を超えると、血管を更に変形させるためには比例する小さな追加力が必要とされる。それゆえ、いったん克服されると、フープ・ストレス成分は有効に除去され、それによって、外部的に血管壁に印加された圧力が血管内を流れる加圧された管腔（例えば、血液）と平衡することができて、血管の断面を本質的に一定に（ゼロ経壁圧に）維持する。この状態は、いったんフープ・ストレスが克服されると、血管内圧が、事実、動脈壁を横断して直接伝送されるという点で重要である。しかしながら、このような平坦化（すなわち、ゼロ経壁圧）状態がその上に存在する圧平プロファイルに沿う点は、個から個で、種から種で、及び単独の個については位置から位置で、変動することになる。したがって、所望経壁圧状態が存在する圧力波形上の点を正確に孤立させるには、追加の技術が必要とされる。
【０１２９】
（下に更に詳細に説明する）１つのこのような技術は、評価される血管の前壁と後壁を検出することであり、この情報は血管の直径の減少を直接決定するために使用される。先に示したように、血管直径の減少は或る種の心収縮期条件と心拡張期条件に直接関係させることができ、かつ前述の血行力学マーカの或る種の生起と相関させることができる。
【０１３０】
管腔と壁検出
本発明の方法論は、都合良く、血管の管腔、及び前壁と後壁の位置を検出するためにＡモード超音波エネルギーを採用する。本明細書に使用されるように、用語「前（ｆｒｏｎｔ）」と「後（ｂａｃｋ）」は、特別な空間的配向（ｏｒｉｅｎｔａｔｉｏｎ）の意味を含むのではなく、超音波源から発射された移動する超音波の波頭が直接伝搬の間中に血管の壁と出合う順序を単に称する。それゆえ、１つのトランスジューサの配向で血行力学性質を測定するとき「前」壁であるといえるものが他のトランスジューサの配向で「後」壁になることがある。本発明に特有なのは、前壁と後壁の検出の手段としての管腔検出のまた使用である。すなわち、他の血管構成要素に比較して管腔による超音波エネルギーの比較的低い後方散乱に因り、管腔の位置は後方散乱のエネルギーから容易に識別されることがあり、かつ前血管壁と後血管壁はこれに対して識別される。
【０１３１】
図１７は、本発明に従って管腔検出と壁検出に後方散乱超音波エネルギーを使用する一般的方法論を例示する。超音波エネルギーは装置（例えば、図１２のシステム１２００の超音波トランスジューサ１２０６）によって発生され、かつ被検体の皮膚を通してその下に横たわる組織内へ直接伝送され、かつ最終的に注目の血管を通る（ステップ１７０２）。本文脈の中で、用語「組織」は、検査される血管の付近に存在する全ての構成要素、いずれもの介在した皮膚、筋肉、骨、静脈ばかりでなく、血管（例えば、動脈）壁及び管腔自体を含む意味に取られる。次に、信号処理がその後の管腔検出と壁検出を援助するために受信された後方散乱エネルギーに遂行される（ステップ１７０４）。
【０１３２】
次いで、処理された信号が、その中で管腔の位置を検出するためにステップ１７０６で分析される。特に、処理された後方散乱エネルギー信号は、評価される種／年齢範囲での考察下の血管の深さと両立する深さの範囲内での後方散乱エネルギーの「最小」の生起を識別するために分析される。この最小の位置は、評価される血管内に存在する管腔と一致すると想定される。
【０１３３】
次に、１つ又は両血管壁が分析の開始点としてステップ１７０６で検出された管腔を使用してステップ１７０８で検出される。前血管壁と後血管壁と関連した信号は、動脈内管腔の識別された位置から前方へかつ後向へ探索することによって検出される。後方向での（トランスジューサへ向けての）探索と関連した後方散乱エネルギーが（下に詳細に説明する）規定された量又は率だけ平均管腔パワー又は血液信号の推定されたレベルを超えるというように、１つ以上の判定基準を満足させるとき、前壁が検出され、かつ線に沿う信号の位置が前壁の位置として解釈される。同様に、前への探索と関連した反射エネルギーが関連後壁判定基準を満足させるとき、後壁が検出されかつ線に沿う信号の位置が後壁の位置として解釈される。壁検出のための種々の方法論を図２９ａと２９ｂを参照して下に説明する。
【０１３４】
図１７ａに示したように、図のＡモード管腔及び壁検出アプローチと関連した信号処理（ステップ１７０４）の１実施の形態は、「生（ｒａｗ）」後方散乱Ａモード信号をベースバンド化（ｂｓｓｅｂａｎｄｉｎｇ）（直角復調）し（ステップ１７２０）、ベースバンド化されたＡモード信号を平方し（ステップ１７２２）、それによって「包絡線平方」データを生じ、かつステップ１７２４によって所定サンプリング・レート（例えば、２ＭＨｚ）にデシメートする。次に、深さ依存利得（例えば、ＴＧＣ）がステップ１７２６によって包絡線平方データに適用されて、伝搬損失に対して調節する。これらの信号処理ステップは、次のパラグラフで更に詳細に説明する。
【０１３５】
図１８は、本明細書で「生（ｒａｗ）」Ａモード・データと称する生（ｌｉｖｉｎｇ）被検体の橈骨動脈と関連したＡモード信号振幅対時間のプロットである。注意するのは、図１８の時間軸は組織を通しての超音波の伝搬速度に基づいて深さに変換されていることである。この変換が遂行されるのは、深さは本発明の目的のために時間よりも容易に解釈されるからである。もっとも、納得するように、時間（又は時間スケールのような他の表現）も等しく使用されて巧く行く。
【０１３６】
図１８に明示した１つの重要な特徴は、深さの関数としてのＡモード信号のレベル又は振幅である。図１８に例示したように、この例では約１．５ｍｍから３．３５ｍｍの明白な範囲があり、そこでは信号レベルが極めて低い。この低い信号振幅は、検査下の血管内を流れる血液（管腔）と関連したエネルギーの低後方散乱に由来する。
【０１３７】
更に、この例に使用された１６ＭＨｚ中心周波数は、図１８に明白に描かれている。中心周波数は動脈壁の位置に関連する情報を輸送しないから、信号処理の第１ステップは、中心周波数を除去することである。このプロセスは、図１７ａにステップ１７２０として識別された、直角復調（又は包絡線検波）として表される。図１８の模範的Ａモード信号と関連した包絡線検出された信号を図１９に示す。
【０１３８】
本実施の形態に使用された直角復調は、一般に、Ａモード信号に正弦関数と余弦関数を乗算すること、及び低域通過フィルタリングを含む。直角復調の目的は、伝送搬送周波数の除去を通して、Ａモード信号をベースバンド化することである。この方法を一般的に図２０に例示する。図２０に示したように、直角復調は、受信された後方散乱「生」Ａモード信号に複素指数関数ｅｘｐ（−２πｆｃｔ）を乗算することからなり、ここにｆｃは信号の伝送中心周波数を表す（ステップ２００４）。これは、複素指数関数とＡモード信号の和周波数と差周波数を表す一連の値を生じる。次に、その結果生じる信号を低域通過フィルタして（ステップ２００６）差（ベースバンド）周波数成分のみを保持する。指数関数は複素数であるから、それぞれ、ＩとＱで指定される同相チャネルと直角相チャネルの両方を生じる。このプロセスを図２１にグラフで描いてある。
【０１３９】
例示の実施の形態では、正弦波の周波数は、Ａモード信号の中心周波数ｆｃに相当する１６ＭＨｚである。正弦波のサンプリング間隔は、６４ＭＨｚの実効レートでサンプリングされるディジタル化Ａモード信号のそれに等しい。正弦波と余弦波による乗算は、次のようである。
【０１４０】
【数２】
2004521660.tif 000003
【０１４１】
ここにａ（ｎＴ）はサンプリング間隔Ｔ＝１／Ｆｓを持つ「生」Ａモード信号を表し、かつＦｓはサンリング・レート（例えば、６４ＭＨｚ）を表す。この乗算は、Ａモード線の全持続時間に遂行される。同相成分と直角相成分、ＸＩ（ｎＴＤ）とＸＱ（ｎＴＤ）は、それぞれ、Ｙｃｏｓ（ｎＴ）、Ｙｓｉｎ（ｎＴ）を低域通フィルタすることによって得られる。ここに、ＴＤは、デジメートされたサンプリング間隔を表す。
【０１４２】
低域通過フィルタリング（図２０のステップ２００６）は、例示の実施の形態では、３７係数からなる有限インパルス応答（ＦＩＲ）フィルタで以て完遂される。フィルタは、８０％帯域幅を持つＡモード・トランスジューサに相当する６ＭＨｚの通過帯域遮断を備えて設計される。これは、サンプリングにおいて４分の１倍だけのデシメーションをするようになっている。すなわち、フィルタ出力は、１６ＭＨｚ（＝ＴＤ）のレートで生成される。周波数応答を図２２に示す。この模範的低域通過フィルタと関連した係数を図２３に示す。
【０１４３】
例示の場合、フィルタ出力は、６４ＭＨｚの実効レートでサンプリングされる元の「生」Ａモード信号と比較して、都合良く、約１６ＭＨｚのレートでサンプリングされるだけで済むことに注意されたい。
【０１４４】
本発明の本実施の形態では、先に説明した複素包絡線の大きさは、後の処理に必要とされる唯一の成分である。もっとも、処理済みの又は未処理のＡモード信号の他の成分も、所望ならば、他の目的に使用される。複素包絡線の大きさを得ることは、低域通過フィルタされた結果の平方の和の平方根を取ることを必要とする。すなわち、結果は複素数であるから、その大きさは、式４に示したように、平方同相成分と平方直角成分との和の平方根に等しい。
【０１４５】
【数３】
2004521660.tif 000004
【０１４６】
しかしながら、包絡線の大きさを計算するのではなく、続く信号処理は、本明細書に「包絡線平方」と称する、非平方根（ｎｏｎ−ｓｑｕａｒｅ−ｒｏｏｔｅｄ）信号に適用されることがある。この後者のアプローチは、平方根演算を遂行する必要を軽減し、それによって結果のアルゴリズム実施をいくらか簡単化する。包絡線平方波形を図２４に示す。
【０１４７】
包絡線平方Ｅｓｑ（ｎＴＤ）は、同相成分と直角成分から導出される。これは、式５に記述したように計算される。
【０１４８】
【数４】
2004521660.tif 000005
【０１４９】
動脈壁の内部の管腔と関連した信号のレベルをより良く査定するために、（１０を底とする）対数が図２５に示したように得られる。管腔と関連した信号２５０２のレベルは、動脈壁からの後方散乱と関連した信号２５０４より大略３０ｄＢ下にあることに注意されたい。
【０１５０】
しかしながら、納得するように、或る種の応用では、管腔検出方法論がＡモード信号の絶対値にもっぱら基づいていたとしたならば、「偽」管腔検出が起こるかもしれない。例えば、絶対値のみに基づく管腔検出判定基準は、管腔がない所に管腔があると又は逆にそれが実在している所にないと誤って検出するかもしれない。このような偽検出は、例えば、トランスジューサと血管との間に介在した組織（静脈、筋肉等）からの後方散乱と反射、測定中の被検体の動き、等々を含む種々のファクタから起こるおそれがある。それゆえ、下に更に詳細に説明するように、本発明の方法論は、いずれか所与の管腔検出の真実性を確認するために管腔によって生じた信号アーチファクトに近い信号アーチファクトを更に検査し、かつ管腔とこれに隣接する管壁の両方の位置を正確に検出する。
【０１５１】
非線形深さ依存利得（図１７ａのステップ１７２６）が音響Ａモード信号が皮膚、組織、及び血管を通して伝搬するに連れてのその信号の減衰を償うために包絡線にまた適用される。この深さ依存利得は、しばしば時間利得補償（Ｔｉｍｅ  Ｇａｉｎ  Ｃｏｍｐｅｎｓａｔｉｏｎ）、つまりＴＧＣと呼ばれ、かつ深さは先に説明したようにこの応用では陰に関係しいる。利得に関する式は、下の式６によって与えられる。
【０１５２】
【数５】
2004521660.tif 000006
【０１５３】
ここに、深さはミリメートルで表されると想定し、かつ利得は包絡線平方に適用される。しかしながら、納得するように、他の利得関数も所によっては又は上の式６のそれと一致してのどちらかでたやすく適用されることがある。このような利得関数は、臨床試験を通してのように実験的に決定されるか、又は数学的誘導、又はＡモード信号の又はこれに関する他のパラメータの逸話的や歴史的サンプリングのような他の手段を介して決定されることがある。
【０１５４】
利得関数のＡモード信号の包絡線平方への適用は、図２６に描いてある。図２６で、実曲線２６０２はＴＧＣを適用されたＡモード信号に相当するのに対して、破曲線２６０４はＴＧＣなしでのＡモード信号に相当することに注意されたい。
【０１５５】
更に納得するように、前述の議論は直角復調、包絡線平方計算、及びフィルタリングの面で行われているが、他の型式の信号処理技術及び条件調整技術も本発明と両立して採用されることがある。
【０１５６】
信号処理の完了した後、図１７の方法１７００の次のステップ１７０６は、選択された血管（例えば、動脈）内の管腔を各Ａモード線について検出する。１模範的実施の形態では、管腔は図２７ａの一般的方法論（以下「プラトー（ｐｌａｔｅａｕ）法」）に従ってステップ１７０６で検出される。図２７ａに例示したように、管腔検出のプラトー法２７００は、一般に、深さの関数としてパワー計量をまず決定する（ステップ２７０２）を含み、パワー計量は前述の包絡線平方関数に基づいている。図２７ｂは、単一Ａモード線についての深さの関数としての後方散乱パワーを例示する。
【０１５７】
次に、パワー計量が下の式７に従って積分される。すなわち、
【０１５８】
【数６】
2004521660.tif 000007
【０１５９】
ここに、Ｐｉｎｔ（ｋ）は包絡線平方列のｋ番サンプルまで積分されたパワーを表し、及びＰ（ｋ）は包絡線平方列のｋ番サンプルを表す。
【０１６０】
プラトー法２７００の第２ステップ２７０４は、管腔内血液からの弱い後方散乱と関連したプラトーの最良推定値を識別することからなる。このステップ２７０４の部分として、積分されたパワー関数が、図２７ｃに描いたように、単位へまず正規化される（ステップ２７０６）。次に、ステップ２７０８で、ゼロから単位までの範囲が所与の数の間隔（例えば、本実施の形態では、８０等間隔、すなわち、０．０１２５の間隔、それであるから第１間隔は０から０．０１２５、第２間隔は０．０１２５から０．０２５０、等々のよう）に分割される。納得するように、他の数の間隔、及び事実、（対数関数のような関数によって決定されるもののような）非等間隔も、所望ならば、また充分に使用される。
【０１６１】
次に、プラトーがステップ２７１０によって検出される。まず、正規化されたパワーの各間隔内に起こるサンプルの数が計算される（ステップ２７１２）。この計算の背後の根本理由は、例えば、血液からの後方散乱に相当して、パワーが低いとき、１つの積分から次へと交差するように多くの値を積分しなければならないことである。それゆえ、多くのサンプルを含む間隔はプラトー（及び管腔）を表示する。図２７ａの例示の実施の形態では、３つの連続する間隔に含まれた点の数からなる試験が適用される（ステップ２７１４）ので、計量は下の式（８）に示すようになる。
【０１６２】
【数７】
2004521660.tif 000008
【０１６３】
ここに、Ｍｃｎｔ（ｋ）はｋ番間隔で終わる３つの連続する間隔内のカウントの和を表し、及びｍ（ｎ）はｎ番間隔内のサンプルの数を表す。この計量は、橈骨動脈の或る１つの横掃引の間のＡモード線数５０について図２７ｄに示してある。間隔番号６０より上の間隔内の比較的大きいピーク２７３０に注意されたい。これらは深い深さで起こる低パワーに一般に相当していて、しばしば管腔の偽検出を生じる。
【０１６４】
深い深さの方の低パワーに因るプラトーが偽管腔検出をトリガしないことを保証するために、追加項が前掲の計量に加えられた。特に、計量は、３つの連続する間隔に後動脈壁に相当する低カウントを持つ間隔を続けなければならいという要件を含む。それゆえ、修正された計量は、式９によって与えられる。
【０１６５】
【数８】
2004521660.tif 000009
【０１６６】
カウントｍ（ｋ＋２）、すなわち、３つの連続間隔のグループに先んじる２つの間隔が式９に含まれていることに注意されたい。単一Ａモード線に対するこの計量のプロットは、本例の正規化され積分されたパワーについて図２７ｅに示してある。
【０１６７】
ピークの検出の結果、管腔に相当する公算の最も大きいそれら３つの連続する間隔の識別が行われる。ピークに相当するそれら３つの間隔に先立つ全ての間隔内のサンプルの合計カウントは、管腔の深さを識別する。しかしながら、１模範的実施の形態では、その後の動脈壁探索の目的のために管腔と関連したＡモード線内のサンプル点が第１間隔からそれら３つの連続する間隔内の第１間隔までの（かつこれを含む）全ての間隔内のサンプルの合計数から導出される。それゆえ、所与のＡモード線についての管腔と関連したサンプルは、式１０によって推定される。
【０１６８】
【数９】
2004521660.tif 000010
【０１６９】
ここに、Ｌは管腔面積に先立つサンプルの数を表し、及びＫは３つの連続する間隔のうちのピークが検出された第３間隔に対する「ポインタ」を表す。
【０１７０】
更に、管腔パワーの平均の推定値がそれら３つの連続する間隔の１つ以上に含まれたパワーから導出されることがある（ステップ２７１６）。例えば、１模範的アプローチでは、それら３つの連続する間隔の第２間隔に相当する包絡線平方のサンプル値を単に平均することによって、平均管腔パワーがこの第２間隔から決定される。管腔後方散乱パワーのこの推定値は、下に更に詳細に説明するように、動脈壁を検出するために検出された管腔の推定された位置と関連してその後使用される。
【０１７１】
注意するのは、積分されたパワー・プロファイル中のプラトーは他の源からも生じることである。特に、人間橈骨動脈の場合、このようなプラトーは、圧平デバイスと注目の血管との間に位置した１つ以上の静脈の存在によって、又は嚢胞（ｃｙｓｔ）や他の成長物の存在によって、また誘発されることがあることである。より普通であるが、静脈に因るプラトーは動脈の圧縮と関連した注目の「真」プラトーと比較して一般に極めて小さく、これは比較的的小さい静脈直径（と壁厚さ）に大部分因っている。更に、本発明の方法論は、比較的薄い壁の静脈を凹ませるために必要なレベルよりも上の圧平圧力で血管直径（及び積分された後方散乱パワー）の測定を実施することによってこの問題点を有効に克服する。静脈のこのような凹みは、大きな血管（例えば、橈骨動脈）の直径に重大に影響するために必要な圧力の充分下の圧力で一般に起こり、それによって、静脈に因る潜在的「偽」プラトーの問題に対する洗練された解決を用意する。
【０１７２】
静脈と対照的に、嚢胞又は成長関連のアーチファクトは、その母集団の非常に小さな分数部分にのみ一般に起こり、かつトランスジューサを動脈をおおう他の位置に再配置するとそれらから起こるいかなる作用も除去する程度に、しばしば空間的に局所化されている。更に、嚢胞性区（ｃｙｓｔｉｃ  ａｒｅａ）は、相当するパルス圧又はパルス圧に関係した運動成分を持っておらず、それゆえ当業者に周知の多くの信号処理技術を使用してたやすく識別しかつ選別する（ｓｃｒｅｅｎ）ことができる。
【０１７３】
第２実施の形態では、図１７のステップ１７０６での管腔の検出は、以下に「間隔法」と称する、Ａモード線に沿う所定距離（例えば、１ｍｍ）内の平均パワーを計算することによって、完遂される。伝搬速度が１５００ｍ／ｓかつ式１１から導出されると想定すると、１ｍｍ間隔は、１．３３μｓのＡモード線に沿う時間間隔に相当する。
【０１７４】
【数１０】
2004521660.tif 000011
ここに、
Ｄ＝深さ、
ｃ＝伝搬の速度
ｔ＝伝搬の時間
【０１７５】
式１１中の２なる係数は、双方向音響波伝搬を勘定に入れる。先に説明した１６ＭＨｚのサンプリング・レートに対しては、１．３３μｓは、デサンプルド（ｄｅ−ｓａｍｐｌｅｄ）Ａモード線に沿う約２２サンプルに相当する。
【０１７６】
パワー計算用１ｍｍ間隔は、例示の実施の形態では、検査下の血管（例えば、橈骨動脈）の直径が成人母集団の本質的に全てで１ｍｍを超えるという、本願の讓受人による観察に基づいて選択されている。それゆえ、パワー計算は、管腔にのみ寄与する信号を構成するように合理的に保証される。しかしながら、納得するように、この間隔は、非成人母集団の測定、他の血管の使用、及び／又は他の種についてのその技術の使用のようなファクタに基づいて調節されることがある。また注意するのは、上に挙げた１ｍｍ間隔は、また都合良く、Ａモード信号パワーの推定値について平均するサンプルの「合理的」数を提供する。
【０１７７】
１６ＭＨｚのサンプリング・レートに対しては、前述の１ｍｍ間隔は、１．３３μｓ、又は約２２サンプルに相当する。時刻ｎＴＤでのパワーの推定値Ｐ１ｍｍ（ｎＴＤ）は、式１２に従って計算される。すなわち、
【０１７８】
【数１１】
2004521660.tif 000012
【０１７９】
ＴＧＣで補正されたＡモード・データの線に沿う１ｍｍ間隔内の音響パワーの模範的プロットは、図２８ａに描いてある。
【０１８０】
橈骨動脈の内部の管腔（例えば、血液）の位置は、１ｍｍ平均パワー計算における第１最小として検出される。この最小に対する探索は、Ａモード線に沿う所定点（例えば、１ｍｍ）で開始し、かつ本例示では７ｍｍに選択された第２所定深さまで続く。１ｍｍから７ｍｍのこの範囲は、橈骨動脈の前壁（すなわち、トランスジューサから放射された音響エネルギーが最初に出会う壁）が成人母集団の本質的に全てにおいてこの間隔内に起こるという観察に基づいて選択された。また注意するのは、前述の最小は図２８ａのプロット内で明確に区別可能であり、それによって、都合良く、良品質（すなわち、低雑音レベルと低他のアーチファクト）のＡモード信号での容易な検出を見込む。電子信号内の最小の検出は、電子技術上周知の多くの技術によって完遂されるといえ、したがって、本明細書には説明しない。
【０１８１】
図２８ｂは、前述の「間隔」法を使用する管腔検出をグラフで例示する。図２８ｂに例示したように、管腔検出は、事前選択された深さで、本場合は１ｍｍで動脈壁の探索を開始することを含む（ステップ２８０２）。最小パワーは、深さ＝１ｍｍに対するパワーで初期化される。次に、最大指定深さ（前述の例では、７ｍｍ）に対する関係がステップ２８０４によって決定される。最大深さを超えたならば、最小は検出されなかったのであり、かつ探索はステップ２８０５で終結され（かつオプショナルに再開始される）。最大指定深さをステップ２８０４で超えていなかったならば、次のパワー・サンプルがステップ２８０６によって分析される。最小値未満であるならば、最小が現在深さのサンプルと関連したパワー・レベルＰにリセットされ、かつカウントがゼロにリセットされる（ステップ２８０８）。逆に、次のパワー・サンプルが現在最小以上であるならば、カウントは１だけ増分され（ステップ２８１０）、かつ整数のカウントが所定カウント制限パラメータ（例示の実施の形態では５）に対するその関係を決定するために分析される（ステップ２８１２）。現在カウントがカウント制限パラメータを超えるならば、第１最小が現在深さで検出されたと想定され、かつ探索は終結される（ステップ２８１４）。検出されないならば、ステップ２８０４から２８１４が適用可能であるに従って再び繰り返される。
【０１８２】
前述の手順の間中、管腔からの後方散乱Ｅｓｑ（ｂｌｏｏｄ）と関連したＡモード信号は、管腔の領域がいったん識別されたならば、容易に識別することができる。推定値は、管腔が検出された深さに中心を置く範囲にわたる包絡線平方のサンプルの単なる平均として計算される。平均は、式１３の模範的関係に従って計算される。
【０１８３】
【数１２】
2004521660.tif 000013
【０１８４】
ここにＮｂｌｏｏｄは、血液検出の深さに相当する包絡線平方の指標を表す。式１３で、平均は包絡線平方の１５サンプルにわたって計算されることに注意されい。もっとも、納得するように、特定応用次第で他の値で置換されることもある。
【０１８５】
例えば、プラトー法又は間隔法のどちらかを使用して管腔が検出された後、前動脈壁と後動脈壁が検出され、かつそれらの空間位置（すなわち、センサ・ユニットに対する深さ）が推定される。特に、前動脈壁は、センサ・ユニットと検出された管腔の深さとの間、特にＡモード信号が管腔と関連した積分されたパワー又は信号レベルへ低下する所の付近に一般に起こると想定される。この想定は、血管の管腔が壁に直ぐ隣接している（すなわち、トランスジューサに関して壁の遠い方の側上）にあるという知識に基礎を置いている。同様に、後動脈壁は、管腔信号の深さよりも深い深さかつ、積分されたパワー・レベル又は信号レベルが管腔と関連したレベルの上へ鋭く立ち上がる所の付近で一般に起こると想定される。動脈の前壁と後壁に相当するＡモードの積分されたパワー・レベル又は信号レベルのこれらの変化は、例えば、ＴＧＣ・Ａモード・データに容易に観測可能である。
【０１８６】
第１模範的実施の形態（図２９ａ）では、壁検出の方法は、（ｉ）動脈壁に対する探索をそこから開始する管腔内の検出された位置、及び（ｉｉ）上の図２７ａに関して説明した検出しきい値をそこから導出する平均管腔パワーの推定値を利用する。Ａモード線に沿う管腔の検出された位置は、後壁に対する探索の開始として役立つ。この位置から開始して（ステップ２９０２）、積分されたパワーは、Ａモード信号の包絡線の平方の連続するサンプルを加算することによって計算される（ステップ２９０４）。定式化は、式１４によって与えられる。
【０１８７】
【数１３】
2004521660.tif 000014
【０１８８】
ここに、ＰｂｋはサンプルＫｓｔで開始してサンプルｍまで後壁の方向における積分されたパワーを表し、及びＰ（ｋ）はＡモード信号の包絡線の平方のｋ番サンプルを表す。Ｐｂｋ（ｍ）が規定しきい値を超えるとき、後壁が検出されたと想定されかつ相当するサンプル値ｍは後壁の深さを表す（ステップ２９０６）。Ａモード線に沿うサンプル数と深さは関係していることに注意されたい。実深さは、式１５に示したようにサンプル数から導出される。
【０１８９】
【数１４】
2004521660.tif 000015
【０１９０】
ここに、ｃは伝搬の速度を表し、Ｎｓａｍｐはサンプル数を表し、及びＦｓは包絡線平方列と関連したサンプリング・レートを表す。
【０１９１】
前壁は、後壁とほとんど同じ仕方で検出される。後壁に対し使用された同じ位置から開始して、積分されたパワーは、前壁の方向にＡモード信号の包絡線の平方の連続するサンプルを加算することによって計算される。それゆえ、サンプルは、サンプルＫｓｔから降順に取られる。積分されたパワーが規定しきい値を超えるとき、前壁が検出されたと想定されかつ相当するサンプル値は前壁の深さを決定する。
【０１９２】
図２９ａの実施の形態では、前述の管腔平均パワーの５００倍に等しい前壁しきい値、及び管腔平均パワーの５０倍に等しい後壁しきい値が利用される。これらの値は、本願の讓受人によって指定深さ依存利得（例えば、ＴＧＣ）を用いて実験的に決定された。しかしながら、納得するように、他のしきい値及び他の型式の判定基準も使用されることがある。例えば、他のパラメータ（例えば、管腔内のクラッタ又は介在した組織内の減衰）の関数として変動する値を有する可変しきい値で置換される。更になお、納得するように、前壁のしきい値と後壁のしきい値の間の関係は、所定関係を考慮に入れることも入れないこともある、例えば、前壁しきい値は後壁しきい値の関数として決定されるか、又はこれに代わりそれから完全に独立である。他の変形も可能である。例えば、下のスコア付け（ｓｃｏｒｉｎｇ）と加重についての議論を見られたい。
【０１９３】
更になお、処理を容易にするために、信号処理技術上周知であるように、Ａモード線の線間平均が管腔検出と壁検出に先立って遂行されることがある。この技術は、管腔検出及び動脈壁深さ推定の両面で、システムの性能と堅牢性を改善するために使用することができる。
【０１９４】
動脈壁を検出する方法の実施の形態（図２９ｂ）では、包絡線平方の振幅Ｅｓｑ（ｎＴＤ）が血液からの後方散乱Ｅｓｑ（ｂｌｏｏｄ）と関連した推定された信号レベルと比較される。これは、先に説明した積分パワー計算アプローチと明確に区別される。前のように、前壁に対する探索は、トランスジューサ素子への方向、すなわち、深さを減らす方向である。包絡線平方アーチファクトの振幅が指定された判定基準を満たす（例えば、推定された信号レベルを超える）とき、この振幅は前壁からの後方散乱に寄与する。この生起の深さは、前壁の位置であると指定され、所望ならば、オプショナル後続の確認プロセスに付される。これに代わりに（又は同時に）、包絡線平方変量の振幅に関係する以外の判定規準が採用されることがあることに注意されたい。
【０１９５】
後壁に対する探索は、前壁について説明したのと類似の仕方で遂行される。特に、探索は、検出された管腔の深さで開始されるが、しかし深さを増す（すなわち、超音波トランスジューサから遠ざかる）方向に進行する。前壁の場合のように、包絡線平方の振幅が指定された判定基準を超える、満たすとき、この振幅は後壁からの後方散乱に寄与する。この生起の深さは、後壁の位置であると想定される。
【０１９６】
したがって、Ａモード・データが「合理的な」品質のものであるとき、簡単な探索アルゴリズムが前壁と後壁とに関連した位置を識別するために信号レベル（包絡線平方）方法論と一緒に使用される。特に、図２９ｂに例示したように、方法２９５０は、管腔信号が検出された深さから出発して（ステップ２９５２）、管腔の信号レベルを（ｉ）トランスジューサへの方向と、（ｉｉ）トランスジューサから遠ざかる方向の両方でＡモード信号レベルと第１比較することを含む。信号レベルは、選択された方向で測定される（ステップ２９５４）。信号レベルが１つ以上の所定判定規準（例えば、上の（ｉ）又は（ｉｉ）の信号レベルが管腔と関連した信号レベルの「ｎ」倍ほどである）を満たすことが発見されるとき、信号は動脈壁に相当すると想定される（ステップ２９５６）。トランスジューサの方へのＡモード線に沿う点であってこの判定規準がそこで満たされる点は前動脈壁の位置として表され、及び同様に、トランスジューサから遠ざかる方へのＡモード線に沿う点であってこの判定規準がそこで満たされる点は後動脈壁の位置として表される。これは、図３０にグラフで描いてあって、平均パワー計算３００２とＡモード包絡線３００８内に存在する前壁アーチファクト、後壁アーチファクト３００４、３００６との間の関係を例示する。
【０１９７】
包絡線平方を計算し、指定された間隔（例えば、１ｍｍ）での相当する平均パワーと管腔信号を決定し、かつ前動脈壁と後動脈壁を検出するプロセスは、本実施の形態では、Ａモード「一線一線（ｌｉｎｅ−ｂｙ−ｌｉｎｅ）」式で遂行される。５０００線からなる模範的Ａモード・データ集合に対してこの仕方で決定された動脈壁位置を図３１に示す。しかしながら、このような一線一線分析は、例えば、１つの線の部分、又は多数の線のグループ化（連続的であろうとそうでなかろと）、又は線平均のような他の分析的アプローチ又は増分で置換されることがある。
【０１９８】
いったん相対的壁位置が決定されると、血管の直径は、後壁と前壁の深さの差を取ることによって簡単に決定されるといえる。例えば、後壁が３．５ｍｍの深さで起こり、かつ前壁が１．５ｍｍの深さで起こるならば、血管の直径は、（３．５ｍｍ−１．５ｍｍ）＝２．０ｍｍで推定することができる。この決定は、１実施の形態では、本明細書に先に説明したように導出される平滑化された壁位置推定値（すなわち、他の信号処理技術を使用して平滑化のために数学的に平均されるもの）を使用して行われる。
【０１９９】
納得するように、壁検出の前述の模範的方法論（例えば、積分パワーと包絡線平方信号レベル）は血管の前壁と後壁の面で説明されたが、前壁又は後壁の検出アプローチは、別々に及び／又は孤立して応用されることがある。例えば、評価される或る種の被検体（又は或る共通の特性を有する被検体のグループ）について前壁アーチファクトが特に弱い、又はそうでなければ非信頼性であることが決定されるならば、後壁検出が補償上より実質的に加重される。同様に、検査下の血管内のクラッタが顕著であるならば、より多く影響された壁が加重の面で選択的にデエンファシスされる（ｄｅ−ｅｍｐｈａｓｉｚｅｄ）。他のこのような修正も特定応用次第でまた採用されることがある。
【０２００】
更に、納得するように、積分されたパワーしきい値アプローチ及び包絡線信号レベル・アプローチは、単独で又は組合わせで使用されることがあり、又は他の方法（例えば、信号処理技術上周知の型式の、いわゆる「ボックスカー（ｂｏｘｃａｒ）」平均）が採用されることがある。例えば、包絡線平均信号レベルと積分されたパワーの並列計算が遂行され、かつ計算された値が応用可能のしきい値判定基準（例えば、ぞれぞれ、４倍管腔信号レベル、及び５００倍平均管腔パワー）と比較される。次いで、これらの比較の結果は、システム及び／又はユーザが後に使用することができる「信頼（ｃｏｎｆｉｄｅｎｃｅ）レベル」計量の計算によるような、壁検出における相対的な信頼性（ｒｅｌｉａｂｉｌｉｔｙ）又は信頼（ｃｏｎｆｉｄｅｎｃｅ）を決定するために使用される。積分されたパワー・アプローチ及び包絡線平方アプローチの両方に対する壁検出スコアが高いならば、結果の信頼計量は高い。壁検出のスコアが１つの技術に対して高くかつ他に対して低いならば、信頼計量の値は減少させられ、かつその後の確認処理が表示される。
【０２０１】
更になお、納得するように、全体として壁検出の文脈の中で、信号処理技術上知られたスコア付けと加重の種々の形式が前述の判定基準を置換して、又はこれと一緒に使用されることがある。このような技術は、都合良く、実臨床使用下でのシステムの堅牢性を増す。例えば、１代替実施の形態では、離散値（例えば、信号レベルの場合には４倍推定値信号レベル、又は積分アプローチの場合には５００倍平均管腔パワー）を満たす又は超える包絡線平方パワー値又は積分されたパワー値を有するＡモード信号は、「１．０」のスコアを割り当てられる。この離散値よりも低い値を有する信号は、窓関数ｗ（ｘ）に対するそれらの関係に基づいてスコアを割り当てられるので、所定しきい値（例えば、２倍の推定された信号レベル、又は２５０倍平均管腔）の下にある値は０．０のスコアを与えられ、及び所定しきい値と設計された判定基準との間の値は関数ｗ（ｘ）に従って０．０と１．０の間のスコアを割り当てられるというようになる。
【０２０２】
なお他の代替実施の形態では、所与のＡモード測定と関連した包絡線平方パワー値又は積分されたパワー値は、その間隔中にサンプリングされたデータの品質に関係している或る外因性又は内因性「品質係数（ｑｕａｌｉｔｙ  ｆｃｔｏｒ）」に基づいて加重されることがある。簡単な例として、評価される被検体が所与の一連のＡモード線の間中に動くならば、データの品質は低下させられ、したがって、このデータに基づいて遂行されたいずれの管腔／壁検出計算も他のサンプルに比較して加重を人工的に下げられる。
【０２０３】
なお他の代替実施の形態では、或る決まった深さに相当するＡモード線によって発生された包絡線平方パワー値又は積分されたパワー値が調節されるというように、「位置的（ｌｏｃａｔｉｏｎａｌ）」加重及び／又は位置的スコア付けが適用されることがある。例えば、先に説明したように、成人における皮膚の表面から１ｍｍ未満に橈骨動脈の前壁を見付ける公算は極端に低く、したがって、（１ｍｍの深さに相当する）最初の「Ｎ」Ａモード線から導出された包絡線平方パワー値又は積分されたパワー値は、大幅に切り下げられるか又は除去さえ（すなわち、ゼロ加重）されることになる。
【０２０４】
同様に、逐次間隔に次ぐ間隔（ｉｎｔｅｒｖａｌ−ｔｏ−ｉｎｔｅｒｖａｌ）処理は、不連続Ａモード線間隔に決定された壁アーチファクトが解決を必要とするあいまい性に従ってマークされるというように、実施されることがある。前壁アーチファクトが最初の１ｍｍ深さ間隔に検出され、かつ再び第３間隔に検出されるならば、（ｉ）最初の間隔に検出されたアーチファクトは前壁でありかつ第３間隔で検出されたアーチファクトは後壁であるか、又は（ｉｉ）最初のアーチファクトは雑音であった、かつ第３間隔アーチファクトは真前壁信号であるかのいずれかを、信頼の或るレベルを以て想定することができる。このようなあいまい性は、成人母集団中の生起の低公算に基づいて第１アーチファクトを除去する前述した「位置的」加重の応用、又はそれらの間隔に対する後続Ａモード信号のサンプリングの使用のような、多くの技術を通して解決することができる。
【０２０５】
本方法論は、必要な所での有意な平滑化／フィルタリングをまた含む。このような平滑化／フィルタリングは、更になお、心収縮期から心拡張期へと起こる時間変動ばかりでなく、呼吸によって引き起こされる拍動間（ｂｅａｔ−ｔｏ−ｂｅａｔ）の変動を除去する。このような平滑化／フィリタリングは、信号処理技術上一般に充分理解されており、したがって、本明細書に更に説明しない。
【０２０６】
納得するように、本明細書に説明した壁検出方法論の正確性は、細い超音波伝送ビーム及び受信ビームを生じる能力に或る点で結び付いている。特に、２方向ビーム応答は、超音波トランスジューサが評価される血管に対してその最適位置の付近にあるとき、その血管の管腔（すなわち、血流量の断面積）内をヒットするためには充分に狭くあるべきである。そうでなければ、血管の側壁からの後方散乱が管腔の検出を妨げ、かつ前壁と後壁の空間位置の推定値の正確性を劣化させる。
【０２０７】
前壁と後壁の位置の評価で他の考慮することは、圧平の間中の組織の相対コンプライアンスと相対圧縮性である。図１７〜３１に関して上に説明したような優勢に静的な壁位置評価の下では、組織内の血管壁の絶対位置と相対位置は、顕著に変動するということはない。しかしながら、（下の図３４〜３７に関して論じる）このような壁位置決定に基づく血管の血行力学パラメータを査定するとき、そこでは血管及びその取り囲む組織の少なくとも部分が本明細書に図６に関して先に説明したように組織の圧平の間中圧縮されるようになっており、前壁と後壁の相対位置と絶対位置は変化する。
【０２０８】
図３２は、本願の讓受人によって得られた実データに基づいて、前述の問題をグラフで例示する。図３２に示したように、血管の後壁３２０２の位置は、ゼロＡモード線（すなわち、走査の開始）で深さ約４ｍｍから６００Ａモード線の後約３ｍｍの深さへ変動する。後壁の深さの変動は、直線並進（ｍｍ）が時間の関数として変動する（図３２ａ）圧平デバイス１２０４の運動によって主として誘発される。圧力トランスジューサ（実圧平面を被検体の組織に接触している）上に記憶されたｍｍＨｇでの圧力を図３２ｂに示す。それゆえ、圧平デバイスと圧力トランスジューサが圧力を血管の上の皮膚の表面に印加するに連れて、いくつかの効果が起こる。すなわち、（ｉ）皮膚、筋肉、及び静脈を含む介在した組織が圧縮する、（ｉｉ）血管の背後の組織（すなわち、後壁の遠い側上）が圧縮する、及び（ｉｉｉ）血管自体が圧縮しかつ、先に説明したようにその固有のフープ・ストレスを克服する。これら３つのファクタは、図３２に示したように前壁と後壁の深さに変化を生じるばかりでなく、血管の直径に変化を生じる。図３３は、圧平の関数として、血管の前壁３３０２と後壁３３０４の位置の変化ばかりでなく、血管（管腔）直径３３０６の変動をグラフで例示する。図３３は本発明の原理をより明確に例示する平滑化データを表すが、実データ又は未平滑化データは本来的に遥かに雑音性であることに注意されたい。
【０２０９】
図３４は、上に説明した前壁と後壁検出方法を使用して血管（管腔）直径の変化の速度を決定する１方法論を例示する。変化の速度は、所与の期間にわたって血管直径の変化を取りかつその同じ間隔中のＡモード線の数で除算することによって簡単に計算される。この後者の数を圧平デバイスの位置の線形変化（図３３ａ）に相関させることもあり、それによっていろいろな種類の伝達関数（すなわち、線形圧平の関数としての血管直径の変化の速度）を発生し、伝達関数は圧縮中の組織と血管の挙動を記述します。
【０２１０】
後方散乱超音波エネルギーを使用して血行力学性質を査定する方法
図３５を参照して前述した後方散乱超音波エネルギー（例えば、Ａモード信号）を使用して血行力学性質を非侵襲的に査定する方法の１模範的代替実施の形態を説明する。納得するように、次の議論と模範的代替実施の形態は人間の橈骨動脈を使用しての動脈血圧の非侵襲的測定の面で行われるが、本明細書に記載した本発明の方法はより一般には、他の血行力学パラメータ、血管、及び種に応用されることがある。
【０２１１】
図３５に示したように、方法３５００は、一般に、（本明細書に図１２に示したそれのような）適当な装置を査定される血管の領域に全体的に位置決めすることを含む（ステップ３５０２）。手頸上の橈骨動脈のような、装置のトランスジューサを置くサイトは、医務者（ｃａｒｅｇｉｖｅｒ）には一般に知られている。次に、ステップ３５０４で、装置が活性化されて、血管及び介在した組織への装置の最適結合を達成するように、トランスジューサを血管に対して横位置決めする。納得するように、横位置決めするための多くの方法が使用されるといえる。例えば、本明細書に先に言及した、又は図４０〜４１に関して本明細書に後に説明する讓受人の同時係属米国特許出願第０９／３４２，５４９号に説明された、観察された反射最小に基づく横位置決めの方法が単独又は他の方法と一緒に採用されることがある。
【０２１２】
「横」位置決めの行為は、縦位置決めの（すなわち、血管の縦軸に沿う）或る成分をまた含むことがあることに注意されたい。というのは、被検体の手頸／前腕上の装置の配置は手頸／前腕内の血管の配向に反して、手頸の肉体的属性によってより多く統制されるからである。特に、橈骨動脈が手頸の骨の軸に完全に平行な方向に走っていない位置にトランスジューサによって測定される点が起こる場合、このような「横」位置決めは、本来的に縦成分をまた充分に含む。更になお、血管に沿う或る種の点は、（例えば、超音波の伝送／反射と潜在的に干渉する静脈や、嚢胞や、他のセグメントの存在に因り）血行力学分析により良く適していることがある。
【０２１３】
次に、初期校正圧平掃引（すなわち、無圧縮から全圧縮へ又は血管の閉塞へ印加圧力を掃引すること）がステップ３５０６で遂行されて終末心拡張期圧プロファイルの傾斜の「平坦化（ｆｌａｔｔｅｎｉｎｇ）を識別し、したがって、圧力波形内の注目の点を局所化するのを援助する。図３５ａに示したように、図３５ａに示したように、終末心拡張期波形３５５４の平坦化領域３５５２は、一般に、急な勾配の２つの領域３５５８と３５６０との間に起こる。もっとも、これは常にたやすく検出可能とは限らない。更になお、注意するのは、多数の平坦領域が単一圧平掃引内に起こることがあることである。
【０２１４】
ステップ３５０８で、上側と下側の圧縮限界（本実施の形態では検査下の血管の直径減少のパーセンテージとして表される）が先行圧平掃引の間中に得られた情報に基づいてセットされる（ステップ３５０６）。直径の減少は、図１７〜３４に関して先に説明した前血管壁と後血管壁の検出技術を応用することによって決定される。例えば、２５％直径減少の下限と５０％直径減少の上限（それぞれ、４４％と７５％の断面積減少）が上側圧縮限界と下側圧縮限界として例示した実施の形態に使用される。もっとも、他の値及び／又は標印（ｉｎｄｉｃｉａ）も使用される。直径の２５％減少の低値は、介在したかつ周囲の組織の圧縮が実質的に起こりかつ血管のフープ・ストレスが克服され始める点に大略等しくするように本願の讓受人によって実験的に決定されている。直径の５０％減少の上側限界は、血管内の終末心拡張期血圧が印加された圧縮力によって克服される点に相関するように実験的に決定されている。
【０２１５】
更に、注意するのは、圧平掃引の間中に、最大拍動圧（ｍａｘｉｍｕｍ  ｐｕｌｓａｔｉｌｅ  ｐｒｅｓｓｕｒｅ）が、平坦化する圧力波形の前述の領域の右手に（すなわち、より後の時刻に／より多量の血管圧縮で）起こることになり、それによって、波形中の注目の点を局所化することを更に助ける。特に、最大拍動動圧（これは技術上周知の多くの手段によってたやすく決定される）よりも多量の圧縮の状態の下で起こる波形中の他のアーチファクトは、差し引くこともでき、それによって、圧力波形中の注目の領域を狭める。したがって、最大拍動圧は、オプショナルにまた監視される（ステップ３５１０）ばかりでなく、動脈圧に上側窓限界をセットするのを援助する。
【０２１６】
ステップ３５０６の初期圧平掃引を完了し、かつ血管を弛緩させた（３５１２）後、血管は印加圧力と血管壁のフープ・ストレスとの間の平衡に大略相当する圧縮のレベルに再び圧平される（ステップ３５１４）。圧縮のこのレベルは、フープ・ストレス成分が除去されているから、終末心拡張期圧が血管に印加された圧縮力の測定値になるという点で、意味深い。先に説明したように、圧縮のこのレベルは、下側限界（例えば、血管直径の２５パーセント減少、つまり断面積の４４パーセント減少）と、上側限界（例えば、５０パーセント減少、つまり断面積の７５パーセント減少）によって制限される。圧平の最適レベルは、経壁圧を平衡させるようになっていて、下側直径減少限界判定基準に関連した圧力と上側直径減少限界判定基準に関連した圧力との間に形成された圧力の窓内に起こることになり、その窓はその中に圧力プロファイルの「平坦化」領域の少なくとも部分を有する。圧平のこの最適レベルで圧力トランスジューサによって測定された圧力は、伝達関数を介しての補正（ステップ３５１８）の基礎として使用される。ステップ３５１８に従って伝達関数を決定しかつ適用する１模範的方法は、図３８と３９を参照して下に説明する。
【０２１７】
納得するように、上側限界判定基準と上側限界判定基準によって作り出された圧力「窓」はかなり広いことがあり、したがって、追加の方法が窓内の最適圧力をより正確に決定するためにオプショナルに適用されることがある。本実施の形態では、これらの追加方法は、図３６に関して下に更に詳細に説明するように最適圧平レベルのアルゴリズムに従う決定（ステップ３５１６）を含む。
【０２１８】
図３６を参照して、最適圧平レベルの決定（ステップ３５１６）の方法論を説明する。図３６に示したように、方法３６００は、一般に、まず、終末心拡張期圧プロファイル内の「平坦化」スポットを識別する（ステップ３６０２）ことを含む。１模範的実施の形態では、圧平掃引にわたっての圧力プロファイルの端から導出されたヒストグラムが発生される（ステップ３６０６）。図３６ａは、ステップ３６０６に従って発生された模範的ヒストグラムである。ヒストグラム中で最大ビン（ｂｉｎ）値は、都合良く、終末心拡張期圧内の平坦スポットに相当する圧力プロファイルの領域で起こり、それによって平坦スポットのたやすい検出が見込まれる（ステップ３６０８）ことに注意されたい。
【０２１９】
次に、ステップ３６０４で、終末心拡張期圧が最大拍動圧で決定される。最大パルス圧（ｐｕｌｓｅ  ｐｒｅｓｓｕｒｅ）は、例示の実施の形態では、パルズ圧データへの高次多項式適合（ｈｉｇｈ  ｏｒｄｅｒ  ｐｏｌｙｎｏｍｉａｌ
ｆｉｔ）を使用して決定され（ステップ３６１０）、次いで、ピークをそれから決定する（ステップ３６１２）。このアプローチは、呼吸作用を含むパルス圧中の雑音を最小限にするのを助ける。次いで最大拍動圧で起こる終末心拡張期圧が記録される（ステップ３６１４）。
【０２２０】
終末心拡張期圧のこれら２つの指標（すなわち、「平坦スポット」と最大拍動圧）は、所望に従って分離して又は互いに一緒に使用することができることに注意されたい。更になお、他の技術が、最適圧平圧力のより正確な決定を血行力学測定の間中維持するために、これら２の指標の１つ又は両方を置換する、又はこれと組み合わせて使用されることがある。
【０２２１】
血行パラメータを査定する方法の代替実施の形態（図３７に示した）では、方法３７００は、血管をおおってトランスジューサを位置決めし（ステップ３７０２）、次いで、その後血管を圧平する（ステップ３７０４）一方、圧力波形を監視する（ステップ３７０６）ことを含む。次いで、終末心拡張期圧波形の屈曲（ｉｎｆｌｅｃｔｉｏｎ）（例えば、勾配の変化）の点がステップ３７０８によって識別される。このような屈曲は、波形中の終末心拡張期圧値に関係した勾配の変化を識別するために適合したアルゴリズムを使用して自動的に識別される、又はこれに代わり手動で、例えば、熟練技術者による波形の視察によって、識別されることがある。前述の限界判定基準（例えば、２５％血管直径減少公称下側限界、及び５０％直径減少上側限界）が波形中の注目の点を更に局所化するためにステップ３７１０で適用される。
【０２２２】
伝達関数
図３８と３９を参照して、拡大縮小（伝達）関数を決定しかつ血行力学パラメータの測定された値に適用する（図３５のステップ３５１８）方法の１模範的実施の形態を説明する。
【０２２３】
納得するように、次の議論は測定された動脈血圧値と使用される拡大縮小関数又は伝達関数の面で行われるが、その方法論は、より一般的には、図３５のアプローチを使用して、又はそうでなく、測定された他の血行力学パラメータ又は波形に適用されることがある。例えば、被検体の循環系へのストレスの圧平が注目の状態（例えば、圧平の間中に圧力測定法によって測定された圧力対動脈直径の或るパーセンテージ減少が終末心拡張期圧と相関する）を生じるならば、本明細書に説明した方法論は、組織コンプライアンス及び他の作用を補正するために必要な拡大縮小を決定するために使用される。
【０２２４】
更になお、図３８の方法論は、一般に、データの母集団に依存する実験に基づく技術を含むが、本発明は、図３８の方法に代わりに又はこれと一緒に他の非実験的アプローチも企図しており、全てのこのような変形及び組合わせは本明細書に添付した特許請求の範囲に入る。したがって、図３８の方法は、単に、性質の上で例証である。
【０２２５】
図３８に従う血行力学パラメータ拡大縮小の方法３８００は、一般に、相関を導出するために実験基準データ（ｅｍｐｉｌｉｃａｌ  ｒｅｆｅｒｅｎｃｅ  ｄａｔａ）を使用することを含み、これらの相関は検査下の被検体から得られた血圧の現場測定に後に応用される。ストレスが実験基準データと現場測定との間で注目の１変量（例えば、終末心拡張期圧）に一致又は「定着（ａｎｃｈｏｒ）」するように、測定される血管に印加される。次いで、他の未定着変量（例えば、心収縮期圧及び平均血圧）に対する導出された補正が、実動脈内圧に事実上等価である補正された波形を得るために、評価される被検体の現場波形の相当する部分に適用される。
【０２２６】
図３８に示したように、方法３８００は、本明細書に先に説明した方法は別として、所望パラメータ（例えば、動脈血圧）の基準値又は基準波形を決定する（ステップ３８０２）ことをまずを含む。１実施の形態では、基準血圧データは、医療技術上周知の形式の「Ａ線」侵襲性カテーテル・デバイスを経由して複数の生被検体から得られる。Ａ線デバイスは、一般には、真動脈血圧（波形）を決定する正確な手段として受け入れられるが、しかし、都合悪く、被検体の循環系内へのカテーテル・プローブの挿入を必要とする。同じ患者についての圧力測定法による圧力データが同時に得られ（ステップ３８０４）、この圧力測定法による測定は、定着変量（ここでは終末心拡張期圧）の圧力測定法による測定値がＡ線によって決定された実終末心拡張期圧と可能な限り密接に整合するというように被検体の血管を圧平することを含む。
【０２２７】
１簡単な実施の形態では、圧平レベルは、圧力測定法の圧力トランスジューサによって測定された（かつ陰極線管（ＣＲＴ）又はこれに類似のもののようなディスプレイ・デバイス上にグラフフで表示された）波形の定着変量がＡ線の同じ変量の値に実用可能な限り密接に整合するまで、漸次高められる。整合は、波形の同時的重なり、又はディスプレイ・デバイス上に表示された２つの数値の差の最小化のような眼に見える手段を介して主としてオペレータによって完遂される。次いで、オペレータは、圧平レベルを可能な整合の最良レベルへ調節する。
【０２２８】
他の実施の形態では、圧平は、Ａ線波形と圧力測定法波形の両方から、注目の定着変量を、周期式に又は連続式に、自動的に計算し、かつそれらをスコア付け計量を使用して可能な限り密接に整合させることを企図するアルゴリズムによって制御される。例えば、離散時間間隔中の終末心拡張期圧どうし間の差に関係しているスコア付け計量が波形間の「適合（ｆｉｔ）」の品質を評価するために使用され、その計量は選択された時間間隔中に計量の値を最小限にするために圧力を変調するように先に説明した圧平モータ制御回路への入力として使用される。このような曲線を適合させる（ｃｕｒｖｅ−ｆｉｔｔｉｎｇ）計量は、信号処理技術における当業者に周知であり、したがって、本明細書に更に説明しない。
【０２２９】
定着変量の整合及びデータの収集の後、データの２つの集合が分析されて、ステップ３８０６によってその個人についての真動脈圧（すなわち、Ａ線値）と圧力測定法によって得られた圧力との間の関係の特性表示を用意する。基準値と圧力測定法のデバイスから導出されたものとの間の差は、なかでも、圧力測定法による測定の間中の介在した組織の圧縮に因る転送損失から生じた差である。特に、例示の実施の形態では、心収縮期圧と平均動脈血圧に対する補正は、２つの波形中の終末心拡張期圧が整合していた事実に基づいて計算される。次いで、これらの補正又は差は、「平均」心収縮期率と「平均」補正率を決定するために実験データがそこから収集された全ての他の被検体からの類似のデータに照らして統計的に分析される（ステップ３８０８）。心収縮期圧について研究された全ての患者にわたってのＡ線と圧力測定法による測定との間の平均誤りは−１３．６％であり、かつ平均圧力についての誤りは−５．６％であったことが、数百の患者の橈骨動脈のサンプリングに基づいて、本願の讓受人によって実験的に決定されている。要するに、これらの数字は、終末心拡張期値を先に説明したように基準値に整合させたように圧力測定法のトランスジューサを圧平させたとき、基準デバイス（すなわち、Ａ線）が心収縮期読取りについて約１４％だけ公称的に高く、かつ平均読取りについて約６％高かったことを表示する。
【０２３０】
ステップ３８０８によって注目のパラメータに対する補正率が得られた後、ステップ３８１０によって本発明の非侵襲的血行力学査定装置が評価される被検体の注目の血管、例えば、橈骨動脈に当てられる。次に、ステップ３８１２で、ストレスが血管に印加される。特に、ストレスは、注目の「定着」変量（本実施の形態では終末心拡張期圧）が、一般に、定着変量の実値に対して或る既知の関係を帯びる所望値を最も達成する点へと血管内の血行力学性質に変化を誘発するために、印加される。本文脈の中で、ストレスのこの印加は、本明細書に先に説明した非侵襲性（圧力測定法の）圧平デバイスを使用して血管を圧平する一方、血管の圧力波形を非侵襲的に測定することを含み、圧平の程度は測定された圧力波形中の定着変量のその実値に対する関係によって決定される。
【０２３１】
次に、ステップ３８１４で、拡大縮小関数又は伝達関数が評価される被検体についての前述のデータから決定される。１実施の形態では、未補正の圧力の入力値と「補正済み」圧力の入力値に基づいて、曲線適合多項式に関する次の３つの同時方程式を解くことによって得られる。すなわち、
【０２３２】
【数１５】
2004521660.tif 000016
【０２３３】
ここに、
Ｐｓｃ、Ｐｄｃ、Ｐｓｃ＝それぞれ、心収縮期、心拡張期、平均「補正済み」（拡大縮小）圧及び
Ｐｓｕ、Ｐｄｕ、Ｐｍｕ＝それぞれ、心収縮期、心拡張期、平均「未補正」（原）圧である。
【０２３４】
定数ｋ１、ｋ２、ｋ３について解くことによって、多項式曲線が発生され、これは、都合良く、圧力の全ての値に対して現被検体について圧力測定法によって得られた波形を拡大縮小する。したがって、これは、入力波形が終末心拡張期値に又は実値と既知の関係を持ついずれかの終末心拡張期点に（本明細書に先に説明した方法又は他の技術によって）整合させられるならば、拡大縮小された出力を用意する。
図３９は、前述のプロセスのグラフ表示である。
【０２３５】
血管位置と位置決め
先に論じたように、装置の超音波トランスジューサと圧力トランスジューサの適正な横位置決めは、動脈血圧を含む血行力学パラメータの測定の強化された正確性の面で有利である。このような横位置きめは、一般に、注目の血管の局所内に装置を初期的に配置する際に断定され、その後、最適結合が達成されるというように、トランスジューサの位置を「微同調」させる。人間の橈骨動脈の場合、動脈の近接位置は、医務者によって知られている、すなわち、手頸の幅よりもいくらか狭いに相当する狭い帯内を、被検体の手頸／前腕の内側表面を縦に下へ発している。しかしながら、他の血管及び／又は他の種については、注目の血管の位置と配向は、容易には識別されないことがある。その近接位置と配向が周知である人間の橈骨動脈の文脈の中でさえも、医務者による装置の不適性な配置は、その位置から得られた血圧測定の一貫性と「堅牢性（ｒｏｂｕｓｔｎｅｓｓ）」に影響するおそれがある。患者（及び他のファクタ）による動きもまた読取りの正確性に影響することがあり、測定装置の周期的再配置／再位置決めを必要とする。
【０２３６】
更になお、例えば、頸動脈（ｃａｒｏｔｉｄ  ａｒｔｅｒｙ）外科（人間の頸に沿って血管をマップする）、大腿動脈（ｆｅｍｏｒａｌ  ａｒｔｅｒｙ）のカテーテル法、及び時間的動脈局所化（ｔｅｍｐｏｒａｌ  ａｒｔｅｒｙ  ｌｏｃａｌｉｚａｔｉｏｎ）を含む他形式の手順があり、そこでは（ｉ）非侵襲的仕方で被検体の組織内の血管を容易にかつ信頼性を以て位置割出しするか、又は（ｉｉ）いったん位置割出しした血管に対して所定位置と配向に監視デバイス又は治療デバイスを維持するかのどちらかができることが大いに望ましい。
【０２３７】
したがって、本発明は、都合良く、注目の血管を位置割出しすること及び１つ以上の選択された装置を血管と所定の関係に維持することの両方に対して非侵襲性技術（と装置）を用意する。特に、１模範的応用では、後方散乱音響エネルギーが分析されて、被検体の、周囲の組織内に埋め込まれている血管を初期的に位置割出しする。他の応用では、後方散乱エネルギーがセンシング装置又は治療装置（例えば、本明細書に先に説明したＮＩＢＰデバイスの圧力トランスジューサ及び／超音波トランスジューサ）を血管に対して最適位置に維持するように使用される。
【０２３８】
１実施の形態では、注目の血管を検出しかつ位置割出しする方法がＡモード線走査に出現した後方散乱超音波エネルギーの低下レベルに基づいて血管と関連した管腔を直接検出することによって完遂される。超音波エネルギーは、周囲の血管壁及び組織によってよりも（血液内に存在する赤血球及び血しようによるというように）管腔内により容易に吸収される。それゆえ、低下後方散乱エネルギーの領域を識別することによって血管が位置割出しされ、及び／又は監視デイスが位置決めされる。このような領域は、一次元（例えば、組織内の深さに関してのみ）、又は血管を横位置つまりトラバース位置及び深さ位置の両方に関して精確に位置割出しすることが望ましい場合のように、多次元で見付けることもできる。
【０２３９】
本発明の他の実施の形態は、後方散乱Ａモード・エネルギーの分析に基づいて血管の前壁と後壁の両方の相対位置を検出する。この壁情報は、血管の有効直径を直接決定するために使用される。次いで、横位置つまりトラバース位置の関数としての血管直径の変動が測定デバイス又は治療デバイスの最適横位置を識別するために使用される。
【０２４０】
図４０を参照して、被検体の組織内の血管を検出する方法の第１実施の形態を詳細に説明する。図４０に示したように、方法４０００は、まず音響（例えば、超音波）エネルギーを発生しかつ所望ビーム幅で一般に血管が位置割出しされている注目の領域内で被検体の組織内へ伝送する（ステップ４００２）。例えば、橈骨動脈の場合、波を伝送する超音波トランスジューサを、超音波で以て血管の少なくとも部分を横断して切るように走査パターン内の注目の領域を横断して物理的に又は電子的に「ステア（ｓｅｅｒ）」することができる。例えば、皮膚の表面にわたって超音波トランスジューサを動かす（本明細書に先に説明したトラバース位置決め装置ような）組立て体が使用されることがある。これに代わり、大きな分散ビームが走査関数（ｓｃａｎ  ｆｕｎｃｔｉｏｎ）の少なくとも部分を取り除くために使用されることがある。ビーム分散は、その後、所望ならば、より良い信号品質のために注目の血管又は領域の識別の際に調節されることがある。なお他の代替実施の形態として、フェーズド・アレー又は他の電子的にステア可能なデバイスがトランスジューサ素子の物理的運動を伴わず組織を走査する比較的狭い分散の伝送／受信ビームを発生するために使用されることがある。
【０２４１】
検査下にある組織は、その組織を「前負荷（ｐｒｅ−ｌｏａｄ）」するために先に説明したように圧平されることがまたある。このような前負荷は、静脈やその他の小さな血管が超音波トランスジューサと注目の血管との間に介在している、又は、人間の橈骨動脈の場合のように、注目の血管の付近に全体的に介在している場合に、有利である。特に、制御された前負荷は小さい静脈を動脈に比較してその小さいフープ・ストレスと低い管腔（血液）圧のゆえに先取的に凹ませ、それによって、下に説明する管腔検出技術の視野からそれを有効に除去する。本願の讓受人は、３０ｍｍＨｇ（４ｋＰａ）の程度での前負荷圧平がトランスジューサと橈骨動脈との間に介在したいずれの静脈をも充分に凹ますために適していると決定している。もっとも、納得するように、特定応用次第で他の値も使用される。
【０２４２】
しかしながら、納得するように、注目の主血管の検出／位置割出しと同時に上に挙げた静脈のような他の血管を検出しかつ位置割出しをするのがまた望ましいことがある。それゆえ、前述の圧平を、管腔検出装置によって眺められたものを「フィルタ」するために所望するに従って選択的に使用することがある。例えば、所与の個人についてＡモード信号に高度のクラッタ及び他の雑音が存在するならば、動脈の回りの静脈からのエネルギー後方散乱の影響を除去するために評価される領域を圧平するのが望ましい。これに代わり、雑音レベルが充分に低いならば、静脈の存在を検出することが望ましく、このような検出は非侵襲的血圧測定と一緒に使用される組織拡大縮小関数又は伝達関数の決定に潜在的に有用である。特に、介在した静脈の存在は、伝達関数の値にいくらか影響することがあり、それでからその存在を検出してかつその影響を補償することが望ましい。再び図４０に帰って、走査される領域内の組織、血管、及び管腔から反射された後方散乱エネルギーが、次に、ステップ４００４で同じ（又は他の）トランスジューサ素子によって受信され、かつステップ４００６で装置内の処理によってＡモード又は他の所望モード性に変換される。後の血管検出を強化するために、ステップ４００８で、「包絡線平方」値の導出を含む（本明細書で先に詳細に説明した）追加信号処理がＡモード線に適用される。最後に、ステップ４０１０で、注目の領域内で血管が存在している位置を識別するために、Ａモード信号が所定計量又は他の判定基準の面で更に分析される。１変形では、計量は、注目の血管内管腔の存在から生じる正規化され積分されたパワー・プロファイル中のプラトーを識別するために適合している。それゆえ、積分された後方散乱パワー中に重要なプラトーがトラバース「掃引」の間中に識別されると、プラトーが検出されるときのトランスジューサの横位置は、血管の横位置に相当する。方法４０００の他の変形では、包絡線平方パラメータの信号レベルは、管腔の存在からまた生じる局所最小を識別するために、また分析される。
【０２４３】
第２実施の形態では、図４０のステップ４０１０での管腔の検出は、本明細書に先に「間隔」法と称した、Ａモード線に沿う所定距離（例えば、１ｍｍ）内の平均パワーを計算することによって完遂される。
【０２４４】
図４０の手順は、いずれかの数の空間的に分布した走査に必要なだけ繰り返されることがあり、このようにして、組織の「スライス」が注目の血管の横位置並びにまたその深さに関する情報を発生するために順序に（又は同時にさえも）深さ（時間）次元で検査されることがあることにまた注意されたい。このような「スライス」は、電子的に溶解又は組み合わされて、断層撮影法で周知であるように測定された包絡線平方値（又は先の実施の形態に関して説明した積分されたパワー）の多次元表現を用意する。本発明は、決定の正確性を向上するために多数のこのような技術の使用を更に企図する。例えば、「プラトー」検出と「間隔」検出の技術は、器械の堅牢性を増すために組み合わせられる。
【０２４５】
更になお、納得するように、管腔（血管）位置情報は、超音波トランスジューサ又は他の測定装置／治療装置を血管に対してどんな所望の配向にも、例えば、固定横オフセットして、（下に説明するように）血管の直径に機能的に関係している横オフセットして等、維持するために使用されることがある。同様に、管腔検出情報は、信号品質、嚢胞性構成要素（ｃｙｓｔｉｃ  ｃｏｍｐｏｎｅｎｔ）又はクラッタの存在、圧平の程度（もしあれば）、等々に基づいてトランスジューサの位置を制御するために使用されることがある。
【０２４６】
図４１を参照して、被検体の組織内の血管を検出しかつ位置割出しする方法の代替実施の形態を説明する。この代替実施の形態では、管腔検出の前述の技術が注目の血管の壁の位置の分析に対する開始点として使用される。それらの壁の位置は、血管の直径の間接的決定を与え、これがトランスジューサと血管との間に最適結合を維持するのを助けるために使用される。
【０２４７】
図４１に示したように、血管に対して測定デバイス又は治療デバイスを位置決めするために壁直径を利用する方法の１模範的実施の形態を説明する。注意するのは、次の議論は人間の橈骨動脈に対して超音波トランスジューサを再位置決めするために使用される（本明細書に先に説明したような）横位置決め組立て体の面で行われるが、その一般的方法論は他の応用にもたやすく適合させられるといえることである。
【０２４８】
方法４１００は、本明細書に先に説明したように、時間の関数として変動する「生」壁直径情報を導出すること（ステップ４１０２）を一般に含む。次に、ステップ４１０４で、ステップ４１０２の壁直径情報が、事前選択された有限時間間隔にわたって積分されて、生データを「平滑化」しかつ呼吸と心周期に関するほとんどの又は全てのアーチファクトを除去する。これは、時間の関数として変動する壁直径計量を生じる（ステップ４１０６）。次いで、壁直径計量は、絶対計量の値と機能的に関係している誤り信号を発生する（ステップ４１０８）ために使用される。次いで、この誤り信号は前述のコントローラ・デバイスに送られ、このデバイスは誤り信号の大きさを可能な限り小さく（又はこれに代えて直径計量を可能な限り大きく）有効に維持するようにトランスジューサ素子の横位置を変動させる（ステップ４１１０）。
【０２４９】
納得するように、上に挙げた壁直径情報は、先に本明細書にその内容を組み入れてある讓受人の同時係属米国特許出願第０９／３４２，５４９号に説明したようなファジィ論理制御含む、多くの異なった制御方式の基礎として使用されることがある。このような変形と応用は、コントローラ技術における当業者に周知であり、したがって、本明細書には更に説明しない。
【０２５０】
最後に注意するのは、上に説明した方法の多くの変形が本発明と両立して利用されるといえることである。特に、或る種のステップは、オプショナルでありかつ所望するに従って遂行されることも削除されることもある。同様に、（例えば、追加データのサンプリングや、処理や、フィルタ動作（ｆｉｌｔｅｒａｔｉｏｎ）や、校正や、数学的分析のような）他のステップが前述の実施の形態に追加されることがある。更に、或る種のステップの遂行の順序は、所望ならば、入れ替えられ、又は並列（又は直列）に遂行される。それゆえ、前述の実施の形態は、本明細書に開示した本発明の広範ないくつもの方法の単に例示である。
【０２５１】
上の詳細な説明は種々の実施の形態に応用された本発明の新規な特徴を図示し、説明し、かつ指摘してきたが、いうまでもなく、例示したデバイス又はプロセスの形式と細部における種々の省略、置換、及び変更は、本発明の精神に反することなく当業者によって行われることがある。前述の説明は、本発明を実施するのに現在思考される最良のモードのものである。この説明は、限定であることを決して意味せず、むしろ本発明の一般原理の例示であると取られるべきである。本発明の範囲は、特許請求の範囲を参照して決定されるべきである。
【図面の簡単な説明】
【図１】
先行技術「最大拍動」理論に従う印加圧力と時間の関数としての動脈の断面形状を例示する複合グラフ図である。
【図２】
血管の壁に狭窄症の変動レベルを形成された血管を例示する図である。
【図３〜３ｅ】
本発明に従って生被検体の循環系内の血行力学パラメータを査定する方法の１模範的実施の形態を例示する論理流れ図である。
【図４ａと４ｂ】
血液速度と血管の有効断面流体面積の減少との間の関係を例示するグラフ図である。
【図５ａ〜図５ｃ】
本発明に従う圧平圧力と、心洞調律と、動脈壁変位との間の関係を例示するグラフ図である。
【図５ｄ】
心拡張期圧と心収縮期圧の両方にについて最高血液速度と流体面積（圧平圧力）のパーセンテージ減少との間の関係を例示するグラフ図である。
【図６】
動脈と、皮膚と、介在した組織及び肉体的構成要素との間の関係を例示する典型的人間手頸の部分の断面図である。
【図７ａ〜７ｂ】
心収縮期状態と心拡張期状態の両方について、図６の断面に対する、それぞれ、模範的伝達関数と伝達分数を例示するグラフ図である。
【図８】
本発明を利用して典型的人間についての測定された、実の、及び選択された（校正された）動脈圧対時間のグラフである。
【図９】
本発明に従って生被検体の循環系の血行力学反応をモデルする方法の１模範的実施の形態を例示する論理流れ図である。
【図１０】
本発明に従って呼吸又は他の周期的誤り源に対して血行力学パラメータ測定を校正する方法の１模範的実施の形態を例示する論理流れ図である。
【図１１】
呼吸の血管への影響を例示する、典型的血管内を流れる血液の速度と運動エネルギ−（ＫＥ）のプロット図である。
【図１２】
本発明に従って生被検体の血管内の血行力学性質を測定する装置の１模範的実施の形態のブロック図である。
【図１３】
人間の橈骨動脈上での図１２の装置の使用を例示するその装置の第２実施の形態の機能ブロック図である。
【図１４ａ〜図１４ｂ】
本発明の方法を組み込んだコンピュータ・プログラムのオブジェクト・コード表現を有する種々の機械読取り可能媒体の斜視図である。
【図１５】
本発明に従ってパラメータ・データを分析する装置の第１実施の形態のブロック図である。
【図１６】
前掲の方法を使用して被検体に治療を施す方法の１模範的実施の形態を例示する論理流れ図である。
【図１７】
本発明に従う管腔と壁検出の一般的方法の１模範的実施の形態を例示する論理ブロック図である。
【図１７ａ】
図１７の方法によるＡモード信号を処理する方法の１模範的実施の形態を例示する論理ブロック図である。
【図１８】
監視される被検体の組織内深さの関数としての「生」Ａモード信号のグラフ図である。
【図１９】
図１８の生Ａモード・データから導出された包絡線のグラフ図である。
【図２０】
本発明に従う直角復調とＡモード信号フィルタリングの一般的方法論を例示する論理流れ図である。
【図２１】
本発明に従うＡモード信号の直角復調の方法を例示するブロック図である。
【図２２】
本発明の低域通過有限インパルス応答（ＦＩＲ）フィルタの１模範的実施の形態の周波数応答のグラフ図である。
【図２３】
本発明の直角復調とフィルタリング方法に使用される（ＦＩＲ）係数の模範的集合をリストするテーブルである。
【図２４】
図１８のＡモード信号から導出された包絡線平方波形のグラフ図である。
【図２５】
図２４の包絡線平方波形の（１０を底とする）対数のグラフ図である。
【図２６】
図２４のＡモード信号の包絡線平方への利得関数（ＴＧＣ）の適用のグラフ図である。
【図２７ａ】
本発明に従う管腔検出の「プラトー」法を例示する論理流れ図である。
【図２７ｂ】
単一Ａモード線について深さの関数として測定された後方散乱パワーを例示するプロット図である。
【図２７ｃ】
図２７ａのプラトー法に従う正規化され積分されたパワー関数のグラフ図である。
【２７ｄ】
図２７ａの方法と一緒に使用されたプラトー検出計量のグラフ図である。
【２７ｅ】
単一Ａモード線に対する正規化されたプラトー検出計量のグラフ図である。
【図２８ａ】
図２６のＴＧＣ補正されたＡモード線に沿うパワー・プロファイル（１ｍｍ間隔）のグラフ図である。
【図２８ｂ】
本発明に従ってＡモード線信号を使用して管腔検出の「間隔」方法論を例示する論理流れ図である。
【図２９ａ】
積分されたパワー計算に基づく前壁と後壁検出の第１模範的方法論を例示する論理流れ図である。
【図２９ｂ】
包絡線平方信号レベル決定に基づく前壁と後壁検出の第２模範的方法論を例示する論理流れ図である。
【図３０】
図２９ｂの方法に従う平均パワー計算と典型的Ａモード包絡線に現れる前壁と後壁アーチファクトとの間の関係のグラフ図である。
【図３１】
図２９ｂの方法を使用する（５０００Ａモード信号線に基づく）時間の関数としての模範的血管の前壁と後壁の位置のグラフ図である。
【図３２】
外部圧平デバイスによる圧平の間中の時間の関数としての血管の後壁の位置のグラフ図である。
【図３２ａ】
典型的橈骨動脈の圧平の間中時間の関数としての圧平デバイスの線形位置（ｍｍ）のグラフ図である。
【３２ｂ】
図３２ａに従う間中に時間の関数として圧平圧力トランスジューサによって測定された圧力プロファイルのグラフ図である。
【図３３】
圧平の関数としての、前及び後動脈壁位置（深さｍｍ）、及び動脈直径（ｍｍ）のグラフ図である。
【図３４】
本発明に従う圧平の関数としての血管直径の変化の速度を決定する方法論のグラフ図である。
【図３５】
本発明に従って管腔検出と壁検出を使用して血行力学パラメータを査定する方法論の１実施の形態を例示する論理流れ図である。
【図３５ａ】
終末心拡張期圧中の「フラット・スポット」を表示する、測定された動脈圧対圧平圧力のグラフである。
【図３６】
本発明に従って血行力学測定の間中に最適圧平圧力を決定する１模範的方法論を例示する論理流れ図である。
【図３６ａ】
本発明に従う模範的終末心拡張期「フラット・スポット」検出アルゴリズムのヒストグラム・プロット図である。
【図３７】
心拡張期圧プロファイル中の屈曲点に基づいて後方散乱超音波エネルギーを使用して血行力学パラメータを査定する方法論の１模範的実施の形態を例示する論理流れ図である。
【図３８】
組織損失を償うために測定された圧力を拡大縮小する方法の１模範的実施の形態を例示する論理流れ図である。
【図３９】
図３８の方法のグラフ図である。
【図４０】
本発明に従う管腔検出を使用して血管を位置割出しする方法を例示する論理流れ図である。
【図４１】
本発明に従う管腔検出と壁検出を使用して血管を位置割出しする方法の１模範的実施の形態を例示する論理流れ図である。
【図１】
2004521660.tif 000017
【図２】
2004521660.tif 000018
【図３】
2004521660.tif 000019
【図３ａ】
2004521660.tif 000020
【図３ｂ】
2004521660.tif 000021
【図３ｃ】
2004521660.tif 000022
【図３ｄ】
2004521660.tif 000023
【図３ｅ】
2004521660.tif 000024
【図４】
2004521660.tif 000025
【図５】
2004521660.tif 000026
【図６】
2004521660.tif 000027
【図７ａ】
2004521660.tif 000028
【図７ｂ】
2004521660.tif 000029
【図８】
2004521660.tif 000030
【図９】
2004521660.tif 000031
【図１０】
2004521660.tif 000032
【図１１】
2004521660.tif 000033
【図１２】
2004521660.tif 000034
【図１３】
2004521660.tif 000035
【図１４ａ】
2004521660.tif 000036
【図１４ｂ】
2004521660.tif 000037
【図１５】
2004521660.tif 000038
【図１６】
2004521660.tif 000039
【図１７】
2004521660.tif 000040
【図１７ａ】
2004521660.tif 000041
【図１８】
2004521660.tif 000042
【図１９】
2004521660.tif 000043
【図２０】
2004521660.tif 000044
【図２１】
2004521660.tif 000045
【図２２】
2004521660.tif 000046
【図２３】
2004521660.tif 000047
【図２４】
2004521660.tif 000048
【図２５】
2004521660.tif 000049
【図２６】
2004521660.tif 000050
【図２７ａ】
2004521660.tif 000051
【図２７ｂ】
2004521660.tif 000052
【図２７ｃ】
2004521660.tif 000053
【図２７ｄ】
2004521660.tif 000054
【図２７ｅ】
2004521660.tif 000055
【図２８ａ】
2004521660.tif 000056
【図２８ｂ】
2004521660.tif 000057
【図２９ａ】
2004521660.tif 000058
【図２９ｂ】
2004521660.tif 000059
【図３０】
2004521660.tif 000060
【図３１】
2004521660.tif 000061
【図３２】
2004521660.tif 000062
【図３２ａ】
2004521660.tif 000063
【図３２ｂ】
2004521660.tif 000064
【図３３】
2004521660.tif 000065
【図３４】
2004521660.tif 000066
【図３５】
2004521660.tif 000067
【図３５ａ】
2004521660.tif 000068
【図３６】
2004521660.tif 000069
【図３７】
2004521660.tif 000070
【図３８】
2004521660.tif 000071
【図３９】
2004521660.tif 000072
【図４０】
2004521660.tif 000073
【図４１】
2004521660.tif 000074
【国際公開パンフレット】
2004521660.tif 000075
2004521660.tif 000076
2004521660.tif 000077
2004521660.tif 000078
2004521660.tif 000079
2004521660.tif 000080
2004521660.tif 000081
2004521660.tif 000082
2004521660.tif 000083
2004521660.tif 000084
2004521660.tif 000085
2004521660.tif 000086
2004521660.tif 000087
2004521660.tif 000088
2004521660.tif 000089
2004521660.tif 000090
2004521660.tif 000091
2004521660.tif 000092
2004521660.tif 000093
2004521660.tif 000094
2004521660.tif 000095
2004521660.tif 000096
2004521660.tif 000097
2004521660.tif 000098
2004521660.tif 000099
2004521660.tif 000100
2004521660.tif 000101
2004521660.tif 000102
2004521660.tif 000103
2004521660.tif 000104
2004521660.tif 000105
2004521660.tif 000106
2004521660.tif 000107
2004521660.tif 000108
2004521660.tif 000109
2004521660.tif 000110
2004521660.tif 000111
2004521660.tif 000112
2004521660.tif 000113
2004521660.tif 000114
2004521660.tif 000115
2004521660.tif 000116
2004521660.tif 000117
2004521660.tif 000118
2004521660.tif 000119
2004521660.tif 000120
2004521660.tif 000121
2004521660.tif 000122
2004521660.tif 000123
2004521660.tif 000124
2004521660.tif 000125
2004521660.tif 000126
2004521660.tif 000127
2004521660.tif 000128
2004521660.tif 000129
2004521660.tif 000130
2004521660.tif 000131
2004521660.tif 000132
2004521660.tif 000133
2004521660.tif 000134
2004521660.tif 000135
2004521660.tif 000136
2004521660.tif 000137
2004521660.tif 000138
2004521660.tif 000139
2004521660.tif 000140
2004521660.tif 000141
2004521660.tif 000142
2004521660.tif 000143
2004521660.tif 000144
2004521660.tif 000145
2004521660.tif 000146
2004521660.tif 000147
2004521660.tif 000148
2004521660.tif 000149
2004521660.tif 000150
2004521660.tif 000151
2004521660.tif 000152
2004521660.tif 000153
2004521660.tif 000154
2004521660.tif 000155
2004521660.tif 000156
2004521660.tif 000157
2004521660.tif 000158
2004521660.tif 000159
2004521660.tif 000160
2004521660.tif 000161
2004521660.tif 000162
2004521660.tif 000163
2004521660.tif 000164
2004521660.tif 000165
2004521660.tif 000166
2004521660.tif 000167
2004521660.tif 000168
2004521660.tif 000169
2004521660.tif 000170
2004521660.tif 000171
2004521660.tif 000172
2004521660.tif 000173
2004521660.tif 000174
2004521660.tif 000175
2004521660.tif 000176
2004521660.tif 000177
2004521660.tif 000178
2004521660.tif 000179
2004521660.tif 000180
2004521660.tif 000181
2004521660.tif 000182
2004521660.tif 000183
2004521660.tif 000184
2004521660.tif 000185
2004521660.tif 000186
2004521660.tif 000187
2004521660.tif 000188
2004521660.tif 000189
2004521660.tif 000190
2004521660.tif 000191
2004521660.tif 000192
2004521660.tif 000193
2004521660.tif 000194
2004521660.tif 000195
2004521660.tif 000196
2004521660.tif 000197
2004521660.tif 000198
2004521660.tif 000199
2004521660.tif 000200
2004521660.tif 000201
2004521660.tif 000202
2004521660.tif 000203
2004521660.tif 000204
2004521660.tif 000205
2004521660.tif 000206
2004521660.tif 000207
2004521660.tif 000208
2004521660.tif 000209
【国際公開パンフレット（コレクトバージョン）】
2004521660.tif 000210
2004521660.tif 000211
2004521660.tif 000212