Schizophrenia classification method and system based on brain wave signal acquisition

Abstract

The invention discloses a schizophrenia classification method and system based on brain wave signal collection, and relates to the technical field of brain wave data processing, and the main scheme is as follows: collecting frequency data and time domain data of brain wave signals; the method comprises the following steps: acquiring environment data and blood oxygen data of an environment when brain wave signals are acquired, calculating a corresponding brain frequency band stability difference index, a brain state stability index, an electroencephalogram signal acquisition comprehensive quality index and a brain health comprehensive evaluation index according to the data, and further calculating a brain function overall health comprehensive index; and presetting a brain function overall state threshold value set, comparing the brain function overall health comprehensive index with a threshold value, and classifying schizophrenia according to the result.

Description

Technical Field

[0001] The present invention relates to the technical field of electroencephalogram data processing, and particularly to a schizophrenia classification method and system based on electroencephalogram signal acquisition. Background Art

[0002] Schizophrenia is a complex neuropsychiatric disorder with diverse clinical manifestations, including symptoms such as hallucinations, delusions, thought disorders, and emotional apathy. In recent years, with the development of electroencephalogram technology, objective diagnostic methods based on electroencephalogram signals have gradually received attention. Electroencephalogram signals can reflect the electrical activities of the brain, have high temporal resolution, and can capture the subtle differences between schizophrenia patients and healthy people.

[0003] Traditional diagnostic methods mainly rely on clinical symptom observation and patient self-report, which have problems such as strong subjectivity and lack of objective indicators. The collected data is relatively single, and there are obvious deficiencies in terms of accuracy and early diagnosis ability as a whole. Summary of the Invention

[0004] (1) Technical Problems to be Solved Aiming at the deficiencies of the prior art, the present invention provides a schizophrenia classification method and system based on electroencephalogram signal acquisition. By collecting the frequency data and time-domain data of electroencephalogram signals; obtaining the environmental data and blood oxygen data of the environment where the electroencephalogram signals are collected, calculating corresponding indicators based on these data, and further calculating the overall brain function health comprehensive index, presetting the overall brain function state threshold set, comparing the overall brain function health comprehensive index with the threshold, and classifying schizophrenia accordingly, which solves the problems of relatively single collected data, strong subjectivity, lack of objective indicators, and obvious deficiencies in terms of accuracy and early diagnosis ability as a whole.

[0005] (2) Technical Solutions To achieve the above object, the present invention is realized through the following technical solutions: A schizophrenia classification method based on electroencephalogram signal acquisition, including the following steps: Step 1: Collect the frequency data and time-domain data of electroencephalogram signals; obtain the environmental data and blood oxygen data of the environment where the electroencephalogram signals are collected; Step 2: Calculate the relative power ratio of electroencephalogram frequency bands , the standard deviation of electroencephalogram frequency band power and the power asymmetry index of electroencephalogram frequency bands ; Calculate the difference index of brain frequency band stability according to the relative power ratio of electroencephalogram frequency bands , the standard deviation of electroencephalogram frequency band power and the power asymmetry index of electroencephalogram frequency bands ; ; Calculate the electroencephalogram (EEG) pulse frequency modulation index based on time-domain data and the coefficient of variation of EEG signal activity ; calculate the brain state stability index based on the EEG pulse frequency modulation index and the coefficient of variation of EEG signal activity ; Calculate the environmental interference index based on environmental data , calculate the signal acquisition comfort adaptation index based on the environmental interference index , and further calculate the comprehensive quality index of EEG signal acquisition ; Calculate the cerebral blood oxygen supply-demand balance index based on blood oxygen data , calculate the cerebral function reserve index based on the cerebral blood oxygen supply-demand balance index , and further calculate the comprehensive evaluation index of brain health ; Calculate the comprehensive index of overall brain function health based on the brain frequency band stability difference index , the brain state stability index , the comprehensive quality index of EEG signal acquisition and the comprehensive evaluation index of brain health ; Step 3: Preset the set of thresholds for the overall state of brain function, compare the comprehensive index of overall brain function health with the set of thresholds for the overall state of brain function, and classify schizophrenia according to the comparison result

[0006] In a preferred embodiment of the above method for classifying schizophrenia based on EEG signal acquisition: The method for calculating the relative power ratio of EEG frequency bands , the standard deviation of EEG frequency band power and the power asymmetry index of EEG frequency bands is as follows The frequency data includes the real-time power of frequency bands ; Calculate the relative power ratio of EEG frequency bands based on the real-time power of frequency bands , and the formula is as follows

[0007] where is the real-time power of the th frequency band , is the serial number corresponding to different frequency bands, and the value range is ​​​​is the number of frequency bands, and the value is a positive integer; The frequency data also includes the average value of the frequency band power ; According to the real-time power of the frequency band and the average value of the frequency band power calculate the standard deviation of the electroencephalogram frequency band power , and the formula is as follows:

[0008] where, is the frequency band power of the th sampling of the th frequency band, is the serial number corresponding to different sampling times, and the value is ; is the number of sampling times, and the value is a positive integer; The frequency data also includes the frequency band power of the right brain region and the frequency band power of the left brain region ; According to the frequency band power of the right brain region and the frequency band power of the left brain region calculate the asymmetry index of the electroencephalogram frequency band power , and the formula is as follows:

[0009] where, is the power of the th frequency band of the right brain region, is the power of the th frequency band of the left brain region.

