backback

Figure 1: Three basic diffusion tensor shapes

setenv FR "2 4 3"

echo "1  1  1  1    3  -4  -3" \
 | emap -i - -fr $FR -amb 0.1 0.1 0.1 -rh -o emap.nrrd

alias SOID soid -sphere -res 20 -r 0.012 -o -

alias OFF2EPS off2eps -i - -fr $FR -or -e emap.nrrd -rh \
  -ur -0.95 0.95 -vr -0.95 0.95 -wd 0.0 0.0 1.5 0.7 0.05

SOID -sc 0.7 0.7 0.7 | OFF2EPS -o demo-sphere.eps 
SOID -sc 0.8 0.18 0.18 | OFF2EPS -o demo-linear.eps
SOID -sc 0.8 0.8 0.18 | OFF2EPS -o demo-planar.eps

rm -f emap.nrrd

foreach SHAPE ( sphere linear planar )
  eps2ppm demo-${SHAPE}.eps 200 \
   | unu project -a 0 -m mean \
   | unu resample -s 256 256 -o demo-${SHAPE}.png
  unu resample -i demo-${SHAPE}.png -s 90 90 -o _demo-${SHAPE}.png
end
unu join -i demo-{sphere,linear,planar}.png -a 0 -o demo3.png
unu join -i _demo-{sphere,linear,planar}.png -a 0 -o _demo3.png