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 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g sqd -sh 1.5 -gsc 0.123 -gr 22 -wd 0.5 0.0 0.0 -vr -0.84 0.81 \
-o barysqdA.eps
tt -n 6 -p 0.45 -0.2 0.35 -ca $CA -v 1 -hack 0 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g sqd -sh 3.0 -gsc 0.123 -gr 22 -wd 0.5 0.0 0.0 -vr -0.84 0.81 \
-o barysqdB.eps
tt -n 6 -p 0.45 -0.2 0.35 -ca $CA -v 1 -hack 0 \
| tend norm -w $WGT \
| tend glyph -sat 0 \
$CAM $PARM \
-nobg -g sqd -sh 6.0 -gsc 0.123 -gr 22 -wd 0.5 0.0 0.0 -vr -0.84 0.81 \
-o barysqdC.eps
rm -f emap.nrrd
unsetenv TEEM_ROOT
unsetenv WGT
unsetenv CAM
unsetenv PARM
unsetenv CA
foreach G ( A B C )
eps2ppm barysqd${G}.eps 300 \
| unu project -a 0 -m mean \
| unu resample -s x0.2488 x0.2488 -o barysqd${G}.png
unu resample -i barysqd${G}.png -s x0.3516 x0.3516 -o _barysqd${G}.png
end
unu join -i barysqd{A,B,C}.png -a 0 -o barysqd.png
unu join -i _barysqd{A,B,C}.png -a 0 -o _barysqd.png
back