# Patent

## Title
Method and MRI device to detect a direction of at least one fiber in a body

## Abstract

Method to extract a diffusion direction from DW-MR signals, in which a magnetic field gradient qi is applied to a body part, a magnetic resonance signal Si is measured and a diffusion direction μj is determined, wherein for a predetermined positive number m of diffusion direction(s) μj, the signals Si are linked to the probability density function of the random vector X of water molecules displacements represented as a mixture of m probability density functions pXj of random vectors Xj characterizing the diffusion along μj, the Xj being decomposed into a first component modelling directional variability of the water molecules displacements with μj as mean direction, and into a second component modelling the length variability of the water molecules displacements along μj.

## Background
This is a non-provisional application claiming the benefit of U.S. Provisional Application No. 61/417,382, filed Nov. 26, 2010.


BACKGROUND OF THE INVENTION
1. Field of the Invention
The invention concerns a method to extract at least one diffusion direction from diffusion-weighted magnetic resonance signals (or DW-MR images for Diffusion-Weighted Magnetic Resonance images) of any part of a body that contains water molecules, comprising
exposing the body to a magnetic resonance imaging process, in which at least one magnetic field gradient is applied to the entire body part, at least one magnetic resonance signal is measured in each voxel of said body part and at least one diffusion direction is determined based on said magnetic resonance signal.
Various methods to extract diffusion directions from DW-MRI data are known in the prior art.
2. Description of Related Art
Two main families of methods are known: on the one hand methods using diffusion tensor imaging DTI according to documents [1], [2] and [3], and on the other hand methods using high angular resolution diffusion imaging NARDI according to documents [4], [5] and [6].
Document [1] U.S. Pat. No. 5,539,310 describes a magnetic resonance spectroscopy system comprising: means for applying a static magnetic field to an object under examination; means for applying magnetic field gradients to the object; means for applying a high frequency pulse to the object; detecting means for detecting a magnetic resonance signal from the object; control means for ordering said means to apply magnetic field gradients, and means for applying a high frequency pulse to operate in a predetermined sequence so as to induce magnetic resonance within the object, wherein said control means includes means for generating a plurality of echo signals corresponding to non-collinear directions of said magnetic field gradients; and means for processing data based on the magnetic resonance signal detected by said detecting means to obtain magnetic resonance information, wherein said means for processing data provides an expression of each said echo signal in terms of a diffusion tensor that summarizes the displacements of water molecules due to the diffusion phenomenon primed by said magnetic field gradients.
Document [1] also provides an estimation of the diffusion tensor.
Document [2] US-2005/0007100 proposes an alternative method for the estimation of the diffusion tensor and furthermore q-space MRI specimen characterization. In q-space MRI, a sample portion is situated in a static magnetic field, typically for a duration long enough to permit some or all spins of one or more species like for example water molecules in the sample portion to align with the static magnetic field. It provides another method to estimate the diffusion tensor from DW-MR images, comprising: obtaining a set of diffusion-weighted magnetic resonance signals associated to a plurality of both magnetic field gradients strengths and directions; and estimating the values of the elements of the tensor D by application of a series of gradient pulses in different directions. According to the same, magnetic resonance methods include modelling magnetic resonance signals obtained from specimens at low and high q-values to obtain parameters associated with specimen structure (wherein brain white matter is a representative anisotropic specimen) and orientation. In evaluation of brain white matter specimens, diffusion within axons is modelled as restricted diffusion along the axon axis according to Neumann's model and as hindered diffusion in the plane perpendicular to the axon axis according to the Gaussian model.
Document [3] U.S. Pat. No. 6,526,305 concerns a method of fiber reconstruction employing data acquired by magnetic resonance imaging. A method of creating an image of brain white matter fibers includes exposing the brain white matter fibers to a magnetic resonance imaging process. The data acquisition from the magnetic resonance imaging includes the acquisition of diffusion-weighted magnetic resonance images that are later employed to calculate an apparent diffusion coefficient in each voxel. The data is introduced into a microprocessor which calculates the six parameters of a three-dimensional diffusion tensor and deduces its three eigenvalues and its three eigenvectors. These eigenvalues along with their associated eigenvectors are subjected to further microprocessing to generate a plurality of images representing a plurality of diffusion properties of the fibers. The process in a preferred embodiment includes that the input of an initiation pixel begins the reconstruction process which involves pixel connecting and ultimately judgement or a decision based upon the tracking reaching termination of the fiber in each direction. If comparison in the computer results in the conclusion that the end of the fiber has been reached (“yes”), that is the end of the process as to that axonal fiber. If not (“no”), the fiber reconstruction process continues. A method of creating an image of individual brain white matter fibers is provided comprising exposing said brain white matter fibers to a DW-MR imaging process, introducing data acquired from said magnetic resonance imaging process into a computer, estimating a diffusion tensor to obtain a plurality of scalar values and a plurality of vectors from said data, employing said values and said vectors to initiate brain fiber tracking at a pixel of interest, and continuing or terminating said fiber tracking based upon a determination regarding whether the fiber is continuous or terminated based on randomness of fiber orientation of adjacent pixels.
One problem of DTI-based methods using estimation of the diffusion tensor according to documents [1], [2] and [3] is that they do not permit to detect multiple fibres when the fibres are crossing. Another inconvenience of DTI is that it provides a had estimation of the diffusion direction when the DW-MR images are too noisy.
Document [4] US-2008/0252291 concerns a High Angular Resolution Diffusion Weighted MRI (HARDI). A magnetic resonance imaging method involves acquisition of magnetic resonance signals with application of diffusion weighting at a plurality of diffusion weighting strengths and a plurality of diffusion directions. The occurrence of one single or several diffusion directions is identified for individual voxels. In this way account is taken of crossing fibres. However, the identification of crossing of fibres is detrimental to the acquisition load and time spent, because it requires a large amount of magnetic resonance (MRI) signals resulting from the application of a huge number of MR gradients.
Document [5] “Von Mises-Fisher Mixture Model of the diffusion ODF”, Tim Mc Graw, Baba C. Vemuri, Bob Yeziersky, Thomas Mareci, concerns high angular resolution diffusion imaging (HARDI) that allows the computation of water molecules displacements probabilities in the whole three-dimensional space. This probability function is referred to as the orientation distribution function ODF which is known as the reverse Fourier transform of the DW-MR image. The ODF over the sphere is obtained by integrating the displacements probability over the radial component. The latter ODF over the sphere is modelled by means of a mixture of von Mises-Fisher distributions.
There is a high loss of information in the method of document [5], because two digital integrations are necessary. This high loss of information can only be compensated by making a greater number of MRI acquisitions, which are then more time-consuming. Furthermore, the ODF is sampled over a set of arbitrary directions and requires a learning algorithm, which makes it not easily and systematically able to be generalised to every clinical cases.
Document [6] “Hyperspherical von Mises-Fisher Mixture (HvMF) Modelling of High Angular Resolution Diffusion MRI”, Ahhir Rhalerao, Carl-Fredrik Westin, MICCAI 2007, Part 1, LNCS 4791, pp 236-243, 2007, describes a mapping of unit vectors onto a 5D hypersphere to model and partition ODFs from HARDI data and makes a link to interpretation of the second order spherical harmonic decompositions of HARDI data.
The problem of techniques in documents [4], [5] and [6] using HARDI is that these methods require the acquisition of a large amount of magnetic resonance (MRI) signals resulting from the application of a huge number of MR gradients. A lower bound of this number is nowadays around 60, leading to a DW-MRI sequence of about 10 to 12 minutes. It is thus not easily applicable to clinical routine where lengthy acquisitions are often not possible on patients suffering from severe pathologies (e.g. psychiatry, multiple sclerosis, Parkinsonian disorders, paediatrics).
BRIEF SUMMARY OF THE INVENTION
A goal of the invention is to provide a method and device to detect diffusion direction(s) in a body, avoiding the inconveniences of the prior art, especially being more robust to noise with respect to DTI-based methods and reducing the acquisition time with respect to HARDI-based methods, and allowing to calculate more than one diffusion direction.
The invention provides a method to extract at least one diffusion direction from diffusion-weighted magnetic resonance signals of any part of a body that contains water molecules, comprising
exposing the body to a magnetic resonance imaging process, in which at least one magnetic field gradient qi is applied to the entire body part, at least one magnetic resonance signal Si is measured in each voxel of said body part and at least one diffusion direction μj is determined based on said magnetic resonance signal Si,
characterized in that
for a predetermined positive number m of diffusion direction(s) μj to be identified in each voxel of said body part with m≧1 said measured magnetic resonance signals Si are linked to the probability density function of the random vector X of water molecules displacements represented as a mixture of m probability density functions pXj of random vectors Xj, each pXj characterizing the diffusion along direction μj, these random vectors Xj being decomposed into a first component modelling directional variability of the water molecules displacements with μj as mean direction, and into a second component that models the length variability of the water molecules displacements along said mean direction μj,
wherein each of the m diffusion direction(s) μj is derived from N magnetic resonance signals Si acquired for respectively a plurality of N non-collinear magnetic field gradients qi applied to said body part, N being a positive integer≧2.
According to an embodiment of the invention, the first component is a 3-dimensional resealed random unit vector with μj as mean direction.
According to an embodiment of the invention, the first component is equal to γ multiplied by Vj, wherein
Vj is a 3-dimensional random unit vector following a Von Mises & Fisher distribution with mean direction μj and concentration parameter κj around said mean direction μj,
γ is the mean distance covered by the water molecules diffusing along any diffusion direction μj during the given diffusion time τ.
According to an embodiment of the invention, the second component Zj is a 3-dimensional random vector whose length varies along the mean direction μj.
According to an embodiment of the invention, Vj and Zj are independent.
According to an embodiment of the invention, the second component is equal Zj, wherein Zj is a 3-dimensional centered Gaussian random vector whose covariance matrix Σj is equal to




  
    
      ∑
      j
    
    ⁢
    
      
    
    ⁢
    
      =
      
        
          
            γ
            2
          
          
            
              κ
              j
            
            +
            1
          
        
        ⁢
        
          (
          
            
              I
              3
            
            +
            
              
                κ
                j
              
              ⁢
              
                μ
                j
              
              ⁢
              
                μ
                j
                ′
              
            
          
          )
        
      
    
  
  ,




