JP6630149B2 - 生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体 - Google Patents

Description

本発明は、例えば断続性ラ音を含む生体音を解析する生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体の技術分野に関する。

この種の装置として、電子聴診器等によって検出される生体の呼吸音について、そこに含まれる副雑音（即ち、正常な呼吸音とは異なる音）の音種を判別するものが知られている。例えば特許文献１では、生体音の音波形を複数の区分に分割し、区分毎に音波形の特徴を特定して音種を判別するという手法が提案されている。

特開２０１３−１２３４９４号公報

上述した特許文献１に記載されている技術では、高周波帯域の周波数成分が多く認められるか否かによって、細かな断続性ラ音と荒い断続性ラ音とを判別している。しかしながら、このような手法では、例えば周囲の環境音の影響等が入る場合において正確に音種を判別することができなくなるおそれがある。即ち、周波数成分だけ判別基準として用いる場合、状況によっては、必ずしも正確な結果が得られないおそれがある。

また、特許文献１に記載されている技術では、複数種類の断続性ラ音（例えば、捻髪音及び水泡音）が同時に発生している場合、それらの両方が発生しているか否かを判別することができない。

本発明が解決しようとする課題には、上記のようなものが一例として挙げられる。本発明は、生体音に含まれる雑音を好適に解析可能な生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体を提供することを課題とする。

上記課題を解決するための生体音解析装置は、生体音の経時的な変化を示す生体音情報を取得する第１取得手段と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得手段と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力手段とを備える。

上記課題を解決するための生体音解析方法は、生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力工程とを備える。

上記課題を解決するためのコンピュータプログラムは、生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力工程とをコンピュータに実行させる。

上記課題を解決するための記録媒体は、上述したコンピュータプログラムが記録されている。

本実施例に係る生体音解析装置の構成を示すブロック図である。 捻髪音を含む生体音の一例を示すスペクトログラム図である。 解析処理後の捻髪音を含む生体音の一例を示すスペクトログラム図である。 水泡音を含む生体音の一例を示すスペクトログラム図である。 解析処理後の水泡音を含む生体音の一例を示すスペクトログラム図である。 本実施例に係る生体音解析装置の動作の流れを示すフローチャートである。 捻髪音の２次元分布を示すモデル図である。 水泡音の２次元分布を示すモデル図である。 捻髪音のパワーの確率分布を示すモデル図である。 水泡音のパワーの確率分布を示すモデル図である。 呼吸音情報のパワー分布比率を区画ごとに示す図表である。

＜１＞
本実施形態に係る生体音解析装置は、生体音の経時的な変化を示す生体音情報を取得する第１取得手段と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得手段と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力手段とを備える。

本実施形態に係る生体音解析装置によれば、その動作時には、先ず第１取得手段によって、生体音の経時的な変化を示す生体音情報が取得される。また、第２取得手段によって、生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報が取得される。なお、「生体音」とは、生体が発する音であり、典型的には呼吸音である。「雑音」とは、本来であれば生体音に含まれるべきではない異音であり、例えば異常呼吸音（副雑音）である断続性ラ音等が挙げられる。また「所定期間」とは、雑音の発生時期及び周波数特性の特徴が顕著に現れる期間であり、例えば１呼吸周期に対応する期間である。

生体音情報及び参照情報が取得されると、それらの情報に基づいて、出力手段から雑音を示す雑音情報が出力される。具体的には、出力手段では生体音情報と参照情報とが比較され、生体音に雑音が含まれているか否かが判定されると共に、雑音と判定された音に関する情報が雑音情報として出力される。なお、出力される雑音情報は複数種類であっても構わない。複数種類の雑音情報を出力するためには、第２手段において、参照情報が複数種類（即ち、雑音毎に）取得されるようにすればよい。

