Print

Print


Tim,

Brilliant! Great colourful images now, thanks.

BTW, I find that BET doesn't work _perfectly_, and the step of making the
binary mask is a bit of a pain. Other software I've tried use a threshold to
mask off non-brain tissues. Could FDT possibly make such a capability
available?

Thanks

J

----- Original Message -----
From: "Tim Behrens" <[log in to unmask]>
To: <[log in to unmask]>
Sent: Sunday, October 17, 2004 4:38 PM
Subject: Re: [FSL] DTIFit


> Hi there -
>
> On Sun, 17 Oct 2004, Jay Ives wrote:
> >
> > Can someone tell me how to fill in the bvecs file correctly? The
> > information I have about the gradient directions in the 30 b direction
> > sequence I am using shows x, y, and z the absolute values of which add
up
> > to one as in the exerpt below.....
> >
> > 1  0.4356 -0.226576 -0.770884 0.004225 0.734449
> > 0  0.564001 0.343396 0.145161 0.613089 0.112896
> > 0  0         0.430336 0.0841        -0.381924 -0.152881
> >
>
>
> No this is not quite right. You need to normalise each gradient direction
> such that its _magnitude_ is equal to 1.
>
> To do this, for each gradient direction, r, compute
>
> r_norm=r/sqrt( r(1)^2+r(2)^2+r(3)^2 );
>
> This is not the same as ensuring that the sum of absolute values in the
> vector is equal to 1.
>
> > I get a good FA map with these, but the V1_RGB looks very smearedn and
the
> > colours are mostly pretty pastels, rather than good strong blue, green,
> > and red.
>
> This is because the display range in fslview is set wrong. If you
> set it as 0->1 you will get nice looking strong primary colours.
>
> I'm not sure about the TCL error in the GUI - we'll look into that next
> week.
>
>
> Hope this helps
>
> Tim
>
>
> >
> > Also, I can run the tools from a command line , but I get an error if I
> > try to use the GUI. Here is the message.....
> >
> > bad cursor spec " watch red white "
> > bad cursor spec " watch red white "
> >     (processing "-cursor" option)
> >     invoked from within
> > "$w configure -cursor { watch red white } "
> >     (procedure "fdt_monitor" line 5)
> >     invoked from within
> > "fdt_monitor $w "${FSLDIR}/bin/eddy_correct $eddy(input) $eddy(output)
> > $eddy(refnum)""
> >     ("eddy_current" arm line 18)
> >     invoked from within
> > "switch -- $tool {
> >  eddy_current {
> >      global eddy
> >
> >      set errorStr ""
> >      if { $eddy(input) == "" } { set errorStr "You need to specify
> > the input ..."
> >     (procedure "fdt:apply" line 7)
> >     invoked from within
> > "fdt:apply .fdt keep"
> >     invoked from within
> > ".fdt.apply invoke"
> >     ("uplevel" body line 1)
> >     invoked from within
> > "uplevel #0 [list $w invoke]"
> >     (procedure "tkButtonUp" line 9)
> >     invoked from within
> > "tkButtonUp .fdt.apply
> > "
> >     (command bound to event)
> >
> > Thanks for this great tool and any assistance you can offer to get me up
> > and running.
> >
> > J
> >
>
>