wherein I3 is the 3×3 identity matrix.

According to an embodiment of the invention, said measured magnetic resonance signals Si is the noisy version of the magnitude of DW-MR complex signal A(qi), said magnitude being denoted by Φ(qi) and linked to the diffusion parameter θ through



  
    
      
    
    ⁢
    
      
        φ
        ⁡
        
          (
          
            
              q
              i
            
            |
            Θ
          
          )
        
      
      =
      
        
          φ
          0
        
        ⁢
        
          
            
              
                (
                
                  
                    ∑
                    
                      j
                      =
                      1
                    
                    M
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      f
                      j
                    
                    ⁢
                    
                      
                        R
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            
                              -
                              
                                q
                                i
                              
                            
                            |
                            
                              θ
                              j
                            
                          
                          ,
                          γ
                        
                        )
                      
                    
                  
                
                )
              
              2
            
            +
            
              
                (
                
                  
                    ∑
                    
                      j
                      =
                      1
                    
                    M
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      f
                      j
                    
                    ⁢
                    
                      
                        I
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            
                              -
                              
                                q
                                i
                              
                            
                            |
                            
                              θ
                              j
                            
                          
                          ,
                          γ
                        
                        )
                      
                    
                  
                
                )
              
              2
            
          
        
      
    
  




  
    
      
    
    ⁢
    with
  




  
    
      
        R
        j
      
      ⁡
      
        (
        
          
            
              -
              
                q
                i
              
            
            |
            
              θ
              j
            
          
          ,
          γ
        
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        sin
                        ⁡
                        
                          (
                          
                            γ
                            ⁢
                            
                              
                                
                                  
                                    
                                    
                                      q
                                      i
                                    
                                    
                                  
                                  2
                                
                                -
                                
                                  
                                    (
                                    
                                      
                                        κ
                                        j
                                      
                                      /
                                      γ
                                    
                                    )
                                  
                                  2
                                
                              
                            
                          
                          )
                        
                      
                      
                        γ
                        ⁢
                        
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            -
                            
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  /
                                  γ
                                
                                )
                              
                              2
                            
                          
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          
                            
                            
                              q
                              i
                            
                            
                          
                          2
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        
                          
                            α
                            j
                          
                          ⁢
                          sinh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          cos
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                        +
                        
                          
                            β
                            j
                          
                          ⁢
                          
                            
                          
                          ⁢
                          cosh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          sin
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                      
                      
                        
                          α
                          j
                          2
                        
                        +
                        
                          β
                          j
                          2
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          (
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            +
                            
                              
                                
                                  κ
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      μ
                                      j
                                      ′
                                    
                                    ⁢
                                    
                                      q
                                      i
                                    
                                  
                                  )
                                
                              
                              2
                            
                          
                          )
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
          
          ⁢
          
            

          
          ⁢
          
            
          
          ⁢
          and
          ⁢
          
            

          
          ⁢
          
            
              I
              j
            
            ⁡
            
              (
              
                
                  
                    -
                    
                      q
                      i
                    
                  
                  |
                  
                    θ
                    j
                  
                
                ,
                γ
              
              )
            
          
        
        =
        
          {
          
            
              
                
                  0
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        
                          
                            α
                            j
                          
                          ⁢
                          cosh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          sin
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                        -
                        
                          
                            β
                            j
                          
                          ⁢
                          
                            
                          
                          ⁢
                          sinh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          cos
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                      
                      
                        
                          α
                          j
                          2
                        
                        +
                        
                          β
                          j
                          2
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          (
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            +
                            
                              
                                
                                  κ
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      μ
                                      j
                                      ′
                                    
                                    ⁢
                                    
                                      q
                                      i
                                    
                                  
                                  )
                                
                              
                              2
                            
                          
                          )
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
          
        
      
    
  



