The present application claims priority under 35 U.S.C. § 119(a) to Korean application number 10-2017-0078414, filed on Jun. 21, 2017, in the Korean Intellectual Property Office, which is incorporated herein by reference in its entirety.
Various embodiments may generally relate to a dementia prediction system and an operating method thereof, and more particularly, to a dementia prediction system using brainwave information and visual information and an operating method thereof.
Technology for measuring various health states of the human body using biometric signals has been studied.
Biometric signals may include, for example, electroencephalogram (EEG brainwave), electromyogram (EMG), electrocardiography (ECG), and the like. EEG among the biometric signals may refer to a waveform obtained by forming an electric field and a magnetic field through flow of current iconized by a stimulus applied to the cerebral cortex between the nerve cells and measuring minute potential change through electronics attached to the scalp. In general, EEG may be distributed in a frequency band of 0 to 100 +Hz and the potential change may be about several tens of μV and thus the amplified potential change may be recorded as EEG.
EEG may be classified into a δ-wave (4 Hz or less), a θ-wave (4-8 Hz), an α-wave (8-12 Hz), a β-wave (12-30 Hz), and a γ-wave (30-50 Hz) according to an activity state of the brain, for example, an oscillating frequency range.
EEG may also be used to diagnose a sleeping and awakened state and an abnormal state of the brain and in recent years, studies on dementia diagnosis through EEG has been actively conducted.
In recent years, with development of digital EEG measuring machines (for example, brainwave sensors), methods for diagnosing dementia through analysis and application of the brainwave have been developed. However, intervention of readers or clinical experts skilled for a reasonable period may be necessary for accurate read-out and subjective judgement criteria even between the skilled clinical experts may be different.
In an embodiment of the present disclosure, a dementia prediction system may be configured to classify a subject into a dementia group and a normal group using parameters obtained by analyzing and operating biometric information and sensory information of the subject.
In an embodiment of the present disclosure, a dementia prediction system may include a bio signal collection module configured to detect brainwave information and visual information of a subject; and a brain aging determination module configured to calculate directional data of a bio signal of the subject from the brainwave information and the visual information and classify the directional data into a dementia group and a normal group using training data of a preset dementia group and a preset normal group.
The brainwave information may include at least one of a peak frequency of an α-wave, an absolute power ratio between the α-wave and a θ-wave, and a power map of the α-wave. The visual information may include information for deposition/non-deposition of an amyloid ingredient in a retina of a human body of the subject.
In an embodiment of the present disclosure, an operating method of a dementia prediction system may include setting training data for a dementia group and a normal group by collecting biometric information and sensory information of the dementia group and collecting biometric information and sensory information of the normal group. The method may further include detecting biometric information and sensory information of a subject. The method may further include generating directional data of sensory information-reflected biometric information by combining and operating the biometric information and the sensory information of the subject. The method may further include classifying the directional data into the dementia group and the normal group using the training data. The method may further include repeatedly performing the detecting of the biometric information and the sensory information of the subject, the generating of the directional data, and the classifying of the directional data.
These and other features, aspects, and embodiments are described below in the section entitled “DETAILED DESCRIPTION”.
The above and other aspects, features and advantages of the subject matter of the present disclosure will be more clearly understood from the following detailed description taken in conjunction with the accompanying drawings, in which:
Various embodiments of the present invention will be described in greater detail with reference to the accompanying drawings. The drawings are schematic illustrations of various embodiments (and intermediate structures). As such, variations from the configurations and shapes of the illustrations as a result, for example, of manufacturing techniques and/or tolerances, are to be expected. Thus, the described embodiments should not be construed as being limited to the particular configurations and shapes illustrated herein but may include deviations in configurations and shapes which do not depart from the spirit and scope of the present invention as defined in the appended claims.
The present invention is described herein with reference to cross-section and/or plan illustrations of idealized embodiments of the present invention. However, embodiments of the present invention should not be construed as limiting the inventive concept. Although a few embodiments of the present invention will be shown and described, it will be appreciated by those of ordinary skill in the art that changes may be made in these embodiments without departing from the principles and spirit of the present invention. In the drawings, the thicknesses of layers and regions may be exaggerated for clarity. In addition, like numbers refer to like element throughout.
Referring to
The bio signal collection module 200 may include a biometric signal detector 210 and a visual information detector 250.
The biometric signal detector 210 may include a brainwave detection sensor configured to detect a brainwave of a subject and a monitoring device. An output of the brainwave detection sensor may be, for example, a result obtained through 8 channels brainwave measuring machines for the brain.
For example, when the 8 channel brainwave measuring machines are used, the brainwave detection sensor may measure the brainwave of the subject using brainwave detection channels (not shown) attached to a brain structure of the subject, for example, 8 portions selected from frontal poles Fp1 and Fp2, frontal portions F3 and F4, central portions C3 and C4, parietal portions P3 and P4, occipital portions O1 and O2, anterior temporal portions F7 and F8, middle temporal portions T3 and T4, posterior temporal portions T5 and T6, a midline-frontal portion Fz, a midline-central portion Cz, and a midline-parietal Pz as illustrated in
For example, the visual information detector 250 may include an eye imaging device. In another example, the visual information detector 250 may include an eye image device and a monitoring device configured to monitor various pieces of visual information detected through the eye imaging device. The visual information detector 250 may detect a check signal of β-amyloid deposited in the retina through the eye image device. The β-amyloid check signal may be a biomarker for sensing a brain aging function.
In general, the deposition of the β-amyloid material in the brain may be a distinct symptom of dementia. However, the deposition of the β-amyloid material in the brain may be founded when the dementia progresses to a certain degree. Since the method for detecting the deposition of the β-amyloid material in the brain is very complicated and is painful to the subject, the early diagnosis of dementia may be difficult.
Unlike other structures of the eye, the retina of the human body is a part of the central nervous system and is known to share many features of the brain. Accordingly, when the retina of the subject is shoot through the eye imaging device, the deposition of the β-amyloid material in the retina may be easily detected. Since the β-amyloid material is earlier deposited in the retina than in the brain, the abnormality of the brain may be early detected. Accordingly, the deposition of the β-amyloid material in the retina may be used as a parameter for determining aging of the brain and thus the test may be easily conducted without pain of the subject unlike a brain test for example, a lumbar puncture or a magnetic resonance imaging (MRI) test.
The eye imaging device may shoot the retina using the hyperspectral imaging technique. The hyperspectral imaging method may be a method for radiating light to the eyeball of the human body, reflecting the light reaching the retina behind the eyeball, and displaying a retinal image using a wavelength of the reflected light.
The deposition/non-deposition of the β-amyloid material in the retina is used as the visual information of the human body in the embodiment, but this is not limited thereto. A variation in widths of the retinal vessels measured through an image sensor of the eye imaging device may be used as the visual information. The widths of the retinal vessels may be values based on the artery and vein of the retina.
The brain aging determination module 300 may include a first analyzer 310, a second analyzer 330, a directional data generator 350, and a machine running unit 370.
The first analyzer 310 may generate parameters necessary for analysis by receiving the brainwave information detected from the biometric signal detector 210.
The first analyzer 310 may generate a waveform feature of the α-wave from various pieces of information input from the biometric signal detector 210. For example, the first analyzer 310 may generate parameters, for example, frequency where the power is maximal within the alpha spectrum (hereinafter, referred to as “peak frequency information of the α-wave”), an absolute power ratio Pα/θ between the α-wave and the θ-wave, a power map of the α-wave, an average, dispersion, kurtosis, and skewness of the α-wave spectrum, and the like.
First, the peak frequency information of the α-wave may be a very important factor in determining a function of the brain. For example, when the brain is aged or is closer to dementia, the peak frequency information of the α-wave may be represented to be slow. When the brain function is good, the peak frequency information of the α-wave may be represented to be high in a range of 8 to 12 Hz which is the frequency range of the α-wave. Since the α-wave is measured in the two portions of the frontal lobe in the steady state as the intrinsic rhythm, signals measured through channels which are in contact with the two portions of the frontal lobe (for example, two portions Fp1 and Fp2 of frontal pole in
It is known that as the aging function of the brain is lowered, the θ-wave power increases in the steady state thus the ratio of θ-wavepower to α-wave power decreases. This is because the cognitive function is injured and thus the brainwave of a low frequency band is relatively activated. Accordingly, the absolute power ratio Pα/θ between the α-wave and the θ-wave in the following Equation 1 may be an another measure in determining the aging of the brain.
P
α/θ=(PαFp1+PαFp2)/(PθFP1+PθFP2) [Equation 1]
Here, PαFp1: absolute power of the α-wave in a first frontal pole,
PαFp2: absolute power of the α-wave in a second frontal pole,
PθFP1: absolute power of the θ-wave in the first frontal pole, and
PθFP2: absolute power of the θ-wave in the second frontal pole.
The absolute power corresponds to the integral of the power spectrum.
The power map of the α-wave may be calculated by an absolute power value of the α-wave obtained through the 8 channels of the biometric signal detector 210. For example, since the absolute power value of the α-wave of a person having a degraded brain function or a dementia patient which is obtained through each channel is reduced as compared with a normal person, such variation amount may be used as a dementia prediction parameter.
As illustrated in
Further, the average and dispersion of the α-wave power spectrum may be used as the dementia prediction parameter. Here, the average of the α-wave power spectrum refers to the average frequency of the power spectrum of the α-wave and the dispersion of the α-wave refers to the extent of spread (for example, a width) of the α-wave power spectrum. For example, as the dementia progresses, average frequency of the α-wave power spectrum may be lowered and the dispersion of the α-wave power spectrum may be widely distributed.
The second analyzer 330 may receive data indicating whether or not the β-amyloid material is deposited in the retina from the eye imaging device of the visual information detector 250 and output a signal of logic 1 in deposition of the β-amyloid material and a signal of logic 0 in non-deposition of the β-amyloid material by processing the received data for the deposition of the β-amyloid material.
The directional data generator 350 may receive the parameter calculated in the first analyzer 310 and the parameter calculated in the second analyzer 330 and form a feature vector for each parameter and a vector space through a combination of the parameters.
For example, the directional data generator 350 may form a feature vector {right arrow over (f)} for dementia prediction and may have following features which are depicted in Equation 2.
{right arrow over (f)}=({circumflex over (α)}, Pα/θ, {right arrow over (Pα)}, avar, Rβ) [Equation 2]
Here, {circumflex over (α)}: peak frequency of the α-wave,
Pα/θ: absolute power ratio between the α-wave and the θ-wave,
{right arrow over (Pα)}: power map of the α-wave,
αmean: average frequency of the α-wave power spectrum,
αavar: dispersion of the α-wave power spectrum, and
Rβ: pigmentation/non-pigmentation of the β-amyloid in the retina.
The directional data generator 350 may induce the feature space (F, F{right arrow over (f)}) for dementia prediction from the following Equation 3.
F=I
{circumflex over (α)}
×I
pα/θ
×I
β
×I
α
×I
α
×I
R
[Equation 3]
Here, I indicates a variable range (distribution) of each parameter and may refer to an effective interval range (effective distribution) for each parameter.
For example, the peak frequency of the α-wave may be corresponding to 8 to 12 Hz and thus the peak frequency distribution I{circumflex over (α)} of the α-wave may be set [8, 12]. The distribution of the absolute power ratio between the α-wave and the θ-wave may be set to [0.1, 3] by considering that the absolute power ratio between the α-wave and the θ-wave in the normal person is 1.5 or more. For example, since the dementia patient has a low peak level of the α-wave and has a high appearance ratio of the θ-wave, the absolute power ratio of the dementia patient may be set between 0.1 and 3. The power map distribution I{circumflex over (P)}
The feature vector of the individual subject may be formed of a space having an effective interval for each parameter, that is, a feature space. The feature space by Equation 3 may be represented with a high-dimensional space, for example, a 13-dimensional space.
The directional data generator 350 may extract the directional data of the biometric information to which the visual information of the subject is reflected by analyzing and operating the biometric signal information and the visual information detected from the subject.
The machine running unit 370 may group similar data entities by receiving the directional data of the biometric information to which the visual information for each subject is reflected. The machine running unit 370 in the embodiment may use supervised learning.
The supervised learning may be a task which derives one function from the training data and may be a series of processes for determining whether the directional data of the subject has a more similar type to which group of data using the training data which is classified into the directional data of the dementia group and the directional data of the normal group on a data base 380.
The training data may be a vector indicating the directivity of the biometric information to which the visual information obtained under the same condition as the directional data measured from the subject is reflected. The training data may be represented to be included in which group through the supervised learning. The term “classification” may refer to a task for representing that the vector of an input entity of the derived function is which type of value. The dementia may be predicted based on the supervised leaning according to the concept of the present disclosure.
The machine running unit 370 in the embodiment may classify the high-dimensional feature space into the normal group and the dementia group and the classification into the normal group and the dementia group will be defined as the following Equation 4.
C={C1,C2} [Equation 4]
Here, C1 refers to the normal group and C2 refers to the dementia group.
When the machine running unit 370 provides ({right arrow over (f)},C) by pairing the normal group of the training data to C1 and pairing the dementia group of the training data to C2, the accurate prediction result may be obtained according to the classified result through a support vector machine (SVM) algorithm.
The SVM algorithm may be typically trained through the supervised learning and may be calculated as the best function which associates the output result with the input data. The SVM algorithm is described in detail in Bernhard Schoelkopf, Alex Smola, “Learning with Kernels” MIT Press, Cambridge, Mass. 2002.
Referring to
The classified group may be formed based on the training data through the machine running process and when target data of an actual subject is input, the target data may be classified to be is included in which group.
For the enhancement and refinement of an analysis signal, it is important to set a gap g between a hyperplane H and support vectors S1 and S2. The hyperplane H may be corresponding to an intermediate value of the gap g between the class C1 and the class C2, the support vector S1 may indicate a boundary line of the class C1 which is closest to the hyperplane H, and the support vector S2 may indicate a boundary line of the class C2 which is closest to the hyperplane H.
As described above, the hyperplane H may be typically set to the intermediate value of the gap g, but the setting of the hyperplane H may be changed by checking the similarity and the degree of correlation within the classified group.
As the number of times that the training data is cumulated is increased, a quality of basic data is increased and thus the machine running unit 370 may accomplish accurate dementia prediction. For example, pieces of directional data of the plurality of subjects may be included in the training data and thus boundary lines of distributions of the classes 1 and 2 C1 and C2 may be more cleared and the probability of error may be further reduced.
For the enhancement and refinement of an analysis signal in the feature vector calculation process, the kurtosis and skewness of the α-wave other than the average Iα
For example, when the waveform of the α-wave has no Gaussian form and is concentrated on an one point or when the waveform of the α-wave has a Gaussian form and the center axis thereof is biased to one side, the kurtosis αkurtosis and skewness αskew of the α-wave power spectrum may be reflected additionally in the feature vector and the vector component of the feature space as in the following Equation 5. The kurtosis and skewness of the α-wave power spectrum may be obtained in the first analyzer 310 through the detection result of the biometric signal detector 210.
{right arrow over (f)}=({circumflex over (α)}, Pα/θ, {right arrow over (Pα)}, αkurtosis, αskew, Rβ)
F=I
{circumflex over (α)}
×I
pα/θ
×I
{right arrow over (P)}
×I
α
×I
α
×I
R
[Equation 5]
The deposition/non-deposition of β-amyloid in the retina is used as the visual information in the embodiment, but a parameter through a pupil light reflection test may be used as the visual information.
Several parameters may be obtained through pupilometry using image analysis technology in a pupillary light reflex (PLR) test such as a standby time to pupil reflection. Accordingly, the subjects may be mapped using feature vectors of Equations 6 and 7.
{right arrow over (f)}=({circumflex over (α)}, Pα/θ, {right arrow over (Pα)}, αkurtosis, αskew, EPLR) [Equation6]
Here, EPLR indicates a variable obtained from the PLR test.
{right arrow over (f)}=({circumflex over (α)}Pα/θ, {right arrow over (Pα)}, αmean, αvar, αkurtosis, αskew, Rβ, EPLR) [Equation 7]
The feature parameters may be added or omitted according to best input for the machine learning algorithm.
An operating method of a dementia prediction system according to an embodiment will be described with reference to
Referring to
The dementia prediction system may detect biometric information and sensory information of a subject under the same condition as the training data (S2).
The dementia prediction system may extract main parameters of the biometric information and the sensory information and generate directional data of the sensory information-reflected biometric information through combination and operation of the main parameters in the brain aging determination module 300 (S3).
The dementia prediction system may classify the directional data into a class 1 (normal group) or a class 2 (dementia group) based on the training data on the database (S4, S5, and S6).
The classification process will be described below in detail. First, the dementia prediction system may determine whether or not the directional data is located within the class 1 C1 (S4). When the directional data of the subject is located within the class 1 C1, the machine running unit 370 may determine the subject to the normal group (S5). For example, the feature vector may be one entity displayed in the feature space and the support vector may correspond to a portion of the feature vectors which determine the hyperplane that separates the two classes.
When the directional data of the subject is not located within the class 1 C10 (S6), as illustrated in
When the directional data of the subject is locate in a gap g2 between the hyperplane H and the class 2 C2 (S63), the machine running unit may determine the subject not to a complete dementia state but to the dementia risk group such as mild cognitive impairment (S64). When the directional data of the subject is not located between the hyperplane H and the feature vector of the class 2 C2 and is corresponding to the class 2 C2, the machine running unit may determine the subject to dementia.
As described above, the error range may be reduced by repeatedly performing operations S2 to S6 and thus more accurate training data may be ensured. For example, the classification-completed directional data of the subject may be continuously reflected to the training data for a reference and thus the more accurate training data for the reference may be formed.
The brainwave detection sensor is used as the bio signal collection module 200 in the embodiment, but this is not limited thereto. The human body signal may be collected through a heartrate variability (HRV) sensing and electromyogram sensor.
According to an embodiment, the dementia group and the normal group may be more accurately classified using the feature vector, which is formed through a combination of the brainwave information as biometric information and the visual information, as the training data.
The subject may be classified into mild cognitive impairment, a dementia potential group (forgetfulness risk group), and a normal group through the feature vector and the feature space obtained from the feature vector, for example, a distance from the hyperplane H and distributions.
The directional data of the subject may be updated as the training data and thus the clearer data reference value may be provided.
The above described embodiments of the present invention are intended to illustrate and not to limit the present invention. Various alternatives and equivalents are possible. The invention is not limited by the embodiments described herein. Nor is the invention limited to any specific type of semiconductor device. Other additions, subtractions, or modifications are obvious in view of the present disclosure and are intended to fall within the scope of the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
10-2017-0078414 | Jun 2017 | KR | national |