Chris<br><br>It's been a while since I've worked with the 0.2 tracking code, but from a quick glance it seems that the behaviour is the following:<br>* If the mask is exited, the first point outside the mask is <i>not</i> included in the streamline.<br>
* If the relevant diffusion metric (FA / FOD amplitude) is below threshold, the first point at which this criterion fails <i>will be</i> included in the streamline.<br>* If the curvature threshold is exceeded, the point at the end of the step that has exceeded this criterion is <i>not</i> included (this applies for DT_STREAM and SD_STREAM only).<br>
* The behaviour on entering an exclusion ROI is irrelevant as the streamline is rejected regardless.<br><br>When we eventually get around to releasing the 0.3 code, this behaviour is made much more explicit in the code design, and can therefore be modified if you prefer specific behaviour for specific termination criteria.<br>
<br>Regarding termination of a streamline if it re-enters the seed region, I think the answer depends on the behaviour you are seeking; do you want the streamlines to be simply terminated if they re-enter the seed region, or would you like them to be rejected as artefact? If the latter, a modification of the filter_tracks command may be more digestible; you could generate a vector specifying whether each point along the streamline is within the seed mask, and reject any streamlines for which there are two discrete streamline segments within the seed mask. Implementing the former may be a little trickier, as you would need to deal with confounds such as multiple seed regions and bi-directional tracking, but it can certainly be done; let me know if you try but get stuck.<br>
<br>I believe there was a streamtrack modification distributed on the discussion board some time ago that wrote each streamline seed to a text file, though I'm having trouble finding it. Unfortunately it's not currently possible to embed the data for each streamline seed within the .tck file, so this needs to be handled externally. Of course, if you are tracking unidirectionally, the seed point will always be the first point on the streamline.<br>
<br>Best regards<br>Rob<br clear="all"><br>--<br><br>Robert Smith<br>PhD Candidate <br><br>The Florey Institute of Neuroscience and Mental Health<br>Melbourne Brain Centre - Austin Campus<br>245 Burgundy Street<br>Heidelberg Vic 3084<br>
Ph: +61 3 9035 7128<br>Fax: +61 3 9035 7301<br><a href="http://www.florey.edu.au" target="_blank">www.florey.edu.au</a><br><span style="font-size:9pt;font-family:"Times New Roman","serif";color:red"></span><br>
<br><br><div class="gmail_quote">On Tue, Nov 6, 2012 at 2:43 AM, Parker, Christopher <span dir="ltr"><<a href="mailto:christopher.parker.10@ucl.ac.uk" target="_blank">christopher.parker.10@ucl.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi Donald,
<div><br>
</div>
<div>I have a few questions about streamtrack. When streamlines terminate (i.e. they enter an exclusion ROI, leave the mask ROI, encounter low FA or high curvature), are the coordinates at the current point (i.e. the coordinate that fulfils the termination
criteria) recorded, or are the streamline coordinates terminated just before this point (the step before termination). I'm guessing for high curvature it is the latter but for the others it is the former.</div>
<div><br>
</div>
<div>I was also wondering if there is any way to specify a termination mask that is also a seed mask using streamtrack, so that streamlines that re-enter the mask are terminated, but are free to initially propagate within the mask. It is not possible with the
command line tools as far as I can see. This would be really useful- maybe I can modify the streamtrack command to make this possible? </div>
<div><br>
</div>
<div>Another thing that would useful is to modify the streamtrack command so that it records the streamline starting coordinates.</div>
<div><br>
</div>
<div>Kind Regards,</div>
<div>Chris</div>
</div>
</div>
<br>_______________________________________________<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>
<br></blockquote></div><br>