wherein qi|Θ means “given the diffusion parameters θ”,

θj={μj,κj}, f or any j=1, . . . , M and Θ={θ1, . . . , θM, γ, f1, . . . , fM-1}.

wherein Φ0 is a given nominal value of the magnetic resonance signal Si when no gradient is applied,

wherein fj is the volume fraction of water molecules that diffuse along diffusion direction μj, with 0≦fj≦1 and





  
    
      
        ∑
        
          j
          =
          1
        
        m
      
      ⁢
      
        
      
      ⁢
      
        f
        j
      
    
    =
    1
  
  ,




wherein μ′j is the transposition or vector μj,




  
    
      Ω
      ij
    
    =
    
      {
      
        
          
            
              -
              
                q
                i
              
            
            ∈
            
                            3
            
          
          :
          
            
              
              
                -
                
                  q
                  i
                
              
              
            
            ≥
            
              
                κ
                j
              
              /
              γ
            
          
        
        ,
        
          
            and
            -
            
              q
              i
            
          
          ∈
          
            
              (
              
                span
                ⁡
                
                  (
                  
                    μ
                    j
                  
                  )
                
              
              )
            
            ⊥
          
        
      
      ]
    
  




  
    
      α
      ij
    
    =
    
      
        
          
            Re
            ⁢
            
              
            
            ⁢
            
              z
              ij
            
          
          +
          
            
            
              z
              ij
            
            
          
        
        2
      
    
  




  
    
      β
      ij
    
    =
    
      
        Im
        ⁢
        
          
        
        ⁢
        
          z
          ij
        
      
      
        
          2
          ⁢
          
            (
            
              
                Re
                ⁢
                
                  
                
                ⁢
                
                  z
                  ij
                
              
              +
              
                
                
                  z
                  ij
                
                
              
            
            )
          
        
      
    
  




  
    
      z
      ij
    
    =
    
      
        κ
        j
        2
      
      -
      
        
          γ
          2
        
        ⁢
        
          
            
            
              q
              i
            
            
          
          2
        
      
      -
      
        2
        ⁢
        i
        ⁢
        
          
        
        ⁢
        
          κ
          j
        
        ⁢
        
          γμ
          j
          ′
        
        ⁢
        
          q
          i
        
      
    
  



(span(μj))L is the space orthogonal to (span(μj)), where (span(μj)) is the vector sub-space of 3 generated by unit vector μj.

It means that any vector of 3 can be decomposed into the sum of a vector belonging to (span(μj)) and of a vector belonging to (span(μj))L.
According to an embodiment of the invention, the concentration parameters are linked to the volume fractions fj by the following relationship

κj=κ·fj 

where κ is a measure of anisotropy in the voxel, with




  κ
  =
  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      
        κ
        j
      
      .
    
  




According to an embodiment of the invention, a multiplicity of said voxels is provided, a magnetic resonance image of said voxels is reconstructed based at least on the values of the diffusion directions {μj}=j=1, . . . , m of the m fibres calculated for each voxel.
The invention provides also an MRI device to carry out the method to extract at least one diffusion direction from diffusion-weighted magnetic resonance signals of any part of a body that contains water molecules as described above, comprising
an emitter-receiver exposing the body part to a diffusion-weighted magnetic resonance imaging process, comprising at least one emitter to apply at least one magnetic field diffusion gradient qi to the entire body part, and at least one receiver to measure at least one magnetic resonance signal Si in each voxel of said body part and a calculator to compute at least one diffusion direction μj based on said magnetic resonance signal Si,
characterized in that
in said calculator, for a predetermined positive number m of diffusion direction(s) μj to be identified in each voxel of said body part with m≧1, said measured magnetic resonance signals Si are linked to the probability density function of the random vector X of water molecules displacements represented as a mixture of m probability density functions pXj of random vectors Xj, each pXj characterizing the diffusion along direction μj, these random vectors Xj being assumed independent and being decomposed into a first component modelling directional variability of the water molecules displacements with as mean direction, and into a second component that models the length variability of the water molecules displacements along said mean direction μj,
wherein in the calculator each of the m diffusion direction(s) μj is derived from N magnetic resonance signals Si measured for respectively a plurality of N non-collinear magnetic field gradients qi applied to said body part, N being a positive integer≧2.



BRIEF DESCRIPTION OF THE DRAWINGS
The invention will be better understood on reading the following description giving solely as an example of embodiment in reference to the drawings in which:
FIG. 1 shows an example of the first component according to the invention, wherein the horizontal and vertical axis are graduated in meters,
FIG. 2 shows an example of the second component according to the invention, wherein the horizontal and vertical axis are graduated in meters,
FIG. 3 shows the result of the first and second components according to FIGS. 1 and 2.
FIG. 4 is a flow chart showing the steps performed according to one embodiment of the invention.
FIG. 5 is a block diagram showing elements according to one embodiment of the invention.



DETAILED DESCRIPTION OF THE INVENTION
The principle of Diffusion-Weighted Magnetic Resonance Imaging (DW-MRI) is to apply sequentially several magnetic field diffusion gradients qi, and to measure, for each of them, their impact on the displacements of water molecules along at least one fibre through the magnetic resonance signal Si, in each voxel of the image. It is provided that N non-collinear magnetic field diffusion gradients qi are successively applied with i being an integer, which indexes each magnetic field diffusion gradient q, and which takes N different values with 1≦i≦N and N≧2, N being a positive integer.
The body to which the magnetic field diffusion gradients are applied is for example the brain of a living human being, in order to detect the orientation of nervous fibres in the brain, constituting the so-called brain white matter, which comprise water molecules and axons. The body can be any part of the living human being provided that there are water molecules, like for example the liver or the heart. The body can be also the body of a living animal.
For each given voxel of the body part, the MR imager gives a complex response signal A(qi) to each magnetic field diffusion gradient qi. The modulus S (qi) of the complex response signal A(qi) is referred to as the measured magnetic resonance signal or MR signal.
The complex measurement A(qi) is the Fourier transform of the probability density function (pdf) pX(x) of the random vector X of water molecules displacements in said voxel according to the following equation (1):




  
    
      A
      ⁡
      
        (
        
          q
          i
        
        )
      
    
    
      A
      ⁡
      
        (
        0
        )
      
    
  
  =
  
    
      
        ∫
        
          R
          3
        
        
          
        
      
      ⁢
      
        
          
            p
            X
          
          ⁡
          
            (
            x
            )
          
        
        ⁢
        exp
        ⁢
        
          {
          
            
              -
              i
            
            ⁢
            
              
            
            ⁢
            
              q
              i
              ′
            
            ⁢
            x
          
          }
        
        ⁢
        
          d
          x
        
      
    
    =
    
      
        φ
        X
      
      ⁡
      
        (
        
          -
          
            q
            i
          
        
        )
      
    
  




