阅读说明：本技术 基于相关性的运动想象的导联选择方法及装置 (Correlation-based motor imagery lead selection method and device ) 是由 金晶 孙浩 于 2021-09-26 设计创作，主要内容包括：本发明涉及一种基于相关性的运动想象的导联选择方法及装置。该运动想象的导联选择方法包括：步骤S1,预处理；步骤S2,数据划分,将用于记录步骤S1所获得的EEG信号的全部导联划分为目标参考导联集合和备选参考导联集合,所述目标参考导联集合包含电极位置C3、C4和Cz的目标导联,所述备选参考导联集合包含非电极位置C3、C4和Cz的备选导联；步骤S3,计算相关性；步骤S4,计算权重；步骤S5,排序选择。通过本发明提供的方案有助于选择出少量有效导联来提升运动想象任务的识别效率。(The invention relates to a correlation-based motor imagery lead selection method and device. The motor imagery lead selection method comprises the following steps: step S1, preprocessing; a step S2 of data division, dividing all leads used for recording the EEG signals obtained in step S1 into a target reference lead set comprising target leads for electrode positions C3, C4 and Cz and an alternative reference lead set comprising alternative leads for non-electrode positions C3, C4 and Cz; step S3, calculating correlation; step S4, calculating weight; step S5, sort the selections. The scheme provided by the invention is beneficial to selecting a small number of effective leads to improve the identification efficiency of a motor imagery task.)

1. A correlation-based motor imagery lead selection method, comprising:

step S1, preprocessing, namely filtering the EEG signals which are originally acquired by using a Butterworth filter to obtain EEG signals of 0.5-30 Hz;

a step S2 of data division, dividing all leads used for recording the EEG signals obtained in step S1 into a target reference lead set comprising target leads for electrode positions C3, C4 and Cz and an alternative reference lead set comprising alternative leads for non-electrode positions C3, C4 and Cz;

step S3, calculating the correlation, traversing each alternative lead, calculating the correlation between each alternative lead and each target lead, and obtaining the correlation between each alternative lead and each target lead;

step S4, calculating weight, calculating corresponding weight according to the correlation, and calculating combined correlation index according to the weight;

and step S5, selecting the alternative leads according to the combined relevance index, selecting the alternative leads with the preset number in the front order, and combining the alternative leads with the target lead to form a target lead set.

2. A lead selection method as in claim 1, further comprising:

step S6, extracting features, namely extracting the features of the target lead set by adopting a common space mode algorithm;

and step S7, classifying, namely classifying by adopting a support vector machine.

3. A lead selection method as claimed in claim 1, wherein in step S3, the target reference lead set Y ═ { Y ═ is defined₁,y₂,y₃}，y₁,y₂,y₃Corresponding to electrode positions C3, C4 and Cz, respectively, each y_iDefining said set of alternative reference leads X ═ X, representing the EEG signal recorded for each lead, with dimensions 1T and T being the number of sample points₁,x₂,…,x_nDimension 1T, T is the number of sampling points, for each of said alternative leads x_iCalculating it from the target lead y_iThe correlation between them;

defining the correlation asThe calculation formula for each component is:

y_prepresenting the p-th target reference lead variable, I representing the mutual information calculation, for each pair of leads (x) for which a mutual information calculation is required_i,y_p) Wherein the sampling point recorded by each lead is T, the formula of mutual information calculation is:

in the formula 2, the first and second groups,is the current x_iAnd y_pThe joint probability density distribution function of (a),is x_iThe edge probability density function of (1).

4. A lead selection method as claimed in claim 3 in which in step S4 weights are calculated for the correlations to obtain a weighted combined mutual information δ_iThe calculation formula is as follows:

where M represents the number of electrode sites, M-3,representing alternative leads x_iCorrelation with the mth target reference lead.

5. A lead selection method according to claim 1, wherein in step S5, if the overall lead number is less than a set value, the set number is 50% of the overall lead number, and if the overall lead number is equal to or greater than the set value, the set number is 23% to 27% of the overall lead number; the overall number of leads is the sum of the target set of reference leads and the alternative reference leads.

