JP2019111347A - Shunt sound analysis device, shunt sound analysis method, computer program and recording medium - Google Patents
Shunt sound analysis device, shunt sound analysis method, computer program and recording medium Download PDFInfo
- Publication number
- JP2019111347A JP2019111347A JP2019002755A JP2019002755A JP2019111347A JP 2019111347 A JP2019111347 A JP 2019111347A JP 2019002755 A JP2019002755 A JP 2019002755A JP 2019002755 A JP2019002755 A JP 2019002755A JP 2019111347 A JP2019111347 A JP 2019111347A
- Authority
- JP
- Japan
- Prior art keywords
- shunt sound
- shunt
- component
- value
- sound
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- External Artificial Organs (AREA)
Abstract
Description
本発明は、被測定者から取得したシャント音を解析するシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体の技術分野に関する。 The present invention relates to the technical fields of a shunt sound analysis device, a shunt sound analysis method, a computer program and a recording medium for analyzing a shunt sound obtained from a subject.
この種の装置として、被測定者から取得したシャント音を解析して、シャント狭窄等に関する医師の診断を支援する装置が知られている。例えば特許文献１では、周波数特性の山の位置が狭窄度合いに対応することを利用して、微細成分を除去した包絡成分のピーク位置に応じた狭窄度合いを提示する技術が記載されている。 As an apparatus of this type, an apparatus is known which analyzes a shunt sound acquired from a subject and assists a doctor's diagnosis of shunt stenosis and the like. For example, Patent Document 1 describes a technique for presenting the degree of stenosis according to the peak position of the envelope component from which the minute component has been removed, using the fact that the position of the mountain in the frequency characteristic corresponds to the degree of stenosis.
上述した特許文献１に記載されている技術では、シャント音の高周波成分（即ち、比較的高い周波数の成分）を利用して血管の狭窄度合いを判定している。しかしながら、狭窄時には必ずしも周波数の高低によらない特徴的な音色（例えば、風切り音等）も発生する。そして特に、このような音色成分は、例えば一時的な血栓や、静脈弁の影響等、狭窄以外の要因でも発生することが知られている。よって、特徴的な音色成分を考慮しない特許文献１の技術では、狭窄度合いを正確に評価できないという技術的問題点が生ずる。 In the technique described in Patent Document 1 described above, the degree of stenosis of a blood vessel is determined using a high frequency component (i.e., a component having a relatively high frequency) of the shunt sound. However, at the time of narrowing, a characteristic tone (for example, wind noise etc.) that does not necessarily depend on the frequency level is also generated. In particular, it is known that such a tonal component also occurs due to factors other than stenosis, such as temporary thrombosis and the influence of a venous valve. Therefore, in the technique of Patent Document 1 in which the characteristic tonal component is not considered, there arises a technical problem that the degree of stenosis can not be accurately evaluated.
また、特徴的な音色成分は、狭窄の度合いによって。周波数特性の包絡成分上に調波性のピークが重畳して現れ、且つ時間的に変動する。このような特性は、特許文献１の周波数特性の山を利用する手法では表現できないため、仮に上述した音色成分を考慮しようと考えた場合であっても、従来の手法では簡単には利用することができない。 Also, the characteristic timbre component depends on the degree of stenosis. The harmonic peak appears superimposed on the envelope component of the frequency characteristic and fluctuates in time. Such a characteristic can not be expressed by the method using the peak of the frequency characteristic of Patent Document 1, so even if it is considered to consider the above-mentioned timbre component, the conventional method should be easily used. I can not
本発明が解決しようとする課題には上記のようなものが一例として挙げられる。本発明は、被測定者から取得したシャント音を解析して、シャント形成部位の狭窄診断を好適に支援することが可能なシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体を提供することを課題とする。 The problems to be solved by the present invention are as mentioned above. The present invention provides a shunt sound analysis device, a shunt sound analysis method, a computer program, and a recording medium capable of suitably assisting in a stenosis diagnosis of a shunt formation site by analyzing a shunt sound obtained from a subject. To be an issue.
上記課題を解決するためのシャント音解析装置は、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得手段と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出手段と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力手段とを備える。 The shunt sound analysis device for solving the above-mentioned subject is produced by the acquisition means which acquires shunt sound information about the shunt sound of a shunt formation part of a person to be measured, and the shunt sound information by narrowing of the blood vessel of the person to be measured. The apparatus may include an extraction unit that extracts a feature component indicating a feature of the shunt sound, and an output unit that outputs evaluation information related to the evaluation of the shunt sound based on the feature component.
上記課題を解決するためのシャント音解析方法は、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得工程と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出工程と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力工程とを備える。 The shunt sound analysis method for solving the above problems includes an acquisition step of acquiring shunt sound information on a shunt sound of a shunt formation site of a subject, and generation of a stenosis of a subject's blood vessel from the shunt sound information. And an output step of outputting evaluation information related to the evaluation of the shunt sound based on the characteristic component.
上記課題を解決するためのコンピュータプログラムは、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得工程と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出工程と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力工程とをコンピュータに実行させる。 A computer program for solving the above problems includes an acquisition step of acquiring shunt sound information on a shunt sound of a shunt formation site of a subject, and the shunt generated by narrowing of a blood vessel of the subject from the shunt sound information. A computer is caused to execute an extraction step of extracting a feature component indicating a feature of a sound and an output step of outputting evaluation information related to the evaluation of the shunt sound based on the feature component.
上記課題を解決するための記録媒体は、上述したコンピュータプログラムが記録されている。 The computer program described above is recorded in a recording medium for solving the above problems.
＜１＞
本実施形態に係るシャント音解析装置は、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得手段と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出手段と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力手段とを備える。
<1>
The shunt sound analysis apparatus according to the present embodiment includes acquisition means for acquiring shunt sound information on a shunt sound of a shunt-formed portion of a subject, and the shunt generated by narrowing of a blood vessel of the subject from the shunt sound information. An extraction unit that extracts a feature component that indicates a feature of a sound, and an output unit that outputs evaluation information related to the evaluation of the shunt sound based on the feature component.
本実施形態に係る第１のシャント音解析装置の動作時には、先ず取得手段により、被測定者のシャント形成部位周辺から、シャント音に関するシャント音情報が取得される。なお、ここでの「シャント音」とは、血液を体外に取り出すためのシャント形成部位周辺において取得される血流音であり、被測定者の脈拍に同期した音である。シャント音の取得は各種センサを用いて行えばよく、その取得方法が特に限定されるものではない。また「シャント音情報」とは、シャント音に関する各種パラメータを含む情報であって、例えば音量や周波数等の時間変化などを含んでいる。 At the time of operation of the first shunt sound analysis device according to the present embodiment, first, the acquisition unit acquires shunt sound information on the shunt sound from the periphery of the shunt formation portion of the subject. Here, the “shunt sound” is a blood flow sound acquired around the shunt formation site for taking blood out of the body, and is a sound synchronized with the pulse of the subject. Acquisition of a shunt sound may be performed using various sensors, and the acquisition method is not particularly limited. The “shunt sound information” is information including various parameters related to the shunt sound, and includes, for example, temporal changes in volume and frequency.
シャント音情報が取得されると、抽出手段により、シャント音情報から特徴成分が抽出される。特徴成分は、被測定者の血管の狭窄によって生じるシャント音の特徴を示す成分であり、例えばシャント音情報を示す時間周波数波形から、周波数特性の包絡成分と、シャント音特有のピッチ成分の周波数範囲を超えた微細成分を除去したものである。 When the shunt sound information is acquired, the extraction unit extracts a feature component from the shunt sound information. The characteristic component is a component indicating the characteristic of the shunt sound caused by the narrowing of the blood vessel of the subject, for example, from the time frequency waveform showing shunt sound information, the envelope component of the frequency characteristic and the frequency range of the pitch component unique to the shunt sound. The removal of fine components exceeding
特徴成分が抽出されると、出力手段により、シャント音の評価に関連する評価情報が出力される。具体的には、出力手段は、抽出された特徴成分に基づいて評価情報を生成し、生成した評価情報を外部のモニタ等に出力する。なお、ここでの「評価情報」は、シャント形成部位の狭窄度合いを評価した結果を示す情報であってもよいし、シャント形成部位の狭窄度合いを評価するための情報であってもよい。より具体的には、狭窄度合いを直接的に示す数値等が評価情報として出力されてもよいし、医師等が狭窄度合いを判断するための一又は複数のパラメータが評価情報として出力されてもよい。 When the feature component is extracted, the output means outputs evaluation information related to the evaluation of the shunt sound. Specifically, the output means generates evaluation information based on the extracted feature component, and outputs the generated evaluation information to an external monitor or the like. Here, the “evaluation information” may be information indicating the result of evaluating the degree of stenosis of the shunt formation site, or may be information for evaluating the degree of stenosis of the shunt formation site. More specifically, a numerical value or the like directly indicating the degree of stenosis may be output as the evaluation information, or one or more parameters for a doctor or the like to determine the degree of stenosis may be output as the evaluation information .
本実施形態に係る評価情報は、上述したように、シャント音情報から抽出された特徴成分に基づくものである。ここで特に、シャント形成部位の血管に狭窄が発生した場合は、シャント音の高周波成分が強くなる傾向にあるが、本願発明者の研究するところによれば、周波数の高低によらない特徴的な音色成分も生じることが判明している。例えば、狭窄発生時には、ヒュー又はビューという擬音で表現される音響的な成分や、シャンシャン又はシャッシャッという擬音で表現される成分が生ずる。 As described above, the evaluation information according to the present embodiment is based on the feature component extracted from the shunt sound information. Here, in particular, when a stenosis occurs in a blood vessel where a shunt is formed, the high frequency component of the shunt sound tends to be strong, but according to the research of the inventor of the present invention, the characteristic is not dependent on the frequency level. It has been found that a tonal component is also generated. For example, when a stenosis occurs, an acoustic component represented by an onomatopoeic sound such as hue or view or a component represented by an onomatopoeic sound such as a shanghai or a shatter occurs.
本実施形態に係るシャント音解析装置によれば、上述した特徴的な音色成分を考慮した評価情報が出力される。よって、例えばシャント音の周波数の高低のみで狭窄を評価する場合と比べると、より正確な評価を実現できる。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。 According to the shunt sound analysis device according to the present embodiment, evaluation information in consideration of the above-described characteristic tone color component is output. Therefore, a more accurate evaluation can be realized, for example, as compared to the case where the narrowing is evaluated only by the frequency of the shunt sound. Therefore, it is possible to preferably support a stenosis diagnosis of the shunt formation site.
＜２＞
本実施形態に係るシャント音解析装置の一態様では、前記抽出手段は、前記シャント音情報に対しケプストラム変換処理を行い、第１の次数以下のケフレンシー及び第２の次数以上のケフレンシーを除去してフーリエ変換することで、前記特徴成分を抽出する。
<2>
In one aspect of the shunt sound analysis device according to the present embodiment, the extraction unit performs cepstrum conversion processing on the shunt sound information to remove quefrance less than the first order and quefrance more than the second order. The feature component is extracted by Fourier transform.
この態様によれば、シャント音情報が取得されると、先ずケプストラム変換処理が行われる。ケプストラム変換処理は、例えばシャント音情報を示す時間周波数波形に対して、逆フーリエ変換を行うことで実現される。 According to this aspect, when shunt sound information is acquired, first, cepstrum conversion processing is performed. The cepstrum conversion process is realized, for example, by performing an inverse Fourier transform on a time frequency waveform indicating shunt sound information.
続いて、ケプストラム変換されたシャント音情報からは、第１の次数以下のケフレンシー及び第２の次数以上のケフレンシーが除去される。ここで、「第２の次数」とは、シャント音特有のピッチ成分の周波数を超えた微細成分を除去するために設定される次数であり、事前のシミュレーション等により最適値が設定されている。また、「第１の次数」とは、シャント音の周波数特性（例えば、周波数ごとのシャント音の大きさを示す波形によって示される特性）の包絡成分を除去するために設定される次数であり、第１の次数と同様に、事前のシミュレーション等により最適値が設定されている。なお、「第１の次数」及び「第２の次数」の最適値は、サンプリング周波数やフレーム長さに応じた値とすればよい。 Subsequently, the quefrency lower than the first order and the quefrance higher than the second order are removed from the cepstrum-converted shunt sound information. Here, the “second order” is an order set to remove a fine component exceeding the frequency of the pitch component unique to the shunt sound, and the optimum value is set by a simulation or the like in advance. The “first order” is an order set to remove the envelope component of the frequency characteristic of the shunt sound (for example, the characteristic indicated by the waveform indicating the magnitude of the shunt sound for each frequency), As in the first order, optimum values are set by prior simulation or the like. The optimum values of the “first order” and the “second order” may be values corresponding to the sampling frequency or the frame length.
第１の次数及び第２の次数に応じたケフレンシーが除去されたシャント音情報は、フーリエ変換されることで、元の周波数波形へと変換される。このようにして得られる波形は、狭窄発生時に生ずる特徴的な音色成分（即ち、特徴成分）を抽出したものとなる。 The shunt sound information from which the quefency according to the first order and the second order has been removed is converted to the original frequency waveform by Fourier transformation. The waveform obtained in this manner is one in which a characteristic tonal component (i.e., a characteristic component) generated at the occurrence of a stenosis is extracted.
以上のように、ケプストラム変換及び所定次数のケフレンシー除去を利用すれば、シャント音情報の特徴成分を好適に抽出することが可能である。 As described above, it is possible to preferably extract the feature component of the shunt sound information by using the cepstral transformation and the quefrency removal of the predetermined order.
＜３＞
本実施形態に係るシャント音解析装置の他の態様では、前記出力手段は、前記特徴成分における時間軸上で隣り合うフレームの相互相関を示す調波性成分値を算出し、前記評価情報として出力する。
<3>
In another aspect of the shunt sound analysis device according to the present embodiment, the output means calculates a harmonic component value indicating a cross-correlation of adjacent frames on the time axis in the feature component, and outputs it as the evaluation information. Do.
この態様によれば、シャント音情報から特徴成分が抽出されると、特徴成分における時間軸上で隣り合うフレーム（例えば、時刻ｎ−１のフレーム及び時刻ｎのフレーム）の相互相関を示す調波性成分値が算出される。このようにして算出される調波性成分値は、シャント音における時間方向に持続的な成分量を示すパラメータである。 According to this aspect, when the feature component is extracted from the shunt sound information, the harmonics indicating the cross-correlation of the adjacent frames (for example, the frame at time n-1 and the frame at time n) in the feature component Sex component values are calculated. The harmonic component value calculated in this manner is a parameter that indicates the amount of component that is persistent in the time direction in the shunt sound.
本願発明者の研究するところによれば、狭窄発生時に生ずる特徴的な音色は、時間方向に連続した成分を有することが判明している。よって、調波性成分値を算出すれば、瞬間的な（即ち、持続しない）特徴成分の影響が軽減され、特徴的な音色成分の大きさ（強さ）に対応する値を得ることができる。従って、調波性成分値を評価情報として出力すれば、狭窄の発生を好適に評価することが可能である。 According to the research of the inventor of the present application, it has been found that the characteristic tone produced at the occurrence of a stenosis has a continuous component in the time direction. Therefore, if the harmonic component value is calculated, the influence of the instantaneous (that is, not sustained) feature component is reduced, and a value corresponding to the size (intensity) of the characteristic tone component can be obtained. . Therefore, if the harmonic component value is output as the evaluation information, it is possible to appropriately evaluate the occurrence of the stenosis.
＜４＞
本実施形態に係るシャント音解析装置の他の態様では、前記出力手段は、前記特徴成分の周波数重心値を算出し、前記評価情報として出力する。
<4>
In another aspect of the shunt sound analysis device according to the present embodiment, the output means calculates a frequency gravity center value of the feature component and outputs it as the evaluation information.
この態様によれば、シャント音情報から特徴成分が抽出されると、特徴成分の周波数重心値が算出される。周波数重心値は、特徴的な音色成分の高さに相当する値であり、例えば狭窄の度合いや、その他の要因によって異なる値として算出される。従って、周波数重心値を評価情報として出力すれば、狭窄の発生を好適に評価することが可能である。 According to this aspect, when the feature component is extracted from the shunt sound information, the frequency centroid value of the feature component is calculated. The frequency centroid value is a value corresponding to the height of the characteristic tone color component, and is calculated as a different value depending on, for example, the degree of narrowing and other factors. Therefore, if the frequency gravity center value is output as the evaluation information, it is possible to preferably evaluate the occurrence of the stenosis.
＜５＞
本実施形態に係るシャント音解析装置の他の態様では、前記出力手段は、前記特徴成分のピーク値から前記特徴成分の平均値を差し引いた音色傾向値を算出し、前記評価情報として出力する。
<5>
In another aspect of the shunt sound analysis device according to the present embodiment, the output means calculates a timbre tendency value obtained by subtracting the average value of the feature components from the peak value of the feature components, and outputs it as the evaluation information.
この態様によれば、シャント音情報から特徴成分が抽出されると、先ず特徴成分のピーク値（所定時間の最大ピーク値）及び平均値（所定時間の平均値）が算出される。そして、特徴成分のピーク値から平均値が差し引かれ、音色傾向値として算出される。 According to this aspect, when the feature component is extracted from the shunt sound information, first, the peak value (maximum peak value of the predetermined time) and the average value (average value of the predetermined time) of the feature component are calculated. Then, the average value is subtracted from the peak value of the feature component, and the tone value is calculated.
音色傾向値は、例えばヒューヒュー等の音響的な成分と、シャッシャッ等の雑音的な成分とのいずれの傾向が強いかを示す値である。具体的には、音色傾向値が高い場合には、音響的な成分の傾向が強いと判定でき、音色傾向値が低い場合には、雑音的な成分の傾向が強いと判定できる。従って音色傾向値を評価情報として出力すれば、狭窄の発生を好適に評価することが可能である。 The timbre tendency value is a value indicating which of an acoustic component such as Hugh and a noise component such as a shatter is more likely. Specifically, when the tone tendency value is high, it can be determined that the tendency of the acoustic component is strong, and when the tone tendency value is low, it can be determined that the tendency of the noise component is strong. Therefore, if the timbre tendency value is output as the evaluation information, it is possible to preferably evaluate the occurrence of the narrowing.
＜６＞
本実施形態に係るシャント音解析方法は、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得工程と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出工程と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力工程とを備える。
<6>
The shunt sound analysis method according to the present embodiment includes an acquisition step of acquiring shunt sound information on a shunt sound of a shunt-formed portion of a subject, and the shunt generated due to a narrowing of a blood vessel of the subject from the shunt sound information. And an output step of outputting evaluation information related to the evaluation of the shunt sound based on the feature component.
本実施形態に係るシャント音解析方法によれば、上述した本実施形態に係るシャント音解析装置と同様に、狭窄によって生じる特量成分に基づいた評価情報が出力される。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。 According to the shunt sound analysis method according to the present embodiment, as in the above-described shunt sound analysis device according to the present embodiment, the evaluation information based on the special component generated due to the narrowing is output. Therefore, it is possible to preferably support a stenosis diagnosis of the shunt formation site.
＜７＞
本実施形態に係るコンピュータプログラムは、被測定者のシャント形成部位のシャント音に関するシャント音情報を取得する取得工程と、前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出工程と、前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力工程とをコンピュータに実行させる。
<7>
The computer program according to the present embodiment includes an acquiring step of acquiring shunt sound information on a shunt sound of a shunt formation site of a subject, and from the shunt sound information, the shunt sound generated by the narrowing of the blood vessel of the subject. A computer is caused to execute an extraction step of extracting a feature component indicating a feature, and an output step of outputting evaluation information related to the evaluation of the shunt sound based on the feature component.
本実施形態に係るコンピュータプログラムによれば、コンピュータに上述した本実施形態に係るシャント音解析方法の各工程を実行させることができる。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。 According to the computer program of the present embodiment, it is possible to cause the computer to execute the respective steps of the shunt sound analysis method of the present embodiment described above. Therefore, it is possible to preferably support a stenosis diagnosis of the shunt formation site.
＜８＞
本実施形態に係る記録媒体は、上述したコンピュータプログラムが記録されている。
<8>
The above-described computer program is recorded in a recording medium according to the present embodiment.
本実施形態に係る記録媒体によれば、記録されたコンピュータプログラムを実行させることで、狭窄によって生じる特量成分に基づいた評価情報を出力させることが可能である。従って、シャント形成部位の狭窄診断を好適に支援することが可能である。 According to the recording medium of the present embodiment, it is possible to output the evaluation information based on the special component caused by the stenosis by executing the recorded computer program. Therefore, it is possible to preferably support a stenosis diagnosis of the shunt formation site.
本実施形態に係るシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体の作用及び他の利得については、以下に示す実施例において、より詳細に説明する。 The functions and other gains of the shunt sound analysis device, the shunt sound analysis method, the computer program, and the recording medium according to the present embodiment will be described in more detail in the examples shown below.
以下では、図面を参照してシャント音解析装置の実施例について詳細に説明する。 Hereinafter, embodiments of a shunt sound analysis apparatus will be described in detail with reference to the drawings.
＜装置構成＞
先ず、図１を参照して、本実施例に係るシャント音解析装置の全体構成について説明する。ここに図１は、実施例に係るシャント音解析装置の全体構成を示すブロック図である。
<Device configuration>
First, the overall configuration of a shunt sound analysis device according to the present embodiment will be described with reference to FIG. FIG. 1 is a block diagram showing the overall configuration of a shunt sound analysis device according to the embodiment.
図１において、本実施例に係るシャント音解析装置は、シャント音入力部１１０と、音声信号解析処理部１２０と、パラメータ演算部１３０と、評価情報演算出力部１４０と、表示部１５０とを備えて構成されている。 In FIG. 1, the shunt sound analysis apparatus according to the present embodiment includes a shunt sound input unit 110, an audio signal analysis processing unit 120, a parameter calculation unit 130, an evaluation information calculation output unit 140, and a display unit 150. Is configured.
シャント音入力部１１０は、「取得手段」の一具体例であり、例えば振動センサ等から入力されたシャント音に対し、所定のサンプリング周波数Ｆｓでアナログデジタル変換を行い、シャント音波形を出力する。シャント音入力部１１０で生成されたシャント音波形は、音声信号解析処理部１２０へと出力される構成となっている。 The shunt sound input unit 110 is one specific example of the “acquisition unit”, and performs analog-to-digital conversion at a predetermined sampling frequency Fs, for example, on a shunt sound input from a vibration sensor or the like, and outputs a shunt sound waveform. The shunt sound waveform generated by the shunt sound input unit 110 is output to the audio signal analysis processing unit 120.
音声信号解析処理部１２０は、シャント音入力部１１０から入力されたシャント音波形に対して各種解析処理（例えば、短時間フーリエ変換や対数変換等）を行い、時間周波数解析波形を出力する。音声信号解析処理部１２０による解析結果は、パラメータ演算部１３０に出力される構成となっている。 The audio signal analysis processing unit 120 performs various analysis processing (for example, short time Fourier transformation, logarithmic conversion, and the like) on the shunt sound waveform input from the shunt sound input unit 110, and outputs a time frequency analysis waveform. The analysis result by the audio signal analysis processing unit 120 is output to the parameter calculation unit 130.
パラメータ演算部１３０は、「抽出手段」の一具体例であり、音声信号解析処理部１２０から入力された時間周波数解析波形に対して各種演算処理（例えば、逆フーリエ変換、ケフレンシーカット、フーリエ変換等）を行うことで、時間特徴成分特性波形を得る。またパラメータ演算部１３０は、時間特徴成分特性波形から、調波性成分値、周波数重心値、及び音色傾向値を算出する。パラメータ演算部１３０によって算出された調波性成分値、周波数重心値、及び音色傾向値は、評価情報演算出力部１４０に夫々出力される構成となっている。 The parameter operation unit 130 is one specific example of the “extraction means”, and performs various operation processing (for example, inverse Fourier transform, querency cut, Fourier transform, etc.) on the time frequency analysis waveform input from the audio signal analysis processing unit 120. Etc., to obtain a time characteristic component characteristic waveform. The parameter calculator 130 also calculates the harmonic component value, the frequency centroid value, and the timbre tendency value from the time feature component characteristic waveform. The harmonic component value, the frequency centroid value, and the timbre tendency value calculated by the parameter calculation unit 130 are output to the evaluation information calculation output unit 140, respectively.
評価情報演算出力部１４０は、「出力手段」の一具体例であり、パラメータ演算部１３０の出力値、及び出力値のピーク値や時間平均値を用いて、被測定者の血管の狭窄度合いを評価するための評価情報を演算し出力する。評価情報演算出力部１４０によって演算された評価情報は、表示部１５０に出力される構成となっている。 The evaluation information calculation output unit 140 is a specific example of “output means”, and uses the output value of the parameter calculation unit 130 and the peak value or time average value of the output value to determine the degree of stenosis of the blood vessel of the subject. Calculates and outputs evaluation information for evaluation. The evaluation information calculated by the evaluation information calculation output unit 140 is output to the display unit 150.
表示部１５０は、例えばモニタ等として構成されており、評価情報演算出力部１４０から入力された評価情報を、例えば医師等の装置使用者に視覚的に提示することが可能に構成されている。 The display unit 150 is configured, for example, as a monitor, and is configured to be able to visually present the evaluation information input from the evaluation information calculation output unit 140 to a device user such as a doctor.
＜動作説明＞
次に、本実施例に係るシャント音解析装置の具体的な動作について説明する。なお、以下では、本実施例に係るシャント音解析装置が有する各部位のうち、本実施例に特有な部位（具体的には、音声信号解析処理部１２０、パラメータ演算１３０、評価情報演算出力部１４０、及び表示部１５０）の動作について詳細に説明する。
<Description of operation>
Next, a specific operation of the shunt sound analysis device according to the present embodiment will be described. In the following, among the portions of the shunt sound analysis device according to the present embodiment, portions unique to the present embodiment (specifically, the audio signal analysis processing unit 120, the parameter calculation 130, the evaluation information calculation output unit The operation of the display unit 140 and the display unit 150 will be described in detail.
＜音声信号解析処理部＞
先ず、図２及び図３を参照して、音声信号解析処理部１２０の動作について詳細に説明する。ここに図２は、シャント音波形の一例を示す波形図である。また図３は、時間周波数波形の一例を示すスペクトログラムである。
<Voice signal analysis processing unit>
First, the operation of the audio signal analysis processing unit 120 will be described in detail with reference to FIGS. 2 and 3. Here, FIG. 2 is a waveform diagram showing an example of a shunt sound wave form. FIG. 3 is a spectrogram showing an example of a time frequency waveform.
図２に示すように、音声信号解析処理部１２０には、シャント音入力部１１０で生成されたシャント音波形が入力される。なお、シャント音波形）は、図を見ても分かるように、被測定者のシャント音形性部位から取得された音声の振幅について、その時間変動を示す波形である。 As shown in FIG. 2, a shunt sound waveform generated by the shunt sound input unit 110 is input to the audio signal analysis processing unit 120. The shunt sound waveform is a waveform showing time fluctuation of the amplitude of the sound acquired from the shunt sound form part of the subject, as can be seen from the figure.
音声信号解析処理部１２０は、シャント音波形の時刻ｎの値Ｘ（n）に対して、長さＮフレーム単位で、Ｎポイント短時間フーリエ変換と、対数変換とを施し、時間周波数解析波形Ｐｌｏｇ［ｎ，ｋ］を算出する。 The audio signal analysis processing unit 120 performs N-point short-time Fourier transformation and logarithmic transformation on the value X (n) at time n of the shunt sound waveform in units of length N frame to obtain a time frequency analysis waveform Plog. Calculate [n, k].
具体的には、音声信号解析処理部１２０は、以下の数式（１）を用いて、短時間フーリエ変換を行う。 Specifically, the audio signal analysis processing unit 120 performs short-time Fourier transform using the following equation (1).
図３に示すように、音声信号解析処理部１２０によって演算された時間周波数解析波形Ｐｌｏｇ［ｎ，ｋ］は、各周波数に対応する音圧の時間変動を示す波形である。 As shown in FIG. 3, the time-frequency analysis waveform Plog [n, k] calculated by the sound signal analysis processing unit 120 is a waveform indicating time fluctuation of sound pressure corresponding to each frequency.
＜パラメータ演算部＞
次に、図４から図９を参照して、パラメータ演算部１３０の動作について詳細に説明する。
<Parameter operation unit>
Next, the operation of the parameter operation unit 130 will be described in detail with reference to FIGS. 4 to 9.
＜特徴成分の抽出＞
先ず、図４から図６を参照して、パラメータ演算部１３０による特徴成分の抽出について詳細に説明する。ここに図４は、時刻ｎにおける周波数波形の一例を示す波形図であり、図５は、時刻ｎにおける特徴成分特性波形の一例を示す波形図である。また図６は、時間特徴成分特性波形の一例を示すスペクトログラムである。
<Extraction of feature components>
First, extraction of feature components by the parameter operation unit 130 will be described in detail with reference to FIGS. 4 to 6. FIG. 4 is a waveform chart showing an example of the frequency waveform at time n, and FIG. 5 is a waveform chart showing an example of the characteristic component characteristic waveform at time n. FIG. 6 is a spectrogram showing an example of a time feature component characteristic waveform.
図４に示すように、時間周波数解析波形Ｐｌｏｇ［ｎ，ｋ］の時刻ｎにおける周波数波形は、周波数毎の音圧を示す波形として表される。パラメータ演算部１３０は、このような時間周波数解析波形Ｐｌｏｇ［ｎ，ｋ］の各フレームに各種演算処理を行い、特徴成分を抽出する。 As shown in FIG. 4, the frequency waveform at time n of the time frequency analysis waveform Plog [n, k] is expressed as a waveform indicating the sound pressure for each frequency. The parameter operation unit 130 performs various operation processing on each frame of such a time frequency analysis waveform Plog [n, k] to extract a feature component.
具体的には、パラメータ演算部１３０は、以下の数式（３）を用いて、時間周波数解析波形Ｐｌｏｇ［ｎ，ｋ］の逆フーリエ変換を行う。 Specifically, the parameter operation unit 130 performs inverse Fourier transform of the time-frequency analysis waveform Plog [n, k] using the following equation (3).
なお、次数ｌｉｆｔ＿ｌｏｗ及びｌｉｆｔ＿ｈｉｇｈの最適値は、サンプリング周波数Ｆｓやフレーム長さＮに応じて定まる値であり、事前のシミュレーション等により最適値が設定されている。ここでの次数ｌｉｆｔ＿ｌｏｗは、「第１の次数」の一具体例であり、次数ｌｉｆｔ＿ｈｉｇｈは、「第２の次数」の一具体例である。 The optimum values of the orders lift_low and lift_high are values determined according to the sampling frequency Fs and the frame length N, and the optimum values are set by simulation in advance. The order lift_low here is one specific example of the “first order”, and the order lift_high is one specific example of the “second order”.
次に、パラメータ演算部１３０は、以下の数式（５）を用いてフーリエ変換を行い、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］を演算する。 Next, the parameter calculation unit 130 performs Fourier transform using the following equation (5) to calculate a time feature component characteristic waveform Fpeaks [n, k].
図６に示すように、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］は、各周波数に対応する音圧の微細成分の時間変動を示す波形である。 As shown in FIG. 6, the time characteristic component characteristic waveform FPeaks [n, k] is a waveform that indicates the time fluctuation of the minute component of the sound pressure corresponding to each frequency.
＜特徴成分パラメータの演算＞
次に、図７及び図８を参照して、パラメータ演算部１３０による特徴成分パラメータの演算について詳細に説明する。ここに図７は、パラメータ演算部の具体的な構成を示すブロック図である。また図８は、調波性成分値、周波数重心値及び音色傾向値の一例を示すタイムチャートである。
<Calculation of feature component parameters>
Next, calculation of feature component parameters by the parameter calculation unit 130 will be described in detail with reference to FIGS. 7 and 8. FIG. 7 is a block diagram showing a specific configuration of the parameter operation unit. FIG. 8 is a time chart showing an example of the harmonic component value, the frequency gravity center value, and the timbre tendency value.
図７において、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］は、パラメータ演算部１３０の調波性成分演算部１３１、周波数重心演算部１３２、及び音色傾向演算部１３３に夫々入力される。 In FIG. 7, the time feature component characteristic waveform Fpeaks [n, k] is input to the harmonic component calculation unit 131, the frequency centroid calculation unit 132, and the timbre tendency calculation unit 133 of the parameter calculation unit 130.
調波性成分演算部１３１は、以下の数式（６）を用いて、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］から、特徴成分パラメータの１つである調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］を算出する。 The harmonic component computing unit 131 uses the following equation (6) to extract the harmonic component value FPeaksValue [n], which is one of the feature component parameters, from the time feature component characteristic waveform FPeaks [n, k]. calculate.
周波数重心演算部１３２は、以下の数式（７）を用いて、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］から、特徴成分パラメータの１つである周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］を算出する。 The frequency centroid calculation unit 132 calculates a frequency centroid value FPeaksCentroid [n], which is one of the feature component parameters, from the time feature component characteristic waveform FPeaks [n, k] using the following equation (7).
音色傾向演算部１３３は、以下の数式（８）を用いて、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］から、特徴成分パラメータの１つである音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］を算出する。 The timbre tendency operation unit 133 calculates the timbre tendency value FPeaksSoundColor [n] which is one of the feature component parameters from the time feature component characteristic waveform FPeaks [n, k] using the following equation (8).
図８に示すように、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］は、特徴的な音色成分の大きさに相当する値として算出される。周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］は、特徴的な音色成分の高さに相当する値として算出される。音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］は、特徴的な音色成分の音色の傾向を示す値として算出される。 As shown in FIG. 8, the harmonic component value FPeaksValue [n] is calculated as a value corresponding to the size of the characteristic tone color component. The frequency centroid value FPeaksCentroid [n] is calculated as a value corresponding to the height of the characteristic tone color component. The timbre tendency value FPeaksSoundColor [n] is calculated as a value indicating the tendency of the timbre of the characteristic timbre component.
＜演算結果を用いた評価＞
次に、図９を参照して、上述した演算によって得られた各種パラメータを用いた狭窄度合いの評価方法について、具体的に説明する。ここに図９は、パラメータ演算部における演算結果を複数のシャント音毎に示す比較一覧図である。
<Evaluation using calculation results>
Next, with reference to FIG. 9, a method for evaluating the degree of stenosis using various parameters obtained by the above-described calculation will be specifically described. FIG. 9 is a comparison list showing the calculation results in the parameter calculation unit for each of a plurality of shunt sounds.
図９において、シャント音Ａは、正常なシャント音の一例である。シャント音Ａは、時間特徴成分特性波形ＦＰｅａｋｓ［ｎ，ｋ］を見ても分かるように、特徴的な特性が現れない。このため、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］、及び音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］のいずれも極めて低い値で安定している。このように、特徴成分パラメータが低い値で安定している場合には、狭窄は発生していない（或いは、極めて狭窄度合いが小さい）と評価することができる。 In FIG. 9, the shunt sound A is an example of a normal shunt sound. As seen from the time feature component characteristic waveform Fpeaks [n, k], the shunt sound A does not exhibit characteristic characteristics. Therefore, all of the harmonic component value FPeaksValue [n], the frequency gravity center value FPeaksCentroid [n], and the timbre tendency value FPeaksSoundColor [n] are stable at extremely low values. As described above, when the feature component parameter is stable at a low value, it can be evaluated that no stenosis has occurred (or the degree of stenosis is extremely small).
シャント音Ｂは、シャッシャッという特徴的な音を含むシャント音の一例である。シャント音Ｂは、シャント音Ａと比べると、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］が高い値として算出されている。また、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］も高い値で推移しており、心拍に同期して高さが揺れ動く様子が分かる。音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］からは、雑音的な傾向が強いことが分かる。この結果、シャント音Ｂには、雑音的な特徴成分が含まれており、狭窄が発生している可能性があると評価できる。 The shunt sound B is an example of a shunt sound including a characteristic sound called shatter. As compared with the shunt sound A, the shunt sound B is calculated as a higher value of the harmonic component value FPeaksValue [n]. Further, the frequency gravity center value FPeaksCentroid [n] also changes at a high value, and it can be seen that the height swings in synchronization with the heartbeat. From the tone tendency value FPeaksSoundColor [n], it can be seen that the noise tendency is strong. As a result, the shunt sound B contains a noise characteristic component, and it can be evaluated that a narrowing may occur.
シャント音Ｃは、ビュービューという特徴的な音を含むシャント音の一例である。シャント音Ｃは、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］が極めて高い値として算出されている。また、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］は極めて高い値で推移している。音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］からは、音響的な傾向が強いことが分かる。この結果、シャント音Ｃには音響的な特徴成分が多く含まれており、狭窄が発生している可能性が高いと評価できる。 The shunt sound C is an example of a shunt sound including a characteristic sound called a view. The shunt sound C is calculated as an extremely high value of the harmonic component value FPeaksValue [n]. In addition, the frequency gravity center value FPeaksCentroid [n] changes at a very high value. The tonal tendency value FPeaksSoundColor [n] indicates that the acoustic tendency is strong. As a result, the shunt sound C contains many acoustic feature components, and it can be evaluated that the possibility of the occurrence of a stenosis is high.
シャント音Ｄは、低い唸りのような特徴的な音を含むシャント音の一例である。シャント音Ｄは、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］が高い値として算出されている。また、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］は断続的に激しく変化している。音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］からは、音響的な傾向が強いことが分かる。この結果、シャント音Ｃには音響的な特徴成分が含まれており、狭窄が発生している可能性があると評価できる。 The shunt sound D is an example of a shunt sound including a characteristic sound such as low whistling. The shunt sound D is calculated as a high value of the harmonic component value FPeaksValue [n]. In addition, the frequency centroid value FPeaksCentroid [n] changes intermittently and sharply. The tonal tendency value FPeaksSoundColor [n] indicates that the acoustic tendency is strong. As a result, the shunt sound C contains an acoustic feature component, and it can be evaluated that a narrowing may occur.
なお、上述した評価方法はあくまで一例であり、同一の演算結果から異なる評価を行っても構わない。 In addition, the evaluation method mentioned above is an example to the last, and you may evaluate different from the same calculation result.
＜評価情報演算出力部＞
次に、評価情報演算出力部１４０の動作について詳細に説明する。
<Evaluation information calculation output unit>
Next, the operation of the evaluation information calculation output unit 140 will be described in detail.
評価情報演算出力部１４０は、パラメータ演算部１３０で算出された調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］、及び音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］の各々を数値化して、狭窄度合いを評価するための（或いは、狭窄度合いを示す）評価情報として出力する。 The evaluation information calculation output unit 140 digitizes each of the harmonic component value FPeaksValue [n], the frequency gravity center value FPeaksCentroid [n], and the timbre tendency value FPeaksSoundColor [n] calculated by the parameter calculation unit 130, The information is output as evaluation information for evaluating the degree (or indicating the degree of stenosis).
具体的には、評価情報演算出力部１４０は、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］を、０ｄＢ〜２０ｄＢの範囲で０〜１００に正規化して出力する。また、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］を、１５０Ｈｚ〜２ｋＨｚの範囲で０〜１００に正規化して出力する。また、評価情報演算出力部１４０は、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］及び周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］の各々について、ピーク値及び時間平均値を演算して出力する。 Specifically, the evaluation information calculation output unit 140 normalizes the harmonic component value FPeaksValue [n] to 0 to 100 in the range of 0 dB to 20 dB and outputs the result. Further, the frequency centroid value FPeaksCentroid [n] is normalized to 0 to 100 in the range of 150 Hz to 2 kHz and output. Further, the evaluation information calculation output unit 140 calculates and outputs a peak value and a time average value for each of the harmonic component value FPeaksValue [n] and the frequency gravity center value FPeaksCentroid [n].
評価情報演算出力部１４０は、音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］を、０ｄＢ〜１０ｄＢの範囲で０〜１００に正規化して出力する。また、評価情報演算出力部１４０は、音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］について、ピーク値及び時間平均値を演算して出力する。 The evaluation information calculation output unit 140 normalizes the timbre tendency value FPeaksSoundColor [n] to 0 to 100 in the range of 0 dB to 10 dB and outputs it. In addition, the evaluation information calculation output unit 140 calculates and outputs a peak value and a time average value for the tone tendency value FPeaksSoundColor [n].
＜表示部＞
次に、図１０から図１３を参照して、表示部１５０の動作について詳細に説明する。ここに図１０から図１３は夫々、表示部における表示例を示す平面図である。
<Display section>
Next, the operation of the display unit 150 will be described in detail with reference to FIGS. 10 to 13. FIGS. 10 to 13 are plan views showing display examples in the display unit.
図１０及び図１１に示すように、表示部１５０の表示領域１５５には、例えば特徴成分の大きさ（即ち、調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］）の平均スコア（即ち、正規化された時間平均値）と、平均スコアについての印象マップが表示される。なお、ここでの平均スコアは、３０以上の場合に狭窄が発生している可能性が高いと評価できるようなパラメータとなっている。 As shown in FIGS. 10 and 11, in the display area 155 of the display unit 150, for example, the average score (ie, normalized time) of the size of the feature component (ie, the harmonic component value FPeaksValue [n]). An impression map for the average score and the average score is displayed. The average score here is a parameter that can be evaluated as having a high possibility of the occurrence of a stenosis in the case of 30 or more.
具体的には、図１０の例では、平均スコア“５”と、特徴成分が殆ど現れていない印象マップが表示されている。この場合、例えば狭窄が発生している可能性は低く、当日の透析に問題はなさそうだと判断できる。一方、図１１の例では、平均スコア“７０”と、特徴成分が強く現れた印象マップが表示されている。この場合、例えば狭窄が発生している可能性が高く、医師に連絡が必要であると判断できる。このような比較的簡易な表示態様は、聴診時などにおいてシャント音の聞こえ方に対する判断が難しい場合に有効である。 Specifically, in the example of FIG. 10, an average score of “5” and an impression map in which a feature component hardly appears are displayed. In this case, for example, the possibility that a stenosis has occurred is low, and it can be determined that there is no problem with dialysis on the day. On the other hand, in the example of FIG. 11, an average score “70” and an impression map in which the feature component appears strongly are displayed. In this case, for example, there is a high possibility that a stenosis has occurred, and it can be determined that the doctor needs to be contacted. Such a relatively simple display mode is effective when it is difficult to determine how to hear the shunt sound at the time of auscultation or the like.
図１２及び図１３に示すように、表示部１５０の表示領域１５５には、例えば調波性成分値ＦＰｅａｋｓＶａｌｕｅ［ｎ］、周波数重心値ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］、及び音色傾向値ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］の各々のピーク値及び平均スコア、並びに解析波形が表示されてもよい。このように評価情報を詳細に表示すれば、より正確で細かな診断（例えば、スコアの上昇が狭窄による影響なのか、或いは他の要因による影響なのかの診断）を行うことも可能となる。 As shown in FIGS. 12 and 13, in the display area 155 of the display unit 150, for example, each of the harmonic component value FPeaksValue [n], the frequency gravity center value FPeaksCentroid [n], and the timbre tendency value FPeaksSoundColor [n]. Peak values and average scores, as well as analysis waveforms may be displayed. By displaying the evaluation information in detail in this manner, it is also possible to make a more accurate and detailed diagnosis (for example, a diagnosis as to whether an increase in score is an influence of a stenosis or an influence of other factors).
具体的には、図１２の例では、特徴的な音色が大きく変動及び断続しているため、狭窄を疑うべきであると判断できる。また、図１３の例では、特徴的な音が見られるが、音が連続しているようなので、他の所見を併せて診断すべきと判断できる。 Specifically, in the example of FIG. 12, it is possible to judge that the narrowing should be suspected because the characteristic timbre is greatly fluctuated and interrupted. Further, in the example of FIG. 13, characteristic sounds can be seen, but since sounds seem to be continuous, it can be determined that other findings should be diagnosed together.
なお、表示部１５０は、上述した表示例以外の態様で表示を行っても構わない。また、装置の使用者が適切な表示態様を適宜選択できるように構成されてもよい。 Note that the display unit 150 may perform display in modes other than the above-described display example. Further, the user of the device may be configured to be able to select an appropriate display mode as appropriate.
以上説明したように、本実施例に係るシャント音解析装置によれば、取得したシャント音情報に基づいて適切な評価情報が出力される。従って、シャント形成部位における狭窄診断を好適に支援することが可能である。 As described above, according to the shunt sound analysis device according to the present embodiment, appropriate evaluation information is output based on the acquired shunt sound information. Therefore, it is possible to preferably support the diagnosis of stenosis at the shunt formation site.
本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴うシャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。 The present invention is not limited to the above-described embodiment, and can be appropriately modified without departing from the scope or spirit of the invention as can be read from the claims and the specification as a whole. An apparatus, a shunt sound analysis method, a computer program and a recording medium are also included in the technical scope of the present invention.
１１０ シャント音入力部
１２０ 音声信号解析処理部
１３０ パラメータ演算部
１３１ 調波性成分演算部
１３２ 周波数重心演算部
１３３ 音色傾向値演算部
１４０ 評価情報演算出力部
１５０ 表示部
１５５ 表示領域
ＦＰｅａｋｓＶａｌｕｅ［ｎ］ 調波性成分値
ＦＰｅａｋｓＣｅｎｔｒｏｉｄ［ｎ］ 周波数重心値
ＦＰｅａｋｓＳｏｕｎｄＣｏｌｏｒ［ｎ］ 音色傾向値
110 Shunt sound input unit 120 Audio signal analysis processing unit 130 Parameter operation unit 131 Harmonic component operation unit 132 Frequency center of gravity operation unit 133 Timbre tendency value operation unit 140 Evaluation information operation output unit 150 Display unit 155 Display area FPeaksValue [n] Wave component value FPeaksCentroid [n] Frequency centroid value FPeaksSoundColor [n] Tone tendency value
Claims (1)
前記シャント音情報から、前記被測定者の血管の狭窄によって生じる前記シャント音の特徴を示す特徴成分を抽出する抽出手段と、
前記特徴成分に基づいて、前記シャント音の評価に関連する評価情報を出力する出力手段と
を備えることを特徴とするシャント音解析装置。 Acquisition means for acquiring shunt sound information on the shunt sound of the shunt-formed portion of the subject;
Extracting means for extracting, from the shunt sound information, a feature component indicating a feature of the shunt sound caused by a narrowing of a blood vessel of the subject;
An output means for outputting evaluation information related to the evaluation of the shunt sound based on the characteristic component.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019002755A JP2019111347A (en) | 2019-01-10 | 2019-01-10 | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019002755A JP2019111347A (en) | 2019-01-10 | 2019-01-10 | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017524292A Division JP6467044B2 (en) | 2015-06-22 | 2015-06-22 | Shunt sound analysis device, shunt sound analysis method, computer program, and recording medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2019111347A true JP2019111347A (en) | 2019-07-11 |
Family
ID=67221880
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019002755A Pending JP2019111347A (en) | 2019-01-10 | 2019-01-10 | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2019111347A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010029434A (en) * | 2008-07-29 | 2010-02-12 | Tom-Medic Co Ltd | Vasoconstriction level determination program, recording medium, information terminal device, vasoconstriction level determination system and vasoconstriction level determination method |
| JP2014008263A (en) * | 2012-06-29 | 2014-01-20 | Univ Of Yamanashi | Shunt constriction diagnostic support system and method, array shaped sound collection sensor device, and successive segmentation self organized map forming device, method and program |
-
2019
- 2019-01-10 JP JP2019002755A patent/JP2019111347A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010029434A (en) * | 2008-07-29 | 2010-02-12 | Tom-Medic Co Ltd | Vasoconstriction level determination program, recording medium, information terminal device, vasoconstriction level determination system and vasoconstriction level determination method |
| JP2014008263A (en) * | 2012-06-29 | 2014-01-20 | Univ Of Yamanashi | Shunt constriction diagnostic support system and method, array shaped sound collection sensor device, and successive segmentation self organized map forming device, method and program |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11116478B2 (en) | Diagnosis of pathologies using infrasonic signatures | |
| EP3678553B1 (en) | Diagnosis of pathologies using infrasonic signatures | |
| US20200051585A1 (en) | Voice activity detection method and apparatus | |
| JP6467044B2 (en) | Shunt sound analysis device, shunt sound analysis method, computer program, and recording medium | |
| WO2018155384A1 (en) | Detection device | |
| Backstrom et al. | Time-domain parameterization of the closing phase of glottal airflow waveform from voices over a large intensity range | |
| JP5077847B2 (en) | Reverberation time estimation apparatus and reverberation time estimation method | |
| Millette et al. | Signal processing of heart signals for the quantification of non-deterministic events | |
| Buder et al. | Quantitative and graphic acoustic analysis of phonatory modulations | |
| JP2019111347A (en) | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium | |
| Majerus et al. | Bruit-enhancing phonoangiogram filter using sub-band autoregressive linear predictive coding | |
| WO2016165220A1 (en) | Traffic safety detection method and apparatus | |
| CN106683665B (en) | Method and system for analyzing musical scale of audio | |
| Naufal et al. | Advancement on automatic blood pressure measurement using auscultatory method: A literature review | |
| Scalassara et al. | Autoregressive decomposition and pole tracking applied to vocal fold nodule signals | |
| WO2016139802A1 (en) | Shunt murmur analysis device, shunt murmur analysis method, computer program, and recording medium | |
| WO2016207950A1 (en) | Shunt murmur analysis device, shunt murmur analysis method, computer program, and recording medium | |
| JP6501917B2 (en) | Body sound analysis device, body sound analysis method, computer program and recording medium | |
| JP6782940B2 (en) | Tongue position / tongue habit judgment device, tongue position / tongue habit judgment method and program | |
| TWI226600B (en) | Nasal detection method and device thereof | |
| JP6630149B2 (en) | Body sound analyzer, body sound analysis method, computer program, and recording medium | |
| RU165114U1 (en) | DEVICE FOR ASSESSING EMOTIONAL HUMAN VOICE BY VOICE | |
| JP6298527B2 (en) | Body sound analysis apparatus, body sound analysis method, computer program, and recording medium | |
| JP7303921B2 (en) | Shunt sound analysis device, shunt sound analysis method, computer program and recording medium | |
| Zhang et al. | Temporal analysis, manipulation, and resynthesis of musical vibrato |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190110 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200204 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20200728 |