This application is based on and claims priority under 35 USC 119 from Japanese Patent Application No. 2022-148335 filed on Sep. 16, 2022, the entire content of which is incorporated herein by reference.
The present disclosure relates to a physiological information processing method and a physiological information processing apparatus. In particular, the present disclosure relates to a physiological information processing method and a physiological information processing apparatus in which information related to an artifact mixed in an electroencephalogram signal is visually presented. Further, the present disclosure further relates to a program that causes a computer to execute the information processing method, and a non-transitory computer-readable storage medium storing the program.
An electroencephalograph configured to noninvasively measures an electroencephalogram indicating an electrical activity of brain neurons of a patient is known. In the electroencephalograph, electroencephalogram signals are obtained through a plurality of electrodes attached to a head of the patient. On the other hand, since the electroencephalogram signals are weak electrical signals of a living body which indicate the electrical activity of the brain neurons, artifacts, which are noises other than the electroencephalogram, are likely to be mixed in the electroencephalogram signals. U.S. Pat. No. 9,055,927B discloses a technique for detecting, separating, and reducing various artifacts mixed in electroencephalogram signals of respective channels. Further, U.S. Pat. No. 9,055,927B discloses a technique of indicating a type of each of the artifacts and re-extracting an electroencephalogram component signal according to the type of artifact to be reconstructed.
At this time, in the electroencephalogram signals from which various artifacts have been reduced by software through a signal process of reducing the artifacts, waveform distortion may occur compared to original electroencephalogram signals before an artifact is reduced. As described above, reliability of the electroencephalogram signals from which the artifacts are reduced through the signal process is lower than reliability of electroencephalogram signals in which the artifacts are not mixed. Therefore, in order to record electroencephalogram signals having a high signal quality, it is essential to eliminate factors of the artifacts when electroencephalogram measurement is started.
On the other hand, the artifacts mixed in the electroencephalogram include various types of artifacts such as artifacts caused by a living body (such as electro cardiogram and eye movement) and artifacts caused by factors other than the living body (such as a contact failure of an electrode). Further, measures for reducing various artifact differ from each other. For this reason, a medical worker who has little experience and is unfamiliar with the artifact reduction measure may find it difficult to record an electroencephalogram signal having a signal quality that is less contaminated with the artifact. From the above viewpoint, there is room for consideration of a user interface to support recording of the electroencephalogram signal having a signal quality that is less contaminated with the artifact.
Aspect of non-limiting embodiments of the present disclosure relates to provide a physiological information processing method and a physiological information processing apparatus that support recording of an electroencephalogram signal having a high signal quality, and relates to provide a program that causes a computer to execute the information processing method, and a computer readable storage medium storing the program.
Aspects of certain non-limiting embodiments of the present disclosure address the features discussed above and/or other features not described above. However, aspects of the non-limiting embodiments are not required to address the above features, and aspects of the non-limiting embodiments of the present disclosure may not address features described above.
According to an aspect of the present disclosure, there is provided a physiological information processing method executed by a computer, the physiological information processing method including:
According to an aspect of the present disclosure, there is provided a physiological information processing apparatus including:
Exemplary embodiment(s) of the present invention will be described in detail based on the following figures, wherein:
Hereinafter, the present embodiment will be described with reference to the drawings.
The processing apparatus 1 may be a medical device (for example, a patient monitor or an electroencephalograph) configured to display physiological information of a subject P (patient), a personal computer, a workstation, a smartphone, a tablet, or a wearable device (for example, an AR glass) worn on a body (for example, an arm or a head) of a medical worker. The processing apparatus 1 is configured to output information indicating electroencephalogram of the subject P.
The controller 2 may include one or more memories and one or more processors. Each memory is configured to store a computer readable instruction (program). For example, the memory may include a read only memory (ROM) in which various programs are stored, a random access memory (RAM) having a plurality of work areas in which various programs executed by the one or more processors are stored, and the like. Each processor may include at least one of, for example, a central processing unit (CPU), a micro processing unit (MPU), and a graphics processing unit (GPU). The CPU may include a plurality of CPU cores. The GPU may include a plurality of GPU cores. The processor may be configured to load a program designated from various programs incorporated in the storage device 3 or ROM into RAM, and configured to execute various processes in cooperation with the RAM. In particular, the processor is configured to load a physiological information processing program for executing a series of processes illustrated in
The storage device 3 is a storage device (storage) such as a hard disk drive (HDD), a solid state drive (SSD), or a flash memory, and stores programs and various data. The physiological information processing program may be incorporated in the storage device 3. Further, the storage device 3 may store electroencephalogram data of the subject P. For example, the electroencephalogram data (electroencephalogram signals) obtained from an electroencephalogram sensor 10 may be stored in the storage device 3 via the sensor interface 8.
The communication unit 5 is configured to connect the processing apparatus 1 to an in-hospital network. Specifically, the communication unit 5 may include various wired connection terminals configured to communicate with a central monitor or a server disposed in the in-hospital network. Further, the communication unit 5 may include a wireless communication module configured to wirelessly communicate with the central monitor or the server. The communication unit 5 may include, for example, a wireless communication module corresponding to a medical telemetry system. The communication unit 5 may include a wireless communication module corresponding to a wireless communication standard such as Wi-Fi (registered trademark) or Bluetooth (registered trademark) and/or a wireless communication module corresponding to a mobile communication system using SIM. The in-hospital network may include, for example, a local area network (LAN) or a wide area network (WAN). The processing apparatus 1 may be connected to Internet via the in-hospital network.
The display 4 is configured to display the physiological information (information related to the electroencephalogram) of the subject P. The display 4 can include, for example, a liquid crystal panel or an organic EL panel. The input operation unit 6 is, for example, a touch panel, a mouse, and/or a keyboard, which are arranged over the display 4. The input operation unit 6 is configured to receive an input operation of the medical worker and configured to generate an operation signal corresponding to the input operation of the medical worker. After the operation signal generated by the input operation unit 6 is transmitted to the controller 2 via the bus 14, the controller 2 is configured to execute a predetermined operation according to the operation signal. The voice output unit 7 can include one or more speakers.
The sensor interface 8 is an interface configured to connect the electroencephalogram sensor 10 to the processing apparatus 1. The sensor interface 8 may include an input terminal to which the electroencephalogram signals output from the electroencephalogram sensor 10 are input. The electroencephalogram sensor 10 is configured to noninvasively measure the electroencephalogram signals indicating an electrical activity of brain neurons of the subject P. The electroencephalogram sensor 10 can include a plurality of electrodes attached to a head of the subject P.
As illustrated in
As a method of synthesizing the electroencephalogram, a homopolar lead method for recording a potential difference between potential of a head and potential of an earlobe, a bipolar lead method for recording a potential difference of the head between two points, or an average reference electrode method is applied.
The sensor interface 8 may include at least a plurality of amplifier circuits and an AD converter. Each of the plurality of amplifier circuits is configured to amplify the electroencephalogram signal output from the corresponding electrode. The AD converter is configured to convert the electroencephalogram signal from an analog signal to a digital signal. The electroencephalogram signal converted into the digital signal is transmitted from the sensor interface 8 to the controller 2.
Next, the physiological information processing method according to the present embodiment will be described below with reference to
As illustrated in
Next, the controller 2 is configured to separate, into a plurality of artifact component signals and an electroencephalogram component signal, each of the plurality of electroencephalogram signals (step S2 in
The independent component analysis is a calculation method for separating a multivariate signal into a plurality of additive components. In the present embodiment, since the 21 types of electroencephalogram signals are simultaneously obtained by 21 electrodes, the 21 types of electroencephalogram signals can be decomposed into components of 21 types of source waveforms S1 to S21 by the independent component analysis (see
In a case where the 21 types of electroencephalogram signals X(t) are set to X(t)=(XFp1(t), . . . , XA2(t))T, and the 21 types of source waveforms S(t) are (S1(t), . . . , S21(t))T, a relationship between X(t) and S(t) in the independent component analysis is expressed as X(t)=H·S(t). Here, H is a mixing matrix of 21 rows×21 columns. Specifically, the relationship between X(t) and S(t) is expressed by the following formula (1).
Thus, the controller 2 can express the electroencephalogram signal X(t) as H·S(t), through the independent component analysis. That is, each of the electroencephalogram signals XFp1 to XA2 can be expressed by 21 types of source signals S1 to S21. Here, the electroencephalogram signal XFp1 of the electrode Fp1 can be expressed as XFp1=h1 1S1+h1 2S2+ . . . h1 21S21. The electroencephalogram signal XA2 of the electrode A2 can be expressed as XA2=h21 1S1+h21 2S2+ . . . h21 21S21. In the following description, the electroencephalogram signals may be collectively referred to simply as electroencephalogram signals X.
Next, the controller 2 is configured to classify the 21 types of source waveforms S1 to S21. In the present embodiment, the source waveform includes a source waveform indicating the electroencephalogram and a source waveform indicating the artifact. The source waveform indicating the artifact includes a source waveform indicating an artifact caused by a living body of the subject P, and a source waveform indicating an artifact caused by a factor other than the living body. The source waveform indicating the artifact caused by the living body includes a source waveform indicating an artifact caused by an electro cardiogram of the subject P, a source waveform indicating an artifact caused by eye movement of the subject P, and a source waveform indicating an artifact caused by an electromyogram of the subject P. On the other hand, the source waveform indicating the artifact caused by the factor other than the living body includes a source waveform indicating an artifact caused by an attachment failure of the electrode, and a source waveform indicating an artifact caused by an electrode lead wire.
As described above, the source waveforms S1 to S21 are classified into any one of i) the source waveform indicating the electroencephalogram, ii) the source waveform indicating the artifact caused by the electro cardiogram, iii) the source waveform indicating the artifact caused by the eye movement, iv) the source waveform indicating the artifact caused by the electromyogram, v) the source waveform indicating the artifact caused by the attachment failure of the electrode, and vi) the source waveform indicating the artifact caused by the electrode lead wire. As an example of a method for classifying the source waveform, a type of the source waveform may be determined based on a correlation coefficient between the source waveform and each of reference waveforms related to various artifacts. For example, in a case where the correlation coefficient between a source waveform S10 and a reference waveform indicating the artifact caused by the electrode lead wire is high, the source waveform S10 may be classified as the source waveform indicating the artifact caused by the electrode lead wire. In this manner, types of the artifacts indicated by the source waveforms are identified based on correlation coefficients between the source waveforms and the reference waveforms for the various artifacts. Further, in a case where a predetermined source waveform does not correlate with any type of artifact-related reference waveform, the predetermined source waveform may be classified as the source waveform indicating the electroencephalogram.
For details of the method for classifying each of the source waveforms, reference is made to the aforementioned U.S. Pat. No. 9,055,927B and the following non-patent literature. The following non-patent literature describes an example of a method for classifying the source waveform in detail. Non-patent Literature: Wallstrom G L, Kass R E, Miller A, Cohn J F, Fox N A (2004). Automatic correction of ocular artifacts in the EEG: a comparison of regression-based and component-based methods. Int J Psychophysiol 53: 105-119.
As described above, after types of the source waveforms S1 to S21 are specified, the controller 2 is configured to separate each of the electroencephalogram signals X into the electroencephalogram component signal B indicating the component of the electroencephalogram and the artifact component signals C1 to C5. As illustrated in
As illustrated in
For example, as illustrated in
Referring back to
In step S4, the controller 2 is configured to determine a signal quality index (SQI) of each of the electroencephalogram signals XFp1 to XA2. The SQI is an index indicating a signal quality of the electroencephalogram signal. In this respect, the controller 2 is configured to determine SQI of each of the electroencephalogram signals X, based on RMS of the electroencephalogram signal X and some of RMS of the artifact component signals C1 to C5 of the electroencephalogram signal X. More specifically, SQI of each of the electroencephalogram signals X is calculated based on the following formula (2). SQI is shown in percentage (%). The higher a value of SQI of the electroencephalogram signal X, the higher the signal quality of the electroencephalogram signal X. Thus, the signal quality of each of electroencephalogram signals can be objectively evaluated through the value of SQI.
Here, RMSx indicates RMS of the electroencephalogram signal X. RMSC5 indicates RMS of the artifact component signal C5 caused by the electromyogram. RMS−C2 indicates RMS of the artifact component signal C2 caused by the electrode lead wire. RMSc1 indicates RMS of the artifact component signal C1 caused by the attachment failure of the electrode. Thus, SQI of each of the electroencephalogram signals X is calculated based on RMSx of each electroencephalogram signal X and RMS (RMSC1, RMSC2, RMSC5) of the artifact component signals C1, C2, and C5 of each of the electroencephalogram signals X. For example, SQI of the electroencephalogram signal XFp1 of the electrode Fp1 is calculated based on RMS of the electroencephalogram signal XFp1 and RMS of the artifact component signals C1, C2, and C5 of the electroencephalogram signal XFp1. In the present embodiment, RMS of the artifact component signals C3 and C4 are not used in the calculation of the SQI, but these RMS may be used in the SQI calculation.
Next, in step S5, the controller 2 is configured to determine a comprehensive SQI of the plurality of electroencephalogram signals X, based on SQI of electroencephalogram signals X. Here, the comprehensive SQI indicates a comprehensive signal quality of the plurality of electroencephalogram signals X. In this respect, in the present embodiment, the controller 2 is configured to determine a minimum value among SQI of the electroencephalogram signals X as the comprehensive SQL. A representative value such as an average value or a median value of each electroencephalogram signal X may be determined as the comprehensive SQI.
After the controller 2 determines the comprehensive SQI, the controller 2 may be configured to change a display content of an SQI indicator 32 indicating the value of the SQI on an electroencephalogram display screen 30 illustrated in
Next, in step S6, the controller 2 is configured to determine whether the comprehensive SQI is 50% or less. In a case where a determination result of step S6 is NO, the processing returns to step S1. On the other hand, in a case where the comprehensive SQI is 50% or less (the determination result of step S6 is YES), the controller 2 is configured to determine whether RMS of each of the artifact component signals C1 to C5 of each of the electroencephalogram signals X (electroencephalogram signals XFp1 to XA2) is greater than a threshold value (step S7). In this example, the controller 2 is configured to determine whether the comprehensive SQI is 50% or less, but in step S6, the controller 2 may be configured to determine whether the comprehensive SQI is X % (X is any value other than 50) or less.
In step S7, the controller 2 is configured to determine whether RMSC1 of the artifact component signal C1 of the electroencephalogram signal X is greater than a threshold value RMSth1 associated with RMSC1. Further, the controller 2 is configured to determine whether the RMSC2 of the artifact component signal C2 of the electroencephalogram signal X is greater than a threshold value RMSth2 associated with RMSC2. The controller 2 is configured to determine whether RMSC3 of the artifact component signal C3 of the electroencephalogram signal X is greater than a threshold value RMSth3 associated with RMSC3. Furthermore, the controller 2 is configured to determine whether RMSC4 of the artifact component signal C4 of the electroencephalogram signal X is greater than a threshold value RMSth4 associated with RMSC4. The controller 2 is configured to determine whether RMSC5 of the artifact component signal C5 of the electroencephalogram signal X is greater than a threshold value RMSth5 associated with RMSC5.
The processing of determining whether RMS of each of the artifact component signals C1 to C5 is greater than the respective one of the threshold values RMSth1 to RMSth5 is performed on each of the electroencephalogram signals XFp1 to XA2. Therefore, the threshold value determination processing related to RMS of the artifact component signals C1 to C5 is performed for each of the 21 types of electroencephalogram signals XFp1 to XA2. That is, 105 times, which is that 21 multiplied by 5, of threshold value determination processing are executed in step S6.
For example, in a case where RMSC1 of the artifact component signal C1 of the electroencephalogram signal XFp1 is greater than the threshold value RMSth1, the controller 2 is configured to determine that the artifact caused by the attachment failure of the electrode is mixed in the electroencephalogram signal XFp1. In a case where RMSC2 of the artifact component signal C2 of the electroencephalogram signal XFp1 is greater than the threshold value RMSth2, the controller 2 is configured to determine that the artifact (AC noise) caused by the electrode lead wire is mixed in the electroencephalogram signal XFp1. In a case where RMSC3 of the artifact component signal C3 of the electroencephalogram signal XFp1 is greater than the threshold value RMSth3, the controller 2 is configured to determine that the artifact (electro cardiogram noise) caused by the electro cardiogram is mixed in the electroencephalogram signal XFp1. In a case where RMSC4 of the artifact component signal C4 of the electroencephalogram signal XFp1 is greater than the threshold value RMSth4, the controller 2 is configured to determine that the artifact (eye movement noise) caused by the eye movement is mixed in the electroencephalogram signal XFp1. In a case where RMSC5 of the artifact component signal C5 of the electroencephalogram signal XFp1 is greater than the threshold value RMSth5, the controller 2 is configured to determine that the artifact (electromyogram noise) caused by the electromyogram is mixed in the electroencephalogram signal XFp1.
In step S8, the controller 2 is configured to visually display information related to the artifact associated with the artifact component signal of which RMS is greater than the threshold value. In this respect, the controller 2 may be configured to cause the electroencephalogram display area 31 of the electroencephalogram display screen 30 to display artifact related information 33 indicating information related to the artifact.
The artifact related information 33 includes information indicating the type of artifact, information indicating an electrode associated with an electroencephalogram signal in which the artifact is mixed, and information indicating a measure for reducing the artifact. For example, in a determination result of step S7, RMSC5 of the artifact component signal C5 of the electroencephalogram signal X associated with each of the electrodes F3, O1, P3, Cz, F4, and Pz is greater than RMSth5. In such a case, as illustrated in
As illustrated in
In a case where RMSC2 of the artifact component signal C2 of the electroencephalogram signal X related to a predetermined electrode is greater than RMSth2, information indicating the AC noise, information indicating an electrode associated with an electroencephalogram signal in which the AC noise is mixed, and information indicating a measure for removing the AC noise are displayed in the artifact related information 33-3. In a case where RMSC1 of the artifact component signal C1 of the electroencephalogram signal X related to the predetermined electrode is greater than RMSth1, information indicating a noise (electrode attachment failure noise) caused by the electrode attachment failure, information indicating an electrode associated with an electroencephalogram signal in which the electrode attachment failure noise is mixed, and information indicating a measure for reducing the electrode attachment failure noise are displayed in the artifact related information 33-4.
In the example shown in
In particular, among artifacts caused by factors other than the living body, the artifact (electrode attachment failure noise) caused by the attachment failure of the electrode may have a higher priority than the artifact (AC noise) caused by the electrode lead wire. Further, among the artifacts caused by the living body, the artifact (electromyogram noise) caused by the electromyogram may have a highest priority, and the artifact (electro cardiogram noise) caused by the electro cardiogram may have a higher priority than the artifact (eye movement noise) caused by the eye movement. That is, a priority order of information display related to each artifact is as follows. Electrode attachment failure noise>AC noise>electromyogram noise>electro cardiogram noise>eye movement noise. For example, in a case in which five noises are simultaneously mixed in the electroencephalogram signals XFp1 to XA2, information related to the electrode attachment failure noise is first displayed, and after the noise is reduced by the medical worker through the measure against the artifact, information related to the AC noise is displayed.
The priority order of the information display related to each artifact is determined in consideration of difficulty of reducing the artifacts and of artifacts that greatly affect the electroencephalogram waveform. For example, since difficulty of reducing the artifact caused by the factor other than the living body is lower than difficulty of reducing the artifact caused by the living body, information related to the artifact caused by the factor other than the living body is preferentially displayed. Furthermore, since an influence of the electrode attachment failure noise on the electroencephalogram waveform is larger than an influence of the AC noise on the electroencephalogram waveform, display of the information related to the electrode attachment failure noise has a higher priority than display of the information related to the AC noise. The priority order of the information display related to each artifact is not particularly limited. The priority order may be appropriately set by the medical worker.
As illustrated in
Specifically, in the example illustrated in
Furthermore, in a case where the artifacts (electromyogram noise, electro cardiogram noise, and eye movement noise) caused by the living body are mixed in the electroencephalogram signals X related to the electrodes Fp1, Fp2, F7, F3, Fz, F4, F8, A1, and A2, information indicating that the artifacts caused by the living body are mixed in the electroencephalogram signals X of the electrodes Fp1, Fp2, F7, F3, Fz, F4, F8, A1, and A2 is displayed on the electrode image 41. More specifically, illustrations of the electrodes Fp1, Fp2, F7, F3, Fz, F4, F8, A1, and A2 are colored in a second color in a state in which the artifacts caused by the living body are associated with the second color different from the first color.
In a case where no artifact is mixed in the electroencephalogram signals X related to the electrodes T3, C3, Cz, C4, T4, T5, P3, Pz, P4, and T6, illustrations of the electrodes T3, C3, Cz, C4, T4, T5, P3, Pz, P4, and T6 are not colored. In this manner, the medical worker can intuitively grasp the type of artifact mixed in the electroencephalogram signal of each electrode by visually recognizing the illustration of each electrode displayed on the electrode image 41.
In the example illustrated in
The electrode display screen 40 illustrated in
In this manner, the controller 2 is configured to repeatedly execute a series of processing of steps S1 to S8 shown in
According to the present embodiment, the controller 2 is configured to determine whether RMS of each of the artifact component signals C1 to C5 of each of the electroencephalogram signals X satisfies the threshold value, and the type of artifact mixed in the electroencephalogram signal X is specified. Thereafter, the artifact related information (particularly, information indicating the type of artifact and information indicating the measure for reducing the artifact) related to the artifact mixed in each of the electroencephalogram signals X is visually presented to the medical worker. In this way, it is possible to take the measure for reducing the artifact from the electroencephalogram signal through the artifact related information even for a medical worker who has little experience in electroencephalogram measurement. As a result, the medical worker can record an electroencephalogram signal having a high signal quality which is less contaminated with the artifact.
In the electrode image 41 illustrated in
Further, according to the present embodiment, in a case where the comprehensive SQI is 50% (an example of a predetermined threshold value) or less, it is determined whether RMS (an example of the effective value) of each of the artifact component signals C1 to C5 of each electroencephalogram signal X satisfies the predetermined threshold value. As described above, only in the case in which the comprehensive SQI is low and a probability that the electroencephalogram signal X is mixed with the artifact is extremely high, a determination processing related to RMS of the artifact component signals C1 to C5 is executed. Accordingly, it is possible to reduce a calculation load and power consumption of a computer configured to execute the determination processing.
In the present embodiment, it is determined whether RMS of each artifact component signal of each electroencephalogram signal is greater than the threshold value in a case where the comprehensive SQI of the electroencephalogram signal X is 50% or less in step S6, but the present embodiment is not limited thereto. For example, the processing of steps S4 to S6 may not be executed in the series of processing of
In order to achieve the processing apparatus 1 according to the present embodiment by software, the physiological information processing program may be incorporated in the storage device 3 or the ROM in advance. Alternatively, the physiological information processing program may be stored in a computer readable storage medium such as a magnetic disk (for example, HDD and a floppy disk), an optical disk (for example, CD-ROM, DVD-ROM, and Blu-ray (registered trademark) disk), a magneto optical disk (for example, MO), a flash memory (for example, a SD card, a USB memory, and SSD). In this case, the physiological information processing program stored in the storage medium may be incorporated in the storage device 3. Further, after the program incorporated in the storage device 3 is loaded onto RAM, the processor may be configured to execute the program loaded on RAM. As described above, the physiological information processing method according to the present embodiment is executed by the processing apparatus 1.
The physiological information processing program may be downloaded from a computer on a communication network via the communication unit 5. In this case, the downloaded program may be incorporated in the storage device 3.
Although the embodiments of the presently disclosed subject matter have been described above, the technical scope of the presently disclosed subject matter should not be construed as being limited to the description of the present embodiments. The present embodiments are merely an example, and it is understood by those skilled in the art that various modifications of the embodiments are possible within the scope of the disclosed subject matters described in the claims. The technical scope of the presently disclosed subject matter should be determined based on the scope of the disclosed subject matters described in the claims and equivalents thereof
Number | Date | Country | Kind |
---|---|---|---|
2022-148335 | Sep 2022 | JP | national |