There’s nothing wrong with the top one: those are SH representations of the DWI signal (not FODs), which has large magnitudes perpendicular to the main fibre orientation(s).
As the bottom one is in fact labelled “FOD.jpg”, it’s the bottom one that has a severe mistake in it: note the orientation of those (supposedly intended to be FODs) is wrong in e.g. the corpus callosum and the corticospinal tract, indicating a problem with the gradient table (a combination of flips of axes and permutations of columns is typically the cause of such a problem).
It’s not wrong, as in: that is definitely the correct output of amp2sh. The question though is what you’re ultimately after (i.e. how you see yourself continuing from this output onwards)… Were you genuinely after an SH representation of the DWI signal? Or are you after FODs? (because the signal is something very different from FODs)
If you’re confused about this, I’d definitely encourage you to read up a bit on the literature (e.g. read Donald’s spherical deconvolution and CSD papers), before proceeding…
What I am after is interpolating the data we have. We started form DSI and used Dsi Studio to interpolate them.
Now we would like to select just the points we want and where we want them. To do this, we passed to SH and then we will choose the points and go back.
Sounds alright then, as long as you take into account that it’s an SH fit of course. Definitely make sure you’re providing amp2sh with a dense enough sampling, and potentially choose the lmax wisely.
I also have some question regarding amp2sh after three-tissue CSD (SS3T)
I want to calculate the the SH coefficient of FOD as one of the local diffusion feature for thalamus clustering
My questions are:
After SS3T, Does the white matter fod (wm_fod.mif) could be directly used as input for amp2sh to calculate the SH coeffcient, however the thalamus tissue consist of white matter and grey matter, so I am not sure if I should use the wm_fod.mif only, only need to combine with gm_fod.mif for the SH coefficient. My knowledge would be the SH coefficient of FOD image is the coefficient of fiber orientation distribution, so the fiber here means ‘white matter fiber’? And what should I do to extract the reliable SH coefficients as a clustering input
I have one dataset only with 34 diffusion MRI volumes, according to the SH coefficient store rules, they should calculate the SH coefficient with lmax order 6 rather the default lmax order 8 due to the limited volume. so another dataset with 54 diffusion MRI volumes, which should calculate the SH coefficient with lmax order 8, if i only want to retrieve the SH coefficient value with spherical order 6 from maximum order 8( to make it comparable across different dataset), can i just directly select the first 28 volumes from the FOD, should i need to re-run CSD with set the lmax order to 6
Any insightful suggestion would be greatly appreciated!