METHOD OF DETERMINING BRAIN ACTIVITY

Abstract

A first aspect provides, in a device for processing biomedical data obtained from a body, a method of determining brain activity. The method comprises obtaining neuron oscillation signal data of a first part of a brain over a period of time. The oscillation signal data is apportioned in timeslots within the period and per time slot, based on the oscillation data for the timeslot, an amplitude value and variation value representing data values in the timeslot are obtained. A correlation value is determined between the amplitude values and the corresponding variation values over at least a substantial part of the first time period; and the correlation value is provided as an output. The correlation value may optionally be used for estimating a ratio between excitation and inhibition activity of the brain under scrutiny.

Claims

1. In a device for processing biomedical data obtained from a body, a A method of determining brain activity, comprising: obtaining neuron oscillation signal data of at least a part of a brain over a period of time; apportioning the neuron oscillation signal data in timeslots within the period; per time slot, based on the neuron oscillation signal data for the timeslot, obtaining an amplitude value and variation value representing data values in the timeslot; determining a correlation value between the amplitude values and the corresponding variation values over at least a substantial part of the time period; and; providing the correlation value as an output.

2. The method according to claim 1, wherein obtaining the variation value comprises determining a variation value of an amplitude envelope or power envelope of the signal data.

3. The method according to claim 1, wherein obtaining the variation value further comprises calculating, on a per data value basis, a cumulative sum based on amplitude envelope data values and determining a variation value of data after calculation of the cumulative sum.

4. The method according to claim 3, further comprising: obtaining a mean or median value of signal data; subtracting the mean or median value from momentary signal values; and summing of the resulting values.

5. The method according to claim 3, wherein obtaining the variation value further comprises detrending the signal data on a per timeslot basis after application of the summation and determining a variation of the detrended signal data.

6. The method according to claim 1, wherein obtaining the variation value comprises applying detrended fluctuation analysis on the signal data per timeslot.

7. The method according to claim 1, wherein obtaining an amplitude value comprises obtaining a mean amplitude per timeslot as the amplitude value.

8. The method according to claim 1, wherein obtaining an amplitude value comprises obtaining a median amplitude per timeslot as the amplitude value.

9. The method according to claim 1, wherein obtaining an amplitude value comprises: obtaining an amplitude envelope for the signal data; and obtaining a mean value or a median value of the amplitude envelope as the amplitude value.

10. The method according to claim 1, wherein the timeslots overlap.

11. The method according to claim 1, wherein the neuron oscillation signal data is obtained by at least one of: electro encephalography; magneto encephalography; single or multi-channel local field potentials; voltage sensitive dies; functional magnetic resonance imaging; and near-infrared spectroscopy.

12. The method according to claim 1, further comprising: if the correlation value is above a first threshold, determining the ratio of excitatory and inhibitory signalling of neurons to be inhibitory dominated; if the correlation value is below a second threshold, determining the ratio of excitatory and inhibitory signalling of neurons to be excitatory dominated; and if the correlation value is between the first and the second threshold, determining the ratio of excitatory and inhibitory signalling of neurons to be balanced.

13. The method according to claim 1, wherein the correlation value is used for providing an estimation of a ratio between excitatory and inhibitory signalling of neurons of at least a part of the brain from which the neuron oscillation signal data originates.

14. A non-transitory computer readable medium comprising computer executable code for programming a computer to enable the computer causing a processor to execute the method according to claim 1.

15. DA device for processing biomedical data obtained from a body, the device comprising: an input module configured to obtain neuron oscillation signal data of a at least a part of a brain over a period of time; a time division module configured to apportion oscillation data in timeslots within the period; a value determination module configured to, per time slot, based on the oscillation data for the timeslot, obtain an amplitude value and variation value representing data values in the timeslot; a correlation module configured to determine a correlation value between the amplitude values and the corresponding variation values over the first time period; and; an output configured to provide the correlation value.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The various aspects and examples thereof will now be discussed in further detail. In the drawings,

[0019] FIG. 1: shows a signal processing device 100 connected to an electro encephalogram signal acquiring device;

[0020] FIG. 2: shows a flowchart; and

[0021] FIG. 3 A through FIG. 3 G show graphs of processed signal data.

DETAILED DESCRIPTION

