(57)【要約】  本発明は、最高血圧（ＳＢＰ）、最低血圧（ＤＢＰ）、および平均血圧（ＭＡＰ）の推定のためのシステムに関する。このシステムは、脈波について、そのエネルギと組み合わせた生理学的モデルを確立し、その後に以前のモデルの値ならびにユーザに関する他の変数（例えば、年齢、性別、身長、体重など）を含む固定長のベクトルを生成する。この固定長のベクトルが、対象の３つの変数の計算のための「ランダムフォレスト」に基づく関数推定器システムの入力として使用される。この関数推定器の主たる利点は、推定すべき関数についていかなる制約も事前に課すことがなく、本発明の場合のような異種のデータにおいてきわめて信頼できる点にある。
2012505679.tif 000002
【特許請求の範囲】
【請求項１】
血圧の非侵襲的測定のためのシステムおよび装置であって、
入力信号についての自己回帰移動平均式（ＡＲＭＡ）の循環機能の確率論的モデル（５）、（６）、および（７）を含み、
入力信号（８）および（９）のＴｅａｇｅｒ−Ｋａｉｓｅｒ演算子についての自己回帰移動平均式（ＡＲＭＡ）の血圧パルスの確率論的モデルを含み、
臨床データ（性別、年齢、身長、ボディ・マス・インデックスなど）およびその関数を含み、
「ランダムフォレスト」（３）に基づく関数推定システム（３）を含んでいる
ことを特徴とする、血圧の非侵襲的測定のためのシステムおよび装置。
【請求項２】
関数推定システム（３）の入力が、患者の以前のモデルおよび情報（性別、年齢、身長、体重、ボディ・マス・インデックス、パルスのリズム、心臓のコヒーレンス、前処理された入力信号のゼロ通過、および前処理された入力信号のゼロ通過の変動性など）を含む固定サイズのベクトルからなることを特徴とする、請求項２に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項３】
システムの入力信号が、光学的、機械的、または音響的に得られて前処理されたプレチスモグラフ波であることを特徴とする、請求項１または２に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項４】
プレチスモグラフ波が、デジタル・パルス・オキシメータによって得られていることを特徴とする、請求項１〜３のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項５】
（３）において推定される関数が、基本パラメータ（ＳＢＰ、ＤＢＰ、およびＭＡＰ）および（３）における推定誤差を減らすための前記パラメータの線形結合であることを特徴とする、請求項１〜４のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項６】
推定されたパラメータ（ＳＢＰ、ＤＢＰ、およびＭＡＰ）の系統誤差および分散を小さくするためにシステム（３）の推定の平均を実行する事後処理システム（４）を含むことを特徴とする、請求項１〜５のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項７】
ＳＢＰ、ＤＢＰ、およびＭＡＰの推定システムが、ＤＳＰ装置によって実現され、例えばＦＰＧＡ形式のマイクロコントローラによって実現されることを特徴とする、請求項１〜６のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項８】
推定されたパラメータ（ＳＢＰ、ＤＢＰ、ＭＡＰ）の分散を少なくする目的のＣＰＵによるデータ処理システムを内部に備えている少なくとも音響、機械、および／または光学カテーテルを含んでいる手動装置を備えていることを特徴とする、請求項１に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項９】
前記ＣＰＵが、ＤＳＰ、ＦＰＧＡ、またはマイクロコントローラ装置によって実現されることを特徴とする、請求項８に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項１０】
保存用メモリを備え、例えばフラッシュ式のメモリを備えていることを特徴とする、請求項６〜８のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項１１】
例えばシリアルポート、ブルートゥース、またはＵＳＢ（１２）によるＰＣへの外部接続手段、及び／又は例えばＷｉＦｉ、Ｚｉｇｂｅｅ、またはＵＷＢによるネットワーク接続手段を備えていることを特徴とする、請求項６〜９のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項１２】
データ視覚化画面を備えていることを特徴とする、請求項６〜１０のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【請求項１３】
制御スイッチ、電池、および／または外部の電源への接続を備えていることを特徴とする、請求項６〜１１のいずれか一項に記載の血圧の非侵襲的測定のためのシステムおよび装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、オシロメトリック法およびコロトコフ法にかかわらず、最高血圧、最低血圧、および平均血圧を非侵襲的に測定するためのシステムに関する。そのようにするために、圧力パルスおよびその瞬間のエネルギの生理学の確率論的モデルが作成され、「ランダムフォレスト」に基づいて関数を近似するためのシステムに組み合わせられる。入力信号は、他の患者の変数に組み合わせられたプレチスモグラフパルスの前処理済みのバージョンである。
【背景技術】
【０００２】
血液の循環の主たる機能は、組織のニーズを満たす（すなわち、組織へと栄養素を運び、廃棄物を運び去り、ホルモンを運び、組織のすべての液体の正しい平衡を維持する）ことである。
【０００３】
組織のニーズに関連した血流束の制御と、血流束に必要な心臓の制御および血圧との間の関係は、きわめて理解が困難であり、例えば評論家または高血圧の患者などといった血流力学的制御を必要とする患者の血液力学的管理について、この分野に多数の文献が存在し、多数の特許も存在している。
【０００４】
心臓が、連続的な方法で大動脈へと血液を送るため、この動脈の圧力は高い（平均で１００ｍｍＨｇ）。心拍はパルスに基づいているため、平均として、動脈の圧力は、１２０ｍｍＨｇの最高血圧（ＳＢＰ）と８０ｍｍＨｇの最低血圧（ＤＢＰ）との間を変動する。血液が体循環を通って流れるとき、その平均圧力（ＭＡＰ）は、心臓の右心耳の大静脈の終わりに達する瞬間の約０ｍｍＨｇまで、次第に減少する。
【０００５】
全身の毛細血管の圧力は、動脈側の極端の付近の３５ｍｍＨｇから、静脈側の極端の付近の１０ｍｍＨｇなどの低いレベルまで、さまざまであるが、血管床の大部分における平均の機能圧は約１７ｍｍＨｇであり、組織の細胞への栄養素の拡散を可能にしつつ少量の血漿が多孔質の毛細血管を通過することができる充分に低い圧力である。
【０００６】
肺循環に関しては、圧力が、やはり大動脈のようにパルス状であるが、最大血圧が２５ｍｍＨｇであって、最低血圧が８ｍｍＨｇであり、肺の平均の動脈圧は１６ｍｍＨｇにすぎない。肺の毛細管圧は、７ｍｍＨｇにすぎない。しかしながら、肺を通過する毎分の総血流束は、体循環と同じである。肺の系におけるこれらの低い圧力は、毛細血管がガスの交換のためだけに血液へと暴露されればよく、心臓へと戻るまでの血液の移動距離が短いため、肺のニーズに適している。これにもとづき、呼吸の機能およびガスの交換が、患者の血液力学、したがって血圧において、きわめて重要であると結論付けることができる。
【０００７】
具体的には、循環機能に３つの基本的な原理が存在する。
【０００８】
１．組織が活動しているとき、待機のときよりも多くの血流束が必要であるため、体のすべての組織の血流束が、組織のニーズに基づいて制御される（すなわち、代謝機能）。
【０００９】
２．心臓の支出（ＣＥ）の制御が、局所の組織の流束のすべての合計によって制御される（すなわち、静脈の戻りの応答が、そのような静脈の戻りの由来である動脈への心臓による送り返しである）。この意味で、心臓は、組織のニーズに応答して状態機械として機能する。しかしながら、心臓の応答は完璧ではなく、必要な量の血液を送らせる特別な神経信号を必要とする。
【００１０】
３．血液は、局所の血流束の制御および／または心臓の支出の制御によって、独立に制御される。例えば、圧力が通常の平均値（１００ｍｍＨｇ）を下回って低下するとき、この圧力を通常の値へと回復させるために、一連の循環の変化をもたらす反射インパルスのカスケードが生成される。神経信号が、心臓の送出圧力を高めるとともに、動脈樹に蓄積するような方法で、全身の細動脈の大部分の全体的収縮によってより多くの血液を心臓に与えるため、大きな静脈リザーバの収縮を強める。
【００１１】
新たな心拍の各々において、血液の新たな波が動脈を満たす。動脈系の膨張性ゆえに、血液は、心臓の収縮期および拡張期の両方において流れる。また、通常の状態においては、動脈樹の容量が、パルスの圧力を、血液が毛細血管に達するときにはほぼ消失して、組織におけるほぼ連続的な血流束（振動がきわめてわずかである）が保証されるように、低下させる。図１が、大動脈の根元のカテーテルによって侵襲的な形態で得られた圧力パルスの典型的な記録を示している。若い健常な成人において、心臓の収縮期における圧力（ＳＢＰ）は、約１２０ｍｍＨｇである一方で、拡張期における圧力（ＤＢＰ）は、約８０ｍｍＨｇである。両方の圧力の間の差が、パルス圧力（ＰＰ）と呼ばれ、通常の状態においては、約４０ｍｍＨｇである。
【００１２】
２つの主たる要因が、パルス圧力を左右する。
１．心臓の収縮期の容積。
２．動脈樹の総膨張性（容量）。
【００１３】
一般に、収縮期の容積が大きいということは、各々の心拍において動脈樹を満たさなければならない血液の量が多く、収縮期および拡張期における圧力の上昇および低下がより大きいことを意味し、より大きなＰＰにつながる。
【００１４】
他方で、ＰＰを、収縮期容積と動脈樹の容量との間の割合として定めることもできる。これら２つの要素のいずれかに影響する循環のあらゆるプロセスが、ＰＰにも影響を及ぼす。
【００１５】
上述にもとづき、心臓が収縮期に大動脈へと血液を送り出すとき、送出の開始においては、血液が慣性ゆえに遠方の血管へと速やかに移動することができないため、動脈の近位側の部分だけが膨張する。しかしながら、中央の大動脈の圧力の上昇が速やかに慣性を上回り、膨張の波面が大動脈のすべてを通って広がる。この現象は、動脈における圧力パルスの伝達として知られている。波面が動脈樹を通って広がるとき、圧力パルスの輪郭は、遠方の血管への伝達の細に減衰させられる（圧力パルスの吸収）。
【００１６】
人間の動脈圧を測定するために、主要な血液の通り道において侵襲的なカテーテルを使用することは、上述のように、重大な事例（例えば、敗血性ショックまたは多臓器不全の患者など、血液力学的に深刻な危険にさらされている患者）を除き、妥当でない。代わりに、動脈圧の測定／監視は、聴診法（コロトコフ音）などの間接的な方法によって行われる。この方法においては、聴診器が肘前の動脈に配置され、ホースが、上腕部に沿った動脈圧で膨張させられる。ホースが、動脈が血液によって膨張したままであるような低い圧力で腕を圧縮するときは、血液が動脈に沿って循環しているにもかかわらず、聴診器で音が聞こえることはない。しかしながら、ホースの圧力が、動脈圧のサイクルの一部において動脈を閉塞させるように充分に高いとき、各々のパルスにおいて音が聞こえる。したがって、この方法においては、最初にホースの圧力が、この圧力が上腕の動脈のＳＢＰよりも高いときに音が聞こえないような方法で、ＳＢＰを上回って高められる。この瞬間において、ホースの圧力の低下が開始され、圧力がＳＢＰ圧力を下回って低下する瞬間に、血液が動脈を通って流れ始め、コロトコフ音が心拍に同期して聞こえ始める。この瞬間に、ＳＢＰが割り出される。ホースの圧力が下がり続けるときに、リズミカルおよび粗い音によってコロトコフの質が変化する。ホースの圧力がＤＢＰに等しくなる瞬間に、音が聞こえなくなり、したがって圧力が割り出される。このシステムは、医学界において基準として考えられ、これについて、この方法によって動脈圧を割り出すためのシステムおよび電子装置を説明する文献および特許が存在している。
【００１７】
圧力波の伝播によって引き起こされる流体柱（通常は水銀）の振動の測定に基づく補完的なシステムも存在している。そのようなシステムは、オシロメトリック法として知られている。この方法を説明する文献、およびこの方法によって動脈圧を割り出すためのシステムおよび装置を記載するいくつかの特許も存在している。これらの方法を、動脈圧の測定を改善するために組み合わせることができる。要約すると、上述の方法は、物理的な形態の圧力を比較する機械的な方法に基づいている。
【００１８】
圧力パルスの伝達は、フォトプレチスモグラフのパルス（ＰＰＧ）に密接に関係している。なぜならば、これらの装置は、通常は近赤外（ＮＩＲ）の波長について、血液のヘモグロビンの光の吸収の変化を測定し、得られる信号が圧力パルスに比例するからである。実際、ＰＰＧを、患者の皮膚上で非侵襲な形態で使用される微小血管レベル（通常は、指または耳たぶ）の血液の体積の変化を測定するための低コストの技法と考えることができる。この技法は、例えばデジタル・パルス・オキシメータおよび血管診断システムなどの市販の医療装置において実現されている（例えば、ＰＰＧで不整脈または期外収縮を信頼できる方法で検出することができる）。
【００１９】
動脈の圧力を間接的に評価するために、いくつかの手段によって得られるＰＰＧの使用に言及するいくつかの特許が存在している。
【００２０】
心臓の拡張期の開始におけるＰＰＧ信号の識別からのＳＢＰの連続的監視のためのシステムが、特許出願ＵＳ１９７４０５２３１９６に確立されている。
【００２１】
循環系の分析モデルおよび決定論によって表わされる血液の体積、心臓の支出などからＳＢＰおよびＤＢＰを推定するためのシステムが、特許出願ＵＳ４４１８７００に記載されている。このシステムは、発明において発見された数学的モデルに反映される各患者についての特有の較正を必要とする（Ｋ定数）。
【００２２】
心電計信号（ＥＣＧ）のピークおよびパルス検出器によって検出されるパルスの間の時間経過に基づくＤＢＰの連続的監視のための方法が、特許出願ＵＳ４０３０４８５に記載されている。この発明は、パルスの伝送時間が動脈の圧力につれて変化するという原理に基づいている。ＳＢＰの測定システムが、最初に血圧計（機械的な方法）によって各々の患者について特定的に較正される。
【００２３】
ＰＰＧ信号に基づくＳＢＰおよびＤＢＰの連続的な監視のための方法が、特許出願ＵＳ５１４０９９０に記載されている。ＳＢＰおよびＤＢＰが、ＰＰＧによって得られる血液の体積、ならびに平均値の前に割り出される患者の動脈圧−血液の体積に関する一定のパラメータＫを使用する各患者に特有の較正期間におけるＳＢＰおよびＤＢＰの指標から割り出される。
【００２４】
耳たぶにおけるＰＰＧ信号のパルスの経過時間から中間動脈圧（ＭＡＰ）を連続的に測定するための方法が、特許出願ＵＳ５２３７９９７に記載されている。ＳＢＰおよびＤＢＰが、耳たぶの血液の体積の密度の測定から得られる。この発明は、従来からの手段（オシロメトリックまたはコロトコフ）による個々の動脈の張力の値の較正を必要とする。
【００２５】
ＥＣＧおよびＰＰＧ信号からＳＢＰおよびＤＢＰを割り出すための方法が、特許出願ＵＳ５８６５７５５およびＵＳ５８５７９７５に記載されている。動脈圧が、各パルスにおけるパルスの到着時間、体積測定の波の形状、および心臓のリズムから得られる。これらの特許は、動脈圧を割り出すために、ＰＰＧパルスの開始と振幅の５０％との間の時間の差によって、ＥＣＧのＲ波とＰＰＧパルスの開始との間の時間差を使用する。
【００２６】
パルスの経過時間と、少なくとも従来からのシステムおよび線形回帰分析による較正の後の心臓のリズムおよびパルスの面積とから、動脈の張力を測定するためのシステムおよび装置が特許出願に記載されている。
【００２７】
体の別々の部分に置かれた２つのＰＰＧセンサに基づいてシステムＮＩＢＰを確立し、これらのセンサによる平均のパルスの経過時間の間の鎖を計算し、心臓によって送られる血液の体積の変化を割り出すための技法が、欧州特許出願ＥＰ０４４３２６７Ａ１に記載されている。このシステムは、従来からのシステムによる較正を必要とする。
【００２８】
ＰＰＧから種々の生理学的パラメータを得るためのシステムが、特許ＵＳ２００４／０２６０１８６Ａ１に記載されている。より具体的には、この特許は、呼吸のリズム、心臓のリズム、心臓のリズムの変動性、血液の体積の変動性、および自律神経系についての情報の推定、ならびに動脈圧の相対変化（絶対変化ではない）の監視を実行する。
【００２９】
フォトプレチスモグラフィ（ＰＰＧ）の原理に基づいてホースを使用せずに動脈の圧力を測定するためのシステムが、特許ＵＳ２００６／００７４３２２Ａ１に記載されている。この特許は、ホースなしで動脈の圧力を測定するためのシステムおよび装置を特許請求しているが、オシロメトリックおよびコロトコフの原理に基づく各ユーザについての較正を必要とする。ひとたびシステムが較正されると、これをＰＰＧの原理によって単独に個人化された方法で使用することができる。
【００３０】
ＰＰＧ波の基本周波数を得、それらから血圧を得るために、心臓血管の波（ＰＰＧから得られる圧力パルス）の形状の調和分析によって動脈樹において見られる血液の体積を得るためのシステムおよび装置が、特許ＵＳ２００７／００３２７３２Ａ１に記載されている。
【００３１】
動脈の圧力を非侵襲な方法で監視（ＮＩＢＰ）するためのシステムであって、膨張可能なホースと、このホースの初期の膨張圧力を割り出すためにパルスオキシメータ（ＳｐＯ２）にて実現されたフォトプレチスモグラフとを備えているシステムが、特許ＵＳ２００８／００４５８４６Ａ１に記載されている。この発明の基本原理は、ＳＢＰよりも高いホースの膨張圧力からＰＰＧ波のパルスが消えることにある。この方法で、ユーザが、ホースの過膨張から保護され、他のより安全な測定を実行することができる。
【００３２】
上述の特許と対照的に、血圧の測定時間を短縮および／または最適化するためにＰＰＧ信号を使用するＮＩＢＰシステムが、特許ＵＳ２００８／００８２００６Ａ１に記載されている。この発明においては、ホースの収縮期間がＰＰＧ信号によって制御される。
【発明の概要】
【発明が解決しようとする課題】
【００３３】
ＰＰＧ信号を個々の較正または血圧の従来からのシステムのグループとの組み合わせにおいて基本的な機能の原理として使用する多数の特許が見られるが、可動の機械部品を含まず、非侵襲で動脈の圧力を臨床的に割り出す（ＮＩＢＰ）ことができ、血圧計による個人の較正を必要とせず、正しい性能および動作のために従来からの手段（あるいは、上述したオシロメトリックまたはコロトコフ法の発展）の支援を必要とせずに機能することができる、より安全で信頼できる連続監視システムを発見するというニーズを、依然として解決すべきである。
【課題を解決するための手段】
【００３４】
すでに述べた本発明の従来技術によれば、パルスが、血圧に依存する減衰および形態の変化を被る。この影響は、ＳＢＰおよびＤＢＰの間の差に依存して変化する。本発明の提案のシステムおよび装置は、パルス（ＰＰＧ）の形状および圧力レベルの間の関数関係の干渉にもとづき、情報が、パルスおよびその統計と患者の血圧状態との間の依存関係から推定される。
【００３５】
ＳＢＰ、ＤＢＰ、およびＭＡＰの推定を実行するための入力情報が、関数推定部の仕事を容易にするために処理される。ＰＰＧ信号は変化する継続時間を有するため、各々の測定について固定長のベクトルを生成するための処理が実行される。このベクトルは、パルスに関する情報（自己相関係数および移動平均）、パルス間の平均距離、その分散、瞬間のエネルギに関する情報、エネルギの変動性、ならびに例えば患者の性別、年齢、体重、身長、臨床情報（ボディ・マス・インデックスまたは同様の測定値）などといった個人の臨床情報を含む。
【００３６】
関数の干渉のためのシステムは、パルスおよび血圧レベルの間の関係にいかなる関数の制約も課されないという意味で、盲目的に働く。ＰＰＧを血圧レベルに関連付ける関数形式が未知であるため、この関数を推測するためのシステムであって、臨床情報およびＰＰＧの波形から導出されるパラメータなどの無関係な入力変数の前で信頼できるシステムが選択される。また、技法は、本発明の従来技術において述べたように、他のパラメータにも関係している。本発明の関数の推定のための好ましいシステムは、例えば決定および回帰ツリー（ＣＡＲＴ）、スプライン、分類器コミッティ、サポート・ベクトル・マシン、およびニューラルネットワークなどの他の「機械学習」システムおよびパターン認識と比べ、「ランダムフォレスト」である。ランダムフォレストは、各ノードにおけるランダム変数の選択によって関数を推定する複数の決定ツリーの並列生成にもとづき、ノードの刈り込みは実行されず、各々のツリーが異なる系統一般化誤差を呈するような方法で学習データベースの無作為な部分集合によって学習される。したがって、各々のツリーの推定の平均を実行するとき、系統誤差が補償され、推定の分散が減少する。
【００３７】
本発明の実行は、２つの異なる段階を含んでいる。第１の段階は、１回だけ実行され、したがって後の較正／個人化を必要としないシステムの学習である。この段階は、性別、体重、年齢、などといった患者の種々のパラメータについての情報と、プレチスモグラフ波の記録とを有するデータベースを得ることからなる。この情報が、決定ツリーのパラメータの推定に使用され、システムに保存される。
【００３８】
第２の段階は、学習の段階において得られた一式のツリーの情報を取り込み、測定の瞬間において患者のプレチスモグラフ波を例えば性別、体重、年齢、などの他の変数とともに記録することからなる。この段階において、システムは、プレチスモグラフパルスの情報を読み取り、これらの処理を実行し、信号を記述する情報を有する固定長のベクトルを生成する。個人に関するさらなる情報がこのベクトルに追加され、対象の変数のいくつかの中間関数を計算する一式の「ランダムフォレスト」が適用される。その後に、対象の変数が、これらの中間関数から計算される。
【図面の簡単な説明】
【００３９】
【図１】本発明の図１は、侵襲カテーテルによって得られたパルス波の形状を示している。
【図２】本発明の図２は、説明されるシステムおよび装置の全体的なブロック図を示している。
【図３】本発明の図３は、デジタル・パルス・オキシメータによって得られたプレチスモグラフ波の形状を示している。
【図４】図４は、本発明において説明される前処理システムの詳細なブロック図を示している。
【図５】図５は、本発明において説明される圧力パルスの生理学の確率論的モデルの確立のためのＡＲフィルタの詳細なブロック図を示している。
【発明を実施するための形態】
【００４０】
本発明は、血圧（最高血圧、最低血圧、および平均血圧）の連続的な監視のためのシステム（図２）からなり、血圧のデータが、プレチスモグラフ信号取得装置（１）（光学的、音響的、または機械的な信号）によって推定され、本発明の好ましい実施例は、パルス・オキシメトリ・システム（ＳｐＯ２）からなる。ＰＰＧ情報が、例えば年齢、性別、身長、体重、などといった個人の他のデータと組み合わせられ、従来技術に提示の循環系の生理学による確率論的モデルを実行するデジタル前処理システム（２）に結び付けられる。このシステムが、圧力パルスの伝達、したがって血圧をランダムな方法で左右するすべてのパラメータをより良好な方法で捕らえる。
【００４１】
得られた確率論的モデルのベクトルが、「ランダムフォレスト」に基づいて関数を近似するデジタルシステム（３）にも結び付けられるが、その主たる機能は、基本パラメータ（ＳＢＰ、ＤＢＰ、およびＭＡＰ）をこれらに関係する種々の別の関数によって推定し、事後処理の工程（４）における推定の誤差を小さくすることにある。システム（４）の主たる機能は、得られるＳＢＰ、ＤＢＰ、およびＭＡＰの推定の系統誤差（バイアス）およびばらつきを小さくすべく、先の工程（３）の関数の平均によってＳＢＰ、ＤＢＰ、およびＭＡＰの最終的な値を推定することにある。システム（２、３、および４）は、ＦＰＧＡ、ＤＳＰ、またはマイクロコントローラなどといったいくつかの装置を備えるＣＰＵによって実現される。
【００４２】
ＰＰＧ曲線を得るためのシステム（１）が、組織の微小血管網の体積の変化を検出するための非侵襲、低コスト、かつ簡単な技法を実現する。このシステムの最も基本的な実施例は、以下を含む少数の光電子部品しか必要としない。
【００４３】
１．組織（例えば、皮膚）を照明するための１つ以上の光源。
２．検出体積内の組織の注入の変化に関係する光強度の小さな変動を測定するための１つ以上のフォトディテクタ。
【００４４】
ＰＰＧは、通常は非侵襲な方法で使用され、赤外または近赤外（ＮＩＲ）の波長で動作する。ＰＰＧにおいて最も認識される波形は、末梢脈拍（図３）であり、各々の心拍に同期している。ＰＰＧによって得られる波動と侵襲カテーテルによって得られるパルスとの間の類似性（図１および３）を認めることが重要である。ＰＰＧによって得られる情報は、きわめて価値のある情報であるため、本発明の主たる入力として考えられる。
【００４５】
ＰＰＧ波は、各々の心拍に同期した血液の体積の変化に関連した生理学的なパルス波（ＡＣ成分）を含んでいる。この成分が、呼吸のリズム、中枢神経系の活動、および体温調節に関係する他の基本的な低周波成分（ＤＣ成分）に重畳している。ＡＣ成分の基本周波数が、心臓のリズムに応じて約１Ｈｚに発見される（図３）。
【００４６】
光と生物学的組織との間の相互作用は複雑であり、散乱、吸収、反射、透過、および蛍光などの光学プロセスを含んでいる。システム（１）について選択される波長が、以下の理由できわめて重要である。
【００４７】
１．水の光学ウインドウ：組織の主たる成分は水である。これは、紫外の波長および赤外帯域内の長い波長をきわめて吸収する。可視光（赤色）またはＮＩＲが組織を通過することを可能にし、これらの波長において血流束またはその体積の測定を可能にするウインドウが、水の吸収スペクトルに存在する。したがって、本発明は、システム（１）のためにＮＩＲの波長を使用する。
【００４８】
２．等吸収の波長：この波長を除き、酸素ヘモグロビン（ＨｂＯ２）と還元ヘモグロビン（Ｈｂ）との間に、吸収の大きな差が存在する。したがって、この波長（すなわち、ＮＩＲ範囲における８０５ｎｍ付近）において実行される測定では、信号が組織の酸素飽和度の変化の影響を受けないと考えられる。
【００４９】
３．組織への進入：所定の放射強度における組織への光の進入の深さは、選択される波長の関数でもある。ＰＰＧにおいては、進入体積が、（１）において使用されるシステムなどの透過システムにおいて約１ｃｍ３である。
【００５０】
ＰＰＧパルス（図３）は、２つの異なる段階を示している。すなわち、パルスの立ち上がりを呈する上昇段階、およびパルスの低下を呈する下降段階である。第１の段階が、心臓の収縮期に関係している一方で、第２の段階は、心臓の拡張期ならびに循環系の末梢において波が被る反射に関係している。下降段階における重拍脈は、動脈硬化でない健康な患者のＰＰＧにおいても一般的に見られる。
【００５１】
本発明の従来技術において説明したとおり、循環樹における圧力パルスＰＰの伝播も考慮に入れなければならない。ＰＰは、循環樹の末梢に向かって移動するにつれてその形状を変化させ、増幅され／減衰させられ、形状および時間的特性の変化を被る。これらの変化は、末梢において動脈が細くなることで引き起こされるＰＰの反射に起因する。さらに、ＰＰパルスの伝播は、位相ひずみの周波数依存によっても影響される。
【００５２】
パルスを生成する生理学的プロセスのために、ＡＲＭＡモデル（自己回帰移動平均モデル）が生成の機構を特徴付け、したがってこのモデルが、ＰＰの表現として考えられている。非線形の相互作用を並行してモデル化するために、Ｔｅａｇｅｒ−Ｋａｉｓｅｒ演算子が使用され、ＡＲ（自己回帰）システム（２）に組み合わせられる。
【００５３】
図１および３に見られるとおり、ＰＰはＰＰＧに類似しており、同様の変化が、血管病理学において観察される（狭窄によって引き起こされるクッション効果または拍動性の変化）。
【００５４】
システムのパルスオキシメータ（１）は、患者の動脈における酸素飽和度（ＳｐＯ２）についての情報を得るためにＰＰＧを使用する。すでに述べたように、ＳｐＯ２を、赤色またはＮＩＲの波長で組織（通常は、指または耳たぶ）を照明することによって得ることができる。通常は、ＳｐＯ２装置は、パラメータを割り出すために両方の波長の間の交換を使用する。両方の波長の大きさが、これらの波長におけるＨｂＯ２およびＨｂの間の吸収の差ゆえに、ＳｐＯ２の変化に影響される。ＳｐＯ２を、大きさの間の比、ＰＰＧ、ならびにＡＣおよびＤＣ成分から得ることができる。
【００５５】
パルスオキシメトリにおいては、組織を透過する光の強度（Ｔ）が、一般にＤＣ信号として知られており、組織の光学特性（すなわち、吸収係数および散乱係数）の関数である。動脈の拍動が、酸素および還元ヘモグロビンの濃度の周期的な変動を引き起こし、さらには吸収係数の周期的な変動をもたらす。
【００５６】
ＰＰＧのＡＣ成分の強度の変動は、以下のようであってもよい。
2012505679.tif 000003
【００５７】
この生理学的な波動の形状は、光の強度の変動に比例し、光の強度そのものは、散乱および吸収係数（それぞれ、
2012505679.tif 000004
）の関数である。
2012505679.tif 000005
の変動を、酸素および還元ヘモグロビンの変動（
2012505679.tif 000006
）の線形な変動として定めることができる:
2012505679.tif 000007
εｏｘおよびεｄｅｏｘは、酸素および還元ヘモグロビンの遮光効率（すなわち、所定の環境において距離単位による散乱および吸収ゆえに失われる光の割合）である。これまでの式にもとづき、動脈の酸素飽和度（ＳｐＯ２）を、下記によって定めることができる。
2012505679.tif 000008
【００５８】
ＡＣ成分に依存したＳｐＯ２の表現を、選択された波長（赤色およびＮＩＲ）に式（Ｉ）および（ＩＩＩ）を直接適用することによって得ることができる。
2012505679.tif 000009
ここで、
2012505679.tif 000010
が成り立つ。
【００５９】
血液中の同期変化に関係しない低周波への影響を補償するためにＤＣ成分でＡＣ成分を正規化する（従来技術を参照）ことで、以下が得られる。
2012505679.tif 000011
【００６０】
このパラメータを（ＩＶ）に含めることで、以下が得られる。
2012505679.tif 000012
ここで、
2012505679.tif 000013
が成り立ち、
ここで、
2012505679.tif 000014
および
2012505679.tif 000015
は、ＲおよびＮＩＲの波長において推定された式（Ｉ）に相当する。
【００６１】
式（ＶＩ）が、ＳｐＯ２への正確な解であるが、ｋを推定することは、
2012505679.tif 000016
を有していないため不可能である。いずれにせよ、ｋおよびＲは、組織の光学的特性の関数であり、ｋをＲの関数として表現することが可能である。より詳しくは、ｋを以下の線形回帰として表わすことができる。
ｋ＝ａＲ＋ｂ    （ＶＩＩ）
（ＶＩＩ）
【００６２】
この線形回帰は、経験的に導出されるが、Ｐの強度において平坦な波（その吸収係数が、以下のように定められる）を仮定している較正係数を意味する。
ｄＰ＝μａＰｄｚ  （ＶＩＩＩ）
ここで、ｄＰは、吸収係数μａを有する一様な環境における無限小ｄｚを通過する光線の強度の微分変化を表わす。したがって、ｚについて積分を行うことによって、Ｂｅｅｒ−Ｌａｍｂｅｒｔの法則が得られる。
2012505679.tif 000017
2012505679.tif 000018
を仮定し、式（ＶＩＩ）はｋ＝１に換算され、これは、本発明において実行されるパルスオキシメトリの測定における好ましい近似である。
【００６３】
得られるシステム（１）のＰＰＧ信号が、本発明のシステム（２）（図４）の励起として使用され、システム（２）の主たる機能は、ＳＢＰ、ＤＢＰ、およびＭＡＰの推定のための循環機能の確率論的モデルを確立することにある。
【００６４】
本発明の先行技術において、圧力パルスＰＰの形状および伝播の両方において重要な役割を有する種々のパラメータを説明した。これらのパラメータは、心臓の支出、心臓のリズム、心臓の同期、呼吸のリズム、代謝機能などに関係している。ＰＰとＰＰＧとの間の密接な関係もすでに説明した。したがって、すでに詳述したパラメータがＰＰの伝播の種類において重要な役割を有するため、それらがＰＰＧも左右すると仮定される。
【００６５】
これを考慮に入れ、本発明の好ましい実施例は、ＡＲＭＡ（ｑ，ｐ）（ｑ（ＭＡ）およびｐ（ＡＲ）の近似のｑによる自己回帰移動平均モデル）の確率論的モデル化のシステム（５）を使用する。
【００６６】
時間的な連なりＰＰＧ（ｎ）、ＰＰＧ（ｎ−１）、・・・、ＰＰＧ（ｎ−Ｍ）を、階差における以下の式が満たされる場合に、ｐ＝Ｍ次のＡＲプロセスとしてモデル化することができる。
2012505679.tif 000019
ここで、係数
2012505679.tif 000020
は、ＡＲとして知られるパラメータであり、
2012505679.tif 000021
は、白色プロセスである。項
2012505679.tif 000022
が、係数
ａｋ
および
2012505679.tif 000023
の内部積であり、ここでｋ＝１、・・・、Ｍである。式（Ｘ）を、
2012505679.tif 000024
のように書き直すことができ、ｖｋ＝−ａｋである。
【００６７】
先の式から、パルス
2012505679.tif 000025
の実際の値が、先の値（
2012505679.tif 000026
）プラス予測誤差の項
2012505679.tif 000027
の有限な線形結合に等しいものとして割り出される。したがって、式（Ｘ）を線形畳み込みとして書き直すことで、以下が得られる。
2012505679.tif 000028
【００６８】
一般性を失うことなく、ａ０＝１を、以下によって定められる予測フィルタのＺ変換として定めることができる。
2012505679.tif 000029
【００６９】
2012505679.tif 000030
をＰＰＧパルスのＺ変換として定めることで、
2012505679.tif 000031
であり、ここで、
2012505679.tif 000032
である。
図５が、システム（１）において得られるパルス
2012505679.tif 000033
のＡＲ成分の分析フィルタを示している。
【００７０】
パルス
2012505679.tif 000034
のｑ＝Ｋ次のＭＡ（移動平均）に関して、これを、ガウス性白色雑音によって励起される不連続な線形パルスの応答として記述することができる。したがって、ＥＤＦとして記述されるフィルタのＭＡ応答は、
2012505679.tif 000035
であり、ここで［ｂ１、ｂ２、□、ｂＫ］は、ＭＡパラメータとして知られる定数であり、
2012505679.tif 000036
は、平均がゼロであって、分散がσ２である白色プロセスである。したがって、式（ＸＩＩ）および（ＸＶＩ）を関連させることで、
2012505679.tif 000037
が得られ、
2012505679.tif 000038
が、ＡＲＭＡ（ｑ，ｐ）モデルの誤差の項である。（ＸＶＩＩ）においてＺ変換を実行し、以下が得られる。
2012505679.tif 000039
ＡＲおよびＭＡベクトルの第１の項を、一般性を失うことなく１に等しくできるため、システム（２）におけるＡＲＭＡ（ｑ，ｐ）フィルタ表現（５）は、下記によって定められる。
2012505679.tif 000040
【００７１】
ここで、Ａ（ｚ）およびＢ（ｚ）は、それぞれの
2012505679.tif 000041
のＡＲおよびＭＡ成分である。本発明の好ましい実施形態は、ｑ＝１およびｐ＝５次のＡＲＭＡモデルを使用するが、［４，１２］の間を含む任意のｐおよびｑ次を使用することができる。
【００７２】
ひとたびＷｏｌｄの分解およびＬｅｖｉｎｓｏｎ−Ｄｕｒｂｉｎの回帰によってＡＲＭＡ（ｑ，ｐ）モデルが計算されると、Ｈ（ｚ）が生成され、入力信号がＨ（ｚ）の逆数によってフィルタ処理される（６）。さらに、残余の統計
2012505679.tif 000042
が、サブシステム（７）において計算される。得られたこれらのサブシステムの情報が、出力の固定長のベクトル
2012505679.tif 000043
に保存される。
【００７３】
本発明の前処理システム（２）は、すでに述べたものと同等のｐ次のＡＲプロセスによってＴｅａｇｅｒ−Ｋａｉｓｅｒ演算子を計算してその出力をモデル化するサブシステム（８）をさらに備えている。
【００７４】
この場合、一般性を失うことなく、ＰＰＧパルスは、一種のＡＭ−ＦＭ変調信号（振幅および周波数の両方において変調されている）と考えられる。
2012505679.tif 000044
2012505679.tif 000045
および
2012505679.tif 000046
は、ＰＰＧの瞬間の振幅および周波数である。決定された信号のＴｅａｇｅｒ−Ｋａｉｓｅｒ演算子は、
2012505679.tif 000047
によって定められ、
2012505679.tif 000048
である。
【００７５】
この演算子が、式（ＸＸ）のＡＭ−ＦＭ変調信号に適用され、ＰＰＧの振動を生み出している源の瞬間のエネルギがもたらされる。これが
2012505679.tif 000049
であり、ここで近似誤差は、ＰＰＧパルスの場合など、瞬間の振幅
2012505679.tif 000050
および瞬間の周波数
2012505679.tif 000051
が
2012505679.tif 000052
という平均値に比べて過度に速くは変化しない場合、大きくない。
【００７６】
2012505679.tif 000053
のｐ次のＡＲプロセスは、図５のものと同等のフィルタ（９）によって実現される。本発明の好ましい実施形態は、ｐ＝５次のＡＲモデルを使用するが、４〜１２の間のｐおよびｑ次において、任意の他のモデルを使用することができる。
【００７７】
ひとたびＡＲＭＡ（ｑ，ｐ）モデル（５、６、および７）に基づく確率論的モデルおよびＡＲＭＡ（ｑ，ｐ）モデルが、Ｔｅａｇｅｒ−Ｋａｉｓｅｒ（８および９）演算子によって計算されると、本発明は、サブシステム（１０）によってＰＰＧから心臓のリズム（ＨＲ）および心臓の同期（すなわち、心臓のリズムの変動性）を計算する。本発明の好ましい実施形態は、信号の自己相関関数によって、２秒〜５分の間でさまざまであってもよいＰＰＧの時間ウインドウにおいて、心臓のリズムを計算する。
【００７８】
前処理システム（２）は、ＰＰＧ信号のゼロ通過ならびにそれらゼロ通過の変動を計算するサブシステム（１１）をさらに備えている。本発明の好ましい実施形態は、２秒〜５分の間でさまざまであってもよいＰＰＧの時間ウインドウにおいて、心臓のリズムを計算する。
【００７９】
最後に、前処理システム（２）は、患者に関係する変数を生成するためのサブシステム（１２）を備えており、それらの変数は、とりわけ、以下の変数である。
１．性別、年齢、体重、身長、患者が何らかの食品を食べたか否か、一日のうちの時刻。
２．ボディ・マス・インデックス。
３．体重を年齢で割ったもの。
４．体重をＨＲで割ったもの。
５．身長をＨＲで割ったもの。
６．ＨＲを年齢で割ったもの。
７．身長を年齢で割ったもの。
８．年齢をボディ・マス・インデックスで割ったもの。
９．ＨＲをボディ・マス・インデックスで割ったもの。
【００８０】
システム（２）を構成するサブシステムにおいて得られたすべてのデータは、固定長の出力ベクトル
2012505679.tif 000054
に保存される。
【００８１】
ひとたび特徴の固定長のベクトル
2012505679.tif 000055
が得られると、ＳＢＰ、ＤＢＰ、およびＭＡＰの推定を、「ランダムフォレスト」に基づいて関数を近似するためのシステム（３）によって実行することができる。本発明において提示される関数推定システムは、ひとたび「ランダムフォレスト」が正しく学習されたならば、いかなる較正も必要としない。
【００８２】
より具体的には、「ランダムフォレスト」は、ツリー構造
2012505679.tif 000056
を有する一式の分類器からなる分類器であり、ここで、
2012505679.tif 000057
は、無作為の独立ベクトルであって、同一に分布しており（ｉ．ｉ．ｄ）、各々のベクトルが、入力Ｖの最も一般的な分離に投票を配置する。この近似は、ユニークツリーに基づく他の分類器と比べ、信頼性において明確な利点を有しており、パルスと血圧レベルとの間の関係にいかなる機能的制約も課さない。
【００８３】
本発明において使用される「ランダムフォレスト」は、指標
2012505679.tif 000058
が数値を有するような方法で乱数ベクトル
2012505679.tif 000059
に依存する決定樹の成長によって生成される。各々のツリーに関するこの乱数ベクトルが、各々のノードの無作為な分布をもたらすと同時に、各々のツリーのためにデータの異なる部分集合を結果として与える学習ベースの無作為なサンプリングについての情報も提供する。この結果にもとづき、本発明において使用される分類器の一般化誤差が、下記の式によって定められる。
2012505679.tif 000060
【００８４】
「ランダムフォレスト」の一般化誤差は、一意の決定樹による一般化誤差よりも小さいため、
2012505679.tif 000061
を定めることで、以下が得られる。
2012505679.tif 000062
【００８５】
各々のツリーが異なる一般化誤差を呈し、ρが、（ＸＸＩＶ）に定められる残余の間の相関を表わす。この事実は、残余（ＸＸＩＶ）の間の小さな関係が大きな推定をもたらすことを意味する。本発明においては、この最小の相関が、サブシステム（２）において学習されているツリーの各々のノードの特徴ベクトルの無作為なサンプリングプロセスから来る。一般化誤差をさらにもっと減らす目的において、本発明は、対象のパラメータ（ＳＢＰ、ＤＢＰ、およびＭＡＰ）ならびにそれらの線形結合の両方を推定する。
【００８６】
「ランダムフォレスト」は、各々のツリーに、ブートストラップシステムによる対称的な変動性の後で、系統誤差を導入するために変更されたＣＡＲＴ式（「分類および回帰ツリー」）の一式の決定樹からなる（両方のランダムプロセスが指標
2012505679.tif 000063
の分析においてパラメータ
2012505679.tif 000064
によってモデル化される）。各々の実施形態における種々の系統誤差が、２つの機構によって導入される。
１．属性の部分集合の各ノードにおける無作為な選択が、各ツリーが異なる方法で振る舞うような方法で類似のノードの間で異なるツリーにおいて形成される仕切りの統計レベルにおいて相当物の確立を不可能にする。
２．ツリーをそれらの最大値まで成長させる。この場合、ツリーが、ルールに基づく探索テーブルと同様の方法で機能する。属性のサンプリングゆえ、それらは異なる構造を有する探索テーブルである。
【００８７】
このプロセスの結果は、各ツリーが異なる統計誤差を呈することである。
【００８８】
さらに、これら２つの修正の各々において、各ツリーがブートストラップ式のサンプルで学習する（すなわち、サンプルが入力データから取られ、結果として入力データの一部が欠ける一方で、残りの部分が繰り返される）。このブートストラップの効果は、平均の推定を行うときに補償される変動性を導入する。
【００８９】
これらの特徴の大域的結果は、系統誤差および誤差の変動性を容易に補償でき、他の種類の関数推定器（ＸＸＶ）を上回る正確さをもたらすシステム（４）である。このシステムにおいて、基本分類器は、（本発明の場合のように）外れ値または異質な種類のデータを有する入力分布において何が強力であるかをレベルに基づいて決定するツリーである。
【００９０】
システム（４）の好ましい実施形態は、ノードレベル（２〜４７の間のレベルを選択することができる）の４７のうちの２つの要素および１００というブートストラップサイズ（サイズを２５〜５００の間で変更できる）の無作為なサンプルを得ることからなる。
【００９１】
本発明による手持ち装置は、データおよび装置の機能のための制御指示を視覚化するための画面を備えることができる。少なくとも音響、機械、および／または光プローブを備え、その信号が、ＤＳＰ、ＦＰＧＡ、またはマイクロコントローラによって実現されるＣＰＵによる事後処理システムによって解釈される。さらに、データおよびシステムの作動プロセスを保存するための作業メモリを備える。
【００９２】
さらに、本発明は、手動装置が、装置の作動および制御のための技術水準によるボタンまたはスイッチパネル、ならびに電池および／または外部の電源へのアクセスを備えることを想定する。
【００９３】
最後に、本発明によって得られた結果を、シリアルポートまたはＵＳＢあるいはネットワーク接続（例えば、ＷｉＦｉまたはブルートゥース）によって、分析のためにＰＣへ送信することができる。
【００９４】
些細な細部の変更による代案が、本明細書に記載のとおりの本発明の技術的範囲に包含されることを、理解されたい。
【図１】
2012505679.tif 000065
【図２】
2012505679.tif 000066
【図３】
2012505679.tif 000067
【図４】
2012505679.tif 000068
【図５】
2012505679.tif 000069
【国際調査報告】
2012505679.tif 000070
2012505679.tif 000071
2012505679.tif 000072
2012505679.tif 000073
2012505679.tif 000074
2012505679.tif 000075