Hi Luis,
there are many issues with all the approaches that you are considering. I do not think you can do a meaningless comparison, for the following reasons:
- It is a bit strange to eliminate data from a good dataset (64 directions) in order to match the poorer quality of a 16 direction dataset. I would try to keep the best quality data and work with those.
- The number of directions alone is not determining the quality of the data. The direction schemes are (or should be) optimized, so that they are uniformly distributed on the sphere. A dataset of 16 uniform directions and another of 16 random directions are not comparable, the latter is likely to give you biased estimates. If you insist eliminating directions from the 64 direction dataset, this should be done in a systematic way that preserves the optimality of the scheme.
- Averaging directions: Are your directions unique? If yes, then you cannot average, you need to have repeats of exactly the same directions to average them. Averaging will increase the SNR for the datasets that you can apply this averaging compared to the datasets that you cannot. This means better resolving power for the former.
- Are the datasets of different resolution? What about the sequence itself, is it the same? Eddy-current compensation schemes have been employed for all acquisitions? Differences in all these factors can cause significant differences in the quality of the data and make them even less comparable.
- Even if we hypothetically ignore all the above issues, and assume that you can reduce all datasets down to 16 directions, they will not be good for probabilistic tractography. Depending on the SNR and the b-value, you will hardly be able to resolve crossings and also the path distribution is likely to be very wide. If however you want to do tractography using such a low number of directions, I would suggest keeping 1 fibre compartment per voxel (i.e. no crossings). Some "easy" tracts (in terms of geometry) should be reconstructed, so it depends on your application and the tract you are after.
Hope this is clear,
Stam
On 29 Mar 2011, at 20:52, Luis Morís wrote:
> Hello everyone,
>
> At the moment I'm working on tractography with DTI images. I'm trying
> to make an study on images which has been done in different machines,
> (I know this is highly discouraged) and some of those images has
> different sequences from 16 to 64 gradient directions. It's not the
> ideal scenario.
>
> My doubt is if there is anyway of reducing the 64 directions images to
> 16 in order to make them comparable. I have thought of simply erasing
> gradients until I have 16, but that could lead to some unbalance. Also
> averaging some images etc...
>
> Is there a known method for acomplishing this? Does averaging DTI
> images and gradients makes any sense?
>
> Cheers and thanks for your help,
>
> Luis.
>
|