[Mrtrix-discussion] calculating recovered HARDI / residuals

[Donald Tournier]

Hi Kerstin,

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...

Get in touch if you get stumped...
Cheers, and good luck!

Donald.


On 10 June 2010 11:04, Kerstin Pannek <k.pannek1 at uq.edu.au> wrote:

> Hi
>
> I would like to calculate the residuals of my HARDI data for bootstrapping.
> 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?
>
> Thanks
> Kerstin
> _______________________________________________
> Mrtrix-discussion mailing list
> Mrtrix-discussion at www.nitrc.org
> http://www.nitrc.org/mailman/listinfo/mrtrix-discussion
>



-- 
Jacques-Donald Tournier (PhD)
Brain Research Institute, Melbourne, Australia
Tel: +61 (0)3 9496 4078
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.nitrc.org/pipermail/mrtrix-discussion/attachments/20100610/83d7ba92/attachment.html