[go: up one dir, main page]

CN118948226A - A thalamus-cortex connectivity detection method and related device - Google Patents

A thalamus-cortex connectivity detection method and related device Download PDF

Info

Publication number
CN118948226A
CN118948226A CN202411448546.3A CN202411448546A CN118948226A CN 118948226 A CN118948226 A CN 118948226A CN 202411448546 A CN202411448546 A CN 202411448546A CN 118948226 A CN118948226 A CN 118948226A
Authority
CN
China
Prior art keywords
target
stimulation
preset
coherence
thalamus
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202411448546.3A
Other languages
Chinese (zh)
Other versions
CN118948226B (en
Inventor
李小俚
关龙舟
何文博
张昊
李继芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangxi Huaheng Jingxing Medical Technology Co ltd
Original Assignee
Jiangxi Huaheng Jingxing Medical Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangxi Huaheng Jingxing Medical Technology Co ltd filed Critical Jiangxi Huaheng Jingxing Medical Technology Co ltd
Priority to CN202411448546.3A priority Critical patent/CN118948226B/en
Publication of CN118948226A publication Critical patent/CN118948226A/en
Application granted granted Critical
Publication of CN118948226B publication Critical patent/CN118948226B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • A61B5/377Electroencephalography [EEG] using evoked responses
    • A61B5/383Somatosensory stimuli, e.g. electric stimulation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4058Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
    • A61B5/4064Evaluating the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00Electrotherapy; Circuits therefor
    • A61N1/18Applying electric currents by contact electrodes
    • A61N1/32Applying electric currents by contact electrodes alternating or intermittent currents
    • A61N1/36Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
    • A61N1/36014External stimulators, e.g. with patch electrodes
    • A61N1/36025External stimulators, e.g. with patch electrodes for treating a mental or cerebral condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N1/00Electrotherapy; Circuits therefor
    • A61N1/18Applying electric currents by contact electrodes
    • A61N1/32Applying electric currents by contact electrodes alternating or intermittent currents
    • A61N1/36Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
    • A61N1/36014External stimulators, e.g. with patch electrodes
    • A61N1/3603Control systems
    • A61N1/36031Control systems using physiological parameters for adjustment
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Neurology (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Psychology (AREA)
  • Neurosurgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Developmental Disabilities (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Social Psychology (AREA)
  • Hospice & Palliative Care (AREA)
  • Child & Adolescent Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本申请公开了一种丘脑‑皮层连接性检测方法及相关装置。通过对目标检测对象实施经颅直流电刺激以及带动目标检测对象进行预设被动运动,进而采集丘脑和皮层对应刺激电极位置的脑电信号,对采集的脑电信号进行相干性分析,得到目标相干度,根据目标相干度可以确定丘脑与皮层的连接性程度。采用本申请实施例,实现了准确地对丘脑‑皮层连接性进行检测。

The present application discloses a thalamus-cortex connectivity detection method and related devices. By implementing transcranial direct current stimulation on the target detection object and driving the target detection object to perform preset passive movements, the EEG signals of the thalamus and cortex corresponding to the stimulation electrode positions are collected, and the collected EEG signals are analyzed for coherence to obtain the target coherence degree, and the degree of connectivity between the thalamus and the cortex can be determined according to the target coherence degree. By using the embodiments of the present application, accurate detection of thalamus-cortex connectivity is achieved.

Description

Thalamus-cortex connectivity detection method and related device
Technical Field
The application relates to the field of biomedical engineering, in particular to a thalamus-cortex connectivity detection method and a related device.
Background
Studies have shown that abnormal connections between thalamus and cortex are associated with a number of neuropsychiatric disorders such as anxiety, insomnia, depression and the like. By detecting the connectivity between the thalamus and the cortex, the connectivity disorder between the thalamus and the cortex can be timely found and quantified, and the method has important significance for early diagnosis and treatment of diseases and deep understanding of brain mechanisms.
The prior art generally uses Electroencephalogram (EEG) and magnetic resonance imaging (Magnetic Resonance Imaging, MRI) techniques to detect connectivity between thalamus and cortex, but this approach requires complex data analysis and is subject to multiple interference factors, resulting in less reliable results of the final detection of thalamus-cortex connectivity.
Therefore, how to accurately detect thalamus-cortex connectivity is a challenge.
Disclosure of Invention
The embodiment of the application provides a thalamus-cortex connectivity detection method and a related device, which can accurately detect thalamus-cortex connectivity.
In a first aspect, an embodiment of the present application provides a thalamus-cortex connectivity detection method, which is applied to a controller of a neural connectivity detection system, wherein the neural connectivity detection system further comprises a transcranial direct current stimulation device, an operation device and an electroencephalogram acquisition device; the operating device comprises a piston, an air compressor and a rubber band, wherein the piston is arranged in an inner hole of the air compressor, and the air compressor is used for driving the piston to move; the rubber band is used for connecting the piston and a preset part of the target detection object; the method comprises the following steps:
Determining a target stimulation parameter for the target test object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulating electrode position is the electrode position corresponding to the cerebral cortex of the target detection object;
The transcranial direct current stimulation device is used for conducting transcranial direct current stimulation on the target detection object for the duration of stimulation according to the target stimulation parameters, and the control device is used for driving the preset part of the target detection object to conduct preset passive movement;
acquiring the electroencephalogram signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively through the electroencephalogram acquisition device to obtain a first electroencephalogram signal and a second electroencephalogram signal;
Performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees;
Determining the degree of connectivity of the thalamus with the cortex based on the plurality of target coherence degrees.
In a second aspect, an embodiment of the present application provides a thalamus-cortex connectivity detection device, which is applied to a controller of a neural connectivity detection system, wherein the neural connectivity detection system further comprises a transcranial direct current stimulation device, an operation device and an electroencephalogram acquisition device; the operating device comprises a piston, an air compressor and a rubber band, wherein the piston is arranged in an inner hole of the air compressor, and the air compressor is used for driving the piston to move; the rubber band is used for connecting the piston and a preset part of the target detection object; the thalamus-cortex connectivity detection device comprises: a determination unit, a coherence analysis unit, wherein,
The determining unit is used for determining target stimulation parameters aiming at the target detection object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulating electrode position is the electrode position corresponding to the cerebral cortex of the target detection object;
The transcranial direct current stimulation device is used for conducting transcranial direct current stimulation on the target detection object for the duration of stimulation according to the target stimulation parameters, and the control device is used for driving the preset part of the target detection object to conduct preset passive movement;
acquiring the electroencephalogram signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively through the electroencephalogram acquisition device to obtain a first electroencephalogram signal and a second electroencephalogram signal;
the coherence analysis unit is used for performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees;
The determining unit is further used for determining the connectivity degree of the thalamus and the cortex according to the target coherence degrees. In a third aspect, an embodiment of the present application provides an electronic device, including: a processor, a memory, a communication interface, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the processor, the programs comprising instructions for performing the steps in the first aspect of the embodiments of the present application.
In a fourth aspect, embodiments of the present application provide a computer-readable storage medium storing a computer program for electronic data exchange, wherein the computer program causes a computer to perform part or all of the steps described in the first aspect of the embodiments of the present application.
In a fifth aspect, embodiments of the present application provide a computer program product, wherein the computer program product comprises a non-transitory computer readable storage medium storing a computer program operable to cause a computer to perform some or all of the steps described in the first aspect of the embodiments of the present application. The computer program product may be a software installation package.
According to the embodiment of the application, transcranial direct current stimulation can be implemented on the target detection object through the transcranial direct current stimulation device, the preset part of the target detection object is driven by the operation device to perform preset passive movement, the electroencephalogram acquisition device is used for respectively acquiring the electroencephalogram signals of the target detection object at the position of the stimulation electrode corresponding to the thalamus and the position of the second stimulation electrode corresponding to the cortex, so that the first electroencephalogram signals and the second electroencephalogram signals are obtained, the coherence analysis is performed on the first electroencephalogram signals and the second electroencephalogram signals, a plurality of target coherence is obtained, the connectivity degree of the thalamus and the cortex is determined according to the plurality of target coherence degrees, and further, the accurate detection on the thalamus-cortex connectivity is realized.
Drawings
In order to more clearly describe the embodiments of the present application or the technical solutions in the background art, the following description will describe the drawings that are required to be used in the embodiments of the present application or the background art.
FIG. 1 is a schematic diagram of a neural connectivity detection system according to an embodiment of the present application;
FIG. 2 is a functional block diagram of a transcranial direct current stimulation device according to an embodiment of the present application;
FIG. 3 is a functional block diagram of an operator provided in accordance with an embodiment of the present application;
FIG. 4 is a schematic flow chart of a method for detecting thalamus-cortex connectivity according to an embodiment of the present application;
Fig. 5 is an application scenario schematic diagram of an operating device provided in an embodiment of the present application;
FIG. 6 is a block diagram showing functional units of a thalamus-cortex connectivity detection device according to an embodiment of the present application;
Fig. 7 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
In order that those skilled in the art will better understand the present application, a technical solution in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present application, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The terms first, second and the like in the description and in the claims and in the above-described figures are used for distinguishing between different objects and not necessarily for describing a sequential or chronological order. Furthermore, the terms "comprise" and "have," as well as any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those listed steps or elements but may include other steps or elements not listed or inherent to such process, method, article, or apparatus.
Reference herein to "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment may be included in at least one embodiment of the application. The appearances of such phrases in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments. Those of skill in the art will explicitly and implicitly appreciate that the embodiments described herein may be combined with other embodiments.
The following describes related content, concepts, meanings, technical problems, technical schemes, beneficial effects and the like related to the embodiment of the application.
Some terms involved in the present application will be explained first:
Thalamus: the thalamus is a relay of sensory information transmission, and almost all sensory signals (except for olfaction) pass through the thalamus before reaching the cortex, where they are initially processed and integrated and then transmitted to the corresponding areas of the cortex.
Cortex: i.e., the cortex, which is responsible for handling advanced neurological functions including thought, consciousness, memory, language, learning, vision, hearing, decision making, and motor control.
Transcranial direct current stimulation: transcranial direct current stimulation (TRANSCRANIAL DIRECT current stimulation, tDCS for short) is a non-invasive neural electrical stimulation technique that applies weak current to the scalp of the brain through positive and negative electrodes to stimulate specific brain regions, regulate neural activity and/or excitability of the cerebral cortex, and in general, anodic (positive) stimulation increases neuronal excitability while cathodic (negative) stimulation inhibits neuronal excitability.
10-20 Standard guide system: the 10-20 standard lead system is a naming and positioning system for standardized electrode placement in electroencephalography, recording electrical activity in different areas of the brain by placing electrodes at specific percentage locations on the scalp, where "10" and "20" represent the percentage of the distance between the electrodes to the widest part of the head and the distance from the front to the back of the head, respectively.
Anti-aliasing filter: the anti-aliasing filter (anti-ALIAS FILTER) is a low-pass filter, which mainly functions to filter out frequency components in the signal above the nyquist frequency before the sampling process of the signal. In digital signal processing, according to the nyquist sampling theorem, when an analog signal is digitized, the sampling frequency must be at least twice the highest frequency of the signal to avoid aliasing. Aliasing refers to the fact that if the sampling frequency is insufficient to capture all frequency components of the signal, frequency components higher than half the sampling frequency may erroneously appear in the low frequency portion of the sampled signal, resulting in loss or distortion of the original signal information.
Nyquist frequency: the nyquist frequency (Nyquist frequency) is the minimum sampling frequency that needs to be defined to prevent aliasing of the signal.
Butterworth filter: the butterworth filter (Butterworth filter) is one of electronic filters that has a flat frequency response at the pass frequency (i.e., the cut-off frequency of the filter) so that it does not cause phase distortion or amplitude variation of the signal at that frequency.
Hanning window: hanning window (Hanning) is a window function commonly used in signal processing, particularly when performing frequency domain analysis, that is used to reduce the weighting function of abrupt effects across samples, thereby reducing spectral leakage.
Referring to fig. 1, fig. 1 is a schematic structural diagram of a nerve connectivity detection system according to an embodiment of the present application, and as shown in the drawing, the nerve connectivity detection system 10 includes a transcranial direct current stimulation device 101, a manipulation device 102, an electroencephalogram acquisition device 103 and a controller 104.
The nerve connectivity detection system 10 is a system for detecting connectivity between thalamus and cortex of a subject, which regulates excitability of cerebral cortex by applying weak direct current to scalp of the subject through the transcranial direct current stimulation device 101, and performs passive exercise on the subject in combination with the manipulation device 102, and records electrical activity of brain of the subject during stimulation and exercise through the electroencephalogram acquisition device 103 and acquires electroencephalogram signals for evaluation and analysis to determine the degree of connectivity between thalamus and cortex of the subject.
The controller 104 is the core control unit of the neuro-connectivity detection system 10, which is primarily responsible for determining stimulation parameters for the subject, including: stimulus intensity, first stimulus electrode position, second stimulus electrode position, and stimulus duration. The subject is then subjected to a pre-stimulus test to ensure that the subject can withstand the stimulus intensity, after which the controller 104 controls the transcranial direct current stimulation device 101 to apply weak direct current to the scalp of the subject for stimulation, and controls the manipulation device 102 to perform passive movements on the subject, while controlling the electroencephalogram acquisition device 103 to record and acquire electroencephalogram signals of the brain of the subject. The degree of connectivity between the thalamus and the cortex is determined by acquiring and evaluating and analyzing the brain electrical signals of the thalamus and the cortex, and determining the synchronicity of the brain electrical signals between the thalamus and the cortex.
In embodiments of the present application, the transcranial direct current stimulation device 101 may apply a weak direct current to the scalp of a subject to modulate the excitability of the cerebral cortex, which may alter the resting membrane potential of neurons, thereby affecting the activity and connectivity of brain regions. The operating device 102 mainly comprises a piston, an air compressor and a rubber band, and is responsible for controlling fingers or limbs of a subject, wherein the piston in the air compressor moves under the action of air pressure, and the connected rubber band drives a preset part of the subject to perform passive movement, so that the natural movement of the subject is simulated. The electroencephalogram acquisition device 103 mainly comprises an electroencephalogram electrode, a lead, a signal amplifier, an analog-to-digital converter, data recording and analysis software and the like, and the electroencephalogram acquisition device provides data support for analyzing connectivity between thalamus and cortex by placing the electrode in a specific brain region and capturing changes of electroencephalogram signals.
Referring to fig. 2, fig. 2 is a functional block diagram of a transcranial direct current stimulation device according to an embodiment of the present application, and as shown in fig. 2, the transcranial direct current stimulation device 101 includes a processing device 201, a storage device 202, an interface device 203, a communication device 204, a display device 205, an input device 206, a speaker 207, and a stimulation device 208.
The transcranial direct current stimulation device 101 is mainly used for applying transcranial direct current stimulation on the scalp of a subject. In the transcranial direct current stimulation device 101, the processing device 201 may be a Central Processing Unit (CPU) or a Micro Controller (MCU) or the like, which is mainly responsible for executing preset program instructions, processing data, and controlling the overall operation of the transcranial direct current stimulation device 101. The storage 202 is used to store operating software of the device, a preset program, and user data, and may be a Read Only Memory (ROM), a Random Access Memory (RAM), a hard disk, and the like. The interface means 203 provides a physical connection, such as a USB interface and a headset interface, allowing the user to connect to an external device, such as a computer or a data recorder. The communication device 204 is responsible for wired or wireless data transmission, including Wi-Fi and bluetooth communications. These communication modes enable the device to exchange data with external systems, facilitating remote monitoring and program updating. The display device 205, which may be a liquid crystal display or a touch screen, is used to provide visual information to the user of the status of the device, stimulation parameters, and operational feedback. The input device 206, which may be a touch screen and a keyboard, is used to provide an intuitive way for the operator to enter instructions and adjust settings so that the operator can conveniently control the stimulation process. The speaker 207 may be a sound feedback providing operation confirmation, such as an alert sound when a user inputs a command or adjusts settings through an input device, or an abnormal situation or a matter requiring user attention is detected, and may be used to sound an alarm or notification sound, or the like. Stimulation device 208 is the core component of transcranial direct current stimulation device 101 that is responsible for generating and delivering electrical stimulation to a target individual, and by precisely controlling the intensity, frequency, and duration of the current, stimulation device 208 can modulate specific neural pathways.
Referring to fig. 3, fig. 3 is a functional block diagram of an operating device according to an embodiment of the present application, and as shown in fig. 3, the operating device 102 includes a piston 301, an air compressor 302, and a rubber band 303.
Wherein the manipulation device 102 is used for controlling limb movement of a subject to assist in acquisition and analysis of brain electrical signals. The manipulator 102 is mainly composed of three parts, namely a piston 301, an air compressor 302 and a rubber band 303, which cooperate together to achieve manipulation control of the finger or limb of the subject.
The piston 301 is a component for transmitting power in the operator 102, which can be moved by the pressure provided by the air compressor 302, thereby simulating natural motion. The piston 301 is connected to the air compressor 302 by a plastic tube, which is capable of telescopic movement when the air compressor 302 provides the pressure of the gas.
The air compressor 302 is a power source of the manipulator 102, and can adjust the air flow through a built-in controllable electric valve, and the opening and closing of the electric valve controls the movement of the piston 301, so as to adjust the movement of the finger or limb of the subject.
The rubber band 303 plays a role in connecting and transmitting motion in the operating device 102, one end of the rubber band is connected with the piston 301, the other end of the rubber band can be fixed on the thumb of a subject, and when the piston 301 moves, the piston 301 can drive the thumb to perform corresponding bending or straightening motion to simulate the actual hand motion.
In the embodiment of the present application, when the manipulating device 102 is used, the elastic band 303 is firstly fixed on the thumb of the subject to ensure the connection stability, then the air compressor 302 is started, and the movement of the piston 301 is controlled by adjusting the opening and closing of the electric valve, so as to drive the thumb to perform the preset passive movement. And furthermore, in the whole passive movement process, the movement state and the brain electrical signals of the subject can be monitored in real time.
Referring to fig. 4, fig. 4 is a schematic flow chart of a thalamus-cortex connectivity detection method according to an embodiment of the present application, wherein the method is applied to a controller of a neural connectivity detection system, and the neural connectivity detection system further comprises a transcranial direct current stimulation device, an operation device and an electroencephalogram acquisition device; the operating device comprises a piston, an air compressor and a rubber band, wherein the piston is arranged in an inner hole of the air compressor, and the air compressor is used for driving the piston to move; the rubber band is used for connecting the piston and a preset part of the target detection object; the method includes, but is not limited to, the steps of:
S401, determining target stimulation parameters aiming at the target detection object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulating electrode position is the electrode position corresponding to the cerebral cortex of the target detection object.
The target detection object is a subject needing to detect connectivity of thalamus and cortex, before the target detection object is stimulated, the stimulation parameters of the transcranial direct current stimulation device need to be determined, and the stimulation parameters are set for the transcranial direct current stimulation device, so that the transcranial direct current stimulation device carries out transcranial direct current stimulation for the target detection object, and the stimulation parameters comprise target stimulation intensity, first stimulation electrode position, second stimulation electrode position and stimulation duration.
In a specific embodiment, the electrode position can be marked on the scalp of the target detection object according to the 10-20 standard lead system, the C 3 position is selected as the corresponding electrode position of the thalamus of the target detection object, the first stimulation electrode position is the C 3 position, the F p2 position is selected as the corresponding electrode position of the cerebral cortex of the target detection object, the second stimulation electrode position is the F p2 position, the target stimulation intensity is further set to be 1mA, and the stimulation duration is set to be 20 minutes.
Further, an EEG cap is provided to the head of the subject, wherein the EEG cap marks the specific location of the electrode placement according to the 10-20 standard lead system, then the anode electrode of the EEG cap is placed at the first stimulation electrode location and the cathode electrode of the EEG cap is placed at the second stimulation electrode location. To optimize the transmission efficiency of the stimulation current, the impedance between the electrode and the scalp of the target detection object is reduced, and the conductive gel may be applied between the electrode and the scalp of the target detection object to reduce the resistance and improve the stimulation effect.
S402, performing transcranial direct current stimulation on the target detection object for the stimulation duration time according to the target stimulation parameters through the transcranial direct current stimulation device, and driving the preset part of the target detection object to perform preset passive movement through the operation device.
The thalamus can receive and process nerve signals from the movement part and transmit the nerve signals to the cerebral cortex for further processing when the target detection object actively activates the preset part. The preset passive movement is preset passive movement, and the preset part of the target detection object is driven by additional equipment to perform preset movement so as to induce communication of nerve signals between thalamus and cerebral cortex of the target detection object.
In a specific embodiment, the stimulation parameters are input into the transcranial direct current stimulation device so that the transcranial direct current stimulation device can perform transcranial direct current stimulation on the target test object for a stimulation duration. Meanwhile, the air compressor is controlled to drive the piston inside the air compressor to regularly stretch and retract so as to drive the preset part of the target detection object to perform preset passive movement.
In one possible embodiment, the preset portion is the right thumb of the target detection object, the right thumb is connected with the piston of the operating device through the rubber band, and the piston can perform all-round passive movement, such as straightening, bending, inward or outward rotation, on the right thumb of the target detection object with the help of the piston.
Further, the thumb of the right hand of the target detection object is subjected to preset passive movement within 20 minutes of stimulation, wherein the movement period of each group is 20 seconds, 7 passive movements are performed within 20 seconds, the duration of each passive movement can be 2-3 seconds, and a rest is performed for 20 seconds after each group of 7 passive movements. During a total of 20 minutes of stimulation, 30 sets of passive exercises were performed, the total number of passive exercises being 210, including in particular 30 passive exercises and 30 rest.
For convenience of understanding, please refer to fig. 5, fig. 5 is a schematic diagram of an application scenario of an operating device according to an embodiment of the present application. As shown in the figure, the manipulating device includes a piston, an air compressor, and a rubber band for connecting the thumb of the right hand of the target object to be detected with the piston, and at the same time, the piston is a movable part inside the air compressor, which is connected with the thumb of the right hand of the target object to be detected through the rubber band. When the operating device is started, the air compressor starts to work, the air compressor is internally inflated, the piston moves outwards, the right thumb of the target detection object is driven to bend, when the air compressor is internally discharged, the piston is driven to move towards the air compressor, the right thumb of the target detection object is driven to straighten, and the piston in the air compressor is driven to regularly stretch and retract, so that the preset part of the target detection object performs preset passive movement.
Optionally, before performing, by the transcranial direct current stimulation device, transcranial direct current stimulation on the target detection object for the duration of the stimulation according to the target stimulation parameter in the step S402, the method may further include the following steps:
a21, acquiring the current stimulation intensity of the transcranial direct current stimulation device at the current moment;
a22, increasing the current stimulation intensity to a preset stimulation intensity in a first preset time period;
a23, continuously stimulating the target detection object according to the preset stimulation intensity in a second preset time period, and monitoring the physiological signal of the target detection object to obtain a target physiological signal; the starting time of the second preset time period is later than the ending time of the first preset time period;
a24, determining the preset stimulus intensity as the target stimulus intensity when the target physiological signal meets a first preset condition; the first preset condition is that the target physiological signal is in a preset safety interval;
A25, when the target physiological signal does not meet the first preset condition, adjusting the preset stimulus intensity based on the target physiological signal to obtain standard stimulus intensity;
a26, determining the standard stimulus intensity as the target stimulus intensity.
It should be noted that, before the transcranial direct current stimulation device performs the transcranial direct current stimulation for the duration of the stimulation on the target test object according to the target stimulation parameters, a pre-stimulation test needs to be performed on the target test object to ensure that the target test object can withstand the stimulation.
The first preset time period is preset, and the transcranial direct current stimulation device slowly improves the stimulation intensity in the first preset time period. The second preset time period is preset, and the fact that the transcranial direct current stimulation device continuously stimulates the target detection object in the second preset time period is indicated, wherein the starting time of the second preset time period is later than the ending time of the first preset time period. The preset condition is that the physiological signal of the target detection object is in a preset safety interval, which means that the target detection object can bear the stimulation of the transcranial direct current stimulation device. The preset safety interval is an acceptable range of physiological parameters of a preset target detection object, and can be adaptively adjusted based on a standard safety interval according to influence possibly generated by nerve connectivity detection. The preset stimulus intensity may be a target stimulus intensity set in advance, or may be a stimulus intensity set in advance according to the condition of the subject.
In a specific embodiment, the target detection object is subjected to a pre-stimulation test by the transcranial direct current stimulation device, so that the current stimulation intensity of the transcranial direct current stimulation device at the current moment can be obtained, and is smaller than the preset stimulation intensity, wherein the preset stimulation intensity can be preset for the transcranial direct current stimulation device, the current stimulation intensity is increased to the preset stimulation intensity in a first preset time period, the target detection object is continuously stimulated according to the preset stimulation intensity in a second preset time period, and meanwhile, the physiological signals of the target detection object are monitored, so that the target physiological signals can be obtained, wherein the target physiological signals comprise, but are not limited to, heart rate, blood pressure, respiratory frequency, body temperature, skin electric activity, electroencephalogram signals and the like, and whether the target detection object can bear the preset stimulation intensity or not is determined by analyzing the target physiological signals.
Further, when the target physiological signal meets the first preset condition, that is, the target physiological signal is in the preset safety interval, the preset stimulation intensity is determined as the target stimulation intensity, and the transcranial direct current stimulation device can perform stimulation of the target stimulation intensity on the target detection object. When the target physiological signal does not meet the first preset condition, the target detection object is not tolerant to the stimulus of the preset stimulus intensity in the pre-test process, the preset stimulus intensity needs to be adjusted based on the target physiological signal to obtain the standard stimulus intensity, the standard stimulus intensity obtained after adjustment is determined to be the target stimulus intensity, and therefore the target detection object can be ensured to bear the target stimulus intensity through the pre-stimulus test to the target detection object, so that connectivity detection interference between thalamus and cortex is reduced, and accuracy of connectivity detection and safety of the target detection object are ensured.
In one possible embodiment, the current stimulation intensity of the transcranial direct current stimulation device at the current moment is 0mA, after the transcranial direct current stimulation device is started, the stimulation intensity is slowly increased to 1mA within 15s, the target detection object is stimulated for 10s based on the stimulation intensity of 1mA, meanwhile, the heart rate change of the target detection object is monitored, if the heart rate change is judged to be always kept within 70-110 times/min, the target detection object can bear the stimulation intensity of 1mA, otherwise, the current stimulation intensity is adjusted, so that the target detection object can bear the adjusted stimulation intensity. After determining the stimulus intensity that the target test subject can withstand, the stimulation of the transcranial direct current stimulation device may be slowly dropped to 0 within 15 seconds in preparation for the neuroconnectivity test.
Optionally, when the target physiological signal includes a heart rate, the step a25 may include adjusting the preset stimulus intensity based on the target physiological signal to obtain a standard stimulus intensity when the target physiological signal does not meet the first preset condition, and may include the following steps:
b21, acquiring a target safe interval and a standard heart rate value of the heart rate;
b22, obtaining a difference value between the heart rate and the standard heart rate value to obtain a heart rate difference value;
b23, determining the adjustment parameters of the preset stimulus intensity according to the heart rate difference value to obtain target adjustment parameters;
b24, when the heart rate is not in the target safety interval, determining the deviation degree of the heart rate based on the target safety interval, and obtaining a target deviation degree;
B25, determining a target adjustment factor corresponding to the target deviation degree;
And B26, adjusting the preset stimulus intensity according to the target adjustment factor and the target adjustment parameter to obtain the standard stimulus intensity.
When the target physiological signal to be detected and analyzed is a heart rate, a target safety interval is set, wherein the target safety interval represents a safety interval of heart rate change of a target detection object when the target detection object receives stimulation, and a standard heart rate value can be set, and the standard heart rate value represents an expected heart rate level of the target detection object after stimulation. Then, according to the deviation between the heart rate variation data and the preset standard, the adjustment parameters and the adjustment factors of the preset stimulus intensity can be further determined so as to adjust the preset stimulus intensity.
In a specific embodiment, an average value of heart rates of the target detection objects detected in the second preset time period may be used as the estimated heart rate, and then, a target safety interval and a standard heart rate value of the heart rate may be obtained. Obtaining a difference value between the heart rate and a standard heart rate value to obtain a heart rate difference value, determining an adjustment parameter of the preset stimulation intensity according to the heart rate difference value to obtain a target adjustment parameter, determining the deviation degree of the heart rate based on the target safety interval when the heart rate is not in the target safety interval to obtain a target deviation degree, determining a target adjustment factor corresponding to the target deviation degree, and adjusting the preset stimulation intensity according to the target adjustment factor and the target adjustment parameter to obtain the standard stimulation intensity.
In one possible embodiment, the preset stimulus intensity is set to 1mA, the target safety interval of the heart rate of the target detection object is set to 70-110 times/min, the standard heart rate value is set to 90 times/min, the average heart rate of the target detection object is detected to be 60 times/min in the continuous stimulus period within 10s, the heart rate difference is 30 times/min, the target adjustment parameter is 6, namely, the target adjustment parameter is the ratio of the heart rate difference to 5 when the heart rate difference is set to 5 times/min. The target adjustment parameter of each unit is set to have an adjustment amplitude of 0.05mA for the stimulus intensity, and the target adjustment parameter is set to 6 at this time, so that the adjustment amplitude is set to 0.3mA at this time. And meanwhile, if the heart rate of the target detection object is detected not to exceed the lower limit of the target safety zone, determining that the difference between the heart rate and the lower limit of the target safety zone is 10 according to the heart rate and the lower limit of the target safety zone, and determining that the target deviation degree is 1/4 if the zone length of the target safety zone is 40, and determining that the target adjustment factor is 1/4. The preset stimulus intensity is adjusted based on the following formula: standard stimulus intensity = preset stimulus intensity + (0.05 x target adjustment parameter x (1 + target adjustment factor)). A standard stimulation intensity of 1.375mA was then obtained.
In one possible embodiment, the preset stimulus intensity is set to 1mA, the target safety interval of the heart rate of the target detection object is set to 70-110 times/min, the standard heart rate value is set to 90 times/min, the average heart rate of the target detection object is detected to be 115 times/min in the continuous stimulus period within 10s, the heart rate difference is 25 times/min, the target adjustment parameter is 5, that is, the ratio of the heart rate difference to 5 is set for every 5 times/min of the heart rate difference. The target adjustment parameter of each unit is set to have an adjustment amplitude of 0.05mA for the stimulus intensity, and the target adjustment parameter is set to be 5 at this time, so that the adjustment amplitude is set to be 0.25mA at this time. And meanwhile, if the heart rate of the target detection object exceeds the upper limit of the target safety zone, determining that the difference between the heart rate and the upper limit of the target safety zone is 5 according to the upper limit of the heart rate and the upper limit of the target safety zone, and determining that the target deviation degree is 1/8 if the zone length of the target safety zone is 40, and determining that the target adjustment factor is 1/8. The preset stimulus intensity is adjusted based on the following formula: standard stimulus intensity = preset stimulus intensity- (0.05 x target adjustment parameter x (1 + target adjustment factor)). A standard stimulation intensity of 0.71875mA was then obtained.
S403, respectively acquiring the brain electrical signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position through the brain electrical acquisition device to obtain a first brain electrical signal and a second brain electrical signal.
The electroencephalogram acquisition device can comprise an electroencephalogram electrode, a wire, a signal amplifier, an analog-to-digital converter, data recording and analysis software and the like, is connected with the stimulation electrode of the EEG cap, and can acquire and record the electroencephalogram signal of the target detection object at the specific stimulation electrode position.
In a specific embodiment, the first stimulating electrode position and the second stimulating electrode position are connected with an electroencephalogram acquisition device, weak electroencephalogram signals are amplified through a signal amplifier of the electroencephalogram acquisition device during the period that the transcranial direct current stimulation device stimulates a target detection object, analog electroencephalogram signals are converted into digital signals through an analog-to-digital converter, and an anti-aliasing filter can be applied to remove signal components with frequencies higher than nyquist frequency before the signals are converted into digital signals by the analog-to-digital converter so as to reduce aliasing. The sampled digital signals comprise brain electrical activity signals in a plurality of frequency bands, wherein the plurality of frequency bands comprise Theta wave bands, alpha wave bands, beta wave bands and the like, the brain electrical activity signals are recorded through data recording and analysis software without limitation, and the first brain electrical signal and the second brain electrical signal can be obtained.
S404, performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees.
The coherence analysis of the first electroencephalogram signal and the second electroencephalogram signal can obtain a plurality of target coherence degrees which can be reflected in the same frequency band, and the synchronism of the activities of the two electroencephalogram signals.
In a specific embodiment, acquiring an electroencephalogram signal at a first stimulation electrode position to obtain a first electroencephalogram signal, acquiring an electroencephalogram signal at a second acquisition electrode position to obtain a second electroencephalogram signal, and performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to evaluate the functional connection strength between thalamus and cortex.
In the embodiment of the application, the collected electroencephalogram signals comprise signals in a plurality of frequency bands, so that the signals in the specific frequency bands can be extracted for coherence analysis, specifically, a band-pass filter is used for preprocessing the first electroencephalogram signals and the second electroencephalogram signals, the signals in the specific frequency bands are extracted, fourier transformation is carried out on the filtered signals, the signals are converted into frequency domain signals from time domain signals, amplitude and phase information on each frequency point are calculated, and the coherence of the first electroencephalogram signals and the second electroencephalogram signals in the specific frequency bands, namely, the synchronicity of amplitude spectrums of the first electroencephalogram signals and the second electroencephalogram signals in the specific frequency bands can be obtained, a plurality of target coherence degrees can be obtained based on the correlation degree of the coherence degrees of thalamus and cortex in the specific frequency bands, and the high coherence degree represents the activity high synchronicity of the thalamus and cortex in the specific frequency bands, namely, the connectivity between the thalamus and cortex is strong.
Optionally, in step S404, coherence analysis is performed on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees, which may include the following steps:
A41, respectively carrying out filtering processing on the first electroencephalogram signal and the second electroencephalogram signal based on a preset filter to obtain a third electroencephalogram signal and a fourth electroencephalogram signal;
A42, discretizing the third electroencephalogram signal and the fourth electroencephalogram signal through a preset window function respectively to obtain a first time domain signal and a second time domain signal;
a43, performing Fourier transform on the first time domain signal and the second time domain signal respectively to obtain a first frequency domain signal and a second frequency domain signal;
A44, determining frequency points in the first frequency domain signal and the second frequency domain signal under the same frequency to obtain a plurality of frequency point groups; each frequency point group comprises a first frequency point set and a second frequency point set; each frequency point set comprises a plurality of frequency points; the first frequency point set is the frequency point set of the first frequency domain signal; the second frequency point set is the frequency point set of the second frequency domain signal;
A45, determining the coherence of the first frequency domain signal and the second frequency domain signal according to the plurality of frequency point groups to obtain a plurality of target coherence degrees.
Specifically, the first electroencephalogram signal and the second electroencephalogram signal are preprocessed by using a preset filter, signals in a specific frequency band are extracted, wherein the specific frequency band can be a signal in an Alpha wave band, the preset filter can be a 2-order filter of Butterworth, the first electroencephalogram signal and the second electroencephalogram signal are filtered by using the 2-order filter of Butterworth, the band-pass filtering range is 8-13Hz, and based on the fact, the electroencephalogram signals in the Alpha wave band can be extracted from the first electroencephalogram signal and the second electroencephalogram signal, and the third electroencephalogram signal and the fourth electroencephalogram signal can be obtained respectively.
Further, a preset window function is adopted to discretize the third electroencephalogram signal and the fourth electroencephalogram signal respectively, so that the third electroencephalogram signal and the fourth electroencephalogram signal are converted into time domain signals suitable for Fourier transform, wherein the preset window function can be a Hanning window function, discretization processing is conducted on the third electroencephalogram signal and the fourth electroencephalogram signal through the Hanning window function respectively, and a first time domain signal and a second time domain signal are obtained, and a calculation formula of the discretization processing is as follows:
In the above formula, ω (n) is the value of the window function at the discrete time point n; n is an index of discrete time points, n=0, 1,2,..n-1, where N is the total number of signal samples; n is the length of the window function, i.e. the total number of signal samples.
Further, fourier transforming the first time domain signal and the second time domain signal to obtain a first frequency domain signal and a second frequency domain signal, where the fourier transforming may use Fast Fourier Transform (FFT) to obtain the first frequency domain signal and the second frequency domain signal, where a formula of the fourier transform is as follows:
Where F (k) is the value of the frequency domain signal at the kth frequency point, k=n=0, 1,2,..n-1, where N is the total number of signal samples; ω (N) is the value of the time domain signal after window function processing at the N-th time point, n=0, 1,2,..; i is an imaginary unit.
Further, frequency points in the first frequency domain signal and the second frequency domain signal under the same frequency are determined to obtain a plurality of frequency point groups, each frequency point group comprises a first frequency point set and a second frequency point set, each frequency point set comprises a plurality of frequency points, wherein the first frequency point set is the frequency point set of the first frequency domain signal, and the second frequency point set is the frequency point set of the second frequency domain signal.
It should be noted that, since the filtering range of the preset filter is 8-13Hz, if the coherence of the third electroencephalogram signal and the fourth electroencephalogram signal is directly calculated in the frequency range, the coherence of the signals may not be accurately reflected due to insufficient frequency resolution, so that the 8-13Hz may be divided into a plurality of frequency sub-bands, such as 8-8.5Hz, 8.5-9Hz, and the like, frequency points in the frequency sub-bands are sampled as the same frequency point, a set of frequency point sets may be obtained, and a plurality of frequency point sets may be obtained based on the plurality of frequency sub-bands.
Specifically, the coherence of the two frequency domain signals on each frequency point can be calculated according to the first frequency point set and the second frequency point set, the correlation analysis of the amplitude spectrums of the first electroencephalogram signal and the second electroencephalogram signal on specific frequencies is determined, a plurality of target coherence degrees are obtained, the coherence degrees of the targets are analyzed and can be accurately reflected in an alpha frequency band, and the synchronicity of the first electroencephalogram signal and the synchronicity of the second electroencephalogram signal are accurately reflected, so that connectivity information about thalamus and cortex is provided.
Optionally, the step a45 of determining the coherence of the first frequency domain signal and the second frequency domain signal according to the plurality of frequency point groups to obtain a plurality of target coherence degrees may include the following steps:
b41, determining a target frequency point group; the target frequency point group is any one of the frequency point groups;
B42, determining the self-power spectral density of the first frequency domain signal according to a first frequency point set in the target frequency point set to obtain a first power spectral density;
b43, determining the self-power spectral density of the second frequency domain signal according to a second frequency point set in the target frequency point set to obtain a second power spectral density;
b44, determining the cross power spectral density of the first frequency domain signal and the second frequency domain signal according to the target frequency point group to obtain a third power spectral density;
And B45, determining target coherence according to the first power spectral density, the second power spectral density and the third power spectral density.
Specifically, any one frequency point group is selected from the plurality of frequency point groups to serve as a target frequency point group, and the power spectrum densities of the first frequency domain signal and the second frequency domain signal are calculated according to the target frequency point group. Determining the self-power spectral density of the first frequency domain signal according to a first frequency point set in the target frequency point set to obtain a first power spectral density, and determining the self-power spectral density of the second frequency domain signal according to a second frequency point set in the target frequency point set to obtain a second power spectral density, wherein the self-power spectral density has the following calculation formula:
In the above formula, P (k) is the power spectral density at the kth frequency point; f 2 is the upper limit of the sub-band of interest; f 1 is the lower limit of the sub-band of interest; f (k) is the value of the frequency domain signal at the kth frequency point, k=n=0, 1, 2.
Further, determining the cross power spectral density of the first frequency domain signal and the second frequency domain signal according to the target frequency point group to obtain a third power spectral density, wherein the calculation formula of the cross power spectral density is as follows:
In the above formula, P xy (k) is the cross-power spectral density of the two signals x and y at the frequency point k, where is the cross-power spectral density of the first frequency domain signal and the second frequency domain signal at the frequency point k; f 2 is the upper limit of the sub-band of interest; f 1 is the lower limit of the sub-band of interest; f x (k) is the value of the first frequency domain signal at the kth frequency bin, k=n=0, 1,2,..n-1, where N is the total number of signal samples; f y (k) is the value of the second frequency domain signal at the kth frequency point.
Then, determining the target coherence according to the first power spectral density, the second power spectral density and the third power spectral density, wherein a calculation formula of the coherence is as follows:
In the above formula, C xy (k) is the coherence degree of the first frequency domain signal and the second frequency domain signal at the frequency point k, the value range is between 0 and 1, when the coherence degree is close to 1, the two signals are highly correlated, and when the coherence degree is close to 0, the two signals are basically uncorrelated; p xy (k) is the cross power spectral density of the first frequency domain signal and the second frequency domain signal, i.e. the third power spectral density; p x (k) is the self-power spectral density of the first frequency domain signal, i.e. the first power spectral density; p y (k) is the self power spectral density of the second frequency domain signal, i.e. the second power spectral density.
In a specific embodiment, the coherence analysis is performed on the first electroencephalogram signal and the second electroencephalogram signal, so that whether the thalamus and the cortex are synchronous in activity on a specific frequency can be determined, and according to the obtained multiple target coherence degrees, the connectivity strength between the thalamus and the cortex can be evaluated.
S405, determining the connectivity degree of thalamus and cortex according to the target coherence degrees.
In a specific embodiment, the overall coherence of the first electroencephalogram signal and the second electroencephalogram signal is determined through a plurality of target coherence degrees, the degree of connectivity between the thalamus and the cortex is analyzed according to the overall coherence degree, if the overall coherence degree is higher than or equal to a preset first threshold value, strong connectivity between the thalamus and the cortex can be determined, if the overall coherence degree is lower than or equal to a preset second threshold value, weak connectivity between the thalamus and the cortex can be determined, and if the overall coherence degree is between the preset second threshold value and the preset first threshold value, connectivity between the thalamus and the cortex needs to be further determined according to the proximity degree of the overall coherence degree relative to the two threshold values. The preset first threshold value and the preset second threshold value are preset threshold values, and can be set according to conditions of different crowds, such as different ages, sexes, health conditions and the like.
Optionally, the determining the connectivity degree of the thalamus and the cortex according to the target coherence degrees in the step S405 may include the following steps:
a501, determining an average value of the plurality of target coherence degrees to obtain a target average value;
a502, when the target average value is greater than or equal to a preset first coherence, determining that connectivity of thalamus and cortex is in a strong connection state;
a503, when the target average value is smaller than or equal to a preset second correlation, determining that connectivity of thalamus and cortex is in a weak connection state;
A504, when the target average value is larger than the preset second coherence degree and the target average value is smaller than the preset first coherence degree, fitting according to the target coherence degrees to obtain a target coherence curve; the horizontal axis of the target coherence curve is frequency and the vertical axis is coherence;
A505, determining target weighted coherence according to the target coherence and the target average value;
A506, determining a dividing line according to the target weighted coherence; the parting line is parallel to the transverse axis;
a507, determining a first area value of a region above the dividing line and a second area value of a region below the dividing line in the target coherence curve;
A508, determining the ratio between the first area value and the second area value to obtain a target ratio;
a509, acquiring a preset ratio threshold;
A510, when the target ratio is greater than or equal to the preset ratio threshold, determining that connectivity of thalamus and cortex is in a strong connection state;
And A511, when the target ratio is smaller than the preset ratio threshold, determining that the connectivity of the thalamus and the cortex is in a weak connection state.
The method comprises the steps of presetting a first coherence degree, presetting a second coherence degree, judging that the connectivity between the thalamus and the cortex is strong connection coherence degree, presetting a threshold value of a preset ratio, and indicating a standard value of a ratio of the number of target coherence degrees larger than the target weighted coherence degree to the number of target coherence degrees smaller than the target weighted coherence degree in a plurality of target coherence degrees, wherein the preset first coherence degree is preset, judging that the connectivity between the thalamus and the cortex is strong connection coherence degree, and the preset second coherence degree is preset, judging that the connectivity between the thalamus and the cortex is weak connection coherence degree.
In a specific embodiment, an average value of a plurality of target coherence degrees is determined, so as to obtain a target average value, when the target average value is greater than or equal to a preset first coherence degree, connectivity between the thalamus and the cortex can be determined to be in a strong connection state, when the target average value is less than or equal to a preset second coherence degree, connectivity between the thalamus and the cortex can be determined to be in a weak connection state, when the target average value is greater than the preset second coherence degree and the target average value is less than the preset first coherence degree, the current plurality of target coherence degrees can be considered to be unstable coherence degrees, overall trend analysis is required to be performed on the plurality of target coherence degrees, whether the plurality of target coherence degrees tend to be in the preset first coherence degrees is judged, if the plurality of target coherence degrees tend to be in the preset first coherence degrees, connectivity between the thalamus and the cortex can be determined to be in the strong connection state, otherwise, the connectivity between the thalamus and the cortex can be determined to be in the weak connection state.
Further, fitting is carried out according to the plurality of target coherence degrees to obtain a target coherence curve, the horizontal axis of the target coherence curve is frequency, the vertical axis of the target coherence curve is coherence degree, the target weighted coherence degree is determined according to the plurality of target coherence degrees and the target average value, and the target weighted coherence degree represents the reference coherence degree in the plurality of target coherence degrees. Determining a dividing line according to the target weighted coherence, wherein the dividing line is parallel to a transverse axis, calculating a first area value of an area above the dividing line and a second area value of an area below the dividing line in the target coherence curve, further determining a ratio between the first area value and the second area value to obtain a target ratio, determining that the connectivity of the thalamus and the cortex is in a strong connection state when the target ratio is greater than or equal to a preset ratio threshold, and determining that the connectivity of the thalamus and the cortex is in a weak connection state when the target ratio is less than the preset ratio threshold.
Optionally, the step a505, determining a target weighted coherence according to the target coherence and the target average value may include the following steps:
B501, determining the median, the maximum value and the minimum value of the plurality of target coherence degrees;
b502, determining a difference value between the maximum value and the minimum value to obtain a first difference value;
B503, determining a difference value between the median and the minimum value to obtain a second difference value;
b504, determining a difference value between the target average value and the minimum value to obtain a third difference value;
b505, determining the ratio of the second difference value to the first difference value to obtain a first ratio;
b506, determining the ratio of the third difference value to the first difference value to obtain a second ratio;
B507, normalizing the first ratio and the second ratio to obtain a first weight and a second weight; the first weight corresponds to the median; the second weight corresponds to the target average value;
and B508, carrying out weighted calculation on the median and the target average value according to the first weight and the second weight to obtain the target weighted coherence.
In a specific embodiment, a median, a maximum value and a minimum value of a plurality of target coherence degrees are determined, a difference value between the maximum value and the minimum value is determined to obtain a first difference value, a difference value between the median and the minimum value is determined to obtain a second difference value, a difference value between a target average value and the minimum value is determined to obtain a third difference value, a ratio of the second difference value to the first difference value is determined to obtain a first ratio value, a ratio of the third difference value to the first difference value is determined to obtain a second ratio value, and the first ratio value and the second ratio value are normalized respectively to obtain a first weight and a second weight, wherein the first weight corresponds to the median, and the second weight corresponds to the average value. And carrying out weighted calculation on the median and the average number according to the first weight and the second weight to obtain the target weighted coherence.
For example, under the condition that the preset first coherence degree is 0.7, the preset second coherence degree is 0.5, the target average value is 0.56 when the target coherence degree is {0.6,0.5,0.5,0.6,0.6}, the target average value is between the preset second coherence degree and the preset first coherence degree, the median of the target coherence degrees is determined to be 0.6, the maximum value is 0.6, the minimum value is 0.5, the first difference value is 0.1, the second difference value is 0.1, the third difference value is 0.06, the first ratio is 1, the second ratio is 0.6, the normalized corresponding value of the first ratio is the first weight, namely 0.625, the normalized corresponding value of the second ratio is the second weight, namely 0.375, and the weighted calculation is performed on the median and the average value according to the first weight and the second weight, so that the target weighted coherence degree is 0.585 can be obtained.
In one possible embodiment, the fitted target coherence curve is segmented based on the target weighted coherence, and a first area value of the region above the segmentation line and a second area value of the region below the segmentation line in the target coherence curve are determined, wherein a target ratio between the first area value and the second area value is greater than a preset ratio threshold, and the connectivity of the thalamus and the cortex is determined to be in a strong connection state.
In summary, by implementing the embodiment of the application, transcranial direct current stimulation can be implemented on a target detection object through a transcranial direct current stimulation device, a preset part of the target detection object is driven by an operating device to perform preset passive movement, electroencephalogram signals of the target detection object at a stimulation electrode position corresponding to thalamus and a second stimulation electrode position corresponding to cortex are respectively acquired through an electroencephalogram acquisition device, so that a first electroencephalogram signal and a second electroencephalogram signal are obtained, coherence analysis is performed on the first electroencephalogram signal and the second electroencephalogram signal, a plurality of target coherence degrees are obtained, the connectivity degree of thalamus and cortex is determined according to the plurality of target coherence degrees, and therefore the accurate detection of thalamus-cortex connectivity is realized.
Referring to fig. 6, fig. 6 is a functional unit block diagram of a thalamus-cortex connectivity detection device provided by an embodiment of the present application, where the thalamus-cortex connectivity detection device 60 is applied to a controller of a neural connectivity detection system, and the neural connectivity detection system further includes a transcranial direct current stimulation device, an operation device, and an electroencephalogram acquisition device; the operating device comprises a piston, an air compressor and a rubber band, wherein the piston is arranged in an inner hole of the air compressor, and the air compressor is used for driving the piston to move; the rubber band is used for connecting the piston and a preset part of the target detection object; the thalamus-cortex connectivity detection device 60 includes: a determination unit 601, a coherence analysis unit 602, wherein,
The determining unit 601 is configured to determine a target stimulation parameter for the target detection object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulating electrode position is the electrode position corresponding to the cerebral cortex of the target detection object;
The transcranial direct current stimulation device is used for conducting transcranial direct current stimulation on the target detection object for the duration of stimulation according to the target stimulation parameters, and the control device is used for driving the preset part of the target detection object to conduct preset passive movement;
acquiring the electroencephalogram signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively through the electroencephalogram acquisition device to obtain a first electroencephalogram signal and a second electroencephalogram signal;
The coherence analysis unit 602 is configured to perform coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees;
The determining unit 601 is further configured to determine a degree of connectivity between the thalamus and the cortex according to the plurality of target coherence degrees.
Optionally, before the transcranial direct current stimulation of the target test object for the duration of the stimulation according to the target stimulation parameter by the transcranial direct current stimulation device, the thalamus-cortex connectivity detection device 60 is further specifically configured to:
Acquiring the current stimulation intensity of the transcranial direct current stimulation device at the current moment;
the current stimulation intensity is increased to a preset stimulation intensity in a first preset time period;
continuously stimulating the target detection object according to the preset stimulation intensity within a second preset time period, and monitoring the physiological signal of the target detection object to obtain a target physiological signal; the starting time of the second preset time period is later than the ending time of the first preset time period;
when the target physiological signal meets a first preset condition, determining the preset stimulus intensity as the target stimulus intensity; the first preset condition is that the target physiological signal is in a preset safety interval;
When the target physiological signal does not meet the first preset condition, adjusting the preset stimulus intensity based on the target physiological signal to obtain standard stimulus intensity;
the standard stimulus intensity is determined as the target stimulus intensity.
Optionally, when the target physiological signal includes a heart rate, the thalamus-cortex connectivity detection device 60 is further specifically configured to, in terms of the adjusting the preset stimulus intensity based on the target physiological signal to obtain a standard stimulus intensity:
Acquiring a target safe interval and a standard heart rate value of the heart rate;
obtaining a difference value between the heart rate and the standard heart rate value to obtain a heart rate difference value;
determining an adjustment parameter of the preset stimulus intensity according to the heart rate difference value to obtain a target adjustment parameter;
When the heart rate is not in the target safety interval, determining the deviation degree of the heart rate based on the target safety interval, and obtaining a target deviation degree;
Determining a target adjustment factor corresponding to the target deviation;
And adjusting the preset stimulus intensity according to the target adjustment factor and the target adjustment parameter to obtain the standard stimulus intensity.
Optionally, in the performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees, the coherence analysis unit 602 is further specifically configured to:
respectively carrying out filtering processing on the first electroencephalogram signal and the second electroencephalogram signal based on a preset filter to obtain a third electroencephalogram signal and a fourth electroencephalogram signal;
discretizing the third electroencephalogram signal and the fourth electroencephalogram signal through a preset window function respectively to obtain a first time domain signal and a second time domain signal;
performing Fourier transform on the first time domain signal and the second time domain signal respectively to obtain a first frequency domain signal and a second frequency domain signal;
determining frequency points in the first frequency domain signal and the second frequency domain signal under the same frequency to obtain a plurality of frequency point groups; each frequency point group comprises a first frequency point set and a second frequency point set; each frequency point set comprises a plurality of frequency points; the first frequency point set is the frequency point set of the first frequency domain signal; the second frequency point set is the frequency point set of the second frequency domain signal;
And determining the coherence of the first frequency domain signal and the second frequency domain signal according to the plurality of frequency point groups to obtain a plurality of target coherence degrees.
Optionally, in the determining the coherence of the first frequency domain signal and the second frequency domain signal according to the plurality of frequency point groups, to obtain a plurality of target coherence degrees, the coherence analysis unit 602 is further specifically configured to:
determining a target frequency point group; the target frequency point group is any one of the frequency point groups;
determining the self-power spectral density of the first frequency domain signal according to a first frequency point set in the target frequency point group to obtain a first power spectral density;
determining the self-power spectral density of the second frequency domain signal according to a second frequency point set in the target frequency point set to obtain a second power spectral density;
Determining cross power spectral densities of the first frequency domain signal and the second frequency domain signal according to the target frequency point group to obtain a third power spectral density;
and determining target coherence according to the first power spectral density, the second power spectral density and the third power spectral density.
Optionally, in the aspect of determining the connectivity degree of the thalamus and the cortex according to the target coherence degrees, the determining unit 601 is further specifically configured to:
determining an average value of the plurality of target coherence degrees to obtain a target average value;
when the target average value is greater than or equal to a preset first coherence, determining that connectivity of thalamus and cortex is in a strong connection state;
when the target average value is smaller than or equal to a preset second correlation degree, determining that connectivity of thalamus and cortex is in a weak connection state;
When the target average value is larger than the preset second coherence degree and the target average value is smaller than the preset first coherence degree, fitting according to the target coherence degrees to obtain a target coherence curve; the horizontal axis of the target coherence curve is frequency and the vertical axis is coherence;
determining a target weighted coherence according to the plurality of target coherence and the target average;
determining a dividing line according to the target weighted coherence; the parting line is parallel to the transverse axis;
determining a first area value of a region above the dividing line and a second area value of a region below the dividing line in the target coherence curve;
determining a ratio between the first area value and the second area value to obtain a target ratio;
Acquiring a preset ratio threshold;
When the target ratio is greater than or equal to the preset ratio threshold, determining that connectivity of thalamus and cortex is in a strong connection state;
and when the target ratio is smaller than the preset ratio threshold, determining that the connectivity of the thalamus and the cortex is in a weak connection state.
Optionally, in the determining the target weighted coherence according to the plurality of target coherence and the target average value, the determining unit 601 is further specifically configured to:
determining a median, a maximum value, and a minimum value of the plurality of target coherence degrees;
Determining a difference value between the maximum value and the minimum value to obtain a first difference value;
determining a difference value between the median and the minimum value to obtain a second difference value;
determining a difference value between the target average value and the minimum value to obtain a third difference value;
determining the ratio of the second difference value to the first difference value to obtain a first ratio;
Determining the ratio of the third difference value to the first difference value to obtain a second ratio;
Normalizing the first ratio and the second ratio respectively to obtain a first weight and a second weight; the first weight corresponds to the median; the second weight corresponds to the target average value;
And carrying out weighted calculation on the median and the target average value according to the first weight and the second weight to obtain target weighted coherence.
The thalamus-cortex connectivity detection device 60 disclosed by the application can be used for implementing transcranial direct current stimulation on a target detection object through a transcranial direct current stimulation device, driving a preset part of the target detection object to perform preset passive movement through an operation device, respectively acquiring the electroencephalogram signals of the target detection object at the corresponding stimulation electrode position of the thalamus and the corresponding second stimulation electrode position of the cortex through an electroencephalogram acquisition device, obtaining a first electroencephalogram signal and a second electroencephalogram signal, performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal, obtaining a plurality of target coherence, determining the connectivity degree of the thalamus and the cortex according to the plurality of target coherence, and further realizing accurate detection on the thalamus-cortex connectivity.
Referring to fig. 7, fig. 7 is a schematic structural diagram of an electronic device according to an embodiment of the present application, where the electronic device may include a processor, a memory, a communication interface, and one or more programs, where the processor, the memory, and the communication interface may be connected to each other through a bus; the one or more programs are stored in the memory and configured to be executed by the processor; the electronic equipment comprises a controller of a nerve connectivity detection system, wherein the nerve connectivity detection system also comprises a transcranial direct current stimulation device, an operating device and an electroencephalogram acquisition device; the operating device comprises a piston, an air compressor and a rubber band, wherein the piston is arranged in an inner hole of the air compressor, and the air compressor is used for driving the piston to move; the rubber band is used for connecting the piston and a preset part of the target detection object; in an embodiment of the present application, the program includes instructions for performing the steps of:
Determining a target stimulation parameter for the target test object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulating electrode position is the electrode position corresponding to the cerebral cortex of the target detection object;
The transcranial direct current stimulation device is used for conducting transcranial direct current stimulation on the target detection object for the duration of stimulation according to the target stimulation parameters, and the control device is used for driving the preset part of the target detection object to conduct preset passive movement;
acquiring the electroencephalogram signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively through the electroencephalogram acquisition device to obtain a first electroencephalogram signal and a second electroencephalogram signal;
Performing coherence analysis on the first electroencephalogram signal and the second electroencephalogram signal to obtain a plurality of target coherence degrees;
Determining the degree of connectivity of the thalamus with the cortex based on the plurality of target coherence degrees.
According to the electronic equipment, transcranial direct current stimulation can be carried out on a target detection object through the transcranial direct current stimulation device, the control device drives the preset part of the target detection object to perform preset passive movement, the electroencephalogram acquisition device respectively acquires the electroencephalogram signals of the target detection object at the position of the stimulation electrode corresponding to the thalamus and the position of the second stimulation electrode corresponding to the cortex, so that the first electroencephalogram signals and the second electroencephalogram signals are obtained, coherence analysis is carried out on the first electroencephalogram signals and the second electroencephalogram signals, a plurality of target coherence degrees are obtained, the connectivity degree of the thalamus and the cortex is determined according to the plurality of target coherence degrees, and therefore the thalamus-cortex connectivity is accurately detected.
The embodiment of the application also provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program for electronic data exchange, and the computer program causes a computer to execute part or all of the steps of any one of the above method embodiments, and the computer includes an electronic device.
Embodiments of the present application also provide a computer program product comprising a non-transitory computer-readable storage medium storing a computer program operable to cause a computer to perform part or all of the steps of any one of the methods described in the method embodiments above. The computer program product may be a software installation package, said computer comprising an electronic device.
Those of ordinary skill in the art will appreciate that implementing all or part of the above-described method embodiments may be accomplished by a computer program to instruct related hardware, the program may be stored in a computer readable storage medium, and the program may include the above-described method embodiments when executed. And the aforementioned storage medium includes: ROM or random access memory RAM, magnetic or optical disk, etc.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied in hardware, or may be embodied in software instructions executed by a processor. The software instructions may be comprised of corresponding software modules that may be stored in RAM, flash memory, ROM, EPROM, electrically Erasable EPROM (EEPROM), registers, hard disk, a removable disk, a compact disk read-only (CD-ROM), or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC. In addition, the ASIC may be located in a terminal device or a management device. The processor and the storage medium may reside as discrete components in a terminal device or management device.
Those skilled in the art will appreciate that in one or more of the examples described above, the functions described in the embodiments of the present application may be implemented, in whole or in part, in software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present application, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital subscriber line (digital subscriber line, DSL)), or wireless (e.g., infrared, wireless, microwave, etc.). The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a magnetic tape), an optical medium (e.g., a digital video disc (digital video disc, DVD)), or a semiconductor medium (e.g., a Solid State Drive (SSD)), or the like.
The respective apparatuses and the respective modules/units included in the products described in the above embodiments may be software modules/units, may be hardware modules/units, or may be partly software modules/units, and partly hardware modules/units. For example, for each device or product applied to or integrated on a chip, each module/unit included in the device or product may be implemented in hardware such as a circuit, or at least some modules/units may be implemented in software program, where the software program runs on a processor integrated inside the chip, and the remaining (if any) part of modules/units may be implemented in hardware such as a circuit; for each device and product applied to or integrated in the chip module, each module/unit contained in the device and product can be realized in a hardware manner such as a circuit, different modules/units can be located in the same component (such as a chip, a circuit module and the like) or different components of the chip module, or at least part of the modules/units can be realized in a software program, the software program runs on a processor integrated in the chip module, and the rest (if any) of the modules/units can be realized in a hardware manner such as a circuit; for each device, product, or application to or integrated with the terminal device, each module/unit included in the device may be implemented in hardware such as a circuit, and different modules/units may be located in the same component (e.g., a chip, a circuit module, etc.) or different components in the terminal device, or at least some modules/units may be implemented in a software program, where the software program runs on a processor integrated within the terminal device, and the remaining (if any) some modules/units may be implemented in hardware such as a circuit.
The foregoing detailed description of the embodiments of the present application further illustrates the purposes, technical solutions and advantageous effects of the embodiments of the present application, and it should be understood that the foregoing description is only a specific implementation of the embodiments of the present application, and is not intended to limit the scope of the embodiments of the present application, and any modifications, equivalent substitutions, improvements, etc. made on the basis of the technical solutions of the embodiments of the present application should be included in the scope of the embodiments of the present application.

Claims (10)

1.一种丘脑-皮层连接性检测方法,其特征在于,应用于神经连接性检测系统的控制器,所述神经连接性检测系统还包括经颅直流电刺激装置、操纵装置、脑电采集装置;所述操纵装置包括活塞、空气压缩机、橡皮筋,所述活塞设置于所述空气压缩机的内部孔洞中,所述空气压缩机用于带动所述活塞进行移动;所述橡皮筋用于连接所述活塞和目标检测对象的预设部位;所述方法包括:1. A thalamus-cortex connectivity detection method, characterized in that the controller is applied to a neural connectivity detection system, wherein the neural connectivity detection system further comprises a transcranial direct current stimulation device, a manipulation device, and an electroencephalogram acquisition device; the manipulation device comprises a piston, an air compressor, and a rubber band, wherein the piston is disposed in an internal hole of the air compressor, and the air compressor is used to drive the piston to move; the rubber band is used to connect the piston and a preset position of a target detection object; the method comprises: 确定针对所述目标检测对象的目标刺激参数;所述目标刺激参数包括:目标刺激强度、第一刺激电极位置、第二刺激电极位置、刺激持续时长;所述第一刺激电极位置为所述目标检测对象的丘脑对应的电极位置;所述第二刺激电极位置为所述目标检测对象的大脑皮层对应的电极位置;Determine target stimulation parameters for the target detection object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulation electrode position is the electrode position corresponding to the cerebral cortex of the target detection object; 通过所述经颅直流电刺激装置根据所述目标刺激参数对所述目标检测对象实施所述刺激持续时长的经颅直流电刺激,并通过所述操纵装置带动所述目标检测对象的所述预设部位进行预设被动运动;Implementing transcranial direct current stimulation of the stimulation duration on the target detection object by the transcranial direct current stimulation device according to the target stimulation parameters, and driving the preset part of the target detection object to perform preset passive movement by the manipulation device; 通过所述脑电采集装置分别采集所述目标检测对象在所述第一刺激电极位置和所述第二刺激电极位置的脑电信号,得到第一脑电信号和第二脑电信号;The EEG acquisition device acquires EEG signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively to obtain a first EEG signal and a second EEG signal; 对所述第一脑电信号和所述第二脑电信号进行相干性分析,得到多个目标相干度;Performing coherence analysis on the first EEG signal and the second EEG signal to obtain multiple target coherences; 根据所述多个目标相干度确定丘脑与皮层的连接性程度。The degree of connectivity between the thalamus and the cortex is determined based on the multiple target coherences. 2.如权利要求1所述的方法,其特征在于,在所述通过所述经颅直流电刺激装置根据所述目标刺激参数对所述目标检测对象实施所述刺激持续时长的经颅直流电刺激之前,所述方法包括:2. The method according to claim 1, characterized in that before the transcranial direct current stimulation device performs the transcranial direct current stimulation of the stimulation duration on the target detection object according to the target stimulation parameters, the method comprises: 获取所述经颅直流电刺激装置当前时刻的当前刺激强度;Obtaining the current stimulation intensity of the transcranial direct current stimulation device at the current moment; 在第一预设时间段内将所述当前刺激强度提升至预设刺激强度;Increasing the current stimulation intensity to a preset stimulation intensity within a first preset time period; 在第二预设时间段内根据所述预设刺激强度持续对所述目标检测对象进行刺激,监测所述目标检测对象的生理信号,得到目标生理信号;所述第二预设时间段的开始时间晚于所述第一预设时间段的结束时间;Continuously stimulating the target detection object according to the preset stimulation intensity within a second preset time period, monitoring the physiological signal of the target detection object, and obtaining a target physiological signal; the start time of the second preset time period is later than the end time of the first preset time period; 在所述目标生理信号满足第一预设条件时,则将所述预设刺激强度确定为所述目标刺激强度;所述第一预设条件为所述目标生理信号处于预设安全区间;When the target physiological signal satisfies a first preset condition, the preset stimulation intensity is determined as the target stimulation intensity; the first preset condition is that the target physiological signal is within a preset safety range; 在所述目标生理信号不满足所述第一预设条件时,基于所述目标生理信号对所述预设刺激强度进行调整,得到标准刺激强度;When the target physiological signal does not meet the first preset condition, adjusting the preset stimulation intensity based on the target physiological signal to obtain a standard stimulation intensity; 将所述标准刺激强度确定为所述目标刺激强度。The standard stimulation intensity is determined as the target stimulation intensity. 3.如权利要求2所述的方法,其特征在于,在所述目标生理信号包括心率时,所述基于所述目标生理信号对所述预设刺激强度进行调整,得到标准刺激强度,包括:3. The method according to claim 2, wherein when the target physiological signal includes a heart rate, adjusting the preset stimulation intensity based on the target physiological signal to obtain a standard stimulation intensity comprises: 获取所述心率的目标安全区间和标准心率值;Obtaining a target safety interval and a standard heart rate value of the heart rate; 获取所述心率与所述标准心率值的差值,得到心率差值;Obtaining a difference between the heart rate and the standard heart rate value to obtain a heart rate difference; 根据所述心率差值确定所述预设刺激强度的调整参数,得到目标调整参数;Determine the adjustment parameter of the preset stimulation intensity according to the heart rate difference to obtain a target adjustment parameter; 在所述心率不处于所述目标安全区间时,确定所述心率基于所述目标安全区间的偏离度,得到目标偏离度;When the heart rate is not within the target safety interval, determining a deviation of the heart rate based on the target safety interval to obtain a target deviation; 确定所述目标偏离度对应的目标调整因子;Determining a target adjustment factor corresponding to the target deviation; 根据所述目标调整因子和所述目标调整参数对所述预设刺激强度进行调整,得到所述标准刺激强度。The preset stimulation intensity is adjusted according to the target adjustment factor and the target adjustment parameter to obtain the standard stimulation intensity. 4.如权利要求1-3任一项所述的方法,其特征在于,所述对所述第一脑电信号和所述第二脑电信号进行相干性分析,得到多个目标相干度,包括:4. The method according to any one of claims 1 to 3, characterized in that the coherence analysis of the first EEG signal and the second EEG signal to obtain multiple target coherences comprises: 基于预设滤波器分别对所述第一脑电信号和所述第二脑电信号进行滤波处理,得到第三脑电信号和第四脑电信号;Based on a preset filter, the first EEG signal and the second EEG signal are filtered to obtain a third EEG signal and a fourth EEG signal; 通过预设窗函数分别对所述第三脑电信号和所述第四脑电信号进行离散化,得到第一时域信号和第二时域信号;Discretizing the third EEG signal and the fourth EEG signal respectively by using a preset window function to obtain a first time domain signal and a second time domain signal; 将所述第一时域信号和所述第二时域信号分别进行傅里叶变换,得到第一频域信号和第二频域信号;Performing Fourier transform on the first time domain signal and the second time domain signal respectively to obtain a first frequency domain signal and a second frequency domain signal; 确定所述第一频域信号和所述第二频域信号中在同一频率下的频率点,得到多个频率点组;每一频率点组包括第一频率点集和第二频率点集;每一频率点集包括多个频率点;所述第一频率点集为所述第一频域信号的频率点集;所述第二频率点集为所述第二频域信号的频率点集;Determine frequency points at the same frequency in the first frequency domain signal and the second frequency domain signal to obtain multiple frequency point groups; each frequency point group includes a first frequency point set and a second frequency point set; each frequency point set includes multiple frequency points; the first frequency point set is the frequency point set of the first frequency domain signal; the second frequency point set is the frequency point set of the second frequency domain signal; 根据所述多个频率点组确定所述第一频域信号和所述第二频域信号的相干性,得到多个目标相干度。The coherence of the first frequency domain signal and the second frequency domain signal is determined according to the multiple frequency point groups to obtain multiple target coherence degrees. 5.如权利要求4所述的方法,其特征在于,所述根据所述多个频率点组确定所述第一频域信号和所述第二频域信号的相干性,得到多个目标相干度,包括:5. The method according to claim 4, wherein determining the coherence of the first frequency domain signal and the second frequency domain signal according to the multiple frequency point groups to obtain multiple target coherences comprises: 确定目标频率点组;所述目标频率点组为所述多个频率点组中任一频率点组;Determine a target frequency point group; the target frequency point group is any frequency point group among the multiple frequency point groups; 根据所述目标频率点组中的第一频率点集确定所述第一频域信号的自功率谱密度,得到第一功率谱密度;Determine the autopower spectral density of the first frequency domain signal according to the first frequency point set in the target frequency point group to obtain a first power spectral density; 根据所述目标频率点组中的第二频率点集确定所述第二频域信号的自功率谱密度,得到第二功率谱密度;Determine the autopower spectral density of the second frequency domain signal according to the second frequency point set in the target frequency point group to obtain a second power spectral density; 根据所述目标频率点组确定所述第一频域信号和所述第二频域信号的互功率谱密度,得到第三功率谱密度;Determine the cross power spectrum density of the first frequency domain signal and the second frequency domain signal according to the target frequency point group to obtain a third power spectrum density; 根据所述第一功率谱密度、所述第二功率谱密度和所述第三功率谱密度确定目标相干度。A target coherence degree is determined according to the first power spectral density, the second power spectral density, and the third power spectral density. 6.如权利要求5所述的方法,其特征在于,所述根据所述多个目标相干度确定丘脑与皮层的连接性程度,包括:6. The method according to claim 5, wherein determining the degree of connectivity between the thalamus and the cortex according to the multiple target coherences comprises: 确定所述多个目标相干度的平均值,得到目标平均值;Determining an average value of the multiple target coherences to obtain a target average value; 在所述目标平均值大于或等于预设第一相干度时,确定丘脑与皮层的连接性处于强连接状态;When the target average value is greater than or equal to a preset first coherence degree, determining that the connectivity between the thalamus and the cortex is in a strong connection state; 在所述目标平均值小于或等于预设第二相干度时,确定丘脑与皮层的连接性处于弱连接状态;When the target average value is less than or equal to a preset second coherence degree, determining that the connectivity between the thalamus and the cortex is in a weak connection state; 在所述目标平均值大于所述预设第二相干度且所述目标平均值小于所述预设第一相干度时,根据所述多个目标相干度进行拟合,得到目标相干曲线;所述目标相干曲线的横轴为频率且纵轴为相干度;When the target average value is greater than the preset second coherence and the target average value is less than the preset first coherence, fitting is performed according to the multiple target coherences to obtain a target coherence curve; the horizontal axis of the target coherence curve is frequency and the vertical axis is coherence; 根据所述多个目标相干度和所述目标平均值确定目标加权相干度;Determining a target weighted coherence according to the multiple target coherences and the target average value; 根据所述目标加权相干度确定分割线;所述分割线平行于所述横轴;Determining a segmentation line according to the target weighted coherence; the segmentation line is parallel to the horizontal axis; 确定所述目标相干曲线中的所述分割线以上区域的第一面积值,以及所述分割线以下区域的第二面积值;Determine a first area value of a region above the dividing line and a second area value of a region below the dividing line in the target coherence curve; 确定所述第一面积值与所述第二面积值之间的比值,得到目标比值;Determine a ratio between the first area value and the second area value to obtain a target ratio; 获取预设比值阈值;Obtaining a preset ratio threshold; 在所述目标比值大于或等于所述预设比值阈值时,确定丘脑与皮层的连接性处于强连接状态;When the target ratio is greater than or equal to the preset ratio threshold, determining that the connectivity between the thalamus and the cortex is in a strong connection state; 在所述目标比值小于所述预设比值阈值时,确定丘脑与皮层的连接性处于弱连接状态。When the target ratio is less than the preset ratio threshold, it is determined that the connectivity between the thalamus and the cortex is in a weak connection state. 7.如权利要求6所述的方法,其特征在于,所述根据所述多个目标相干度和所述目标平均值确定目标加权相干度,包括:7. The method according to claim 6, wherein determining the target weighted coherence according to the plurality of target coherences and the target average value comprises: 确定所述多个目标相干度的中位数、最大值和最小值;Determining the median, maximum, and minimum values of the plurality of target coherences; 确定所述最大值和所述最小值之间的差值,得到第一差值;Determine a difference between the maximum value and the minimum value to obtain a first difference; 确定所述中位数与所述最小值的差值,得到第二差值;Determine the difference between the median and the minimum value to obtain a second difference; 确定所述目标平均值与所述最小值的差值,得到第三差值;Determine a difference between the target average value and the minimum value to obtain a third difference value; 确定所述第二差值与所述第一差值的比值,得到第一比值;determining a ratio of the second difference to the first difference to obtain a first ratio; 确定所述第三差值与所述第一差值的比值,得到第二比值;Determine a ratio of the third difference to the first difference to obtain a second ratio; 分别对所述第一比值和所述第二比值进行归一化,得到第一权重和第二权重;所述第一权重与所述中位数对应;所述第二权重与所述目标平均值对应;Normalizing the first ratio and the second ratio respectively to obtain a first weight and a second weight; the first weight corresponds to the median; and the second weight corresponds to the target average; 根据所述第一权重和所述第二权重对所述中位数和所述目标平均值进行加权计算,得到目标加权相干度。The median and the target average are weightedly calculated according to the first weight and the second weight to obtain a target weighted coherence. 8.一种丘脑-皮层连接性检测装置,其特征在于,应用于神经连接性检测系统的控制器,所述神经连接性检测系统还包括经颅直流电刺激装置、操纵装置、脑电采集装置;所述操纵装置包括活塞、空气压缩机、橡皮筋,所述活塞设置于所述空气压缩机的内部孔洞中,所述空气压缩机用于带动所述活塞进行移动;所述橡皮筋用于连接所述活塞和目标检测对象的预设部位;所述丘脑-皮层连接性检测装置包括:确定单元、相干性分析单元,其中,8. A thalamus-cortex connectivity detection device, characterized in that it is a controller applied to a neural connectivity detection system, wherein the neural connectivity detection system further comprises a transcranial direct current stimulation device, a manipulation device, and an electroencephalogram acquisition device; the manipulation device comprises a piston, an air compressor, and a rubber band, wherein the piston is arranged in an internal hole of the air compressor, and the air compressor is used to drive the piston to move; the rubber band is used to connect the piston and a preset position of a target detection object; the thalamus-cortex connectivity detection device comprises: a determination unit, a coherence analysis unit, wherein, 所述确定单元,用于确定针对所述目标检测对象的目标刺激参数;所述目标刺激参数包括:目标刺激强度、第一刺激电极位置、第二刺激电极位置、刺激持续时长;所述第一刺激电极位置为所述目标检测对象的丘脑对应的电极位置;所述第二刺激电极位置为所述目标检测对象的大脑皮层对应的电极位置;The determination unit is used to determine target stimulation parameters for the target detection object; the target stimulation parameters include: target stimulation intensity, first stimulation electrode position, second stimulation electrode position, stimulation duration; the first stimulation electrode position is the electrode position corresponding to the thalamus of the target detection object; the second stimulation electrode position is the electrode position corresponding to the cerebral cortex of the target detection object; 通过所述经颅直流电刺激装置根据所述目标刺激参数对所述目标检测对象实施所述刺激持续时长的经颅直流电刺激,并通过所述操纵装置带动所述目标检测对象的所述预设部位进行预设被动运动;Implementing transcranial direct current stimulation of the stimulation duration on the target detection object by the transcranial direct current stimulation device according to the target stimulation parameters, and driving the preset part of the target detection object to perform preset passive movement by the manipulation device; 通过所述脑电采集装置分别采集所述目标检测对象在所述第一刺激电极位置和所述第二刺激电极位置的脑电信号,得到第一脑电信号和第二脑电信号;The EEG acquisition device acquires EEG signals of the target detection object at the first stimulation electrode position and the second stimulation electrode position respectively to obtain a first EEG signal and a second EEG signal; 所述相干性分析单元,用于对所述第一脑电信号和所述第二脑电信号进行相干性分析,得到多个目标相干度;The coherence analysis unit is used to perform coherence analysis on the first EEG signal and the second EEG signal to obtain multiple target coherences; 所述确定单元,还用于根据所述多个目标相干度确定丘脑与皮层的连接性程度。The determination unit is further used to determine the degree of connectivity between the thalamus and the cortex based on the multiple target coherences. 9.一种电子设备,其特征在于,包括:处理器、存储器、通信接口以及一个或多个程序;所述一个或多个程序被存储在所述存储器中,并且被配置成由所述处理器执行,所述程序包括用于执行如权利要求1-7任一项所述的方法中的步骤的指令。9. An electronic device, characterized in that it comprises: a processor, a memory, a communication interface and one or more programs; the one or more programs are stored in the memory and are configured to be executed by the processor, and the program includes instructions for executing the steps in the method as described in any one of claims 1-7. 10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被处理器执行时使所述处理器执行如权利要求1-7任一项所述的方法。10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program, wherein the computer program comprises program instructions, and when the program instructions are executed by a processor, the processor is caused to execute the method according to any one of claims 1 to 7.
CN202411448546.3A 2024-10-17 2024-10-17 A thalamus-cortex connectivity detection method and related device Active CN118948226B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411448546.3A CN118948226B (en) 2024-10-17 2024-10-17 A thalamus-cortex connectivity detection method and related device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411448546.3A CN118948226B (en) 2024-10-17 2024-10-17 A thalamus-cortex connectivity detection method and related device

Publications (2)

Publication Number Publication Date
CN118948226A true CN118948226A (en) 2024-11-15
CN118948226B CN118948226B (en) 2025-02-25

Family

ID=93400067

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411448546.3A Active CN118948226B (en) 2024-10-17 2024-10-17 A thalamus-cortex connectivity detection method and related device

Country Status (1)

Country Link
CN (1) CN118948226B (en)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080009458A (en) * 2006-07-24 2008-01-29 중앙대학교 산학협력단 Emotion Recognition Device Using Neural Network
US20080103547A1 (en) * 2004-09-21 2008-05-01 University Of Florida Research Foundation, Inc. Multiple lead method for deep brain stimulation
WO2009023162A1 (en) * 2007-08-11 2009-02-19 The Johns Hopkins University Focal noninvasive stimulation of the sensory cortex of a subject with cerebral palsy
US20120116475A1 (en) * 2010-11-10 2012-05-10 Medtronic, Inc. Arousal state modulation with electrical stimulation
CN102488514A (en) * 2011-12-09 2012-06-13 天津大学 Method for analyzing relativity between electroencephalograph and myoelectricity based on autonomous and stimulation movement modalities
US20150133812A1 (en) * 2012-01-09 2015-05-14 Richard Christopher DeCharms Methods and systems for quantitative measurement of mental states
KR20160026889A (en) * 2013-06-28 2016-03-09 코치 유니버시티 An electro-stimulation device
WO2020049514A1 (en) * 2018-09-06 2020-03-12 Alpha Omega Neuro Technologies Ltd. Therapeutic space assessment
CN111494799A (en) * 2020-03-20 2020-08-07 上海脑晟智能科技有限公司 Intelligent nerve pressure reduction system
CN113662564A (en) * 2021-09-24 2021-11-19 电子科技大学 An event-evoked potential extraction device based on eigenmode function coherence analysis
US20220379132A1 (en) * 2021-05-25 2022-12-01 Wave Neuroscience, Inc. Brain stimulation apparatus, method, and medium
WO2023201863A1 (en) * 2022-04-18 2023-10-26 北京理工大学 Method and system for closed-loop regulation of neural network by adaptive alternating current stimulation
CN117084693A (en) * 2023-10-07 2023-11-21 吉林大学 An encephalopathy detection system based on multi-domain feature fusion and attention mechanism
US20240017084A1 (en) * 2022-07-15 2024-01-18 The Florida State University Research Foundation, Inc. fNIRS for Dosing Transcranial Magnetic Brain Stimulation

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080103547A1 (en) * 2004-09-21 2008-05-01 University Of Florida Research Foundation, Inc. Multiple lead method for deep brain stimulation
KR20080009458A (en) * 2006-07-24 2008-01-29 중앙대학교 산학협력단 Emotion Recognition Device Using Neural Network
WO2009023162A1 (en) * 2007-08-11 2009-02-19 The Johns Hopkins University Focal noninvasive stimulation of the sensory cortex of a subject with cerebral palsy
US20120116475A1 (en) * 2010-11-10 2012-05-10 Medtronic, Inc. Arousal state modulation with electrical stimulation
CN102488514A (en) * 2011-12-09 2012-06-13 天津大学 Method for analyzing relativity between electroencephalograph and myoelectricity based on autonomous and stimulation movement modalities
US20150133812A1 (en) * 2012-01-09 2015-05-14 Richard Christopher DeCharms Methods and systems for quantitative measurement of mental states
KR20160026889A (en) * 2013-06-28 2016-03-09 코치 유니버시티 An electro-stimulation device
WO2020049514A1 (en) * 2018-09-06 2020-03-12 Alpha Omega Neuro Technologies Ltd. Therapeutic space assessment
CN111494799A (en) * 2020-03-20 2020-08-07 上海脑晟智能科技有限公司 Intelligent nerve pressure reduction system
US20220379132A1 (en) * 2021-05-25 2022-12-01 Wave Neuroscience, Inc. Brain stimulation apparatus, method, and medium
CN113662564A (en) * 2021-09-24 2021-11-19 电子科技大学 An event-evoked potential extraction device based on eigenmode function coherence analysis
WO2023201863A1 (en) * 2022-04-18 2023-10-26 北京理工大学 Method and system for closed-loop regulation of neural network by adaptive alternating current stimulation
US20240017084A1 (en) * 2022-07-15 2024-01-18 The Florida State University Research Foundation, Inc. fNIRS for Dosing Transcranial Magnetic Brain Stimulation
CN117084693A (en) * 2023-10-07 2023-11-21 吉林大学 An encephalopathy detection system based on multi-domain feature fusion and attention mechanism

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
胡琰茹;刘晓莉;乔德才;: "一次性力竭运动过程中大鼠"黑质-丘脑-皮层"通路神经元相干性及递质动态变化", 中国应用生理学杂志, no. 03, 28 May 2017 (2017-05-28), pages 14 - 21 *
靳静娜;金芳;王欣;刘志朋;殷涛;: "低频重复经颅磁刺激刺激初级运动皮层对大脑功能连接影响的研究", 生物医学工程学杂志, no. 04, 25 August 2017 (2017-08-25), pages 15 - 21 *

Also Published As

Publication number Publication date
CN118948226B (en) 2025-02-25

Similar Documents

Publication Publication Date Title
Chen et al. Adaptive real-time identification of motor unit discharges from non-stationary high-density surface electromyographic signals
US20210267465A1 (en) Systems and methods for detecting strokes
EP3509694B1 (en) System and method for generating electromagnetic treatment protocols
CN117503163B (en) Method and device for detecting abnormal discharge of human brain, storage medium and electronic equipment
CN113577559B (en) Closed-loop deep brain stimulation device, system and equipment based on multiple signals
WO2016075324A1 (en) Scoring method based on improved signals analysis
CN104571504B (en) A kind of online brain-machine interface method based on Imaginary Movement
CN117462146A (en) Human brain abnormal discharge detection method, device, storage medium and electronic equipment
CN112826510A (en) Electrophysiological Signal Acquisition System
CN115192907A (en) Real-time biofeedback percutaneous vagus nerve electronic acupuncture device
Angrisani et al. Instrumentation and measurements for non-invasive EEG-based brain-computer interface
CN109284009B (en) System and method for improving auditory steady-state response brain-computer interface performance
CN118416390A (en) Sleep control method, device, system and equipment
CN117839068A (en) An automatic electrical stimulation system with adjustable phase based on real-time EEG signals
CN116474264A (en) Myoelectric signal closed-loop control method and device based on functional electrical stimulation
CN118948226B (en) A thalamus-cortex connectivity detection method and related device
CN117339103A (en) Stomach function regulating system based on stomach rhythm and stomach brain coupling
CN115177272A (en) System and method for monitoring spinal cord and nerve in intervertebral foramen mirror operation
Ab Rahman et al. Classification of EEG signal for body earthing application
Atchuthan et al. Classification of noxious and non-noxious event-related potentials from S1 in pigs using a convolutional neural network
Soni et al. Automated Epilepsy Detection System Based On Tertiary Wavelet Model (TWM) Techniques
CN120571158B (en) Acupoint stimulation current decision method and system based on integrated fingerstall electrode
TWI846436B (en) System and method for assessing movement disorders in parkinson's disease using brain activity
CN114288634B (en) Fitness action self-recognition and alarm system based on EMG signal collection
Ghosh et al. Real-Time Brain Signals Monitoring System Using Single-Channel EEG for BCI Applications

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant