JP4389239B2 - 聴診信号における律動を推定するための方法および装置 - Google Patents

Description

【０００１】
本発明は、聴診信号を処理するための装置に関する。この装置は、聴診信号を受信すると共にこの聴診信号の包絡線を表す包絡線信号を供給するための手段と、聴診信号の推定律動を表す信号を算出するための推定装置とを備える。
【０００２】
更に、本発明は、聴診信号を処理する方法、コンピュータ読み取り可能媒体に関し、具体的には、心臓学において用いるのに適した電子聴診器に関する。
【０００３】
近年、医師は、心血管の疾病の診断のためにすばらしい多くの機器を自由に用いられるようになった。かかる機器の１つは、周知の聴診器であり、心臓および隣接する大きな血管から発する音を検出するために用いられる。聴診器の導入によって、聴診は、心臓の機能的完全性の貴重な情報を臨床家に与える技術へと発展した。今日では、臨床的な研究によって得られる組み込まれた知識を用いることによって、聴診信号内に存在する情報を電子的に処理することができる。従って、臨床医を、つまらない作業から解放することができる。かかる作業の１つは、聴診信号に基づいて鼓動する心臓の心拍数を推定することであろう。
【０００４】
聴診信号は、血液をポンプ作用で送出する過程で心臓の収縮および拡張から起こる特徴的な音、すなわち、血液の加速および減速からの振動によって生じる音を含む。鼓動する心臓を聴診する場合、それは、聴かれている鼓動振動数ではなく、特徴的な音である。更に、正常な状態の健康な中年の人々では、特徴的な音は、多くの場合、それぞれ心臓の収縮および拡張に関連するＳ１およびＳ２と示され、時間的に十分に大きい間隔によって分離されて、これらの特徴的な音を個別の音として認知することを可能とする。すなわち、特徴的な音は、典型的に、成人では１５０ないし２００ミリ秒以上によって、子供では約５０ミリ秒によって分離されている。通常、特徴的な音は、約１５Ｈｚおよび１００Ｈｚ間に分布し３０ないし４０Ｈｚにピークがある連続周波数スペクトルを有する。
【０００５】
しかしながら、特徴的な音には、心臓内の異なる種類の振動から発生する音が重なる場合がある。これらの音は、いわゆる雑音であり、血液の乱れまたは構造の粗動のような機構によって起こる恐れがある。雑音は、２つの主要なグループに分けられる。すなわち、第１のグループは、ランダムな「雑音性の」音から成り、第２のグループは、連続的な狭帯域周波数スペクトルを有する「音楽的な」音から成る。人において、「音楽的な」音は、通常、心臓の鼓動振動数に関連することが、経験によって示されている。雑音の周波数スペクトルは、５００ないし２０００Ｈｚまでの分布を有する場合がある。
【０００６】
また、特徴的な音は、ある時点またはランダムな時点で現れる心臓の偽の収縮によって妨害されることがある。更に、特徴的な音は、ランダムな時点で現れることもあり得る。この現象は、不整脈として既知である。
【０００７】
雑音があればそれを含む、鼓動する心臓の音を、例えば患者の胸部上に配置した変換器によって聴診する場合、この変換器によって低周波数のランブル雑音が拾われる。このランブル雑音は、心拍の周波数のすぐ下からＤＣ方向に主に分布する周波数内容を有する場合がある。
【０００８】
深刻な場合には、これらの妨害作用が、聴診信号内の特徴的な音を完全に破壊するか、または、少なくとも、従来技術の技法を用いて近似的な心拍数でさえ推定することが不可能となり得る程度まで、聴診信号を妨害する可能性がある。
【０００９】
しかしながら、非電子聴診器を用いて訓練した人では、特徴的な心音という形態の情報は、繰り返しごとに比較的容易に認知することができ、従って、計数は比較的容易である。これは、人による聴診信号の知覚が、妨害作用またはアーチファクトを減じることにおいて極めて優れており、心拍数を推定するための情報を提供する聴診信号内の特徴に注目するという事実による。
【００１０】
米国特許第４，９７２，８４１号は、推定した脈拍数を表示するための脈拍数ディスプレイを有する電子聴診器を開示する。この聴診器では、変換器が、いわゆるコロトコフ音を電気信号に変換する。コロトコフ音は、血液が血管を流れることが可能となった際にへこんだ血管の壁が振動することによって生成される音である。電気信号は増幅されてカウンタに供給される。カウンタでは、単位時間当たりの検出脈拍数を算出し、この結果をデジタルの脈拍数として表示する。
【００１１】
しかしながら、この特許は、測定音が予測される通常の心音と大きく異なる場合がある、人から発せられる心拍数の推定には関連しない。この特許は、心拍数に関連する聴診信号内の情報と取り違えるかまたは混同される恐れのあるアーチファクトを抑制するように構成された手段を何ら開示していない。
【００１２】
米国特許第４，４３６，０９６号は、デジタル式の心拍数聴診器を開示する。変換器が、人の患者の心音を表す電気信号を生成し、この電気信号を受信するためにフィルタが接続されている。フィルタは、狭帯域バンド・パス・フィルタであり、その通過帯域は、特徴的な心音の周波数、例えば３３Ｈｚを中心とする。フィルタは、これによって、信号対雑音比を向上させるように構成されており、患者の衣服の上から変換器を使用可能とする。濾波された信号はパルスに変換され、このパルスはカウンタによって計数され、デジタルの心拍数指示として表示される。
【００１３】
この場合も、この特許は、測定音が予測される通常の心音とは大きく異なる場合がある、人から発せられる心拍数の推定には関連しない。また、この特許は、心拍数に関連する聴診信号内の情報と取り違えるまたは混同される恐れのあるアーチファクトを抑制するように構成された手段を何ら開示しない。
【００１４】
国際出願ＷＯ９７／０００４５号は、鼓動する心臓によって生成される音波の特徴を調べるための聴診器を開示する。この聴診器は、閾値波を確立するための手段を備えており、この閾値波を超えると、周期的な基準複合体が現れて識別される。更に、この聴診器は、隣接する基準複合体上の同様の時点の間で経過する時間を心拍レートの指示として測定する時間測定装置も備えている。
【００１５】
この聴診器では、心臓の収縮および拡張から発する心音の間に実質的にノイズが存在しないという仮定に基づいて、心拍数を推定する。従って、この聴診器は、心臓の雑音または不整脈のある人の心拍数の推定に用いるには適していない。
【００１６】
心拍数を推定するための他の装置は、聴診信号の包絡線を決定するための手段を用いる。この包絡線は、聴診信号の絶対値を取り、この絶対値信号を帯域濾波することによって、決定される。
【００１７】
結果として、従来技術には、ある程度の雑音および／または不整脈が聴診信号に存在する場合に心拍数が誤って推定されるという問題が伴う。
【００１８】
このため、本発明の主な目的は、心臓の実際の血液のポンプ式送出による収縮および拡張から発する心拍数を推定する、より精密かつロバストな心拍数推定を提供することである。
【００１９】
これが達成されるのは、最初の段落で述べた装置において、推定装置が、包絡線信号の一部を選択するように構成された手段を備えており、前記一部が、第１信号として包絡線信号における一連のサイクルの少なくとも１つを含み、この手段が、更に、第１信号と包絡線信号との間の一致度を算出するように構成されていることを特徴とする場合である。
【００２０】
結果として、心拍数に関連する情報と取り違えるまたは混同される恐れのある聴診信号内のアーチファクトまたはランダム・ノイズを抑制することができる。これによって、心拍を非常に明確に規定した信号が得られる。
【００２１】
聴診信号を処理するための好都合な装置では、推定装置は、相互相関関数を算出する。この相互相関関数は、一致度が高い大きな信号値を有すると共に一致度が低い小さな信号値を有する明確に規定された信号を供給する。
【００２２】
バイアスをかけた聴診信号の一部が、バイアスをかけた聴診信号の一連のサイクルの１つを表す場合、この一部は、いくつかの面で理想的と見なされる基準サイクルを表すことができる。
【００２３】
聴診信号を処理するための別の好都合な装置では、推定装置は、自己相関関数を算出する。自己相関関数は、心拍数が所与の時間間隔に渡って比較的安定している場合に、特に有利である。
【００２４】
聴診信号の品質は、聴診信号の一致度を表す信号において、以下の３つの項目の少なくとも１つを検証することによって、確認することができる。
【００２５】
ａ）位置決めされた極値間の時間差は、所定の制限内でなければならない。
【００２６】
ｂ）この時間差の平均に対する最小および最大の時間差の比は、所定の制限内でなければならない。
【００２７】
ｃ）極値の位置決めにおける相関の結果の大きさは、所定の制限内でなければならない。
【００２８】
このように、律動推定の統計データに直接関連する好適な制限によって、好ましい品質レベルを規定することができる。
【００２９】
バイアス・プロセッサは、聴診信号のＡ特性周波数重み付けバージョンまたは聴診信号の近似的なＡ特性周波数重み付けバージョンを算出するためのフィルタを含むことができる。シミュレーションによって、聴診信号のＡ特性周波数重み付けバージョンは、様々な状況において安定して心拍数を推定するために十分な情報を含むことが示されている。
【００３０】
所定の周波数範囲における聴診信号の二重微分に対応する近似によって、聴診信号のＡ特性周波数重み付けバージョンを算出する場合、濾波プロセスに必要な計算上の労力は最小となる。
【００３１】
バイアス・プロセッサは、包絡線信号を濾波するための適応バンド・パス・フィルタを備えていても良く、前記バンド・パス・フィルタは、少なくとも選択可能な高い通過帯域および低い通過帯域を有し、更に、第１の周波数帯域における包絡線信号内容の第２の周波数帯域における包絡線信号内容に対する比を算出すると共に算出した比に応じて通過帯域を選択するコントローラを備える。このため、比較的広い範囲の周波数に渡って聴診信号内の律動を推定することができる一方で、信号対雑音比の大きな劣化が回避される、すなわち、バンド・パス・フィルタの比較的狭い通過帯域を維持することができる。
【００３２】
聴診信号は、サンプル・レートで到着するサンプルを含むことができ、この装置は、サンプル・レートに対応するレートで動作する同期プロセッサを備えることができ、更に、要求によって開始する時間間隔で動作する非同期プロセッサも備える。これによって、リアル・タイムで聴診信号の新たなサンプルを取得することができる一方で、聴診信号の前のサンプル値に作用する相関関数を算出する。
【００３３】
あるいは、本発明の目的が達成されるのは、この装置において、手段が、Ａ特性周波数重み付けまたは近似的なＡ特性周波数重み付けに対応する周波数応答を有するフィルタを備え、このフィルタは、少なくとも対象の周波数範囲について算出された二重微分として具現化されることを特徴とする場合である。
【００３４】
結果として、鼓動する心臓の特徴的な音および心拍数に関連する雑音に基づいて、鼓動する心臓の心拍数を推定することができる。このため、特徴的な音が心拍数を推定するための十分な情報を与えない場合に、心拍数に関連する雑音に基づいて心拍数を推定することができる。更に、低周波数ランブル雑音は軽減され、高周波数成分は強調され、これによって、時間に対する特徴的な音の相対的な位置決めの明確に規定された情報が得られる。この明確に規定された情報は、心拍数を電子的に推定するために適している。
【００３５】
二重微分に対応した周波数応答を有するフィルタによって聴診信号を濾波する場合、聴診のＡ特性周波数重み付けバージョンを算出するために必要な計算上の労力は軽減する。
【００３６】
対象の周波数範囲は、２０００Ｈｚ未満の周波数範囲に及ぶ場合がある。特徴的な心音および心拍数に関連する雑音は、この周波数範囲に位置する周波数成分を含み、鼓動する心臓の律動を推定するために、この範囲において十分な情報を与える。
【００３７】
更に、本発明は、聴診信号における律動を推定する方法にも関する。
【００３８】
本発明を、好適な実施形態に関連付けて、図面を参照しながら、以下で更に十分に説明する。
【００３９】
図１は、近似的なＡ特性周波数重み付け聴診信号を算出するバイアス・フィルタによって聴診信号を処理するためのプロセッサのブロック図を示す。聴診信号は、音響信号を表す信号であり、鼓動する心臓の音を表す聴診信号を供給するように位置付けられたある種の変換器によってピックアップすることができる。以下では、聴診信号をＶＣＧと示し、これは、デジタルまたはアナログの信号とすることができる。
【００４０】
このブロック図は、バイアス・プロセッサ１０６を含む。ＶＣＧ信号は、バイアス・フィルタ１０１において濾波される。このフィルタの周波数応答は、ＶＣＧ信号の近似的なＡ特性周波数重み付けに対応する。Ａ特性周波数重み付けは、可聴周波数範囲すなわち約２０ないし２００００Ｈｚにおける信号振幅の人による知覚をモデル化する、人の耳の線形モデルである。Ａ特性周波数重み付けは、通常、２０００Ｈｚの中心周波数の二次バンド・パス・フィルタとしてモデル化される。Ａ特性周波数重み付けモデルは、時間に対するＶＣＧ信号の二重微分として近似され、これによって、ＶＣＧ信号のＡ特性周波数重み付けバージョンすなわちＡ−ＶＣＧが与えられる。この近似は、約２０ないし２００Ｈｚ間の周波数について有効である。このため、バイアス・フィルタ１０１は、二次高域フィルタの傾きとして見ることができる。このフィルタは、これによって、聴診信号内に存在し得る低周波数のランブル信号を減じる。更に、フィルタは、特徴的な心音の周波数スペクトルの上部、および最終的に雑音を強調する。しかしながら、Ａ特性周波数重み付け信号Ａ−ＶＣＧも、特徴的な心音にも最終的に雑音にも関連しない源から発するノイズを強調し、これによって、信号対雑音比を大きく劣化させる場合がある。従って、Ａ特性周波数重み付け信号Ａ−ＶＣＧは、フィルタ１０２によって低域濾波される。低域濾波された信号は、次いで、入力信号として整流器１０３に供給される。整流器１０３は、入力信号の絶対値を算出し、これによって出力信号を供給する。整流器１０３からの出力信号は、フィルタ１０４によって帯域濾波される。この帯域濾波された信号は、Ａ特性周波数重み付け信号Ａ−ＶＣＧの包絡線を構成する。このバンド・パス・フィルタは、予想心拍数の範囲、例えば０．５から３．３ＨＺに対応する通過帯域を有する。更に、バンド・パス・フィルタの拒絶帯によって、いかなるＤＣ成分も、基本的な心拍周波数を超えるいかなる周波数成分すなわち比較的高い周波数成分も、減じられる。
【００４１】
反復性の明確に規定されたサイクルとして存在する帯域濾波された信号内の情報は、これによって、心拍数を表す。心拍数は、推定装置１０５によって推定することができる。推定装置は、帯域濾波された信号の明確に規定された極値のタイミングを決定するピーク・ロケータの形態とすれば良い。あるいは、推定装置は、推定心拍数を表すある周波数を有する方形信号を発生する閾値またはゼロ交差検出器の形態とすれば良い。推定装置が閾値検出器から成る場合、バイアスされた聴診信号および閾値は、双方とも比較器に入力される。閾値は、一定の値または適応的な閾値の形態とすることができる。適応的な閾値は、バイアスされた聴診信号を濾波するロー・パス・フィルタによって供給することができる。
【００４２】
図２は、特徴的な心音、雑音、およびランブル雑音の予想周波数スペクトルを示す。このスペクトルは、粗い近似であり、例示の目的のためのみのものである。スペクトルの大きさアンペアは、任意の単位であり、周波数Ｆ（Ｈｚ）の関数として示す。鼓動する心臓の収縮および拡張に関連する特徴的な心音の周波数スペクトル２０１は、概略的に示されている。文献では、収縮に関連する音をＳ１と示し、拡張に関連する音をＳ２と示す。Ｓ１およびＳ２の周波数内容の周波数スペクトルは、主として２０および５０Ｈｚ間に位置付けられているが、上は１００ないし２００Ｈｚまで、下は１０Ｈｚ未満まで分布する場合がある。
【００４３】
しかしながら、雑音の周波数内容は、通常、より高い周波数に分布する。これは、周波数応答２０２によって示されている。「音楽的な」音から成る雑音の周波数内容は、比較的狭い帯域の周波数範囲２０６に分布している。この比較的狭い帯域は、心拍数の推定に使用可能な情報を含む。通常、雑音は、心拍数に関連する情報を含まない比較的広い帯域２０５にも関連する。それらの周波数範囲の絶対数値の定義およびこれらの帯域における聴診信号の大きさは、対象の雑音によって異なる。
【００４４】
低周波数ランブル雑音の周波数スペクトル２０４は、通常、特徴的な心音の周波数スペクトル２０１に重なる。
【００４５】
伝達関数２０３は、バイアスをかけた聴診信号を算出するためのフィルタに好適な伝達関数を示す。フィルタは、これによって、特徴的な心音の周波数スペクトルの上部の周波数成分の振幅および、雑音の周波数スペクトルの周波数成分の振幅を強調することができる。更に、聴診信号内に存在し得る低周波数ランブル雑音は軽減される。伝達関数２０３は、デケード当たり高域４０ｄＢの傾きを含む。この傾きは、近似的なＡ特性周波数重み付け関数の一部とすることができ、この場合、近似は、２０００Ｈｚの中心周波数を有する二次バンド・パス・フィルタとして実施される。
【００４６】
図３は、相関器によって聴診信号の鼓動周波数を算出するためのプロセッサのブロック図を示す。この実施形態では、受信される聴診信号ＶＣＧは、デジタル信号であり、サンプル周波数ｆｓで到着するデジタル・サンプルを含む。例えば、ｆｓは４０００Ｈｚに等しい。
【００４７】
この実施形態は、サンプル周波数と同期したレートでＶＣＧサンプルを処理するバイアス・プロセッサ３０１と、バイアス・プロセッサに対して非同期で動作することができ、要求された場合に心拍数を表す信号ＨＲを供給する推定装置３０２とを含む。推定装置３０２は、バイアス・プロセッサ３０１から、バッファ３０３を介して、サンプルの形態のデータを受信する。バッファ３０３は、循環バッファまたは先入れ先出し（ＦＩＦＯ）メモリの形態とすれば良い。
【００４８】
バイアス・プロセッサ３０１では、整流器３０４によって、デジタル聴診信号ＶＣＧを受信し、ＶＣＧ信号の各絶対値を計算する。次いで、聴診信号の包絡線を算出するために、フィルタ３０５によって、絶対値を帯域濾波する。このバンド・パス・フィルタは、予想心拍数の範囲、例えば０．５ないし３．３Ｈｚに対応する通過帯域を有する。更に、バンド・パス・フィルタの拒絶帯によって、いかなるＤＣ成分も、また、基本的な心拍周波数を超える周波数成分すなわち比較的高い周波数成分も軽減される。
【００４９】
バンド・パス・フィルタからのサンプルは、同じ時点でバッファ３０３に供給される。バッファは、例えば５秒の信号シーケンスに対応した多数のサンプル点のための容量を有する循環バッファとすることができる。バンド・パス・フィルタから新たなサンプルが到着すると、バッファ内の最も古いサンプルは重ね書きされる。推定装置３０２が開始すると、バッファからサンプルを取得する。推定装置は、例えばデジタル表示の更新用の情報を与えるために、ある一定の時点または不規則な時点で開始することができる。
【００５０】
バンド・パス・フィルタによってバッファ３０３を介して供給された包絡線信号は、円滑であり、反復性サイクルを含むが、特定の時間間隔の間に１つ以上のサイクルが存在するか否かを判定することが難しい場合がある。しかしながら、特定の時間間隔内の正確なサイクル数を判定するために、包絡線信号の情報内容を強化することができる。これは、包絡線信号自体に対する包絡線信号の一致度を算出することによって行われる。これは、包絡線信号の大部分のサイクルが、血液をポンプ式に送出する過程における実際の相互関連する心臓の収縮および拡張を表すという仮定に基づいている。包絡線信号の一致度は、相関器３０６が包絡線信号の自己相関を計算することによって算出される。
【００５１】
好適な実施形態では、相関器は、最大の振幅を有する包絡線信号内のサイクルを選択するための手段を備える。これは、聴診信号内の最大の振幅が、アーチファクトを含む可能性が最も低い通常の心臓の収縮／拡張によって生じるという有望な仮定のもとで行われる。選択したサイクルの範囲は、ゼロ未満の最も近い極小を探すことによって見出される。これによって、ゼロを超える大きさの極小になることが回避される。この最小値の検索は、最大の振幅の位置に対して前方向および後方向の双方に適用される。選択したサイクルを格納し、包絡線信号の残りと相互相関させる。あるいは、信号自体に対する信号の一致度を算出することができる畳み込みまたはいずれかの他の関数を使用可能である。
【００５２】
図４は、相関器およびＡ特性周波数重み付けを含むバイアス・フィルタによって聴診信号の鼓動周波数を算出するためのプロセッサのブロック図を示す。この好適な実施形態では、全体の性能を向上させるために、Ａ特性周波数重み付けを含むバイアス・フィルタおよび相関を結合している。すなわち、心臓の特徴的な音が大きく歪められていても、鼓動する心臓の律動を推定することが可能である。
【００５３】
聴診信号ＶＣＧは、サンプル周波数ｆｓで到着するサンプル値を含むデジタル信号の形態である。例えば、ｆｓ＝４０００Ｈｚである。聴診信号は分岐されている。第１の分岐は、聴診信号を、整流器または絶対値検出器、ＡＢＳ４０６に供給する。聴診信号の絶対値は、更に、１００Ｈｚを超える周波数成分を減じるように約３Ｈｚの遮断周波数を有するロー・パス・フィルタ４０７に供給される。第２の分岐は、聴診信号を、聴診信号の二重微分を計算するフィルタ４０１に供給する。フィルタ４０１の周波数応答は、２０００Ｈｚ未満の周波数範囲における二次高域傾きまたはＡ特性周波数重み付け近似に相当する。聴診信号の二重微分は、ロー・パス・フィルタ４０２に供給され、ここで、フィルタ４０１において微分によって強調された比較的高い周波数ノイズ成分を減衰させる。フィルタ４０２は、約１００Ｈｚの遮断周波数を有するものとすれば良い。フィルタ４０２によって低域濾波された信号は、整流器または絶対値検出器、ＡＢＳ４０３に供給される。聴診信号の絶対値は、更に、ロー・パス・フィルタ４０７に供給され、１００Ｈｚを超える周波数成分が減じられる。
【００５４】
デシメータＤＥＣ４０５および４０８は、１００Ｈｚを超える周波数成分が減じられているので、より低いサンプル・レートｆｓＲを発生する。これは、例えば２００Ｈｚとすれば良い。このため、聴診信号の後の信号処理において必要な計算上の労力を軽減することができる。
【００５５】
バンド・パス・フィルタＢＰＦ＃１およびＢＰＦ＃２から成るバンド・パス・フィルタ４１１は、デシメートされた聴診信号を受信し、帯域濾波された信号を算出する。バンド・パス・フィルタは、例えば０．５ないし３．３Ｈｚの予想心拍数範囲に対応する通過帯域を有するものとすれば良い。
【００５６】
ＤＣ成分および２５Ｈｚを超える周波数成分が減じられることを注記しておく。バンド・パス・フィルタの中心周波数は、律動推定装置の性能を向上させるために、適応的に調節することができる。これについては、後に論じる。
【００５７】
このように帯域濾波された信号は、低いサンプル・レートでの聴診信号および二重微分された聴診信号の包絡線を表す。
【００５８】
デシメータＤＥＣ４１２および４１３は、２５Ｈｚを超える周波数成分が減じられているので、更に低いサンプル・レートを発生する。これは、例えば５０Ｈｚとすれば良い。このため、聴診信号の後の信号処理において必要な計算上の労力を更に軽減することができる。
【００５９】
バンド・パス・フィルタからのサンプルを、同じ時点で、バッファ４１４および４１５に供給する。バッファは、例えば５秒の信号シーケンスに対応した多数のサンプル点のための容量を有する循環バッファとすることができる。推定装置４１０が開始すると、バッファからサンプルを取得する。推定装置は、例えばデジタル表示の更新用の情報を供給するために、ある一定の時点または不規則の時点で開始することができる。
【００６０】
バッファ４１４および４１５から取得されたサンプルは、オートスケール手段４１６および４１７によってそれぞれスケーリングされる。オートスケール手段は、ある数のスケーリングされたサンプルにおけるピーク値が１．０の値を取ることを保証することによって、受信サンプル値をスケーリングすることができる。例えば、サンプル数は、バッファ４１４および４１５に格納された信号シーケンスに対応し、これは１．０に等しいピーク値に正規化された信号シーケンスである。スケーリングされたサンプルは、加算器４１８において、等しい重みを加算される。なぜなら、２つの信号内の情報は等しく重要であるからである。加算の結果は、和信号として供給される。あるいは、２つの信号に、異なる重みを加算しても良い。
【００６１】
相関器Ｘ−Ｃｏｒｒ４１９によって、和信号の一致度を算出する。好適な実施形態では、相関器は、最大の振幅を有する包絡線信号内のサイクルを選択するための手段を備えている。選択されたサイクルの時間範囲は、ゼロ未満の最も近い極小を検索することによって見出される。これによって、ゼロを超える大きさを有する極小になることが回避される。この最小値の検索は、最大の振幅の位置に対して前方向および後ろ方向の双方に適用される。選択したサイクルを格納し、包絡線信号の残りと相互相関を取る。
【００６２】
ピーク・ロケータ４２０は、相互相関の結果を表す信号における極値を識別する。これは、例えば、相互相関の結果の一次導関数を算出し変曲を識別することによって実行可能である。好ましくは、ピーク・ロケータは、相互相関の結果を表す信号における最大値を識別する。
【００６３】
好ましくは、品質推定装置４２１によって、極値の位置決めを処理する。推定装置は、相互相関の結果が、聴診信号の律動がある場合にこれを推定することに適しているか否かを決定する。品質推定装置４２１は、品質を推定するプロセスにおいて３つのステップを含むことができる。
【００６４】
１．２つの極値間の時間差の平均は、所定の制限内でなければならない。
【００６５】
２．この時間差の平均に対する最小および最大の時間差の比は、所定の制限内でなければならない。
【００６６】
３．極値の位置決めにおける相互相関の結果の大きさは、所定の制限内でなければならない。
【００６７】
上述の項目１ないし３の１つ以上に当てはまらない場合、聴診信号内の律動を推定することが不可能であると判定される。あるいは、上述の項目に当てはまる場合は、聴診信号内の律動を推定することができる。律動の推定は、単位時間当たりの相互相関の結果における極値の数を計数し、これを１分当たりの鼓動数として表示することによって推定することができる。
【００６８】
代替的な実施形態では、相関器３０６または４１９は、自己相関関数を算出するように構成することができる（すなわち、バッファ３０３または４１４、４１５の容量に対応する長さを有する２つの同一の信号を互いに相関させる）。自己相関関数の算出から得られる信号を用いて、聴診信号の律動を推定し、聴診信号の品質を調べることができる。聴診信号内の極大を検索することによって、聴診信号の一次周波数を見出す。この検索は、０．３秒（２００ｂｐｍ）から開始し、一次周波数は、固定閾値を超える第１の極大として規定される（自己相関における２秒の制限前に極大が見出されない場合、一次周波数は決定されない）。一次周波数を承認する前に、以下の項目を確認しなければならない。
【００６９】
１．バッファ信号の最後の２秒の平均レベルは固定閾値を超えなければならない。この確認は、間違った種類の信号（または存在しない信号）に基づいた推定を防ぐために行われる。
【００７０】
２．第１のピークの２倍の時間位置決めされた自己相関におけるピーク（±７．５％）が存在しなければならない。このチェックは、バッファ信号が周期的か否かを確認するために行われる。また、このピークは、固定の閾値を超えなければならない。
【００７１】
３．Ｎ個のピークが、第１のピーク（ｐｅａｋ１）の左部に見出された場合、等しく離間した位置において見出される値の平均＝（ｐｅａｋ１）／Ｎ＋１→（Ｎ・ｐｅａｋ１）／（Ｎ＋１）は、固定の閾値未満でなければならない。このチェックは、包絡線が、（Ｎ＋１）・（一次周波数）に等しい鼓動周波数を有する成分を含む場合に、低い心拍数推定値を防ぐために行われる。
【００７２】
４．この推定は、ある時間間隔に渡って安定していなければならない。一次周波数が、毎分４０の鼓動（ｂｍｐ）未満であると判定された場合、この推定は、８倍安定していなければならない（±３０％）＝１．６秒。一次周波数が、４０ｂｍｐ以上であると判定された場合、この推定は、４倍安定していなければならない（±３０％）＝０．８秒。このチェックは、推定の妥当性を増すために行われる。
【００７３】
他の妥当性基準が実行可能であるが、上述の基準は十分に機能すると証明されている。
【００７４】
図５は、適応バンド・パス・フィルタを示す。このフィルタは、図４に示したバンド・パス・フィルタ４１１の代わりに用いることができる。適応フィルタは、比較的広い周波数範囲に渡って聴診信号における律動の推定に用いることができる一方で、信号対雑音比を大幅に劣化させることを回避する、すなわち、バンド・パス・フィルタの比較的狭い通過帯域を維持することができる。
【００７５】
バンド・パス・フィルタ５１２は、２つの入力信号Ｉ１およびＩ２を受信し、２つの帯域濾波された出力信号Ｏ１およびＯ２を供給する。バンド・パス・フィルタ５１２は、このため、２つのバンド・パス・フィルタＢＰＦ＃１およびＢＰＦ＃２を備える。フィルタＢＰＦ＃１およびＢＰＦ＃２の各々は、比較的低い中心周波数を有する通過帯域ＬＢおよび比較的高い中心周波数を有する通過帯域ＨＢを備える。入力ＬＢまたはＨＢのいずれかにイネーブル信号を供給することによって、通過帯域の１つを選択することができる。この通過帯域の選択は、適応コントローラ５１３によって実行される。
【００７６】
高域フィルタ５０１は、バンド・パス・フィルタ５１２に入力される信号の１つに対応する入力信号を受信する。高域フィルタ５０１は、約０．３Ｈｚの遮断周波数を有し、これによって、適応通過帯域選択に用いられない極めて低い周波数成分を減じる。フィルタ５０１は、高域濾波した信号をロー・パス・フィルタ５０２に供給する。ロー・パス・フィルタは、約１Ｈｚの遮断周波数を有する。
【００７７】
バッファ５０３および５０４は、それぞれ、約３秒の持続時間を有する信号シーケンスに対応するサンプル値を格納することができる。平均値検出器５０５および５０６によって、バッファに格納されたシーケンスの予想値または平均値が算出される。平均値検出器５０５によって算出された平均値は、バンド・パス・フィルタ５１２に入力される全信号の平均値を表す。平均値検出器５０６によって算出される平均値は、バンド・パス・フィルタ５１２に入力される信号の低周波数部分すなわち約１Ｈｚ未満の平均値を表す。
【００７８】
このため、除算装置５０７は、比較的低周波数である信号の分数Ａ／Ｂを算出する。初期状態では、閾値Ｔｈ．は、値０．４５に等しい。この値を、加算器５０８によって、この分数Ａ／Ｂから減算する。判断状態５０９では、周波数比Ａ／Ｂが閾値より大きいか否かを確認する。これが当てはまる場合（Ｔ）、閾値Ｔｈ．を０．３５に等しくセットし、比較的低い中心周波数を有する通過帯域ＬＢを選択する。これは、聴診信号がゆっくり鼓動している心臓から発するという仮定に基づいて行われる。このため、比較的低い中心周波数を有する通過帯域ＬＢを選択することによって、より高い、余分な周波数成分を減じることができる。好ましくは、比較的低い中心周波数は約０．７Ｈｚである。
【００７９】
あるいは、周波数比Ａ／Ｂが閾値Ｔｈ．未満である場合（Ｆ）、閾値Ｔｈ．を値０．５５に等しくセットし、比較的高い中心周波数を有する通過帯域ＨＢを選択する。これは、聴診信号において低周波数内容が乏しいため、聴診は極めて速く鼓動している心臓、例えば子供から発する可能性が高いという仮定に基づいて行われる。比較的高い中心周波数を有する通過帯域ＨＢを選択することによって、聴診信号の律動が正確に推定される可能性が最大となる。好ましくは、比較的高い中心周波数は約１．２Ｈｚである。
【００８０】
図６は、相関器およびＡ特性周波数重み付けを備えるバイアス・フィルタを用いた実施形態に関連する信号を示す。聴診信号ＶＣＧおよびＶＣＧ信号の対応する帯域濾波された包絡線を、時間の関数として示す。更に、Ａ特性周波数重み付けまたは二重微分した聴診信号ｄ^２ＶＣＧ／ｄｔ^２および対応する帯域濾波した包絡線を示す。グラフの解釈を考慮すると、帯域濾波した信号の１つのみを用いて心拍数を推定可能であることは明らかである。しかしながら、上述のように、信号を加算し、相関器に対する入力として用いると、より強力な推定が得られる。これは、信号「ＸＣ−Ｉ」および「ＸＣ−Ｏ」によって示されている。
【００８１】
図７は、聴診信号の包絡線における情報が不十分である場合の例を示す。聴診信号ＶＣＧおよびＶＣＧ信号の対応する帯域濾波した包絡線を、時間の関数として示す。更に、Ａ特性周波数重み付けまたは二重微分した聴診信号ｄ^２ＶＣＧ／ｄｔ^２および対応する帯域濾波した包絡線を示す。グラフの解釈を考慮すると、二重微分した聴診信号ｄ^２ＶＣＧ／ｄｔ^２によって最も強力な推定が得られることは明らかである。しかしながら、聴診信号の大部分について、信号「ＸＣ−Ｏ」が、最も強力な全体的推定を与える。
【００８２】
図８は、Ａ特性周波数重み付け聴診信号の包絡線における情報が不十分である場合の例を示す。聴診信号ＶＣＧおよびＶＣＧ信号の対応する帯域濾波した包絡線を、時間の関数として示す。更に、Ａ特性周波数重み付けまたは二重微分した聴診信号ｄ^２ＶＣＧ／ｄｔ^２および対応する帯域濾波した包絡線を示す。グラフの解釈を考慮すると、帯域濾波したＶＣＧ包絡線信号によって最も強力な推定が得られることは明らかである。しかしながら、この場合も、信号「ＸＣ−Ｏ」が、最も強力な全体的推定を与える。
【００８３】
本発明は、聴診器またはいずれかの他の機器もしくは装置の一部として具現化することができる。あるいは、本発明は、コンピュータのメモリ内にロードすることができ、それから実行することができるコンピュータ・プログラムまたはコンピュータ・プログラムの一部として具現化することができる。コンピュータ・プログラムは、いずれかのデータ記憶装置またはデータ伝送媒体によって配布することができる。この記憶媒体は、磁気テープ、光ディスク、コンパクト・ディスク（ＣＤまたはＣＤ−ＲＯＭ）、ミニ・ディスク、ハード・ディスク、フロッピ・ディスク、強誘電体メモリ、電気的消去可能プログラマブル読み出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュ・メモリ、ＥＰＲＯＭ、読み出し専用メモリ（ＲＯＭ）、スタティック・ランダム・アクセス・メモリ（ＳＲＡＭ）、ダイナミック・ランダム・アクセス・メモリ（ＤＲＡＭ）、強磁性体メモリ、光記憶装置、電荷結合デバイス、スマート・カード等とすることができる。伝送媒体は、例えばローカル・エリア・ネットワーク（ＬＡＮ）、ワイド・エリア・ネットワーク（ＷＡＮ）のようなネットワーク、またはこれらの何らかの組み合わせ、例えばインターネットとすることができる。ネットワークは、有線または無線通信リンクから成るものとすることができる。ネットワークを介して、本発明のソフトウエア実施形態（すなわちプログラム）、またはその部分を、ネットワークを介してプログラムを転送することによって、配布することができる。
【図面の簡単な説明】
【図１】 近似的なＡ特性周波数重み付け聴診信号を算出するバイアス・フィルタによって聴診信号を処理するためのプロセッサのブロック図を示す。
【図２】 特徴的な心音、雑音、およびランブル雑音の予想周波数スペクトルを示す。
【図３】 相関器によって聴診信号の鼓動周波数を算出するためのプロセッサのブロック図を示す。
【図４】 相関器およびＡ特性周波数重み付けを含むバイアス・フィルタによって聴診信号の鼓動周波数を算出するためのプロセッサのブロック図を示す。
【図５】 適応バンド・パス・フィルタを示す。
【図６】 相関器およびＡ特性周波数重み付けを含むバイアス・フィルタを用いた実施形態に関連する信号を示す。
【図７】 聴診信号の包絡線内の情報が不十分である例を示す。
【図８】 Ａ特性周波数重み付け聴診信号の包絡線内の情報が不十分である例を示す。