本実施形態では特に、雑音情報を出力するために用いられる参照情報が、雑音の所定期間内における発生時期及び周波数特性を含む情報として取得されている。このため、生体音に含まれる雑音を、所定期間内における発生時期及び周波数特性の観点から精度よく判別することができる。具体的には、所定の周波数範囲において、所定期間内における一の期間でのみで発生するような雑音等を正確に判別することが可能となる。生体音に含まれる雑音の中には、発生時期及び周波数特性に顕著な特徴を有するものがあることが一般的に知られている。よって、雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得することで、極めて好適に雑音情報を出力することが可能となる。

なお、生体音に含まれる雑音は、周波数特性のみに基づいて判別することもできると考えられる。しかしながら、周波数特性のみを基準とする判別方法では、例えば周囲の環境音の影響等が入る場合において、正確な判別が行えなくなる可能性がある。即ち、判別基準が周波数特性のみであるが故に、予期せぬノイズ等により判別制度が大きく低下してしまうおそれがある。しかるに本実施形態では、既に説明したように、所定期間内における発生時期及び周波数特性を含む参照情報を利用して生体音に含まれる雑音を示す情報を出力することができる。従って、生体音を好適に解析して、正確な雑音情報を得ることが可能である。

＜２＞
本実施形態に係る生体音解析装置の一態様では、前記生体音情報に基づいて、前記生体音に断続性ラ音が含まれているか否かを判定する判定手段を更に備え、前記出力手段は、前記断続性ラ音が含まれていると判定された前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する。

この態様によれば、先ず判定手段によって生体音に断続性ラ音が含まれているか否かが判定される。なお、判定手段は、既存の技術を利用して生体音に断続性ラ音が含まれているか否かを判定することができる。言い換えれば、判定手段による断続性ラ音の判定手法は、特に限定されるものではない。

判定手段によって断続性ラ音が判定されると、断続性ラ音が含まれていると判定された生体音情報及び参照情報に基づいて、雑音情報が出力される。即ち、この態様では、雑音情報を出力する対象となる生体音が、断続性ラ音である雑音を含む生体音に限られる。このため、出力すべき雑音情報が予め限定されることになり、より正確に雑音情報を出力することが可能となる。

＜３＞
本実施形態に係る生体音解析装置の他の態様では、前記生体音情報を前記所定期間に対応づける時間軸情報を取得する第３取得手段を更に備え、前記出力手段は、前記所定期間に対応づけられた前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する。

この態様によれば、第３取得手段によって、生体音情報を所定期間に対応づける時間軸情報が取得される。これにより、生体音情報と所定期間とを好適に対応づけることが可能となり、参照情報（即ち、雑音の所定期間内における発生時期及び周波数特性）に基づき、より正確な雑音情報を出力することが可能となる。

＜４＞
上述した第３取得手段を備える態様では、前記所定期間は、１呼吸周期に対応する期間であり、前記第３取得手段は、吸気相及び呼気相を含む呼吸相を前記時間軸情報として取得してもよい。

この場合、１呼吸周期における発生時期及び周波数特性に特徴のある雑音に関して、極めて正確な雑音情報を出力することが可能となる。例えば、吸気相においてのみ高周波数で発生するような雑音や、吸気相及び呼気相の両方で発生するが、その多くが低周波数で発生する雑音等を示す雑音情報を、正確に判別して出力することができる。

＜５＞
上述した１呼吸周期に対応する期間を所定期間とする態様では、前記第２取得手段は、前記発生時期が前記吸気相であること、前記周波数特性が第１周波数よりも高いことを、捻髪音に対応する前記参照情報として取得してもよい。

断続性ラ音の一種である捻髪音は、呼吸相の吸気相において第１周波数よりも高い音として発生する。よって、これに対応する参照情報を取得することで、生体音に含まれる捻髪音の雑音情報を極めて高い精度で出力することが可能となる。

＜６＞
或いは１呼吸周期に対応する期間を所定期間とする態様では、前記第２取得手段は、前記発生時期が前記吸気相及び呼吸相の両方に及ぶこと、前記周波数特性が第２周波数よりも低いことを、水泡音に対応する前記参照情報として取得してもよい。