[0022] FIG. 1 shows an electro encephalogram signal acquiring device 150, a signal processing device 100 and a display screen 160. The electro encephalogram signal acquiring device 150 comprises electrodes 152 that may be placed on the head of a subject under scrutiny and a human being in particular. The electro encephalogram signal acquiring device 150 further comprises a signal conditioning circuit 154. The signal conditioning circuit 154 is arranged for signal conditioning of signals received from the electrodes.

[0023] The conditioning may comprise amplification, filtering, transformation, other, or a combination thereof. The signal conditioning circuit 154 may also be arranged for providing a signal to the electrodes 152, which electrodes 152 may modify the provided signal. The modified signal is received by the signal conditioning circuit 154 and may be conditioned. the electrodes 152 may be provided internally or externally of the person under scrutiny.

[0024] The signal processing device comprises an input module 132, an optional filter module 134, an analogue to digital converter 136, a processing module 102, a storage module 104 and an output module 138. The input module 132 is arranged for receiving a signal from the electro encephalogram signal acquiring device 150. As such, it may comprise a plug or a socket. The filter module 134 may be arranged for filtering the received signal. The filter module 134 is arranged to pass through AC signal components of a wide frequency range (e.g., 0-20.000 Hz), preferably with no high-pass filtering. Optionally, the filter module 134 may be a bandpass filter in a more narrow frequency range, for example 0.1-500 Hz or any other band. Optionally, a DC signal may be passed through as well.

[0025] Such filtering may comprise low pass, band pass or high pass filtering, other filtering, or a combination thereof. The filtering may be active or passive filtering. The analogue to digital converter 136 is arranged for sampling the filtered signal and transforming it to the digital domainwith respect to time and amplitude. Alternatively, the analogue to digital converter 136 is provided upstream of the filter module 134. In that case, the filter module 134 is a digital filter; otherwise an analogue filter.

[0026] The processing module 102 comprises various submodules for processing the digitised signal. It is noted that, optionally, the signal processing may also be arranged for processing a signal in the analogue domain, but the digital domain is preferred. For processing the signal, the processing module 102 comprises a timeslot submodule 112 for assigning one or more consecutive samples to a timeslot or, worded differently, a time window.

[0027] The processing module 102 further comprises an amplitude analysis submodule 114 arranged for determining an amplitude envelope of a signal. Alternatively or additionally, a power envelope may be obtained. It is noted that power of a signal is directly related to the amplitude of the signal. Therefore, in this document, unless specified otherwise, a reference to amplitude of the acquired signal is equivalent to a reference to power of the acquired signal, either processed or not.

[0028] Furthermore, the processing module 102 comprises a signal integration submodule 116 as a summing submodule. In a preferred embodiment, the signal integration module is arranged to sum a demeaned signal. Such demeaned signal is obtained by subtracting from an actual signal value a mean or median value calculated over a particular timeframe. The processing module 102 also comprises a normalising submodule 118 for normalising a signal. Additionally, the processing module 102 comprises a detrending submodule 120 for detrending a signal.

[0029] The processing module 102 further comprises a variation determining submodule 122 for determining variation of a signal within a particular amount of time, like within a timeslot. The variation of a signal may be determined as the standard deviation, but determining the variation by other means may be envisaged as well. The processing module 102 also comprises a correlation submodule 124 for determining correlation between two datasets.

[0030] Furthermore, the processing module 102 comprises a comparison module for comparing a calculated or otherwise obtained value with a pre-determined value that may be stored in the storage module 104.

[0031] It is noted that in various examples, the processing module 102 may comprise one or more of the processing submodules discussed above, depending on various demands in particular scenarios; not all submodules are required for implementation of the various aspects.

[0032] The storage module 104 is arranged for storing results of processing signals by the processing module 102. Furthermore, the storage module 104 may, as for example a non-transitory medium, comprise computer executable code for programming the processing module 102 to carry out various methods for signal processing. And the storage module 104 may hold various values used for signal processing.

[0033] The screen 160 is connected to an output module 138 comprised by the signal processing device 100. As such, the screen 160 may be provided with processed data signals, either in values, in graphical representations, other or a combination thereof.