Claims (29)

1. 聴診信号を処理するための装置であって：
   聴診信号を受信し、前記聴診信号の包絡線を表す包絡線信号を供給するための処理手段（１０６；３０１；４０９）と；
   前記聴診信号の推定律動を表す信号を算出するための推定装置（１０５；３０２；４１０）と；
   を備える装置において、
   前記推定装置は、第１信号として前記包絡線信号における一連のサイクルの少なくとも１つを含む前記包絡線信号の一部を選択し、前記第１信号と前記包絡線信号との間の一致度を算出するように構成された手段を備えることを特徴とする装置。
2. 請求項１に記載の装置において、前記一致度の算出は、前記第１信号と前記包絡線信号の相互相関関数の計算として具現化されることを特徴とする装置。
3. 請求項１または２に記載の装置において、前記第１信号は、前記包絡線信号の一連のサイクルの１つを表すように選択されることを特徴とする装置。
4. 請求項２に記載の装置において、前記推定装置が備える前記手段は、前記相互相関関数の計算を行うための前記第１信号と前記包絡線信号として、同じサンプル・シーケンスを選択し、これによって、前記包絡線信号に対する自己相関関数を計算するように構成されていることを特徴とする装置。
5. 請求項１ないし４のいずれか１項に記載の装置であって、前記聴診信号の品質を検証するために、前記聴診信号の前記一致度を表す信号において、以下の３つの項目：
   ａ）前記一致度を表す信号の、位置が検出された極値間の時間差は、所定の制限内でなければならない；
   ｂ）前記時間差の平均値に対する時間差の最小値および最大値の比は、所定の制限内でなければならない；
   ｃ）前記一致度を表す信号の前記極値の位置における相関の結果の大きさは、所定の制限内でなければならない；
   のうち少なくとも１つを確認する手段（４２１）を、前記推定装置が備えることを特徴とする装置。