断続性ラ音の一種である水泡音は、呼吸相の吸気相及び呼気相の両方において、第２周波数よりも低い音として発生する。よって、これに対応する参照情報を取得することで、生体音に含まれる水泡音の雑音情報を極めて高い精度で出力することができる。

＜７＞
本実施形態に係る生体音解析装置の他の態様では、前記参照情報は、前記発生時期及び前記周波数特性の２軸上で表された前記雑音情報の分布情報を含み、前記出力手段は、前記分布情報及び前記生体音情報の一致度合いに基づいて、前記雑音情報を出力する。

この態様によれば、参照情報には、発生時期及び周波数特性の２軸上で表された雑音情報の分布情報が含まれている。このため、分布情報と生体音情報の一致度合いを利用すれば、生体音情報がどの程度雑音情報に近い傾向を有しているのかを容易に知ることができる。この結果、例えば一致度合いが所定の閾値を超えているか否かを判定することにより、好適に参照情報に対応する雑音か否かを判別することができる。従って、より正確な雑音情報を出力することが可能となる。

＜８＞
本実施形態に係る生体音解析方法は、生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力工程とを備える。

本実施形態に係る生体音解析方法によれば、上述した本実施形態に係る生体音解析装置と同様に、正確な雑音情報を出力することができる。

なお、本実施形態に係る生体音解析方法においても、上述した本実施形態に係る生体音解析装置における各種態様と同様の各種態様を採ることが可能である。

＜９＞
本実施形態に係るコンピュータプログラムは、生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、前記生体音に含まれ得る雑音の所定期間内における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、前記生体音情報及び前記参照情報に基づいて、前記生体音に含まれる前記雑音を示す雑音情報を出力する出力工程とをコンピュータに実行させる。

本実施形態に係るコンピュータプログラムによれば、上述した本実施形態に係る生体音解析方法と同様の処理をコンピュータに実行させることができるため、正確な雑音情報を出力することができる。

なお、本実施形態に係るコンピュータプログラムにおいても、上述した本実施形態に係る生体音解析装置における各種態様と同様の各種態様を採ることが可能である。

＜１０＞
本実施形態に係る記録媒体は、上述したコンピュータプログラムが記録されている。

本実施形態に係る記録媒体によれば、上述したコンピュータプログラムをコンピュータにより実行させることにより、正確な雑音情報を出力することができる。

本実施形態に係る生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体の作用及び他の利得については、以下に示す実施例において、より詳細に説明する。

以下では、生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体の実施例について、図面を参照しながら詳細に説明する。なお、以下では、呼吸音の解析を行う生体音解析装置を例に挙げて説明する。

＜装置構成＞
先ず、本実施例に係る生体音解析装置の構成について、図１を参照して説明する。ここに図１は、本実施例に係る生体音解析装置の構成を示すブロック図である。

図１において、本実施例に係る生体音解析装置は、呼吸音取得部１１０と、参照情報取得部１２０と、呼吸相取得部１３０と、処理部２００と、結果表示部３００とを備えて構成されている。

呼吸音取得部１１０は、生体の呼吸音を呼吸音信号として取得可能に構成されたセンサである。呼吸音取得部１１０は、例えばＥＣＭ（Electret Condenser Microphone）やピエゾを利用したマイク、振動センサ等で構成されている。また、呼吸音取得部１１０は、生体の呼吸音を呼吸音信号として取得可能に構成されたセンサだけでなく、センサからの呼吸音信号を取得するものを含んでもよい。呼吸音取得部１１０で取得された呼吸音信号は、時間周波数解析部２１０及び断続性ラ音判定部２２０に出力される構成となっている。なお、呼吸音取得部１１０は、「第１取得手段」の一具体例である。

参照情報取得部１２０は、生体音に含まれる雑音を判別するための参照情報を記憶するデータベースを含んで構成されている。参照情報取得部１２０は、記憶された参照情報を適宜取得して、捻髪音傾向算出部２３０及び水泡音傾向算出部２４０に出力可能に構成されている。本実施例に係る参照情報取得部１２０は特に、断続性ラ音である捻髪音及び水泡音に関する参照情報を取得可能とされている。なお、参照情報取得部１２０は、「第２取得手段」の一具体例である。

呼吸相取得部１３０は、生体の呼吸相情報（即ち、吸気相及び呼気相に関する時間軸情報）を取得可能に構成されたセンサである。呼吸相取得部１３０は、呼吸音取得部１１０と一体的に構成されていてもよい。呼吸相取得部１３０で取得された呼吸相情報は、捻髪音傾向算出部２３０及び水泡音傾向算出部２４０に出力される構成となっている。なお、呼吸相取得部１３０は、「第３取得手段」の一具体例である。

処理部２００は、複数の演算回路やメモリ等を含んで構成されている。処理部２００は、時間周波数解析部２１０と、断続性ラ音判定部２２０と、捻髪音傾向算出部２３０と、水泡音傾向算出部２４０と、捻髪音判定部２５０と、水泡音判定部２６０とを備えて構成されている。

時間周波数解析部２１０は、呼吸音取得部１１０で取得された呼吸音情報に対して時間周波数解析処理を実行する。具体的には、時間周波数解析部２１０は、ＦＦＴ（Fast Fourier Transform：高速フーリエ変換）処理等を実行可能に構成されている。時間周波数解析部２１０の解析結果は、捻髪音傾向算出部２３０及び水泡音傾向算出部２４０に出力される構成となっている。

断続性ラ音判定部２２０は、呼吸音取得部１１０で取得された呼吸音情報に基づいて、呼吸音に断続性ラ音が含まれているか否かを判定する。なお、断続性ラ音判定部２２０における具体的な判定方法については、既存の技術を利用することができるため、ここでの詳細な説明は省略する。断続性ラ音判定部２２０における判定結果は、捻髪音傾向算出部２３０及び水泡音傾向算出部２４０に出力される構成となっている。

捻髪音傾向算出部２３０は、時間周波数解析部２１０から入力される時間周波数解析された呼吸音信号、参照情報取得部１２０から入力される捻髪音に関する参照情報、呼吸相取得部１３０から入力される呼吸相情報、及び断続性ラ音判定部２２０から入力される判定結果に基づいて、生体音に含まれる音が捻髪音にどの程度近いものであるかを示す捻髪音傾向ＦＣを算出する。なお、捻髪音傾向ＦＣの算出方法については、後の動作説明において詳述する。

水泡音傾向算出部２４０は、時間周波数解析部２１０から入力される時間周波数解析された呼吸音信号、参照情報取得部１２０から入力される水泡音に関する参照情報、呼吸相取得部１３０から入力される呼吸相情報、及び断続性ラ音判定部２２０から入力される判定結果に基づいて、生体音に含まれる音が水泡音にどの程度近いものであるかを示す水泡音傾向ＣＣを算出する。なお、水泡音傾向ＣＣの算出方法については、後の動作説明において詳述する。

捻髪音判定部２５０は、捻髪音傾向算出部２３０で算出された捻髪音傾向ＦＣに基づいて、生体音に捻髪音が含まれているか否かを判定する。具体的には、捻髪音判定部２５０は、捻髪音傾向ＦＣが所定の閾値を超えているか否かによって、生体音に捻髪音が含まれているか否かを判定する。捻髪音判定部２５０の判定結果は、結果表示部３００に出力される構成となっている。

水泡音判定部２６０は、水泡音傾向算出部２４０で算出された水泡音傾向ＣＣに基づいて、生体音に水泡音が含まれているか否かを判定する。具体的には、水泡音判定部２６０は、水泡音傾向ＣＣが所定の閾値を超えているか否かによって、生体音に水泡音が含まれているか否かを判定する。水泡音判定部２６０の判定結果は、結果表示部３００に出力される構成となっている。