[0010] In the above preferred scheme of a schizophrenia classification method based on electroencephalogram signal acquisition: The method for calculating the brain frequency band stability difference index is as follows: According to the relative power ratio of the electroencephalogram frequency band , the standard deviation of the electroencephalogram frequency band power and the asymmetry index of the electroencephalogram frequency band power calculate the brain frequency band stability difference index , and the formula is as follows:

[0011] In the above preferred scheme of a schizophrenia classification method based on electroencephalogram signal acquisition: The method for calculating the electroencephalogram pulse frequency modulation index and the coefficient of variation of the electroencephalogram signal activity is as follows: The time-domain data includes the instantaneous frequency of the electroencephalogram pulse and the average instantaneous frequency ; According to the instantaneous frequency of the electroencephalogram (EEG) pulse and the average instantaneous frequency calculate the EEG pulse frequency modulation index , and the formula is as follows:

[0012] where, is the instantaneous frequency of the th EEG signal pulse, is the influence coefficient corresponding to the th EEG signal pulse, is the serial number corresponding to different EEG signals, and the value range is ; is the number of EEG signals, and the value range is positive integers; The time-domain data also includes the EEG wave amplitude and the average value of the EEG wave amplitude ; According to the EEG wave amplitude and the average value of the EEG wave amplitude calculate the coefficient of variation of EEG signal activity , and the formula is as follows:

[0013] where, is the EEG wave amplitude at the th time point, is the serial number corresponding to different time points, and the value range is ; is the number of all time points, and the value range is positive integers.

[0014] In a preferred scheme of the above schizophrenia classification method based on EEG signal acquisition: The method for calculating the brain state stability index is as follows: According to the EEG pulse frequency modulation index and the coefficient of variation of EEG signal activity calculate the brain state stability index , and the formula is as follows:

[0015] In a preferred scheme of the above schizophrenia classification method based on EEG signal acquisition: The method for calculating the comprehensive quality index of EEG signal acquisition is as follows: The environmental data includes the actual temperature value , the rated temperature value 、Actual humidity value 、Rated humidity value and electromagnetic radiation intensity value ; According to the actual temperature value 、Rated temperature value 、Actual humidity value 、Rated humidity value and electromagnetic radiation intensity value calculate the environmental interference index , and the formula is as follows:

[0016] Among them, is the weight coefficient of the temperature deviation value, and the value range is 0.1~0.3; is the weight coefficient of the humidity deviation value, and the value range is 0.2~0.3; is the weight coefficient of the electromagnetic radiation intensity value, and the value range is 0.5~0.6; and ; is the normalization constant; According to the environmental interference index calculate the signal acquisition comfort adaptation index , and the formula is as follows:

[0017] According to the environmental interference index and the signal acquisition comfort adaptation index calculate the comprehensive quality index of EEG signal acquisition , and the formula is as follows:

[0018] In the preferred scheme of the above schizophrenia classification method based on EEG signal acquisition: The method for calculating the comprehensive brain health assessment index is: Blood oxygen data includes arterial oxygen saturation value 、Hemoglobin value 、Venous oxygen saturation value and cardiac output value ; According to the arterial oxygen saturation value 、Hemoglobin value 、Venous oxygen saturation value and cardiac output value calculate the brain blood oxygen supply-demand balance index , and the formula is as follows:

[0019] According to the cerebral blood oxygen supply-demand balance index Calculate the cerebral function reserve index , and the formula is as follows:

[0020] Wherein, is the influence coefficient of the cerebral blood oxygen supply-demand balance index , is the normalization constant, and the value is a positive integer; According to the cerebral blood oxygen supply-demand balance index and the cerebral function reserve index Calculate the comprehensive cerebral health assessment index , and the formula is as follows:

[0021] In a preferred embodiment of the above-mentioned schizophrenia classification method based on electroencephalogram signal acquisition: The method for calculating the comprehensive cerebral function overall health index is as follows: According to the cerebral frequency band stability difference index , the cerebral state stability index , the comprehensive electroencephalogram signal acquisition quality index and the comprehensive cerebral health assessment index Calculate the comprehensive cerebral function overall health index , and the formula is as follows:

[0022] In a preferred embodiment of the above-mentioned schizophrenia classification method based on electroencephalogram signal acquisition: The criteria for classifying schizophrenia are: The preset cerebral function overall state threshold set includes a risk-free threshold , a low-risk threshold and a severe-risk threshold ; Compare the comprehensive cerebral function overall health index with the cerebral function overall state threshold set, and classify schizophrenia according to the comparison result. The criteria are as follows:

[0023] The present invention also discloses a schizophrenia classification system based on electroencephalogram signal acquisition, which is characterized in that it includes: A data acquisition module, configured to collect frequency data and time-domain data of electroencephalogram signals; obtain environmental data and blood oxygen data of the environment where the electroencephalogram signals are collected; An index calculation module, configured to calculate the relative power ratio of electroencephalogram frequency bands according to the frequency data and the standard deviation of the brain wave frequency band power and the asymmetry index of the brain wave frequency band power ; calculating the difference index of the brain frequency band stability according to the relative power ratio of the brain wave frequency band and the standard deviation of the brain wave frequency band power and the asymmetry index of the brain wave frequency band power ; calculating the difference index of the brain frequency band stability ; Calculating the electroencephalogram pulse frequency modulation index and the coefficient of variation of the electroencephalogram signal activity according to the time domain data ; calculating the brain state stability index according to the electroencephalogram pulse frequency modulation index and the coefficient of variation of the electroencephalogram signal activity ; calculating the brain state stability index ; ; Calculating the environmental interference index according to the environmental data , calculating the signal acquisition comfort adaptation index according to the environmental interference index , and further calculating the comprehensive quality index of the electroencephalogram signal acquisition ; ; Calculating the brain blood oxygen supply-demand balance index according to the blood oxygen data , calculating the brain function reserve index according to the brain blood oxygen supply-demand balance index , and further calculating the comprehensive evaluation index of brain health ; ; Calculating the comprehensive index of the overall brain function health according to the difference index of the brain frequency band stability , the brain state stability index , the comprehensive quality index of the electroencephalogram signal acquisition and the comprehensive evaluation index of brain health ; ; A classification module, configured to preset a set of brain function overall state thresholds, compare the comprehensive index of the overall brain function health with the set of brain function overall state thresholds, and classify schizophrenia according to the comparison result