6. 請求項１ないし５のいずれか１項に記載の装置において、前記処理手段（１０６；３０１；４０９）は、前記聴診信号のＡ特性周波数重み付けバージョン、または前記聴診信号の近似的なＡ特性周波数重み付けバージョンを算出するためのフィルタ（１０１；４０１）を備えることを特徴とする装置。
7. 請求項６に記載の装置において、前記聴診信号の前記Ａ特性周波数重み付けバージョンは、前記聴診信号の二重微分に対応する近似によって算出されることを特徴とする装置。
8. 請求項１ないし７のいずれか１項に記載の装置において、前記処理手段（４０９）は、前記包絡線信号を濾波するための適応バンド・パス・フィルタ（５１２）を備えており、前記適応バンド・パス・フィルタは、少なくとも、それぞれ選択可能な高い通過帯域および低い通過帯域を有し、前記適応バンド・パス・フィルタは、前記適応バンド・パス・フィルタに入力される信号が相対的に低周波数である場合は前記低い通過帯域を選択すると共に前記適応バンド・パス・フィルタに入力される信号が相対的に高周波数である場合は前記高い通過帯域を選択するコントローラ（５１３）を備えることを特徴とする装置。
9. 請求項１ないし８のいずれか１項に記載の装置において、前記聴診信号は、サンプル・レートで到着するサンプルを含み、前記処理手段は、前記サンプル・レートに対応するレートで動作する前記聴診信号の前記包絡線を供給するための同期プロセッサ（３０１；４０９）を備え、前記推定装置は、要求によって開始する時点で動作する前記一致度算出を実行するための非同期プロセッサ（３０２；４１０）も備えることを特徴とする装置。
10. 請求項１ないし９のいずれか１項に記載の装置において、前記包絡線信号は、前記聴診信号の絶対値を表す信号を低域濾波するための第１分岐の手段（４０６、４０７）によって供給されると共に、前記包絡線信号は、前記聴診信号のＡ特性周波数重み付けバージョンまたは前記聴診信号の近似的なＡ特性周波数重み付けバージョンを低域濾波することによって得られる信号の絶対値を表す信号を低域濾波するための第２分岐の手段（４０１、４０２、４０３、４０４）によって供給されることを特徴とする装置。
11. 請求項１０に記載の装置において、前記包絡線信号は、前記第１および第２分岐によって供給される前記信号の和を算出するための手段（４１８）によって供給されることを特徴とする装置。
12. 請求項１ないし１１のいずれか１項に記載の装置と、聴診に利用される音を検出して当該装置に聴診信号を供給するための手段とを備えることを特徴とする聴診器。
13. 聴診信号における律動を推定するための装置であって：
    聴診信号を受信し、前記聴診信号の包絡線を表す包絡線信号を供給するための処理手段（１０６；３０１；４０９）と；
    前記聴診信号の推定律動を表す信号を算出するための推定装置（１０５；３０２；４１０）と；
    を備える装置において、
    前記処理手段（１０６；３０１；４０９）は、Ａ特性周波数重み付けまたは近似的なＡ特性周波数重み付けに対応する周波数応答を有するフィルタ（１０１；４０１）を備えており、前記フィルタは、少なくとも対象の周波数範囲について算出された二重微分として具現化されることを特徴とする装置。
