CN106845137A - A kind of sacred disease analysis method based on brain network structure - Google Patents
A kind of sacred disease analysis method based on brain network structure Download PDFInfo
- Publication number
- CN106845137A CN106845137A CN201710107488.1A CN201710107488A CN106845137A CN 106845137 A CN106845137 A CN 106845137A CN 201710107488 A CN201710107488 A CN 201710107488A CN 106845137 A CN106845137 A CN 106845137A
- Authority
- CN
- China
- Prior art keywords
- network
- node
- brain
- shortest path
- metric
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
A kind of sacred disease analysis method based on brain network structure, comprises the following steps:Diffusion tensor imaging data are pre-processed by step one;Step 2, node and side are two basic elements of brain network;Definition on encephalomere point, 246 cortex interested and infracortical region, a node of each Regional Representative's network are divided into according to human brain group of networks collection of illustrative plates by whole brain;For the side of network, using the fiber tracking in FSL instruments, and constructed by sparse threshold method and have no right two-value network;Step 3, calculates the whole world and Local Area Network measurement, and calculate the TG-AUC AUC of each network metric;AUC to each network metric and each network metric performs nonparametric displacement test determination in network attribute with the presence or absence of significant group difference, while by existing experimental data, coming judgement sample disease condition and information.Auxiliary diagnosis and the analysis for effectively realizing sacred disease of the invention.
Description
Technical field
It is especially a kind of to be based on brain net the present invention relates to the medical imaging under computer graphics, Nervous System Anatomy field
The sacred disease analysis method of network structure.
Background technology
With the development in epoch, the progress of Medical Imaging Technology, diffusion tensor imaging is in the research of Neuscience
Increasing influence power is accounted for, it is that this epoch is indispensable to possess advanced neuroimaging technology;Diffusion tensor imaging skill
Art is a kind of method of emerging description brain structure;At present, diffusion tensor imaging is just widely used in psychiatric department
In disease and the supplementary means of diagnosis, it might even be possible to for the formulation of pre-operative surgical scheme, it may be said that its tribute in medical domain
Offer the advantage that can not be substituted;So having great meaning for brain science to the research based on diffusion tensor imaging
Justice.
Now, the brain network built with diffusion tensor imaging data has many important topological properties, such as " small generation
Boundary " attribute, modular institutional framework;On the other hand, researcher has found that many mental illnesses (divide by such as stages alzheimer's disease and spirit
Split disease etc.) it is relevant with the change in topology of brain structure, brain function Network Abnormal, these researchs are not only to understand neuropsychiatric disease
Pathomechanism provides New Century Planned Textbook, can also be evaluated for the early prevention of disease and treatment and provide brain network video mark.
The content of the invention
In order to overcome the shortcomings of that the auxiliary diagnosis of sacred disease cannot be realized in the prior art and analyze, the present invention is combined and expanded
Tensor imaging technology is dissipated, the analysis based on brain network structure of a kind of auxiliary diagnosis for effectively realizing sacred disease and analysis is proposed
Method.
In order to solve the above-mentioned technical problem following technical scheme is provided:
A kind of sacred disease analysis method based on brain network structure, comprises the following steps:
Diffusion tensor imaging data are pre-processed by step one
1.1), the head in DTI data is corrected to the rigid body translation of b0 images by using each diffusion weighted images
Motion and vortex distortion;
1.2), brain area is extracted using brain extracting tool in FSL;
1.3), it is registered in often with brigadier's normed space data and each participant's b0 image using the linear image of FMRIB
Individual participant T1 data;
Step 2, sets up brain network
Node and side are two basic elements of brain network;Definition on encephalomere point, will according to human brain group of networks collection of illustrative plates
Whole brain is divided into 246 cortex interested and infracortical region, a node of each Regional Representative's network;For net
The side of network, using the fiber tracking in FSL instruments, and is constructed by sparse threshold method and haves no right two-value network;
Step 3, calculates parameters, carries out data analysis
According to the definition of brain network, the whole world and Local Area Network measurement, including cluster coefficient are calculatedShortest path length L,
Worldlet degree σ, component efficiency Eloc, global efficiency Eglob, centrad Nbc, worldlet degree σ, and calculate each network metric
TG-AUC AUC;AUC to each network metric and each network metric performs nonparametric displacement test determination in network
Whether there is significant group difference in attribute, while by existing experimental data, coming judgement sample disease condition and information.
Further, in the step 3, various brain network parameters are calculated, including following parameter:
3.1 cluster coefficients:What cluster coefficient was weighed is the grouping of the world economy degree of network, is the important parameter for measuring network:
CiRepresent the cluster coefficient of node i, eiThe number on the side actually connected between representing the nodes neighbors,Table
Show possible maximum connection side number;
Network node quantity is represented,Represent whole network;
3.2 component efficiencies:
E (i) represents the component efficiency of arbitrary node i, GiRefer to the subgraph that the neighbours of node i are constituted,Table constitutes subgraph
Quantity, i, j are arbitrary node, l in networkjkRepresent node j, the shortest path length between k;
3.3 shortest path lengths:The information that shortest path features a certain node in network reaches the optimal of another node
Path, by shortest path can quickly transmission information, so as to save system resource;
3.4 global efficiencies:Usual shortest path length will carry out computing in some connected graph, because if in network
There is disconnected node can cause the shortest path length value between the two nodes to be infinite, therefore propose:
3.5 centrads:Centrad is a statistical indicator for portraying nodes effect and status;
NbcI () is the centrad of any point i in network, G represents whole network, and i, j, k are arbitrary nodes in network,
δjkIt is the quantity of all shortest paths from node j to node k, δjkI () is the quantity by node i in these shortest paths;
3.6 worldlet degree, by worldlet degree, we can weigh the power of " worldlet " attribute, and δ, λ are calculated first;δ
The ratio of real network and random network cluster coefficients is represented, λ represents the ratio of real network and random network shortest path length
Value;
CrealRepresent real network cluster coefficient, CrandomThe cluster coefficient of table random network.
LrealRepresent real network cluster coefficient, LrandomThe cluster coefficient of table random network.
Worldlet degree σ is defined as follows:
Further, in the step 3, nonparametric displacement test mode:The average value of each network metric is calculated first
Between group difference;For each network metric, randomly all values are re-assigned in two groups, and recalculate two
Mean difference between randomization group;The randomization program repeat 10,000 time, and each distribution preceding 5% be used as have probability
Error is the critical value of the one tailed test of 0.05 null hypothesis.
Beneficial effects of the present invention are:Effectively realize auxiliary diagnosis and the analysis of sacred disease.
Specific implementation process
The invention will be further described below.
A kind of sacred disease analysis method based on brain network structure, it is characterised in that:Comprise the following steps:
Diffusion tensor imaging data are pre-processed by step one:
1.1), the head in DTI data is corrected to the rigid body translation of b0 images by using each diffusion weighted images
Motion and vortex distortion;
1.2), brain area is extracted using brain extracting tool in FSL;
1.3), it is registered in often with brigadier's normed space data and each participant's b0 image using the linear image of FMRIB
Individual participant T1 data.
Step 2, sets up brain network:
Node and side are two basic elements of brain network;Definition on encephalomere point, will according to human brain group of networks collection of illustrative plates
Whole brain is divided into 246 cortex interested and infracortical region, a node of each Regional Representative's network;For net
The side of network, using the fiber tracking in FSL instruments, and constructs two-value network by sparse threshold method;
Step 3, calculates parameters, carries out data analysis
According to the definition of brain network, the whole world and Local Area Network measurement, including cluster coefficient are calculatedShortest path length L,
Worldlet degree σ, component efficiency Eloc, global efficiency Eglob, centrad Nbc, worldlet degree σ, and calculate each network metric
TG-AUC AUC;AUC to each network metric and each network metric performs nonparametric displacement test determination in network
Whether there is significant group difference in attribute, while by existing experimental data, coming judgement sample disease condition and information.
Body calculating parameter is as follows:
3.1 cluster coefficients:What cluster coefficient was weighed is the grouping of the world economy degree of network, is the important parameter for measuring network:
CiRepresent the cluster coefficient of node i, eiThe number on the side actually connected between representing the nodes neighbors,Table
Show possible maximum connection side number;
Network node quantity is represented,Represent whole network;
3.2 component efficiencies:
E (i) represents the component efficiency of arbitrary node i, GiRefer to the subgraph that the neighbours of node i are constituted,Table constitutes subgraph
Quantity, i, j are arbitrary node, l in networkjkNode j is represented, (side number minimum leads to the shortest path length between k
Road);
3.3 shortest path lengths:The information that shortest path features a certain node in network reaches the optimal of another node
Path, by shortest path can quickly transmission information, so as to save system resource;
3.4 global efficiencies:Usual shortest path length will carry out computing in some connected graph, because if in network
There is disconnected node can cause the shortest path length value between the two nodes to be infinite, therefore propose:
3.5 centrads:Centrad is a statistical indicator for portraying nodes effect and status;
NbcI () is the centrad of any point i in network, G represents whole network, and i, j, k are arbitrary nodes in network,
δjkIt is the quantity of all shortest paths from node j to node k, δjkI () is the quantity by node i in these shortest paths;
3.6 worldlet degree, by worldlet degree, we can weigh the power of " worldlet " attribute, and δ, λ are calculated first;δ
The ratio of real network and random network cluster coefficients is represented, λ represents the ratio of real network and random network shortest path length
Value;
CrealRepresent real network cluster coefficient, CrandomThe cluster coefficient of table random network.
LrealRepresent real network cluster coefficient, LrandomThe cluster coefficient of table random network.
Worldlet degree σ is defined as follows:
Data are analyzed:Calculate the group difference between the average value of each network metric;In order to test what is observed
, for each network metric, randomly be re-assigned to all values in two groups by the possible occurrent null hypothesis of group difference,
And recalculate the mean difference between two randomization groups;The randomization program is repeated 10,000 time, and before each distribution
5% critical value for being used as the one tailed test with the null hypothesis that probable error is 0.05.
Claims (3)
1. a kind of sacred disease analysis method based on brain network structure, it is characterised in that:Comprise the following steps:
Diffusion tensor imaging data are pre-processed by step one
1.1), the head movement in DTI data is corrected to the rigid body translation of b0 images by using each diffusion weighted images
Distorted with vortex;
1.2), brain area is extracted using brain extracting tool in FSL;
1.3), each is registered in brigadier's normed space data and each participant's b0 image using the linear image of FMRIB to join
With person's T1 data;
Step 2, sets up brain network
Node and side are two basic elements of brain network;Definition on encephalomere point, will be whole according to human brain group of networks collection of illustrative plates
Brain is divided into 246 cortex interested and infracortical region, a node of each Regional Representative's network;For network
Side, using the fiber tracking in FSL instruments, and is constructed by sparse threshold method and haves no right two-value network;
Step 3, calculates parameters, carries out data analysis
According to the definition of brain network, the whole world and Local Area Network measurement, including cluster coefficient C, shortest path length L, small generation are calculated
Boundary degree σ, component efficiency Eloc, global efficiency Eglob, centrad Nbc, worldlet degree σ, and calculate the curve of each network metric
Lower area AUC;AUC to each network metric and each network metric performs nonparametric displacement test determination in network attribute
In whether there is significant group difference, while by existing experimental data, coming judgement sample disease condition and information.
2. a kind of sacred disease analysis method based on brain network structure as claimed in claim 1, it is characterised in that:The step
In rapid three, various brain network parameters are calculated, including following parameter:
3.1 cluster coefficients:What cluster coefficient was weighed is the grouping of the world economy degree of network, is the important parameter for measuring network:
CiRepresent the cluster coefficient of node i, eiThe number on the side actually connected between representing the nodes neighbors,Express possibility
Maximum connection side number;
N represents network node quantity, and V represents whole network;
3.2 component efficiencies:
E (i) represents the component efficiency of arbitrary node i, GiRefer to the subgraph that the neighbours of node i are constituted,Table constitutes the number of subgraph
Amount, i, j are arbitrary node, l in networkjkRepresent node j, the shortest path length between k;
3.3 shortest path lengths:Shortest path features the optimal path of information another node of arrival of a certain node in network,
By shortest path can quickly transmission information, so as to save system resource;
3.4 global efficiencies:Usual shortest path length will carry out computing in some connected graph, because if existing in network
Disconnected node can cause the shortest path length value between the two nodes to be infinite, therefore propose:
3.5 centrads:Centrad is a statistical indicator for portraying nodes effect and status;
NbcI () is the centrad of any point i in network, G represents whole network, and i, j, k are arbitrary node, δ in networkjkIt is
The quantity of all shortest paths from node j to node k, δjkI () is the quantity by node i in these shortest paths;
3.6 worldlet degree, by worldlet degree, we can weigh the power of " worldlet " attribute, and δ, λ are calculated first;δ is represented
The ratio of real network and random network cluster coefficients, λ represents the ratio of real network and random network shortest path length;
CrealRepresent real network cluster coefficient, CrandomThe cluster coefficient of table random network.
LrealRepresent real network cluster coefficient, LrandomThe cluster coefficient of table random network.
Worldlet degree σ is defined as follows:
3. a kind of sacred disease analysis method based on brain network structure as claimed in claim 1 or 2, it is characterised in that:Institute
State in step 3, nonparametric replaces test mode:The group difference between the average value of each network metric is calculated first;For every
, randomly be re-assigned to all values in two groups by individual network metric, and recalculates average between two randomization groups
Difference;The randomization program repeat 10,000 time, and each distribution preceding 5% be used as with the null hypothesis that probable error is 0.05
One tailed test critical value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710107488.1A CN106845137A (en) | 2017-02-27 | 2017-02-27 | A kind of sacred disease analysis method based on brain network structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710107488.1A CN106845137A (en) | 2017-02-27 | 2017-02-27 | A kind of sacred disease analysis method based on brain network structure |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106845137A true CN106845137A (en) | 2017-06-13 |
Family
ID=59133503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710107488.1A Pending CN106845137A (en) | 2017-02-27 | 2017-02-27 | A kind of sacred disease analysis method based on brain network structure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106845137A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107680677A (en) * | 2017-10-11 | 2018-02-09 | 四川大学 | Neuropsychiatric disease sorting technique based on brain network analysis |
CN109920550A (en) * | 2018-12-25 | 2019-06-21 | 天津大学 | A method of teenager's lafora's disease is studied based on dMRI |
CN110298479A (en) * | 2019-05-20 | 2019-10-01 | 北京航空航天大学 | A kind of brain volume atrophy prediction technique based on brain function network |
CN111627553A (en) * | 2020-05-26 | 2020-09-04 | 四川大学华西医院 | Method for constructing individualized prediction model of first-onset schizophrenia |
CN112274145A (en) * | 2019-07-22 | 2021-01-29 | 苏州布芮恩智能科技有限公司 | Method and device for processing near-infrared brain function imaging data and storage medium |
CN113425312A (en) * | 2021-07-30 | 2021-09-24 | 清华大学 | Electroencephalogram data processing method and device |
CN117690537A (en) * | 2024-02-04 | 2024-03-12 | 中日友好医院(中日友好临床医学研究所) | Cross-modal method and device for QSM (quad slit) and brain atrophy and brain connection group association |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931281A (en) * | 2016-04-14 | 2016-09-07 | 中国人民解放军国防科学技术大学 | Method for quantitatively describing cerebral function network based on network characteristic entropy |
CN106021949A (en) * | 2016-05-31 | 2016-10-12 | 常州大学 | Function connection analysis method of brain default network |
CN106447023A (en) * | 2016-08-18 | 2017-02-22 | 常州大学 | Brain default network directed connection analysis method based on motif structures |
-
2017
- 2017-02-27 CN CN201710107488.1A patent/CN106845137A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931281A (en) * | 2016-04-14 | 2016-09-07 | 中国人民解放军国防科学技术大学 | Method for quantitatively describing cerebral function network based on network characteristic entropy |
CN106021949A (en) * | 2016-05-31 | 2016-10-12 | 常州大学 | Function connection analysis method of brain default network |
CN106447023A (en) * | 2016-08-18 | 2017-02-22 | 常州大学 | Brain default network directed connection analysis method based on motif structures |
Non-Patent Citations (6)
Title |
---|
AARONALEXANDER-BLOCH等: "The discovery of population differences in network community structure: New methods and applications to brain functional networks in schizophrenia", 《NEUROIMAGE》 * |
SOPHIE ACHARD等: "A Resilient, Low-Frequency, Small-World Human Brain Functional Network with Highly Connected Association Cortical Hubs", 《JOURNAL OF NEUROSCIENCE》 * |
刘薇: "磁共振弥散张量成像中的失真伪影校正及其相关算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
张方风: "复杂脑网络研究:现状与挑战", 《上海理工大学学报》 * |
王志江: "脑功能网络的拓扑结构及任务效应的功能磁共振成像研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
田丽霞: "基于图论的复杂脑网络分析", 《北京生物医学工程》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107680677A (en) * | 2017-10-11 | 2018-02-09 | 四川大学 | Neuropsychiatric disease sorting technique based on brain network analysis |
CN109920550A (en) * | 2018-12-25 | 2019-06-21 | 天津大学 | A method of teenager's lafora's disease is studied based on dMRI |
CN110298479A (en) * | 2019-05-20 | 2019-10-01 | 北京航空航天大学 | A kind of brain volume atrophy prediction technique based on brain function network |
CN110298479B (en) * | 2019-05-20 | 2021-09-03 | 北京航空航天大学 | Brain volume atrophy prediction method based on brain function network |
CN112274145A (en) * | 2019-07-22 | 2021-01-29 | 苏州布芮恩智能科技有限公司 | Method and device for processing near-infrared brain function imaging data and storage medium |
CN111627553A (en) * | 2020-05-26 | 2020-09-04 | 四川大学华西医院 | Method for constructing individualized prediction model of first-onset schizophrenia |
CN113425312A (en) * | 2021-07-30 | 2021-09-24 | 清华大学 | Electroencephalogram data processing method and device |
CN117690537A (en) * | 2024-02-04 | 2024-03-12 | 中日友好医院(中日友好临床医学研究所) | Cross-modal method and device for QSM (quad slit) and brain atrophy and brain connection group association |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106845137A (en) | A kind of sacred disease analysis method based on brain network structure | |
McKenzie et al. | Distributions of cherries for two models of trees | |
Green et al. | Effect of the number of scale points on chi‐square fit indices in confirmatory factor analysis | |
Cheng et al. | Optimization of seed density in DTI tractography for structural networks | |
CN107358022A (en) | A kind of Modularity analysis method of cerebral function network | |
CN105574191B (en) | Online community network multi-source point information source tracing system and method | |
Kennedy et al. | Why data coherence and quality is critical for understanding interareal cortical networks | |
CN105893637A (en) | Link prediction method in large-scale microblog heterogeneous information network | |
CN112365962B (en) | Medical resource allocation method and equipment based on lung CT image and clinical information | |
DE112012005087T5 (en) | Method, unit and computer program for visualizing risk estimates in event sequences | |
CN111710415A (en) | Whole brain oriented network analysis method based on Granger neuropathy | |
CN105228185A (en) | A kind of method for Fuzzy Redundancy node identities in identification communication network | |
Li et al. | Asymptotic analysis of endemic equilibrium to a brucellosis model | |
CN107358534A (en) | The unbiased data collecting system and acquisition method of social networks | |
CN103984794B (en) | Fixed-outline integrated circuit floorplanning method for soft module | |
Hu et al. | Hypernetwork models based on random hypergraphs | |
CN108470251A (en) | Community based on Average Mutual divides quality evaluating method and system | |
Mancini et al. | Network attack simulations in Alzheimer's disease: The link between network tolerance and neurodegeneration | |
Chen et al. | Predicting resting-state functional connectivity with efficient structural connectivity | |
JP2010191890A (en) | Method and program for inversely estimating od traffic volume between large zones using centroid | |
Preciado et al. | Moment-based spectral analysis of random graphs with given expected degrees | |
CN113706459B (en) | Detection and simulation repair device for abnormal brain area of autism patient | |
Xia et al. | Connectomic insights into topologically centralized network edges and relevant motifs in the human brain | |
Morcrette | Fully analyzing an algebraic Pólya urn model | |
CN104468196B (en) | Virtual network method for diagnosing faults and device based on evidence screening |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170613 |