[0024] (III) Beneficial effects The present invention provides a schizophrenia classification method and system based on electroencephalogram signal acquisition, having the following beneficial effects (1) By collecting the frequency data and time domain data of the electroencephalogram signal; obtaining the environmental data and blood oxygen data of the environment where the electroencephalogram signal is collected, it is possible to comprehensively obtain information related to brain activities, realizing multi-dimensional brain function evaluation, and providing a rich data basis for subsequent analysis and evaluation

[0025] (2) By processing the frequency data, time-domain data, environmental data, and blood oxygen data of the collected electroencephalogram (EEG) signals, the brain frequency band stability difference index, brain state stability index, comprehensive quality index of EEG signal collection, and comprehensive brain health assessment index are calculated respectively. Then, the overall comprehensive brain function health index is obtained. These indexes provide a scientific basis for comprehensively and objectively evaluating the brain function state and avoid the limitations of a single index.

[0026] (3) By presetting a set of thresholds for the overall brain function state and comparing the overall comprehensive brain function health index with them, based on objective quantitative indexes, the uncertainty of subjective judgment is avoided. This provides a method based on objective data for the early diagnosis of schizophrenia, helps improve the accuracy of diagnosis. In this way, schizophrenia patients can be more accurately identified and classified, which is helpful for early diagnosis, treatment, and disease monitoring. Description of the Drawings

[0027] Figure 1 It is a schematic diagram of the working steps of a schizophrenia classification method based on EEG signal collection according to the present invention. Detailed Embodiment

[0028] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0029] Please refer to Figure 1 , the present invention provides a schizophrenia classification method based on EEG signal collection, including the following steps: Step 1: Collect the frequency data and time-domain data of the EEG signals; obtain the environmental data and blood oxygen data of the environment where the EEG signals are collected.

[0030] When in use, in combination with the content of Step 1: By collecting the above relevant data, information related to brain activities can be comprehensively obtained, realizing multi-dimensional brain function evaluation, and providing a rich data basis for subsequent analysis and evaluation.

[0031] Step 2: Calculate the relative power ratio of EEG frequency bands , the standard deviation of EEG frequency band power and the power asymmetry index of EEG frequency bands ; according to the relative power ratio of EEG frequency bands , the standard deviation of EEG frequency band power and the power asymmetry index of EEG frequency bands Calculate the difference index of brain frequency band stability .

[0032] Calculate the electroencephalogram (EEG) pulse frequency modulation index and the coefficient of variation of EEG signal activity according to the time-domain data ; Calculate the brain state stability index according to the EEG pulse frequency modulation index and the coefficient of variation of EEG signal activity ; Calculate the brain state stability index .

[0033] Calculate the environmental interference index according to the environmental data , calculate the signal acquisition comfort adaptation index according to the environmental interference index , and further calculate the comprehensive quality index of EEG signal acquisition . .

[0034] Calculate the brain blood oxygen supply-demand balance index according to the blood oxygen data , calculate the brain function reserve index according to the brain blood oxygen supply-demand balance index , and further calculate the comprehensive brain health assessment index . .

[0035] Calculate the comprehensive index of overall brain function health according to the difference index of brain frequency band stability , the brain state stability index , the comprehensive quality index of EEG signal acquisition and the comprehensive brain health assessment index ; .

[0036] Step 201: The method for calculating the relative power ratio of brain wave frequency bands , the standard deviation of brain wave frequency band power and the power asymmetry index of brain wave frequency bands is as follows: The frequency data includes the real-time power of frequency bands .

[0037] It should be noted that the real-time power of frequency bands represents the real-time power of the th frequency band. The brain wave signal is collected by an electroencephalogram device, and the brain wave signal is decomposed into and other n different frequency bands. The power of the signal is calculated within each frequency band. For example, the Fourier transform is performed on the brain wave signal to obtain the energy of the signal, and then the energy is divided by the number of sample points within the frequency band to obtain the real-time power of the frequency band .

[0038] According to the real-time power of frequency bands Calculating the relative power ratio of brain wave frequency bands , based on the following formula:

[0039] Wherein, is the real-time power of the th frequency band, is the serial number corresponding to different frequency bands, and the value range is ; is the number of frequency bands, and the value range is a positive integer.

[0040] It should be noted that in this formula: is the sum of the real-time powers of all frequency bands, representing the total energy of the entire brain wave, represents the energy magnitude of the brain wave in this frequency band, calculates the proportion of the power of this frequency band in the total power. By summing up this proportion for all frequency bands, the relative power ratio of brain wave frequency bands is obtained. The relative power ratio of brain wave frequency bands can reflect the relative distribution of brain wave powers in different frequency bands.

