Hi Marina,<div><br></div><div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">How the command streamtrack choose the seed point?<br></blockquote>
<div><br></div><div>Seed points are generated at random by uniform sampling from the mask. The mask is defined as any location in 3D space where the image intensity after interpolation is greater than 0.5. You can disable the mask interpolation step with the -nomaskinterp option, so that seeds are only placed strictly within the voxels of interest. Unfortunately there is no option currently to ensure regular seed placement. </div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Morever, while the command is running, I notice that there are generated streamlines but only some of them are selected, and the number of selected streamlines corresponds to the value given at the -number parameter.<br>
Which is the criterium by which streamlines are selected?<br></blockquote><div><br></div><div>Streamlines are selected if they satisfy all constraints. By default, there is a minimum length constraint, which I guess is the only reason you would currently get different numbers of generated vs. selected tracks. You can change that using the -minlength option. In general, additional constraints can be placed using the -include & -exclude options. </div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Is there an option to set the number of generated streamlines instead of the number of selected ones?<br>
I've tried to run the command without any options but the numbers of generated tracks and selected ones do not match.<br></blockquote><div><br></div><div>There is no explicit option, but same result can be achieved very easily by setting -maxnum and -number to the same value (by default, maxnum is 100xnumber). </div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
My aim is to have a certain number N of tracks beginning from each seed point (voxel with intensity equal to 1) of the mask.<br></blockquote><div><br></div><div>This is a common request, but it's still not possible to do very easily. If the number of streamlines generated is large, you can certainly set the number of tracks to num_voxels x num_tracks_per_voxels, and get a very good approximation to uniform sampling of your ROI (you can find out the number of voxels in the ROI using the mrstats command). If you need to save the tracking results for each voxel separately, then you'll need to break up your ROI into individual voxel ROIs, and launch one streamtrack command per voxel. This can easily be done with a bit of scripting, but there is no easy way within MRtrix to break up the ROI into individual voxel ROIs. Maybe this is something that I could provide, since it would be trivial to code up. I'll have a think about that...</div>
<div> </div><div>Hope that helps,</div><div><br></div><div>Donald.</div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Thank you for your helpfulness,<br>
Marina_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org">Mrtrix-discussion@www.nitrc.org</a><br>
<a href="http://www.nitrc.org/mailman/listinfo/mrtrix-discussion" target="_blank">http://www.nitrc.org/mailman/listinfo/mrtrix-discussion</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Jacques-Donald Tournier (PhD)<br>Brain Research Institute, Melbourne, Australia<br>Tel: +61 (0)3 9035 7033<br>
</div>