EP3445239A2 - Method and system for estimating a location of an epileptogenic zone of a mammalian brain - Google Patents
Method and system for estimating a location of an epileptogenic zone of a mammalian brainInfo
- Publication number
- EP3445239A2 EP3445239A2 EP17720722.2A EP17720722A EP3445239A2 EP 3445239 A2 EP3445239 A2 EP 3445239A2 EP 17720722 A EP17720722 A EP 17720722A EP 3445239 A2 EP3445239 A2 EP 3445239A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- zone
- brain
- nodes
- model
- location
- 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
Links
- 230000002397 epileptogenic effect Effects 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 77
- 210000005171 mammalian brain Anatomy 0.000 title claims abstract description 22
- 210000004556 brain Anatomy 0.000 claims description 130
- 230000001537 neural effect Effects 0.000 claims description 59
- 230000000694 effects Effects 0.000 claims description 40
- 206010001497 Agitation Diseases 0.000 claims description 33
- 230000003902 lesion Effects 0.000 claims description 28
- 230000008878 coupling Effects 0.000 claims description 23
- 238000010168 coupling process Methods 0.000 claims description 23
- 238000005859 coupling reaction Methods 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 18
- 238000004088 simulation Methods 0.000 claims description 18
- 238000002610 neuroimaging Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 3
- 230000036962 time dependent Effects 0.000 claims 1
- 206010010904 Convulsion Diseases 0.000 description 65
- 238000002595 magnetic resonance imaging Methods 0.000 description 25
- 238000003062 neural network model Methods 0.000 description 15
- 230000001709 ictal effect Effects 0.000 description 13
- 238000002513 implantation Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 8
- 230000006399 behavior Effects 0.000 description 7
- 230000007177 brain activity Effects 0.000 description 7
- 230000002109 interictal effect Effects 0.000 description 7
- 238000013459 approach Methods 0.000 description 6
- 238000002271 resection Methods 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 210000001320 hippocampus Anatomy 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 208000015134 congenital hypothalamic hamartoma syndrome Diseases 0.000 description 4
- 238000005265 energy consumption Methods 0.000 description 4
- 230000003190 augmentative effect Effects 0.000 description 3
- 230000001054 cortical effect Effects 0.000 description 3
- 230000001934 delay Effects 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 238000013515 script Methods 0.000 description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 description 3
- 230000000287 tissue oxygenation Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 2
- 210000003016 hypothalamus Anatomy 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 230000036244 malformation Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002503 metabolic effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000007115 recruitment Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 210000003625 skull Anatomy 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 210000003478 temporal lobe Anatomy 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 206010012177 Deja vu Diseases 0.000 description 1
- 206010018852 Haematoma Diseases 0.000 description 1
- 208000004547 Hallucinations Diseases 0.000 description 1
- 208000002927 Hamartoma Diseases 0.000 description 1
- 206010048911 Lissencephaly Diseases 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 206010048910 Pachygyria Diseases 0.000 description 1
- 206010033546 Pallor Diseases 0.000 description 1
- 206010033557 Palpitations Diseases 0.000 description 1
- 208000037158 Partial Epilepsies Diseases 0.000 description 1
- 206010061334 Partial seizures Diseases 0.000 description 1
- 206010035039 Piloerection Diseases 0.000 description 1
- 230000007321 biological mechanism Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000090 biomarker Substances 0.000 description 1
- 230000006931 brain damage Effects 0.000 description 1
- 231100000874 brain damage Toxicity 0.000 description 1
- 208000029028 brain injury Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 230000001055 chewing effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 206010013395 disorientation Diseases 0.000 description 1
- 206010015037 epilepsy Diseases 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 201000007186 focal epilepsy Diseases 0.000 description 1
- 210000001652 frontal lobe Anatomy 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000007917 intracranial administration Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000003957 neurotransmitter release Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 210000000869 occipital lobe Anatomy 0.000 description 1
- 210000001152 parietal lobe Anatomy 0.000 description 1
- 230000005371 pilomotor reflex Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 201000009570 retrograde amnesia Diseases 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012421 spiking Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000002739 subcortical effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000000472 traumatic effect Effects 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
- A61B5/4094—Diagnosing or monitoring seizure diseases, e.g. epilepsy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/291—Bioelectric electrodes therefor specially adapted for particular uses for electroencephalography [EEG]
- A61B5/293—Invasive
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6846—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive
- A61B5/6867—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be brought in contact with an internal body part, i.e. invasive specially adapted to be attached or implanted in a specific body part
- A61B5/6868—Brain
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/291—Bioelectric electrodes therefor specially adapted for particular uses for electroencephalography [EEG]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
Definitions
- This application discloses a method for estimating a location of an epileptogenic zone of a mammalian brain optimizing the placement of stereotactic electrodes, a system for performing said method, and a computer-readable medium containing instructions for performing such a method, as outlined by the claims and the content of this description.
- Personalized medicine proposes the customization of healthcare with medical decisions, practices, and products being tailored to the individual patient. Inter-individual variability between different patients has clear effects upon the responsiveness to treatment approaches.
- Patients suffering from brain seizures may receive drug treatment.
- a possible treatment for patients having, for example, a drug-resistant form of seizures is the surgical resection of the epileptogenic zone.
- An epileptogenic zone is a localized region or network within the patient's brain where seizures arise, before recruiting secondary brain networks called the propagation zone.
- the propagation zone includes areas of the patient's brain that are affected by, for example, a seizure arising in the epileptogenic zone, the seizure then spreading into the areas encompassing the propagation zone.
- SEEG stereotactic electroencephalograms
- MRI structural magnetic resonance imaging
- EEG classical electroencephalograms
- MEG magnetoencephalograms
- CT computed tomography
- SPECT single-photon emission computed tomography
- the implantation of electrodes into the mammalian brain for recording SEEGs is a high-risk procedure and should only be performed when the professional tasked with the procedure is quite certain about where said electrodes are most likely to correctly identify the epileptogenic zone.
- a method for estimating a location of an epileptogenic zone of a mammalian brain may include receiving a structural skeleton model of a mammalian brain, wherein the structural skeleton model comprises a plurality of nodes and is based on non-invasive neuroimag- ing data and wherein connectivity information of the brain between different nodes is extracted from the non-invasive neuroimaging data.
- a coupled brain network model is provided by populating each node of the structural skeleton model with a neural population model, where- in a neural population model corresponding to a node is coupled to further neural population models corresponding to further nodes according to the connectivity information extracted from the non-invasive neuroimaging data.
- a propagation zone is predicted in the brain network model based on an evolution of a simulated seizure starting from the first estimate of the location of the epileptogenic zone of the coupled brain network model.
- any estimate of the location of the epileptogenic zone includes only a real subset of all the nodes of the brain network model. Furthermore, it is often assumed that the predicted propagation zone also does not include all nodes of the brain network model, but only includes a real subset of the nodes.
- a structural skeleton model may be provided by importing structural and dMRI data of a mammalian brain via an input/output interface into a memory unit or a central processing unit of a computer system.
- data may include raw data, pre-processed data, or processed data.
- Raw structural or dMRI data of the patient's brain may be processed, for example, by using Scripts
- Structural and dMRI data may be used for building up a structural skeleton model specific to the particular patient's brain. Thus, the structural and dMRI data need to be taken from the specific brain, in which the epileptogenic zone is to be identified.
- a structural skeleton model derived from structural and dMRI data is divided into a plurality of voxels, wherein a single node may be assigned to a single voxel or a plurality of voxels.
- a typical amount of nodes in a structural skeleton model of the human brain may add up to thousands of nodes.
- Each node may represent a volume of the mammalian brain either located on the cortical surface, within one of a plurality of cortical regions or within a plurality of subcortical regions.
- the parcellation may be performed, for example, using the Desikan-Killiany atlas.
- tractography may be performed and a connectivity matrix connecting the different nodes of the structural skeleton model may be obtained, for example by summing track counts over each region of the parcellation and by normalizing the values so that the maximum value of the connectivity matrix is 1.
- the method further provides a coupled brain network model such as a large- scale brain network model.
- a coupled brain network model such as a large- scale brain network model.
- These large-scale brain network models are known to provide insights into the mechanisms underlying the emergence of the resting-state network dynamics.
- a coupled brain network model is provided by populating each node of the structural skeleton model with a neural population model or neural mass model.
- different nodes can be provided with different neural population models.
- nodes of the structural skeleton model are populated with the same neural population model; furthermore different nodes may include different parameters for the neural population model, for example, at least a subset of nodes includes a value of an excitability parameter which allows the neural population model to reach an ictal state, while other nodes include a value of the excitability parameter which does not push the neural population model into an ictal state without external input from other nodes.
- Neural population models are often represented by a system of coupled differential equations.
- One such example is the Epileptor as described, for example, in Jirsa et al., "On the nature of seizure dynamics", Brain, vol.
- the neural population model can be chosen to exhibit a bistable behaviour, such as a saddle-node bifurcation or subcritical Hopf bifurcation, to enter the ictal state.
- the same neural population model can also show a bistable behaviour, such as a saddle homoclinic bifurcation or fold limit cycle bifurcation, to exit the ictal state, i.e. to enter a non-ictal state.
- Each neural population model is coupled to further neural population models representing a different node in the structural skeleton model.
- the coupled neural population models are coupled so as to represent the connectivity information extracted from the non-invasive neuroim- aging data.
- the neural population model may include state variables evolving over time and parameters, which govern the behaviour of the state variables over time.
- a neural population model parameter may govern a degree of excitability of the model, i.e. an increase (or decrease) leads to a model more prone to exhibit autonomous seizure (or spiking) activity.
- Other parameters may govern the strength of a connection between two neural population models placed at different, yet connected nodes.
- the neural population model may include two or more coupled differential equations.
- the system of coupled differential equations may include a fast and a slow subsystem. This means that the two differential equations operate on different time scales.
- a biological correlate of a fast subsystem may be local neuronal activi- ty in the respective brain region.
- the slow subsystem may represent modulations of connections between different brain regions and their corresponding biological mechanisms such as neurotransmitter release, metabolic processes, or other, slower mechanisms forwarding signals through different brain networks.
- a programmer inputting an estimated location of an epileptogenic zone by learned experience
- a computer system inputting an estimated location of an epileptogenic zone based on algorithms
- Such an estimate may include at least one of the plurality of nodes.
- the neural population model is provided with a first set of parameters, the set of parameters known to be able to simulate or "trigger" a seizure.
- Triggering a seizure is, in some examples of the method, equivalent to a parameter set of a neural population model or a neural population model triggering seizures autonomously, i.e., the neural population model does not require external input into the neural population model to exhibit fast oscillations representing an ictal state occurring during a seizure.
- the nodes not identified as the first estimate of the location of the epileptogenic zone may be provided with a first set of parameters such that the neural population model does not trigger a seizure autonomously.
- the neural population models of areas in a non-epileptogenic zone require external input to exhibit fast oscillations representing an ictal state of a seizure. External input is understood as input from other nodes.
- a propagation zone may arise, i.e. a simulation of the brain network model may result in a propagation pattern of activity representing a location of a propagation zone in the brain network model.
- an epileptogenic zone exhibiting autonomous seizures may acquire further nodes through the coupling of the neural population models, which may show seizure activity in response to an autonomous seizure in the epileptogenic zone.
- Such a time evolution may be simulated by, for example, a system including a central processing unit, a memory unit, an input/output unit, or a computer-readable medium containing instructions for performing the method described above in connection with the system as outlined described above.
- the method repeatedly estimates an epileptogenic zone and predicts a resulting propagation zone. If the resulting propagation zone does not match an observed propagation zone, or if a clinician, using his experience comes to the conclusion that the resulting propagation zone is unlikely, a new hypothesis of an epileptogenic zone is chosen and a new, resulting propagation zone is determined.
- parameters of the neural population models may be changed. A change in the parameter values may lead to a different predicted propagation zone than the propagation zone determined using a first set of parameter values, and the new propagation zone may resemble previously recorded patient seizure or interictal data.
- the neural population model includes a parameter representing an excitability of the neural population model.
- Neural population models belonging to a node or a plurality of nodes of the epileptogenic zone may be assigned parameter values, such that the neural population models located in the estimated epileptogenic zone exhibit autonomous seizure activity.
- the model may also exhibit interictal and preictal activity.
- Other nodes are also assigned parameter values of the parameter representing a degree of excitability, however, the parameter values are chosen such that the neural model shows less excitability.
- the distribution of the parameter values is based on the distance from the epileptogenic zone.
- a distance can hereby be understood as a spatial distance or a time distance, i.e. the further away a node is from the epileptogenic zone, the lesser the excitability is chosen to be. Whether different nodes are connected may also influence the distribution of the parameter values of the parameter representing excitability.
- the state variables of neural population models may represent biological parameters such as energy consumption, tissue oxygenation, or extracellular ion concentrations, or as previously mentioned a degree of excitability
- the time evolution of simulation data gathered from the coupled brain network model may show a propagation zone that may or may not resemble seizure activity in a real patient, in particular in the patient from which the structural skeleton model was derived.
- a further or new prospective epileptogenic zone can be chosen and a corresponding propagation zone can be predicted.
- An exemplary procedure is outlined in this application.
- the probability of identifying the correct epileptogenic zone is highly increased in contrast to empirical or heuristic approaches.
- a second estimate of the epileptogenic zone replaces the first estimate of the location of the epileptogenic zone if the simulated propagation zone of the first estimate differs from an observed propagation zone.
- a wrongly identified epileptogenic zone can be falsified by simulation.
- a simulated activity or signal may be compared to a real, recorded brain signal by standard statistical methods. If the simulated activity based on a first epileptogenic zone resembles the recorded signal well, or better than simulated activity from other epileptogenic zone estimates, the first epileptogenic zone can be kept and the other zones can be discarded.
- the propagation zone prediction may include the prediction of electrical activity data by employing a forward model.
- Forward models are used to map, for example, electrical activity or tissue oxygenation occurring within the brain volume at different times in different places to surface potentials, which may be measured by invasive or noninvasive surface electrodes during the recording of a SEEG, EEG or MEG. Since the analysis of SEEG and EEG data, in particular the comparison of different EEG and SEEG recordings, is easily accomplished by a plurality of commercial and non-commercial tools, simulated seizures can easily be compared to recorded seizure patterns. In particular, the implantation of electrodes can be postponed until a number of potential epileptogenic zones have been ruled out by the method as suggested in this application. Furthermore, the implantation scheme can be optimized with regard to best coverage and minimal invasiveness. This significantly improves the chances of performing a successful surgical resection.
- the proposed position of the SEEG electrodes can be mapped to a three-dimensional model of the brain, to help a practicioner during the implantation.
- the location of the proposed position can be shown in an augmented reality device, such as augmented reality glasses.
- the position of the SEEG electrodes is forwarded to the augmented reality software which includes a three dimensional model of the patient's brain.
- the coupled brain network model may be further adapted by including a parameter representing a structural anomaly, such as an MRI lesion or a malformation such as pachygyria or an hematoma, in at least one node.
- a structural anomaly such as an MRI lesion or a malformation such as pachygyria or an hematoma
- structural anamolies e.g. MRI lesions
- MRI lesions are known to have an effect on epileptic seizures and are a very good indicator for the location of the epileptogenic zone.
- structural anamolies e.g. MRI lesions
- An MRI lesion or other structural anomalies are important structural data further informing and improving the model.
- an MRI lesion may be identified by characteristic dark patterns by a software or a practicioner when analyzing patient specific neuroimaging data, such as MRI images. While the structural skeleton model may result in showing that two nodes are connect- ed with each other, an MRI lesion may add additional information to the brain network model in that the strength of the connection between nodes located in the area of the MRI lesions and nodes neighboring those nodes needs to be modified. However, in other examples the structural anomalies may represent a change in the dynamics of the node itself. By including MRI lesion information, the quality of the prediction and subsequent correct identification of the epileptogenic zone is greatly enhanced.
- the parameter values of the parameter representing the structural anomaly are distributed throughout the different nodes of the brain network model, such that the distribution resembles the occurrence of a structural anomaly in the mammalian brain.
- Fig. 1 is a schematic overview of a system for estimating the location of an epileptogenic zone
- Fig. 2 illustrates the relation between different ways of understanding neural activity
- Figs. 3 are illustrations of a structural skeleton model and a coupled brain network model
- Figs. 4 show a flow chart of an exemplary method for finding an epileptogenic zone using the brain network model.
- Fig. 1 shows a schematic overview of a system configured for estimating the epileptogenic zone in a patient.
- System 10 includes a central processing unit (CPU) 20, a memory unit 30, an input/output interface 40, and several input/output devices coupled to the interface 40.
- the system can be a computer or a distributed system, in which the different units may include further components to function independent of each other, but are coupled via different data connections.
- the CPU may include a multiprocessor and/or multi-core processor for fast algorithm and data execution.
- the memory unit 30 may include primary storage linked to the CPU(s), random access memory (RAM), volatile memory as well as hard disks, flash memory units, or EEPROM units for storing data such as a structural skeleton 32, the coupled brain network model 34, various implementations of a neural population model 36 to placed in different nodes of the structural skeleton, and/or instructions 38 for executing different algorithms when processing and/or comparing the data.
- the input/output interface 40 can include several interfaces such as an input for a keyboard, a wired or wireless data connection for uploading or downloading data into the memory unit, or interfaces for connecting a display, such as a monitor, or a keyboard for inputting commands.
- the interface of Fig. 1 is coupled with the internet 60, a display 70, a keyboard 80 and a flash memory reader 90. Additionally, the interface 40 may be coupled to an EEG displaying unit 100.
- Fig. 2 illustrates how different brain signals can be recorded and how different imaging methods are connected to each other.
- Fig. 2 shows an exemplary scalp EEG recording 210 including different channels, which include different signals, such as ictal periods 212 or 214 and non-ictal or interictal period 216.
- EEG recordings have a high temporal resolution, but provide little information on the spatial distribution of recorded electrical signals throughout the brain.
- Neuroimaging methods such as MRI (represented by longitudinal section through a skull 220) or methods based on MRI provide a much better spatial resolution.
- the temporal resolution of MRI based methods is limited as the signal recorded in MRI is not the electrical brain activity, but a signal representing the energy consumption.
- Modern techniques may link the seizure activity 212 or 214 to increased metabolic activity in area 222 or 224, respectively.
- MRI techniques include dMRI techniques, which allow a non-invasive extraction of connectivity information of a patient.
- the connectivity information is displayed in matrices to establish a time-space structure of the coupling (for further information, see for example Jirsa, V.K. "Neural field dynamics with local and global connectivity and time delay.” Philos. Trans. A. Math. Phys. Eng. Sci. 367, 1131-1143, which is incorporated herewith in its entirety).
- the actual brain activity which underlies the EEG recording 210 and MRI 220, is based on the electrical activity of different brain areas, which are connected with each other.
- the brain activity can be visualized by different nodes representing different areas of the brain, which are connected to each other according to the connectivity information.
- the resulting network (represented by reference numeral 230) dis- plays localized, electric activity.
- the electric (and metabolic) energy needed for generating the electric brain activity can be seen in the MRI images as energy consumption. Furthermore, the localized, electric activity can be recorded by EEG electrodes.
- Techniques for mapping localized electric brain activity in different brain areas to EEG-like signals are generally known as forward models 240. Techniques for mapping EEG signals to localized, electric activity in different brain areas include inverse models 250 and usually much harder to describe than forward models or solutions (see for example, Jirsa, V.K. et al "Spatiotemporal forward solution of the EEG and MEG using network model- ing", IEEE Trans. Med. Imaging 21, 493-504.).
- the forward models map activity in the simulated three-dimensional physical brain (represented by the nodes in different positions in three-dimensional space) to the surface at which the EEG signals are to be recorded.
- Fig. 3A and Fig. 3B illustrate the structural brain skeleton and the coupled brain network model.
- the structural brain skeleton of Fig. 3A is an individualized or patient-specific representation of large-scale brain connectivity.
- the structural skeleton model of the patient represents the patient's individual pattern of connections between different brain areas.
- the structural skeleton model may include information on the time it takes electric activity of one brain area to effect electric activity in another area.
- a structural skeleton often involves a panellation of a patient's brain in voxels.
- the number of voxels (see an exemplary voxel 299) of a typical structural skeleton model may include between 100 and 100000 voxels, preferably between 5000 and 50000 voxels.
- the structural skeleton may include a representation of the patient's cortical surface only, or include further brain areas such as parts of the temporal gyrus, the frontal gyrus, the frontal lobe, the temporal pole, the occipital lobe, the parietal lobe.
- the dMRI data is used to extract the connections between different voxels of different brain areas.
- the neural network model will include a sys- tern of coupled differential or difference equations.
- the neural network model used is the Epileptor as described in Jirsa 2014. The Epileptor includes five state variables acting on three different time scales. On the fastest time scale, state variables xl and yl account for fast discharges during a seizure.
- the permittivity state variable z accounts for slow processes such as a variation in extracellular ion concentrations, energy consumption, and/or tissue oxygenation.
- the system exhibits fast oscillations during the ictal state through the variables xl and yl.
- Autonomous switching between interictal and ictal states is realized via the permittivity variable z through saddle-node and homoclinic bifurcation mechanisms (i.e. bistable behaviour) for the seizure onset and offset, respectively.
- the switching is accompanied by a direct current (DC) shift, which has been recorded in vitro and in vivo (see for example, Jirsa 2014).
- DC direct current
- x2 -yl + x2 - x2 3 + 12 + 0.0002#(xl) - 0.3(z - 3.5)
- y2 —(-y2 + f2(xl, x2))
- the parameter xO controls the tissue excitability, and is epileptogenic, i.e. is triggering seizures autonomously for a critical value xO >—2.05, otherwise the CodeBox Computerdienste GmbH
- 177PCT 0962 14 tissue is healthy.
- the five dimensional epileptor can be reduced to a two-dimensional system:
- the fast variable of the neural network model represents fast discharges
- the slow variable represents the switching between ictal and interictal states through a bifurcation of the dynamic system of the neural network model.
- the different nodes of the brain network model are coupled by permittivity coupling, i.e. the neural network models of connected nodes are coupled in their slow variables.
- K y includes the connectome Q (a value representing the connectivity between different voxels or nodes), and a scaling factor G, wherein ⁇ ⁇ ⁇ ⁇ x G.
- the values for Ci j are determined from dMRI or DTI data.
- the index j runs over all nodes connected to node i.
- the equation fort he slow variable then reads, for example
- the permittivity coupling from node i to node j can also be chosen to include a signal transmission delay to account for real brain transmission delays.
- the delay is introduced by modifying the coupling term such that the input from node j to node i is delayed by a time tdeiay, > e. Kjj(xl](t)-xlj(t-tdeiay)-
- the brain network model can be used for estimating the location of an epileptogenic zone using the brain network model.
- said nodes neural network model i.e.
- the neural network models of a first number of nodes are provided with a parameter set, which allows the neural network models to exhibit autonomous seizures.
- the parameter set may include values for K, j , II, 12, xO and optionally others.
- the parameter values of a parameter representing excitability can be distributed by a Gaussian distribution, i.e. the parameter values being highest in the epileptogenic zone and being smaller the further away the corresponding node is from the epileptogenic zone. Changing the distribution scheme may consequently result in different propagation zones. Distance between different nodes can also be based on the strength of connections between different nodes in the brain network model.
- the SEEG data are windowed and Fourier transformed to obtain estimates of their spectral density over time. Then SEEG power above 10 Hz is summed to capture the temporal variation of the fast activity. These time series are corrected to a pre-ictal baseline, log-transformed and linearly detrended over the time window encompassing the seizure. Contacts are selected, which present greater high-frequency activity than their neighbors on the same electrode. Given that, contrary to M/EEG, the SEEG lead field is very sparse, three nodes per contact are used in the network model. Other nodes are not recruited and rest at their fixed points. The effect is approximated by the fitting through a constant sum over the corresponding elements of the structural connectivity matrix. Next, one may use an observation model that incorporates an SEEG forward model, under the assumption that the variable describes fluctuations in the log power of high frequency activity, predicting sensor log power, with normally distributed observation error.
- Uninformative priors are placed on the hidden states' initial conditions, while their evolution follows an Euler-Maruyama discretization of the corresponding stochastic differential equations with linear additive normally distributed noise. Uninformative priors are also placed on the excitability parameter per node, observation baseline power, scale and noise. Finally, the length of the seizure is also allowed to freely vary to match that of a given recorded seizure. Structural connectivity specifies a prior on the connectivity used in the generative method. This model is implemented using Stan, a software for Bayesian inference, which implements both Hamiltonian Monte-Carlo and automatic variational inference algorithms for generic differential probability models (see Hoffman, M.D., Gelman, A.
- an Euler- Maruyama integration scheme with an exemplary integration step of 0.05 may be used.
- Additive white Gaussian noise is introduced in the variables and with mean 0 and variance 0.0025 (see Jirsa 2014).
- Other variables experienced only little or no noise due to their high sensitivity.
- 256 time steps are equivalent to one second of real time to obtain realistic frequency ranges, seizure lengths, and matched intracranial EEG sampling frequency.
- a rectangular function in time was applied on the z variable of the stimulated region (amplitude: 0.5, length : 2 s).
- the result of the above fitting procedure is a parameter value for each node, thereby establishing a distribution (spatial, temporal or spatio-temporal) for each node of the brain network model).
- a map of the distribution may also be referred to as a heat map.
- the choice of the first estimate of an epileptogenic zone may be provided by a clinician, followed by an iterative process between clinician and the method disclosed herein or a fully automated approach, resulting in a refinement of the epileptogenic zone via a second, third or fourth estimate, which may be subsequently provided; however, the systems and methods as disclosed in this application may also include an automated choice simulation of different first, second and further epileptogenic zones to arrive at a preferred epileptogenic zone, the propagation zone, of which resembles previously recorded brain activity best.
- the relationship between the brain network model, an estimate of an epileptogenic zone, a propagation zone, and other regions is illustrated in Fig. 3B.
- the (simplified due to a reduced number of nodes) brain network model 300 includes a plurality of nodes 310. Each node includes a neural network model such as the epileptor. The different nodes are coupled to further nodes according to the connectivity information, for example the connectome Cy.
- the coupling between nodes i and j can be uni- or bi-directional (Ci j ⁇ 0,C j i ⁇ 0), represented by uni-directional or bi-directional arrows, respectively.
- An epileptogenic zone is first estimated to include two nodes 312 and 314, i.e.
- the neural network models of each node of the estimated epileptogenic zone can exhibit autonomous seizures.
- the resulting propagation zone includes nodes which are coupled to the nodes 312 and 314, for example nodes 322, 324 and 326.
- the nodes of the propagation zone while not able to exhibit seizures autonomously, are captured by the nodes of the epileptogenic zone, i.e. receive input from the nodes of the epileptogenic zone through the permittivity coupling and are driven to exhibit seizure behavior as well.
- the remaining nodes of the brain network are not captured by the epilepto- genie zone in the sense that the nodes do not receive sufficient input from other nodes to be driven into a seizure state.
- Further information may include structural anomalies.
- an MRI lesion may be included by modifying the local connectivity of the area involved. This may be achieved by changing the scalar factor G to G /oco// with G ⁇ Gi ocatl to represent the structural anomaly.
- the local connectivity topology of the area i.e. the nodes affected by the structural anomaly, is not changed but scaled up by the factor 6 /0C0/ .
- anomalies which may be modeled can include pachygaria, hamartoma and others.
- Including a value Gi oca i introduces a lesion map into the brain network model.
- structural anomalies may be identified by a practicioner or a software from neuroimaging data, such as MRI images.
- the areas of the patient's brain to be modeled which appear to show a structural anomaly are mapped to the nodes representing the respective area in the network model. While the area including the identified anomaly can in some example include a single node, said area can encompass several nodes, i.e. a subset A of nodes.
- the nodes included in the subset A can be identified through streamlines from DTI measurements, for example.
- the DTI measurement cannot be used for introducing a lesion map, as the DTI streamlines do not indicate the strength of a coupling between areas, but merely the existence of a coupling between areas.
- G ioC ai is a scalar factor to the global G value introduced earlier.
- each subset is assigned a value G ioca i > 1, i.e. the local value of G for the nodes of the subsets, A, B, ... is larger than the global value of G .
- different subsets A and B ca n be assigned different values of G loca to represent different anomalies or different anomaly effects.
- the values of the G local of all nodes i, i.e. Gi ocaXi i completes the initial lesion map.
- This lesion map is based on structural anomaly data and can, in some embodiments, be also used in the fitting or estimation scheme discussed in this application, i.e. the initial lesion map can be changed to acquire a better estimate of the epileptogenic zone of the specific patient's brain.
- the personalized brain network model of a patient may therefore include a connectome represented by C i; - (derived from DTI), a lesion map represented by G local i (derived from structural anomalies or malformations of the patient a scaling the global G value for each node) and a heat map or distribution of epileptogenicity or excitability represented by e.g. a map of values xO for the epileptor neural population model.
- the lesion map can remain the initial lesion map and only the distribution epileptogenicity or excitability is changed to find the most likely epileptogenic zone for the patient, or, in other examples, the initial lesion map is changed while the distribution of excitability (i.e.
- both the lesion map and the heat map are changed to find the most likely candidate for the epileptogenic zone and thereby an adequate candidate area for a implantation site of one or more SEEG electrodes or a lesion site for surgical resectioning.
- the brain network model When using the brain network model to find the propagation zone of an estimate of an epileptogenic zone, it has proven helpful, to translate the brain network model activity, using a forward model, into SEEG, fM RI, or EEG activi- ty. Since SEEG, fM RI and EEG data (i.e. real, not simulated patient data) can be easily obtained from a patient, the translated brain network model activity can be easily compared to rea l EEG patient data to access whether the estimate of the epileptogenic zone provides brain activity which is similar to measured brain data, or whether a new epileptogenic zone must be chosen to arrive at a better fit between the simulated brain network model activity and the measured brain data.
- SEEG, fM RI and EEG data i.e. real, not simulated patient data
- the systems and methods discussed in this application may also encompass the uploading, reception or generation of a patient's brain network model based on a patient's structural skeleton model as indicated in Fig. 4.
- a patient's brain network model 400 has been generated or uploaded into a system such as the system of Fig. 1 and real brain data 410 (such as EEG 412, fMRI 414, SPECT and/or SEEG data), from the patient including a period of a simple seizure (i.e. a seizure not capturing a propagation zone) or a complex seizure (i.e.
- an estimation method 420 is initiated which includes: choosing a first estimate of an epileptogenic zone 430; a subsequent simulation of a resulting propagation zone 440 and an optional translation 450 via a forward model into simulated EEG, fMRI, SPECT or SEEG activity 460; a comparison (step 470)between the measured, real brain data 410 and the simu- lated activity 460.
- the estimation method may optionally include using prior information 435 to obtain an estimate of a first epileptogenic zone including at least one node such as a clinician's estimate of the first epileptogenic zone or an automated starting choice.
- the epileptogenic zone may be varied after generating the simulated activity 460 for a specific epileptogenic zone and steps 430 through 470 may be repeated with a new estimate of a location of the epileptogenic zone.
- the varied epileptogenic zone estimate may include neighboring nodes, additional nodes, or a subset of nodes from the previous or first estimate of the epileptogenic zone.
- the estimation method may then include a decision engine 480 (applying statistical methods as known in the prior art to compare time-series; the decision engine may be a software module for deciding which simulated time-series fits the real data best) to find the best match 485 between simulated brain network model activity and real brain data and output the epileptogenic zone of the best matching data 490 as a suggestion for a brain area of surgical resection or a suggestion for a pre-surgical implantation site of electrodes for recording
- a decision engine 480 applying statistical methods as known in the prior art to compare time-series; the decision engine may be a software module for deciding which simulated time-series fits the real data best) to find the best match 485 between simulated brain network model activity and real brain data and output the epileptogenic zone of the best matching data 490 as a suggestion for a brain area of surgical resection or a suggestion for a pre-surgical implantation site of electrodes for recording
- SEEG data to minimize the traumatic impact for a patient.
- a clinician obtains non-invasive brain images of the patient including MRI, DTI, EEG, MEG.
- the clinician interprets all data and makes his/her initial hypothesis on the location of the EZ and a proposal of an implantation scheme for the SEEG electrodes.
- a first virtual patient brain model i.e. brain network model
- the Clinician's first hypothesis on EZ is used as a first estimate of the location of the epileptogenic zone.
- the computer sys- tern generates simulated imaging data (EEG, MEG), which will be analysed using the clinical standard visualization (visual inspection of EEG time series by expert eye) and analysis tools (biomarkers such as Functional Connectivity, Epileptogenicity index, H2).
- EEG simulated imaging data
- MEG simulated imaging data
- analysis tools biomarkers such as Functional Connectivity, Epileptogenicity index, H2
- EZ hypothesis i.e. by changing the spatial distribution of excitability values, as well as including MRI lesion information in the brain network model
- the clinician will attempt to converge model behavior with empirical patient data. This iterative process is close to the spirit of contemporary patient management staff meetings in hospitals, where individual patient cases are discussed.
- model parameters in the brain network model will be fitted automatically against empirical EEG/MEG data.
- the computer system Once the updated location of the epileptogenic zone is deemed suffi- cient, the computer system generates forward solutions of SEEG data and proposes an optimal implantation scheme for SEEG electrodes.
- the method outlined by Fig. 4 and its corresponding description can be performed using non-invasive real brain data, such as EEG data, or invasive real brain data, such as SEEG data.
- non-invasive real brain data such as EEG data
- SEEG data invasive real brain data
- This epileptogenic zone can than be used as a candidate for surgical resection.
- an exemplary implementation of a method embodiment A patient was diagnosed with bi-temporal epilepsy and experienced simple and secondary generalized seizures, which were accompanied by deja-vu hallucinations, associated with palpitations, horripilation, and frisson sensations. Fixed gaze, chewing up and pallor were observed during the seizure. In the post-critic period the patient showed temporal disorientation, repetition of the same questions and retrograde amnesia during one week.
- the MRI examination revealed a hypothalamic hamartoma.
- Surface EEG recordings revealed interictal spikes and indicated a bias towards the left hemisphere. Based on the presurgical evaluation, seven SEEG electrodes were implanted in the left hemisphere, and two in the right hemisphere.
- One electrode was implanted in the hypothalamic hamartoma. During two weeks of continuous SEEG recordings, 6 simple seizures localized in the right hippocampus, and two complex seizures starting in the right hippocampus and then recruiting the left hippocampus, the left temporal lobe and the hypothalamic hamartoma were recorded.
- the large-scale connectivity of the patient was reconstructed, in particular the weight and tract length matrices generating a structural skeleton model.
- the tract length matrix divided by the signal transmission speed defines the time delays, thereby establishing the space-time structure of the cou- pling and allowing a full virtualization of the patient's brain model.
- a hypothalamic hamartoma was included in the model by changing its effect on local connectivity through factor K, j .
- each voxel of the structural skeleton was replaced by a node.
- Each node of the resulting network was set with an epileptor as described above.
- the nodes were connected via permittivity coupling, which a cti ng on a slow time scale and allowing the spread of the seizure though the network by recruiting regions not in the epileptogenic zone.
- Each node was set with a different excitability parameter xO.
- the value of the excitability was set heterogeneously across the network:
- the epileptogenic zone was set with values of excitability for initiating autonomous triggering of seizures
- the propagation zone was set with values of excitability below a value for autonomous triggering of seizures, however, the values of excitability were higher than in areas which were not considered to be part of the propagation zone.
- a systematic parameter space exploration was performed by varying the following parameters: (i) the global coupling strength G, which is a scalar factor multiplying the whole connectivity matrix, (ii) the local coupling strength G_hyp of the hypothalamus, which is a scalar factor multiplying the contribution of the hypothalamus to the connectivity matrix, (iii) the excitability values of the right hippocampus, (iv) the excitability values of the regions not recruited in the propagation zone (other regions). The excitability values of the other regions in the epileptogenic zone and the propagation zone were fixed (see Table 1).
- the brain network model was used for simulations over a period of 20 seizures and the forward solution for the SEEG electrodes was computed. Simple seizures and complex seizures were generated with similar regions recruited compared to the real SEEG recordings.
- the left hippocampus was stimulated and a propagation pattern in the left temporal lobe, similar to the SEEG recordings, was observed.
- a parameter value distribution for a parameter representing an MRI lesion can also be included in the brain network model.
- the simulation of the brain network model can be very time-intensive, methods for simplifying the calculations without loosing too much of the brain network model's dynamic behaviour will be discussed in the following.
- the coupling term KJx, - is weak compared to the other terms of the IN dimensional system, and in particular the term in x 0 .
- the difference coupling leads to weak coupling terms in a linear stability analysis.
- the weakness depends on the global coupling strength factor G and the normalization of the connectivity matrix, however one can show that this approximation holds over a larger range than for other neural network models.
- the leading eigenvectors are derived and the eigenvectors only depend on the topological connection to the epileptogenic zone.
- the linearization procedure can be easily implemented in a set of instructions for running simulations of a linearized brain network model. Since linear problems can be much more easily be computed than non-linear problems, the simulation time can be greatly reduced.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Neurology (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Neurosurgery (AREA)
- Physiology (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP16166489.1A EP3235427A1 (en) | 2016-04-21 | 2016-04-21 | Method and system for estimating a location of an epileptogenic zone of a mammalian brain |
PCT/EP2017/059532 WO2017182637A2 (en) | 2016-04-21 | 2017-04-21 | Method and system for estimating a location of an epileptogenic zone of a mammalian brain |
Publications (2)
Publication Number | Publication Date |
---|---|
EP3445239A2 true EP3445239A2 (en) | 2019-02-27 |
EP3445239B1 EP3445239B1 (en) | 2024-01-31 |
Family
ID=55970776
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP16166489.1A Withdrawn EP3235427A1 (en) | 2016-04-21 | 2016-04-21 | Method and system for estimating a location of an epileptogenic zone of a mammalian brain |
EP17720722.2A Active EP3445239B1 (en) | 2016-04-21 | 2017-04-21 | Method and system for estimating a location of an epileptogenic zone of a mammalian brain |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP16166489.1A Withdrawn EP3235427A1 (en) | 2016-04-21 | 2016-04-21 | Method and system for estimating a location of an epileptogenic zone of a mammalian brain |
Country Status (3)
Country | Link |
---|---|
US (1) | US20190223779A1 (en) |
EP (2) | EP3235427A1 (en) |
WO (1) | WO2017182637A2 (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018015778A1 (en) | 2016-07-18 | 2018-01-25 | Universite d'Aix-Marseille (AMU) | A method of modulating epileptogenicity in a patient's brain |
US20180235704A1 (en) * | 2017-02-21 | 2018-08-23 | The Cleveland Clinic Foundation | Configuring a stereo-electroencephalography procedure |
US11266355B2 (en) * | 2017-05-19 | 2022-03-08 | Cerner Innovation, Inc. | Early warning system and method for predicting patient deterioration |
EP3770918A1 (en) * | 2019-07-22 | 2021-01-27 | Universite d'Aix-Marseille (AMU) | Method for determining an onset time and an excitability of a brain region |
CN110931123B (en) * | 2019-12-06 | 2024-01-05 | 深圳市德力凯医疗设备股份有限公司 | Matrix type brain network and construction method thereof |
CN110859600A (en) * | 2019-12-06 | 2020-03-06 | 深圳市德力凯医疗设备股份有限公司 | Method for generating electroencephalogram signal, storage medium and electronic equipment |
CN111462887B (en) * | 2020-03-31 | 2023-08-29 | 首都医科大学宣武医院 | Wearable epileptic digital assistant system |
CN114926396B (en) * | 2022-04-13 | 2023-06-20 | 四川大学华西医院 | Mental disorder magnetic resonance image preliminary screening model construction method |
CN114869300B (en) * | 2022-07-08 | 2022-09-06 | 首都医科大学附属北京天坛医院 | Epileptic zone positioning device and method based on electroencephalogram, electronic device and storage medium |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8560073B2 (en) * | 2009-03-23 | 2013-10-15 | Flint Hills Scientific, Llc | System and apparatus for automated quantitative assessment, optimization and logging of the effects of a therapy |
GB0919739D0 (en) * | 2009-11-12 | 2009-12-30 | Univ Gent | Methods and systems for channel selection |
EP2950715A4 (en) * | 2013-01-31 | 2016-11-30 | Univ California | System and method for modeling brain dynamics in normal and diseased states |
-
2016
- 2016-04-21 EP EP16166489.1A patent/EP3235427A1/en not_active Withdrawn
-
2017
- 2017-04-21 US US16/094,701 patent/US20190223779A1/en active Pending
- 2017-04-21 WO PCT/EP2017/059532 patent/WO2017182637A2/en active Application Filing
- 2017-04-21 EP EP17720722.2A patent/EP3445239B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP3235427A1 (en) | 2017-10-25 |
EP3445239B1 (en) | 2024-01-31 |
WO2017182637A2 (en) | 2017-10-26 |
WO2017182637A3 (en) | 2017-12-14 |
US20190223779A1 (en) | 2019-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3445239B1 (en) | Method and system for estimating a location of an epileptogenic zone of a mammalian brain | |
Jirsa et al. | The virtual epileptic patient: individualized whole-brain models of epilepsy spread | |
US20220039736A1 (en) | Method of modulating epileptogenicity in a patient's brain | |
Samuelsson et al. | Spatial fidelity of MEG/EEG source estimates: A general evaluation approach | |
KR20200086316A (en) | Efficacy and/or therapeutic parameter recommendations using individual patient data and therapeutic brain network maps | |
US20110028827A1 (en) | Spatiotemporal pattern classification of brain states | |
JP2018508321A (en) | System and method for measuring neurovascular reactivity to brain stimulation | |
US10588561B1 (en) | Noninvasive system and method for mapping epileptic networks and surgical planning | |
EP2499585B1 (en) | Methods and systems for channel selection | |
US20150325139A1 (en) | Apparatus and method for supporting rehabilitation of brain-damaged patient | |
US20210346096A1 (en) | Methods and apparatus for electromagnetic source imaging using deep neural networks | |
Chen et al. | Functional magnetic resonance imaging for imaging neural activity in the human brain: the annual progress | |
Storti et al. | Exploring the epileptic brain network using time-variant effective connectivity and graph theory | |
US7461045B1 (en) | Optimization of spatio-temporal pattern processing for seizure warning and prediction | |
Chang et al. | Assessing recurrent interactions in cortical networks: modeling EEG response to transcranial magnetic stimulation | |
CN110392549B (en) | Systems, methods and media for determining brain stimulation that elicits a desired behavior | |
Proix et al. | Individual structural connectivity defines propagation networks in partial epilepsy | |
US20240006050A1 (en) | Locating an epileptogenic zone for surgical planning | |
Li et al. | Neural fragility of the intracranial eeg network decreases after surgical resection of the epileptogenic zone | |
Wang et al. | Virtual Epileptic Patient (VEP): Data-driven probabilistic personalized brain modeling in drug-resistant epilepsy | |
Mesejo et al. | A differential evolution-based approach for fitting a nonlinear biophysical model to fMRI BOLD data | |
Shokooh et al. | Identification of global and local states during seizures using quantitative functional connectivity and recurrence plot analysis | |
Sorrentino et al. | On the topochronic map of the human brain dynamics | |
Dunstan et al. | Global nonlinear approach for mapping parameters of neural mass models | |
Eiber et al. | A ‘Total Unique Variation Analysis’ for Brain-Machine Interfaces |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: UNKNOWN |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20181109 |
|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20210401 |
|
RAP3 | Party data changed (applicant data changed or rights of an application transferred) |
Owner name: CODEBOX GMBH |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: A61B 5/00 20060101ALI20230531BHEP Ipc: A61B 5/055 20060101AFI20230531BHEP |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
INTG | Intention to grant announced |
Effective date: 20230822 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20231219 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D Ref country code: CH Ref legal event code: EP |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602017078861 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20240219 Year of fee payment: 8 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG9D |