setenv CAM "-fr 0 0 -30 -up 0 1 0 -rh"
echo "1 0.86 0.86 0.86 2.0 -2.0 -4" \
| emap -i - -amb 0.1 0.1 0.1 $CAM -o tmp0.nrrd
unu 2op ^ tmp0.nrrd 1.5 | unu 2op - 1 - | unu 2op x - 3.14159 | unu 1op cos \
| unu 2op + - 1 | unu 2op / - 1.8 -o tmp1.nrrd
unu 3op lerp 0.5 tmp0.nrrd tmp1.nrrd -o emap.nrrd
rm -f tmp0.nrrd tmp1.nrrd
setenv PARM "-psc 150 -ur -1.13 1.13 -vr -0.98 0.9 -or -atr 0 -emap emap.nrrd"
setenv CA 0.9
setenv WGT "2 1.5 0.4"
tt -n 6 -p 0.45 -0.2 0.35 -ca $CA -v 1 -hack 0.06 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g box -gsc 0.110 -gr 20 -wd 0.5 0.2 0.0 \
-o barybox.eps
tt -n 6 -p 0.45 -0.2 0.35 -ca $CA -v 1 -hack 0.05 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g cyl -gsc 0.115 -gr 100 -wd 0.5 0.2 0.0 \
-o barycyl.eps
tt -n 6 -p 0.45 -0.2 0.35 -ca $CA -v 1 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g sph -gsc 0.123 -gr 22 -wd 0.5 0.0 0.0 \
-o barysph.eps
rm -f emap.nrrd
unsetenv TEEM_ROOT
unsetenv WGT
unsetenv CAM
unsetenv PARM
unsetenv CA
foreach SHP ( box cyl sph )
eps2ppm bary${SHP}.eps 300 \
| unu project -a 0 -m mean \
| unu resample -s x0.217 x0.217 -o bary${SHP}.png
unu resample -i bary${SHP}.png -s x0.3516 x0.3516 -o _bary${SHP}.png
end
unu join -i bary{box,cyl,sph}.png -a 0 -o bary3.png
unu join -i _bary{box,cyl,sph}.png -a 0 -o _bary3.png
back