[0041] The frequency data also includes the average value of the frequency band power .

[0042] It should be noted that the average value of the frequency band power represents the average value of the real-time powers of the frequency bands sampled each time for all frequency bands .

[0043] Based on the real-time power of the frequency band and the average value of the frequency band power calculate the standard deviation of the brain wave frequency band power , based on the following formula:

[0044] Wherein, is the power of the th frequency band at the th sampling, is the serial number corresponding to different sampling times, and the value range is ; is the number of sampling times, and the value range is a positive integer.

[0045] It should be noted that in this formula: The numerator calculates the square of the deviation between the power of the frequency band at each sampling point and the average value of the frequency band power. By accumulating the squared deviations of all sampling points, the total sum of squared deviations is obtained, Denote the total number of data points. By dividing the total sum of squared deviations by the number of data points, the average degree of deviation is obtained, and then taking the square root, which is the power standard deviation.

[0046] The frequency data also includes the power of the brain wave bands in the right brain region and the power of the brain wave bands in the left brain region .

[0047] It should be noted that the power of the brain wave bands in the right brain region represents the power of the th brain wave band in the right brain region. By collecting the electroencephalogram signals through an electroencephalogram device, the electroencephalogram signals are decomposed into and other n different brain wave bands, and the signal power of the right brain region is calculated respectively using Fourier transform within each brain wave band; the power of the brain wave bands in the left brain region represents the power of the th brain wave band in the left brain region. By collecting the electroencephalogram signals through an electroencephalogram device, the electroencephalogram signals are decomposed into and other n different brain wave bands, and the signal power of the left brain region is calculated respectively using Fourier transform within each brain wave band.

[0048] According to the power of the brain wave bands in the right brain region and the power of the brain wave bands in the left brain region calculate the brain wave band power asymmetry index , and the formula is as follows:

[0049] where is the power of the th brain wave band in the right brain region, is the power of the th brain wave band in the left brain region.

[0050] It should be noted that in this formula: calculates the relative difference in power between the left and right brain regions for each brain wave band. By summing up such relative differences for all brain wave bands, the total power asymmetry degree is obtained. is a normalization factor, aiming to map the calculation result to a specific range so that the final asymmetry index has a more reasonable numerical range and physical meaning.

[0051] Step 202: Calculate the brain wave band stability difference index The method is as follows: According to the relative power ratio of the brain wave bands , the power standard deviation of the brain wave bands and the brain wave band power asymmetry index calculate the brain wave band stability difference index , and the formula is as follows:

[0052] It should be noted that in this formula: the numerator reflects the stability of power. When is smaller, the value is larger. It reflects the symmetry of the power between the left and right brain regions. When is smaller, the value is larger; Multiply the three parts and divide by for normalization to obtain the brain frequency band stability difference index .

[0053] Step 203: Calculate the electroencephalogram (EEG) pulse frequency modulation index and the coefficient of variation of EEG signal activity The method is as follows: The time-domain data includes the instantaneous frequency of the EEG pulse and the average instantaneous frequency .

[0054] It should be noted that the instantaneous frequency of the EEG pulse represents the instantaneous frequency of the th EEG signal pulse. The EEG wave signal is collected by an EEG device, and the obtained signal is decomposed in the time and frequency domains to obtain the frequency value of each pulse at a specific moment, that is, the instantaneous frequency; The average instantaneous frequency represents the average instantaneous frequency of all EEG signal pulses and is obtained by averaging the instantaneous frequencies.

[0055] Calculate the EEG pulse frequency modulation index and the average instantaneous frequency according to the instantaneous frequency of the EEG pulse. The formula is as follows:

[0056] Among them, is the instantaneous frequency of the th EEG signal pulse, is the influence coefficient corresponding to the th EEG signal pulse. By analyzing the EEG signal data, the least squares method is used to fit the relationship between the pulse frequency and the modulation index, and then the influence coefficient is determined. is the serial number corresponding to different EEG signals, and the value is ; is the number of EEG signals, and the value is a positive integer.

[0057] It should be noted that in this formula: The sum of the deviations between the instantaneous frequencies of all EEG signal pulses and the average instantaneous frequency is calculated, multiplied by the corresponding influence coefficient, and divided by to obtain the EEG pulse frequency modulation index , The larger the value, the greater the fluctuation of the EEG pulse frequency relative to the average instantaneous frequency, that is, the more significant the frequency modulation; The smaller the value, the relatively more stable the EEG pulse frequency and the lower the degree of frequency modulation.

[0058] The time-domain data also includes the EEG wave amplitude and the average value of the EEG wave amplitude .

[0059] It should be noted that the EEG wave amplitude represents the EEG wave amplitude at the th time point, which is obtained by collecting EEG signals through an electroencephalogram device and then measuring the amplitudes of EEG waves at different time points; the average value of the EEG wave amplitude represents the average value of the amplitudes within all time points, which is obtained by averaging the amplitudes of EEG waves at all time points.

[0060] According to the EEG wave amplitude and the average value of the EEG wave amplitude calculate the coefficient of variation of EEG signal activity , and the formula is as follows:

[0061] Among them, is the EEG wave amplitude at the th time point, is the serial number corresponding to different time points, and the value range is ; is the number of all time points, and the value range is a positive integer.

[0062] It should be noted that in this formula: is to normalize the total sum of squared deviations, is the degree of freedom, which is used to correct the deviation of the sample variance. The denominator is the average value of the EEG wave amplitude. By dividing the standard deviation by the average value, the coefficient of variation is obtained. The coefficient of variation can eliminate the influence of dimensions and facilitate comparing the dispersion degrees of different data sets.