6. A lead selection method as claimed in claim 1 in which in step S1, the originally acquired EEG signal is filtered using a butterworth filter of order five.

7. A correlation-based motor imagery lead selection device, comprising

The preprocessing module is suitable for filtering the EEG signals which are originally acquired by using a Butterworth filter to obtain EEG signals of 0.5-30 Hz;

a data partitioning module adapted to partition all leads of EEG signals obtained after processing by the pre-processing module into a target reference lead set comprising target leads of electrode positions C3, C4 and Cz and a set of alternative reference leads comprising alternative leads of non-electrode positions C3, C4 and Cz;

the correlation calculation module is suitable for traversing each alternative lead and calculating the correlation between the alternative lead and each target lead respectively to obtain the correlation between each alternative lead and each target lead;

the weight calculation module is suitable for calculating corresponding weight according to the correlation and calculating a combined correlation index according to the weight;

a sorting module adapted to sort the alternative leads according to the combined relevance indicator, select a set number of the alternative leads sorted top, and combine them with the target lead to form a target lead set.

8. A lead selection device as claimed in claim 7, further comprising an evaluation module adapted to apply a co-spatial mode algorithm to feature extraction and classification of the target lead set using a support vector machine.

9. A computer readable medium having stored thereon computer instructions, wherein the computer instructions when executed perform the steps of the lead selection method of any one of claims 1 to 6.

10. A lead selection system comprising a memory and a processor, the memory having stored thereon computer instructions executable on the processor, wherein the processor, when executing the computer instructions, performs the steps of the lead selection method of any one of claims 1 to 6.

Technical Field

The invention relates to the field of brain-Computer Interface (BCI for short) and medical rehabilitation, in particular to a method and a device for selecting leads of electroencephalograms (EEG) based on motor imagery.

Background

The brain-computer interface technology is a direct path established without depending on conventional muscles and nerves, and has been applied to various aspects of medical treatment, military, assistance for disabled persons, life entertainment, and the like. The brain-computer interface system identifies the thinking state or the movement intention of the subject by analyzing the acquired electroencephalogram information, and further converts the thinking state or the movement intention into a control command for controlling the external equipment. The motor imagery brain-computer interface system is a spontaneous brain-computer interface system which does not need to rely on external stimulation, and a user can autonomously decide the starting and the ending of a motor imagery task and flexibly control external equipment.

In order to comprehensively record the brain electrical information, the conventional brain-computer interface system usually adopts a relatively large-scale dense lead distribution form to acquire EEG signals, but the acquisition range is relatively wide, and possibly excessive artifact information and noise information are recorded at the same time, and the artifact information and the noise information will have adverse effects on the analysis of the signals, so that the performance of the motor imagery brain-computer interface system is influenced. Leads containing useful information can be selected by conducting lead selection, redundant leads are effectively eliminated, and further the influence of artifact and noise information on the performance of the motor imagery brain-computer interface system can be avoided. Lead selection algorithms can be divided into two categories: 1) filtering selection, wherein the method is mainly characterized in that a target lead combination is rapidly selected according to an evaluation index by constructing the evaluation index of a lead, and the method has the main defect that an accurate lead combination cannot be obtained; 2) the classifier method mainly determines the optimal lead combination according to the classification accuracy rate obtained by the classifier by traversing all possible lead combinations, and is longer in required calculation time and not applicable under the condition of larger number of leads.

Disclosure of Invention

The invention aims to provide a correlation-based motor imagery lead selection method and device, which can improve the identification efficiency of a motor imagery task by selecting a small number of effective leads.

In order to solve the technical problem, the invention provides a correlation-based motor imagery lead selection method, which comprises the following steps:

step S1, preprocessing, namely filtering the EEG signals which are originally acquired by using a Butterworth filter to obtain EEG signals of 0.5-30 Hz;