14. 請求項１３に記載の装置において、前記対象の周波数範囲は２０００Ｈｚ未満の周波数範囲であることを特徴とする装置。
15. 聴診信号を処理する方法であって：
    聴診信号を受信し、前記聴診信号の包絡線を表す包絡線信号を供給するステップと；
    前記聴診信号の推定律動を表す信号を算出するステップと；
    を備える方法において、更に：
    第１信号として前記包絡線信号における一連のサイクルの１つを少なくとも含む前記包絡線信号の一部を選択するステップと、
    前記第１信号と前記包絡線信号との間の一致度を算出するステップと；
    を備えることを特徴とする方法。
16. 請求項１５に記載の方法において、前記一致度を算出するステップは、前記第１信号と前記包絡線信号の相互相関関数によって算出されることを特徴とする方法。
17. 請求項１５または１６に記載の方法において、前記包絡線の前記一部は、前記包絡線信号の一連のサイクルの１つを表すように選択されることを特徴とする方法。
18. 請求項１６に記載の方法であって、前記相互相関関数の計算を行うための前記第１信号と前記包絡線信号として、同じサンプル・シーケンスを選択し、これによって、前記包絡線信号に対する自己相関関数を計算するステップを備えることを特徴とする方法。
19. 請求項１５ないし１８のいずれか１項に記載の方法であって、前記聴診信号の前記一致度を表す信号において、以下の３つの項目：
    ａ）前記一致度を表す信号の、位置が検出された極値間の時間差は、所定の制限内でなければならない；
    ｂ）前記時間差の平均値に対する時間差の最小値および最大値の比は、所定の制限内でなければならない；
    ｃ）前記一致度を表す信号の前記極値の位置における相関の結果の大きさは、所定の制限内でなければならない；
    のうち少なくとも１つを確認することによって、前記聴診信号の品質を検証するステップを備えることを特徴とする方法。