[0063] Step 204: Calculate the brain state stability index The method is as follows: According to the EEG pulse frequency modulation index and the coefficient of variation of EEG signal activity calculate the brain state stability index , the formula is as follows:

[0064] It should be noted that in this formula: the numerator is 1, representing the reference value, and the denominator comprehensively considers the numerical values of electroencephalogram (EEG) pulse frequency modulation and EEG signal activity variation. To avoid the situation where the denominator is 0, and at the same time, a certain offset is made to the calculation result. The two are divided to obtain the brain state stability index . The smaller the index value, the smaller the comprehensive influence of EEG pulse frequency modulation and EEG signal activity variation, and the more stable the brain state; the larger the index value, the more unstable the brain state.

[0065] Step 205: Calculate the comprehensive quality index of EEG signal acquisition The method is as follows: The environmental data includes the actual temperature value , the rated temperature value , the actual humidity value , the rated humidity value and the electromagnetic radiation intensity value .

[0066] It should be noted that the actual temperature value represents the actual temperature measurement value in the environment where the EEG signal is collected, which is directly measured in the environment by a thermometer; the rated temperature value represents the ideal temperature value specified by the standard for the environment where the EEG signal is collected, which is obtained from the design requirements; the actual humidity value represents the actual humidity measurement value in the environment where the EEG signal is collected, which is directly measured in the environment by a hygrometer; the rated humidity value represents the ideal humidity value specified by the standard for the environment where the EEG signal is collected, which is obtained from the design requirements; the electromagnetic radiation intensity value represents the electromagnetic radiation intensity measurement value in the environment where the EEG signal is collected, which is measured in the environment by an electromagnetic radiation measuring instrument.

[0067] According to the actual temperature value , the rated temperature value , the actual humidity value , the rated humidity value and the electromagnetic radiation intensity value calculate the environmental interference index , and the formula is as follows:

[0068] Among them, is the weight coefficient of the temperature deviation value. According to the temperature deviation value, the environmental interference index The influence degree is determined, and the value ranges from 0.1 to 0.3; is the weight coefficient of the humidity deviation value. According to the influence degree of the humidity deviation value on the environmental interference index it is determined, and the value ranges from 0.2 to 0.3; is the weight coefficient of the electromagnetic radiation intensity value. According to the influence degree of the electromagnetic radiation intensity value on the environmental interference index it is determined, and the value ranges from 0.5 to 0.6; and ; is the normalization constant, which is used to normalize the calculation result to ensure that the environmental interference index is within a reasonable range.

[0069] It should be noted that in this formula: represents the deviation between the actual temperature and the rated temperature, and the influence of the temperature deviation on the environmental interference index is reflected by multiplying the temperature deviation by its weight coefficient. represents the deviation between the actual humidity and the rated humidity, and the influence of the humidity deviation on the environmental interference index is reflected by multiplying the humidity deviation by its weight coefficient. The influence of electromagnetic radiation on the environmental interference index is reflected by multiplying the electromagnetic radiation intensity by its weight coefficient. The sum of the three parts is divided by to obtain the environmental interference index .

[0070] According to the environmental interference index calculate the signal acquisition comfort adaptation index , and the formula is as follows:

[0071] It should be noted that in this formula: By dividing the numerator 100 by the denominator , the signal acquisition comfort adaptation index is obtained. When the value of the environmental interference index is larger, the value of the denominator is larger, resulting in the value being smaller, which indicates that the environmental interference is larger and the comfort of signal acquisition is lower; when the value of the environmental interference index is smaller, the value of the denominator is smaller, resulting in the value being larger, which indicates that the environmental interference is smaller and the comfort of signal acquisition is higher.

[0072] According to the environmental interference index and the signal acquisition comfort adaptation index calculate the comprehensive quality index of EEG signal acquisition , and the formula is as follows:

[0073] It should be noted that in this formula: by dividing the signal acquisition comfort adaptation index by , the comprehensive quality index of EEG signal acquisition is obtained. Adding 1 to the environmental interference index is to avoid the situation of the denominator being 0 and ensure the stability of the calculation; the larger the value of the index, the higher the comprehensive quality of EEG signal acquisition in the current environment; the smaller the value of the index, the lower the acquisition quality.

[0074] Step 206: Calculate the comprehensive brain health assessment index The method is as follows: The blood oxygen data includes the arterial oxygen saturation value , hemoglobin value , venous oxygen saturation value and cardiac output value .

[0075] It should be noted that the oxygen saturation value represents the saturation of oxygen in arterial blood, usually expressed as a percentage, and is obtained by measuring the proportion of oxyhemoglobin in total hemoglobin in arterial blood with a pulse oximeter; the hemoglobin value represents the content of hemoglobin in the blood, usually expressed in grams per liter, and is obtained by measuring the hemoglobin content in a blood sample using a hemoglobin analyzer; the venous oxygen saturation value represents the saturation of oxygen in venous blood, usually expressed as a percentage, and is obtained by measuring with a venous blood oxygen measurement device; the cardiac output value represents the amount of blood pumped out by the heart per minute, usually expressed in liters per minute, and is directly measured by echocardiogram.

[0076] According to the arterial oxygen saturation value , hemoglobin value , venous oxygen saturation value and cardiac output value calculate the brain blood oxygen supply-demand balance index , and the formula is as follows:

[0077] It should be noted that in this formula: The difference reflects the oxygen consumption situation during blood circulation, represents the maximum amount of oxygen that the blood can carry under the current hemoglobin content and cardiac output, is a constant related to the ability of hemoglobin to bind oxygen. Usually, each gram of hemoglobin can bind at most milliliters of oxygen; by dividing the oxygen content difference between arterial blood and venous blood by the maximum amount of oxygen that the blood can carry , the cerebral blood oxygen supply-demand balance index is obtained , and this index can reflect the balance between the supply and demand of oxygen in the brain during blood circulation.

[0078] According to the cerebral blood oxygen supply-demand balance index the cerebral function reserve index is calculated , and the formula is as follows:

[0079] where is the influence coefficient of the cerebral blood oxygen supply-demand balance index . By collecting data on multiple groups of cerebral blood oxygen supply-demand balance indexes and the corresponding cerebral function reserve conditions, using the relationship model between GOI and cerebral function reserve, applying interference, collecting data, and training the model with the data, the influence coefficient is obtained. is the normalization constant, which takes a positive integer value and is used to normalize the calculation result to ensure that the cerebral function reserve index is within a reasonable range.

[0080] It should be noted that in this formula: by dividing the numerator 1 by the denominator , a fractional value is obtained, and then subtracting this fractional value from 1 to obtain the cerebral function reserve index ; when the value of the cerebral blood oxygen supply-demand balance index is larger, the value of the denominator is smaller, resulting in having a larger value, which indicates better cerebral function reserve; when the value of the cerebral blood oxygen supply-demand balance index is smaller, the value of the denominator is larger, resulting in having a smaller value, which indicates poorer cerebral function reserve.

[0081] According to the cerebral blood oxygen supply-demand balance index and the cerebral function reserve index the comprehensive cerebral health assessment index is calculated , and the formula is as follows:

[0082] It should be noted that in this formula: by multiplying the cerebral blood oxygen supply-demand balance index and the cerebral function reserve index , the comprehensive cerebral health assessment index is obtained. When and are both larger, has a larger value, which indicates that the brain is in a better state in terms of both blood oxygen supply-demand balance and function reserve, and the overall cerebral health condition is good; when or One of them is smaller, and the value is smaller, indicating that there are problems in the balance between blood oxygen supply and demand or functional reserve in the brain, and the overall brain health condition is poor.

[0083] Step 207: Calculate the comprehensive index of overall brain function health The method is as follows: According to the brain frequency band stability difference index , the brain state stability index , the comprehensive quality index of electroencephalogram signal acquisition and the comprehensive evaluation index of brain health calculate the comprehensive index of overall brain function health , and the formula is as follows:

[0084] It should be noted that in this formula: by multiplying these four indexes and then taking the fourth root, the comprehensive index of overall brain function health is obtained. When the value is larger, it indicates that the brain is in a better state in terms of frequency band stability, state stability, signal acquisition quality and overall health condition, and the overall brain function health condition is good; when the value is smaller, it indicates that there are problems in the corresponding aspects, and the overall brain function health condition is poor.

[0085] When in use, combine the content of steps 201 to 207: By processing the frequency data, time domain data, environmental data and blood oxygen data of the collected electroencephalogram signals, calculate the brain frequency band stability difference index, the brain state stability index, the comprehensive quality index of electroencephalogram signal acquisition and the comprehensive evaluation index of brain health respectively, and then obtain the comprehensive index of overall brain function health. These indexes provide a scientific basis for comprehensively and objectively evaluating the brain function state and avoid the limitations of single indexes.

[0086] Step 3: Preset the set of overall brain function state thresholds, and compare the comprehensive index of overall brain function health with the set of overall brain function state thresholds, and classify schizophrenia according to the comparison results.

[0087] Step 301: The method for presetting the set of overall brain function state thresholds is as follows: Collect historical data of electroencephalogram signals of a large number of schizophrenia patients, and calculate the comprehensive index of overall brain function health of different patients through the historical data, screen out the lower 20% of the comprehensive index of overall brain function health and take its average value as the risk-free threshold ; Screen out the top 50% of the overall brain function health composite indicators, calculate their mean and standard deviation, and subtract twice the standard deviation from the mean as the low-risk threshold ; Add twice the standard deviation to the mean as the severe-risk threshold ; Among them, .

[0088] Step 302: The criteria for classifying schizophrenia are as follows: The set of overall brain function state thresholds includes the risk-free threshold , the low-risk threshold and the severe-risk threshold ; Compare the overall brain function health composite indicator with the set of overall brain function state thresholds, and classify schizophrenia according to the comparison results. The criteria are as follows:

[0089] When in use, combine the content of Steps 301 to 302: By presetting the set of overall brain function state thresholds and comparing the overall brain function health composite indicator with it, based on objective quantitative indicators, the uncertainty of subjective judgment is avoided, the accuracy of diagnosis is improved, and it is helpful for early diagnosis, treatment, and condition monitoring.

[0090] On the other hand, the present invention also discloses a schizophrenia classification system based on electroencephalogram signal acquisition, including: A data acquisition module for acquiring the frequency data and time-domain data of electroencephalogram signals; obtaining the environmental data and blood oxygen data of the environment where the electroencephalogram signals are acquired; An index calculation module for calculating the relative power ratio of electroencephalogram frequency bands , the standard deviation of electroencephalogram frequency band power and the power asymmetry index of electroencephalogram frequency bands ; calculating the difference index of brain frequency band stability according to the relative power ratio of electroencephalogram frequency bands , the standard deviation of electroencephalogram frequency band power and the power asymmetry index of electroencephalogram frequency bands ; ; Calculating the electroencephalogram pulse frequency modulation index and the coefficient of variation of electroencephalogram signal activity according to the time-domain data ; calculating the brain state stability index and the coefficient of variation of electroencephalogram signal activity according to the electroencephalogram pulse frequency modulation index ; ; Calculating the environmental interference index according to the environmental data , calculate the signal acquisition comfort adaptation index according to the environmental interference index , and further calculate the comprehensive quality index of EEG signal acquisition ; ; Calculate the cerebral blood oxygen supply-demand balance index according to the blood oxygen data , calculate the cerebral functional reserve index according to the cerebral blood oxygen supply-demand balance index , and further calculate the comprehensive evaluation index of brain health ; ; According to the cerebral frequency band stability difference index , cerebral state stability index , the comprehensive quality index of EEG signal acquisition and the comprehensive evaluation index of brain health calculate the comprehensive index of overall brain function health ; A classification module, configured to preset a set of thresholds for the overall state of brain function, and compare the comprehensive index of overall brain function health with the set of thresholds for the overall state of brain function, and classify schizophrenia according to the comparison result.

[0091] The above embodiments can be implemented in whole or in part by software, hardware, firmware, or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. Those of ordinary skill in the art will realize that the units and algorithm steps of the examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware or a combination of computer software and electronic hardware. Whether these functions are executed in hardware or software depends on the specific application and design constraints of the technical solution.

[0092] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0093] The above is only the specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed in the present application, and all should be covered by the protection scope of the present application.

Claims

1. A method for classifying schizophrenia based on electroencephalogram signal acquisition, characterized in that: Including the following steps: Step 1: Collect the frequency data and time-domain data of the electroencephalogram (EEG) signals; obtain the environmental data and blood oxygen data of the environment where the EEG signals are collected; Step 2: Calculate the relative power ratio of the EEG frequency band based on the frequency data , the standard deviation of the EEG frequency band power and the power asymmetry index of the EEG frequency band ; Calculate the brain frequency band stability difference index based on the relative power ratio of the EEG frequency band , the standard deviation of the EEG frequency band power and the power asymmetry index of the EEG frequency band ; ; Calculating the electroencephalogram pulse frequency modulation index based on time-domain data and the coefficient of variation of electroencephalogram signal activity ; According to the electroencephalogram pulse frequency modulation index and the coefficient of variation of electroencephalogram signal activity calculate the brain state stability index ; Calculate the environmental interference index based on environmental data , calculate the signal acquisition comfort adaptation index based on the environmental interference index , and further calculate the comprehensive quality index of EEG signal acquisition ; ; Calculating the cerebral blood oxygen supply-demand balance index based on blood oxygen data , calculating the cerebral functional reserve index based on the cerebral blood oxygen supply-demand balance index , and further calculating the comprehensive cerebral health assessment index ; ; According to the brain frequency band stability difference index , the brain state stability index , the comprehensive quality index of electroencephalogram signal acquisition and the comprehensive brain health assessment index Calculate the overall comprehensive index of brain function health ; Step 3: Preset the overall brain function status threshold set, and compare the comprehensive index of overall brain function health with the overall brain function status threshold set, and classify schizophrenia according to the comparison results.

2. The schizophrenia classification method based on electroencephalogram signal acquisition according to claim 1, characterized in that: Calculating the relative power ratio of electroencephalogram frequency bands , the standard deviation of electroencephalogram frequency band power and the power asymmetry index of electroencephalogram frequency bands The method is as follows: The frequency data includes the real-time power of the frequency band ; Real-time power according to frequency band Calculate the relative power ratio of electroencephalogram frequency bands , and the formula is as follows: ； Among them, is the real-time power of the th frequency band, is the serial number corresponding to different frequency bands, and the value range is ; is the number of frequency bands, and the value is a positive integer; The frequency data also includes the average value of the band power ; According to the real-time power of the frequency band and the average value of the frequency band power Calculate the standard deviation of the electroencephalogram frequency band power , and the formula is as follows: ； Among them, is the th sub-sampled band power of the th frequency band, is the serial number corresponding to different sampling times, and the value range is ; is the number of sampling times, and the value is a positive integer; The frequency data also includes the band power of the right brain region and the band power of the left brain region ; According to the frequency band power of the right brain region and the frequency band power of the left brain region Calculate the electroencephalogram frequency band power asymmetry index , and the basis formula is as follows: ； Among them, is the power of the th frequency band in the right brain region, is the power of the th frequency band in the left brain region.

3. The schizophrenia classification method based on electroencephalogram signal acquisition according to claim 2, characterized in that: Method for calculating brain frequency band stability difference index is as follows: According to the relative power ratio of brain wave frequency bands , the standard deviation of brain wave frequency band power and the power asymmetry index of brain wave frequency band calculate the brain frequency band stability difference index , and the basis formula is as follows:

4. A schizophrenia classification method based on electroencephalogram signal acquisition according to claim 3, characterized in that: Method for calculating electroencephalogram pulse frequency modulation index and coefficient of variation of electroencephalogram signal activity is as follows: The time-domain data includes the instantaneous frequency of the EEG pulse and the average instantaneous frequency ; According to the instantaneous frequency of the electroencephalogram (EEG) pulse and the average instantaneous frequency calculate the EEG pulse frequency modulation index , and the formula is as follows: ； Among them, is the instantaneous frequency of the th electroencephalogram signal pulse, is the influence coefficient corresponding to the th electroencephalogram signal pulse, is the serial number corresponding to different electroencephalogram signals, and the value range is ; is the number of electroencephalogram signals, and the value is a positive integer; The time-domain data also includes the brain wave amplitude and the average value of the brain wave amplitude ; According to the brain wave amplitude and the average value of the brain wave amplitude calculate the coefficient of variation of the electroencephalogram (EEG) signal activity , and the formula is as follows: ； Among them, is the amplitude of the brain wave at the th time point, is the serial number corresponding to different time points, and the value range is ; is the number of all time points, and the value is a positive integer.

5. A schizophrenia classification method based on electroencephalogram signal acquisition according to claim 4, characterized in that: Method for calculating brain state stability index is as follows: According to the electroencephalogram pulse frequency modulation index and the coefficient of variation of electroencephalogram signal activity calculate the brain state stability index , and the formula is as follows:

6. The schizophrenia classification method based on electroencephalogram signal acquisition according to claim 5, wherein: Method for calculating comprehensive quality index of electroencephalogram signal acquisition is as follows: The environmental data includes the actual temperature value , the rated temperature value , the actual humidity value , the rated humidity value and the electromagnetic radiation intensity value ; According to the actual temperature value , rated temperature value , actual humidity value , rated humidity value and electromagnetic radiation intensity value calculate the environmental interference index , and the formula used is as follows: ； Among them, is the weight coefficient of the temperature deviation value, and its value ranges from 0.1 to 0.3; is the weight coefficient of the humidity deviation value, and its value ranges from 0.2 to 0.3; is the weight coefficient of the electromagnetic radiation intensity value, and its value ranges from 0.5 to 0.6; and ; is the normalization constant; According to the environmental interference index Calculate the signal acquisition comfort adaptation index , and the basis formula is as follows: ； According to the environmental interference index and the signal acquisition comfort adaptation index calculate the comprehensive quality index of EEG signal acquisition , and the formula is as follows:

7. A schizophrenia classification method based on electroencephalogram signal acquisition according to claim 6, characterized in that: Method for calculating comprehensive brain health assessment index is as follows: The blood oxygen data includes arterial blood oxygen saturation value , hemoglobin value , venous blood oxygen saturation value and cardiac output value ; According to the arterial oxygen saturation value , hemoglobin value , venous oxygen saturation value and cardiac output value Calculate the cerebral oxygen supply-demand balance index , and the formula is as follows: ； According to the cerebral oxygen supply-demand balance index calculate the cerebral functional reserve index , and the formula is as follows: ； Among them, is the influence coefficient of the cerebral oxygen supply-demand balance index, is the normalization constant, and its value is a positive integer; According to the cerebral blood oxygen supply-demand balance index and the cerebral functional reserve index calculate the comprehensive cerebral health assessment index , and the formula is as follows:

8. A schizophrenia classification method based on electroencephalogram signal acquisition according to claim 7, characterized in that: Method for calculating an overall comprehensive index of brain function health is as follows: According to the brain frequency band stability difference index , brain state stability index , comprehensive quality index for electroencephalogram signal acquisition and comprehensive brain health assessment index Calculate the comprehensive index of overall brain function health , and the formula is as follows:

9. A schizophrenia classification method based on electroencephalogram signal acquisition according to claim 8, characterized in that: The criteria for classifying schizophrenia are as follows: The preset overall brain function state threshold set includes a risk-free threshold , a low-risk threshold and a severe-risk threshold ; among which, ; Compare the overall comprehensive index of brain function health with the set of thresholds of the overall state of brain function, and classify schizophrenia according to the comparison results. The criteria are as follows:

10. A schizophrenia classification system based on electroencephalogram signal acquisition, characterized in that: Including: A data acquisition module, configured to collect the frequency data and time-domain data of the EEG signals; Obtain the environmental data and blood oxygen data of the environment where the EEG signals are collected; An index calculation module, configured to calculate the relative power ratio of electroencephalogram frequency bands based on frequency data , the standard deviation of the electroencephalogram frequency band power and the power asymmetry index of the electroencephalogram frequency band ; based on the relative power ratio of the electroencephalogram frequency band , the standard deviation of the electroencephalogram frequency band power and the power asymmetry index of the electroencephalogram frequency band calculate the brain frequency band stability difference index ; Calculating the electroencephalogram (EEG) pulse frequency modulation index based on time-domain data and the coefficient of variation of EEG signal activity ; calculating the brain state stability index based on the EEG pulse frequency modulation index and the coefficient of variation of EEG signal activity ; ; Calculate the environmental interference index based on environmental data , calculate the signal acquisition comfort adaptation index based on the environmental interference index , and further calculate the comprehensive quality index of EEG signal acquisition ; ; Calculating the cerebral blood oxygen supply-demand balance index based on blood oxygen data , calculating the cerebral functional reserve index based on the cerebral blood oxygen supply-demand balance index , and further calculating the comprehensive cerebral health assessment index ; ; According to the brain frequency band stability difference index , the brain state stability index , the comprehensive quality index of electroencephalogram signal acquisition and the comprehensive brain health assessment index calculate the overall comprehensive index of brain function health ; A classification module, configured to preset a set of overall brain function state thresholds and compare the comprehensive index of overall brain function health with the set of overall brain function state thresholds, and classify schizophrenia according to the comparison result.

Images