a step S2 of data division, dividing all leads used for recording the EEG signals obtained in step S1 into a target reference lead set comprising target leads for electrode positions C3, C4 and Cz and an alternative reference lead set comprising alternative leads for non-electrode positions C3, C4 and Cz;

step S3, calculating the correlation, traversing each alternative lead, calculating the correlation between each alternative lead and each target lead, and obtaining the correlation between each alternative lead and each target lead;

step S4, calculating weight, calculating corresponding weight according to the correlation, and calculating combined correlation index according to the weight;

and step S5, selecting the alternative leads according to the combined relevance index, selecting the alternative leads with the preset number in the front order, and combining the alternative leads with the target lead to form a target lead set.

According to an embodiment of the invention, the lead selection method further comprises:

step S6, extracting features, namely extracting the features of the target lead set by adopting a common space mode algorithm;

and step S7, classifying, namely classifying by adopting a support vector machine.

According to an embodiment of the present invention, in step S3, the target reference lead set Y ═ { Y ═ is defined₁,y₂,y₃}，y₁,y₂,y₃Corresponding to electrode positions C3, C4 and Cz, respectively, each y_iDefining said set of alternative reference leads X = { X } representing EEG signals recorded for each lead with dimensions 1T, T being the number of sample points₁,x₂,…,x_nDimension 1T, T is the number of sampling points, for each of said alternative leads x_iCalculating it from the target lead y_iThe correlation between them;

defining the correlation asEach one of which isThe calculation formula of the components is as follows:

y_prepresenting the p-th target reference lead variable, I representing the mutual information calculation, for each pair of leads (x) for which a mutual information calculation is required_i,y_p) Wherein the sampling point recorded by each lead is T, the formula of mutual information calculation is:

in the formula 2, the first and second groups,is the current x_iAnd y_pThe joint probability density distribution function of (a),is x_iThe edge probability density function of (1).

According to an embodiment of the present invention, in step S4, weights are calculated for the correlations to obtain a weighted combined mutual information δ_iThe calculation formula is as follows:

where M represents the number of electrode sites, M-3,representing alternative leads x_iCorrelation with the mth target reference lead.

According to an embodiment of the present invention, in step S5, if the number of global leads is less than the set value, the set number is 50% of the number of global leads, and if the number of global leads is greater than or equal to the set value, the set number is 23% -27% of the number of global leads; the overall number of leads is the sum of the target set of reference leads and the alternative reference leads.

In step S1, the originally acquired EEG signal is filtered using a butterworth filter of order five, according to one embodiment of the present invention.

The invention also provides a lead selection device based on the correlation motor imagery, which comprises:

the preprocessing module is suitable for filtering the EEG signals which are originally acquired by using a Butterworth filter to obtain EEG signals of 0.5-30 Hz;

a data partitioning module adapted to partition all leads of EEG signals obtained after processing by the pre-processing module into a target reference lead set comprising target leads of electrode positions C3, C4 and Cz and a set of alternative reference leads comprising alternative leads of non-electrode positions C3, C4 and Cz;

the correlation calculation module is suitable for traversing each alternative lead and calculating the correlation between the alternative lead and each target lead respectively to obtain the correlation between each alternative lead and each target lead;

the weight calculation module is suitable for calculating corresponding weight according to the correlation and calculating a combined correlation index according to the weight;

a sorting module adapted to sort the alternative leads according to the combined relevance indicator, select a set number of the alternative leads sorted top, and combine them with the target lead to form a target lead set.

According to an embodiment of the present invention, the lead selection apparatus further comprises an evaluation module, the evaluation module is adapted to perform feature extraction on the target lead set by using a common space mode algorithm and perform classification by using a support vector machine.

The present invention also provides a computer readable medium having stored thereon computer instructions which, when executed, perform the steps of the aforementioned lead selection method.

The invention also provides a lead selection system comprising a memory and a processor, the memory having stored thereon computer instructions executable on the processor, the processor executing the computer instructions to perform the steps of the lead selection method as previously described.