20. 請求項１５ないし１９のいずれか１項に記載の方法であって、更に、前記聴診信号のＡ特性周波数重み付けバージョン、または前記聴診信号の近似的なＡ特性周波数重み付けバージョンを算出するステップを備えることを特徴とする方法。
21. 請求項２０に記載の方法において、前記聴診信号の前記Ａ特性周波数重み付けバージョンは、前記聴診信号の二重微分に対応する近似によって算出されることを特徴とする方法。
22. 請求項１５ないし２１のいずれか１項に記載の方法であって、少なくとも、それぞれ選択可能な高い通過帯域および低い通過帯域を有する適応バンド・パス・フィルタによって前記包絡線信号を濾波するステップをさらに備え、前記適応バンド・パス・フィルタは、前記適応バンド・パス・フィルタに入力される信号が相対的に低周波数である場合は前記低い通過帯域を選択すると共に前記適応バンド・パス・フィルタに入力される信号が相対的に高周波数である場合は前記高い通過帯域を選択するように制御されることを特徴とする方法。
23. 請求項１５ないし２２のいずれか１項に記載の方法において、前記聴診信号は、サンプル・レートで到着するサンプルを含み、前記包絡線信号を供給するステップは、前記サンプル・レートに対応するレートで実行している前記聴診信号の前記包絡線を供給するための同期ステップを備え、前記聴診信号の推定律動を表す信号を算出するステップは、要求によって開始する時点で動作する前記一致度算出を実行するための非同期ステップも備えることを特徴とする方法。
24. 請求項１５ないし２３のいずれか１項に記載の方法において、前記包絡線信号は、前記聴診信号の絶対値を表す信号を低域濾波する第１ステップによって供給されると共に、前記包絡線信号は、前記聴診信号のＡ特性周波数重み付けバージョンまたは前記聴診信号の近似的なＡ特性周波数重み付けバージョンを低域濾波することによって得られる信号の絶対値を表す信号を低域濾波する第２ステップによって供給されることを特徴とする方法。
25. 請求項２４に記載の方法であって、前記第１および第２ステップによって供給される前記信号を加算するステップを備えることを特徴とする方法。
26. 記録されたプログラムを有するコンピュータ読み取り可能媒体において、前記プログラムは、実行されると、請求項１５ないし２５のいずれか１項に記載の方法を前記コンピュータに実行させるコンピュータ読み取り可能媒体。
27. 聴診信号における律動を推定するための方法であって：
    聴診信号を受信し、前記聴診信号の包絡線を表す包絡線信号を供給するステップと；
    前記聴診信号の推定律動を表す信号を算出するステップと；
    を備える方法において、
    前記推定律動を表す前記信号は、Ａ特性周波数重み付けまたは近似的なＡ特性周波数重み付けに対応する周波数応答を有するフィルタ（１０１；４０１）によって前記聴診信号を濾波することによって算出され、前記フィルタは、少なくとも対象の周波数範囲についての二重微分の計算として具現化されることを特徴とする方法。
28. 請求項２７に記載の方法において、前記対象の周波数範囲は２０００Ｈｚ未満の周波数範囲であることを特徴とする方法。
29. 記録されたプログラムを有するコンピュータ読み取り可能媒体において、前記プログラムは、実行されると、請求項２７または２８に記載の方法を前記コンピュータに実行させるコンピュータ読み取り可能媒体。

Images