(57)【要約】  呼吸音解析装置は、呼吸音のスペクトルを、呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解手段（１２０）と、第１分解手段によって分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解手段（１３０）と、第１分解手段及び第２分解手段によって分解されたスペクトルに基づいて、呼吸音に含まれる分解されたスペクトルの各々の割合に関する情報を出力する出力手段（１４０，１５０）とを備える。この呼吸音解析装置によれば、呼吸音に含まれる複数の音種を好適に分解することが可能である。
2016002004.tif 000002
【特許請求の範囲】
【請求項１】
呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解手段と、
前記第１分解手段によって分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解手段と、
前記第１分解手段及び前記第２分解手段によって分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力手段と
を備えることを特徴とする呼吸音解析装置。
【請求項２】
前記第１分解手段は、
前記呼吸音のスペクトルの第１所定特徴に対応する周波数に関する情報を取得する周波数情報取得手段と、
前記複数の基準スペクトルを、前記周波数に関する情報に応じてシフトさせ、周波数シフト基準スペクトルを取得するシフト手段と、
前記周波数シフト基準スペクトルに基づいて、前記呼吸音のスペクトルを分解する第１分解処理手段と
を有することを特徴とする請求項１に記載の呼吸音解析装置。
【請求項３】
前記第２分解手段は、
前記第１分解手段によって分解された少なくとも一部のスペクトルの第２所定特徴に対応する周波数を取得する周波数取得手段と、
前記第２所定特徴に対応する周波数が時系列的に連続するか否かに応じて、前記第１分解手段によって分解された少なくとも一部のスペクトルを分解する第２分解処理手段と
を有することを特徴とする請求項１又は２に記載の呼吸音解析装置。
【請求項４】
前記第２分解手段は、前記第２分解処理手段によって、前記第２所定特徴に対応する周波数が時系列的に連続するとして分解されたスペクトルを、前記第２所定特徴に対応する周波数と所定の閾値との関係に応じて更に分解する第３分解処理手段を有することを特徴とする請求項３に記載の呼吸音解析装置。
【請求項５】
前記第２分解手段は、前記第１分解手段によって分解された少なくとも一部のスペクトルの振幅値を取得する振幅値取得手段と、
前記振幅値が時系列的に連続するか否かに応じて、前記第１分解手段によって分解された少なくとも一部のスペクトルを分解する第４分解処理手段と
を有することを特徴とする請求項１から４のいずれか一項に記載の呼吸音解析装置。
【請求項６】
呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、
前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、
前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程と
を備えることを特徴とする呼吸音解析方法。
【請求項７】
呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、
前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、
前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程と
をコンピュータに実行させることを特徴とするコンピュータプログラム。
【請求項８】
請求項７に記載のコンピュータプログラムが記録されていることを特徴とする記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の音種を含む呼吸音を解析する呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の技術分野に関する。
【背景技術】
【０００２】
この種の装置として、電子聴診器等によって検出される生体の呼吸音について、正常音種と異常音種とを判別するものが知られている。例えば特許文献１では、副雑音の変換信号の時間／周波数表現に基づいて副雑音を識別するという技術が提案されている。特許文献２では、データベースに記憶されている異常音の音情報の中から最も類似する音情報を検索するという技術が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特表２００４−５３１３０９号公報
【特許文献２】国際公開２０１０／０４４４４２号
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、上述した特許文献１及び２に記載されているような技術では、正常音種と異常音種との分離を十分に行えないという技術的問題点がある。具体的には、例えば異常音種である笛声音、類鼾音、捻髪音、及び水泡音等を夫々分離することができない。【０００５】
本発明が解決しようとする課題には、上記のようなものが一例として挙げられる。本発明は、呼吸音に含まれる複数の音種を好適に分解することが可能な呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体を提供することを課題とする。
【課題を解決するための手段】
【０００６】
上記課題を解決するための呼吸音解析装置は、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解手段と、前記第１分解手段によって分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解手段と、前記第１分解手段及び前記第２分解手段によって分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力手段とを備える。【０００７】
上記課題を解決するための呼吸音解析方法は、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程とを備える。【０００８】
上記課題を解決するためのコンピュータプログラムは、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程とをコンピュータに実行させる。【０００９】
上記課題を解決するための記録媒体は、上述したコンピュータプログラムが記録されている。
【図面の簡単な説明】
【００１０】
【図１】本実施例に係る呼吸音解析装置の全体構成を示すブロック図である。
【図２】本実施例に係る呼吸音解析装置の動作を示すフローチャートである。
【図３】捻髪音を含む呼吸音の周波数解析結果を示すスペクトログラム図である。
【図４】笛声音を含む呼吸音の周波数解析結果を示すスペクトログラム図である。
【図５】捻髪音を含む呼吸音の所定タイミングにおけるスペクトルを示すグラフである。
【図６】捻髪音を含む呼吸音のスペクトルの近似方法を示す概念図である。
【図７】笛声音を含む呼吸音の所定タイミングにおけるスペクトルを示すグラフである。
【図８】笛声音を含む呼吸音のスペクトルの近似方法を示す概念図である。
【図９】周波数解析方法の一例を示すグラフである。
【図１０】周波数解析結果の一例を示す図である。
【図１１】スペクトルのピーク検出結果を示す概念図である。
【図１２】正常肺胞呼吸音基底を示すグラフである。
【図１３】捻髪音基底を示すグラフである。
【図１４】連続性ラ音基底を示すグラフである。
【図１５】ホワイトノイズ基底を示すグラフである。
【図１６】周波数シフトされた連続性ラ音基底を示すグラフである。
【図１７】スペクトルと、基底及び結合係数との関係を示す図である。
【図１８】観測されたスペクトル及び近似に用いられる基底の一例を示す図である。
【図１９】スペクトルを示す各基底及び結合係数を示す図である。
【図２０】ピーク周波数の時間的連続性を利用した分解方法を示す概念図（その１）である。
【図２１】ピーク周波数の時間的連続性を利用した分解方法を示す概念図（その２）である。
【図２２】第１実施例に係る笛声音と類鼾音との分解に用いる閾値を示すグラフである。
【図２３】第２実施例に係る笛声音と類鼾音との分解に用いる閾値の初期値を示すグラフである。
【図２４】第２実施例に係る笛声音と類鼾音との分解に用いる閾値の調整後の値を示すグラフ（その１）である。
【図２５】第２実施例に係る笛声音と類鼾音との分解に用いる閾値の調整後の値を示すグラフ（その２）である。
【図２６】笛声音を含む呼吸音のスペクトログラム図である。
【図２７】笛声音のピーク周波数及びピーク数を示すグラフである。
【図２８】類鼾音を含む呼吸音のスペクトログラム図である。
【図２９】類鼾音のピーク周波数及びピーク数を示すグラフである。
【図３０】振幅値の時間的連続性を利用した分解方法を示す概念図（その１）である。
【図３１】振幅値の時間的連続性を利用した分解方法を示す概念図（その２）である。
【図３２】分解処理の各段階において得られるスペクトログラムを、順を追って示す概念図（その１）である。
【図３３】分解処理の各段階において得られるスペクトログラムを、順を追って示す概念図（その２）である。
【図３４】表示部における表示例を示す平面図（その１）である。
【図３５】表示部における表示例を示す平面図（その２）である。
【発明を実施するための形態】
【００１１】
＜１＞
本実施形態に係る呼吸音解析装置は、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解手段と、前記第１分解手段によって分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解手段と、前記第１分解手段及び前記第２分解手段によって分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力手段とを備える。
【００１２】
本実施形態に係る呼吸音解析装置によれば、その動作時には、先ず第１分解手段によって呼吸音のスペクトルが分解される。第１分解手段では、呼吸音を分類する基準となる複数の基準スペクトルに基づいて、呼吸音のスペクトルが分解される。なお、ここでの「基準スペクトル」とは、呼吸音に含まれる複数の音種（例えば、正常呼吸音や、連続性ラ音、捻髪音等）を分類するために、各音種に応じて予め設定されたスペクトルである。基準スペクトルは、例えば各音種に特有の形状を有するスペクトルとして設定されている。
【００１３】
基準スペクトルを用いることで、各基準スペクトルに対応する音種が夫々どのような割合で呼吸音のスペクトルに含まれているかを知ることができる。即ち、呼吸音のスペクトルを基準スペクトルに対応する各音種に分解することができる。
【００１４】
第１分解手段による分解が行われた後には、第２分解手段によって、第１分解手段によって分解された少なくとも一部のスペクトルが分解される。即ち、第１分解手段によって分解されたスペクトルは、少なくとも部分的に更に分解されることになる。第２分解手段では、第１分解手段によって分解されたスペクトルが、所定の時系列特徴に基づいて分解される。ここで「所定の時系列特徴」とは、各スペクトルの時系列変化を判別するための特徴であり、例えば各音種に特有のピーク周波数や振幅値の時系列的な連続性を判別するためのものとして予め設定されている。なお、所定の時系列特徴は、複数種類設定されていてもよい。
【００１５】
所定の時系列特徴を用いれば、スペクトルが所定の時系列特徴を有する音種を含むものであるか否かを判別することができる。この結果、スペクトルを、所定の時系列特徴を有する音種のスペクトルと、所定の時系列特徴を有しない音種のスペクトルとに分解できる。
【００１６】
第２分解手段による分解が行われた後には、出力手段によって、第１分解手段及び第２分解手段によって分解されたスペクトルに基づいて、呼吸音に含まれる分解されたスペクトルの各々の割合に関する情報が出力される。即ち、第１分解手段及び第２分解手段による分解の結果として最終的に得られた各スペクトルについて、その各々が呼吸音に含まれる割合に関する情報が出力される。
【００１７】
ここで仮に、基準スペクトルを用いた分解（即ち、第１分解手段による分解）のみを利用するだけでも、呼吸音のスペクトルをある程度は分解することができる。しかしながら、基準スペクトルを用いた分解では、例えば基準スペクトルが互いに類似するような音種を分解したり、或いは基準スペクトルの設定が困難な音種を分解したりすることが難しい。よって、基準スペクトルを用いた分解だけでは、呼吸音のスペクトルを十分なまでに分解することができないおそれがある。
【００１８】
しかるに本発明では特に、上述したように、基準スペクトルを用いて分解された少なくとも一部のスペクトルは、所定の時系列特徴を用いた分解（即ち、第２分解手段による分解）によって更に分解される。これにより、基準スペクトルを用いた分解だけでは対応できない音種についても、好適に分解することが可能となる。具体的には、所定の時系列特徴を用いた分解では、例えば連続性ラ音とその他の音とを分解することができる。また、正常呼吸音と水泡音とを分解することもできる。
【００１９】
以上説明したように、本実施形態に係る呼吸音解析装置によれば、呼吸音に含まれるスペクトルを好適に分解して、呼吸音に含まれる各音種の割合に関する情報を出力することが可能である。
【００２０】
＜２＞
本実施形態に係る呼吸音解析装置の一態様では、前記第１分解手段は、前記呼吸音のスペクトルの第１所定特徴に対応する周波数に関する情報を取得する周波数情報取得手段と、前記複数の基準スペクトルを、前記周波数に関する情報に応じてシフトさせ、周波数シフト基準スペクトルを取得するシフト手段と、前記周波数シフト基準スペクトルに基づいて、前記呼吸音のスペクトルを分解する第１分解処理手段とを有する。
【００２１】
この態様によれば、第１分解手段による分解時には、先ず周波数情報取得手段によって、呼吸音のスペクトルの第１所定特徴に対応する周波数に関する情報が取得される。なお、ここでの「第１所定特徴」とは、呼吸音のスペクトルに含まれる音種に応じて特定の周波数に発生する特徴を意味しており、例えば周波数解析された信号に現れるピーク等である。また「周波数に関する情報」とは、周波数を直接的に示す情報に限定されず、その周波数を間接的に導き出すことができるような情報を含む趣旨である。周波数情報取得手段では、例えば呼吸音を示す信号に対して、高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）等による周波数解析が実行され、解析結果の極大値（即ち、ピーク）に対応する周波数に関する情報が取得される。
【００２２】
周波数に関する情報が取得されると、シフト手段により、呼吸音を分類する基準となる複数の基準スペクトルが周波数に関する情報に応じてシフトされ、周波数シフト基準スペクトルが取得される。基準スペクトルは、例えば呼吸音から取得された第１所定特徴であるピーク位置等に応じて周波数シフトされ、周波数シフト基準スペクトルとされる。
【００２３】
周波数シフト基準スペクトルが取得されると、第１分解処理手段により、周波数シフト基準スペクトルに基づくスペクトルの分解が実行される。即ち、呼吸音のスペクトルが、複数の周波数シフト基準スペクトルの各々に対応する音種に分解される。より具体的には、例えば呼吸音のスペクトルに対して、複数の周波数シフト基準スペクトルを基底とする演算が実行されることで、呼吸音のスペクトルに含まれる周波数シフト基準スペクトルの各々の割合が結合係数として算出される。なお、周波数シフト基準スペクトルの割合を算出する際には、非負の近似法（即ち、結合係数が負にならない近似法）が用いられてもよい。非負の近似法としては、例えば非負値行列因子分解（ＮＭＦ：Non-negative Matrix Factorization）が挙げられる。
【００２４】
以上説明したように、本態様に係る第１分解手段によれば、分解対象である呼吸音に応じてシフトされた周波数シフト基準スペクトルが用いられるため、より好適に呼吸音のスペクトルを分解することができる。
【００２５】
＜３＞
本実施形態に係る呼吸音解析装置の他の態様では、前記第２分解手段は、前記第１分解手段によって分解された少なくとも一部のスペクトルの第２所定特徴に対応する周波数を取得する周波数取得手段と、前記第２所定特徴に対応する周波数が時系列的に連続するか否かに応じて、前記第１分解手段によって分解された少なくとも一部のスペクトルを分解する第２分解処理手段とを有する。
【００２６】
この態様によれば、第２分解手段による分解時には、先ず周波数取得手段によって、第１分解手段により分解された少なくとも一部のスペクトルから第２所定特徴に対応する周波数が取得される。なお、ここでの「第２所定特徴」とは、既出の第１所定特徴と同様に、呼吸音成分に含まれる音種に応じて特定の周波数に発生する特徴を意味しており、例えば周波数解析された信号に現れるピーク等である。第２所定特徴は、第１所定特徴と同一の特徴であってもよいし、異なる特徴であってもよい。ちなみに、第２所定特徴に対応する周波数の取得は、後述する時間的な連続性を判定するために、連続して複数回行われる。このようにして取得された複数の周波数は、例えばバッファ等の記憶手段に一時的に記憶される。
【００２７】
周波数の取得後には、第２分解処理手段により、第１分解手段によって分解された少なくとも一部のスペクトルの分解が実行される。第２分解処理手段では、第２所定特徴に対応する周波数が時系列的に連続するか否かに応じてスペクトルが分解される。なお、ここでの「時系列的に連続する」とは、取得された周波数の時系列上での変化に連続性が認められる状態を指し、例えば時間的に連続して取得された２つの周波数が、所定の周波数範囲に入っているか否かによって判定できる。
【００２８】
第２所定特徴に対応する周波数が時系列的に連続するか否かを判定できれば、第２所定特徴が時間的な連続性を有する音種（例えば、連続性ラ音等）と、時間的な連続性を有さない音種（例えば、正常呼吸音や水泡音等）との分解を実現できる。よって、第１分解手段における基準スペクトルを用いた分解では分解できない、或いは分解することが難しい音種についても、第２分解手段において好適に分解することが可能である。
【００２９】
＜４＞
上述した第２分解処理手段を備える態様では、前記第２分解手段は、前記第２分解処理手段によって、前記第２所定特徴に対応する周波数が時系列的に連続するとして分解されたスペクトルを、前記第２所定特徴に対応する周波数と所定の閾値との関係に応じて更に分解する第３分解処理手段を有するように構成されてもよい。
【００３０】
この場合、第２分解処理手段によって分解されたスペクトルのうち、第２所定特徴に対応する周波数が時系列的に連続するとして分解されたスペクトル（言い換えれば、第２所定特徴に対応する周波数が時系列的に連続しないとして分解されたスペクトル以外のスペクトル）について、第３分解処理手段による分解が更に実行される。
【００３１】
第３分解処理手段では、第２所定特徴に対応する周波数と所定の閾値との関係に応じて分解が実行される。なお、ここでの「所定の閾値」とは、時系列的に連続するとして分解されたスペクトルを分解するための閾値であり、時系列的に連続するとして分解されたスペクトルに含まれる異なる２以上の音種を互いに判別できる値として設定されている。所定の閾値を用いれば、例えば時系列的に連続するとして分解されたスペクトルを、周波数が所定の閾値以上となるスペクトルと、周波数が所定の閾値未満となるスペクトルとに分解することができる。より具体的には、連続性ラ音に対応するスペクトルを、笛声音に対応するスペクトルと、類鼾音に対応するスペクトルとに分解することができる。
【００３２】
なお、所定の閾値は、複数設定されていてもよい。また、所定の閾値は、取得された周波数に応じて可変な値とされてもよい。
【００３３】
＜５＞
本実施形態に係る呼吸音解析装置の他の態様では、前記第２分解手段は、前記第１分解手段によって分解された少なくとも一部のスペクトルの振幅値を取得する振幅値取得手段と、前記振幅値が時系列的に連続するか否かに応じて、前記第１分解手段によって分解された少なくとも一部のスペクトルを分解する第４分解処理手段とを有する。
【００３４】
この態様によれば、第２分解手段による分解時には、先ず振幅値取得手段によって、第１分解手段により分解された少なくとも一部のスペクトルの振幅値が取得される。振幅値の取得は、後述する時間的な連続性を判定するために、連続して複数回行われる。このようにして取得された複数の振幅値は、例えばバッファ等の記憶手段に一時的に記憶される。
【００３５】
振幅値の取得後には、第４分解処理手段により、第１分解手段によって分解された少なくとも一部のスペクトルの分解が実行される。第４分解処理手段では、取得された振幅値が時系列的に連続するか否かに応じてスペクトルが分解される。なお、ここでの「時系列的に連続する」とは、取得された振幅値の時系列上での変化に連続性が認められる状態を指し、例えば時間的に連続して取得された２つの振幅値が、所定の振幅値範囲に入っているか否かによって判定できる。
【００３６】
振幅値が時系列的に連続するか否かを判定できれば、振幅値が時間的な連続性を有する音種（例えば、正常呼吸音等）と、時間的な連続性を有さない音種（例えば、水泡音等）との分解を実現できる。よって、第１分解手段における基準スペクトルを用いた分解では分解できない、或いは分解することが難しい音種についても、第２分解手段において好適に分解することが可能である。
【００３７】
なお、第４分解処理手段に加えて、上述した第２分解処理手段（即ち、周波数が時系列的に連続しているか否かに応じて分解を実行する手段）が備えられている場合には、第２分解処理手段と第４分解処理手段は、典型的には第１分解手段によって分解されたスペクトルのうち相異なる一部を夫々分解するように設定される。ただし、第２分解処理手段で分解されたスペクトルを第４分解処理手段で更に分解するようにしてもよいし、逆に第４分解処理手段で分解されたスペクトルを第２分解処理手段で更に分解するようにしてもよい。
【００３８】
＜６＞
本実施形態に係る呼吸音解析方法は、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程とを備える。
【００３９】
本実施形態に係る呼吸音解析方法によれば、上述した本実施形態に係る呼吸音解析装置と同様に、呼吸音に含まれるスペクトルを好適に分解して、呼吸音に含まれる各音種の割合に関する情報を出力することが可能である。
【００４０】
なお、本実施形態に係る呼吸音解析方法においても、上述した本実施形態に係る呼吸音解析装置における各種態様と同様の各種態様を採ることが可能である。
【００４１】
＜７＞
本実施形態に係るコンピュータプログラムは、呼吸音のスペクトルを、前記呼吸音を分類する基準となる複数の基準スペクトルに基づいて分解する第１分解工程と、前記第１分解工程において分解された少なくとも一部のスペクトルを、所定の時系列特徴に基づいて分解する第２分解工程と、前記第１分解工程及び前記第２分解工程において分解されたスペクトルに基づいて、前記呼吸音に含まれる前記分解されたスペクトルの各々の割合に関する情報を出力する出力工程とをコンピュータに実行させる。
【００４２】
本実施形態に係るコンピュータプログラムによれば、上述した本実施形態に係る呼吸音解析方法と同様の処理をコンピュータに実行させることができるため、呼吸音に含まれるスペクトルを好適に分解して、呼吸音に含まれる各音種の割合に関する情報を出力することが可能である。
【００４３】
なお、本実施形態に係るコンピュータプログラムにおいても、上述した本実施形態に係る呼吸音解析装置における各種態様と同様の各種態様を採ることが可能である。
【００４４】
＜８＞
本実施形態に係る記録媒体は、上述したコンピュータプログラムが記録されている。
【００４５】
本実施形態に係る記録媒体によれば、上述したコンピュータプログラムをコンピュータにより実行させることにより、呼吸音に含まれるスペクトルを好適に分解して、呼吸音に含まれる各音種の割合に関する情報を出力することが可能である。
【００４６】
本実施形態に係る呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の作用及び他の利得については、以下に示す実施例において、より詳細に説明する。
【実施例】
【００４７】
以下では、図面を参照して呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の実施例について詳細に説明する。【００４８】
＜全体構成＞
先ず、本実施例に係る呼吸音解析装置の全体構成について、図１を参照して説明する。ここに図１は、本実施例に係る呼吸音解析装置の全体構成を示すブロック図である。【００４９】
図１において、本実施例に係る呼吸音解析装置は、主な構成要素として、生体音取得部１１０と、第１分解部１２０と、第２分解部１３０と、成分量算出部１４０と、結果出力部１５０とを備えて構成されている。【００５０】
生体音取得部１１０は、生体の呼吸音を取得可能なセンサ等として構成される。生体音取得部１１０は、例えばＥＣＭ（Electret Condenser Microphone）やピエゾを利用したマイク、振動センサ等で構成されている。生体音取得部１１０で取得された呼吸音は、第１分解部１２０に出力される構成となっている。【００５１】
第１分解部１２０は、「第１分解手段」の一具体例であり、複数の演算回路やメモリ等を含んで構成されている。具体的には、第１分解部１２０は、周波数解析部１２１、ピーク周波数検出部１２２、基底集合生成部１２３、混合モデル基準データベース記憶部１２４、及び基底混合比算出部１２５を備えている。第１分解部１２０は、生体音取得部１１０で取得された呼吸音を、基底集合を利用して複数の音種に対応する成分へと分解する。第１分解部１２０での分解結果は、第２分解部１３０へと出力される構成となっている。第１分解部１２０に備えられる各部の動作については、後に詳述する。【００５２】
第２分解部１３０は、「第２分解手段」の一具体例であり、複数の演算回路やメモリ等を含んで構成されている。具体的には、第２分解部１３０は、ピーク周波数連続性判定部１３１、周波数記憶部１３２、結合部１３３、振幅値連続性判定部１３４、及び振幅値記憶部１３５を備えている。第２分解部１３０は、第１分解部１２０で分解された複数の音種に対応する成分を、時系列的特徴に基づいて更に分解する。第２分解部１３０での分解結果は、成分量算出部１４０へと出力される構成となっている。第２分解部１３０に備えられる各部の動作については、後に詳述する。【００５３】
成分量算出部１４０は、第１分解部１２０及び第２分解部１３０での分解結果に基づいて、生体音取得部１１０で取得された呼吸音に含まれる各音種の成分量を算出するものとして構成されている。成分量算出部１４０で算出された成分量を示す情報は、結果出力部１５０に出力される構成となっている。【００５４】
結果出力部１５０は、成分量算出部１４０で算出された各音種の成分量を、例えばディスプレイ等の画像や映像を表示可能な機器、或いはスピーカ等の音声を出力可能な機器に出力するものとして構成されている。【００５５】
＜動作説明＞
次に、本実施例に係る呼吸音解析装置の動作について、図１に加えて図２を参照して説明する。ここに図２は、本実施例に係る呼吸音解析装置の動作を示すフローチャートである。ここでは、本実施例に係る呼吸音解析装置が実行する処理の全体的な流れを把握するための簡単な説明を行う。各処理の詳細については、後述する。【００５６】
図１及び図２において、本実施例に係る呼吸音解析装置の動作時には、先ず生体音取得部１１０において呼吸音が取得される（ステップＳ１０１）。【００５７】
呼吸音信号が取得されると、周波数解析部１２１において周波数解析（例えば、高速フーリエ変換）が実行される（ステップＳ１０２）。また、ピーク周波数検出部１２２において、周波数解析結果を用いたピーク（極大値）の検出が実行される。【００５８】
続いて、基底集合生成部１２３において基底集合が生成される（ステップＳ１０３）。具体的には、基底集合生成部１２３は、混合モデル基準データベース記憶部１２４に記憶されている基底を用いて基底集合を生成する。この際、基底集合生成部１２３は、周波数解析結果から得られたピーク位置（即ち、対応する周波数）に基づいて、基底をシフトさせる。【００５９】
基底集合が生成されると、基底混合比算出部１２５において、周波数解析結果及び基底集合に基づく結合係数（即ち、各基底に対応する音種の成分量に対応する値）の算出が実行され（ステップＳ１０４）、結合係数に応じた信号強度が算出される（ステップＳ１０５）。また、基底混合比算出部１２５では更に、算出結果に基づいて音種パターン判定が行われ、分解された各成分が音種毎に異なる部位へと出力される（ステップＳ１０６）。【００６０】
具体的には、比較的なだらかな形状で表されるスペクトルを特徴とする成分は捻髪音成分と判定され、成分量算出部１４０に出力される（ステップＳ１０７）。即ち、捻髪音と判定された成分については、第２分解部１３０での分解は行われない。【００６１】
また、鋭いピークを有するスペクトルを特徴とする成分は主に連続性ラ音を含む成分として判定され、ピーク周波数連続性判定部１３１に出力される。ピーク周波数連続性判定部１３１では、ピーク周波数が時系列的に連続しているか否かが判定される（ステップＳ１０８）。ピーク周波数は、所定の時間間隔で取得された上で周波数記憶部１３２に記憶されており、ピーク周波数連続性判定部１３１は、周波数記憶部１３２に記憶されている過去のピーク周波数を利用して時系列的な連続性を判定する。【００６２】
ここで、ピーク周波数が時系列的に連続していると判定された成分は（ステップＳ１０８：ＹＥＳ）、連続性ラ音成分として更に周波数が所定閾値以上であるか否か判定される（ステップＳ１０９）。ピーク周波数が所定閾値以上である成分は（ステップＳ１０９：ＹＥＳ）、笛声音成分と判定され、成分量算出部１４０に出力される（ステップＳ１１０）。一方、ピーク周波数が所定閾値未満である成分は（ステップＳ１０９：ＮＯ）、類鼾音成分と判定され、成分量算出部１４０に出力される（ステップＳ１１１）。なお、ピーク周波数が時系列的に連続していないと判定された成分は（ステップＳ１０８：ＮＯ）、正常呼吸音及び水泡音を含む成分として、結合部１３３へと出力される。【００６３】
捻髪音成分及び連続性ラ音を含む成分以外の成分として判定された成分は、正常呼吸音及び水泡音を含む成分として、結合部１３３に出力される。この正常呼吸音及び水泡音を含む成分は、結合部１３３において、ピーク周波数連続性判定部１３１においてピーク周波数が時系列的に連続していないと判定された成分と結合される。結合された正常呼吸音及び水泡音を含む成分は、振幅値連続性判定部１３４に出力される。【００６４】
振幅値連続性判定部１３４では、スペクトルの振幅値が時系列的に連続しているか否かが判定される（ステップＳ１１２）。スペクトルの振幅値は、所定の時間間隔で取得された上で振幅値記憶部１３５に記憶されており、振幅値連続性判定部１３４は、振幅値記憶部１３５に記憶されている過去の振幅値を利用して時系列的な連続性を判定する。【００６５】
ここで、振幅値が時系列的に連続していると判定された成分は（ステップＳ１１２：ＹＥＳ）、正常呼吸音成分と判定され、成分量算出部１４０に出力される（ステップＳ１１３）。一方、振幅値が時系列的に連続していないと判定された成分は（ステップＳ１１２：ＮＯ）、水泡音と判定され、成分量算出部１４０に出力される（ステップＳ１１４）。【００６６】
続いて、成分量算出部１４０では、第１分解部１２０及び第２分解部１３０での分解結果に基づいて、信号強度が算出される（ステップＳ１１５）。信号強度が算出されると、結果出力部１５０において、信号強度（即ち、呼吸音に含まれる各音種の成分量）を示す画像データ等が生成され、外部のディスプレイ等において解析結果として表示される（ステップＳ１１６）。【００６７】
その後、解析処理を継続するか否かの判定が実行される（ステップＳ１１７）。解析処理を継続すると判定された場合（ステップＳ１１７：ＹＥＳ）、ステップＳ１０１からの処理が再び実行される。解析処理を継続しないと判定された場合（ステップＳ１１７：ＮＯ）、一連の処理は終了する。【００６８】
＜呼吸音信号の具体例＞
次に、本実施例に係る呼吸音解析装置で解析される呼吸音信号の具体例について、図３及び図４を参照して説明する。ここに図３は、捻髪音を含む呼吸音の周波数解析結果を示すスペクトログラム図であり、図４は、笛声音を含む呼吸音の周波数解析結果を示すスペクトログラム図である。【００６９】
図３に示す例では、正常呼吸音に対応するスペクトログラムパターンに加えて、異常呼吸音の１つである捻髪音に対応するスペクトログラムパターンが観測されている。捻髪音に対応するスペクトログラムパターンは、図中の拡大部分に示すように、菱形に近い形状である。【００７０】
図４に示す例では、正常呼吸音に対応するスペクトログラムパターンに加えて、異常呼吸音の１つである笛声音に対応するスペクトログラムパターンが観測されている。笛声音に対応するスペクトログラムパターンは、図中の拡大部分に示すように、白鳥の首のような形状である。【００７１】
このように、異常呼吸音には複数の音種が存在し、その音種によって異なる形状のスペクトログラムパターンとして観測される。ただし、図を見ても分かるように、正常呼吸音及び異常呼吸音は互いに混じり合った状態で検出される。本実施例に係る呼吸音解析装置は、このように混じり合った複数の音種を分解するための処理を実行する。【００７２】
＜呼吸音信号の近似方法＞
次に、本実施例に係る呼吸音解析装置の第１分解部１２０による分解処理（即ち、基底集合を利用した分解）の概略について、図５から図８を参照して簡単に説明する。ここに図５は、捻髪音を含む呼吸音の所定タイミングにおけるスペクトルを示すグラフであり、図６は、捻髪音を含む呼吸音のスペクトルの近似方法を示す概念図である。また図７は、笛声音を含む呼吸音の所定タイミングにおけるスペクトルを示すグラフであり、図８は、笛声音を含む呼吸音のスペクトルの近似方法を示す概念図である。【００７３】
図５において、捻髪音を含む呼吸音信号（図３参照）について、捻髪音に対応するスペクトログラムパターンが強く現れているタイミングでスペクトルを抽出すると、図に示すような結果が得られる。このスペクトルは、正常呼吸音と捻髪音とを含んでいると考えられる。【００７４】
図６において、正常呼吸音に対応するスペクトル及び捻髪音に対応するスペクトルは、予め実験等により推定できる。このため、予め推定したパターンを利用すれば、上述したスペクトルについて、正常呼吸音に対応する成分と捻髪音に対応する成分とがどのような割合で含まれているかを知ることができる。【００７５】
図７において、笛声音を含む呼吸音信号（図４参照）について、笛声音に対応するスペクトログラムパターンが強く現れているタイミングでスペクトルを抽出すると、図に示すような結果が得られる。このスペクトルは、正常呼吸音と笛声音とを含んでいると考えられる。【００７６】
図８において、上述した正常呼吸音及び捻髪音と同様に、笛声音に対応するスペクトルについても、予め実験等により推定できる。このため、予め推定したパターンを利用すれば、上述したスペクトルについて、正常呼吸音に対応する成分と笛声音に対応する成分とがどのような割合で含まれているかを知ることができる。【００７７】
以下では、このような解析を実現するための各処理について、より具体的に説明する。【００７８】
＜周波数解析＞
呼吸音信号の周波数解析及び解析結果におけるピークの検出について、図９から図１１を参照して詳細に説明する。ここに図９は、周波数解析方法の一例を示すグラフであり、図１０は、周波数解析結果の一例を示す図である。また図１１は、スペクトルのピーク検出結果を示す概念図である。【００７９】
図９において、取得された呼吸音信号に対しては、先ず周波数解析が実行される。周波数は、高速フーリエ変換等の既存の技術を利用して行うことができる。本実施例では、周波数毎の振幅値（即ち、振幅スペクトル）を周波数解析結果として用いている。なお、データ取得時のサンプリング周波数、窓サイズ、窓関数（例えば、ハニング窓等）については、適宜決定すればよい。【００８０】
図１０に示すように、周波数解析結果は、ｎ個の値で構成されるものとして得られる。なお、「ｎ」は、周波数解析における窓サイズ等によって決まる値である。【００８１】
図１１において、周波数解析によって得られたスペクトルについては、ピークの検出が実行される。図に示す例では、１００Ｈｚ、１３０Ｈｚ、１８０Ｈｚ，及び３２０Ｈｚの位置にピークｐ１〜ｐ４が検出されている。なお、ピークの検出処理については、どの周波数にピークが存在するかだけわかればよいため、簡易的な処理でも構わない。ただし、小さなピークでも取りこぼしがないよう、ピーク検出のパラメータ設定されていることが好ましい。【００８２】
本実施例では、極大値を取る点を求め、更にその点の２階微分値の小さいもの（即ち、絶対値が大きいもの）から順に最大Ｎ個（Ｎは所定の値）を検出している。極大値は、差分の符号が正から負に切り替わる点から求められる。２階微分値は差分の差分で近似する。この値が所定の閾値（負の値）より小さいものを、小さいものから順に最大Ｎ個選び、その位置を記憶する。【００８３】
＜基底集合の生成＞
次に、基底集合の生成について、図１２から図１６を参照して詳細に説明する。ここに図１２は、正常肺胞呼吸音基底を示すグラフである。また図１３は、捻髪音基底を示すグラフであり、図１４は、連続性ラ音基底を示すグラフであり、図１５は、ホワイトノイズ基底を示すグラフである。図１６は、周波数シフトされた連続性ラ音基底を示すグラフである。【００８４】
図１２から図１５に示すように、各音種に対応する基底（「基準スペクトル」の一具体例）は、特有の形状を有している。なお、各基底は周波数解析結果と同じｎ個の数値（即ち、周波数ごとの振幅値）で構成されている。なお、各基底は、周波数毎の振幅値を示す線と周波数軸とで囲まれた面積が所定の値（例えば１）になるように正規化されている。【００８５】
ちなみに、ここでは正常肺胞呼吸音基底、捻髪音基底、連続性ラ音基底、ホワイトノイズ基底の４つの基底を示しているが、１つの基底しかない場合でも分解処理を実行することができる。また、ここで挙げた基底以外の基底を用いることもできる。なお、ここで挙げた呼吸音に対応する基底に代えて、例えば心拍音や腸音に対応する基底を用いれば、心拍音や腸音の分解を実行することも可能である。【００８６】
図１６において、上述した基底のうち連続性ラ音に対応する基底は、周波数解析の結果から検出されたピーク位置に合わせて周波数シフトされる。ここでは、図１１で示したピークｐ１〜ｐ４の各々に合わせて、連続性ラ音基底を周波数シフトさせた例を示している。なお、連続性ラ音に対応する基底以外の基底を周波数シフトさせてもよい。【００８７】
以上の結果、基底集合は、正常肺胞呼吸音基底、捻髪音基底、ピーク検出個数分の連続性ラ音基底、及びホワイトノイズ基底の集合として生成される。【００８８】
＜結合係数の算出＞
次に、結合係数の算出について、図１７から図１９を参照して詳細に説明する。ここに図１７は、スペクトルと、基底及び結合係数との関係を示す図であり、図１８は、観測されたスペクトル及び近似に用いられる基底の一例を示す図である。また図１９は、非負値行列因子分解による近似結果を示す図である。【００８９】
解析対象であるスペクトルｙ、基底ｈ（ｆ）、及び結合係数ｕの関係は、以下の数式（１）で表すことができる。【００９０】
【数１】
2016002004.tif 000003
図１７に示すように、スペクトルｙ及び各基底ｈ（ｆ）は、ｎ個の値を有している。他方、結合係数は、ｍ個の値を有している。なお、「ｍ」は、基底集合に含まれる基底の数である。【００９１】
本実施例に係る呼吸音解析装置では、非負値行列因子分解を利用して基底集合に含まれる各基底の結合係数を算出する。具体的には、以下の数式（２）で示される最適化基準関数Ｄを最小化するｕ（ただし、ｕの各成分値は非負）を求めればよい。【００９２】
【数２】
2016002004.tif 000004
なお、一般的な非負値行列因子分解は、基底スペクトルの集合を表す基底行列と、結合係数を表すアクティベーション行列を共に算出する手法であるが、本実施例においては、基底行列を固定して結合係数のみを算出している。【００９３】
ちなみに、結合係数を算出するための手段として、非負値行列因子分解以外の近似法を用いてもよい。ただし、この場合においても非負であるという条件が望まれる。以下では、非負の近似法を用いる理由について、具体例を挙げて説明する
図１８に示すように、観測されたスペクトルを、基底Ａ〜Ｄの４つの基底で近似して結合係数を算出する場合を考える。なお、非負であることを条件とした場合の期待する結合係数ｕは、基底Ａに対応するものが１、基底Ｂに対応するものが１、基底Ｃに対応するものが０、基底Ｄに対応するもの０である。即ち、非負であることを条件とした場合、観測されたスペクトルは、基底Ａに１を乗じたものと、基底Ｂに１を乗じたものとを足し合わせたスペクトルとして近似される。【００９４】
一方、非負であることを条件としない場合に算出され得る結合係数ｕは、基底Ａに対応するものが０、基底Ｂに対応するものが０、基底Ｃに対応するものが１、基底Ｄに対応するものが−０．５である。即ち、非負であることを条件としない場合、観測されたスペクトルは、基底Ｃに１を乗じたものと、基底Ｄに−０．５を乗じたものとを足し合わせたスペクトルとして近似される。【００９５】
上述した２つの例を比較した場合、非負であることを条件とする場合よりも、非負であることを条件としない場合の方が高い近似精度を得られることがある。しかしながら、ここでの結合係数ｕはスペクトルごとの成分量を表すものであるため、非負の値として得られなければならない。言い換えれば、結合係数ｕが負の値で得られた場合には、成分量としての解釈ができない。これに対し、非負の条件を課して近似を行えば、成分量に対応する結合係数ｕを算出することができる。【００９６】
図１９において、本実施例に係る生体解析装置では、上述したように、正常肺胞呼吸音基底、捻髪音基底、４つの連続性ラ音基底、及びホワイトノイズ基底からなる基底集合を用いて結合係数ｕを算出するため、結合係数ｕは、ｕ１からｕ７の７個の値を有するものとして算出される。【００９７】
ここで、正常肺胞呼吸音基底に対応する結合係数ｕ１は、呼吸音に対する正常肺胞呼吸音の割合を示す値であると言える。同様に、捻髪音基底に対応する結合係数ｕ２、ホワイトノイズ基底に対応する結合係数ｕ３、１００Ｈｚにシフトした連続性ラ音基底に対応する結合係数ｕ４、１３０Ｈｚにシフトした連続性ラ音基底に対応する結合係数ｕ５、１８０Ｈｚにシフトした連続性ラ音基底に対応する結合係数ｕ６、及び３２０Ｈｚにシフトした連続性ラ音基底に対応する結合係数ｕ７の各々についても、呼吸音に対する各音種の割合を示す値であると言える。従って、結合係数ｕから各音種の信号強度を算出することができる。【００９８】
以上のように、本実施例では、各音種に対応する複数の基底を利用して呼吸音に含まれる複数の音種を分解する。【００９９】
＜ピーク周波数の時間的連続性を利用した分解＞
次に、ピーク周波数連続性判定部１３１において行われる分解処理について、図２０及び図２１を参照して具体的に説明する。ここに図２０及び図２１は夫々、ピーク周波数の時間的連続性を利用した分解方法を示す概念図である。【０１００】
ピーク周波数連続性判定部１３１は、上述したように、第１分解部１２０において連続性ラ音を含む成分として分解された一つ以上の成分を、連続性ラ音とそれ以外の音（例えば、正常呼吸音及び水泡音等）に分解する。具体的には、呼吸音信号の周波数解析結果から検出されるピーク周波数が、所定の範囲内で変動している場合に連続性ラ音であると判定する。【０１０１】
図２０に示すように、連続性ラ音である笛声音や類鼾音は、時間軸上で連続して検出されるピークの位置が所定の範囲内に収まるように変動する。言い換えれば、ピーク周波数が時間的に連続性を有するように変化する。よって、連続するピーク位置が所定の範囲内にある場合には、その音が連続性ラ音であると判別できる。【０１０２】
他方、図２１に示すように、連続性ラ音以外の音は、時間軸上で連続して検出されるピークの位置が所定の範囲内に収まらないように変動する。言い換えれば、ピーク周波数が時間的な連続性を有さず離散的に変化する。よって、連続するピーク位置が所定の範囲内でない場合には、その音が連続性ラ音でないと判別できる。【０１０３】
なお、連続性ラ音の判定には、複数回の判定結果を用いることもできる。具体的には、時間軸上で連続して検出されるピークの位置が所定の範囲内に収まるように変動している回数が所定回数以上継続した場合に、その音が連続性ラ音であると判定するようにしてもよい。【０１０４】
＜ピーク周波数に対する閾値を利用した分解＞
次に、ピーク周波数の時間的連続性を利用した分解の後に実行される、ピーク周波数に対する閾値を利用した分解について説明する。なお、以下では、３つの異なる実施例について説明する。【０１０５】
＜第１実施例＞
先ず、第１実施例に係る閾値を利用した分解方法について、図２２を参照して説明する。ここに図２２は、第１実施例に係る笛声音と類鼾音との分解に用いる閾値を示すグラフである。【０１０６】
第１実施例に係る分解方法では、ピーク周波数の時間的連続性を利用した分解の結果として得られる一つ以上の連続性ラ音成分を、所定閾値との比較により笛声音と類鼾音とに分解する。ここで、笛声音は高音性連続性ラ音、類鼾音は低音性連続性ラ音と呼ばれるように、笛声音と類鼾音とは音の高さ（即ち、周波数）で判別することが可能である。しかしながら、笛声音及び類鼾音は、ピーク周波数が時間的に変化する。このため、ピーク周波数に対する単一の閾値（即ち、値が変動しない一つの閾値）を利用して笛声音及び類鼾音を判定しようとすると、時間の経過により、判定結果が変化してしまうことがある。例えば、ピーク周波数が判定閾値を跨ぐように変化してしまうと、それまでは正確に判定されていたものが、誤った音種として判定されることになってしまう。このため第１実施例では、ピーク周波数に応じて判定閾値を変動させる。【０１０７】
図２２に示すように、第１実施例に係る分解処理では、笛声音と判定する割合及び類鼾音と判定する割合がピーク周波数に応じてなめらかに変化するように閾値が変動する。例えば、ピーク周波数が２００Ｈｚの場合には、笛声音が７％含まれ、類鼾音が９３％含まれると判定する。ピーク周波数が２５０Ｈｚの場合には、笛声音が５０％含まれ、類鼾音が５０％含まれると判定する。ピーク周波数が２８０Ｈｚの場合には、笛声音が７８％含まれ、類鼾音が２２％含まれると判定する。なお、ここでの具体的な数値はあくまで一例であり、異なる値を設定してもよい。また、測定対象である生体の性別、年齢、身長、体重等によって異なる変動特性を有するようにしてもよい。【０１０８】
上述した変動する閾値を利用することで、ピーク周波数の変動に起因する誤判定を好適に防止することができる。即ち、第１実施例に係る分解処理では、笛声音及び類鼾音を判定するための閾値がピーク周波数に応じて適切な値にとなるよう変動するため、例えば変動しない単一の閾値を用いる場合と比較して、より正確な分解が行える。【０１０９】
＜第２実施例＞
次に、第２実施例に係る分解処理について、図２３から図２５を参照して説明する。ここに図２３は、第２実施例に係る笛声音と類鼾音との分解に用いる閾値の初期値を示すグラフである。また図２４及び図２５は夫々、第２実施例に係る笛声音と類鼾音との分解に用いる閾値の調整後の値を示すグラフである。【０１１０】
図２３に示すように、第２実施例に係る分解処理では、閾値である２５０Ｈｚを境にして判定結果が変化するものとして設定されている。具体的には、ピーク周波数が２５０Ｈｚ以上である場合には、連続性ラ音は笛声音成分を１００％含んでおり、類鼾音は含んでいないと判定される。一方、ピーク周波数が２５０Ｈｚ未満である場合には、連続性ラ音は類鼾音成分を１００％含んでおり、笛声音は含んでいないと判定される。【０１１１】
図２４に示すように、第２実施例に係る分解処理では、直前の判定において笛声音成分を１００％含むものであると判定された場合、閾値が２５０Ｈｚから２２０Ｈｚへと低くされる。よって、笛声音成分を１００％含むものとして判定され易くなる。具体的には、ピーク周波数が２３０Ｈｚの場合を考えると、初期の閾値（図２３参照）によれば類鼾音と判定されることになるが、調整後の閾値（図２４参照）によれば笛声音と判定される。【０１１２】
図２５に示すように、第２実施例に係る分解処理では、直前の判定において類鼾音成分を１００％含むものであると判定された場合、閾値が２５０Ｈｚから２８０Ｈｚへと高くされる。よって、類鼾音成分を１００％含むものとして判定され易くなる。具体的には、ピーク周波数が２７０Ｈｚの場合を考えると、初期の閾値（図２３参照）によれば笛声音と判定されることになるが、調整後の閾値（図２５参照）によれば類鼾音と判定される。【０１１３】
上述したように閾値を調整すれば、ピーク周波数の変動に起因する誤判定を好適に防止することができる。即ち、第２実施例に係る分解処理では、笛声音及び類鼾音を判定するための閾値が過去の判定結果に基づいて適切なものへと調整されるため、例えば調整されない単一の閾値を用いる場合と比較して、より正確な判定が行える。【０１１４】
なお、閾値の調整は、直前の判定結果だけによらず、複数回の過去の判定結果に基づいて行われてもよい。また、複数回の過去の判定結果を利用する場合は、各判定結果に対して重み付けを行ってもよい。例えば、過去の判定結果であるほど影響が小さくなるように重み付けを行ってもよい。また、調整する閾値の初期値として、第１実施例のなめらかな閾値を用いてもよい（図２２参照）。【０１１５】
＜第３実施例＞
次に、第３実施例に係る分解処理について、図２６から図２９を参照して説明する。ここに図２６は、笛声音を含む呼吸音のスペクトログラム図であり、図２７は、笛声音のピーク周波数及びピーク数を示すグラフである。また図２８は、類鼾音を含む呼吸音のスペクトログラム図であり、図２９は、類鼾音のピーク周波数及びピーク数を示すグラフである。【０１１６】
図２６において、笛声音を含む呼吸音は、所定のピークを有するスペクトル波形として検出される。ここからピーク周波数Ｆ及びピーク数Ｎを検出するには、先ずスペクトル波形の単一時間（即ち、図中の白枠で囲った領域）に対応する周波数−振幅グラフを作成する。【０１１７】
図２７に示すグラフから、笛声音のピーク周波数Ｆ１及びピーク数Ｎ１が検出できる。なお、笛声音のピーク周波数の分布は、１８０〜９００Ｈｚ程度であることが分かっている。また、図を見ても分かるように、笛声音のピーク数Ｎ１は１個である。【０１１８】
図２８において、類鼾音を含む呼吸音は、笛声音とは異なる所定のピークを有するスペクトル波形として検出される。ここからピーク周波数Ｆ及びピーク数Ｎを検出するには、同様にスペクトル波形の単一時間に対応する周波数−振幅グラフを作成する。【０１１９】
図２９に示すグラフから、類鼾音のピーク周波数Ｆ２及びピーク数Ｎ２が検出できる。なお、類鼾音のピーク周波数の分布は、１００〜２６０Ｈｚ程度であることが分かっている。即ち、類鼾音のピーク周波数Ｆ２は、笛声音のピーク周波数Ｆ１よりも低い領域に分布していることになる。また、図を見ても分かるように、類鼾音のピーク数Ｎ２は例えば、３個である。即ち、類鼾音のピーク数Ｎ２は、笛声音のピーク数Ｎ１のように１つでなく、複数である。【０１２０】
第３実施例に係る分解処理では、上述した笛声音及び類鼾音の特性の違いを利用して判定が行われる。具体的には、ピーク周波数Ｆ及びピーク数Ｎの各々に基づいて、笛声音と類鼾音とが分解される。このようにすれば、例えばピーク周波数Ｆだけを利用して笛声音と類鼾音とを分解する場合と比べて、より正確な分解が行える。【０１２１】
＜振幅値の時間的連続性を利用した分解＞
次に、振幅値連続性判定部１３４において行われる分解処理について、図３０及び図３１を参照して具体的に説明する。ここに図３０及び図３１は夫々、振幅値の時間的連続性を利用した分解方法を示す概念図である。【０１２２】
振幅値連続性判定部１３４は、上述したように、第１分解部１２０及びピーク周波数連続性判定部１３１において正常呼吸音及び水泡音を含む成分として分解された一つ以上の成分を、正常呼吸音と水泡音とに分解する。具体的には、呼吸音信号の周波数解析結果から検出されるスペクトルの振幅値が、所定の範囲内で変動している場合に正常呼吸音であると判定し、所定の範囲内で変動していない場合に水泡音であると判定する。【０１２３】
図３０に示すように、正常呼吸音は、時間軸上で連続して検出される振幅値が所定の範囲内に収まるように変動する。言い換えれば、振幅値が時間的に連続性を有するように変化する。よって、連続する振幅値が所定の範囲内にある場合には、その音が正常呼吸音であると判別できる。【０１２４】
他方、図３１に示すように、水泡音は、時間軸上で連続して検出される振幅値が所定の範囲内に収まらないように変動する。言い換えれば、振幅値が時間的な連続性を有さず離散的に変化する。よって、連続する振幅値が所定の範囲内でない場合には、その音が水泡音であると判別できる。【０１２５】
なお、正常呼吸音と水泡音との分解処理には、複数回の判定結果を用いることもできる。具体的には、時間軸上で連続して検出される振幅値が所定の範囲内に収まるように変動している回数が所定回数以上継続した場合に、その音が正常呼吸音であると判定し、それ以外の場合には水泡音であると判定するようにしてもよい。【０１２６】
＜分解処理の具体例＞
次に、本実施例に係る呼吸音解析装置における分解処理について、図３２及び図３３を参照して、より具体的に説明する。ここに図３２及び図３３は夫々、分解処理の各段階において得られるスペクトログラムを、順を追って示す概念図である。【０１２７】
図３２に示す例では、捻髪音を多く含み、笛声音及び類鼾音は殆ど含まれない呼吸音を解析対象としている。【０１２８】
取得された呼吸音波形をＳＴＦＴにより解析することで、解析処理間隔毎のスペクトルが得られる。周波数解析後スペクトルは、基底集合を利用した分解によって、正常・水泡音代表成分（即ち、正常呼吸音及び水泡音の殆どを含む成分）、ピーク型スペクトル成分（即ち、主として連続性ラ音を含む成分）、捻髪音成分、ノイズ成分（図示省略）に分解される。【０１２９】
ピーク型スペクトラム成分は、ピーク周波数の連続性を利用した分解により、連続性ラ音成分と、正常・水泡音成分の一部とに分解される。連続性ラ音成分は、周波数閾値を利用した分解により、類鼾音成分と笛声音成分とに分解される。他方、正常・水泡音成分の一部は、基底集合を利用して分解された正常・水泡音代表成分と結合された後、振幅値の連続性を利用した分解によって、正常呼吸音成分と水泡音成分とに分解される。【０１３０】
ここで、分解後のスペクトログラムを見ても分かるように、捻髪音成分は多く抽出される一方で、笛声音成分及び類鼾音成分は殆ど抽出されていない。この結果から、呼吸音に含まれる各音種を正確に分解できていることが分かる。【０１３１】
一方、図３３に示す例では、笛声音及び類鼾音を多く含み、捻髪音は殆ど含まれない呼吸音を解析対象としている。【０１３２】
この場合の分解後のスペクトログラムを見てみると、笛声音成分及び類鼾音成分は多く抽出される一方で、捻髪音成分は殆ど抽出されていない。この結果からも、図３２の例と同様に、呼吸音に含まれる各音種を正確に分解できていることが分かる。【０１３３】
＜解析結果の表示例＞
次に、解析結果の具体的な表示例について、図３４及び図３５を参照して詳細に説明する。ここに図３４及び図３５は夫々、表示部における表示例を示す平面図である。なお、図３４の表示例は、図３２に示す解析結果（即ち、捻髪音を多く含む呼吸音の解析結果）に対応しており、図３５の表示例は、図３３に示す解析結果（即ち、笛声音及び類鼾音を多く含む呼吸音の解析結果）に対応している。【０１３４】
図３４及び図３５に示すように、本実施形態に係る呼吸音解析装置に接続される表示部の表示領域２００には、解析結果が複数の画像として表示される。具体的には、領域２００ａには、取得された呼吸音の波形が表示されている。領域２００ｂには、取得された呼吸音のスペクトルが表示されている。領域２００ｃには、取得された呼吸音のスペクトログラムが表示されている。領域２００ｄには、分解された各音種（ここでは、正常呼吸音、類鼾音、笛声音、捻髪音、水泡音の５音種）の成分量の時系列変化を表すグラフが表示されている。領域２００ｅには、分解された各音種の割合がレーダーチャートとして表示されている。【０１３５】
このような表示を行うことにより、解析結果を視覚的に分かりやすく提供することが可能である。即ち、ユーザに対して、呼吸音に含まれている各成分の割合を直感的に知らせることができる。なお、このような解析結果の表示態様はあくまで一例であり、他の表示態様で解析結果が表示されてもよい。例えば、分解された各音種の割合は、棒グラフや円グラフとして表示されてもよいし、数値化して表示されてもよい。【０１３６】
以上説明したように、本実施例に係る呼吸音解析装置によれば、呼吸音に含まれるスペクトルを好適に分解して、呼吸音に含まれる各音種の割合に関する情報を出力することが可能である。【０１３７】
本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。
【符号の説明】
【０１３８】
１１０  生体音取得部
１２０  第１分解部
１２１  周波数解析部
１２２  ピーク周波数検出部
１２３  基底集合生成部
１２４  混合モデル基準データベース記憶部
１２５  基底混合比算出部
１３０  第２分解部
１３１  ピーク周波数連続性判定部
１３２  周波数記憶部
１３３  結合部
１３４  振幅値連続性判定部
１３５  振幅値記憶部
１４０  成分量算出部
１５０  結果出力部
２００  表示領域
ｙ  スペクトル
ｈ（ｆ）  基底
ｕ  結合係数
【図１】
2016002004.tif 000005
【図２】
2016002004.tif 000006
【図３】
2016002004.tif 000007
【図４】
2016002004.tif 000008
【図５】
2016002004.tif 000009
【図６】
2016002004.tif 000010
【図７】
2016002004.tif 000011
【図８】
2016002004.tif 000012
【図９】
2016002004.tif 000013
【図１０】
2016002004.tif 000014
【図１１】
2016002004.tif 000015
【図１２】
2016002004.tif 000016
【図１３】
2016002004.tif 000017
【図１４】
2016002004.tif 000018
【図１５】
2016002004.tif 000019
【図１６】
2016002004.tif 000020
【図１７】
2016002004.tif 000021
【図１８】
2016002004.tif 000022
【図１９】
2016002004.tif 000023
【図２０】
2016002004.tif 000024
【図２１】
2016002004.tif 000025
【図２２】
2016002004.tif 000026
【図２３】
2016002004.tif 000027
【図２４】
2016002004.tif 000028
【図２５】
2016002004.tif 000029
【図２６】
2016002004.tif 000030
【図２７】
2016002004.tif 000031
【図２８】
2016002004.tif 000032
【図２９】
2016002004.tif 000033
【図３０】
2016002004.tif 000034
【図３１】
2016002004.tif 000035
【図３２】
2016002004.tif 000036
【図３３】
2016002004.tif 000037
【図３４】
2016002004.tif 000038
【図３５】
2016002004.tif 000039
【国際調査報告】
2016002004.tif 000040
2016002004.tif 000041
2016002004.tif 000042
2016002004.tif 000043