[0034] FIG. 2 shows a flowchart 200 depicting a procedure for signal processing that may be carried out by the signal processing device 100 as shown by FIG. 1. The various parts of the flowchart are briefly summarised in the list below. It is noted that whereas parts of the flowchart 200 may be referred to as steps, it is noted that processing parts indicated as one step may in other examples be distributed over multiple steps and the other way around. The various parts are summarised as follows: [0035] 202 start [0036] 204 obtain signal [0037] 206 sample signal [0038] 207 bandpass filter signal [0039] 208 assign samples to timeslots [0040] 210 determine amplitude envelope [0041] 212 sum values [0042] 214 normalise values [0043] 216 detrend values [0044] 218 determine variation of the values [0045] 220 determine mean amplitude [0046] 222 correlate variation and mean amplitude for slots [0047] 224 higher than first threshold? [0048] 226 lower than second threshold? [0049] 228 balance excitation and inhibition [0050] 230 end [0051] 240 inhibition dominant [0052] 242 excitation dominant

[0053] The process starts in a terminator 202 and proceeds to step 204 in which the signal processing device 100, or, in another example, the signal processing module 102 receives a signal comprising data representing brain activity. Such signal may be obtained from the electro encephalogram signal acquiring device 150, but also data from magneto encephalography or invasively recorded local field potentials may be used. Additionally or alternatively, data is acquired by means of single or multi-channel local field potentials, voltage sensitive dies, functional magnetic resonance imaging or near-infrared spectroscopy. The received signal is processed by means of one or more different operations. Each operation may be executed by an applicable submodule of the signal processing module 102. Not every operation is necessary and various operations may be omitted or combined. Furthermore, other operations may be added. The operations may be executed as shown by the flowchart 200 or in another order without departing from a general inventive concept.

[0054] The process continues to step 206, in which the received signal is sampled and preferably also transformed to the digital domain with respect to the amplitude. The signal thus obtained is optionally provided to a further bandpass filter in step 207. The filtering of the further bandpass filter is preferably obtained by filtering in bands that are narrower than an optional filtering band of the bandpass filter 134. The pass band of the further bandpass filter may be adjustable providing for example a band from 1 Hz to 4 Hz, from 4 Hz to 8 Hz, from 8 Hz to 13 Hz and so on. The output of this operation is provided in FIG. 3 A. Consecutive values of the digitised signal values are assigned to timeslots in step 208. Hence, one timeslot preferably comprises multiple signal values. Consecutive timeslots may overlap, such that at last some signal values may be assigned to two or more consecutive timeslots.

[0055] Of the various signal values, the amplitude envelope is subsequently determined in step 210. The output of this operation is provided in FIG. 3 B. The amplitude envelope is determined over all signal values or at least a substantial part thereof, irrespective of the timeslots. The amplitude envelope may be determined by means of an absolute Hilbert transform, but also, additionally or alternatively, by means of other methods. Additionally or alternatively, an envelope of the power of the signal is determined.

[0056] Subsequently, the process splits in two branches. In the left branch, the signal envelope is integrated in step 212. Alternatively or additionally, signal values are summed in another way. Such summing may be preceded by a demeaning step as discussed above. The demeaning and/or the summing may be processed over the whole time period or over a timeslot. The output of this operation is provided in FIG. 3 C. Such integration may be executed by simple arithmetic adding or specific integration algorithms. In another example, the signal is provided to a low-pass filter.

[0057] The values of the integrated signal are, on a per-timeslot basis, normalised in step 214. The normalisation may be executed based on an average value per timeslot. The output of this operation is provided in FIG. 3 D. As this is preferably the first operation that is executed on a per-timeslot basis, the signal values may alternatively, also be assigned to timeslots between step 214 and step 216.

[0058] The signal values thus normalised are, on a per-timeslot basis, detrended in step 216. Various algorithms are available for detrending a signal. The output of this operation is provided in FIG. 3 E.

[0059] Of the integrated, normalised and detrended values of the determined signal envelope, the variation is calculated in step 218. The output of this operation is provided in FIG. 3 F. Such variation may be the standard deviation. Alternatively or additionally, also other algorithms for determining a variation of a parameter may be employed.

[0060] In the right branch, a mean or median amplitude or amplitude envelope is calculated for values assigned to a particular timeslot in step 220. Depending on the application, either operationmean or medianmay be calculated. In further examples, outliers may be removed from a dataset before calculating the mean or median value. Alternatively or additionally the mean or median of the power or power envelope may be calculated per timeslot.