The invention provides a correlation-based motor imagery lead selection method and device, a computer readable medium and a lead selection system, which are used for calculating correlation weights among leads and constructing a target lead set according to weight sequencing, so that effective leads are selected to improve the identification efficiency of a motor imagery task.

Drawings

In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below, wherein:

figure 1 shows an electrode placement position profile under an international standard 10/20EEG measurement;

FIG. 2 shows a flow chart of a correlation-based motor imagery lead selection method according to an embodiment of the invention.

Fig. 3 shows a schematic structural diagram of a lead selection device based on correlation motor imagery according to an embodiment of the present invention.

Detailed Description

In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below.

In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those specifically described herein, and thus the present invention is not limited to the specific embodiments disclosed below.

As used in this application and the appended claims, the terms "a," "an," "the," and/or "the" are not intended to be inclusive in the singular, but rather are intended to be inclusive in the plural unless the context clearly dictates otherwise. In general, the terms "comprises" and "comprising" merely indicate that steps and elements are included which are explicitly identified, that the steps and elements do not form an exclusive list, and that a method or apparatus may include other steps or elements.

In describing the embodiments of the present application in detail, the cross-sectional views illustrating the structure of the device are not enlarged partially in a general scale for convenience of illustration, and the schematic drawings are only examples, which should not limit the scope of the present application. In addition, the three-dimensional dimensions of length, width and depth should be included in the actual fabrication.

For convenience in description, spatial relational terms such as "below," "beneath," "below," "under," "over," "upper," and the like may be used herein to describe one element or feature's relationship to another element or feature as illustrated in the figures. It will be understood that these terms of spatial relationship are intended to encompass other orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "below" or "beneath" other elements or features would then be oriented "above" the other elements or features. Thus, the exemplary words "below" and "beneath" can encompass both an orientation of up and down. The device may have other orientations (rotated 90 degrees or at other orientations) and the spatial relationship descriptors used herein should be interpreted accordingly. Further, it will also be understood that when a layer is referred to as being "between" two layers, it can be the only layer between the two layers, or one or more intervening layers may also be present.

In the context of this application, a structure described as having a first feature "on" a second feature may include embodiments in which the first and second features are formed in direct contact, and may also include embodiments in which additional features are formed in between the first and second features, such that the first and second features may not be in direct contact.

The motor imagery EEG signal is characterized in that specific cortex of a brain is activated when a subject performs limb motor imagery, and at the moment, EEG signals with certain rules can be recorded on the scalp of the motor cortex by using electrodes.

The central nervous system of humans comprises the spinal cord and the brain. The brain is divided into three main areas: brainstem, cerebellum and brain. The brainstem is an extension of the spinal cord that connects the cerebral cortex, spinal cord and cerebellum, and it is also the control center for motor reflexes. Within the nervous system, there are ascending nerve bundles from the spinal cord responsible for transmitting information to the brainstem (e.g., pain or touch), and descending nerve bundles connecting the brain and cerebellum, etc., to motor neurons, thereby controlling skeletal muscle movement.

Motor imagery EEG signals recorded by the placement of electrodes on the human scalp when the subject imagines a particular limb movement, the recorded brain potentials are in fact a superposition of collective discharges of neuronal populations within the brain. Brain activity is associated with different events and stimulation types, e.g. sensory stimulation may cause changes in neuronal collective activity that can be observed in EEG signals.

Figure 1 shows an electrode placement position profile under an international standard 10/20EEG measurement. There are many studies that use electrode signals at electrode positions C3, C4, CZ and their vicinity as data sources when performing data analysis in terms of motor imagery. Since C3, C4 and CZ contain most of the motor imagery related information corresponding to the left hand, right hand and foot.

FIG. 2 shows a flow chart of a correlation-based motor imagery lead selection method according to an embodiment of the invention. As shown, a correlation-based motor imagery lead selection method includes:

step S1, preprocessing, namely filtering the EEG signals which are originally acquired by using a Butterworth filter to obtain EEG signals of 0.5-30 Hz;

