AthenaEPI / dmipy

The open source toolbox for reproducible diffusion MRI-based microstructure estimation
MIT License
96 stars 30 forks source link

A little confuse about the explaination for noddi parameter's maps #62

Closed foxet closed 5 years ago

foxet commented 5 years ago

The explanation for parameter's maps followed is from example_noddiwatson.ipynb

"Above we see the familiar maps. Remember that partial_volume_0/1/2 represent fractions of models in the order they were given to the MultiCompartmentMicrostructureModel object. In this case, 0 represents the fraction of Ball (CSF), and the 1 represents the Stick (intra-axonal) and 2 represents the zeppelin (extra-axonal).

In the bottom-right ODI map you can see that the concentration is highest in the corpus callosum and in coherent white matter bundles. Finally, the πœ†βˆ₯ and πœ†iso maps are fixed to the values we gave them."

_ 图片

I think the partial_volume_0 is from G1ball, represents fraction of CSF

I have two questions: **1. what's the parameters for partial_volume_1 and partial_volume_2 in the explaination

Does it refers to 'partial_volume_1' and 'SD1WatsonDistributed_1_partial_volume_0'? or refers to: 'SD1WatsonDistributed_1_partial_volume_0' and 'SD1WatsonDistributed_1_SD1Watson_1_mu'

  1. Is the map for Stick (intra-axonal) parameter equal to NDI (neurite density index ) map?**

Here's the watson parameters: NODDI_mod.parameter_names

['G1Ball_1_lambda_iso',
'SD1WatsonDistributed_1_SD1Watson_1_mu', 'SD1WatsonDistributed_1_partial_volume_0', 'SD1WatsonDistributed_1_SD1Watson_1_odi', 'partial_volume_0', 'partial_volume_1']

It would be appreciated that if you could help, : )

rutgerfick commented 5 years ago

Hey!

I updated the example where i explicitly compute the stick and zeppelin maps from the parameters: https://github.com/AthenaEPI/dmipy/blob/master/examples/example_noddi_watson.ipynb

Answers to your question

  1. partial_volume_0 is for the ball, and partial_volume_1 is for the Watson Bundle. SD1WatsonDistributed_1_partial_volume_0 is for the Stick INSIDE the watson bundle, and (1-SD1WatsonDistributed_1_partial_volume_0) is for the zeppelin (so the complement). To the get signal contribution of the Stick and Zeppelin to the total signal, you have to multiply them by the fraction of the watson bundle (partial_volume_1).

  2. So, the stick in the model represents the intra-axonal/neurite fraction, and the zeppelin the extra-axonal/neurite fraction. The total neurite / white matter fraction is thus the sum of these two, which is just partial_volume_1. So, it depends what you want, and what you define NDI to be :-)

foxet commented 5 years ago

very clear now, thankyou rutgerfick