以上のように、処理部２００は、呼吸音取得部１１０で取得された呼吸音情報、参照情報取得部１２０で取得された参照情報、呼吸相取得部１３０で取得された呼吸相情報に基づいて、生体音に捻髪音又は水泡音が含まれているか否かを判定することが可能とされている。また、処理部２００は、生体音に捻髪音又は水泡音が含まれているか否かだけでなく、捻髪音及び水泡音の強度等を出力可能に構成されてもよい。処理部２００は、「出力手段」の一具体例である。

結果表示部４００は、例えば液晶モニタ等のディスプレイとして構成されており、処理部２００から出力される各種情報を画像データとして表示する。

＜捻髪音と水泡音の特性＞
次に、本実施例に係る生体音解析装置において判別される捻髪音及び水泡音の特性について、図２から図５を参照して詳細に説明する。ここに図２は、捻髪音を含む生体音の一例を示すスペクトログラム図であり、図３は、解析処理後の捻髪音を含む生体音の一例を示すスペクトログラム図である。また図４は、水泡音を含む生体音の一例を示すスペクトログラム図であり、図５は、解析処理後の水泡音を含む生体音の一例を示すスペクトログラム図である。

なお、図３に示すスペクトログラムは、図２に示すスペクトログラムに対してＣＭＮ（Cepstral Mean Normalization）処理と、リフタリング処理（具体的には、ケプストラムの高次ケフレンシ―成分のカットとを実行したものである。同様に、図４に示すスペクトログラムは、図３に示すスペクトログラムに対してＣＭＮ処理と、リフタリング処理とを実行したものである。ＣＭＮ処理及びリフタリング処理を実行することで、捻髪音及び水泡音の特徴が強調され、その特性がより明確に判別できるようになる。なお、ＣＭＮ処理及びリフタリング処理については既存の技術であるため、ここでのより詳細な説明は省略する。

図２及び図３を見ると、捻髪音は、吸気相において強く検出され、呼気相においては殆ど検出されないことが分かる。また、捻髪音は、比較的高い周波数帯で発生していることも分かる。よって、捻髪音は、吸気相において比較的高めの周波数帯に現れる特性を有していると言える。

一方、図４及び図５を見ると、水泡音は、吸気相及び呼気相の両方において検出されていることが分かる。また、水泡音は、比較的低い周波数帯で発生していることも分かる。よって、水泡音は、吸気相及び呼気相の両方において比較的低めの周波数帯に現れる特性を有していると言える。

以上のように、捻髪音及び水泡音を含む生体音には、呼吸相における発生時期と周波数特性において明確な違いが存在している。よって、この違いを利用すれば、捻髪音及び水泡音の存在を個別に判定することが可能である。本実施形態に係る生体音解析装置は、以下に詳述する処理を実行することにより、生体音に含まれる捻髪音及び水泡音を判別する。

＜動作説明＞
次に、本実施例に係る生体音解析装置の動作について、図６を参照して説明する。ここに図６は、本実施例に係る生体音解析装置の動作の流れを示すフローチャートである。

図６において、本実施例に係る生体音解析装置の動作時には、先ず呼吸音取得部１１０において、生体の呼吸音を示す呼吸音信号が取得される（ステップＳ１０１）。呼吸音信号が取得されると、断続性ラ音判定部２２０において、生体音に断続性ラ音が含まれているか否かの判定が行われる（ステップＳ１０２）。判定の結果、生体音に断続性ラ音が含まれていないと判定された場合（ステップＳ１０３：ＹＥＳ）、断続性ラ音判定部２２０からは、生体音には捻髪音及び水泡音が含まれないという判定結果が出力される（ステップＳ１０４）。一方、生体音に断続性ラ音が含まれていると判定された場合（ステップＳ１０３：ＹＥＳ）、生体音に含まれる断続性ラ音が捻髪音であるのか、或いは水泡音であるのかを判別するための処理（即ち、ステップＳ１０５以降の処理）が実行される。