step S2, data division, dividing all leads used for recording the EEG signals obtained in step S1 into a target reference lead set comprising target leads for electrode positions C3, C4 and Cz and an alternative reference lead set comprising alternative leads for non-electrode positions C3, C4 and Cz, i.e. the leads remaining except for electrode positions C3, C4 and Cz;

step S3, calculating the correlation, namely calculating the correlation between the alternative leads and the target lead, traversing each alternative lead, and calculating the correlation between the alternative lead and each target lead in the target reference lead set respectively to obtain the correlation between each alternative lead and each target lead;

step S4, calculating weight, calculating corresponding weight according to the correlation, and calculating combined correlation index according to the weight;

and step S5, sorting and selecting, namely sorting the alternative leads according to the combined relevance index, selecting the alternative leads with the set number in the front sorting, and combining the alternative leads with the target leads to form a target lead set. The combined relevance index is equivalent to a weighted value, the combined relevance index is sorted in size, the weighted value in the top sorting is large, and the relevance is higher.

Preferably, the lead selection method further comprises:

and step S6, extracting features, namely extracting the features of the target lead set by adopting a common space mode algorithm, and taking the result obtained after spatial feature filtering in the common space mode as a feature vector. The real nerve activity on the cerebral cortex is collected through the multi-channel electrode, and the result obtained after the filtering through the common space mode is proved to be an important signal source of the identification information.

And step S7, classifying, namely classifying by adopting a support vector machine. And classifying the extracted features through a support vector machine so as to accurately identify the movement tendency.

Preferably, in step S3, a target reference lead set Y ═ Y is defined₁，y₂，y₃}，y₁，y₂，y₃Corresponding to electrode positions C3, C4 and Cz, respectively, each y_iRepresents eachEEG signal recorded by one lead, dimension is 1 × T, T is number of sampling points, and alternative reference lead set X ═ X is defined₁，x₂，...，x_nDimension is 1 × T, T is the number of sample points, x for each alternative lead_iCalculating its and target lead y_iThe correlation between them;

defining a correlation asThe calculation formula for each component is:

y_pthe p-th target reference lead variable is represented, and I represents the mutual information calculation. In the process of evaluating lead similarity by using mutual information, the mutual information between each alternative lead and all target leads is calculated, so that the relation between the alternative leads and the target lead set can be evaluated in an all-around manner. For each pair of leads (x) for which mutual information needs to be calculated in equation 1_i，y_p) Wherein the sampling point recorded by each lead is T,the formula of mutual information calculation is:

in the formula 2, the first and second groups,is the current x_iAnd y_pThe joint probability density distribution function of (a),is x_iThe edge probability density function of (1).

Preferably, in step S4, correlation is performedThe weight is calculated to obtain a weighted combined mutual information delta_iThe calculation formula is as follows:

where M represents the number of electrode positions, electrode positions C3, C4 and Cz are included in the present invention, so that M is 3,representing alternative leads x_iCorrelation with the mth target reference lead. The proposed weighted combination mutual information can effectively integrate the measured mutual information to achieve the omnibearing evaluation of the alternative leads, and the evaluation of the alternative leads by the weighted combination mutual information can effectively select the leads containing important information.

Preferably, in step S5, if the number of global leads is less than the set value, the number is set to 50% of the number of global leads, and if the number of global leads is greater than or equal to the set value, the number is set to 23% to 27% of the number of global leads; the overall number of leads is the sum of the target reference lead set and the alternative reference leads. By way of example and not limitation, in the present embodiment, the set value is 32.

Preferably, in step S1, the originally acquired EEG signal is filtered using a butterworth filter of order five.

Fig. 3 shows a schematic structural diagram of a lead selection device based on correlation motor imagery according to an embodiment of the present invention. As shown, a correlation-based motor imagery lead selection apparatus 300 includes a pre-processing module 301, a data partitioning module 302, a correlation calculation module 303, a weight calculation module 304, and an ordering module 305.

The preprocessing module 301 is adapted to filter the originally acquired EEG signal by using a butterworth filter to obtain an EEG signal of 0.5 to 30 Hz.

The data partitioning module 302 is adapted to partition all leads of the EEG signals obtained after processing by the pre-processing module 301 into a target reference lead set comprising target leads for electrode positions C3, C4 and Cz, and a set of alternative reference leads comprising alternative leads for non-electrode positions C3, C4 and Cz.

The correlation calculation module 303 is adapted to traverse each candidate lead and calculate its correlation with each target lead separately, resulting in a correlation of each candidate lead with the target lead.

The weight calculation module 304 is adapted to calculate corresponding weights according to the correlations, and to calculate a combined correlation index according to the weights.

The ranking module 305 is adapted to rank the alternative leads according to the combined relevance indicator, selecting a set number of alternative leads ranked top, to be combined with the target lead into a target lead set.

Preferably, lead selection apparatus 300 further comprises an evaluation module 306. The evaluation module 306 is adapted to perform feature extraction on the target lead set using a common spatial mode algorithm and classify using a support vector machine.

The present invention also provides a computer readable medium having stored thereon computer instructions which, when executed, perform the steps of the aforementioned lead selection method.

The invention also provides a lead selection system, which comprises a memory and a processor, wherein the memory is stored with computer instructions capable of running on the processor, and the processor executes the computer instructions to execute the steps of the lead selection method.

Some aspects of the correlation-based motor imagery lead selection method of the present invention may be performed entirely by hardware, entirely by software (including firmware, resident software, micro-code, etc.), or by a combination of hardware and software. The above hardware or software may be referred to as "data block," module, "" engine, "" unit, "" component, "or" system. The processor may be one or more Application Specific Integrated Circuits (ASICs), Digital Signal Processors (DSPs), digital signal processing devices (DAPDs), Programmable Logic Devices (PLDs), Field Programmable Gate Arrays (FPGAs), processors, controllers, microcontrollers, microprocessors, or a combination thereof. Furthermore, aspects of the present application may be represented as a computer product, including computer readable program code, embodied in one or more computer readable media. For example, computer-readable media may include, but are not limited to, magnetic storage devices (e.g., hard disk, floppy disk, magnetic strips … …), optical disks (e.g., Compact Disk (CD), Digital Versatile Disk (DVD) … …), smart cards, and flash memory devices (e.g., card, stick, key drive … …).

The computer readable medium may comprise a propagated data signal with the computer program code embodied therein, for example, on a baseband or as part of a carrier wave. The propagated signal may take any of a variety of forms, including electromagnetic, optical, and the like, or any suitable combination. The computer readable medium can be any computer readable medium that can communicate, propagate, or transport the program for use by or in connection with an instruction execution system, apparatus, or device. Program code on a computer readable medium may be propagated over any suitable medium, including radio, electrical cable, fiber optic cable, radio frequency signals, or the like, or any combination of the preceding.

According to the correlation-based motor imagery lead selection method and device, the computer readable medium and the lead selection system, a few lead sets with important information can be selected from the high-density lead sets, the screened lead sets can cover the important information for identifying the motor imagery task, and the influence of artifact information and noise on a brain-computer interface system is effectively avoided. The invention aims to measure the correlation among leads of an EEG data acquisition system of a motor imagery brain-computer interface, select leads of electrode positions C3, C4 and Cz as a target reference lead set, use the rest as alternative reference lead sets, construct a new joint lead evaluation index according to the correlation measurement index between each alternative lead and the target lead, sort the leads according to the new lead evaluation index, and select a set number of leads to form the target lead set. The lead selection process mainly comprises the following steps: data preprocessing, partitioning of target reference lead set and alternative reference lead set, calculation of correlation and mutual information thereof, calculation of alternative lead weighted combination mutual information and sorting selection process of alternative leads.

Although the present invention has been described with reference to the present specific embodiments, it will be appreciated by those skilled in the art that the above embodiments are merely illustrative of the present invention, and various equivalent changes and substitutions may be made without departing from the spirit of the invention, and therefore, it is intended that all changes and modifications to the above embodiments within the spirit and scope of the present invention be covered by the appended claims.