(57)【要約】  本発明の実施形態は、末梢脈波形が波形を構成する別個の異なる位相にフィルタリングされる、改善された変換方法を提供する。別個の位相は、再結合されて推定された動脈内伝達関数を提供する前に変換される。たとえば、一実施形態では、末梢脈波形は第１ハイパスフィルタによってフィルタリングされ、末梢脈波形のコピーは、カットオフ周波数の異なる第２ハイパスフィルタによってフィルタリングされる。そして、２つのフィルタリング波形を、たとえば元の波形に再度加算することによってさらに処理することができ、その後、それら波形は、時分割式に合わせて多重化されて最終波形を提供する。たとえば、収縮期に対応する第１フィルタリング波形の部分を、拡張期に対応する第２フィルタリング波形の部分と結合して最終波形を生成することができ、心臓のそれぞれの時相の特徴を抽出するように、それぞれのフィルタカットオフ周波数を選択することができる。
2014534023.tif 000002
【特許請求の範囲】
【請求項１】
オシロメトリック測定から導出される末梢脈波形を処理して対応する動脈内波形を推定する方法であって、
前記末梢脈波形をそれぞれのフィルタを用いてフィルタリングして複数の波形を提供するステップと、
前記複数の波形を時分割多重化して、前記推定された動脈内波形を表す結合波形を提供するステップと、
を含む方法。
【請求項２】
前記末梢脈波形が、第１周波数でフィルタリングされて第１フィルタリング波形を提供し、第２周波数でフィルタリングされて第２フィルタリング波形を提供し、
前記第１周波数が心周期の第１時相に関連し、前記第２周波数が前記心周期の第２時相に関連する、請求項１に記載の方法。
【請求項３】
前記複数の波形のうちの少なくとも１つが、結合される前に変換される、請求項１または２に記載の方法。
【請求項４】
前記複数の波形のうちの少なくとも１つが、結合される前に増幅される、請求項２または３に記載の方法。
【請求項５】
前記第１時相が、前記心周期の収縮期である、請求項２〜４のいずれか一項に記載の方法。
【請求項６】
前記第１時相が、前記心周期の収縮期を開始する急な立上りである、請求項２〜４のいずれか一項に記載の方法。
【請求項７】
前記第２時相が、前記心周期の拡張期である、請求項２〜６のいずれか一項に記載の方法。
【請求項８】
前記第１周波数が、前記心周期の収縮期の最高周波数に対応する、請求項２〜７のいずれか一項に記載の方法。
【請求項９】
前記第１周波数が、３Ｈｚと２０Ｈｚとの間の周波数である、請求項２〜８のいずれか一項に記載の方法。
【請求項１０】
前記第１周波数がおよそ３．２Ｈｚである、請求項２〜９のいずれか一項に記載の方法。
【請求項１１】
第２周波数が、前記心周期の拡張期の最高周波数に対応する、請求項２〜１０のいずれか一項に記載の方法。
【請求項１２】
前記第２周波数が、１．５Ｈｚと２０Ｈｚとの間の周波数である、請求項２〜１１のいずれか一項に記載の方法。
【請求項１３】
前記第２周波数がおよそ１．７５Ｈｚである、請求項２〜１２のいずれか一項に記載の方法。
【請求項１４】
前記末梢脈波形が、第１ハイパスフィルタによってフィルタリングされて前記第１フィルタリング波形を提供し、
第２ハイパスフィルタによってフィルタリングされて前記第２フィルタリング波形を提供する、請求項２〜１３のいずれか一項に記載の方法。
【請求項１５】
前記第１ハイパスフィルタが、前記第２フィルタより高次である、請求項１４に記載の方法。
【請求項１６】
前記第１ハイパスフィルタが７次フィルタである、請求項１４または１５に記載の方法。
【請求項１７】
前記第２ハイパスフィルタが４次フィルタである、請求項１４〜１６のいずれか一項に記載の方法。
【請求項１８】
平滑化フィルタを用いて前記結合波形を平滑化するステップをさらに含む、請求項１〜１７のいずれか一項に記載の方法。
【請求項１９】
前記平滑化フィルタが２０Ｈｚを超える周波数を減衰させる、請求項１８に記載の方法。
【請求項２０】
前記平滑化フィルタが、カットオフ周波数がおよそ６Ｈｚであるローパスフィルタである、請求項１８または１９に記載の方法。
【請求項２１】
前記平滑化フィルタが、前記第１ハイパスフィルタおよび前記第２ハイパスフィルタより低次である、請求項１４に従属する請求項１８〜２０のいずれか一項に記載の方法。
【請求項２２】
前記平滑化フィルタが３次フィルタである、請求項１４に従属する請求項１８〜２１のいずれか一項に記載の方法。
【請求項２３】
前記結合波形が全心周期である、請求項１〜２２のいずれか一項に記載の方法。
【請求項２４】
前記末梢脈波形が対象のものであり、
前記対象の末梢収縮期血圧および拡張期血圧が得られ、
前記方法が、前記推定された動脈内波形、前記末梢収縮期血圧および前記拡張期血圧を使用して真の平均動脈内圧を計算するステップをさらに含む、請求項２３に記載の方法。
【請求項２５】
前記平均動脈内圧を使用して推定された中心圧力波形を計算するステップをさらに含む、請求項２４に記載の方法。
【請求項２６】
前記対象がヒトである、請求項２４または２５に記載の方法。
【請求項２７】
前記複数のフィルタのうちの少なくとも１つがバターワースフィルタである、請求項１〜２６のいずれか一項に記載の方法。
【請求項２８】
オシロメトリック測定から導出される末梢脈波形を処理して対応する動脈内波形を推定する装置であって、
前記末梢脈波形を受け取る入力と、
前記末梢脈波形を第１周波数でフィルタリングして第１フィルタリング波形を提供するように構成された第１フィルタと、
末梢脈波形を第２周波数でフィルタリングして第２フィルタリング波形を提供するように構成された第２フィルタと、
時分割マルチプレクサと、
前記推定された動脈内波形を出力する出力と、
を具備し、
前記時分割マルチプレクサが、前記第１フィルタリング波形および前記第２フィルタリング波形を結合して、前記推定された動脈内波形を生成する、装置。
【請求項２９】
前記第１周波数が３Ｈｚと２０Ｈｚとの間の周波数である、請求項２８に記載の装置。
【請求項３０】
前記第１周波数がおよそ３．２Ｈｚである、請求項２８に記載の装置。
【請求項３１】
前記第２周波数が１．５Ｈｚと２０Ｈｚとの間の周波数である、請求項２８に記載の装置。
【請求項３２】
前記第２周波数がおよそ１．７５Ｈｚである、請求項２８に記載の装置。
【請求項３３】
前記第１フィルタおよび前記第２フィルタがハイパスフィルタである、請求項２８〜３２のいずれか一項に記載の装置。
【請求項３４】
前記第１フィルタが前記第２フィルタより高次である、請求項２８〜３３のいずれか一項に記載の装置。
【請求項３５】
前記出力波形をフィルタリングする平滑化フィルタをさらに具備する、請求項２８〜３４のいずれか一項に記載の装置。
【請求項３６】
前記平滑化フィルタが、２０Ｈｚを超える周波数を減衰させる、請求項３５に記載の装置。
【請求項３７】
前記平滑化フィルタが、カットオフ周波数がおよそ６Ｈｚであるローパスフィルタである、請求項３５または３６に記載の装置。
【請求項３８】
前記推定された動脈内波形が、前記第１フィルタリング波形、前記第２フィルタリング波形およびフィルタリングされていない末梢脈波形の時分割多重を含む、請求項２８〜３７のいずれか一項に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、対象の血圧を測定し推定する方法および装置に関し、より詳細には、対象の動脈内血圧波形を推定する方法および装置に関する。
【背景技術】
【０００２】
血圧（ＢＰ：ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）は、循環血液が血管の壁に及ぼす力を指し、主な生命徴候のうちの１つを構成する。循環血液の圧力は、血液が動脈、細動脈、毛細血管および静脈を通って移動するに従って低下する。血圧という用語は、一般に、動脈圧、すなわち相対的に大きい動脈内の圧力を指す。動脈は、心臓から血液を運び出す血管である。
【０００３】
動脈内の血圧は、概して振動して変化し、波形（時間に対する圧力のグラフ）として表示することが可能である。動脈内のピーク圧力は、収縮期血圧（ＳＢＰ：ｓｙｓｔｏｌｉｃ  ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）として知られ、心周期の開始の近くで発生する。心周期の休止期に発生する、動脈内の最低圧力は、拡張期血圧（ＤＢＰ：ｄｉａｓｔｏｌｉｃ  ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）として知られている。心周期を通しての平均圧力は、平均動脈圧（ＭＡＰ：ｍｅａｎ  ａｒｔｅｒｉａｌ  ｐｒｅｓｓｕｒｅ）として知られており、脈圧（ＰＰ：ｐｕｌｓｅ  ｐｒｅｓｓｕｒｅ）は、収縮期圧と拡張期圧との差である。
【０００４】
既存の血圧モニタは、末梢動脈、たとえば腕における収縮期血圧および拡張期血圧の測定を行う。しかしながら、上腕動脈、橈骨動脈または指動脈で測定される収縮期血圧は、大動脈基部における中心収縮期血圧（ｃＳＢＰ：ｃｅｎｔｒａｌ  ｓｙｓｔｏｌｉｃ  ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）を超え、それは、収縮期血圧は、反射圧力波の結果として上肢の末梢動脈に沿った伝播により心臓の近くの中心収縮期血圧を超えて増幅されるためである、ということが長く認識されてきた。これにより、末梢脈圧と比較した中心脈圧の差ももたらされる。心拡張中の圧力の変化は低速であるため、拡張期血圧は、中心部位および末梢部位において同様である。平均動脈圧もまた、中心部位および末梢部位で同様である。ｃＳＢＰは、心臓に対する負荷をより適切に表すことができるので、末梢収縮期血圧（ｐＳＢＰ：ｐｅｒｉｐｈｅｒａｌ  ｓｙｓｔｏｌｉｃ  ｂｌｏｏｄ  ｐｒｅｓｓｕｒｅ）よりも心疾患により密接に関連することが予測される。
【０００５】
血圧は、通常、オシロメトリック法によって上腕において上腕の周囲で膨らませたカフを用いて測定される。カフ内の圧力を、腕内の収縮期血圧（末梢収縮期血圧、ｐＳＢＰ）を超える圧力まで膨張させ、その後、腕における拡張期血圧（ＤＢＰ）未満の圧力まで低速で収縮させる。収縮中のいずれの時点においても、１心周期中の平均カフ圧（ＭＣＰ）がｐＳＢＰ未満に低下すると、カフ圧は、ＭＣＰの周囲でわずかな量で振動する。ｐＳＢＰおよびＤＢＰを、カフ内の圧力振動の振幅から推定することができる。
【０００６】
平均動脈血圧（ＭＡＰ）およびＤＢＰは、大動脈と腕の導管動脈との間でほとんど差がないが、上肢で測定されるｐＳＢＰは、反射による増幅をもたらす上肢に沿った伝播により、大動脈の中心収縮期血圧（ｃＳＢＰ）を超えて増幅されることが知られている。ｃＳＢＰは、ｐＳＢＰより心血管イベントのリスクのより適切な推定を可能にすると考えられる。さらに、ｐＳＢＰには同様の影響を与えるがｃＳＢＰには特異な影響を与える降圧薬計画を比較した場合、低いｃＳＢＰほど、転帰の改善に関係した。したがって、ｃＳＢＰの非侵襲性推定方法が必要とされている。
【０００７】
最も一般的に採用される方法は、アプラネーショントノメトリ（ａｐｐｌａｎａｔｉｏｎ  ｔｏｎｏｍｅｔｒｙ）（橈骨動脈の上に圧力センサを保持してそれを下にある骨に対して緩やかに押し付ける）によって、橈骨動脈からの末梢血圧波形を測定することを含む。そして、この圧力波形を、ｐＳＢＰおよびＤＢＰのオシロメトリック法による測定値から校正することができる。そして、この末梢カフ圧波形に対して「一般化圧伝達関数」（ＧＴＦ：ｇｅｎｅｒａｌｉｓｅｄ  ｔｒａｎｓｆｅｒ  ｆｕｎｃｔｉｏｎ）を適用して、その末梢カフ圧波形を、ｃＳＢＰを推定することができる中心波形に変換する。ＧＴＦを、周波数領域において高速フーリエ変換を使用して、または時間領域においてパラメトリック関数を使用して導出することができる。ＧＴＦは、中心波形の所与の形状または周波数成分に対して、上腕が、年齢および他の個体内特性とは無関係に波形に対して比較的一定の影響を及ぼすという事実を利用する。アプラネーショントノメトリは、訓練された観察者によって行われる必要があり、行うために数分間かかり、比較的高価な機器を必要とし、オシロメトリック法または血圧の他の測定法が必要である。
【０００８】
トノメトリに基づく方法の進歩により、血圧のオシロメトリック測定で使用される上腕カフから血圧波形を直接導出するようになる。カフが（通常、ＤＢＰとｐＳＢＰとの間の圧力まで）膨張すると、カフから記録された圧力波形は、腕内の動脈内の圧力波形（およびトノメトリによって取得された圧力波形）に幾分か類似している。カフを、最高収縮期圧（ｓｕｐｒａｓｙｓｔｏｌｉｃ  ｐｒｅｓｓｕｒｅ）（ｐＳＢＰを超える）まで膨張させることができるが、これは、患者には不快であり、また、本方法が、血圧の定期的なオシロメトリック測定中に適用することがより困難であることも意味する。こうした測定中、目的は、ｐＳＢＰからＤＢＰへの収縮中、カフ圧振動に関する情報を最大化するために、カフ圧がｐＳＢＰを超える時間を最小限にすることである。
【０００９】
トノメータから得られた末梢波形を中心波形に変換するために使用されるＧＴＦは、カフ波形を中心波形に変換するために使用することはできない。しかしながら、代替的なＧＴＦを使用して、カフ波形を中心波形に変換することができる。ＧＴＦの特徴は、カフ波形の取得中のカフ内の平均圧力によって決まる。カフ波形は動脈内波形またはトノメータから得られた波形とは異なるため、ｐＳＢＰおよびＤＢＰから校正することができない。カフ波形にＧＴＦを適用することから得られる中心波形を、中心部位および末梢部位におけるＭＡＰおよびＤＢＰが等しいためにＭＡＰおよびＤＰＢから校正することができる。ＭＡＰ値およびＤＢＰ値を、オシロメトリック法によって得ることができる。
【００１０】
したがって、上記方法により、オシロメトリック血圧の従来の測定中に（または測定の直前あるいは直後に）、患者または観察者に対して従来の血圧の測定以上の面倒をかけない方法で、血圧カフからｃＳＢＰを求めることができる。この血圧カフ手法の第１の不都合は、ＭＡＰおよびＤＢＰから波形を校正する要件である。目下、オシロメトリック法によって導出されるようなＭＡＰの精度を確認するために合意された規格はない。第２の不都合は、ＭＡＰおよびＤＢＰからの校正は、オシロメトリック法によるＤＢＰの推定の誤差によってより著しく影響を受ける可能性がある、ということである。逆に、ｐＳＢＰおよびＤＢＰによる校正が行われる場合、ｐＳＢＰの推定の誤差により、ｃＳＢＰに同様の誤差が発生する。したがって、オシロメトリック法が幾分かの誤差を受け易い場合であっても、推定されたｃＳＢＰとｐＳＢＰとの間の差は、真のｃＳＢＰとｐＳＢＰとの間の差の正確な尺度であり続ける。
【００１１】
中心波形に変換する前の、ｐＳＢＰおよびＤＢＰによるカフから得られた波形の校正は、同様の手続きがトノメータから得られた波形に適用される場合に比較して相対的に不正確である。これは、トノメータから得られた波形とは異なり、カフ波形の形状および形態は、真の動脈内圧力波形のものと一定の関係を有していないためである。カフ波形は、ＭＣＰと動脈内圧力との差と、心周期の時相とによって決まる量によって歪められる。ｐＳ
ＢＰを超える圧力までカフを膨張させることによって動脈内圧力により近似させることができるが、これには上述した不都合がある。
【発明の概要】
【発明が解決しようとする課題】
【００１２】
オシロメトリック法を用いて測定される、すなわちカフ圧力装置を用いて得られる末梢脈波形を、対応する動脈内波形に変換する一方法は、一般的伝達関数を使用する。本方法は、完全には正確ではなく、したがって、末梢脈波形を動脈内波形に変換する方法の精度を向上させるという動機がある。
【課題を解決するための手段】
【００１３】
本発明の実施形態は、末梢脈波形が波形を構成する別個の異なる位相にフィルタリングされる、改善された変換方法を提供する。別個の位相は、再結合されて推定された動脈内伝達関数を提供する前に変換される。たとえば、一実施形態では、末梢脈波形は第１ハイパスフィルタによってフィルタリングされ、末梢脈波形のコピーは、カットオフ周波数の異なる第２ハイパスフィルタによってフィルタリングされる。そして、２つのフィルタリング波形を、たとえば元の波形に再度加算することによってさらに処理することができ、その後、それら波形は、時分割式に合わせて多重化されて最終波形を提供する。たとえば、収縮期に対応する第１フィルタリング波形の部分を、拡張期に対応する第２フィルタリング波形の部分と結合して、最終波形を生成することができ、心臓のそれぞれの時相の特徴を抽出するように、それぞれのフィルタカットオフ周波数を選択することができる。
【００１４】
本発明の第１態様によれば、オシロメトリック測定から導出される末梢脈波形を処理して対応する動脈内波形を推定する方法が提供され、本方法は、末梢脈波形をそれぞれのフィルタを用いてフィルタリングして複数の波形を提供するステップと、複数の波形を時分割多重化して、推定された動脈内波形を表す結合波形を提供するステップとを含む。
【００１５】
オシロメトリック測定を使用することにより、必要な技能レベルが低い単純な装置を用いてデータを収集することができ、それにより、装置は、訓練されていない人が使用するのに対してかつ自動患者家庭モニタリングに対して適したものとなる。時間領域において波形を多重化することにより、個々の波形を相互干渉なしに再結合することができる。
【００１６】
末梢脈波形を、第１周波数でフィルタリングして第１フィルタリング波形を提供し、第２周波数でフィルタリングして第２フィルタリング波形を提供することができ、第１周波数は心周期の第１時相に関連し、第２周波数は心周期の第２時相に関連する。
【００１７】
複数の波形のうちの少なくとも１つを、結合する前に変換し、任意選択的に、結合する前に増幅することができる。フィルタリング波形は、その対応する動脈内波形を最も適切に再構成するために増幅を必要とする場合がある。増幅は正の場合も負の場合もある。
【００１８】
第１時相は、心周期の収縮期、または心周期の収縮期を開始する急な立上りであり得る。第２時相は、心周期の拡張期であり得る。本発明の上記態様を使用して心周期の所定成分を隔離することにより、別個の成分を別個に変換して対応する動脈内波形を推定することができる。
【００１９】
第１周波数は、３Ｈｚと２０Ｈｚとの間の周波数、好ましくは３．２Ｈｚであり得る。これらの周波数は、心周期の収縮期に関連する。第２周波数は、１．５Ｈｚと２０Ｈｚとの間の周波数、好ましくはおよそ１．７５Ｈｚであり得る。これらの周波数は、心周期の拡張期に関連する。
【００２０】
末梢脈波形は、第１ハイパスフィルタによってフィルタリングされて第１フィルタリング波形を提供し、第２ハイパスフィルタによってフィルタリングされて第２フィルタリング波形を提供する。所定の周波数をフィルタリングするために個々のフィルタが使用される場合、その個々のフィルタをより容易に同調させることができる。
【００２１】
第１フィルタは第２フィルタより高次であり得る。この構成により、複数の周波数のフィルタリングが可能になるとともに、回路またはソフトウェアを提供する構成の複雑性も最小限になる。第１フィルタを７次フィルタとすることができ、第２フィルタを４次フィルタとすることができる。フィルタのこうした次数は、本発明を実施するために特に有利であることが分かった。
【００２２】
結合波形を、平滑化フィルタを用いて平滑化することができる。好ましくは、平滑化フィルタは、２０Ｈｚを超える周波数を減衰させる。平滑化フィルタは、カットオフ周波数がおよそ６Ｈｚであるローパスフィルタであり得る。平滑化フィルタは、第１フィルタリング波形および第２フィルタリング波形が結合された時に導入される、結合波形の望ましくない部分を除去する。最も重要な周波数は２０Ｈｚ未満であり、６Ｈｚが、フィルタに対して最適な応答を提供する好ましいカットオフ周波数であることが分かった。
【００２３】
平滑化フィルタは、第１フィルタおよび第２フィルタより低次であり得る。これにより、本発明のより単純な実施が可能になるとともに、推定された動脈内波形の精度が維持される。３次フィルタが、好適な周波数応答を提供することにより、平滑化フィルタとして有用であることが示された。
【００２４】
結合波形は全心周期であり得る。これにより、波形からより多くの情報を得ることができる。
【００２５】
動脈内心臓波形を使用して、正確な動脈内平均動脈血圧（ＭＡＰａｖｗａｖｅ）を推定することができる。末梢血圧測定を使用するカフ測定装置を使用して対象から取得されたオシロメトリック測定値のみを使用して、正確なＭＡＰａｖｗａｖｅを達成することができる。
【００２６】
ＭＡＰａｖｗａｖｅの推定には、中心波形を校正するためにｐＳＢＰおよびＤＢＰが必要である。これは、（本発明によって行われるように）カフ波形を動脈内波形に変換することによって可能である。そして、この動脈内波形を、ｐＳＢＰおよびＤＢＰから校正し、この動脈内波形の平均としてＭＡＰａｖｗａｖｅを得るために使用することができる。これは、ＭＡＰを推定するオシロメトリック法ではなく、オシロメトリック法によって得られる収縮期圧および拡張期圧に依存する。そして、元の末梢カフ波形を中心動脈波形に変換してＭＡＰａｖｗａｖｅおよびＤＢＰから校正することができる。
【００２７】
本発明の第２態様によれば、オシロメトリック測定から導出される末梢脈波形を処理して対応する動脈内波形を推定する装置が提供され、本装置は、末梢脈波形を受け取る入力と、末梢脈波形を第１周波数でフィルタリングして第１フィルタリング波形を提供するように構成された第１フィルタと、末梢脈波形を第２周波数でフィルタリングして第２フィルタリング波形を提供するように構成された第２フィルタと、時分割マルチプレクサと、推定された動脈内波形を出力する出力とを備え、時分割マルチプレクサは、第１フィルタリング波形および第２フィルタリング波形を結合して、推定された動脈内波形を生成する。
【００２８】
第１周波数は、３Ｈｚと２０Ｈｚとの間であり得るが、好ましくはおよそ３．２Ｈｚであり得る。第２周波数は１．５Ｈｚと２０Ｈｚとの間であり得るが、好ましくはおよそ１
．７５Ｈｚであり得る。第１フィルタおよび第２フィルタをハイパスフィルタとすることができ、第１フィルタを第２フィルタより高次とすることができる。本装置は、出力波形をフィルタリングする平滑化フィルタをさらに備えることができる。いかなる平滑化フィルタも、２０Ｈｚを超える周波数を減衰させることができ、特に、平滑化フィルタは、カットオフ周波数がおよそ６Ｈｚであるローパスフィルタであり得る。これらの特徴の利点は、上記対応する方法の特徴の利点に対応する。
【００２９】
時分割マルチプレクサはまた、フィルタリングされていない末梢脈波形を第１フィルタリング波形および第２フィルタリング波形と結合することができる。これにより、元の末梢脈波形の一部を使用して推定された動脈内波形のより適切な再構成が可能になる。
【００３０】
フィルタは、本発明の少なくとも１つの態様では、好適な位相および周波数応答を提供するためバターワースフィルタであるものとして明示的に挙げられている。本発明のいずれかの態様または実施形態のフィルタのうちの１つまたはすべてを置き換えるように、または本発明の実施形態のうちのいずれかのフィルタのうちのいずれかに加えて、他のタイプのフィルタを使用することができる。こうした他のタイプのフィルタとしては、パッシブフィルタおよびアクティブフィルタを含むアナログフィルタ、および電子部品を使用して実装されるデジタルフィルタまたはマイクロプロセッサを使用して実装されるソフトウェアフィルタが挙げられる。他のフィルタの例としては、チェビシェフ型フィルタ、ベッセル型フィルタおよび楕円型フィルタが挙げられる。フィルタの選択された次数は、特に好適であることがわかったが、フィルタの他の実数を使用することができる。
【００３１】
ここで、本発明の実施形態を、添付図面を参照して単に例として説明する。
【図面の簡単な説明】
【００３２】
【図１】本発明の一実施形態のブロック図を示す。
【図２】本発明の実施形態による信号を結合する方法を示す。
【図３】本発明の別の実施形態のブロック図を示す。
【図４】本発明のさらなる実施形態のブロック図を示す。
【図５】従来技術の方法の結果の線形単回帰を示す。
【図６】図５のデータセットの平均および標準偏差のグラフを示す。
【図７】従来技術の方法と比較して本発明の実施形態に従って実行される方法例の結果の線形単回帰グラフを示す。
【図８】図８のデータセットの平均および標準偏差のグラフを示す。
【発明を実施するための形態】
【００３３】
ここで、本発明の実施形態について説明する。
【００３４】
図１は、入力１として末梢脈波形を受け入れ、出力２として推定された対応する動脈内波形を提供する信号処理ブロック図を示す。入力１を受け入れる３つのフィルタＦ１、Ｆ２、Ｆｎが示されている。フィルタのうちの１つＦｎおよびその関連する構成要素Ａｎ、８は任意選択的であり、これは、それらが破線を用いて描かれていることによって示されている。各フィルタＦ１、Ｆ２、Ｆｎは、末梢脈波形の所定位相に関連する異なる周波数に同調される。各フィルタの出力信号ｆ１、ｆ２、ｆｎは増幅され、増幅されたフィルタ出力信号の各々は入力信号に加算される。増幅器Ａ１、Ａ２、Ａｎが、その関連するフィルタ出力ｆ１、ｆ２、ｆｎを、１を超える、１に等しいまたは１未満の増幅率で増幅することができる。増幅されたフィルタ出力は、構成要素６、７、８を用いて入力信号と結合される。結合構成要素の出力は、出力２において推定された動脈内波形を形成するように再構成される。出力信号は入力信号に対応するが、脈波形の所定の位相に対応する入力波形の各成分は、フィルタ処理および増幅プロセスによって変換される。
【００３５】
図１の任意選択的な構成要素および接続を、３つのフィルタリング経路を提供するように含めることができる。本発明のこの実施形態はまた、図示する任意選択的なフィルタ経路に加えて他のフィルタリング経路も含むことができる。図示するものおよび図示するものに対する追加であるものとして記載されるものとを含む各フィルタ経路は、信号処理回路への末梢位相波形入力の所定の位相をフィルタリングし増幅するように配置されている。増幅器Ａ１、Ａ２、Ａｎは、フィルタＦ１、Ｆ２、Ｆｎに対して別個であるものとして示されているが、実際にはフィルタと結合されることが可能であり、それにより、各フィルタＦ１、Ｆ２、Ｆｎは、増幅率も提供するように配置される。
【００３６】
出力２信号処理回路において追加のフィルタ（図示せず）を設けることができる。この追加のフィルタは、出力２に提供される再構成された信号を平滑化するように配置される。
【００３７】
図１の、フィルタリングされかつ増幅された信号を、混合により、または多重化により再構成することができる。図２（ａ）〜（ｅ）は、こうした変調プロセスがいかに作用するかを示す。図２（ａ）は、図１の入力１に提供される入力信号例を示す。（ａ）〜（ｅ）のグラフの各々の垂直軸は、入力信号の大きさを表し、水平軸は時間を表す。
【００３８】
図２（ａ）は、３つの位相ｐ１、ｐ２、ｐ３に分割された入力信号を示す。各位相ｐ１、ｐ２、ｐ３は、末梢脈波形の位相に対応する。各位相ｐ１、ｐ２、ｐ３は、その特定の位相に対して特有の他の位相と比較して異なる周波数特性を有する。
【００３９】
図２（ｂ）、図２（ｃ）および図２（ｄ）は、図１の結合モジュール６、７、８からのそれぞれの信号出力を表す。結合構成要素６の出力を図２（ｂ）によって示し、結合構成要素７の出力を図２（ｃ）によって示し、結合構成要素８の出力を図２（ｄ）によって示す。各フィルタおよび増幅器対のフィルタリング特性および増幅特性が異なるため、出力は、各結合構成要素６、７、８に対して異なる。結合された出力信号は、再構成ブロック９によって再構成され、再構成ブロック９は、この実施形態では、出力を多重化することにより推定された動脈内波形を生成する。再構成ブロック９の出力２を図２（ｅ）に示す。図２（ｅ）の信号は、図２（ａ）〜（ｄ）の信号と同様に３つの位相ａ１、ｂ２、ｎ３を有し、それにより各位相は、異なるフィルタ−増幅器対によって提供される。この多重化プロセスを無期限に繰り返して、図１の信号処理回路の出力２に連続波を提供することができる。図１の信号処理装置の経路の数を増加させることができ、刺激された動脈内波形に再構成される前に別個に処理することができる、末梢脈波形の別個の位相の数は、対応して増加する。
【００４０】
再構成ブロック９の出力を受け取るフィルタを含めることができる。このフィルタは、多重化信号出力の間のいかなる不連続性も平滑化する。
【００４１】
図３は、図１に示すものに類似する信号処理回路の実施形態を示す。図１のものに対応する機能を含むブロックおよび特徴には、同一の番号を付している。この実施形態では、増幅器Ａ１、Ａ２、Ａｎの出力は、再構成ブロック９により元の入力１と結合される。図１の回路とは異なり、図９の回路は、入力１を再構成ブロック９に結合する追加の接続を有している。
【００４２】
図４は、より複雑なフィルタＦを利用する、本発明の実施形態による信号処理回路の別の実施形態を示す。フィルタＦは、入力１から入力信号を受け取り、後に増幅される複数の信号ｆ１、ｆ２、ｆｎを出力する。この実施形態では、単一フィルタのみが必要である。フィルタＦは、ソフトウェアによってプログラムされる電子フィルタであり得る。フィ
ルタリングされた各信号は、別個の増幅器Ａ１、Ａ２、Ａｎによって増幅されているように示されているが、フィルタリングされた信号ｆ１、ｆ２、ｆｎの増幅を、フィルタＦによって提供することができる。そして、フィルタリングされかつ増幅された信号と元のフィルタリングされていない信号とは、再構成ブロック９によって再構成される。
【００４３】
図５および図６は、カフ波形から推定されたｃＳＢＰとトノメトリ波形から推定されたｃＳＢＰとの比較を示し、カフ波形は、（１）ＧＴＦを使用して最初にこの波形を局所末梢動脈内波形に変換し、（２）局所末梢動脈内波形をｐＳＢＰおよびＤＢＰのオシロメトリック測定値を用いて校正し、（３）この校正された局所末梢動脈内波形からＭＡＰａｖｗａｖｅを取得し、（４）ＭＡＰａｖｗａｖｅおよびＤＢＰを使用して、元のカフ波形を変換することから得られる中心血圧波形を校正することによって、校正される。
【００４４】
図６は、カフから推定されたｃＳＢＰとトノメトリ法から推定されたｃＳＢＰとの間の２．２（±５．２）ｍｍＨｇの平均差（±ＳＤ）を示す。図の上方の破線および下方の破線は、平均差±（２Ｄの２倍）、すなわち±１０．４を表す。
【００４５】
この手法は、代替的な手法より精度が低いことが示されている。これは、心周期の時相内で異なるカフ波形の歪みによる可能性がある。動脈内波形の部位（たとえば中心対末梢）による変化は、主に収縮期中に発生し、拡張期では主に指数関数的減衰として残り、動脈内波形からのカフ波形の偏差もまた、拡張期中に発生する。拡張期中のカフ波形のこの予測できない偏差は、特に拡張期の再構成に対し、ＧＴＦを使用する場合のカフ波形の動脈内波形への変換の精度に影響を与える。
【００４６】
図７および図８は、本発明の実施形態のようにカフ波形を局所動脈内波形に変換する好ましい方法による方法の例の結果を示す。本方法は、対象から収集されたカフ圧力測定値を使用して、カフ圧力波形から収縮期波形位相および拡張期波形位相を抽出することによって動脈脈波形を再構成する。
【００４７】
図７は、カフから推定されたｃＳＢＰとトノメトリ法から推定されたｃＳＢＰとの間の１．６２（±４．４７）ｍｍＨｇの平均差（±ＳＤ）を示す。図の上方の破線および下方の破線は、平均差±（ＳＤの２倍）、すなわち±８．９４を表す。
【００４８】
特に、図７および図８は、（１）後述する時分割多重方法を使用してカフ波形を局所末梢動脈内波形に変換し、（２）局所末梢動脈内波形をｐＳＢＰおよびＤＢＰのオシロメトリック測定値で校正し、（３）この校正された局所末梢動脈内波形からＭＡＰａｖｗａｖｅを取得し、（４）ＭＡＰａｖｗａｖｅおよびＤＢＰを使用して、元のカフ波形を変換することから得られる中心血圧波形を校正することによる、カフ波形から推定されたｃＳＢＰとトノメトリ波形から推定されたｃＳＢＰとの比較を示す。
【００４９】
ここで、実行される具体的な処理例と得られる結果とに関連して、本発明の具体的な実施形態について説明する。
【００５０】
方法
収縮期中に発生する末梢動脈内脈波の急な立上りを再構成するために必要な、収縮期中に発生する高周波特徴を抽出するために、カフ圧波形Ａに第１フィルタが適用される。フィルタリングの結果は波形Ｂである。
【００５１】
波形Ｂは、増幅され波形Ａに加算されて、波形Ｃが得られる。波形Ｃは、カフ波形で喪失される動脈内波形の収縮期の高周波成分を組み込むように再構成されたカフ波形である。
【００５２】
末梢動脈内脈波の拡張期部分を再構成するために必要な高周波特徴を抽出するために、カフ圧波形Ａに第２フィルタが適用される。フィルタリングの結果は波形Ｄである。
【００５３】
波形Ｄは、増幅され波形Ａに加算されて、波形Ｅが得られる。波形Ｅは、カフ波形で喪失される動脈内波形の拡張期の指数関数的減衰およびあり得る拡張期ピークを組み込むように再構成されたカフ波形である。
【００５４】
脈波形ＣおよびＥはまとめられ、再構成された動脈脈波が得られる。脈波Ｃの収縮期部分を使用して、末梢動脈内脈波の収縮期部分が再構成され、脈波Ｅの拡張期部分を使用して、末梢動脈内脈波の拡張期部分が再構成される。
【００５５】
そして、再構成された動脈脈波の不連続性のいかなる点も除去するように、第３フィルタが適用される。第３フィルタは、好ましくはゼロ位相フィルタである。
【実施例】
【００５６】
以下の実施例を、上記方法に従って行った。実施例方法および序文で記載したトノメトリ法を、１００の対象から収集されたデータに対して行った。両方法を使用して各対象から３つの測定値が得られた。
【００５７】
実施例は、ＭＡＰにＭＡＰとｐＳＢＰとの差の５０％を足したものに等しい単一ＭＣＰを使用し、心周期の２つの時相、すなわち収縮期および拡張期に対して逐次フィルタリングを適用した。
【００５８】
カットオフ周波数が３．２Ｈｚである７次バターワースフィルタを第１フィルタとして使用して、収縮期中の最高周波数を特定した。収縮期の最高周波数成分を、元の脈波に戻される前に１．７５倍増幅した。
【００５９】
カットオフ周波数が１．７５Ｈｚである４次のバターワースフィルタを第２フィルタとして使用して、拡張期中の最高周波数を特定した。拡張期の最高周波数成分を元の脈波に加算した。
【００６０】
２つの別個のフィルタリング操作からもたらされる２つの脈波の収縮期位相および拡張期位相を、単一脈波にまとめた。
【００６１】
第３フィルタは、カットオフ周波数が６Ｈｚである３次のバターワースローパスフィルタであった。
【００６２】
結果
図７および図８において、白い点は、実施例のフィルタ特性を導出するために使用されたデータのセットを表す。黒い菱形は、本方法を検証するために使用されたデータのセットを表す（このデータのセットは、フィルタ特性を導出するために使用されたデータのセットを含み、したがって、図の菱形のうちのいくつかの上には対応する円がプロットされている）。
【００６３】
カフ圧波形から推定されるｃＳＢＰとトノメトリ法を介して得られるｃＳＢＰとの間に、両タイプの波形がｐＳＢＰおよびＤＢＰのオシロメトリック値を使用して校正される時に、密な関係を見ることができる。
【００６４】
上述した実施形態に対して、追加として、置換としてまたは削除としてでも、さまざま
な変更を行うことができ、それらは添付の特許請求の範囲の範囲内にあるように意図されている。たとえば、１つの代替実施形態では、末梢脈波形を、いかなるフィルタ処理または他の信号処理も適用される前に、収縮期位相および拡張期位相等、その異なる位相に時分割することができる。そして、それぞれの分割された位相信号を、必要に応じて、たとえば適切にフィルタリングし、その後、多重式に再結合することによって処理することができる。
【図１】
2014534023.tif 000003
【図２】
2014534023.tif 000004
【図３】
2014534023.tif 000005
【図４】
2014534023.tif 000006
【図５】
2014534023.tif 000007
【図６】
2014534023.tif 000008
【図７】
2014534023.tif 000009
【図８】
2014534023.tif 000010
【国際調査報告】
2014534023.tif 000011
2014534023.tif 000012
2014534023.tif 000013