生体音に断続性ラ音が含まれている場合には、先ず時間周波数解析部２２０において、呼吸音信号の時間周波数解析処理が行われる（ステップＳ１０５）。また、参照情報取得部１２０において、捻髪音及び水泡音に関する参照情報が取得される（ステップＳ１０６）。更に、呼吸相取得部１３０において、生体音の呼吸相情報が取得される（ステップＳ１０７）。

上述した各情報が取得されると、捻髪音傾向算出部２３０において、捻髪音傾向ＦＣが算出される（ステップＳ１０８）。算出された捻髪音傾向ＦＣは、捻髪音判定部２５０において閾値Ａよりも大きいか否かが判定される（ステップＳ１０９）。そして、捻髪音傾向ＦＣが閾値Ａよりも大きいと判定された場合（ステップＳ１０９：ＹＥＳ）、生体音には捻髪音が含まれていると判定される（ステップＳ１１０）。一方、捻髪音傾向ＦＣが閾値Ａよりも大きくないと判定された場合（ステップＳ１０９：ＮＯ）、生体音には捻髪音が含まれていないと判定される（ステップＳ１１１）。

同様に、水泡音傾向算出部２４０においては、水泡音傾向ＣＣが算出される（ステップＳ１１２）。算出された水泡音傾向ＣＣは、水泡音判定部２６０において閾値Ｂよりも大きいか否かが判定される（ステップＳ１１３）。そして、水泡音傾向ＣＣが閾値Ｂよりも大きいと判定された場合（ステップＳ１１３：ＹＥＳ）、生体音には水泡音が含まれていると判定される（ステップＳ１１４）。一方、水泡音傾向ＣＣが閾値Ｂよりも大きくないと判定された場合（ステップＳ１１３：ＮＯ）、生体音には水泡音が含まれていないと判定される（ステップＳ１１５）。

以上説明した、ステップＳ１０４、Ｓ１１４及びＳ１１５の判定結果は、結果表示部３００に出力される（ステップＳ１１６）。これにより、結果表示部３００では、呼吸音に捻髪音及び水泡音が含まれているか否かが画像データとして表示される。

＜傾向算出方法＞
次に、上述した捻髪音傾向ＦＣ及び水泡音傾向ＣＣの算出方法について、図７から図１１を参照して具体的に説明する。ここに図７は、捻髪音の２次元分布を示すモデル図であり、図８は、水泡音の２次元分布を示すモデル図である。また図９は、捻髪音のパワーの確率分布を示すモデル図であり、図１０は、水泡音のパワーの確率分布を示すモデル図である。図１１は、呼吸音情報のパワー分布比率を区画ごとに示す図表である。

図７において、捻髪音は吸気相において比較的高めの周波数帯に現れる特性を有しているため、この特性を呼吸相における発生時期及び周波数特性の２軸上で表すと、図のような２次元モデルが得られる。捻髪音傾向ＦＣは、このような２次元モデルを参照情報として用いることで算出することができる。

図８において、水泡音は吸気相及び呼気相の両方において比較的低めの周波数帯に現れる特性を有しているため、この特性を呼吸相における発生時期及び周波数特性の２軸上で表すと、図のような２次元モデルが得られる。水泡音傾向ＣＣは、このような２次元モデルを参照情報として用いることで算出することができる。

図９において、捻髪音傾向ＦＣの算出時には、先ず捻髪音に対応する２次元モデル（図７参照）のパワー分布を区画ごとに数値化する。具体的には、図に示すように、高周波数帯（７５０Ｈｚ〜４ｋＨＺ）の吸気相数値Ｑ_ＦＣ（１）が「０．７」となり、呼気相の数値Ｑ_ＦＣ（２）が「０．１」となる。また、低周波数帯（０Ｈｚ〜７５０Ｈｚ）の吸気相の数値Ｑ_ＦＣ（３）が「０．１」となり、呼気相の数値Ｑ_ＦＣ（４）が「０．１」となる。

