<div dir="ltr"><div>Hi,<br><br></div>Thank you for your answer.<br><div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word"><div><br>
</div>
<div>We'll take a look into making the necessary adjustments to Camino to handle the spatially varying b-value; probably a bit fiddly, but doable. In the meantime, I think you can get a reasonable feel for how well the different processing techniques will
work by just ignoring the spatial variation - not a long term solution though.</div></div></blockquote><div><br></div><div>That sounds good, I am looking forward to the update. <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
<div><br>
</div>
<div>I see no reason why shredder won't work with longs. Its three arguments are just in bytes, so can't you just use a <chunksize> of 8?</div></div></blockquote><div><br></div><div>I am sorry I am traveling at the moment so I can't reproduce the error remotely. However I remember that calling the shredder with long arguments simply resulted in <null> output. <br>
</div><div><br></div><div>NB. I used an older version of <a href="http://web4.cs.ucl.ac.uk/research/medic/camino/pmwiki/pmwiki.php?n=Tutorials.ParallelPASTutorial"><parallel processing tutorial></a> as a basis for my scripting. A nice aspect about the tutorial is that the data is split in several slices with the shredder function, which is convenient for displaying the reconstruction. That way the <space> (sometimes also the offset) argument necessarily has to be >long<. Naturally this could be overcome by shredding the diffusion data differently. However 'cat'ing them back together would require some more tweaking.<br>
<br></div><div>With kind regards,<br><br> Emil<br></div><div><br><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">
<div><br>
</div>
<div>Danny</div><div><div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<div>
<div>On 29 Nov 2013, at 13:58, E.H.J. Nijhuis <<a href="mailto:emil.nijhuis@gmail.com" target="_blank">emil.nijhuis@gmail.com</a>> wrote:</div>
<br>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>Dear Danny and Philip,<br>
<br>
thank you very much for your quick and very helpful answer.<br>
<br>
</div>
I am aiming to produce streamlines which I can use for a whole brain connectivity analysis for the publicly available HCP data (multi shell data with b-values of about 0, 1000, 2000 and 3000). With the dataset it was possible to estimate FA values in camino
but I ran into difficulties when processing the data with the PAS approach. <br>
<br>
The first thing which was problematic when processing a HCP dataset was that the shredder function is not able to accept longs (raw hcp.bfloats are rather big >4GB). There was no easy and clean fix to the issue because the function uses DataInputStream.read()
which cannot accept longs either. I overcame the issue by downcasting longs to int somewhere posthoc in the source, but that is certainly not robust enough for an actual fix.<br>
<br>
</div>
The error considering the lag of b=0 images was with the mesd function which I called with the following options:<br>
<div><br>
mesd -schemefile hcp.scheme -filter PAS 1.4 -fastmesd -mepointset 16 -bgmask hcp.Bchar -maskdatatype ubyte < hcp_SliceX.Bfloat > hcp_SliceX.Bdouble<br>
<br>
</div>
<div>This resulted in the warning "WARNING: No b=0 data, cannot normalize", while the output was essentially blank. I guess if I set the lowest b-values to 0 it will work, but given your elaborate answer I now feel reluctant to process the data any further
because of the caveats with the spatially varying b-value and gradient directions.<br>
<br>
</div>
<div>Which brings me to my next set of questions:<br>
<br>
</div>
<div>1. Are you planning to bring out an update which will be able to process the HCP dataset (or more generally datasets with spatially varying gradients)?<br>
</div>
<div>2. If the goal is to produce streamlines for a full brain connectivity analysis: what would be your preferred way of processing for a multi-shell dataset? (I know it is rather impossible to answer this question properly, but I would be very interested
in your and others users opinion and advice).<br>
<br>
</div>
<div>Thank you a lot for your consideration and expertise.<br>
</div>
<div><br>
</div>
<div>With kind regards,<br>
<br>
Emil<br>
<div><br>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Mon, Nov 25, 2013 at 5:28 PM, Alexander, Daniel <span dir="ltr">
<<a href="mailto:d.alexander@ucl.ac.uk" target="_blank">d.alexander@ucl.ac.uk</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 style="word-wrap:break-word">What are you trying to do? Most model fitting routines should work fine without a b=0, but I guess some of the other reconstruction algorithms may demand a b=0 image. Assuming the low b-value images in the HCP data are b=0,
ie set the entries in the schemefile to zero, I would expect to give very similar results to accounting for the slightly non-zero b.
<div><br>
</div>
<div>When you say "DSI acquisition", do you really mean DSI, ie a grid of locations in q-space? HCP use a three-shell HARDI protocol, which is quite different. Or do you have some special data? </div>
<div><br>
</div>
<div>I don't know if anyone has tried processing the standard HCP data using Camino. Most Camino algorithms should run fine on data acquired with their protocol, although currently Camino doesn't account for the spatially varying b-value and gradient directions
in the HCP data. I am not sure how much of an issue that is in terms of results, but you should at least be able to get an output to look at.</div>
<div><br>
</div>
<div>Danny</div>
<div><br>
</div>
<div>PS: Please write to <a href="mailto:camino-users@www.nitrc.org" target="_blank">camino-users@www.nitrc.org</a> for Camino questions rather than direct to me. You will get faster and better answers!</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>
<div>On 25 Nov 2013, at 10:06, "E.H.J. Nijhuis" <<a href="mailto:emil.nijhuis@gmail.com" target="_blank">emil.nijhuis@gmail.com</a>></div>
<div> wrote:</div>
<br>
<blockquote type="cite">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>Dear Daniel,<br>
<br>
</div>
I recently have attempted to process the diffusion dataset of the human connectome project.<br>
<br>
</div>
Unfortunately I noticed that camino was unable to process this dataset due to the lag of actual b=0 values (see attached the sample fsl and corresponding camino tables). The small b-value variations are due certain hardware aspects (see Sotiropoulos et al.,
Neuroimage 2013). <br>
<br>
</div>
<div>The lag of actual b0 values results in the following error messages during reconstruction:<br>
WARNING: No b=0 data, cannot normalize<br>
<br>
</div>
Could you please advice me how to best handle the HCP data? Simply round the 'b0'-values? Also, since it is an DSI aquisition protocol, would you still recommend me to use camino?<br>
<br>
</div>
<div>Thank you a lot for your expertise and consideration.<br>
</div>
<div><br>
</div>
Yours sincerely,<br>
<br>
Emil<br clear="all">
<div><br>
-- <br>
<div dir="ltr"><span style="font-family:arial,helvetica,sans-serif">----------------------------------------------<br>
Postdoctoral Researcher - MR Techniques in Brain Function <br>
<b style="font-weight:normal">Donders Institute: Centre for Cognitive Neuroimaging</b><br>
<br>
tel <a href="tel:%2B31%20%280%2924%2036%2010885" value="+31243610885" target="_blank">
+31 (0)24 36 10885</a><br>
fax <a href="tel:%2B31%20%280%2924%2036%2010989" value="+31243610989" target="_blank">
+31 (0)24 36 10989</a><br>
web <a href="http://www.ru.nl/donders/" target="_blank">http://www.ru.nl/donders/</a><br>
<br>
Postal address: Trigon 204, P.O. Box 9101, 6500HB Nijmegen, The Netherlands<br>
Visiting address: Trigon, Kapittelweg 29, 6525EN Nijmegen, The Netherlands </span>
<br>
</div>
</div>
</div>
</div>
<span><bvals></span><span><bvecs></span><span><hcp100408.scheme></span></blockquote>
</div>
<br>
</div>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
<div dir="ltr"><span style="font-family:arial,helvetica,sans-serif">----------------------------------------------<br>
Postdoctoral Researcher - MR Techniques in Brain Function <br>
<b style="font-weight:normal">Donders Institute: Centre for Cognitive Neuroimaging</b><br>
<br>
tel <a href="tel:%2B31%20%280%2924%2036%2010885" value="+31243610885" target="_blank">
+31 (0)24 36 10885</a><br>
fax <a href="tel:%2B31%20%280%2924%2036%2010989" value="+31243610989" target="_blank">
+31 (0)24 36 10989</a><br>
web <a href="http://www.ru.nl/donders/" target="_blank">http://www.ru.nl/donders/</a><br>
<br>
Postal address: Trigon 204, P.O. Box 9101, 6500HB Nijmegen, The Netherlands<br>
Visiting address: Trigon, Kapittelweg 29, 6525EN Nijmegen, The Netherlands </span>
<br>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div></div></div>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><span style="font-family:arial,helvetica,sans-serif">----------------------------------------------<br>Postdoctoral Researcher - MR Techniques in Brain Function <br>
<b style="font-weight:normal">Donders Institute: Centre for Cognitive Neuroimaging</b><br><br>tel <a href="tel:%2B31%20%280%2924%2036%2010885" value="+31243610885" target="_blank">+31 (0)24 36 10885</a><br>fax <a href="tel:%2B31%20%280%2924%2036%2010989" value="+31243610989" target="_blank">+31 (0)24 36 10989</a><br>
web <a href="http://www.ru.nl/donders/" target="_blank">http://www.ru.nl/donders/</a><br>
<br>Postal address: Trigon 204, P.O. Box 9101, 6500HB Nijmegen, The Netherlands<br>Visiting address: Trigon, Kapittelweg 29, 6525EN Nijmegen, The Netherlands
</span><br></div>
</div></div></div>