Disclosure of Invention
Aiming at the problems, the invention provides a personalized hearing loss modeling method which can ensure that a hearing aid is matched with hearing experts as few as possible and make up the defect that the existing hearing aid technology simply depends on audiogram to classify the hearing loss.
In order to achieve the technical purpose and achieve the technical effect, the invention is realized by the following technical scheme:
a personalized hearing loss modeling method, comprising:
acquiring audiogram of a large number of samples of hearing-impaired patients and corresponding hearing aid insertion gain;
dividing samples of the hearing-impaired patients into three types of moderate hearing loss, severe hearing loss and extremely severe hearing loss according to hearing loss degrees;
step (C), classifying the hearing aid insertion gains of various hearing-impaired patient samples by using an unsupervised clustering method aiming at the classified hearing-impaired patient samples with moderate hearing loss, severe hearing loss and extremely severe hearing loss;
step (D), calculating the average value of audiogram curves corresponding to the insertion gains of the hearing aids in each category to represent the hearing loss of each individual;
and (E) calculating the distance between the audiogram to be classified and each individual hearing loss and classifying according to the minimum distance between the audiogram to be classified and each individual hearing loss.
Preferably, the hearing aid insertion gain G in step (a) has 133 dimensions, respectively, of insertion gain at 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz, 500Hz, 630Hz, 800Hz, 1000Hz, 1250Hz, 1600Hz, 2000Hz, 2500Hz, 3150Hz, 4000Hz, 5000Hz, 6300Hz, and 8000Hz at the input sound pressure level of 50dB, 55dB, 60dB, 65dB, 70dB, 80dB, and 90dB SPL input.
Preferably, the step (B) of classifying the samples of the hearing-impaired patients according to hearing loss degree comprises the following steps:
B1) for each audiogram HLiCalculating the sample g of hearing-impaired patientiAverage hearing loss μ HL ofi:
In the formula, HLi(phi) denotes the sample g of the hearing impaired patientiHearing loss at frequency point phi;
B2) according to hearing impaired patient sample giTo perform a classification of moderate hearing loss, severe hearing loss and very severe hearing loss, wherein: g1For a medium hearing loss set, G2For severe hearing impairment set, G3For the very severe hearing loss set:
G1={gi|μHLi∈(40dB,60dB]}
G2={gi|μHLi∈(60dB,80dB]}
G3={gi|μHLi>80dB}。
preferably, the step (C) specifically includes the steps of:
C1) setting the maximum iteration number as N and the classification number as k, and constructing a sample set G of the insertion gain according to the classification of the step (B), wherein the sample set G is { G ═ G1,g2,g3,…,gnN represents the number of samples;
C2) randomly selecting a point from the input sample set as a first cluster center mu1;
C3) For each point g
iCalculating its minimum distance from the selected cluster center
r represents the number of cluster centers already present, wherein,
denotes g
iAnd cluster center mu
iThe distance of (d);
C4) selection DiThe largest point is taken as the new cluster center mur+1;
C5) Repeating steps C3) and C4) until k cluster centers [ mu ] are selected1,μ2,μ3,…,μk};
C6) Calculate each sample giTo respective cluster center vector muj(j ═ 1,2, …, k) distance dij:
G is prepared from
iIs classified as
ijClass c corresponding to the minimum
j;
C7) For class cjCalculating the mean value of audiogram of all samples as the new cluster center muj;
C8) If all k cluster center vectors mujIf no change occurs or the iteration number N is reached, executing the step C9), otherwise, repeatedly executing the steps C6) to C8);
C9) for each classified subset Gm, the number of samples is
Calculating the contour coefficient of Gm
Wherein s is
iIs the coefficient of the sample;
C10) calculating the contour coefficients under all the k values, and keeping the corresponding classification number k and the corresponding clustering center { mu ] under the condition that the contour coefficients are maximum1,μ2,μ3,…,μk}。
Preferably, 2. ltoreq. k.ltoreq.6.
Preferably, coefficients of samples
Wherein, a
iRepresents g
iAverage distance to homogeneous samples; b
iRepresents g
iAverage distance to all samples in the class closest thereto.
The invention has the beneficial effects that:
compared with the prior art, the invention constructs an individual hearing loss modeling method based on hearing aid gain, classifies hearing-impaired patients with different hearing losses, and selects a representative audiogram to represent the difference of the hearing-impaired patients with different categories so as to overcome the current situation that the existing hearing aid is very dependent on the fitting of hearing experts, wherein:
1) the classification method based on the gain is more in line with the actual compensation effect of the hearing-impaired patients, and can better reflect the individual difference of the hearing-impaired patients;
2) the hearing loss compensation can be better realized by classifying the patients with different hearing losses, so that the hearing aid can be matched with hearing experts as few as possible;
3) the invention can enhance the interpretability of the audiogram, fully explain or optimally utilize a large amount of clinical information contained in the audiogram, and enable non-experts to decide whether to refer a patient to an audiologist, a hearing instrument expert or a physician according to needs.
Detailed Description
The present invention will be better understood and implemented by those skilled in the art by the following detailed description of the technical solution of the present invention with reference to the accompanying drawings and specific examples, which are not intended to limit the present invention.
A personalized hearing loss modeling method, comprising:
and (A) acquiring audiogram of a large number of samples of hearing-impaired patients and corresponding hearing aid insertion gain.
The step is to acquire a sample data source, wherein the more the number of the acquired data sources is, the higher the classification precision of the method is, and the suggested number is not less than (500- & ltSUB & gt 1000- & gt).
Preferably, the hearing aid insertion gain G in step (a) has 133 dimensions, i.e. insertion gain at 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz, 500Hz, 630Hz, 800Hz, 1000Hz, 1250Hz, 1600Hz, 2000Hz, 2500Hz, 3150Hz, 4000Hz, 5000Hz, 6300Hz, and 8000Hz at the respective 19 frequency points at the input sound pressure level of 50dB, 55dB, 60dB, 65dB, 70dB, 80dB, and 90dB SPL, respectively. The hearing aid insertion gain G has 133 dimensions, and includes:
1) 50dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
2) 55dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
3) 60dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
4) 65dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
5) 70dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
6) 80dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level;
7) 90dBSPL (sound pressure level) insertion gain at 19 frequency points below the input sound pressure level.
And (B) dividing the samples of the hearing-impaired patients into three types of moderate hearing loss, severe hearing loss and extremely severe hearing loss according to the hearing loss degree.
Preferably, the step (B) of classifying the sample of the hearing-impaired patient according to hearing loss degree comprises the following steps:
B1) for each audiogram HLiCalculating the sample g of hearing-impaired patientiAverage hearing loss μ HL ofi:
In the formula, HLi(phi) denotes the sample g of the hearing impaired patientiHearing loss at frequency point phi;
B2) according to hearing impaired patient sample giTo perform a classification of moderate hearing loss, severe hearing loss and very severe hearing loss, wherein: g1For a medium hearing loss set, G2For severe hearing impairment set, G3For the very severe hearing loss set:
G1={gi|μHLi∈(40dB,60dB]}
G2={gi|μHLi∈(60dB,80dB]}
G3={gi|μHLi>80dB}。
and (C) classifying the hearing aid insertion gains of the hearing-impaired patient samples of the classified moderate hearing loss, severe hearing loss and extremely severe hearing loss respectively by using an unsupervised clustering method (such as K-Means, mean shift, DBSCAN and the like).
Preferably, the step (C) specifically includes the steps of:
C1) setting the maximum iteration number as N and the classification number as k, and constructing a sample set G of the insertion gain according to the classification of the step (B), wherein the sample set G is { G ═ G1,g2,g3,…,gnN represents the number of samples, wherein preferably, k is more than or equal to 2 and less than or equal to 6, namely, the hearing aid insertion gain classification number of each type of samples of the hearing-impaired patients is between 2 and 6, including 2 and 6, in each classification of moderate hearing loss, severe hearing loss or extremely severe hearing loss.
C2) Randomly selecting a point from the input sample set as a first cluster center mu1;
C3) For each point g
iCalculating its minimum distance from the selected cluster center
r represents the number of cluster centers already present, wherein,
denotes g
iAnd cluster center mu
iThe distance of (d);
C4) selection DiThe largest point is taken as the new cluster center mur+1;
C5) Repeating steps C3) and C4) until k cluster centers [ mu ] are selected1,μ2,μ3,…,μk};
C6) Calculate each sample giTo respective cluster center vector muj(j ═ 1,2, …, k) distance dij:
G is prepared from
iIs classified as
ijClass c corresponding to the minimum
j;
C7) For class cjCalculating the mean value of audiogram of all samples as the new cluster center muj;
C8) If all k cluster center vectors mujIf no change occurs or the iteration number N is reached, executing the step C9), otherwise, repeatedly executing the steps C6) to C8);
C9) for each classified subset Gm, the number of samples is
Calculating the contour coefficient of Gm
Wherein s is
iIs the coefficient of the sample;
wherein coefficients of the samples
Wherein, a
iRepresents g
iAverage distance to homogeneous samples; b
iRepresents g
iAverage distance to all samples in the class closest thereto.
C10) Calculating the contour coefficients under all the k values, and keeping the corresponding classification number k and the corresponding clustering center { mu ] under the condition that the contour coefficients are maximum1,μ2,μ3,…,μk}。
And (D) calculating the average value of audiogram curves corresponding to the insertion gains of the hearing aids in each category to represent the hearing loss of each individual.
And (E) calculating the distance between the audiogram to be classified and each individual hearing loss and classifying according to the minimum distance between the audiogram to be classified and each individual hearing loss.
The following description is provided in connection with specific examples in which experimental data is from the national health and nutrition survey (NHANES), a portion of which assesses the hearing status of a subject by pure tone audiometry, and the NHANES data set contains a number of essentially pure tone audiograms. Experiments searched audiograms obtained from 1999 to 2016, which were obtained using a conventional audiometer with an otoacoustic or plug-in earphone using a standard pure tone audiometric protocol, to obtain a data set of 21436 audiograms for participants aged 12 to 85 years (mean: 39 ± 21 years). The data contains air conduction thresholds at 7 test frequencies: 500Hz, 1000Hz, 2000Hz, 3000Hz, 4000Hz, 6000Hz and 8000 Hz.
For efficient evaluation, we preprocessed the data:
1) deleting the incomplete audiogram with at least one missing threshold;
2) discard data below moderate hearing loss, i.e., audiogram with average thresholds below 40dB HL at 500Hz, 1000Hz, 2000Hz, and 4000 Hz.
The screened data comprises 1578 moderate hearing loss samples, 356 severe hearing loss samples and 63 extremely severe hearing loss samples. Wherein: the medium hearing loss samples were classified into 3 classes, each consisting of 494, 393, and 691 samples; severe hearing loss samples were classified into 3 classes, each consisting of 135, 122, and 99 samples; the very severe hearing loss samples were classified into 2 classes, each consisting of 22 and 41 samples, for a total of 8 classes. The contour coefficients of the classification results of moderate, severe and extremely severe hearing loss samples can reach 0.345, 0.380 and 0.466 respectively. The results of the classification are shown in fig. 1 to 3:
as can be seen from FIG. 1, the discrimination of different categories in the high-frequency part is higher, and the characteristics of high-frequency weight loss of hearing-impaired patients are met. The whole hearing-impaired people are divided into 8 classes, which are lower than the traditional classification method, and the realization possibility is provided for the class switching of the non-fitting hearing aid. Moreover, the time for the expert to check is far longer than that of the method of the invention, and the checking time of the invention is almost zero.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.