(57)【要約】【課題】被測定者から取得したシャント音を解析して、シャント狭窄診断を好適に支援する。【解決手段】シャント音解析装置は、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得手段（１１０）と、取得部が取得したシャント音情報に基づいて、シャント音の聴感上の断続度合いを示す情報を出力する出力手段（１６０）とを備える。出力される断続度合いを示す情報は、シャント形成部位における狭窄発生可能性を定量的に示すものである。従って、シャント音解析装置によれば、狭窄診断を好適に支援することが可能である。【選択図】図１
2020199293.tif 000002
【特許請求の範囲】
【請求項１】
被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得手段と、
前記取得部が取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力手段と
を備えることを特徴とするシャント音解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被測定者から取得したシャント音を解析するシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体の技術分野に関する。
【背景技術】
【０００２】
この種の装置として、被測定者から取得したシャント音を解析して、シャント狭窄等に関する医師の診断を支援する装置が知られている。例えば特許文献１では、アレイ状採音センサを被測定者の腕に固定して複数箇所からシャント音を取得し、取得したサンプルデータと、予め用意した多種多様なシャント狭窄音を含む指標サンプルとについて、ＳＴＭＥＭ法等による解析を行うシャント音解析装置が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−８２６３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、上述した特許文献１に記載されているようなシャント音解析装置では、解析を行うに際し、予め膨大な量の指標サンプルを用意することが要求されてしまうだけでなく、被測定者固有の情報等が考慮されないために、正確な解析結果が得られないおそれがあるという技術的問題点が生ずる。また、解析処理も複雑であり、処理負荷が極めて大きくなってしまうという技術的問題点も生ずる。
【０００５】
本発明が解決しようとする課題には上記のようなものが一例として挙げられる。本発明は、被測定者から取得したシャント音を解析して、シャント狭窄診断を好適に支援することが可能なシャント音解析装置を提供することを課題とする。
【課題を解決するための手段】
【０００６】
上記課題を解決するためのシャント音解析装置は、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得手段と、前記取得部が取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力手段とを備える。
【０００７】
上記課題を解決するためのシャント音解析方法は、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得工程と、前記取得工程で取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力工程とを備える。
【０００８】
上記課題を解決するためのコンピュータプログラムは、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得工程と、前記取得工程で取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力工程とをコンピュータに実行させる。
【０００９】
上記課題を解決するための記録媒体は、上述したコンピュータプログラムが記録されている。
【図面の簡単な説明】
【００１０】
【図１】実施例に係るシャント音解析装置の全体構成を示すブロック図である。
【図２】時間周波数波形の一例を示すスペクトログラム（その１）である。
【図３】音量解析波形の一例を示すグラフである。
【図４】時間周波数波形の一例を示すスペクトログラム（その２）である。
【図５】周波数重心解析波形の一例を示すグラフである。
【図６】個人差入力部の動作を示すフローチャートである。
【図７】被測定者の過去から蓄積されたシャント音の最小値の分布を示すヒストグラムである。
【図８】被測定者の過去から蓄積されたシャント音の最大値の分布を示すヒストグラムである。
【図９】正常時における１心拍区間の時間周波数波形の一例を示すスペクトログラムである。
【図１０】正常時における１心拍区間の周波数重心解析波形の一例を示すグラフである。
【図１１】狭窄発生時における１心拍区間の時間周波数波形の一例を示すスペクトログラムである。
【図１２】狭窄発生時における１心拍区間の周波数重心解析波形の一例を示すグラフである。
【図１３】正常時における１心拍区間の音量変化量の一例を示すグラフである。
【図１４】狭窄発生時における１心拍区間の音量変化量の一例を示すグラフである。
【図１５】正常時における１心拍区間の音量最小値の一例を示すグラフである。
【図１６】狭窄発生時における１心拍区間の音量最小値の一例を示すグラフである。
【図１７】正常時における１心拍区間の音量減衰率の一例を示すグラフである。
【図１８】狭窄発生時における１心拍区間の音量減衰率の一例を示すグラフである。
【図１９】歪み発生時の音量解析波形の一例を示すグラフである。
【図２０】歪み発生時の音量解析微分波形の一例を示すグラフである。
【発明を実施するための形態】
【００１１】
＜１＞
本実施形態に係るシャント音解析装置は、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得手段と、前記取得部が取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力手段とを備える。
【００１２】
本実施形態に係るシャント音解析装置の動作時には、先ず取得手段により、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報が取得される。なお、ここでの「シャント音」とは、血液を体外に取り出すためのシャント形成部位周辺において取得される血流音であり、被測定者の脈拍に同期した音である。シャント音の取得は各種センサを用いて行えばよく、その取得方法が特に限定されるものではない。また「シャント音情報」とは、シャント音に関する各種パラメータを含む情報であって、例えば音量や周波数等の時間変化などを含んでいる。
【００１３】
シャント音情報が取得されると、出力手段において各種解析が実行され、シャント音の聴感上の断続度合いを示す情報が出力される。ここでの「断続度合いを示す情報」とは、単にシャント音が途切れる程度を示す情報を意味するのではなく、狭窄を診断する際の指標とされる断続感が生じる原因となる各種パラメータを含む複合的な情報、或いは断続感に関連する複数のパラメータが統合された情報を意味している。なお、断続度合いを示す情報は、例えば数値化した状態、或いはグラフやチャートによって可視化された状態で出力される。
【００１４】
本願発明者の研究するところによれば、シャント形成部位に狭窄が生じていない場合には、シャント音は大きく連続した低い音となる傾向があるのに対し、狭窄が進行すると、シャント音が小さくなり、音の途切れや狭窄部位の乱流による高い音の成分が単独で又は複合的に表れてくることが判明している。よって、シャント音の断続度合いを示す情報を利用すれば、シャント形成部位に狭窄が生じているか否か、或いは狭窄の程度を容易且つ的確に診断することが可能となる。具体的には、従来から狭窄診断に利用されているエコー診断等を利用せずとも、容易に狭窄診断を行うことができる。また、シャント音の断続度合いを示す情報により、狭窄を診断する医師の技量や経験に左右されない定量的な診断が可能となる。
【００１５】
以上説明したように、本実施形態に係るシャント音解析装置によれば、シャント形成部位の狭窄診断を好適に支援することが可能である。
【００１６】
＜２＞
本実施形態に係るシャント音解析装置の一態様では、前記出力手段は、前記シャント音情報に基づいて、時間経過に対する前記シャント音の音量の変化を示す音量情報を導出する音量情報導出手段を有する。
【００１７】
この態様によれば、音量情報導出手段によって導出された音量情報に基づいて、シャント音の断続度合いを示す情報を出力することが可能となる。ここで特に、本願発明者の研究するところによれば、シャント音の断続感は、時間経過に対するシャント音の音量の変化に大きく依存していることが判明している。よって、導出された音量情報を利用すれば、シャント音の断続度合いを示す情報として、より適切な情報を出力することが可能となる。
【００１８】
＜３＞
上述した音量情報導出手段を備える態様では、前記出力手段は、前記音量情報から前記シャント音の一周期に相当する一周期音量情報を抽出する抽出手段を有してもよい。
【００１９】
この場合、音量情報が導出されると、そこからシャント音の一周期（言い換えれば、被測定者の脈動の一周期）に相当する一周期音量情報が抽出される。このため、シャント音の一周期における音量変化（例えば、収縮期から拡張期にかけての音量変化）を好適に解析できる。従って、シャント音の断続度合いを示す情報として、より適切な情報を出力することが可能となる。
【００２０】
＜４＞
上述した抽出手段を備える態様では、前記出力手段は、前記一周期音量情報における前記シャント音の音量の最大値及び最小値の差分と、正常時におけるシャント音の最大値及び最小値の差分とを対比して、前記断続度合いを数値化する第１数値化手段を有してもよい。
【００２１】
この場合、先ず一周期音量情報におけるシャント音の音量の最大値（例えば、収縮期の音量）から最小値（例えば、拡張期の音量）が減算され、一周期における音量変化量がされる。そして、算出された音量変化量が、予め記憶された正常時における音量変化量と対比されることで、断続度合いが数値化される。
【００２２】
本願発明者の研究するところによれば、シャント音の断続度合いが高まると、一周期における音量変化量が大きくなる傾向があることが判明している。よって、算出された音量変化量が、正常時の音量変化量より大きい場合には、断続度合いが高い数値として出力されればよい。一方で、算出された音量変化量が、正常時の音量変化量より小さい場合には、断続度合いが低い数値として出力されればよい。
【００２３】
＜５＞
或いは抽出手段を備える態様では、前記出力手段は、前記一周期音量情報における前記シャント音の音量の最小値と、正常時におけるシャント音の最小値とを対比して、前記断続度合いを数値化する第２数値化手段を有してもよい。
【００２４】
この場合、先ず一周期音量情報におけるシャント音の音量の最小値（例えば、拡張末期の音量）が抽出される。そして、抽出された最小値が、予め記憶された正常時における最小値と対比されることで、断続度合いが数値化される。
【００２５】
本願発明者の研究するところによれば、シャント音の断続度合いが高まると、一周期における音量の最小値が小さくなる傾向があることが判明している。よって、算出された最小値が、正常時の最小値より大きい場合には、断続度合いが低い数値として出力されればよい。一方で、算出された最小値が、正常時の最小値より小さい場合には、断続度合いが高い数値として出力されればよい。
【００２６】
＜６＞
或いは抽出手段を備える態様では、前記出力手段は、前記一周期音量情報における前記シャント音の音量の最大値から最小値に至る傾きに基づいて、前記断続度合いを数値化する第３数値化手段を有してもよい。
【００２７】
この場合、先ず一周期音量情報におけるシャント音の音量の最大値から最小値に至る傾き（言い換えれば、音量が最大値となってから最小値に至るまでの減衰率）が算出される。そして、算出された傾きに基づいて、断続度合いが数値化される。
【００２８】
本願発明者の研究するところによれば、シャント音の断続度合いが高まると、収縮期から拡張期にかけての音量の減衰が急激になる傾向があることが判明している。よって、算出された傾きの変動が大きいほど、断続度合いが高い数値として出力されればよい。
【００２９】
＜７＞
或いは抽出手段を備える態様では、前記出力手段は、前記一周期音量情報における前記シャント音の音量の極大値及び極小値の差分と、正常時におけるシャント音の極大値及び極小値の差分とを対比して、前記断続度合いを数値化する第４数値化手段を有してもよい。
【００３０】
この場合、先ず一周期音量情報におけるシャント音の音量の極大値及び極小値が抽出される。極大値及び極小値は、例えば微分波形を平滑化したもの利用して抽出することができる。続いて、抽出された極大値から極小値が減算され、極大値から極小値までの音量変化量がされる。そして、算出された音量変化量が、予め記憶された正常時における音量変化量と対比されることで、断続度合いが数値化される。
【００３１】
本願発明者の研究するところによれば、シャント音の断続度合いが高まると、極大値から極小値までの音量変化量が大きくなる（具体的には、収縮期から拡張期にかけての音量の変化に歪みが生じる）傾向があることが判明している。よって、算出された音量変化量が、正常時の音量変化量より大きい場合には、断続度合いが高い数値として出力されればよい。一方で、算出された音量変化量が、正常時の音量変化量より小さい場合には、断続度合いが低い数値として出力されればよい。
【００３２】
＜８＞
本実施形態に係るシャント音解析装置の他の態様では、前記出力手段は、前記シャント音情報に基づいて、周波数ごとの時間経過に伴う前記シャント音の音量を示す分布情報を導出する分布情報導出手段を有する。
【００３３】
この態様によれば、分布情報導出手段によって導出された分布情報に基づいて、シャント音の断続度合いを示す情報を出力することが可能となる。ここで特に、本願発明者の研究するところによれば、シャント音の断続感は、周波数ごとのシャント音の音量の分布に大きく依存していることが判明している。よって、導出された分布情報を利用すれば、シャント音の断続度合いを示す情報として、より適切な情報を出力することが可能となる。
【００３４】
＜９＞
上述した分布情報導出手段を備える態様では、前記出力手段は、前記分布情報が示す時間経過に伴う周波数重心の変化量と、正常時における時間経過に伴う周波数重心の変化量とを対比して、前記断続度合いを数値化する第５数値化手段を有してもよい。
【００３５】
この場合、先ず分布情報から時間経過に伴う周波数重心の変化量が算出される。周波数重心の変化量は、例えば最大重心周波数から最小重心周波数を減算することで算出できる。そして、算出された周波数重心の変化量が、予め記憶された正常時における周波数重心の変化量と対比されることで、断続度合いが数値化される。
【００３６】
本願発明者の研究するところによれば、シャント音の断続度合いが高まると、周波数重心の変化量が大きくなる（具体的には、高周波成分が増加する）傾向があることが判明している。よって、算出された周波数重心の変化量が、正常時の周波数重心の変化量より大きい場合には、断続度合いが高い数値として出力されればよい。一方で、算出された周波数重心の変化量が、正常時の周波数重心の変化量より小さい場合には、断続度合いが低い数値として出力されればよい。
【００３７】
＜１０＞
本実施形態に係るシャント音解析装置の他の態様では、前記出力手段は、前記第１から第５数値化手段を少なくとも２種類有しており、前記第１から第５数値化手段が数値化した前記断続度合いの各々を統合的に判定して、前記断続度合いを示す情報を出力する。
【００３８】
この態様によれば、互いに異なる指標を利用して数値化された断続度合いが統合的に判定されるため、より正確な断続度合いを示す情報を出力することができる。統合的な判定する際には、例えば各数値化手段で数値化された断続度合いを正規化した上で、所定の重み付けを行えばよい。なお、出力される断続度合いを示す情報としては、例えば聴感に合うことを重視した１つの値、波形の断続性を重視した１つの値、或いは聴感に合った断続傾向と波形的な断続傾向を示す２つの値等が挙げられる。
【００３９】
＜１１＞
本実施形態に係るシャント音解析方法は、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得工程と、前記取得工程で取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力工程とを備える。
【００４０】
本実施形態に係るシャント音解析方法によれば、上述した本実施形態に係るシャント音解析装置と同様に、シャント音の聴感上の断続度合いを示す情報が出力される。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。
【００４１】
＜１２＞
本実施形態に係るコンピュータプログラムは、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報を取得する取得工程と、前記取得工程で取得した前記シャント音情報に基づいて、前記シャント音の聴感上の断続度合いを示す情報を出力する出力工程とをコンピュータに実行させる。
【００４２】
本実施形態に係るコンピュータプログラムによれば、コンピュータに上述した本実施形態に係るシャント音解析方法の各工程を実行させることができる。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。
【００４３】
＜１３＞
本実施形態に係る記録媒体は、上述したコンピュータプログラムが記録されている。
【００４４】
本実施形態に係る記録媒体によれば、記録されたコンピュータプログラムを実行させることで、シャント音の聴感上の断続度合いを示す情報が出力させることが可能である。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。
【００４５】
本実施形態に係るシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体の作用及び他の利得については、以下に示す実施例において、より詳細に説明する。
【実施例】
【００４６】
以下では、図面を参照してシャント音解析装置の実施例について詳細に説明する。
【００４７】
＜装置構成＞
先ず、図１を参照して、本実施例に係るシャント音解析装置の全体構成について説明する。ここに図１は、実施例に係るシャント音解析装置の全体構成を示すブロック図である。
【００４８】
図１において、本実施例に係るシャント音解析装置は、音声信号入力部１１０と、音声信号解析部１２０と、個体差入力部１３０と、大局的特徴抽出部１４０と、局所的特徴抽出部１５０と、統合判定部１６０と、表示部１７０とを備えて構成されている。
【００４９】
音声信号入力部１１０は、例えば振動センサ等により構成されており、被測定者のシャント形成部位からシャント音を検出する。音声信号入力部１１０で検出されたシャント音は、音声信号として音声信号解析部１２０へと出力される構成となっている。音声信号入力部１１０は、「取得部」の一具体例である。
【００５０】
音声信号解析部１２０は、音声信号入力部１１０から入力された音声信号に対して時間周波数解析を行う。音声信号解析部１２０による解析結果は、大局的特徴抽出部１４０及び局所的特徴抽出部１５０の各々に出力される構成となっている。
【００５１】
個体差入力部１３０は、予め被測定者から取得した音量波形に基づいて、被測定者固有のシャント音を示すパラメータを算出する。個体差入力部１３０で算出された被測定者固有のパラメータは、大局的特徴抽出部１４０及び局所的特徴抽出部１５０の各々に出力される構成となっている。
【００５２】
大局的特徴抽出部１４０は、シャント音の断続感を生む要因のうち、大局的特徴（具体的には、周波数重心の変化量、１心拍区間の音量変化量、及び１心拍区間の音量最小値）を抽出し、各特徴に基づいて断続度合いを示す情報を算出する。なお、大局的特徴抽出部１４０は、大局的特徴を抽出する際に、個体差入力部１３０から入力される被測定者固有の情報を考慮する。大局的特徴抽出部１４０において算出された断続度合いを示す情報は、それぞれ統合判定部１６０に出力される構成となっている。
【００５３】
局所的特徴抽出部１５０は、シャント音の断続感を生む要因のうち、局所的特徴（具体的には、１心拍区間の音量減衰率、及び音量変化の歪み）を抽出し、各特徴に基づいて断続度合いを示す情報を算出する。なお、局所的特徴抽出部１５０は、局所的特徴を抽出する際に、個体差入力部１３０から入力される被測定者固有の情報を考慮する。局所的特徴抽出部１５０において算出された断続度合いを示す情報は、それぞれ統合判定部１６０に出力される構成となっている。
【００５４】
統合判定部１６０は、大局的特徴抽出部１４０及び局所的特徴抽出部１５０から入力されたシャント音の断続度合いを示す情報を統合的に判定し、シャント形成部位における狭窄可能性を示す情報を算出する。統合判定部１６０において算出された狭窄可能性を示す情報は、表示部１７０に出力される構成となっている。統合判定部１６０は、「出力部」の一具体例である。
【００５５】
表示部１７０は、例えば液晶ディスプレイ等として構成されており、統合判定部１６０から出力された狭窄可能性を示す情報（シャント音の断続度合いを示す情報）を、例えば医師等の使用者に視覚的に提示することが可能に構成されている。また、表示部１７０は、上述した、周波数重心の変化量、１心拍区間の音量変化量、１心拍区間の音量最小値、１心拍区間の音量減衰率、及び音量変化の歪み、の何れかを視覚的に提示するようにしてもよい。
【００５６】
＜動作説明＞
次に、本実施例に係るシャント音解析装置の動作について説明する。なお、以下では、本実施例に係るシャント音解析装置が有する各部位のうち、本実施例に特有な部位（即ち、音声信号解析部１２０、個体差入力部１３０、大局的特徴抽出部１４０、局所的特徴抽出部１５０、及び統合判定部１６０）の動作について詳細に説明する。
【００５７】
＜音声信号解析部＞
先ず、図２から図５を参照して、音声信号解析部１２０の動作について具体的に説明する。ここに図２は、時間周波数波形の一例を示すスペクトログラム（その１）であり、図３は、音量解析波形の一例を示すグラフである。また図４は、時間周波数波形の一例を示すスペクトログラム（その２）であり、図５は、周波数重心解析波形の一例を示すグラフである。
【００５８】
図２において、音声信号解析部１２０は、音声信号入力部１１０から入力された音声信号の時間周波数解析を行い、時間周波数波形を取得する。時間周波数波形は、シャント音の周波数ごとのパワーを時系列で示すものである。音声信号解析部１２０は、時間周波数波形からシャント音の音量を算出する。具体的には、音声信号解析部１２０は、以下の数式（１）を用いて、シャント音の音量ｙ（ｔ）を算出する。
【００５９】
【数１】
2020199293.tif 000003
なお、ｆ（ｎ）は周波数であり、ｐ（ｎ）は周波数ｆ（ｎ）におけるパワーである。
【００６０】
図３に示すように、音声信号解析部１２０が算出した音量解析波形（即ち、シャント音の音量の時間変化を示す波形）は、被測定者の脈動に応じた周期的な波形となる。音量解析波形は、１心拍区間ごとに分割された後、大局的特徴抽出部１４０及び局所的特徴抽出部１５０に出力され、それぞれ大局的特徴及び局所的特徴の抽出に利用される。
【００６１】
図４において、音声信号解析部１２０は更に、時間周波数波形からシャント音の周波数重心を算出する。具体的には、音声信号解析部１２０は、以下の数式（２）を用いて、シャント音の周波数重心ｇ（ｔ）を算出する。
【００６２】
【数２】
2020199293.tif 000004
図５に示すように、音声信号解析部１２０が算出した周波数重心解析波形（即ち、シャント音の周波数重心の時間変化を示す波形）は、音量解析波形と同様に、被測定者の脈動に応じた周期的な波形となる。周波数重心解析波形は、大局的特徴抽出部１４０に出力され、大局的特徴（具体的には、周波数重心の変化）の抽出に利用される。
【００６３】
＜個体差入力部＞
次に、図６から図８を参照して、個体差入力部１３０の動作について具体的に説明する。ここに図６は、個人差入力部１３０の動作を示すフローチャートである。また図７は、被測定者の過去から蓄積されたシャント音の最小値の分布を示すヒストグラムであり、図８は、被測定者の過去から蓄積されたシャント音の最大値の分布を示すヒストグラムである。
【００６４】
図６において、個体差入力部１３０は、予め取得した音声信号から被測定者固有のパラメータを算出する。具体的には、個体差入力部１３０では、音声信号が取得されると（ステップＳ１０１）、音量解析によって音量解析波形が取得される（ステップＳ１０２）。音量解析波形は１心拍区間単位で分割され（ステップＳ１０３）、各区間における最大値及び最小値が検出される（ステップＳ１０４）。検出された最大値及び最小値は、メモリに記憶され（ステップＳ１０５）、そのヒストグラムから最も頻度の高かった値が、その被測定者の固有のシャント音の音量最大値及び音量最小値として決定される（ステップＳ１０６）。
【００６５】
図７及び図８において、最小値及び最大値の検出頻度が図に示すような分布となった場合、被測定者の固有のシャント音の音量最小値は約１３７と決定される。同様に、被測定者の固有のシャント音の音量最大値は約４２５と決定される。
【００６６】
このようにして決定された被測定者の固有の音量最大値及び音量最小値は、音量解析波形として得られるシャント音の音量の個人差補正に用いられる。補正音量Ｙ（ｔ）は、音量解析によって得られた音量ｙ（ｔ）の平均値をｙ＿ａｖｅ、最大値をｙ＿ｍａｘ、最小値をｙ＿ｍｉｎとし、被測定者の固有の音量の最大値をＹ＿ｍａｘ、最小値をＹ＿ｍｉｎとすると、以下の数式（３）を用いて算出できる。
【００６７】
Ｙ（ｔ）＝ｙ＿ａｖｅ＋（ｙ（ｔ）−ｙ＿ａｖｅ）×（ｙ＿ｍａｘ−ｙ＿ｍｉｎ）／（Ｙ＿ｍａｘ−Ｙ＿ｍｉｎ）  ・・・（３）
なお、ここでの説明は省略するが、周波数重心についても被測定者ごとの固有の最大値及び最小値が決定されてもよい。周波数重心の個人差補正についても、上述した音量の個人差補正と同様の方法で行うことができる。
【００６８】
＜大局的特徴抽出部＞
次に、図９から図１６を参照して、大局的特徴抽出部１４０の動作について具体的に説明する。なお、以下では、大局的特徴抽出部１４０で抽出可能な複数の特徴について、特徴別に断続度合いを算出する方法を説明する。
【００６９】
＜１心拍区間の周波数重心の変化量＞
先ず、図９から図１２を参照して、１心拍区間の周波数重心の変化量に基づく断続度合いの算出方法について説明する。ここに図９は、正常時における１心拍区間の時間周波数波形の一例を示すスペクトログラムであり、図１０は、正常時における１心拍区間の周波数重心解析波形の一例を示すグラフである。また図１１は、狭窄発生時における１心拍区間の時間周波数波形の一例を示すスペクトログラムであり、図１２は、狭窄発生時における１心拍区間の周波数重心解析波形の一例を示すグラフである。
【００７０】
図９及び図１０において、正常時のシャント音は、低い周波数を多く含む音として検出される。図１０を見ると分かるように、正常時のシャント音の周波数重心の最小値は約６１０Ｈｚであり、最大値は約７００Ｈｚである。よって、正常時におけるシャント音の周波数重心の変化量（即ち、最大値−最小値）は約９０Ｈｚとなる。
【００７１】
図１１及び図１２において、狭窄発生時（即ち、断続感が強い）シャント音は、正常時と比べると、高周波成分を多く含む音として検出される。図１２を見ると分かるように、狭窄発生時のシャント音の周波数重心の最小値は約５９０Ｈｚであり、最大値は約７４０Ｈｚである。よって、狭窄発生時におけるシャント音の周波数重心の変化量（即ち、最大値−最小値）は約１５０Ｈｚとなる。
【００７２】
以上の結果から分かるように、正常時のシャント音と狭窄発生時のシャント音とでは、１心拍区間の周波数重心の変化量に明確な差が生ずる。よって、周波数重心解析によって得られた周波数重心の変化量を正常時の周波数重心の変化量と対比すれば、狭窄発生可能性（即ち、シャント音の断続度合い）を好適に算出できる。
【００７３】
＜１心拍区間の音量変化量＞
次に、図１３及び図１４を参照して、１心拍区間の音量変化量に基づく断続度合いの算出方法について説明する。ここに図１３は、正常時における１心拍区間の音量変化量の一例を示すグラフである。また図１４は、狭窄発生時における１心拍区間の音量変化量の一例を示すグラフである。
【００７４】
図１３において、正常時のシャント音では、１心拍区間における音量の最大値と最小値との差分が小さい。即ち、正常時のシャント音は、収縮期から拡張期にかけての音量変化量が比較的小さい。
【００７５】
図１４において、狭窄発生時のシャント音では、１心拍区間における音量の最大値と最小値との差分が、正常時と比べて大きい。即ち、狭窄発生時のシャント音は、収縮期から拡張期にかけての音量変化量が比較的大きい。
【００７６】
従って、音量解析波形から得られた１心拍区間における音量変化量が、正常時の１心拍区間における音量変化量と比べてどの程度大きいかによって、狭窄発生可能性（即ち、シャント音の断続度合い）を好適に算出できる。
【００７７】
＜１心拍区間の音量最小値＞
次に、図１５及び図１６を参照して、１心拍区間の音量最小値に基づく断続度合いの算出方法について説明する。ここに図１５は正常時における１心拍区間の音量最小値の一例を示すグラフである。また図１６は。狭窄発生時における１心拍区間の音量最小値の一例を示すグラフである。
【００７８】
図１５において、正常時のシャント音では、１心拍区間における音量の最小値が比較的大きい。具体的には、正常時のシャント音は、音量が小さくなる拡張期においても、比較的音量が大きいまま維持される。
【００７９】
図１６において、狭窄発生時のシャント音では、１心拍区間における音量の最小値が、正常時と比べて大きい。具体的には、狭窄発生時のシャント音は、収縮期の音量については正常時に近い値まで大きくなるが、拡張期の音量が大きく減衰してしまう。
【００８０】
従って、音量解析波形から得られた１心拍区間における音量最小値が、正常時の１心拍区間における音量最小値と比べてどの程度小さいかによって、狭窄発生可能性（即ち、シャント音の断続度合い）を好適に算出できる。
【００８１】
＜局所的特徴抽出部＞
次に、図１７から図２０を参照して、局所的特徴抽出部１５０の動作について具体的に説明する。なお、以下では、局所的特徴抽出部１５０で抽出可能な複数の特徴について、特徴別に断続度合いを算出する方法を説明する。
【００８２】
＜１心拍区間の音量減衰率＞
先ず、図１７及び図１８を参照して、１心拍区間の音量減衰率に基づく断続度合いの算出方法について説明する。ここに図１７は、正常時における１心拍区間の音量減衰率の一例を示すグラフである。また図１８は、狭窄発生時における１心拍区間の音量減衰率の一例を示すグラフである。
【００８３】
図１７において、正常時のシャント音では、１心拍区間における音量の減衰が直線的である。具体的には、正常時のシャント音では、音量の最大値から最小値までの傾き（減衰率）が一定である。
【００８４】
図１８において、狭窄発生時のシャント音では、１心拍区間における音量の減衰が急激な部分が生じている。具体的には、狭窄発生時のシャント音では、収縮期末期から拡張期初期にかけて音量が著しく減少し、その後は比較的緩やかに減少する。このため、音量の最大値から最小値までの傾き（減衰率）は一定とならない。
【００８５】
従って、音量解析波形から得られた１心拍区間における音量減衰率を利用すれば、狭窄発生可能性（即ち、シャント音の断続度合い）を好適に算出できる。減衰率ｄｒは、以下の数式（４）を用いて算出することができる。
【００８６】
【数３】
2020199293.tif 000005
なお、Ｎは１心拍区間のデータ数であり、ｍａｘ及びｍｉｎは、それぞれ１心拍区間内の音量の最大値及び最小値である。減衰率ｄｒは、図１８中に破線で示す直線的な減衰率である０．５を基準として、そこからの乖離が大きいほど小さい値として算出される。よって、算出された減衰率ｄｒが０．５よりどの程度小さいかによって、狭窄発生可能性（即ち、シャント音の断続度合い）を好適に算出できる。
【００８７】
＜音量変化の歪み＞
次に、図１９及び図２０を参照して、収縮期から拡張期にかけて発生する音量変化の歪みに基づく断続度合いの算出方法について説明する。ここに図１９は、歪み発生時の音量解析波形の一例を示すグラフである。また図２０は、歪み発生時の音量解析微分波形の一例を示すグラフである。
【００８８】
図１９において、狭窄発生時のシャント音では、収縮期末期から拡張期初期にかけて、音量の変化に歪みが生じることがある。具体的には、図中の破線で囲った部分のように、音量が急激に低下した後、一時的に上昇することがある。よって、音量解析波形における極大値の存在によって、狭窄発生可能性（即ち、シャント音の断続度合い）を算出することができる。
【００８９】
図２０において、音量解析波形における極大値は、音量解析波形の微分波形を平滑化したものを用いることで容易に検出できる。具体的には、音量解析波形の微分波形の極大値と極小値との差分を算出すること好適に検出することができる。
【００９０】
＜統合判定部＞
最後に、統合判定部１６０の動作について具体的に説明する。以下では、大局的特徴である周波数重心の変化量として抽出された特徴量をρ１、１心拍区間の音量変化量として抽出された特徴量をρ２、１心拍区間の音量最小値として抽出された特徴量をρ３、局所的特徴である１心拍区間の音量減衰率として抽出された特徴量をρ４、音量変化の歪みとして抽出された特徴量をρ５として説明する。なお、各直量ρ１〜ρ５は、それぞれ正規化されている（例えば、断続度合いが最も強い場合が０、最も弱い場合（即ち、正常時）が１００となるように調整されている）ものとする。
【００９１】
統合判定部１６０は、各特徴量ρ１〜ρ５に対して重み付けを行って断続度を算出する。具体的には、断続度は、各特徴量ρ１〜ρ５に対応する重みをω１〜ω５として、以下の数式（５）を用いて算出される。
【００９２】
断続度＝ω１×ρ１＋ω２×ρ２＋ω３×ρ３＋ω４×ρ４＋ω５×ρ５  ・・・（５）
ここで、断続度を聴感に合うことを重視した１つの数値として出力する場合には、重みω１〜ω５についてω１、ω２、ω３＞＞ω４、ω５の関係が成立するようにすればよい。このようにして得られる断続度は、例えば人工透析現場において当日の透析が可能か否かを判定するのに適している。
【００９３】
また、断続度を波形の断続性を重視した１つの数値として出力する場合には、重みω１〜ω５についてω４、ω５＞ω２＞ω１、ω３の関係が成立するようにすればよい。このようにして得られる断続度は、例えばエコー診断との比較に適している。
【００９４】
或いは、断続度を聴感に合うことを重視した数値及び波形の断続性を重視した数値の２つの値として出力する場合には、以下の数式（６）及び（７）を利用して別々に算出すればよい。
【００９５】
聴感による断続度＝ω１×ρ１＋ω２×ρ２＋ω３×ρ３  ・・・（６）
波形による断続度＝ω４×ρ４＋ω５×ρ５  ・・・（７）
このようにして得られる断続度は、例えば聴感と波形との両方の経時変化から断続（即ち、狭窄）の進行具合を判定する場合に適している。
【００９６】
なお、統合判定部１６０による統合的な判定を行わずに、各特徴量ρ１〜ρ５を単独で出力するようにしてもよい。この場合、表示部１７０には、例えば各特徴量ρ１〜ρ５に応じた数値が別々に表示されたり、レーダーチャートとして表示されたりする。これにより、ユーザである医師等は、各特徴量ρ１〜ρ５の各々を利用して複数の視点から狭窄診断を行うことができる。或いは、各特徴量ρ１〜ρ５の各々を用いて、医師等が自ら統合的な判断をして、狭窄診断を行うことができる。
【００９７】
以上説明したように、本実施例に係るシャント音解析装置によれば、シャント音の解析結果として、断続度合いを示す情報が出力される。従って、シャント形成部位における狭窄診断を好適に支援することが可能である。
【００９８】
本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴うシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。
【符号の説明】
【００９９】
１１０  音声信号入力部
１２０  音声信号解析部
１３０  個体差入力部
１４０  大局的特徴抽出部
１５０  局所的特徴抽出部
１６０  統合判定部
１７０  表示部
【図１】
2020199293.tif 000006
【図２】
2020199293.tif 000007
【図３】
2020199293.tif 000008
【図４】
2020199293.tif 000009
【図５】
2020199293.tif 000010
【図６】
2020199293.tif 000011
【図７】
2020199293.tif 000012
【図８】
2020199293.tif 000013
【図９】
2020199293.tif 000014
【図１０】
2020199293.tif 000015
【図１１】
2020199293.tif 000016
【図１２】
2020199293.tif 000017
【図１３】
2020199293.tif 000018
【図１４】
2020199293.tif 000019
【図１５】
2020199293.tif 000020
【図１６】
2020199293.tif 000021
【図１７】
2020199293.tif 000022
【図１８】
2020199293.tif 000023
【図１９】
2020199293.tif 000024
【図２０】
2020199293.tif 000025