setenv CAM "-up 0 0 1 -rh -fr 5 6 5"
echo "1 0.7 0.7 0.7 4 -2.0 -4" \
| emap -i - -amb 0.3 0.3 0.3 $CAM -o tmp0.nrrd
unu 2op ^ tmp0.nrrd 1.55 | 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.7 tmp0.nrrd tmp1.nrrd \
| unu 3op clamp 0 - 1 \
| unu 3op lerp 0.6 1 - -o emap.nrrd
rm -f tmp0.nrrd tmp1.nrrd
alias SOID soid -res 16 -r 0.0
alias OFF2EPS off2eps -i - -or -e emap.nrrd -rh \
$CAM -ur -0.9 0.9 -vr -0.58 0.58 -wd 0 0 1.5 0.4 0.4
SOID -sh 3 -os 1.55 -sc 0.5983 0.2683 0.1333 -o - | OFF2EPS -o seam-A.eps
SOID -sh 3 -os 1.55 -sc 0.5883 0.2783 0.1333 -o - | OFF2EPS -o seam-B.eps
SOID -sh 3 -os 1.55 -sc 0.5783 0.2883 0.1333 -o - | OFF2EPS -o seam-C.eps
SOID -sh 3 -os 1.55 -sc 0.5683 0.2983 0.1333 -o - | OFF2EPS -o seam-D.eps
rm -f emap.nrrd
unsetenv TEEM_ROOT
unsetenv CAM
unsetenv PARM
unsetenv L
unset SH
foreach X ( A B C D )
eps2ppm seam-${X}.eps 300 \
| unu project -a 0 -m mean \
| unu resample -s x0.207 x0.207 -o seam-${X}.png
end
unu join -i seam-{A,B,C,D}.png -a 0 -o seam.png
unu resample -i seam.png -s x0.25 x0.25 -o _seam.png
back