<div dir="ltr">Hi David.<br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>
<br></div><div>Luis, I'm guessing when you say tract-specific, you are averaging the AFD within a tract-of-interest? Since AFD is proportional to the intra-axonal volume, ideally you need to be summing the AFD to account for differences in fibre bundle <i>width</i> (as the modulation does in voxel-wise analysis). <br>
</div></div></blockquote><div>This is where I still have questions and value your input. What I want in the end is a volume in which each voxel has the AFD for a specific bundle. That is, I could have as many of these AFD volumes as dissected bundles. Let's imagine just one. If I end up with such an AFD volume, I imagine I could then simply average all the voxels that contain the bundle and get a tract measurement. I could also sample these AFDs along the tracks and get a spatially varying AFD. <br>
<br></div><div>My question has to do with the summing of AFDs. Do you mean that if 50 streamlines traverse a given voxel, then the AFD for that voxel should be the sum of 50 AFDs? In theory, these 50 FOD lobes should be very similar (give or take some interpolation). I still do not understand why summing 50 AFDs would be better than the mean of them. Perhaps I did not understand correctly, and perhaps I have not expressed myself properly in the sense that I would like a voxel-wise AFD for a given tract.<br>
<br>Helen: Once I figure out these details, I will be happy to post my script.<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">
<br><div><br>
</div><div>The next major release of MRtrix will contain SIFT and some AFD tools. However, FOD registration probably won't be included until a few months after the initial release.</div></div></blockquote><div dir="ltr">
Music to my ears!<br><br></div><div><br></div><div dir="ltr"><br></div><div>Luis<br></div><div dir="ltr"><div><br></div><div><br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div class="gmail_extra"><div><div><br><br><div class="gmail_quote">On 27 March 2014 08:06, Luis Concha <span dir="ltr"><<a href="mailto:lconcha@unam.mx" target="_blank">lconcha@unam.mx</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Dorian and fellow mrtrixers. <br><br>Yes, it is very much like HMOA, which is actuallly just like AFD (*,**) . There is, however, one difference: In the approach I mentioned, the tracts are first computed using the tractography algorithms supplied by streamtrack, and in the paper by Dell'Acqua, HMOA is used to drive the tractography algorithm itself, as well as to select tract-specific HMOA a posteriori (Figure 9 in Dell'Acqua). <br>
<br></div><div>My main concern is that I use find_SH_peaks in order to simplify the identification of the corresponding (parallel) FOD peak to the tract segment. Given that find_SH_peaks output can be affected by the number of directions used as seeds for peak finding, is there any min number that should be used? (I used 120).<br>
<br><br></div><div>I would be happy to share my scripts to obtain AFD/HMOA, if anyone is interested.<br></div><div>
<br><br></div>(*) Marco Catani, during the ISMRM Diffusion Workshop last year, said that it was one of the rewiewers who asked them to change AFD to HMOA. <br>(**) One key point in Raffelt's paper on AFD is the implementation of a voxel-based framework for the analysis of this measurement for group analyses, which I am hoping will be available in mrtrix soon.<span><font color="#888888"><br>
<br>
<div class="gmail_extra">Luis<br></div></font></span><div><div><div class="gmail_extra"><br><br clear="all"><div><div dir="ltr"><br></div></div>
<br><br><div class="gmail_quote">On Wed, Mar 26, 2014 at 2:17 PM, Dorian P. <span dir="ltr"><<a href="mailto:alb.net@gmail.com" target="_blank">alb.net@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">Hello everybody, hello Dr. Concha.<div><br></div><div>As we are on the topic, it seems to me the idea of Dr. Concha is similar to the HMOA measure of Dell'Acqua et al. <a href="http://www.ncbi.nlm.nih.gov/pubmed/22488973" target="_blank">http://www.ncbi.nlm.nih.gov/pubmed/22488973</a></div>
<div><br></div><div>Am I correct thinking they are similar, if not identical measures?</div><div><br></div><div>I am looking for ways to obtain this measure in mrtrix. Anyone can give some hints or share the steps to do this? Will mrtrix plan to output the amplitude of FODs for tracts in next releases?</div>
<div><br></div><div><br></div><div>Thank you</div><div>Dorian</div><div>Thomas Jefferson University</div><div><div><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Wed, Mar 26, 2014 at 3:58 PM, Luis Concha <span dir="ltr"><<a href="mailto:lconcha@unam.mx" target="_blank">lconcha@unam.mx</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><br></div>Hi Thijs.<br><br></div>It is great to know that the new release of mrtrix will have tools for AFD. In the meantime, I have played around with some scripts to obtain such a measurement in a tract-specific way and I would appreciate your input and that of others in this list. The scripts do this:<br>
<br></div>1. Perform tractography and select a set of streamlines for analysis.<br></div>2. For each streamline segment (i.e., the straight line between two adjacent steps), find the FOD lobe that is most parallel (dot product between streamline segment and FOD peaks).<br>
</div>3. Once the FOD lobe has been identified, extract its amplitude (tract-specific AFD). <br></div><br></div>I am assuming that if DWI signal is normalized (as you mentioned), then FOD amplitudes should be comparable between subjects. Also, since there is no spatial normalization of images, there is no need to modulate the AFD to compensate for the spatial transformation.<br>
<br></div>If you find any flaws with this, please let me know. Oh, and if you could give us a teaser of what new features we will see in the new release of mrtrix, that would be fantastic!<br><div><br><br><div><br></div>
</div>
</div><div class="gmail_extra"><br clear="all"><div><div dir="ltr">Dr. Luis Concha<br>Instituto de Neurobiología<br>Laboratorio C-13<br>UNAM, Campus Juriquilla<br>Boulervard Juriquilla 3001<br>Juriquilla, Querétaro.<br>C.P. 76230<br>
México<br>Tel <a href="tel:%28442%29%202%2038%2010%2054" value="+14422381054" target="_blank">(442) 2 38 10 54</a><br>Fax <a href="tel:%28442%29%202%2038%2010%2046" value="+14422381046" target="_blank">(442) 2 38 10 46</a><br>
<a href="http://personal.inb.unam.mx/lconcha/" target="_blank">http://personal.inb.unam.mx/lconcha/</a><br></div></div><div><div>
<br><br><div class="gmail_quote">On Wed, Mar 26, 2014 at 1:47 PM, Thijs Dhollander <span dir="ltr"><<a href="mailto:thijs.dhollander@gmail.com" target="_blank">thijs.dhollander@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr"><div><div><div><div><div>Hi Helen,<br><br></div>Opinions on FA/ADC quantification can be wildly varying... :-) There's no problem calculating and reporting their values, but interpretation is a big issue: they can be quite sensitive in many cases, but they are always far from specific. To avoid starting a long rant on how we should abandon these measures as quickly as possible, I'd rather refer to this nice article: <a href="http://www.sciencedirect.com/science/article/pii/S1053811912007306" target="_blank">http://www.sciencedirect.com/science/article/pii/S1053811912007306</a> ; it contains some nice and clear passages explaining some of the relevant issues to be wary of, as well as further references that are worth a good read.<br>
</div>So the question remains: what's better then? There's quite some alternatives, but I personally particularly like apparent fibre density (AFD). A nice and in-depth article on it can be found here: <a href="http://www.sciencedirect.com/science/article/pii/S1053811911012092" target="_blank">http://www.sciencedirect.com/science/article/pii/S1053811911012092</a> . Last year's ISMRM had an abstract on some further improvements on analysing AFD: <a href="http://cds.ismrm.org/protected/13MPresentations/abstracts/0841.pdf" target="_blank">http://cds.ismrm.org/protected/13MPresentations/abstracts/0841.pdf</a> (if you or anyone in your lab attended, they can access it; the accompanying oral presentation can be found at <a href="http://cds.ismrm.org/protected/13MPresentations/0841/" target="_blank">http://cds.ismrm.org/protected/13MPresentations/0841/</a> ). The next big version of MRtrix (expected to be released somewhere around this year's ISMRM in May... ;-)) will contain much more tools that provide you with the means to perform such analyses.<br>
</div>In the mean time, you can already easily obtain an (orientationally averaged) AFD map by simply calculating the average of your (normalized) DWI images from a single shell acquisition. The same contrast (up to a fixed constant factor) can also be obtained by opening up the result from an SH fit of the signal (or even a CSD outcome) in mrview, and looking at the first coefficient (i.e. the 0th order).<br>
<br></div><div>Hope this already helps somewhat!<br></div><div><br></div>Cheers,<br>Thijs<br><br></div>PS: now that I've mentioned the upcoming ISMRM; I'd like to encourage everyone on this list that still has many questions (or not) to surely attend and look for each other -- discussions on these new and better ways of quantifying are always much more lively in person :-)<br>
<div><div><div><div><div><div><div class="gmail_extra"><div><br clear="all"><div><div dir="ltr"><span style="color:rgb(68,68,68)"><b><font size="4">Thijs Dhollander</font><br><font><a href="mailto:thijs.dhollander@gmail.com" target="_blank">thijs.dhollander@gmail.com</a><br>
Tel. <a href="tel:%2B32%20475%2036%2044%2027" value="+32475364427" target="_blank">+32 475 36 44 27</a></font></b><br><font size="1">Medical Image Computing (MIC), ESAT-PSI, Department of Electrical Engineering, KU Leuven</font></span></div>
</div>
<br><br></div><div><div><div class="gmail_quote">On Wed, Mar 26, 2014 at 4:52 PM, Helen Carlson <span dir="ltr"><<a href="mailto:Helen.Carlson@albertahealthservices.ca" target="_blank">Helen.Carlson@albertahealthservices.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div link="blue" vlink="purple" lang="EN-CA"><div><p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)">Hello Thijs<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)">Thank-you for your reply to my question. I would like to quantify changes in white matter over time, specifically pre/post treatment. Is there a better way to do this than using the mean FA and ADC? Or is this measurement appropriate? <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)">Thanks again<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)">Helen.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:"Calibri","sans-serif";color:rgb(31,73,125)"><u></u> <u></u></span></p><div style="border-width:1pt medium medium;border-style:solid none none;border-color:rgb(181,196,223) -moz-use-text-color -moz-use-text-color;padding:3pt 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:10pt;font-family:"Tahoma","sans-serif"" lang="EN-US">From:</span></b><span style="font-size:10pt;font-family:"Tahoma","sans-serif"" lang="EN-US"> <a href="mailto:mrtrix-discussion-bounces@www.nitrc.org" target="_blank">mrtrix-discussion-bounces@www.nitrc.org</a> [mailto:<a href="mailto:mrtrix-discussion-bounces@www.nitrc.org" target="_blank">mrtrix-discussion-bounces@www.nitrc.org</a>] <b>On Behalf Of </b>Thijs Dhollander<br>
<b>Sent:</b> March 23, 2014 9:47<br><b>To:</b> <a href="mailto:mrtrix-discussion@www.nitrc.org" target="_blank">mrtrix-discussion@www.nitrc.org</a><br><b>Subject:</b> Re: [Mrtrix-discussion] FA and ADC values for a tract<u></u><u></u></span></p>
</div><p class="MsoNormal"><u></u> <u></u></p><div><div><div><p class="MsoNormal" style="margin-bottom:12pt">Hi Helen,<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12pt">What you're seeing is (probably) correct: FAs larger than 1 and ADCs below 0 can appear due to the tensor estimation in MRtrix being a simple (unconstrained) linear least squares one. The tensors are not explicitely constrained to be positive definite, and thus negative eigenvalues might appear. When one or more eigenvalues of a given tensor are negative, the FA can go beyond 1, while the ADC can go below 0.<u></u><u></u></p>
</div><p class="MsoNormal">Cheers,<br>Thijs<u></u><u></u></p><div><div><div><div><p class="MsoNormal"><br clear="all"><u></u><u></u></p><div><div><p class="MsoNormal"><b><span style="font-size:13.5pt;color:rgb(68,68,68)">Thijs Dhollander</span><span style="color:rgb(68,68,68)"><br>
<a href="mailto:thijs.dhollander@gmail.com" target="_blank">thijs.dhollander@gmail.com</a><br>Tel. <a href="tel:%2B32%20475%2036%2044%2027" value="+32475364427" target="_blank">+32 475 36 44 27</a></span></b><span style="color:rgb(68,68,68)"><br>
</span><span style="font-size:7.5pt;color:rgb(68,68,68)">Medical Image Computing (MIC), ESAT-PSI, Department of Electrical Engineering, KU Leuven</span><u></u><u></u></p>
</div></div><p class="MsoNormal" style="margin-bottom:12pt"><u></u> <u></u></p><div><p class="MsoNormal">On Fri, Mar 21, 2014 at 10:10 PM, Helen Carlson <<a href="mailto:Helen.Carlson@albertahealthservices.ca" target="_blank">Helen.Carlson@albertahealthservices.ca</a>> wrote:<u></u><u></u></p>
<div><div><p class="MsoNormal">Hello all,<u></u><u></u></p><p class="MsoNormal">I would like to determine the average FA and ADC values for a given track (in this case, the arcuate). So I followed previous directions posted here and have used the following commands on my arcuate which successfully creates a probability mask for FA and ADC:<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">tracks2prob –template fa.mif Arcuate.tck FA_mask.mif<u></u><u></u></p><p class="MsoNormal">tracks2prob –template adc.mif Arcuate.tck ADC_mask.mif<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">I then used mrstats to get the average, max, min etc of the tracks by using these commands:<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">
mrstats fa.mif –mask FA_mask.mif<u></u><u></u></p><p class="MsoNormal">mrstats adc.mif –mask ADC_mask.mif<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">Which also seemed to work except that I was expecting my FA to range from 0.1 to 1.0 (approximately). Instead, the FA values have a min of 0 and a max of 1.13. In addition, my ADC values have a negative minimum (-0.0024). It seems these values are out of range and do not make sense but I do not know what I am doing wrong.<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal">Thanks in advance for any help<u></u><u></u></p><p class="MsoNormal">Helen<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"><b><span style="color:rgb(54,95,145)">Helen Carlson, Ph.D</span></b><u></u><u></u></p><p class="MsoNormal"><span style="font-size:10pt;color:rgb(54,95,145)">Neuroimaging Research Associate</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:10pt;color:rgb(54,95,145)">Calgary Pediatric Stroke Program & Neuropsychology</span><u></u><u></u></p><p class="MsoNormal"><span style="font-size:10pt;color:rgb(54,95,145)">Alberta Children's Hospital</span><u></u><u></u></p>
</div><p class="MsoNormal"><u></u> <u></u></p><div class="MsoNormal" style="text-align:center" align="center"><hr align="center" size="2" width="100%"></div><p class="MsoNormal"><span style="font-size:7.5pt;font-family:"Arial","sans-serif";color:gray">This message and any attached documents are only for the use of the intended recipient(s), are confidential and may contain privileged information. Any unauthorized review, use, retransmission, or other disclosure is strictly prohibited. If you have received this message in error, please notify the sender immediately, and then delete the original message. Thank you.</span><u></u><u></u></p>
</div><p class="MsoNormal" style="margin-bottom:12pt"><br>_______________________________________________<br>Mrtrix-discussion mailing list<br><a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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><u></u><u></u></p></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div>
</div></div></div></div><br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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></div></div></div></div></div></div></div></div></div></div>
<br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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></div></div></div>
<br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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></div></div></div></div>
<br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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></div></div></div></div>
<br>_______________________________________________<br>
Mrtrix-discussion mailing list<br>
<a href="mailto:Mrtrix-discussion@www.nitrc.org" target="_blank">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><br clear="all"><div><br></div></div></div><span><font color="#888888">-- <br><div dir="ltr"><div><b><font color="#ff6600">David Raffelt (PhD)</font></b></div><div><font color="#ff6600">Postdoctoral Fellow</font></div>
<div><br></div>
<div>The Florey Institute of Neuroscience and Mental Health</div><div>Melbourne Brain Centre - Austin Campus</div><div>245 Burgundy Street</div><div>Heidelberg Vic 3084<div>Ph: <a value="+61390357024">+61 3 9035 7024</a></div>
</div><div><a value="+61390357024">www.florey.edu.au</a></div></div>
</font></span></div>
</blockquote></div><br></div></div>