[Mrtrix-discussion] Number of tracks

[Kerstin Pannek]

Hi Luis

Are you doing tensor tractography? You should be able to use constrained spherical deconvolution with your data. I have some 30 direction data with b=1000, and use csdeconv with -lmax 4.

Medinria uses deterministic tractography (correct me if I'm wrong). With deterministic tractography, there is only one possible streamline for every seed point. With probabilistic tractography, there is a large number of possible streamlines for every seed point. Therefore, you need to calculate a much larger number of streamlines using probabilistic tractography.

It very much depends on your application if -number 100000 is enough. The easiest way to find out would be to do your processing all the way for one participant; then repeat the processing, say, 10 times, starting from streamtrack. Have a look at the metric you are interested in, and how variable it is across those 10 trials. If your results are fairly tight, 100000 is enough; if there is huge variation, you will need to increase the number of streamlines.

Btw, FA initcutoff 0.6 and cutoff 0.3 seems fairly high to me. The MRtrix defaults are 0.2 (init) and 0.1; this works very well in my experience. I'm not sure what the cutoffs are for other software.

Cheers
Kerstin




On 12/05/2011, at 12:39 AM, Luis Morís wrote:

> Hi again,
> 
> Thanks for your quick answer René.
> 
> My data is 31 directions  with b=1000. The point is that no matter how
> much fibers I put in the -number parameter, it always finds enough
> fibers even a couple millions fibres. Previously I was working with
> Medinria and I had around 30000 fibres because seeding was dependant
> on resolution, one seed per voxel.
> 
> My parameters are FA cutoff 0.3 and initcutoff 0.6 min and max length
> 10 and 200 mm respectively, at the moment I'm doing a whole brain
> seeding. My problem is how to set this number of tracks. If I use a
> -number of let's say a hundred thousand, will I have a nice
> distribution of tracks that will be representantive of the brain's
> track population?
> 
> 
> Thanks in advance,
> 
> Luis.
> 
> On Wed, May 11, 2011 at 4:18 PM, René Besseling
> <r.m.h.besseling at gmail.com> wrote:
>> Hi Luis,
>> 
>> Just type streamtrack -h; the tracks that are selected are those of
>> sufficient length (-minlength option) and sufficiently high FA or FOD
>> amplitude (-cutoff value). Note that the -maxnum option sets the
>> maximum allowed number of trials. If your data is of low quality, the
>> number of tracks generated can be lower than what you set by -number
>> if not enough valid tracks can be constructed within -maxnum trials.
>> Furthermore, tractography start points (seed points) are selected at
>> random within the ROI you specify using -seed, so if you use the same
>> streamtrack command twice, the result will be different (due to
>> different random seed points). Information about a streamtrack result
>> can be obtained by typing track_info your_file_name.tck
>> 
>> Best regards,
>> 
>> René
>> 
>> On Wed, May 11, 2011 at 3:54 PM, Luis Morís
>> <luis.moris.fernandez at gmail.com> wrote:
>>> Hello everyone,
>>> 
>>> I write because I have started to use mrtrix, and I've found a doubt.
>>> To be able to do tractography you have to use the parameter -number
>>> that will set the exact amount of tracks that will be calculated.
>>> 
>>> My doubt is how does the program decide which tracks to include or
>>> not. If I set it to 1000 and I run it again with 1000000. What
>>> difference do I find between the two set of tracks? Is it purely
>>> random, depends on the length, the FA...
>>> 
>>> Thanks a lot in advance,
>>> 
>>> Luis.
>>> _______________________________________________
>>> Mrtrix-discussion mailing list
>>> Mrtrix-discussion at www.nitrc.org
>>> http://www.nitrc.org/mailman/listinfo/mrtrix-discussion
>>> 
>> 
> _______________________________________________
> Mrtrix-discussion mailing list
> Mrtrix-discussion at www.nitrc.org
> http://www.nitrc.org/mailman/listinfo/mrtrix-discussion