Hi Kerstin,<div><br></div><div>There is no code to do that currently, but it really shouldn't be too hard. The simplest thing to do would be to copy and modify cmd/sdeconv.cpp (the simple linear spherical deconvolution), to produce the modeled DW signals. You'll need to invert the operation, so that the SH coefficients file is the input, and the DWI file is the output. You can remove a bunch of options that you won't need (lmax, mask, filter & normalise) and all the associated code that depends on those options. You'll be determining lmax from the CSD file, using DWI::SH::LforN (SH.dim(3)), and instead of converting from amplitudes to SH, you'll be doing the opposite operation - at line 213, use SHT.SH2A (sigs, res). There's obviously quite a lot more work to it, but the rest should be *relatively* trivial...</div>
<div><br></div><div>Get in touch if you get stumped...</div><div>Cheers, and good luck!</div><div><br></div><div>Donald.</div><div><br><br><div class="gmail_quote">On 10 June 2010 11:04, Kerstin Pannek <span dir="ltr"><<a href="mailto:k.pannek1@uq.edu.au">k.pannek1@uq.edu.au</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi<br>
<br>
I would like to calculate the residuals of my HARDI data for bootstrapping.<br>
Is there a way in mrtrix to recover the HARDI signal from the csd, or could you let me know how to bring csd, response and gradient directions together to get the HARDI signal?<br>
<br>
Thanks<br>
Kerstin<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>
</blockquote></div><br><br clear="all"><br>-- <br>Jacques-Donald Tournier (PhD)<br>Brain Research Institute, Melbourne, Australia<br>Tel: +61 (0)3 9496 4078<br>
</div>