捻髪音傾向ＦＣは、上記Ｑ_ＦＣ（ｉ）を利用して、下記の数式（１）から算出することができる。

なお、Ｐ（ｉ）は、実際に計測された呼吸音信号の区画ごとのパワーを示す数値である。

このようにして算出される捻髪音傾向ＦＣは、対数尤度であり、負の値として算出される。このため、捻髪音傾向ＦＣが大きいほど（即ち、０に近い値であるほど）、生体音がモデルに近い（即ち、捻髪音に近い）と判定できる。

一方、図１０において、水泡音傾向ＣＣの算出時には、先ず水泡音に対応する２次元モデル（図８参照）のパワー分布を区画ごとに数値化する。具体的には、図に示すように、高周波数帯（７５０Ｈｚ〜４ｋＨＺ）の吸気相の数値Ｑ_ＣＣ（１）が「０．１」となり、呼気相の数値Ｑ_ＣＣ（２）が「０．１」となる。また、低周波数帯（０Ｈｚ〜７５０Ｈｚ）の吸気相の数値Ｑ_ＣＣ（３）が「０．４」となり、呼気相の数値Ｑ_ＣＣ（４）が「０．４」となる。

水泡音傾向ＣＣは、上記Ｑ_ＣＣ（ｉ）を利用して、下記の数式（２）から算出することができる。

このようにして算出される水泡音傾向ＣＣは、捻髪音傾向ＦＣと同様に対数尤度であり、負の値として算出される。このため、水泡音傾向ＣＣの絶対値が大きいほど（即ち、０に近い値であるほど）、呼吸音がモデルに近い（即ち、捻髪音に近い）と判定できる。

以下では、上述した算出処理及びその算出結果を用いた判定処理について、具体例を挙げて説明する。

図１１において、例えば図に示すような比率でパワーが分布する呼吸音信号が得られたとする。この呼吸音信号について、上述した数式（１）を利用して捻髪音傾向ＦＣを算出すると、その結果は以下のようになる。

ＦＣ＝０．８２×ｌｏｇ_１０（０．７）＋０．００×ｌｏｇ１０（０．１）＋０．１６×ｌｏｇ_１０（０．１）＋×０．０２ｌｏｇ１０（０．１）
＝−０．３０７０２
この時、捻髪音傾向ＦＣに対する閾値Ａが−０．５であったとすると、ＦＣ＞閾値Ａの関係が成立する。よって、この場合の呼吸音信号には、捻髪音が含まれると判定されることになる。

＜実施例の効果＞
最後に、本実施例に係る生体音解析装置によって得られる技術的効果について詳細に説明する。

本実施例に係る生体音解析装置によれば、図７から図１１で説明したように、呼吸相における発生時期及び周波数特性の２軸上で表された２次元モデル（即ち、参照情報取得部１２０から取得される参照情報）を利用して、断続性ラ音である捻髪音及び水泡音が判定される。

ここで仮に、周波数特性のみで捻髪音及び水泡音を判定しようとすると、周波数が高いか低いかでしか判定できないため、呼気相で発生している高周波数帯の成分を捻髪音として判定してしまうおそれがある。また、吸気相でしか発生していない低周波数帯の成分を水泡音と判定してしまうおそれがある。このように、周波数特性だけを利用する場合には、捻髪音及び水泡音を正確に判別できないおそれがある。特に、周囲の環境音等の影響を受けてしまうと、判別精度は大きく低下してしまう。

これに対し、本実施例では、呼吸相のどこで発生しているのかを考慮して、捻髪音であるのか、或いは水泡音であるのかを判定することができる。従って、周波数特性のみを利用して判定する場合と比べると、極めて正確に捻髪音及び水泡音を判定することができる。また、捻髪音傾向ＦＣ及び水泡音傾向ＣＣを別々に算出して判定しているため、呼吸音に捻髪音のみが含まれている場合、水泡音のみが含まれている場合、或いは捻髪音及び水泡音の両方が含まれている場合を、それぞれ好適に判別することが可能である。