wherein φX is the characteristic function of the 3-dimensional random vector X of water molecules displacements.

According to the invention, m fibre(s) direction(s) are searched in the voxel or in each of the voxels, wherein m is a pre-determined positive integer equal to one or greater than one, m≧1.
It is considered that each fibre j (with j being an integer between 1 and m) extends along an associated diffusion direction μj in the considered voxel, i.e. dimensional random unit vector in the voxel. Diffusion directions μj are assumed to coincide with fibre directions.
The invention provides means to find the diffusion/fibre directions {μj}j=1, . . . , m directly from the acquired signals A{qi)}i=1, . . . , N.
According to the invention, the random vector Xj of the water molecules displacements along each fibre j is as follows:

Xj=γ·Vj+Zj 

The first component γ·Vj of the random vector Xj of water molecules displacements i a 3-dimensional random vector having direction μj as mean direction.
Vj is a 3-dimensional random unit vector distributed as a Von Mises & Fisher random unit vector having diffusion direction μj as mean direction and concentration parameter κj around said mean direction μj.
γ is the mean distance covered by water molecules diffusing along any diffusion direction μj during the given diffusion time τ.
The second component Zj of the random vector Xj of water molecules displacements is a 3-dimensional random vector whose length varies along the diffusion direction μj, as shown for example on FIG. 2, in which the points are simulations of Zj.
In an embodiment, N≧4 m. This allows to estimate not only the diffusion direction μj but also the concentration parameter κj and the mean distance γ.
The second component is equal to Zj, wherein Zj is a 3-dimensional centered Gaussian random vector whose covariance matrix Σj is equal to




  
    
      ∑
      j
    
    ⁢
    
      
    
    ⁢
    
      =
      
        
          
            γ
            2
          
          
            
              κ
              j
            
            +
            1
          
        
        ⁢
        
          (
          
            
              I
              3
            
            +
            
              
                κ
                j
              
              ⁢
              
                μ
                j
              
              ⁢
              
                μ
                j
                ′
              
            
          
          )
        
      
    
  
  ,




wherein I3 is the 3×3 identity matrix, and Vj and Zj are independent.

The displacement vector Xj of the water molecules is the sum of two independent components:
the first component γ·Vj takes into account the directional variability of water molecules displacements by estimating the diffusion/fibre direction μj on the sphere given the covered mean diffusion distance as radius;
the second component Zj takes into account the length variability of water molecules displacements by estimating the covered mean diffusion distance γ given the diffusion/fibre direction μj as direction along which water molecules diffuse.
The overall estimation is performed by the estimation of these two parts of the above equation simultaneously.
According to an embodiment, the role of κj is to provide a measure of anisotropy along the diffusion direction μj. Parameter κj evaluates the concentration of the data around the mean direction μj. The greater the value of κj, the higher the concentration of the data around the mean direction μj, as shown for example on FIG. 1, in which the simulated displacement direction γ. Vj of water molecules is shown around the mean direction μj. Several κj's, but not all of them, equal to 0, means that the number m of diffusion directions, i.e. the number m of fibres is over estimated, and there is no need to consider so many diffusion directions {μj}j=1, . . . , m. Therefore, the expected water molecules displacement Xj towards these directions μj for which κj=0 must be and is indeed null. All the κj's equal to 0 means that the region is isotropic and the proposed model is coherent with this situation in the sense described below: the MRI signal is invariant to changes of directions of gradients qi, as there are no preferred directions anymore. Finally, when κj>0 for every j, the expected water molecules displacement Xj towards any direction is equal to γ·μj as expected, i.e., the mean covered distance multiplied by the mean direction.
The sum of the first component of FIG. 1 and of the second component of FIG. 2 gives the displacement vector Xj of FIG. 3.
In an improvement, the concentration parameters are linked to the volume fractions fj by the following relationship

κj=κ·fj 

where κ is a measure of anisotropy in the voxel, with




  κ
  =
  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      κ
      j
    
  




With this improvement, sensitivity to noise is furthermore reduced.
These parameters {μj}j=1, . . . , m, {κj}j=1, . . . , m and γ are part of the definition of the probability density function (pdf) pX of water molecules displacements Xj due to the diffusion process. The parameters are extracted as follows.
In the context of multi-fibres detection, it is convenient to propose a mixture model for pX. It consists in modelling pX according to




  
    
      p
      X
    
    ⁡
    
      (
      x
      )
    
  
  =
  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      
        f
        j
      
      ⁢
      
        
          p
          Xj
        
        ⁡
        
          (
          x
          )
        
      
    
  




wherein fj is the volume fraction of water molecules that diffuse along diffusion direction μj or along fibre j with 0≦fj≦1 and




  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      f
      j
    
  
  =
  1.




    
    
        This assumption induces the following relationship between the MRI signal and the mixed densities pXj:
    
    






  
    
      A
      ⁡
      
        (
        
          
            q
            i
          
          |
          Θ
        
        )
      
    
    
      A
      ⁡
      
        (
        0
        )
      
    
  
  =
  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      
        f
        j
      
      ⁢
      
        
          φ
          Xj
        
        ⁡
        
          (
          
            
              
                -
                
                  q
                  i
                
              
              |
              
                θ
                j
              
            
            ,
            γ
          
          )
        
      
    
  




where {φXj}j=1, . . . , m are the characteristic functions of the 3-dimensional random vectors Xj, assumed independent and having probability density function pXj,
where

θj={μj,κj}, f or any j=1, . . . , M and Θ={θ1, . . . , θM, γ, f1, . . . , fM-1}.

Denote by Rj and Ij the real and imaginary parts of φXj, respectively. Then,





  
    
      Re
      ⁡
      
        [
        
          
            A
            ⁡
            
              (
              
                
                  q
                  i
                
                |
                Θ
              
              )
            
          
          
            A
            ⁡
            
              (
              0
              )
            
          
        
        ]
      
    
    =
    
      
        ∑
        
          j
          =
          1
        
        M
      
      ⁢
      
        
      
      ⁢
      
        
          f
          j
        
        ⁢
        
          
            R
            j
          
          ⁡
          
            (
            
              
                
                  -
                  
                    q
                    i
                  
                
                |
                
                  θ
                  j
                
              
              ,
              γ
            
            )
          
        
      
    
  
  ,
  
    

  
  ⁢
  
    
      Im
      ⁡
      
        [
        
          
            A
            ⁡
            
              (
              
                
                  q
                  i
                
                |
                Θ
              
              )
            
          
          
            A
            ⁡
            
              (
              0
              )
            
          
        
        ]
      
    
    =
    
      
        ∑
        
          j
          =
          1
        
        M
      
      ⁢
      
        
      
      ⁢
      
        
          f
          j
        
        ⁢
        
          
            I
            j
          
          ⁡
          
            (
            
              
                
                  -
                  
                    q
                    i
                  
                
                |
                
                  θ
                  j
                
              
              ,
              γ
            
            )
          
        
      
    
  
  ,




The derivation of the model is then given by the following equations:




  
    
      R
      j
    
    ⁡
    
      (
      
        
          
            -
            
              q
              i
            
          
          |
          
            θ
            j
          
        
        ,
        γ
      
      )
    
  
  =
  
    {
    
      
        
          
            
              
                
                  
                    
                      κ
                      j
                    
                    
                      sinh
                      ⁢
                      
                        
                      
                      ⁢
                      
                        κ
                        j
                      
                    
                  
                  ⁢
                  
                    
                      sin
                      ⁡
                      
                        (
                        
                          γ
                          ⁢
                          
                            
                              
                                
                                  
                                  
                                    q
                                    i
                                  
                                  
                                
                                2
                              
                              -
                              
                                
                                  (
                                  
                                    
                                      κ
                                      j
                                    
                                    /
                                    γ
                                  
                                  )
                                
                                2
                              
                            
                          
                        
                        )
                      
                    
                    
                      γ
                      ⁢
                      
                        
                          
                            
                              
                              
                                q
                                i
                              
                              
                            
                            2
                          
                          -
                          
                            
                              (
                              
                                
                                  κ
                                  j
                                
                                /
                                γ
                              
                              )
                            
                            2
                          
                        
                      
                    
                  
                  ⁢
                  
                    e
                    
                      
                        -
                        
                          
                            γ
                            2
                          
                          
                            2
                            ⁢
                            
                              (
                              
                                
                                  κ
                                  j
                                
                                +
                                1
                              
                              )
                            
                          
                        
                      
                      ⁢
                      
                        
                          
                          
                            q
                            i
                          
                          
                        
                        2
                      
                    
                  
                
                ,
              
            
            
              
                
                  
                    q
                    i
                  
                  ∉
                  
                    Ω
                    ij
                  
                
                ,
              
            
          
          
            
              
                
                  
                    κ
                    j
                  
                  
                    sinh
                    ⁢
                    
                      
                    
                    ⁢
                    
                      κ
                      j
                    
                  
                
                ⁢
                
                  
                    
                      
                        α
                        j
                      
                      ⁢
                      sinh
                      ⁢
                      
                        
                      
                      ⁢
                      
                        α
                        j
                      
                      ⁢
                      cos
                      ⁢
                      
                        
                      
                      ⁢
                      
                        β
                        j
                      
                    
                    +
                    
                      
                        β
                        j
                      
                      ⁢
                      
                        
                      
                      ⁢
                      cosh
                      ⁢
                      
                        
                      
                      ⁢
                      
                        α
                        j
                      
                      ⁢
                      sin
                      ⁢
                      
                        
                      
                      ⁢
                      
                        β
                        j
                      
                    
                  
                  
                    
                      α
                      j
                      2
                    
                    +
                    
                      β
                      j
                      2
                    
                  
                
                ⁢
                
                  e
                  
                    
                      
                        -
                        
                          
                            γ
                            2
                          
                          
                            2
                            ⁢
                            
                              (
                              
                                
                                  κ
                                  j
                                
                                +
                                1
                              
                              )
                            
                          
                        
                      
                      ⁢
                      
                        (
                        
                          
                            
                              
                              
                                q
                                i
                              
                              
                            
                            2
                          
                          +
                          
                            
                              
                                κ
                                j
                              
                              ⁡
                              
                                (
                                
                                  
                                    μ
                                    j
                                    ′
                                  
                                  ⁢
                                  
                                    q
                                    i
                                  
                                
                                )
                              
                            
                            2
                          
                        
                        )
                      
                    
                    ,
                  
                
              
            
            
              
                
                  
                    q
                    i
                  
                  ∉
                  
                    Ω
                    ij
                  
                
                ,
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          

        
        ⁢
        
          
            I
            j
          
          ⁡
          
            (
            
              
                
                  -
                  
                    q
                    i
                  
                
                |
                
                  θ
                  j
                
              
              ,
              γ
            
            )
          
        
      
      =
      
        {
        
          
            
              
                0
                ,
              
            
            
              
                
                  
                    q
                    i
                  
                  ∉
                  
                    Ω
                    ij
                  
                
                ,
              
            
          
          
            
              
                
                  
                    
                      κ
                      j
                    
                    
                      sinh
                      ⁢
                      
                        
                      
                      ⁢
                      
                        κ
                        j
                      
                    
                  
                  ⁢
                  
                    
                      
                        
                          α
                          j
                        
                        ⁢
                        cosh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          α
                          j
                        
                        ⁢
                        sin
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          j
                        
                      
                      -
                      
                        
                          β
                          j
                        
                        ⁢
                        
                          
                        
                        ⁢
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          α
                          j
                        
                        ⁢
                        cos
                        ⁢
                        
                          
                        
                        ⁢
                        
                          β
                          j
                        
                      
                    
                    
                      
                        α
                        j
                        2
                      
                      +
                      
                        β
                        j
                        2
                      
                    
                  
                  ⁢
                  
                    e
                    
                      
                        -
                        
                          
                            γ
                            2
                          
                          
                            2
                            ⁢
                            
                              (
                              
                                
                                  κ
                                  j
                                
                                +
                                1
                              
                              )
                            
                          
                        
                      
                      ⁢
                      
                        (
                        
                          
                            
                              
                              
                                q
                                i
                              
                              
                            
                            2
                          
                          +
                          
                            
                              
                                κ
                                j
                              
                              ⁡
                              
                                (
                                
                                  
                                    μ
                                    j
                                    ′
                                  
                                  ⁢
                                  
                                    q
                                    i
                                  
                                
                                )
                              
                            
                            2
                          
                        
                        )
                      
                    
                  
                
                ,
              
            
            
              
                
                  
                    q
                    i
                  
                  ∉
                  
                    Ω
                    ij
                  
                
                ,
              
            
          
        
      
    
  




The modulus of





  
    A
    ⁡
    
      (
      
        q
        i
      
      )
    
  
  
    A
    ⁡
    
      (
      0
      )
    
  




is given by





  
    
      
        Φ
        i
      
      ⁡
      
        (
        
          q
          i
        
        )
      
    
    
      Φ
      0
    
  
  ,




with equation (E) below:





  
    φ
    ⁡
    
      (
      
        
          q
          i
        
        |
        Θ
      
      )
    
  
  =
  
    
      φ
      0
    
    ⁢
    
      
        
          
            (
            
              
                ∑
                
                  j
                  =
                  1
                
                M
              
              ⁢
              
                
              
              ⁢
              
                
                  f
                  j
                
                ⁢
                
                  
                    R
                    j
                  
                  ⁡
                  
                    (
                    
                      
                        
                          -
                          
                            q
                            i
                          
                        
                        |
                        
                          θ
                          j
                        
                      
                      ,
                      γ
                    
                    )
                  
                
              
            
            )
          
          2
        
        +
        
          
            (
            
              
                ∑
                
                  j
                  =
                  1
                
                M
              
              ⁢
              
                
              
              ⁢
              
                
                  f
                  j
                
                ⁢
                
                  
                    I
                    j
                  
                  ⁡
                  
                    (
                    
                      
                        
                          -
                          
                            q
                            i
                          
                        
                        |
                        
                          θ
                          j
                        
                      
                      ,
                      γ
                    
                    )
                  
                
              
            
            )
          
          2
        
      
    
  




wherein

Φ0 is the given nominal magnetic resonance signal when no MR gradient is applied.

In the equations, R or Re is the real part and I or Im is the imaginary part, and the operator μ′j is the transposition of vector μj.
The above model for the pdf of the random vector X of water molecules displacements is a coherent and realistic representation of the diffusion process in terms of μj, κj, fj and γ. Its characteristic function can also be analytically derived.
By coherent, we intend the following behaviour. Suppose that m=1, i.e, there is one fiber direction at concerned voxel. If κj>>0 or κj>>0, i.e, if the data are highly concentrated around the diffusion direction, it characterizes an anisotropic region and, thus, the expected value of the random vector X of water molecules is equal to γμ. On the contrary, if κj=0 or κ=0, i.e. if the data are uniformly distributed on the sphere, it characterizes an isotropic region and, thus, the expected value of the random vector X of water molecules is null and the MRI signal variation is the same no matter what magnetic field gradient direction is applied.
By a realistic representation of the diffusion process, it is understood that a fundamental difference between the above-mentioned DTI model of document [1] and the invention is that DTI implies symmetry of the probability density function of the random vector X of water molecules, whereas the invention build it asymmetric. In each voxel, the DTI allows the water molecules to diffuse indifferently in a direction and in its opposite direction. On the contrary, the invention allows only the diffusion in one direction and not in its opposite direction which has low probability to occur. Hence, the invention not only covers the real case of one fibre covered by DTI, but also the case of fibre crossings much better than DTI. For example in case of a fibre bifurcating, the water molecules should be allowed only to diffuse along 2 directions and not also along the 2 opposite directions where DTI allows them to move which would involve that the water molecules could go against a “wall”.
Then the values of the m diffusion/fibre direction(s) μj are calculated in each voxel from the N magnetic resonance signals Si acquired for N different magnetic field diffusion gradients qi applied to each voxel. The application of at least 4 m different magnetic field diffusion gradients qi enables to solve the underlying minimization problem to find the m diffusion/fibre direction(s) μj along with concentration parameter κj and the covered mean distance γ in each voxel of the brain.
The above measured signal {Si(qi)}j=1, . . . , m is affected by noise.
The following noise model is minimized




  
    
      
        
          S
          ⁡
          
            (
            
              
                
                  q
                  i
                
                |
                Θ
              
              ,
              σ
            
            )
          
        
        =
        
        ⁢
        
          Rice
          ⁡
          
            (
            
              
                φ
                ⁡
                
                  (
                  
                    
                      q
                      i
                    
                    |
                    Θ
                  
                  )
                
              
              ,
              σ
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            Rice
            ⁡
            
              (
              
                
                  
                    S
                    0
                  
                  ⁢
                  
                    
                      
                        
                          (
                          
                            
                              ∑
                              
                                j
                                =
                                1
                              
                              M
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              
                                f
                                j
                              
                              ⁢
                              
                                
                                  R
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      
                                        -
                                        
                                          q
                                          i
                                        
                                      
                                      |
                                      
                                        θ
                                        j
                                      
                                    
                                    ,
                                    γ
                                  
                                  )
                                
                              
                            
                          
                          )
                        
                        2
                      
                      +
                      
                        
                          (
                          
                            
                              ∑
                              
                                j
                                =
                                1
                              
                              M
                            
                            ⁢
                            
                              
                            
                            ⁢
                            
                              
                                f
                                j
                              
                              ⁢
                              
                                
                                  I
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      
                                        -
                                        
                                          q
                                          i
                                        
                                      
                                      |
                                      
                                        θ
                                        j
                                      
                                    
                                    ,
                                    γ
                                  
                                  )
                                
                              
                            
                          
                          )
                        
                        2
                      
                    
                  
                
                ,
                
                  

                
                ⁢
                
                  
                
                ⁢
                σ
              
              )
            
          
          .
        
      
    
  




One way to approximate the Rician law to normality is due to Salvador, R., Pena, A., Menon, D. K., Carpenter, T. A., Pickard, J. D., and Bullmore. F. T. (2005), Formal characterization and extension of the linearised diffusion tensor model. Human Brain Mapping, 24:144-155. They show that the magnitude DW-MR images' logarithm log(S(q|θ,σ)), which is theoretically log-Rician, can be considered Gaussian with mean parameter log Φ(q|θ,σ)) and standard deviation σ/Φ(g|θ), for any SNR≧3, wherein SNR=Φ(q|θ)/σ. An efficient solution to estimate the diffusion parameters is then to perform a weighted least squares (WLS) optimization, since the variance of the approximating Gaussian distribution varies from one image to another one. These weights w(q) need to be estimated. Once notice that the variance of the approximately Gaussian-distributed random variable Φ(q|θ,σ)log(S(qθ,σ)) is equal to σ2, constant for all the magnitude DW-MR images, and assuming to know good prior estimates {circumflex over (Θ)} of the diffusion parameters θ, we set w(q)=Φ(q|{circumflex over (Θ)})2. Good prior estimate {circumflex over (Θ)} can be found by performing a direct Least Squares optimization as suggested by above-mentioned Salvador et al.'s (2005), i.e,





  
    Θ
    ^
  
  =
  
    arg
    ⁢
    
      
    
    ⁢
    
      
        min
        Θ
      
      ⁢
      
        
          ∑
          
            i
            =
            1
          
          N
        
        ⁢
        
          
        
        ⁢
        
          
            (
            
              
                log
                ⁢
                
                  
                
                ⁢
                
                  S
                  ⁡
                  
                    (
                    
                      q
                      i
                    
                    )
                  
                
              
              -
              
                log
                ⁢
                
                  
                
                ⁢
                
                  φ
                  ⁡
                  
                    (
                    
                      
                        q
                        i
                      
                      |
                      Θ
                    
                    )
                  
                
              
            
            )
          
          2
        
      
    
  




to estimate the diffusion parameters θ.
The problem of the choice of the number of fibers m is a so-called model selection problem. As explained above, if the chosen value of m is higher than the true value, the model tells us to reduce it by estimating some of the κj's around zero. There are also several methods in the literature to overcome such a problem; for instance, the minimization of the Akaike Information Criterion (AIC) is one of them.
The whole method enables to use only 4 m different q to calculate each diffusion direction μj. Of course, more than 4 different qi can be used to calculate each diffusion direction μj. Then 4 m different qi are sufficient to calculate the m diffusion direction μj. Of course, more than 4 m different qi can be used to calculate the m diffusion directions μj.
Thus the number N of magnetic field diffusion gradient qi to apply for the DW-MRI sequence is reduced in view of the state of the art.
If it is assumed that a clinical limit of an MRI examination carried out on a human being is that the time spent by the human being in the MRI device for the DW-MRI sequence should not exceed three minutes, it allows the application of around N=30 magnetic field diffusion gradient qi successively and thus allows the detection of m=4 fibre-crossings in each voxel, letting 13 acquisitions to confer to the estimates robustness to noise.
In case of using κj=κ·fj, there are (3 m+2) parameters to estimate. Then, in an example to estimate m=4 fibres, it gives N−1−m=30−1−(3*4+2)=15 acquisitions for robustness to noise.
An MRI image is reconstructed based on the values of diffusion directions {μj}j=1, . . . , m of the m fibres calculated in each voxel by the above process. In each voxel, the fibre(s) j, for j going from 1 to m, extend along the diffusion directions {μj}j=1, . . . , m. Then the image can be displayed, and/or stored, wherein a multiplicity of voxels for which the above calculation process has been carried out enables to reconstruct multiple fibres in each voxel, and to detect fibre-crossings in case of different μj in the same voxel.
Scalar maps of γ and κ defined as the sum of the κj as described above can also be made.
The method is automatically executed by a computer calculating for each given number m of fibre(s) the m diffusion directions {μj}j=1, . . . , m of the m fibres in each voxel of the image, and possibly also the m parameters {κj}j=1, . . . , m and parameter γ, parameters fj and/or κ defined as the sum of the κj as described above based on the given applied magnetic field diffusion gradients {qi}i=1, . . . , N and the resulting acquired magnetic resonance signals {Si}i=1, . . . , N wherein at least an emitter to apply the magnetic field diffusion gradients {qi}=i=1, . . . , N and a receiver to receive the acquired magnetic resonance signals {Si}i=1, . . . , N are provided.

## Claims

1. A method using an MRI device to extract m diffusion direction(s) μj from N diffusion-weighted magnetic resonance signals Si of any part of a body that contains water molecules, wherein m is a predetermined positive number m of diffusion direction(s) μj to be identified in each voxel of said body part with m≧1,
comprising
exposing by an emitter-receiver the body to a magnetic resonance imaging process, in which N magnetic field diffusion gradients qi, which are non-collinear, are applied by at least one emitter to the entire body part with N being a positive integer ≧2 and the N diffusion-weighted magnetic resonance signals Si are measured by at least one receiver in each voxel of said body part,
wherein using a calculator a first component modelling directional variability of water molecules displacements with diffusion direction μj as mean direction and a second component that models a length variability of the water molecules displacements along said diffusion direction μj are estimated from the N diffusion-weighted magnetic resonance signals Si for each of the in diffusion direction(s) μj,
wherein a sum of the first component and of the second component is calculated in the calculator for each of the m diffusion direction(s) μj,
wherein, for each of the m diffusion direction(s) μj, the diffusion direction μj is extracted in the calculator from the sum being a random vector Xj of the water molecules displacements having a probability density function pXj characterizing the diffusion along the diffusion direction μj and for respectively the N magnetic field gradients qi applied to said body part.

2. Method according to claim 1, characterized in that the first component is a 3-dimensional resealed random unit vector with diffusion direction μj as mean direction.

3. Method according to claim 1 or 2, wherein the first component is equal to γ multiplied by Vj, wherein
Vj is a 3-dimensional random unit vector following a Von Mises & Fisher distribution with mean direction μj and concentration parameter κj around said mean direction μj,
γ is the mean distance covered by the water molecules diffusing along any diffusion direction μj during a given diffusion time τ.

4. Method according to claim 1, wherein the second component Zj is a 3-dimensional random vector whose length varies along the diffusion direction μj.

5. Method according to claim 3, wherein Vj and the second component Zj are independent.

6. Method according to claim 5 wherein the second component is equal to Zj, wherein Zj is a 3-dimensional centered Gaussian random vector whose covariance matrix Σj is equal to




  
    
      ∑
      j
    
    ⁢
    
      
    
    ⁢
    
      =
      
        
          
            γ
            2
          
          
            
              κ
              j
            
            +
            1
          
        
        ⁢
        
          (
          
            
              I
              3
            
            +
            
              
                κ
                j
              
              ⁢
              
                μ
                j
              
              ⁢
              
                μ
                j
                ′
              
            
          
          )
        
      
    
  
  ,



wherein I3 is the 3×3 identity matrix,
and wherein μj is a transposition of a vector corresponding to diffusion direction μj.

7. Method according to claim 3, wherein the N diffusion-weighted magnetic resonance signals Si are the noisy version of a magnitude of DW-MR complex signals A(qi), said magnitude being denoted by Φ(qi) and linked to a diffusion parameter θ through



  
    
      
    
    ⁢
    
      
        φ
        ⁡
        
          (
          
            
              q
              i
            
            |
            Θ
          
          )
        
      
      =
      
        
          φ
          0
        
        ⁢
        
          
            
              
                (
                
                  
                    ∑
                    
                      j
                      =
                      1
                    
                    M
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      f
                      j
                    
                    ⁢
                    
                      
                        R
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            
                              -
                              
                                q
                                i
                              
                            
                            |
                            
                              θ
                              j
                            
                          
                          ,
                          γ
                        
                        )
                      
                    
                  
                
                )
              
              2
            
            +
            
              
                (
                
                  
                    ∑
                    
                      j
                      =
                      1
                    
                    M
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      f
                      j
                    
                    ⁢
                    
                      
                        I
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            
                              -
                              
                                q
                                i
                              
                            
                            |
                            
                              θ
                              j
                            
                          
                          ,
                          γ
                        
                        )
                      
                    
                  
                
                )
              
              2
            
          
        
      
    
  




  
    
      
    
    ⁢
    with
  




  
    
      
        R
        j
      
      ⁡
      
        (
        
          
            
              -
              
                q
                i
              
            
            |
            
              θ
              j
            
          
          ,
          γ
        
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        sin
                        ⁡
                        
                          (
                          
                            γ
                            ⁢
                            
                              
                                
                                  
                                    
                                    
                                      q
                                      i
                                    
                                    
                                  
                                  2
                                
                                -
                                
                                  
                                    (
                                    
                                      
                                        κ
                                        j
                                      
                                      /
                                      γ
                                    
                                    )
                                  
                                  2
                                
                              
                            
                          
                          )
                        
                      
                      
                        γ
                        ⁢
                        
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            -
                            
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  /
                                  γ
                                
                                )
                              
                              2
                            
                          
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          
                            
                            
                              q
                              i
                            
                            
                          
                          2
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        
                          
                            α
                            j
                          
                          ⁢
                          sinh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          cos
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                        +
                        
                          
                            β
                            j
                          
                          ⁢
                          
                            
                          
                          ⁢
                          cosh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          sin
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                      
                      
                        
                          α
                          j
                          2
                        
                        +
                        
                          β
                          j
                          2
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          (
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            +
                            
                              
                                
                                  κ
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      μ
                                      j
                                      ′
                                    
                                    ⁢
                                    
                                      q
                                      i
                                    
                                  
                                  )
                                
                              
                              2
                            
                          
                          )
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
          
          ⁢
          
            

          
          ⁢
          
            
          
          ⁢
          and
          ⁢
          
            

          
          ⁢
          
            
              I
              j
            
            ⁡
            
              (
              
                
                  
                    -
                    
                      q
                      i
                    
                  
                  |
                  
                    θ
                    j
                  
                
                ,
                γ
              
              )
            
          
        
        =
        
          {
          
            
              
                
                  0
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
            
              
                
                  
                    
                      
                        κ
                        j
                      
                      
                        sinh
                        ⁢
                        
                          
                        
                        ⁢
                        
                          κ
                          j
                        
                      
                    
                    ⁢
                    
                      
                        
                          
                            α
                            j
                          
                          ⁢
                          cosh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          sin
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                        -
                        
                          
                            β
                            j
                          
                          ⁢
                          
                            
                          
                          ⁢
                          sinh
                          ⁢
                          
                            
                          
                          ⁢
                          
                            α
                            j
                          
                          ⁢
                          cos
                          ⁢
                          
                            
                          
                          ⁢
                          
                            β
                            j
                          
                        
                      
                      
                        
                          α
                          j
                          2
                        
                        +
                        
                          β
                          j
                          2
                        
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          
                            
                              γ
                              2
                            
                            
                              2
                              ⁢
                              
                                (
                                
                                  
                                    κ
                                    j
                                  
                                  +
                                  1
                                
                                )
                              
                            
                          
                        
                        ⁢
                        
                          (
                          
                            
                              
                                
                                
                                  q
                                  i
                                
                                
                              
                              2
                            
                            +
                            
                              
                                
                                  κ
                                  j
                                
                                ⁡
                                
                                  (
                                  
                                    
                                      μ
                                      j
                                      ′
                                    
                                    ⁢
                                    
                                      q
                                      i
                                    
                                  
                                  )
                                
                              
                              2
                            
                          
                          )
                        
                      
                    
                  
                  ,
                
              
              
                
                  
                    
                      q
                      i
                    
                    ∉
                    
                      Ω
                      ij
                    
                  
                  ,
                
              
            
          
        
      
    
  



wherein qi|Θ means “given the diffusion parameters θ”,


θj={μj,κj}, for any j=1, . . . , M and Θ={θ1, . . . , θM, γ, f1, . . . , fM-1},


wherein Φ0 is a given nominal value of the diffusion-weighted magnetic resonance signal Si when no magnetic field diffusion gradient qi is applied,
wherein fj is the volume fraction of water molecules that diffuse along diffusion direction μj, with 0≦fj≦1 and




  
    
      
        ∑
        
          j
          =
          1
        
        m
      
      ⁢
      
        
      
      ⁢
      
        f
        j
      
    
    =
    1
  
  ,




wherein μ′j is a transposition of a vector corresponding to diffusion direction μj,



  
    
      Ω
      ij
    
    =
    
      {
      
        
          
            
              -
              
                q
                i
              
            
            ∈
            
                            3
            
          
          :
          
            
              
              
                -
                
                  q
                  i
                
              
              
            
            ≥
            
              
                κ
                j
              
              /
              γ
            
          
        
        ,
        
          

        
        ⁢
        
          
            and
            ⁢
            
              

            
            -
            
              q
              i
            
          
          ∈
          
            
              (
              
                span
                ⁡
                
                  (
                  
                    μ
                    j
                  
                  )
                
              
              )
            
            ⊥
          
        
      
      ]
    
  




  
    
      α
      ij
    
    =
    
      
        
          
            Re
            ⁢
            
              
            
            ⁢
            
              z
              ij
            
          
          +
          
            
            
              z
              ij
            
            
          
        
        2
      
    
  




  
    
      β
      ij
    
    =
    
      
        Im
        ⁢
        
          
        
        ⁢
        
          z
          ij
        
      
      
        
          2
          ⁢
          
            (
            
              
                Re
                ⁢
                
                  
                
                ⁢
                
                  z
                  ij
                
              
              +
              
                
                
                  z
                  ij
                
                
              
            
            )
          
        
      
    
  




  
    
      z
      ij
    
    =
    
      
        κ
        j
        2
      
      -
      
        
          γ
          2
        
        ⁢
        
          
            
            
              q
              i
            
            
          
          2
        
      
      -
      
        2
        ⁢
        i
        ⁢
        
          
        
        ⁢
        
          κ
          j
        
        ⁢
        γ
        ⁢
        
          
        
        ⁢
        
          μ
          j
          ′
        
        ⁢
        
          q
          i
        
      
    
  



(span(μj))⊥ is a space orthogonal to (span(μj)), where (span(μj)) is a vector sub-space of generated by unit vector μj.

8. Method according to claim 7, wherein the concentration parameters κj are linked to the volume fractions fj by the following relationship


κj=κ·fj 

where κ is a measure of anisotropy in the voxel, with





  κ
  =
  
    
      ∑
      
        j
        =
        1
      
      m
    
    ⁢
    
      
    
    ⁢
    
      
        κ
        j
      
      .

9. Method according to claim 1 or claim 2, wherein a multiplicity of said voxels is provided, a magnetic resonance image of said voxels is reconstructed by the calculator based at least on the values of the diffusion directions μj calculated for each voxel.

10. An MRI device to extract m diffusion direction(s) μj from N diffusion-weighted magnetic resonance signals Sj of any part of a body that contains water molecules, wherein m is a predetermined positive number of diffusion direction(s) μj to be identified in each voxel of said body part with m≧1,
comprising
an emitter-receiver for exposing the body to a magnetic resonance imaging process, the emitter-receiver comprising at least one emitter for applying N magnetic field diffusion gradients qi, which are non-collinear, to the entire body part with N being a positive integer ≧2 and at least one receiver for measuring the N diffusion-weighted magnetic resonance signals Si in each voxel of said body part,
a calculator having means to estimate, for each of the m diffusion direction(s) μj, a first component modelling directional variability of water molecules displacements with diffusion direction μj as mean direction and a second component that models a length variability of the water molecules displacements along said diffusion direction μj from the N diffusion-weighted magnetic resonance signals Si,
wherein the calculator has means to calculate a sum of the first component and of the second component for each of the m diffusion direction(s) μj,
wherein the calculator has means to extract, for each of the in diffusion direction(s) μj, the diffusion direction μj from the sum being a random vector Xj of the water molecules displacements having a probability density function pXj characterizing the diffusion along the diffusion direction μj and for respectively the N magnetic field diffusion gradients qi applied to said body part.