1. Technical Field
The present disclosure relates to a bioacoustic sound testing device, more particularly, to a bioacoustic sound testing device for supporting the estimation of the state of a living body by measuring bioacoustic sound.
2. Background Art
The treatment of respiratory diseases, such as asthma, often takes a long time. For the treatment of asthma, it is basically necessary to take anti-inflammatory agent for controlling the inflammation of the respiratory tract every day and also take bronchodilator for controlling the stenosis of the respiratory tract at the time of an attack, thereby maintaining the respiratory tract in a properly controlled state for a long time. A medical compound of anti-inflammatory agent and bronchodilator is also used as a daily medicine in some cases.
A comprehensive judgment based on the experience of a medical specialist is required for the judgment of the state of the respiratory tract, and support using objective indexes obtained from testing devices is extensively desired to make a comprehensive and appropriate judgment. As a testing device of this kind, for example, the bioacoustic sound testing device described in WO2012/060107 is available.
However, in WO2012/060107, further improvements are required.
In one general aspect, the techniques are disclosed here that a bioacoustic sound testing device comprising:
a band power calculator that calculates a band power in a predetermined period for each of a plurality of predetermined frequency bands on the basis of a bioacoustic sound;
an index value calculator that calculates at least one index value on the basis of the reference power and the band powers of the plurality of predetermined frequency bands.
These general and specific aspects may be implemented using a system, a method, and a computer program, and any combination of systems, methods, and computer programs.
With the present disclosure, an index value properly reflecting the state of the respiratory tract can be obtained while the number of the sensors for measuring the bioacoustic sound is reduced. Furthermore, since the number of the sensors can be reduced, the operation of the apparatus is simplified.
The inventors found that the techniques described in WO2012/060107 cause following problems.
It is known that respiratory sound becomes stronger as respiratory airflow is higher and becomes weaker as physical constitution is larger, and the index value of the standard of the respiratory sound in which the influence of these is reduced is known. In the bioacoustic sound testing device described in WO2012/060107, respiratory sound is measured at two points on the cervical region and the chest, and the ratio of the powers at the cervical region and the chest is used to reduce the influence of the respiratory airflow and the physical constitution. Furthermore, one of a plurality of index values is obtained by predicting the state in a high band from that in the band reflecting respiratory airflow at the cervical region, and a predictive residue is calculated, whereby the influence of the respiratory airflow is reduced.
However, the bioacoustic sound testing device described in WO2012/060107 described above has a problem in that the operation thereof becomes complicated because two sensors are used for the cervical region and the chest. In the comparison between the cervical region and the chest, the testing at the chest close to the peripheral respiratory tract is more suited for the detection of the small change in the state of the respiratory tract, whereby it is assumed that the state of the respiratory tract can be estimated by using only the sensor for the chest; however, since the respiration sound at the chest is weaker, the influence due to factors other than the state of the respiratory tract that is desired to be known is significant, whereby an index value properly reflecting the state of the respiratory tract is not obtained. Hence, the two sensors for the chest and the cervical region are required to obtain the index value properly reflecting the state of the respiratory tract.
Therefore, the inventors considered following improvements.
(1) In a first aspect of the present disclosure, there is provided a bioacoustic sound testing device according to the present disclosure includes a band power calculator that calculates a band power in a predetermined period for each of a plurality of predetermined frequency bands on the basis of a bioacoustic sound; a reference power calculator that calculates a reference power on the basis of the band power of the band at least partially including any of frequencies ranging from 250 Hz or more to 1050 Hz or less and the band power of the band at least partially including any of frequencies ranging from 1250 Hz or more to 1550 Hz or less among from the band powers of the plurality of predetermined frequency bands; and an index value calculator that calculates at least one index value on the basis of the reference power and the band powers of the plurality of predetermined frequency bands.
With the above-mentioned configuration, an index value properly reflecting a state of a respiratory tract is obtained.
(2) In the above aspect, exemplarily, the index value calculator includes: a predicted value calculator that calculates a predicted value for each of the bands on the basis of at least the reference power; and a corrected power calculator that calculates a corrected power for each of the bands by subtracting the predicted value for each of the bands from the band power for each of the bands. The index value calculator calculates an index value on the basis of the corrected power.
(3) In the above aspect, exemplarily, the bands where correction is performed include 700 Hz and 1400 Hz.
(4) In the above aspect, exemplarily, the index value calculator performs correction for one band on the basis of each of two or more kinds of reference powers and sets each corrected power as the index value.
(5) In the above aspect, exemplarily, the two or more kinds of reference powers include a reference power based on weighted addition of respiratory airflow reflecting power and specific change reflecting power and a reference power based on weighted subtraction between the respiratory airflow reflecting power and the specific change reflecting power.
(6) In the above aspect, exemplarily, the device further includes a display that displays the index value.
(7) In the above aspect, exemplarily, the display plots two or more kinds of index values along a time axis.
(8) In the above aspect, exemplarily, the display displays two kinds of index values as a two-dimensional map and displays the two kinds of index values along time series by connecting measurement points.
(9) In the above aspect, exemplarily, the index value calculator calculates the index value on the basis of a linear sum of two or more corrected powers.
(10) In the above aspect, exemplarily, the predicted value calculator predicts a standard value of the band power corresponding to the reference power on the basis of a measurement sample in a standard state.
(11) In the above aspect, exemplarily, the band power calculator includes: a frequency analyzer that performs frequency analysis for a measured bioacoustic sound; a period designator that designates a period to be analyzed among from a plurality of respiration periods; and a period representative band power calculator that calculates a representative band power in the period designated by the period designator for the plurality of predetermined frequency bands on the basis of the bioacoustic sound subjected to frequency analysis using the frequency analyzer.
(12) In the above aspect, exemplarily, a bioacoustic sound measurer for measuring the bioacoustic sound measures the bioacoustic sound at one point on a chest of a living body.
(12) In the above aspect, exemplarily, the band reflecting the specific change is a band in which a difference between a power after a medication of bronchodilator and a power after the medication of anti-inflammatory agent is large.
(13) In the above aspect, exemplarily, the reference power is an average value of the powers of continuous wide bands including the band reflecting the respiratory airflow and the band reflecting the specific change.
(14) In the above aspect, exemplarily, the reference power is an average value of the powers of the bands including frequencies ranging from 150 Hz or more to 1550 Hz or less.
A preferred embodiment according to the present disclosure will be described below in detail referring to the drawings.
The frequency analyzer 202 performs frequency analysis using, for example, fast Fourier transform. The period designator 203 designates a period to be used for the analysis performed using the period representative band power calculator 204 described later. The period to be designated may merely be designated under the conditions that inhalation is performed, that the level of the sound is appropriate and that noise is low and so on, and the professional knowledge of a medical doctor is not required. Extraction of the period can be made automatically by executing a program so that the above-mentioned conditions are satisfied. Furthermore, the respiration period may be designated manually through GUI (Graphical User Interface) or the like.
The period representative band power calculator 204 calculates the power representing the period designated by the period designator 203, for example, for each of 15 bands (bands having central frequencies of 100 Hz, 200 Hz, 300 Hz, . . . , 1300 Hz, 1400 Hz and 1500 Hz and having a bandwidth of 100 Hz). For example, powers in analysis frame unit for frequency analysis are averaged throughout analysis frames included in the period. In this case, since the width of each band is 100 Hz, the band having a central frequency of 100 Hz ranges from 50 to 150 Hz, the band having a central frequency of 200 Hz ranges from 150 to 250 Hz, the band having a central frequency of 300 Hz ranges from 250 to 350 Hz, . . . . , the band having a central frequency of 1300 Hz ranges from 1250 to 1350 Hz, the band having a central frequency of 1400 Hz ranges from 1350 to 1450 Hz, and the band having a central frequency of 1500 Hz ranges from 1450 to 1550 Hz. The “15” frequency bands are just taken as examples and the number of the frequency bands may be larger or smaller than 15. Moreover, the bandwidth is not limited to 100 Hz but may be variable for each band. Also, in a similar way, the period representative band power calculator 204 calculates the period representative band powers of eight bands having central frequencies of 300 Hz, 400 Hz, 500 Hz, 600 Hz, 700 Hz, 800 Hz, 900 Hz and 1000 Hz. These eight bands are bands reflecting the speed of respiration. The “eight” frequency bands are just taken as examples and the number of the frequency bands may be larger or smaller than eight.
Still further, in a similar way, the period representative band power calculator 204 calculates the period representative band powers of three bands having central frequencies of 1300 Hz, 1400 Hz and 1500 Hz. These three bands are bands reflecting a specific change in the living body. More specifically, the difference between the power after the medication of bronchodilator and the power after the medication of anti-inflammatory agent is large in the bands. The “three” frequency bands reflecting the specific change in the living body are just taken as examples and the number of the frequency bands may be larger or smaller than three.
After calculating the band powers of the 15 bands at the central frequencies of 100 to 1500 Hz, the period representative band power calculator 204 outputs the calculated band powers of the 15 bands to the corrected power calculator 503 of the index value calculator 5. Furthermore, after calculating the band powers of the eight bands at the central frequencies of 300 to 1000 Hz, the period representative band power calculator 204 outputs the calculated band powers of the eight bands to the respiratory airflow reflecting power calculator 3 included in the reference power calculator 8. Moreover, after calculating the band powers of the three bands at the central frequencies of 1300 to 1500 Hz, the period representative band power calculator 204 outputs the calculated band powers of the three bands to the specific change reflecting power calculator 4 included in the reference power calculator 8.
The respiratory airflow reflecting power calculator 3 calculates respiratory airflow reflecting power by averaging the band powers of the eight bands serving as bands reflecting the respiratory airflow. In this embodiment, although the bands reflecting the respiratory airflow in the respiratory airflow reflecting power calculator 3 have frequencies ranging from 250 Hz or more to 1050 Hz or less (from a central frequency of 300 to 1000 Hz), the bands may at least partially include any of frequencies ranging from 250 or more to 1050 Hz or less. The specific change reflecting power calculator 4 calculates specific change reflecting power by averaging the band powers of the three bands serving as bands reflecting the specific change other than the respiratory airflow. In this embodiment, although the bands reflecting the specific change in the specific change reflecting power calculator 4 have frequencies ranging from 1250 or more to 1550 Hz or less (from a central frequency of 1300 to 1500 Hz), the bands may at least partially include any of frequencies ranging from 1250 or more to 1550 Hz or less. Furthermore, operation other than averaging may also be used as the operation for calculating the respiratory airflow reflecting power and the specific change reflecting power.
The index value calculator 5 is used to calculate at least one index value from the respiratory airflow reflecting power calculated by the respiratory airflow reflecting power calculator 3, the specific change reflecting power calculated by the specific change reflecting power calculator 4 and the band power calculated by the band power calculator 2, and includes a prediction coefficient storage 501 for holding a prediction coefficient for each band to predict a predicted value for each band; a predicted value calculator 502 for calculating a predicted value for each band on the basis of the respiratory airflow reflecting power, the specific change reflecting power and the prediction coefficient for each band; the corrected power calculator 503 for calculating corrected power for each band by subtracting the predicted value for each band from the representative band power for each band; and an index value calculator 504 for calculating the index value on the basis of the corrected power for each band. Personal profile (for example, physical constitution, age, and medication situation) may be input to the predicted value calculator 502. Although respiratory sound is generally influenced by the level of the power of each band depending on the respiratory airflow, the degree of the influence differs for each band, and the influence of factors other than the respiratory airflow also differs for each band. In this embodiment, the power in a band susceptible to the influence of the respiratory airflow or the like is set to reference power and a standard power value corresponding to the reference power is predicted so that the desired information on the state of the respiratory tract remains in the difference (predictive residue) between the predicted value and the measured value and so that the influence of the other factors is eliminated. The reference power may merely be calculated on the basis of the respiratory airflow reflecting power calculated by the respiratory airflow reflecting power calculator 3 and the specific change reflecting power calculated by the specific change reflecting power calculator 4. In other words, the reference power is calculated by the reference power calculator 8.
The display 6 has a displaying means, such as a liquid crystal display device or an organic EL display device, and displays the index values in a visible mode. In the case that the index values are displayed together with the samples of the past index values stored, the comparison among the index values can be performed. For example, the display plots two or more index values along a time axis. Furthermore, the display displays two index values as a two-dimensional map and displays the index values along time series by connecting measurement points.
Next, the operation of the bioacoustic sound testing device 1 according to this embodiment will be described.
After calculating the band power of the bioacoustic sound, the band power calculator 2 outputs the band powers of the bands having central frequencies of 300 to 1000 Hz and strongly reflecting the respiratory airflow among the 15 bands and the band powers of the bands having central frequencies of 1300 to 1500 Hz and strongly reflecting the specific change in the living body (for example, the change due to respiratory tract stenosis or the differences between respiratory tract inflammation and respiratory tract stenosis) to the reference power calculator 8. On the basis of these band powers, the reference power calculator 8 calculates the reference power (at step S2). In other words, the reference power calculator 8 executes a reference power calculation step. The calculated reference power is then output to the index value calculator 5.
More specifically, the respiratory airflow reflecting power calculator 3 calculates the respiratory airflow reflecting power on the basis of the band powers of the bands having central frequencies of 300 to 1000 Hz. In other words, the respiratory airflow reflecting power calculator 3 calculates the respiratory airflow reflecting power by averaging the respiratory airflow values of the bands reflecting the respiratory airflow among the 15 bands. The calculated respiratory airflow reflecting power is then output to the index value calculator 5.
Furthermore, the specific change reflecting power calculator 4 calculates the specific change reflecting power on the basis of the band powers of the bands having central frequencies of 1300 to 1500 Hz. In other words, the specific change reflecting power calculator 4 calculates the specific change reflecting power by averaging the specific change values of the bands reflecting the specific change among the 15 bands. The calculated specific change reflecting power is then output to the index value calculator 5.
On the basis of the reference power (the respiratory airflow reflecting power and the specific change reflecting power) calculated by the reference power calculator 8 and the band power calculated by the band power calculator 2, the index value calculator 5 calculates one or more index values (at step S3). In other words, the index value calculator 5 sequentially executes the step of reading the prediction coefficient for each band to predict the predicted value for each band (prediction coefficient reading step); the step of calculating the predicted value for each band on the basis of the respiratory airflow reflecting power, the specific change reflecting power and the prediction coefficient for each band (predicted value calculation step); the step of calculating the corrected power for each band by subtracting the predicted value for each band from the representative band power for each band (corrected power value calculation step); and the step of calculating at least one index value on the basis of the corrected power for each band (index value calculation step). The index value calculated by the index value calculator 5 is obtained as a single corrected power value or as the linear sum of a plurality of corrected powers.
It may be possible that a plurality of inhalation periods are designated at step S1 and that a representative value (for example, a median value) is determined from among the index values for the designated periods at step S3.
Next, the features of the bioacoustic sound testing device 1 according to this embodiment will be described in detail.
In the case that the power of the band in which the difference between the respiratory sound improved by daily controller medication and the respiratory sound obtained after the respiratory tract was dilated by only bronchodilator medication is large is set to the specific change reflecting power and that the power obtained by linear operation in the direction of adding the respiratory airflow reflecting power to the specific change reflecting power is set to the reference power, the index value emphasizing the change in the respiratory sound due to the improvement in chronic inflammation can be calculated. Conversely, in the case that the power obtained by linear operation in the direction of subtracting the specific change reflecting power from the respiratory airflow reflecting power is set to the reflecting power, the index value emphasizing the change in the respiratory sound due to the improvement in respiratory tract stenosis can be calculated.
A method for obtaining the prediction coefficient will herein be described. Prior to the use of the bioacoustic sound testing device according to this embodiment, it is assumed that numerous respiratory sound samples of medical cases not diagnosed with asthma by medical doctors or medical cases diagnosed with asthma but being good in the state of the respiratory tract have been collected and adopted as the standard respiratory sound. The standard value (the value of the band power value in the standard state) can be obtained by analyzing the relationship between the power at the single band of 700 Hz and the reference power based on the plurality of bands. Various methods, such as regression analysis, are available as analysis methods, and a parameter representing physical constitution, such as height, and a parameter representing age may be used additionally in addition to the reference power. Furthermore, it may also be possible that the reference power is divided into a plurality of zones depending on the level thereof, that representative values based on the least-squares method are determined for the respective zones, and that the representative values are connected. These numerical values representing the relationship between the reference power and the standard value of the band power, such as the coefficients by regression analysis and the representative values for the respective zones, are referred to as prediction coefficients.
Since the power of a single band usually becomes larger as the reference power becomes larger, a monotonous increase relationship is obtained. Although the relationship is generally linear, the relationship is non-linear if it is observed minutely. First, an example in the case that the respiratory airflow reflecting power is set to the reference power is shown. Referring to
The reference power of the correction (a): base_pow=Fpow
where base_pow represents the reference power and Fpow represents the respiratory airflow reflecting power.
Next, a case in which two kinds of reference powers are used will be described as an example. Two or more kinds of reference powers can be used for the correction. When a case (“white circle P3”) is considered in which the expression of the reference power is modified and the reference power has increased, as the reference power increases, the predicted value also increases, but the power of a single band remains unchanged, whereby the predictive residue S decreases. Conversely, the predictive residue S increases when the reference power decreases. Expressions in which the specific change reflecting power is added or subtracted to increase or decrease the reference power are the expressions of the reference power of the correction (b) and the reference power of the correction (c). However, it is necessary to prepare a predicted value corresponding to the expression of the reference power after the modification. The correction (b) is used to reduce the influence of respiratory airflow and the influence of respiratory tract stenosis, and the correction (c) is used to reduce the influence of respiratory airflow and to emphasize the influence of respiratory tract stenosis.
The reference power of the correction (b): base_pow=Fpow+0.3×Npow
The reference power of the correction (c): base_pow=Fpow−0.2×Npow
where base_pow represents the reference power, Fpow represents the respiratory airflow reflecting power, and Npow represents the specific change reflecting power.
The band in which the specific change reflecting power is decreased by the effect of bronchodilator (the improvement of respiratory tract stenosis) is the band in which the specific change reflecting power is increased by the effect of anti-inflammatory agent (the improvement of respiratory tract inflammation); hence, in the correction (b) in which the weighted addition of the specific change reflecting power to the respiratory airflow reflecting power is performed, the power after the correction tends to increase due to the improvement of respiratory tract stenosis and the power after the correction tends to decrease due to the improvement of respiratory tract inflammation. Furthermore, in the correction (c) in which the weighted subtraction of the specific change reflecting power from the respiratory airflow reflecting power is performed, the power after the correction tends to decrease due to the improvement of respiratory tract stenosis and the power after the correction tends to increase due to the improvement of respiratory tract inflammation.
The weights of the calculation expressions of the correction (b) and the correction (c) were adjusted to the extent that the change due to β2 stimulant becomes 0 in average in the correction (b) and to the extent that the change due to anti-inflammatory agent becomes 0 in average in the correction (c). The weights of the correction (b) and the correction (c) were set, for example, to “0.3” for the correction (b) and to “0.2” for the correction (c). Hence, it is assumed that the index value (b) obtained by the correction (b) mainly represents the degree of respiratory tract inflammation and that the index value (c) obtained by the correction (c) mainly represents the degree of respiratory tract stenosis.
As described above, the correction can be performed by preparing a plurality of reference power calculation expressions and by obtaining standard values corresponding thereto. In addition, standard values are also obtained similarly beforehand for each single band other than 700 Hz and correction is performed. Also in the correction (b), the degree of respiratory tract stenosis is expressed in a band in which the power increases due to the deterioration of respiratory tract stenosis, for example, the 1400 Hz band. In other words, the standard values exist in the number of individual bands for each definition of the reference power. However, for bands that are not used in the end, operation is not required to be performed and parameters are not required to be stored.
As a modification of the correction (b), the average value of the powers of wide bands (for example, 300 to 1500 Hz and 200 to 1500 Hz) including the respiratory airflow reflecting power band and the specific change reflecting power band may be set as the reference power. In addition, in the case that noise tends to be included in a band, the band may be eliminated or interpolation may be performed by using the powers of the adjacent bands.
Furthermore, the effects obtained by the correction (c) and the correction (b) are also obtained respectively when the power value at 700 Hz after the correction and the power value at 1400 Hz after the correction are subjected to weighted addition or weighted subtraction.
The index value can be used properly depending on the contents of medication.
Furthermore, the “15” frequency bands are just taken as examples and the frequency bands may be defined otherwise.
When the correction is performed similarly at 1400 Hz as in the example at 700 Hz shown in
As the daily control state of the respiratory tract, it is desired to know whether the control of inflammation is successful, instead of knowing whether the respiratory tract has been dilated using bronchodilator. Bronchodilator literally has the action of dilating the respiratory tract and changes the physical characteristics of the respiratory tract serving as an acoustic tube, whereby the respiratory sound is also changed. However, the control state of the inflammation and the control state of the stenosis are mixed in the state of the respiratory tract appearing in the respiratory sound. In this respect, since the influence of respiratory tract stenosis can be eliminated in the present disclosure, it is possible to know whether the control of the inflammation is successful.
As described above, the bioacoustic sound testing device 1 according to this embodiment calculates the band power in the 15 frequency bands having a bandwidth of 100 Hz and ranging from 100 to 1500 Hz (central frequencies) on the basis of the bioacoustic sound obtained at one point on the chest; calculates the band powers of the eight bands ranging from 300 to 1000 Hz reflecting the respiratory airflow and the band powers of the three bands ranging from 1300 to 1500 Hz reflecting the specific change of the living body, from among the 15 bands ranging from 100 to 1500 Hz; calculates the respiratory airflow reflecting power by averaging the band powers of the bands reflecting the respiratory airflow; calculates the specific change reflecting power by averaging the band powers of the bands reflecting the specific change; obtains the reference power by the linear sum of the calculated respiratory airflow reflecting power and the calculated specific change reflecting power; calculates the difference between the predicted value corresponding to the reference power and the band power for each band as the corrected power; and calculates at least one index value on the basis of the corrected power, whereby an index value properly reflecting the state of the respiratory tract can be obtained while the number of the sensors for measuring the bioacoustic sound is reduced. Furthermore, since the number of the sensors can be reduced, the operation of the apparatus is simplified.
Although linear operations have been described as examples of the operations in this embodiment, the operations may include non-linear operations, such as logarithmic and exponential operations.
Furthermore, the bioacoustic sound testing device 1 according to this embodiment can be achieved using a proprietary circuit configuration or can be achieved using a configuration including a computer. It is also possible to divide processing between the terminal and the cloud.
Moreover, the program shown in
Although the present disclosure has been described in detail referring to the specific embodiments, it is obvious to those skilled in the art that the present disclosure can be changed and modified variously without departing from the spirit and scope of the present disclosure.
The present disclosure has an effect capable of obtaining an index value properly reflecting the state of the respiratory tract while the number of the sensors for measuring bioacoustic sound is reduced, and is applicable to medical apparatuses for supporting the estimation of the state of a living body.
Number | Date | Country | Kind |
---|---|---|---|
2013-155811 | Jul 2013 | JP | national |
This application is a continuation of PCT application No. PCT/JP2014/003884, which was filed on Jul. 23, 2014 based on Japanese Patent Application (No. 2013-155811) filed on Jul. 26, 2013, the contents of which are incorporated herein by reference. Also, all the references cited herein are incorporated as a whole.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/JP2014/003884 | Jul 2014 | US |
Child | 14677272 | US |