setenv FR "2.5 4 1.8"
echo "1 1 1 1 3 -4 -3" \
| emap -i - -fr $FR -amb 1 1 1 -rh -o emap.nrrd
alias SOID soid -sc 0.75 0.75 0.75 -r 0.0 -res 22 -o -
alias OFF2EPS off2eps -i - -fr $FR -or -e emap.nrrd \
-rh -ur -0.99 0.99 -vr -0.99 0.99 -wd 0.0 0.0 2.0 1.3 1.3 \
-ws 370 -concave -nobg
echo superquad 0
SOID -AB 0.25 0.25 | OFF2EPS -o super00.eps
SOID -AB 0.50 0.25 | OFF2EPS -o super10.eps
SOID -AB 1.00 0.25 | OFF2EPS -o super20.eps
SOID -AB 2.00 0.25 | OFF2EPS -o super30.eps
SOID -AB 4.00 0.25 | OFF2EPS -o super40.eps
echo superquad 1
SOID -AB 0.25 0.50 | OFF2EPS -o super01.eps
SOID -AB 0.50 0.50 | OFF2EPS -o super11.eps
SOID -AB 1.00 0.50 | OFF2EPS -o super21.eps
SOID -AB 2.00 0.50 | OFF2EPS -o super31.eps
SOID -AB 4.00 0.50 | OFF2EPS -o super41.eps
echo superquad 2
SOID -AB 0.25 1.00 | OFF2EPS -o super02.eps
SOID -AB 0.50 1.00 | OFF2EPS -o super12.eps
SOID -AB 1.00 1.00 | OFF2EPS -o super22.eps
SOID -AB 2.00 1.00 | OFF2EPS -o super32.eps
SOID -AB 4.00 1.00 | OFF2EPS -o super42.eps
echo superquad 3
SOID -AB 0.25 2.00 | OFF2EPS -o super03.eps
SOID -AB 0.50 2.00 | OFF2EPS -o super13.eps
SOID -AB 1.00 2.00 | OFF2EPS -o super23.eps
SOID -AB 2.00 2.00 | OFF2EPS -o super33.eps
SOID -AB 4.00 2.00 | OFF2EPS -o super43.eps
echo superquad 4
SOID -AB 0.25 4.00 | OFF2EPS -o super04.eps
SOID -AB 0.50 4.00 | OFF2EPS -o super14.eps
SOID -AB 1.00 4.00 | OFF2EPS -o super24.eps
SOID -AB 2.00 4.00 | OFF2EPS -o super34.eps
SOID -AB 4.00 4.00 | OFF2EPS -o super44.eps
rm -f emap.nrrd
foreach IJEPS ( super??.eps )
set IJ=${IJEPS:r}
echo ========== $IJ ===============
eps2ppm ${IJ}.eps 100 \
| unu project -a 0 -m mean \
| unu resample -s 128 128 -o ${IJ}.png
end
unu join -i super??.png -a 2 \
| unu axsplit -a 2 -s 5 5 \
| unu flip -a 2 \
| unu permute -p 0 3 1 2 \
| unu axmerge -a 0 | unu axmerge -a 1 -o super.png
unu resample -i super.png -s 120 120 -o _super.png
rm -f super??.png
back