CN115005802A - Method, system and device for positioning onset part of brain network disease - Google Patents
Method, system and device for positioning onset part of brain network disease Download PDFInfo
- Publication number
- CN115005802A CN115005802A CN202210858565.8A CN202210858565A CN115005802A CN 115005802 A CN115005802 A CN 115005802A CN 202210858565 A CN202210858565 A CN 202210858565A CN 115005802 A CN115005802 A CN 115005802A
- Authority
- CN
- China
- Prior art keywords
- space
- brain
- pet
- data
- metabolic
- 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
- 210000004556 brain Anatomy 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 44
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 38
- 201000010099 disease Diseases 0.000 title claims abstract description 37
- 238000002600 positron emission tomography Methods 0.000 claims abstract description 69
- 230000002503 metabolic effect Effects 0.000 claims abstract description 56
- 210000004884 grey matter Anatomy 0.000 claims abstract description 17
- 238000007619 statistical method Methods 0.000 claims abstract description 14
- 230000000877 morphologic effect Effects 0.000 claims abstract description 13
- 230000001360 synchronised effect Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 4
- 238000002271 resection Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000004807 localization Effects 0.000 claims description 3
- 238000000692 Student's t-test Methods 0.000 claims description 2
- 230000003902 lesion Effects 0.000 claims description 2
- 238000012353 t test Methods 0.000 claims description 2
- 230000001717 pathogenic effect Effects 0.000 claims 1
- 206010015037 epilepsy Diseases 0.000 description 26
- 206010010904 Convulsion Diseases 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000000977 initiatory effect Effects 0.000 description 4
- 238000001356 surgical procedure Methods 0.000 description 4
- 230000001037 epileptic effect Effects 0.000 description 3
- 230000002397 epileptogenic effect Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 210000004885 white matter Anatomy 0.000 description 3
- 208000019901 Anxiety disease Diseases 0.000 description 2
- 206010003805 Autism Diseases 0.000 description 2
- 208000020706 Autistic disease Diseases 0.000 description 2
- 206010021034 Hypometabolism Diseases 0.000 description 2
- 230000036506 anxiety Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000001054 cortical effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 201000000980 schizophrenia Diseases 0.000 description 2
- 238000010176 18-FDG-positron emission tomography Methods 0.000 description 1
- 206010020674 Hypermetabolism Diseases 0.000 description 1
- 238000012274 Preoperative evaluation Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 230000007585 cortical function Effects 0.000 description 1
- 238000007428 craniotomy Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- -1 depression Diseases 0.000 description 1
- 238000000537 electroencephalography Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000001061 forehead Anatomy 0.000 description 1
- 150000002303 glucose derivatives Chemical class 0.000 description 1
- 230000004153 glucose metabolism Effects 0.000 description 1
- 230000002806 hypometabolic effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003068 molecular probe Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 210000003478 temporal lobe Anatomy 0.000 description 1
- 201000008914 temporal lobe epilepsy Diseases 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Images
Classifications
-
- 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/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/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
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/501—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
-
- 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
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Neurology (AREA)
- Biomedical Technology (AREA)
- Neurosurgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Physiology (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The embodiment of the specification provides a method, a system and a device for positioning an onset part of a brain network disease, wherein the method comprises the following steps: projecting PET metabolic data of a subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people; in the same structure space, eliminating the morphological difference among the subjects, and carrying out the same structure registration among the subjects; specifically extracting PET (positron emission tomography) metabolic data of gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain PET metabolic statistical data; determining regions of deviation from normal based on the PET metabolic statistics, and locating the regions of deviation from normal as foci.
Description
Technical Field
The present invention relates to the field of computer technologies, and in particular, to a method, a system, and an apparatus for positioning an onset of brain network disease.
Background
In the prior art, the magnetic resonance negative epilepsy, namely the epilepsy with definite brain structure change can not be found by the magnetic resonance and the conventional means, and the seizure-free rate of the epilepsy after 10 years is only 59 percent in the latest literature. At present, the latest method for realizing positioning based on the comparison of PET and a healthy person template in the world mainly depends on a Statistical Parametric Mapping (SPM) technical system, and the positioning efficiency of the system is different from 40 to 69 percent according to the positioning difficulty of the grouped cases. Currently, magnetic resonance negative epilepsy cannot be effectively localized by a single examination method on a global scale, and therefore, the effect of surgical diagnosis and treatment is severely limited.
The non-temporal lobe derived magnetic resonance negative epilepsy technology can not find structural focus due to routine examination, can not distinguish from normal cortex under microscope in craniotomy, and does not consider surgical treatment in the past. However, as technology has developed, the proportion and number of magnetic resonance negative epilepsy procedures performed has increased over the last decade. Recent clinical studies have shown that 90% of negative epileptic patients of focal origin can benefit from surgical procedures. By 10 years post-surgery, however, only 59% of patients reach a seizure-free state.
The improvement of the operation curative effect mainly depends on the improvement of the positioning accuracy. However, the traditional method is difficult and serious: firstly, noninvasive electroencephalography (EEG) collected from scalp can only find that the surface of brain is more than 6cm 2 The cooperative discharge of the brain cannot collect the brain electricity of the deep part of the brain (temporal fundus, orbital forehead, deep side of the longitudinal columns of the brain, and the like), and because the origin of the epilepsy lacks specificity, the missed diagnosis risk is also large; secondly, because the encephalic brain electricity is invasive, the coverage range is limited, the earliest detected electrical activity is probably a conducted signal rather than an initial signal, and if the noninvasive means is inaccurate, the accuracy of the invasive means cannot be further ensured; and thirdly, as the name suggests, the magnetic resonance negative epilepsy is the most accurate structural scanning method, and the magnetic resonance can not find clear focus. The above aspects illustrate that magnetic resonance negative epilepsy is challenging to locate by means in the prior art literature. At present, no single checking method can realize the full accuracy of epilepsy locationismAnd covering, and the preoperative evaluation depends on the experience of a doctor, integrates multiple examinations, discusses authenticity and judges the position of the epilepsy empirically.
Research shows that epilepsy is mainly the pathological change of cortical gray matter, and when a technical system of the traditional SPM is compared with a healthy person, the epilepsy not only contains cortical gray matter, but also contains white matter, so that the statistical accuracy is reduced. Moreover, when the SPM is used for comparing the testees, the different tested sulci brains with different walks cannot be completely projected to a uniform space, so that the accuracy of comparison among the testees is reduced.
Currently, F-FDG-PET can realize the observation of the metabolic activity of the cerebral cortex through the uptake of a glucose analog labeled radionuclide. Clinically, early localization of epilepsy is currently discussed primarily by relying on the flesh eye to look for areas of decreased glucose metabolism during the interval between seizures. However, it has limitations in that: there are physiological hypometabolism areas in the brain, i.e., hypometabolism areas that exist in healthy people; ② in patients, hypo-metabolic regions may also be present in the brain areas to which epilepsy spreads; ③ hypermetabolism may also occur in the focus of epilepsy, because this part of the brain is higher in the healthy people, even if slightly reduced, it may be higher than other brain areas in the brain. Therefore, the accurate judgment of the single mode of the epileptogenic focus cannot be realized by reading the picture clinically by naked eyes.
In the image post-processing, based on the method of SPM, the first-line method is to transform the brain PET data of the patient and the healthy brain into a unified template for comparison through deformation in the original state of a 3D matrix.
As mentioned above, the latest first-line technology in the prior art is mainly the post-processing technology under the SPM system. Because different individuals have different sulci gyrus depths and bending degrees, the method system can only unify the widths, the lengths and the like of different brains, can not ensure the regular alignment of the same gyrus of different individuals on spatial positions, has clear registration errors generally, and obstructs the observation of tiny focuses.
Epilepsy is a gray-colored pathology. In the traditional method, adjacent cerebrospinal fluid and white matter structures are simultaneously brought into the statistics, so that the number of mixed factors brought into the data is increased when the data is compared and analyzed. For various reasons, the accuracy of the traditional first-line method for epileptogenic focus is only 40% -69%. Among them, the accuracy rate of the simple temporal lobe epilepsy is about 69%, and the accuracy rate of the magnetic resonance negative epilepsy (the epilepsy with the highest positioning difficulty in the report) in which the cortical function silent region originates and the PET original data can not be sufficiently positioned is only about 40%.
Disclosure of Invention
The invention aims to provide a method, a system and a device for positioning an attack initiation part of a brain network disease, and aims to solve the problems in the prior art.
The invention provides a method for positioning an onset part of brain network diseases, which comprises the following steps:
projecting PET metabolic data of a subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people;
in the same structure space, eliminating the morphological difference among the subjects, and carrying out the same-original structure registration among the subjects;
specifically extracting PET (positron emission tomography) metabolic data of gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain PET metabolic statistical data;
determining a region deviating from normal based on the PET metabolic statistics, and locating the region deviating from normal as a focus.
The invention provides a system for positioning the onset part of brain network diseases, which comprises:
a projection module for projecting the PET metabolic data of a subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people;
the registration module is used for eliminating the morphological difference among the subjects in the same structural space and carrying out the same-original structural registration among the subjects;
the extraction module is used for specifically extracting the PET metabolic data of the gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain the PET metabolic statistical data;
and the positioning module is used for determining a region deviating from normal based on the PET metabolic statistical data and positioning the region deviating from normal as a focus.
The embodiment of the present invention further provides a device for positioning an onset part of brain network diseases, including: the computer program is stored on the memory and can run on the processor, and when being executed by the processor, the computer program realizes the steps of the method for positioning the onset part of the brain network disease.
The embodiment of the invention also provides a computer readable storage medium, wherein an implementation program for information transmission is stored on the computer readable storage medium, and the program is executed by a processor to implement the steps of the method for positioning the onset part of the brain network disease.
By adopting the embodiment of the invention, the accurate projection of the brain metabolic data of a single patient and healthy people in a unified space is realized, the projection difference among the subjects is fully eliminated, meanwhile, the gray matter range is defined based on 3-dimensional reconstruction and segmentation, the PET metabolic signal in the middle of the gray matter is accurately extracted, the comparison among samples is carried out, and the positioning accuracy is improved.
Drawings
In order to more clearly illustrate one or more embodiments or prior art solutions of the present specification, the drawings that are needed in the description of the embodiments or prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the present specification, and that other drawings can be obtained by those skilled in the art without inventive exercise.
FIG. 1 is a flowchart of a method for locating an onset of brain network disease according to an embodiment of the present invention;
FIG. 2 is a graphical representation of PET values at 3D brain surface vertices for an embodiment of the present invention;
FIG. 3 is a schematic diagram of the effect of locating actual problem cases according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a system for locating the onset of brain network disease according to an embodiment of the present invention;
fig. 5 is a schematic diagram of a device for locating an onset site of brain network disease according to an embodiment of the present invention.
Detailed Description
In order to make those skilled in the art better understand the technical solutions in one or more embodiments of the present disclosure, the technical solutions in one or more embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in one or more embodiments of the present disclosure, and it is obvious that the described embodiments are only a part of the embodiments of the present disclosure, and not all embodiments. All other embodiments that can be derived by a person skilled in the art from one or more of the embodiments described herein without making any inventive step shall fall within the scope of protection of this document.
Method embodiment
According to an embodiment of the present invention, a method for positioning an onset start site of a brain network disease is provided, fig. 1 is a flowchart of the method for positioning the onset start site of the brain network disease according to the embodiment of the present invention, as shown in fig. 1, the method for positioning the onset start site of the brain network disease according to the embodiment of the present invention specifically includes:
step 101, projecting PET metabolic data of a subject to the same structural space, wherein the subject specifically comprises: brain network disease patients and healthy people; wherein, in the embodiment of the invention, the brain network diseases are: magnetic resonance negative epilepsy, and may also include other brain network disorders, such as parkinson, depression, schizophrenia, autism, anxiety, and the like.
Step 101 specifically includes:
acquiring PET metabolic data of a subject and 3D MR T1 weighted data, namely structural space data, of a magnetic resonance negative epileptic patient, registering the PET metabolic data into a T12D space, and performing brain structure reconstruction and grey-white substance segmentation of a T13D space through a T13D reconstruction algorithm of FreeSenfer technical standard based on the 3D MR T1 weighted data, wherein the same structural space comprises a T12D space and a T13D space.
102, eliminating morphological differences among the subjects in the same structural space, and carrying out original structure registration among the subjects; step 102 specifically includes:
based on the brain structure information in the T12D space and the T13D space, morphological differences of the inter-subject brain structures are eliminated through a synchronous volume and surface CVSR algorithm, and the homomorphic structure registration of the inter-subject brain structures is carried out.
103, specifically extracting PET (positron emission tomography) metabolic data of gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain PET metabolic statistical data; step 103 specifically comprises:
and interpolating the PET metabolic data into the registration coordinates of the corresponding brain structures based on the registration coordinates among the subjects of the brain structures in the same space, and performing peak-to-peak statistics on the PET metabolic data values of the same space structure at the level of the peak of the 3D model surface by adopting a preset statistical method to obtain PET metabolic statistical data. In an embodiment of the present invention, the predetermined statistical method specifically includes one of the following: z-transformation and t inspection;
and 104, determining a deviation area from normal based on the PET metabolic statistical data, and positioning the deviation area from normal as a focus. Step 104 specifically includes: based on the PET metabolism statistical data, determining a deviation area from normal according to a position corresponding to a lowest z value of the whole brain or a position corresponding to a value with statistical difference in a T distribution diagram of a T13D space, and positioning the deviation area from normal as a focus.
Preferably, in the embodiment of the present invention, the brain space after locating the focus can be used as a standard space of a population level for scientific research or observation of rules and attributes of population level focus distribution; in addition, the brain space after the focus is positioned can be projected to the individual space of the patient through the inverse algorithm of the synchronous volume and surface CVSR algorithm, and is guided into a surgical navigation and robot system for guiding surgical resection.
In conclusion, according to the technical scheme of the embodiment of the invention, the position matching between the brain of the patient and the brain homomorphism of the healthy person is realized by eliminating the brain morphological difference among the subjects, so that a basis is provided for further quantitative statistics; the embodiment of the invention provides a method for accurately extracting the metabolic value in grey matter to count the degree of deviation of a patient from a normal population; according to the embodiment of the invention, the most obvious change part is accurately found through the quantification method, and the surgical treatment is guided.
The following describes the above technical solution of the embodiment of the present invention in detail by taking magnetic resonance negative epilepsy as an example with reference to the accompanying drawings.
According to the embodiment of the invention, PET metabolic data of a patient and PET metabolic data of healthy people are projected to the same space, on the premise of fully eliminating morphological differences among subjects, gray matter signals are specifically extracted, statistical comparison is carried out on epilepsy which is a gray matter disease, and a region deviating from normal is found, namely an epileptogenic focus. The method comprises the following concrete steps:
step 1, pretreatment.
Step 101, space division:
the embodiment of the invention requires that the patient has 3D MR T1 weighted data at the same time, and the weighted data can be acquired by synchronous PET/MR data acquisition and can also be acquired by PET/CT and common magnetic resonance acquisition respectively; registering the PET into T12D space;
based on the 3D MR T1 data, brain structure reconstruction and grey-white substance segmentation were performed by the FreSurfer's Standard T13D reconstruction algorithm.
Step 102, registration among subjects, unifying structural spaces of patients and healthy people:
inter-subject Registration is further performed by a Combined Volumetric and Surface Registration algorithm. The algorithm simultaneously refers to grey white matter information in a 2D space and structure information reconstructed in a 3D space, so that the matching degree of the same structure alignment among the testees can reach more than 95 percent;
step 103, extracting gray matter median PET metabolic information:
as shown in fig. 2, the PET signal is interpolated into corresponding structure coordinates based on the inter-subject registration coordinates of the structure, and vertex-to-vertex statistics is performed on the PET signal values of the same spatial structure at the level of the vertices of the 3D model plane (corresponding to the pixels of the 2D plane). Statistical methods include, but are not limited to: z-transform (statistical standard deviation of patients from normal human-mean), t-test of single patient against healthy population, etc.;
and 2, performing positioning analysis.
Specifically, seizure-inducing property is defined by the lowest z-value of the whole brain, or the position corresponding to a statistically different value in the 3D t distribution map (mathematically meaning the most significant position from the normal population). The positioning effect is shown in fig. 3.
The brain space after focus positioning can continue to be a standard space of a crowd level and is used for scientific research and observation of the rule and attribute of the crowd level seizure focus distribution;
or projected to the individual space of the patient through the inverse algorithm in step 102, and introduced into the surgical navigation and robotic system to guide the surgical resection.
In the embodiment of the invention, based on grouped 42 magnetic resonance negative epilepsy originated from the cortex of the functional silent zone (the case with the highest positioning difficulty all over the world), the SPM method system (first-line method) in the prior art is adopted, so that the epilepsy positioning rate can be only 47%, and the positioning rate can reach 87% by the technical scheme of the embodiment of the invention, so that the problem of magnetic resonance negative epilepsy operation evaluation can be fully solved.
In conclusion, by means of the technical scheme of the embodiment of the invention, the problem of accurate alignment of homologous brain regions of different subjects is fully solved, so that the unified structures of patients and healthy people can be perfectly aligned and overlapped, and the uniformity of space during statistics is ensured; meanwhile, only signals of the grey matter structure are extracted for statistics, and the attributes of the epileptic grey matter disease are better met. By adopting the technical scheme of the embodiment of the invention, even if the epilepsy is in a magnetic resonance negative area, a PET non-specific area and a functional silent area, the positioning accuracy rate reaches 87%. The positioning accuracy is obviously improved.
It should be noted that based on the technical scheme of the embodiment of the present invention, a norm can be added according to different PET molecular probes to find therapeutic targets for other brain network diseases such as parkinson, depression, schizophrenia, autism, anxiety and the like.
System embodiment
According to an embodiment of the present invention, a system for locating a brain network disease attack initiation site is provided, fig. 4 is a schematic view of the system for locating a brain network disease attack initiation site according to an embodiment of the present invention, as shown in fig. 4, the system for locating a brain network disease attack initiation site according to an embodiment of the present invention specifically includes:
a projection module 40 for projecting the PET metabolic data of the subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people; the projection module 40 is specifically configured to:
acquiring PET metabolic data of a subject and 3D MR T1 weighted data, namely structural space data, of a brain network disease patient, registering the PET metabolic data into a T12D space, and performing brain structural reconstruction and grey-white substance segmentation of a T13D space through a T13D reconstruction algorithm of FreeSenfer technical standard based on the 3D MR T1 weighted data, wherein the same structural space comprises a T12D space and a T13D space.
A registration module 42, configured to eliminate morphological differences between subjects in the same structural space, and perform orthotopic structural registration between subjects; the registration module 42 is specifically configured to:
based on the brain structure information in the T12D space and the T13D space, morphological differences of the inter-subject brain structures are eliminated through a synchronous volume and surface CVSR algorithm, and the homomorphic structure registration of the inter-subject brain structures is carried out.
An extraction module 44, configured to specifically extract the gray PET metabolic data based on the registration coordinates between the subjects in the same structural space, and perform statistical analysis to obtain PET metabolic statistical data; the extraction module 44 is specifically configured to: and interpolating the PET metabolic data into the registration coordinates of the corresponding brain structures based on the registration coordinates among the subjects of the brain structures in the same space, and performing peak-to-peak statistics on the PET metabolic data values of the same space structure at the level of the peak of the 3D model surface by adopting a preset statistical method to obtain PET metabolic statistical data. The predetermined statistical method specifically comprises one of: z-transformation and t inspection;
a localization module 46 for determining regions deviating from normal based on the PET metabolic statistics, and localizing the regions deviating from normal as foci. The positioning module 46 is specifically configured to: based on the PET metabolic statistics, determining an off-normal region according to a position corresponding to a lowest z value of the whole brain or a position corresponding to a statistically different value in a T distribution diagram of a T13D space, and positioning the off-normal region as a lesion.
In the embodiment of the present invention, the method further includes:
the first processing module is used for taking the brain space after the focus is caused to be positioned as a standard space of a crowd level and is used for scientific research or observing the rule and the attribute of the distribution of the focus caused to be caused to the crowd level;
and the second processing module is used for projecting the brain space after the focus is positioned to the individual space of the patient through the inverse algorithm of the synchronous volume and surface CVSR algorithm, and guiding the brain space to an operation navigation and robot system for guiding the surgical resection.
The embodiment of the present invention is a system embodiment corresponding to the above method embodiment, and specific operations of each module may be understood with reference to the description of the method embodiment, which is not described herein again.
Apparatus embodiment one
The embodiment of the present invention provides a device for positioning an onset part of brain network diseases, as shown in fig. 5, including: a memory 50, a processor 52 and a computer program stored on the memory 50 and executable on the processor 52, which computer program, when executed by the processor 52, performs the steps as described in the method embodiments.
Device embodiment II
An embodiment of the present invention provides a computer-readable storage medium, on which an implementation program for information transmission is stored, and when the program is executed by the processor 52, the steps described in the method embodiment are implemented.
The computer-readable storage medium of this embodiment includes, but is not limited to: ROM, RAM, magnetic or optical disks, and the like.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for locating the onset of brain network disease, comprising:
projecting PET metabolic data of a subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people;
in the same structure space, eliminating the morphological difference among the subjects, and carrying out the same-original structure registration among the subjects;
specifically extracting PET (positron emission tomography) metabolic data of gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain PET metabolic statistical data;
determining regions of deviation from normal based on the PET metabolic statistics, and locating the regions of deviation from normal as foci.
2. The method of claim 1, wherein projecting the subject's PET metabolic data to the same structural space specifically comprises:
acquiring PET metabolic data of a subject and 3D MR T1 weighted data, namely structural space data, of a brain network disease patient, registering the PET metabolic data into a T12D space, and performing brain structure reconstruction and grey white segmentation of a T13D space through a T13D reconstruction algorithm of FreeKufer technical standard based on the 3D MR T1 weighted data, wherein the same structural space comprises a T12D space and a T13D space.
3. The method according to claim 2, wherein inter-subject morphological differences are eliminated in the same structural space, and performing inter-subject orthotopic structural registration specifically comprises:
based on the brain structure information in the T12D space and the T13D space, morphological differences of the inter-subject brain structures are eliminated through a synchronous volume and surface CVSR algorithm, and the homomorphic structure registration of the inter-subject brain structures is carried out.
4. The method according to claim 1, wherein PET metabolic data of gray matter are specifically extracted and statistically analyzed based on the inter-subject registration coordinates of the same structural space, and the obtaining of PET metabolic statistical data specifically comprises:
and interpolating the PET metabolic data into the registration coordinates of the corresponding brain structures based on the registration coordinates among the subjects of the brain structures in the same space, and performing peak-to-peak statistics on the PET metabolic data values of the same space structure at the level of the peak of the 3D model surface by adopting a preset statistical method to obtain PET metabolic statistical data.
5. The method according to claim 4, wherein the predetermined statistical method comprises in particular one of: z-transformation and t-test.
6. The method of claim 1, wherein determining regions of departure from normality based on the PET metabolic statistics, the locating regions of departure from normality as a pathogenic cooktop body comprises:
based on the PET metabolic statistics, determining an off-normal region according to a position corresponding to a lowest z value of the whole brain or a position corresponding to a statistically different value in a T distribution diagram of a T13D space, and positioning the off-normal region as a lesion.
7. The method of claim 1, further comprising:
taking the brain space after the focus is positioned as a standard space of the crowd level, and using the brain space for scientific research or observing the rule and the attribute of the distribution of the focus of the crowd level;
the brain space after the focus is positioned is projected to the individual space of the patient through the inverse algorithm of the synchronous volume and surface CVSR algorithm, and is guided into a surgical navigation and robot system for guiding surgical resection.
8. A system for locating the onset of brain network disease, comprising:
a projection module for projecting the PET metabolic data of a subject to the same structural space, wherein the subject comprises in particular: brain network disease patients and healthy people;
the registration module is used for eliminating the morphological difference among the subjects in the same structural space and carrying out the same-original structural registration among the subjects;
the extraction module is used for specifically extracting the PET metabolic data of the gray matter based on the registration coordinates between the subjects in the same structural space, and performing statistical analysis to obtain the PET metabolic statistical data;
and the positioning module is used for determining a region deviating from the normal based on the PET metabolic statistical data and positioning the region deviating from the normal as a focus.
9. A device for locating the onset of brain network disease, comprising: memory, processor and computer program stored on the memory and executable on the processor, which computer program, when executed by the processor, carries out the steps of the method of brain network seizure onset localization according to any one of claims 1 to 7.
10. A computer-readable storage medium, wherein a program for implementing information transmission is stored on the computer-readable storage medium, and when the program is executed by a processor, the program implements the steps of the method for locating the onset of brain network disease onset according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210858565.8A CN115005802B (en) | 2022-07-21 | 2022-07-21 | Method, system and device for positioning initial part of brain network disease onset |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210858565.8A CN115005802B (en) | 2022-07-21 | 2022-07-21 | Method, system and device for positioning initial part of brain network disease onset |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115005802A true CN115005802A (en) | 2022-09-06 |
CN115005802B CN115005802B (en) | 2023-04-28 |
Family
ID=83082478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210858565.8A Active CN115005802B (en) | 2022-07-21 | 2022-07-21 | Method, system and device for positioning initial part of brain network disease onset |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115005802B (en) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100004527A1 (en) * | 2008-07-01 | 2010-01-07 | Dale Anders M | Identifying White Matter Fiber Tracts Using Magnetic Resonance Imaging (MRI) |
US20130223714A1 (en) * | 2012-02-28 | 2013-08-29 | Albert Einstein College Of Medicine Of Yeshiva University | Methods for quantitative assessment of volumetric image from a subject and uses therof |
US20140270424A1 (en) * | 2013-03-15 | 2014-09-18 | Mim Software Inc. | Population-guided deformable registration |
CN107374705A (en) * | 2017-06-20 | 2017-11-24 | 上海交通大学 | A kind of lung puncture location of operation method under X-ray and preoperative CT guiding |
WO2018098213A1 (en) * | 2016-11-23 | 2018-05-31 | Wake Forest University Health Sciences | Medical image analysis using mechanical deformation information |
WO2018106713A1 (en) * | 2016-12-06 | 2018-06-14 | Darmiyan, Inc. | Methods and systems for identifying brain disorders |
CN109952597A (en) * | 2016-11-16 | 2019-06-28 | 索尼公司 | Brain registration between patient |
CN110013263A (en) * | 2019-04-04 | 2019-07-16 | 首都医科大学宣武医院 | One kind generating the method and system of standardized uptake value (SUV) based on medical image data |
CN112927187A (en) * | 2021-01-27 | 2021-06-08 | 张凯 | Method for automatically identifying and positioning focal cortical dysplasia epileptic focus |
CN113112476A (en) * | 2021-04-14 | 2021-07-13 | 中国人民解放军北部战区总医院 | Method and system for identifying epileptogenic focus of temporal lobe epilepsy caused by hippocampus sclerosis and/or predicting pathological typing of temporal lobe epilepsy |
WO2022020391A2 (en) * | 2020-07-20 | 2022-01-27 | The Trustees Of Columbia University In The City Of New York | Precision targeted retromer therapeutics for the treatment of neurodegenerative diseases and disorders |
CN114334130A (en) * | 2021-12-25 | 2022-04-12 | 浙江大学 | PET molecular image computer-aided diagnosis system based on brain symmetry |
-
2022
- 2022-07-21 CN CN202210858565.8A patent/CN115005802B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100004527A1 (en) * | 2008-07-01 | 2010-01-07 | Dale Anders M | Identifying White Matter Fiber Tracts Using Magnetic Resonance Imaging (MRI) |
US20130223714A1 (en) * | 2012-02-28 | 2013-08-29 | Albert Einstein College Of Medicine Of Yeshiva University | Methods for quantitative assessment of volumetric image from a subject and uses therof |
US20140270424A1 (en) * | 2013-03-15 | 2014-09-18 | Mim Software Inc. | Population-guided deformable registration |
CN109952597A (en) * | 2016-11-16 | 2019-06-28 | 索尼公司 | Brain registration between patient |
WO2018098213A1 (en) * | 2016-11-23 | 2018-05-31 | Wake Forest University Health Sciences | Medical image analysis using mechanical deformation information |
WO2018106713A1 (en) * | 2016-12-06 | 2018-06-14 | Darmiyan, Inc. | Methods and systems for identifying brain disorders |
CN107374705A (en) * | 2017-06-20 | 2017-11-24 | 上海交通大学 | A kind of lung puncture location of operation method under X-ray and preoperative CT guiding |
CN110013263A (en) * | 2019-04-04 | 2019-07-16 | 首都医科大学宣武医院 | One kind generating the method and system of standardized uptake value (SUV) based on medical image data |
WO2022020391A2 (en) * | 2020-07-20 | 2022-01-27 | The Trustees Of Columbia University In The City Of New York | Precision targeted retromer therapeutics for the treatment of neurodegenerative diseases and disorders |
CN112927187A (en) * | 2021-01-27 | 2021-06-08 | 张凯 | Method for automatically identifying and positioning focal cortical dysplasia epileptic focus |
CN113112476A (en) * | 2021-04-14 | 2021-07-13 | 中国人民解放军北部战区总医院 | Method and system for identifying epileptogenic focus of temporal lobe epilepsy caused by hippocampus sclerosis and/or predicting pathological typing of temporal lobe epilepsy |
CN114334130A (en) * | 2021-12-25 | 2022-04-12 | 浙江大学 | PET molecular image computer-aided diagnosis system based on brain symmetry |
Non-Patent Citations (8)
Title |
---|
ANAND A. JOSHI等: "Surface-Constrained Volumetric Brain Registration Using Harmonic Mappings" * |
KUN GUO等: "Assessment of localization accuracy and postsurgical prediction of simultaneous 18F-FDG PET/MRI in refractory epilepsy patients", 《EUROPEAN RADIOLOGY》 * |
严序等: "基于大脑模板的MRI扫描自动定位方法" * |
吴小伟: "颞叶癫痫患者脑灰质与白质的基于体素的形态测量学研究", 《中国优秀硕士论文全文数据库医药卫生科技辑》 * |
孙强: "一体化PET/MR在难治性癫痫致痫灶诊断中的价值", 《中国优秀硕士论文全文数据库医药卫生科技辑》 * |
罗述谦等: "医学图像配准技术" * |
赵春雷: "基于统计参数图方法对颞叶癫痫PET定侧诊断价值的研究", 《中国博士论文全文数据库医药卫生科技辑》 * |
郭六生: "基于特征点的图像配准技术研究", 《中国优秀硕士论文全文数据库信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115005802B (en) | 2023-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109965895B (en) | Method for constructing epileptic focus positioning classifier based on brain image fusion characteristics | |
Stufflebeam et al. | Localization of focal epileptic discharges using functional connectivity magnetic resonance imaging | |
CN110689536B (en) | Brain grey matter and white matter tracking method and device based on multi-mode magnetic resonance image | |
Van Mierlo et al. | Ictal EEG source localization in focal epilepsy: Review and future perspectives | |
US11205267B2 (en) | Method for localizing implanted intracranial electrode | |
WO2006060727A2 (en) | System and method for eeg imaging of cerebral activity using electrode sets | |
WO2023280086A1 (en) | Target determination method and apparatus, electronic device, storage medium, and neuromodulation device | |
CN105708462A (en) | Data processing method based on rfMRI (resting-state functional magnetic resonance imaging) of idiopathic epilepsy | |
Baroumand et al. | Automated ictal EEG source imaging: A retrospective, blinded clinical validation study | |
CN113112476B (en) | Method and system for identifying epileptogenic focus and/or predicting pathological typing of epileptogenic focus | |
CN110537915A (en) | Corticospinal tract fiber tracking method based on FMRI and DTI fusion | |
CN115005802A (en) | Method, system and device for positioning onset part of brain network disease | |
Liao et al. | Tracing the deformed midline on brain CT | |
Buch et al. | Evaluating the accuracy of perfusion/metabolism (SPET/PET) ratio in seizure localization | |
van de Velden et al. | Effects of inverse methods and spike phases on interictal high-density EEG source reconstruction | |
US11972564B2 (en) | Recording medium, information processing device, information processing method, trained model generation method, and correlation image output device | |
Qu et al. | Interictal electrophysiological source imaging based on realistic epilepsy head model in presurgical evaluation: A prospective study | |
RU2691306C1 (en) | Method for assessing the involvement of disturbed liquor dynamics in the development of a diffuse white matter disorder in cerebral microangiopathy | |
Jeong et al. | Localization of function‐specific segments of the primary motor pathway in children with Sturge‐Weber syndrome: a multimodal imaging analysis | |
CN115312181B (en) | Prognosis judging method, system and device for brain network disease operation | |
CN110192859B (en) | Method for measuring human brain lymph-like pathway clearing efficiency by magnetic resonance | |
Ji et al. | Diffusion tensor imaging combined with nerve fiber bundle tracing in acute cerebral infarction | |
Moreno | A Retrospective Chart Review on Young Adult Temporal Lobe Epileptic Patients with and without Mesial Temporal Sclerosis Using Radiomic Feature Analysis | |
Ivanovska et al. | Epidemiological Image Analysis: Insights, Perspectives, and Challenges | |
Leiberg et al. | Effects of anterior temporal lobe resection on cortical morphology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |