CN110533758B - Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation - Google Patents
Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation Download PDFInfo
- Publication number
- CN110533758B CN110533758B CN201910680865.XA CN201910680865A CN110533758B CN 110533758 B CN110533758 B CN 110533758B CN 201910680865 A CN201910680865 A CN 201910680865A CN 110533758 B CN110533758 B CN 110533758B
- Authority
- CN
- China
- Prior art keywords
- ftd
- fiber
- voxel
- voxels
- continuity
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Graphics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Materials For Medical Uses (AREA)
Abstract
A brain fiber asymmetric reconstruction method based on a hydrodynamic differential equation, which regards fiber bundles as flow bundles and models fiber reconstruction by introducing related concepts in hydrodynamic; describing the continuity of the fiber within the voxel in terms of a divergence concept in the fluid mechanics; describing the spatial consistency of the fiber bundles by extending the concept of divergence between voxels; and calculating the microstructure reconstruction of the asymmetric brain fiber through the optimization of the two constraint conditions. The experiment is respectively compared with the current popular fiber reconstruction method in the aspects of simulated magnetic resonance data and actual clinical data, and the experimental result proves that the method provided by the invention improves the accuracy of fiber reconstruction.
Description
Technical neighborhood
The invention relates to medical imaging and neuroanatomical neighborhood under computer graphics, in particular to a brain fiber microstructure reconstruction method.
Background
The brain is a comprehensive organ for controlling the complex functional activities of human beings such as logic thinking, learning and memory, movement, emotion and the like, and the exploration of the working mechanism of the human brain is the front edge and hot spot of the current scientific research. The brain white matter fiber reconstruction technology is a unique non-invasive method for imaging the fiber tissue microstructure in the living brain at present by imaging the brain fiber space microstructure information with anatomical significance, and has become an important technical means for brain science researches such as brain cognitive mechanism exploration, nerve disease pathological analysis, brain operation navigation and the like.
Reconstruction of brain fiber microstructure is a fundamental step in brain fiber imaging, providing accurate fiber direction estimation for fiber bundle tracking. The traditional reconstruction method usually only depends on the monosomic information to reconstruct the microstructure of the voxels, and the model type is basically a centrosymmetric model, so that the accuracy of fiber direction reconstruction is limited. It is a hot spot of research to propose new, more accurate brain fiber microstructure reconstruction models.
Disclosure of Invention
In order to solve the problems of low dependence and low precision of the existing brain fiber microstructure reconstruction method on single element information and central symmetry of a model, the invention provides an asymmetric brain fiber microstructure reconstruction method based on a hydrodynamic differential equation by combining neighborhood information of voxels.
The technical scheme adopted for solving the technical problems is as follows:
a method for asymmetrically reconstructing brain fibers based on a hydrodynamic differential equation, comprising the steps of:
step one, asymmetric fiber trajectory distribution (fiber trajectory distribution, FTD) function:
considering the fiber bundles as a stream, the stream consists of a series of sets of streamlines s= { S i I=1, ····, n, any point in space (x, y, z) is the flow field at that point, the fiber direction at any point within a voxel is represented by a flow field:
FTD is represented as a flow field distribution throughout the voxel, approximated by a set of ternary quadratic polynomials:
υ(x,y,z)=AC(x,y,z) (2)
wherein the coefficient matrix a is defined as follows:
c (x, y, z) = [ x ] denoted by C 2 ,y 2 ,z 2 ,xy,xz,yz,x,y,z,1] T ;
Combining the spatial continuity of the fiber bundles on the neighborhood, wherein the process is as follows:
2.1 continuity constraints for intra-voxel FTD
Assuming that the diffusion displacement of water molecules in the same fiber bundle remains continuous, using the continuous incompressible fluid theory, the spatial continuity of the fiber trajectory is described by introducing the concept of divergence of the fiber flow on the diffusion tensor vector field:
when the FTD does not belong to the start or end region of the nerve fiber bundle, divΩ satisfies the following formula:
divΩ=0 (5)
the simultaneous formulas (2), (3), (4) and (5) are obtained:
2.2 spatial continuity constraints between FTD voxels while equation (6) ensures that the FTD within a voxel meets continuity, but does not indicate inter-voxel fiber bundle continuity, the same fiber bundle should be consistent between adjacent voxels, i.e., the corresponding FTD should meet inter-voxel continuity, an FTD inter-voxel consistency function is proposed to characterize fiber bundle continuity between adjacent voxels, assuming that the voxel is a unit cube, N c =(c 0 ,c 1 ,…,c 5 ) Representing six voxels adjacent to the center voxel c, denoted by A c =(A 0 ,A 1 ,…,A 5 ) FTD coefficients representing neighborhood voxels, for adjacent voxels c and c traversed by the stream i The continuity function divΓ (x, y, z) of any point of its intersection is defined as follows:
wherein Γ is the adjacent voxel connecting surface through which the flow passes, S is the flow, and v i C and c are respectively i FTD of (c) i Mapping to the same coordinate system as c:
wherein a is jk Anda and A respectively i The continuity function div ψ over the entire plane Γ is the area integral of divΓ (x, y, z):
each voxel has six adjacent voxels, and the continuity function among the voxels is the sum of the continuity functions of six connecting surfaces:
the combined formula (2) -formula (9) is obtained:
wherein a is jk Elements in FTD coefficient matrix a representing the center voxel,FTD coefficient matrix A representing the ith voxel of the neighborhood i Elements of (a) and (b);
step three, calculating FTD
Calculating the FTD by minimizing intra-voxel and inter-voxel continuity functions aims at making the fiber trajectory distribution most closely to the fiber direction distribution function (fiber orientation distribution, FOD), and the coefficient matrix a in the FTD can be calculated by optimizing the cost function as follows:
where Φ (v (x, y, z)) is the probability of FOD at point (x, y, z), we take the 26 neighborhood c= [ C ] of the central voxel for simplicity of calculation 1 ,c 2 ,…,c 26 ]Peak p= [ P ] of middle FOD 1 ,p 2 ,…,p 26 ]As an approximation of Φ (v (x, y, z)), the formula (11) is simplified as follows:
after the flow field coefficient A is obtained, the FTD of the voxel is obtained.
The beneficial effects of the invention are as follows: the accuracy of the fiber reconstruction is improved.
Drawings
FIG. 1 is a schematic diagram of a center voxel and a neighborhood voxel.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1, a method for asymmetrically reconstructing brain fibers based on a hydrodynamic differential equation includes the steps of:
step one, asymmetric fiber trajectory distribution (fiber trajectory distribution, FTD) function:
considering the fiber bundles as a stream, the stream consists of a series of sets of streamlines s= { S i I=1, ····, n, any point in space (x, y, z) is the flow field at that point, the fiber direction at any point within a voxel is represented by a flow field:
FTD is represented as a flow field distribution throughout the voxel, approximated by a set of ternary quadratic polynomials:
υ(x,y,z)=AC(x,y,z) (2)
wherein the coefficient matrix a is defined as follows:
c (x, y, z) = [ x ] denoted by C 2 ,y 2 ,z 2 ,xy,xz,yz,x,y,z,1] T ;
Combining the spatial continuity of the fiber bundles on the neighborhood,
since during dMRI data acquisition, the brain nerve fibers can be seen as unchanged, corresponding to related concepts in fluid mechanics, i.e. the stream is seen as an incompressible constant stream; the process is as follows:
2.1 continuity constraints for intra-voxel FTD
Assuming that the diffusion displacement of water molecules in the same fiber bundle remains continuous, using the continuous incompressible fluid theory, the spatial continuity of the fiber trajectory is described by introducing the concept of divergence of the fiber flow on the diffusion tensor vector field:
when the FTD does not belong to the start or end region of the nerve fiber bundle, divΩ satisfies the following formula:
divΩ=0 (5)
the simultaneous formulas (2), (3), (4) and (5) are obtained:
2.2 spatial continuity constraints between FTD voxels while equation (6) ensures that the FTD within a voxel meets continuity, but does not indicate inter-voxel fiber bundle continuity, the same fiber bundle should be consistent between adjacent voxels, i.e., the corresponding FTD should meet inter-voxel continuity, an FTD inter-voxel consistency function is proposed to characterize fiber bundle continuity between adjacent voxels, assuming that the voxel is a unit cube, N c =(c 0 ,c 1 ,…,c 5 ) Representing six voxels adjacent to the center voxel c, denoted by A c =(A 0 ,A 1 ,…,A 5 ) FTD coefficients representing neighborhood voxels, for adjacent voxels c and c traversed by the stream i The continuity function divΓ (x, y, z) of any point of its intersection) The definition is as follows:
wherein Γ is the adjacent voxel connecting surface through which the flow passes, S is the flow, and v i C and c are respectively i FTD of (c) i Mapping to the same coordinate system as c:
wherein a is jk Anda and A respectively i The continuity function div ψ over the entire plane Γ is the area integral of divΓ (x, y, z):
each voxel has six adjacent voxels, and the continuity function among the voxels is the sum of the continuity functions of six connecting surfaces:
the combined formula (2) -formula (9) is obtained:
wherein a is jk Elements in FTD coefficient matrix a representing the center voxel,FTD coefficient matrix A representing the ith voxel of the neighborhood i Elements of (a) and (b);
step three, calculating FTD
Calculating the FTD by minimizing intra-voxel and inter-voxel continuity functions aims at making the fiber trajectory distribution most closely to the fiber direction distribution function (fiber orientation distribution, FOD), and the coefficient matrix a in the FTD can be calculated by optimizing the cost function as follows:
where Φ (v (x, y, z)) is the probability of FOD at point (x, y, z), we take the 26 neighborhood c= [ C ] of the central voxel for simplicity of calculation 1 ,c 2 ,…,c 26 ]Peak p= [ P ] of middle FOD 1 ,p 2 ,…,p 26 ]As an approximation of Φ (v (x, y, z)), the formula (11) is simplified as follows:
after the flow field coefficient A is obtained, the FTD of the voxel is obtained.
Claims (1)
1. A method for asymmetric reconstruction of brain fibers based on hydrodynamic differential equations, the method comprising the steps of:
step one, an asymmetric fiber track distribution FTD function:
considering the fiber bundles as a stream, the stream consists of a series of sets of streamlines s= { S i I=1, ····, n, any point in space (x, y, z) is the flow field at that point, the fiber direction at any point within a voxel is represented by a flow field:
FTD is represented as a flow field distribution throughout the voxel, approximated by a set of ternary quadratic polynomials:
υ(x,y,z)=AC(x,y,z) (2)
wherein the coefficient matrix a is defined as follows:
c (x, y, z) = [ x ] denoted by C 2 ,y 2 ,z 2 ,xy,xz,yz,x,y,z,1] T ;
Combining the spatial continuity of the fiber bundles on the neighborhood, wherein the process is as follows:
2.1 continuity constraints for intra-voxel FTD
Assuming that the diffusion displacement of water molecules in the same fiber bundle remains continuous, using the continuous incompressible fluid theory, the spatial continuity of the fiber trajectory is described by introducing the concept of divergence of the fiber flow on the diffusion tensor vector field:
when the FTD does not belong to the start or end region of the nerve fiber bundle, divΩ satisfies the following formula:
divΩ=0 (5)
the simultaneous formulas (2), (3), (4) and (5) are obtained:
2.2 spatial continuity constraints between FTD voxels
Although formula (6) can ensure that the FTD in the voxels meets the continuity, but can not represent the fiber bundle continuity among the voxels, the same fiber bundle is consistent among the adjacent voxels, namely the corresponding FTD should meet the continuity among the voxels, a consistency function among the FTD voxels is proposed to characterize the continuity of the fiber bundle among the adjacent voxels, and the voxels are assumed to be unit cubes, and N is c =(c 0 ,c 1 ,…,c 5 ) Representing that the central voxels c are adjacentSix voxels, denoted A c =(A 0 ,A 1 ,…,A 5 ) FTD coefficients representing neighborhood voxels, for adjacent voxels c and c traversed by the stream i The continuity function divΓ (x, y, z) of any point of its intersection is defined as follows:
wherein Γ is the adjacent voxel connecting surface through which the flow passes, S is the flow, and v i C and c are respectively i FTD of (c) i Mapping to the same coordinate system as c:
wherein a is jk Anda and A respectively i I=0, …,5; j=1, 2,3; k=0, …,9, the continuity function div ψ over the whole plane Γ is the area integral of divΓ (x, y, z):
each voxel has six adjacent voxels, and the continuity function among the voxels is the sum of the continuity functions of six connecting surfaces:
the combined formula (2) -formula (9) is obtained:
wherein a is jk Elements in FTD coefficient matrix a representing the center voxel,FTD coefficient matrix A representing the ith voxel of the neighborhood i Elements of (a) and (b);
step three, calculating FTD
Calculating FTD by minimizing intra-voxel and inter-voxel continuity functions aims at making the fiber track distribution most closely to the fiber direction distribution function FOD, and the coefficient matrix a in FTD is calculated by optimizing the cost function:
where Φ (v (x, y, z)) is the probability that FOD is at point (x, y, z), the 26 neighborhood c= [ C ] of the center voxel is taken 1 ,c 2 ,…,c 26 ]Peak p= [ P ] of middle FOD 1 ,p 2 ,…,p 26 ]As an approximation of Φ (v (x, y, z)), the formula (11) is simplified as follows:
after the flow field coefficient A is obtained, the FTD of the voxel is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910680865.XA CN110533758B (en) | 2019-07-26 | 2019-07-26 | Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910680865.XA CN110533758B (en) | 2019-07-26 | 2019-07-26 | Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110533758A CN110533758A (en) | 2019-12-03 |
CN110533758B true CN110533758B (en) | 2023-06-27 |
Family
ID=68661810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910680865.XA Active CN110533758B (en) | 2019-07-26 | 2019-07-26 | Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110533758B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112489220A (en) * | 2020-10-23 | 2021-03-12 | 浙江工业大学 | Nerve fiber continuous tracking method based on flow field distribution |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7657071B2 (en) * | 2005-04-11 | 2010-02-02 | Siemens Medical Solutions Usa, Inc. | Systems, devices, and methods for bundle segmentation in diffusion tensor magnetic resonance imaging |
US7511495B2 (en) * | 2005-04-25 | 2009-03-31 | University Of Utah | Systems and methods for image reconstruction of sensitivity encoded MRI data |
US7623725B2 (en) * | 2005-10-14 | 2009-11-24 | Hewlett-Packard Development Company, L.P. | Method and system for denoising pairs of mutually interfering signals |
WO2012125829A2 (en) * | 2011-03-15 | 2012-09-20 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | Directional diffusion fiber tracking |
-
2019
- 2019-07-26 CN CN201910680865.XA patent/CN110533758B/en active Active
Non-Patent Citations (1)
Title |
---|
曹喆文.基于邻域非对称结构成像的神经纤维聚类算法研究.《中国优秀硕士学位论文全文数据库基础科学辑》.2019,(第undefined期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN110533758A (en) | 2019-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cerliani et al. | Probabilistic tractography recovers a rostrocaudal trajectory of connectivity variability in the human insular cortex | |
Forkel et al. | The anatomy of fronto-occipital connections from early blunt dissections to contemporary tractography | |
Pullan et al. | Mathematically modelling the electrical activity of the heart: from cell to body surface and back again | |
Oosterhof et al. | Surface-based information mapping reveals crossmodal vision–action representations in human parietal and occipitotemporal cortex | |
Nowicka et al. | Transcallosal transfer of information and functional asymmetry of the human brain | |
Miller et al. | Diffeomorphometry and geodesic positioning systems for human anatomy | |
Miller et al. | The emerging discipline of computational functional anatomy | |
Makris et al. | Human cerebellum: surface-assisted cortical parcellation and volumetry with magnetic resonance imaging | |
CN104392019A (en) | High-order diffusion tensor mixed sparse imaging method for brain white matter fiber tracking | |
Zhang et al. | Three-dimensional reconstruction of medical images based on 3D slicer | |
CN107145702A (en) | Cardiogenic embolic type of stroke Risk Forecast System and method based on medical image | |
Wei et al. | Edge sparse basis network: a deep learning framework for EEG source localization | |
CN110533758B (en) | Brain fiber asymmetric reconstruction method based on hydrodynamic differential equation | |
Sun et al. | The forward EEG solutions can be computed using artificial neural networks | |
Johnson et al. | Computational medicine: Bioelectric field problems | |
Schwartz et al. | Evolution of cortical geometry and its link to function, behaviour and ecology | |
Tay et al. | Towards patient-specific cardiovascular modeling system using the immersed boundary technique | |
Hyde et al. | Evaluation of numerical techniques for solving the current injection problem in biological tissues | |
Zhu et al. | Tracking the adaptation and compensation processes of patients’ brain arterial network to an evolving glioblastoma | |
Deng et al. | Bionic design model for co-creative product innovation based on deep generative and BID | |
Song et al. | The effect of deep brain structure modeling on transcranial direct current stimulation-induced electric fields: An in-silico study | |
Guo et al. | Analysis of the elastic stress for the bifurcated region of blood vessel | |
Kermi et al. | 3D-Computerized facial reconstructions from 3D-MRI of human heads using deformable model approach | |
Ren et al. | Modeling hierarchical spatial and temporal patterns of naturalistic fMRI volume via volumetric deep belief network with neural architecture search | |
Tricoche et al. | Visual analysis of bioelectric fields |
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 |