[0061] In step 222, a correlation between the variation per timeslot and the mean value calculated in step 220 is determined over the consecutive timeslots in which the obtained signal is divided. The output of this operation is provided in FIG. 3 G. The correlation value is subsequently compared to a first threshold value in step 224. The value of the threshold will be discussed below. If the correlation is lower than the first threshold, the procedure branches to the right and it is determined in step 240 that the ratio of excitation and inhibition signalling is inhibition dominant.

[0062] If the correlation value is lower than the first threshold, the procedure proceeds to step 226. In step 226, the correlation value is compared to a second threshold value. If the correlation value is lower than the second threshold, the process branches to the right and it is determined in step 242 that the ratio of excitation and inhibition is excitation dominant.

[0063] If the correlation value is higher than the first threshold and lower than the second threshold, the process continues to step 228 in which it is determined that the ratio of excitation and inhibition is balanced.

[0064] In an alternative implementation, the comparison step 224 and the comparison step 226 are combined in one step. In such single comparison step, the correlation value is compared to one single third threshold value.

[0065] It has been determined that the correlation between the mean or median of the power or power envelope per slot and the variation of the signal and more in particular of the signal power or the signal power envelope (or signal amplitude envelop) provides an estimation of the ratio of excitation and inhibition in a part of the brain of which the signal has been obtained. For obtaining such correlation value, one or more operations on the obtained signal data may be required, as discussed above. With 1 minus the correlation being below 1, the ratio of excitation and inhibition is inhibition dominant. And with 1 minus the correlation being above 1, the ratio of excitation and inhibition is excitation dominant. Hence, such third threshold may be set at 1; if 1 minus the correlation has been determined to be substantially equal to 1, the ratio of excitation and inhibition is balanced. In such case, there is no correlation between the two datasets.

[0066] In case of using two threshold values, as depicted by FIG. 2, the lower first threshold may be set at 0.99, 0.95, 0.9, 0.85 or 0.8 and the higher second threshold may be set at 1.01, 1.05, 1.1, 1.15 or 1.2. Alternatively or additionally, not the value of 1 minus the correlation is assessed, but the correlation may be assessed. In such case, the lower first threshold may be set at 0.01, 0.05, 0.1, 0.15 or 0.2 and the higher second threshold may be set at 0.01, 0.05, 0.1, 0.15 or 0.2.

[0067] Having determined the ratio of excitation and inhibition signalling in a population of neurons, the process ends in a terminator 230.

[0068] In an alternative implementation, further factors may be taken into account for determining whether the ratio of excitation and inhibition is either inhibition or excitation dominatedor balanced. For example, the DFA value may be taken into account. If the DFA is found to be less than a pre-determined DFA value and the received signal has no major oscillations, i.e. the median amplitude or power envelope is relative low, it may be determined that the ratio of excitation and inhibition is inhibition dominant. For such determination, the pre-determined DFA value may be set at any value between 0.5 and 0.65 and preferably at 0.6. If the received signal is determined to oscillate significantly, the ratio of excitation and inhibition is determined to be excitation dominant. Alternatively, as an option, if the signal is smaller than the pre-determined DFA value, the ratio of excitation to inhibition is not estimated.

[0069] In the description above, it will be understood that when an element such as layer, region or substrate is referred to as being on or onto another element, the element is either directly on the other element, or intervening elements may also be present. Also, it will be understood that the values given in the description above, are given by way of example and that other values may be possible and/or may be strived for.

[0070] Furthermore, the invention may also be embodied with less components than provided in the embodiments described here, wherein one component carries out multiple functions. Just as well may the invention be embodied using more elements than depicted in the Figures, wherein functions carried out by one component in the embodiment provided are distributed over multiple components.

[0071] It is to be noted that the figures are only schematic representations of embodiments of the invention that are given by way of non-limiting examples. For the purpose of clarity and a concise description, features are described herein as part of the same or separate embodiments, however, it will be appreciated that the scope of the invention may include embodiments having combinations of all or some of the features described. The word comprising does not exclude the presence of other features or steps than those listed in a claim. Furthermore, the words a and an shall not be construed as limited to only one, but instead are used to mean at least one, and do not exclude a plurality.

[0072] A person skilled in the art will readily appreciate that various parameters and values thereof disclosed in the description may be modified and that various embodiments disclosed and/or claimed may be combined without departing from the scope of the invention.

[0073] It is stipulated that the reference signs in the claims do not limit the scope of the claims, but are merely inserted to enhance the legibility of the claims.