本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う生体音解析装置及び生体音解析方法、並びにコンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。

１１０ 呼吸音取得部
１２０ 参照情報取得部
１３０ 呼吸相取得部
２００ 処理部
２１０ 時間周波数解析部
２２０ 断続性ラ音判定部
２３０ 捻髪音傾向算出部
２４０ 水泡音傾向算出部
２５０ 捻髪音判定部
２６０ 水泡音判定部
３００ 結果表示部
ＦＣ 捻髪音傾向
ＣＣ 水泡音傾向

Claims (10)

1. 生体音の経時的な変化を示す生体音情報を取得する第１取得手段と、
   前記生体音に含まれ得る雑音の１呼吸周期に対応する期間における発生時期及び周波数特性を含む参照情報を取得する第２取得手段と、
   前記生体音情報及び前記参照情報に基づいて、前記雑音に含まれる断続性ラ音を示す情報を出力する出力手段と
   を備えることを特徴とする生体音解析装置。
2. 前記生体音情報に基づいて、前記生体音に前記断続性ラ音が含まれているか否かを判定する判定手段を更に備え、
   前記出力手段は、前記断続性ラ音が含まれていると判定された前記生体音情報及び前記参照情報に基づいて、前記雑音に含まれる前記断続性ラ音を示す情報を出力する
   ことを特徴とする請求項１に記載の生体音解析装置。
3. 前記生体音情報を前記１呼吸周期に対応する期間に対応づける時間軸情報を取得する第３取得手段を更に備え、
   前記出力手段は、前記１呼吸周期に対応する期間に対応づけられた前記生体音情報及び前記参照情報に基づいて、前記雑音に含まれる前記断続性ラ音を示す情報を出力する
   ことを特徴とする請求項１又は２に記載の生体音解析装置。
4. 前記第３取得手段は、吸気相及び呼気相を含む呼吸相を前記時間軸情報として取得する
   ことを特徴とする請求項３に記載の生体音解析装置。
5. 前記第２取得手段は、前記発生時期が前記吸気相であること、前記周波数特性が第１周波数よりも高いことを、捻髪音に対応する前記参照情報として取得することを特徴とする請求項４に記載の生体音解析装置。
6. 前記第２取得手段は、前記発生時期が前記吸気相及び呼気相の両方に及ぶこと、前記周波数特性が第２周波数よりも低いことを、水泡音に対応する前記参照情報として取得することを特徴とする請求項４又は５に記載の生体音解析装置。
7. 前記参照情報は、前記発生時期及び前記周波数特性の２軸上で表された前記断続性ラ音を示す情報の分布情報を含み、
   前記出力手段は、前記分布情報及び前記生体音情報の一致度合いに基づいて、前記断続性ラ音を示す情報を出力する
   ことを特徴とする請求項１から６のいずれか一項に記載の生体音解析装置。
8. 生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、
   前記生体音に含まれ得る雑音の１呼吸周期に対応する期間における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、
   前記生体音情報及び前記参照情報に基づいて、前記雑音に含まれる断続性ラ音を示す情報を出力する出力工程と
   を備えることを特徴とする生体音解析方法。
9. 生体音の経時的な変化を示す生体音情報を取得する第１取得工程と、
   前記生体音に含まれ得る雑音の１呼吸周期に対応する期間における発生時期及び周波数特性を含む参照情報を取得する第２取得工程と、
   前記生体音情報及び前記参照情報に基づいて、前記雑音に含まれる断続性ラ音を示す情報を出力する出力工程と
   をコンピュータに実行させることを特徴とするコンピュータプログラム。
10. 請求項９に記載のコンピュータプログラムが記録されていることを特徴とする記録媒体。