-------------------------------------------------------------------------
The purpose of this function is to help visualize tractographic
output of 3dTrackID, in particular the fully probabilistic mode
('-mode PROB') that doesn't output a set of 1D tracts for
viewing. Here, local orientation of WM is viewed on a surface that
encloses the tracked region, and the brightness is modulated by
the FA.
There are two halves to this program, each of which can be run
separately. First, V1 and FA values are used to compute RGB ->
HSL coloration. An smoothed isosurface surrounding the tracked
region is computed, and the coloration is mapped onto the surface.
Second, both AFNI and SUMA are opened up with 'talking' turned on,
and the data sets are visualized: the thresholded FA>0.2 map in
AFNI, and the RGB colorized surface in SUMA.
-------------------------------------------------------------------------
OUTPUTS:
1) PREFIX_RGB_HUE.nii.gz :an HSL coloration volume file with four
bricks from the V1 and FA volumes:
[0] Hue
[1] Saturation
[2] Luminosity
[3] Brightness
2) PREFIX_RGB_iso.ply :a slightly smoothed isosurface file made by
IsoSurface
3) PREFIX_RGB_iso.spec :a spec file made by quickspec. Useful
description, huh?
4) PREFIX_RGB.niml.dset :a projection of appropriate coloration onto
the surface
... and a set of AFNI+SUMA commands will also open up viewers and
drive them with appropriate over/underlays and some
probably-useful parameter settings.
-------------------------------------------------------------------------
RUNNING:
@fat_tract_colorize -in_fa FILE_FA -in_v1 FILE_V1 \
-in_tracts FILE_TR -prefix PREFIX \
{ -in_ulay FILE_UL } \
{ -no_view } { -only_view }
where:
-in_fa FILE_FA :FA values of the DT fitting, which can be used to
modulate the brightness of the RGB coloration.
-in_v1 FILE_V1 :first eigenvector of the DT fitting, such as by
3dDWItoDT. The volume is supposed to be a unit
vector with 3 components. The magnitudes of the
components are each between [0, 1], so that
(|x|, |y|, |z|) -> gets mapped to (R, G, B).
-in_tracts FILE_TR :the INDIMAP or PAIRMAP file output by 3dTrackID,
specifying the subbrick as well, if there are >1
in it (you likely need to put the subbrick in
quotes, like NAME_INDIMAP+orig'[0]').
-prefix PREFIX :prefix of all output files.
-in_ulay FILE_UL :optional ability load in a separate data set to
underlay in both the AFNI SUMA viewers (as
'-vol ...' slices in SUMA). For example, you
might want to to load in an anatomical
volume. Default is to use the FA data set.
-no_view :switch to turn off the auto-running of AFNI_SUMA
commands to view the output immediately
-only_view :switch to *only* view the data with AFNI+SUMA.
This assumes that you have run the command at least
once previously, so that there be data to view.
-------------------------------------------------------------------------
EXAMPLE:
# at the moment, relative paths are not allowed for the FA volume
# (because of how it is loaded into the AFNI GUI for viewing as
# part of this command
@fat_tract_colorize -in_fa DT_FA+orig. \
-in_v1 DT_V1+orig. \
-in_tracts o.NETS_AND_000_PAIRMAP+orig'[0]' \
-prefix RGB
-------------------------------------------------------------------------
