The invention pertains to an angular resolved diffusion weighted magnetic resonance imaging method. In angular resolved diffusion magnetic resonance imaging method, magnetic resonance signals are acquired that are diffusion weighted. The diffusion weighting is effected by way of diffusion magnetic gradient fields. These diffusion weighted magnetic resonance signals are also spatially encoded by way of encoding magnetic gradient fields such as read gradient fields and phase encoding gradients. In particular in diffusion tensor imaging (DTI), diffusion weighting is performed for several spatial directions. From the diffusion weighted magnetic resonance signals and on the basis of a tensor analysis, local principal diffusion directions are derived for individual voxels. The diffusion process is a stochastic process of the population of nuclear (proton) spins and the tensor analysis derives the main diffusion directions into which the diffusive motion of the individual spins. These main directions correspond to the directions of the eigenvectors of the diffusion tensor and the main direction relating to the largest eigenvalue is the principal diffusion direction. This principal diffusion direction represents the direction in which diffusion mainly takes place in the voxel at issue. Information of diffusion directions and the apparent diffusion coefficients is useful to extract the directional fibre structure in neurological systems such as the human or animal brain and spinal cord.
Angular resolved diffusion magnetic resonance imaging is known from the paper ‘Characterization of anisotropy in high angular resolution diffusion-weighted MRI’ by L. R. Frank in MRM47 (2002)1083-1099.
The cited paper mentions the problem that multiple principal diffusion directions may appear in a single voxel and that characterisation of diffusion in such voxels becomes problematic. The known magnetic resonance imaging method applies methods of group theory to this problem to show that the measurements can be decomposed into irreducible representations of the rotation group in which isotropic, single fibre, multiple fibre components are separable direct sum subspaces. Multiple fibres passing through a single voxel are represented in a decomposition based on spherical harmonics and the state of a voxel with several fibres passing through it can be expressed as a direct sum of irreducible representations of the rotation group.
An object of the invention is to provide a high angular resolution diffusion-weighted magnetic resonance imaging method which requires less computational effort than the known magnetic resonance imaging method.
This object is achieved by the magnetic resonance imaging method of the invention comprising
According to the invention contributions from different principal diffusion directions (fibre directions) in individual voxels are distinguished on the basis of diffusion weighted magnetic resonance signals for several values of the diffusion weighting. The invention is based on the insight that the dependence of the signal level of the diffusion weighted magnetic resonance signals on the diffusion weighting is different in voxels where there is only one single principal diffusion direction as compared to voxels where there is a superposition of contributions from several principal diffusion directions. Even, the way the signal level of the diffusion weighted magnetic resonance signals depends on the applied diffusion weighting reflects the number of principal diffusion directions that occur in the voxel at issue. Hence, on the basis of a variation with the diffusion strength of the apparent diffusion coefficient voxels having a single diffusion direction are distinguished from voxels having several diffusion directions. That is, voxels through which fibres pass at different directions can be identified. Accordingly, in the further analysis of the object dataset account can be taken of voxels in which contributions of several principal diffusion directions occur. Notably, in these voxels a decomposition of contribution from the respective principal diffusion directions carried out.
These and other aspects of the invention will be further elaborated with reference to the embodiments defined in the dependent Claims.
From the magnetic resonance signals at several diffusion weightings, respective values of the apparent diffusion coefficients for individual voxels are computed. From these diffusion weighted dependent values, the contributions for separate principal diffusion directions can be computed for the voxel at issue. Hence, contributions to the apparent diffusion coefficient from various fibres passing through the voxel at issue are obtained. Accordingly, the local directional structure of fibres can be better resolved, even if several fibres are crossing at the voxel at issue.
It appears that in practice diffusion strengths for several principal diffusion directions do not vary substantially at the scale of an individual voxel. According to one aspect of the invention the apparent diffusion coefficient can be accurately decomposed into contributions to the respective identified principal diffusion directions. This decomposition can be made on the basis of the assumption that diffusion strengths are equal for the respective principal diffusion direction in the voxel at issue. This assumption appears often to be quite accurate, e.g. because an individual voxel usually pertains to a single type of tissue. This assumption substantially reduces the computational effort to compute the contributions to the apparent diffusion coefficient. The ultimate accuracy of the resolution of the local directional structure of the directional structure of the fibres is hardly affected.
The invention also pertains to a method of analysis of an object dataset as defined in claim 4. The method of analysis of an object dataset of the invention achieves to analyse the directional fibre structure separately from the acquisition of the magnetic resonance signals. That is, the patient can be scanned to acquired the magnetic resonance signal data and these data are later analysed to analyse the directional structure. This analysis can at option also be performed at a different location.
The invention also pertains to a computer programme as defined in claim 5. The computer programme of the invention can be installed in a general purpose workstation so as to enable the workstation to perform the method of analysis of the object dataset of the invention. This workstation may be separate from the magnetic resonance imaging system which acquires the diffusion weighted magnetic resonance signals.
The invention further relates to an magnetic resonance imaging system as defined in claim 6. The magnetic resonance imaging system of the invention comprises an image processing unit which carries out the method of the invention. Notably, in the image processing unit the computer programme of the invention is installed.
These and other aspects of the invention will be elucidated with reference to the embodiments described hereinafter and with reference to the accompanying drawing wherein
The degeneracy check 31 and the decomposition 32 of the apparent diffusion coefficient into components for the respective principal diffusion directions is based on the following considerations. In general the measured magnetic resonance signal intensity in multi-fibre situation is given as:
Here S(v) is the measured signal at voxel position v, S0 is the measured signal when no diffusion sensitation is applied, Dk is the 3×3 diffusion matrix for the k-th fibre, b represents the diffusion strength of the diffusion sensitive pulse sequence and fk is the volume fraction of the k-th fibre in that voxel. Because several measurement are made for b-values and diffusion directions, the quantities Dk and fk can be obtained, e.g. on the basis of a model fitting method.
Number | Date | Country | Kind |
---|---|---|---|
04100442.5 | Feb 2004 | EP | regional |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB05/50402 | 1/31/2005 | WO | 00 | 7/31/2006 |