echo "-1 1 -1 1" \
 | unu reshape -s 2 2 \
 | unu resample -s 25 25 -k tent -c node -o x

echo "1 1 -1 -1" \
 | unu reshape -s 2 2 \
 | unu resample -s 25 25 -k tent -c node -o y

unu 2op x -1 y \
 | unu join -i - x -a 0 -incr -o vl

unu 2op x -1 x \
 | unu join -i y - -a 0 -incr -o vr

unu join -i x y -a 0 -incr -o za

unu join -i x y -a 0 -incr | unu 2op x - -1 -o zb

unu join -i y x -a 0 -incr -o zc

unu join -i y x -a 0 -incr | unu 2op x - -1 -o zd

unu 1op exists -i x \
 | unu pad -min -6 -6 -max M+6 M+6 -b pad -v 0 \
 | unu resample -s x1 x1 -k cubic:5,1,0 \
 | unu resample -s 25 25 -k tent \
 | unu quantize -b 8 \
 | unu 2op / - 255 -t float \
 | unu axinsert -a 0 \
 | unu pad -min 0 0 0 -max 1 M M -o m

foreach D ( vl vr za zb zc zd )
  unu 2op x $D m -o $D
end

unu join -i vl vr -a 2 -o 0
unu join -i za zb -a 2 -o 1
unu join -i zc zd -a 2 -o 2
unu join -i 0 1 2 -a 1 \
 | unu resample -s = 75 50 -k cubic:1,0 -b wrap \
 | unu 2op x - 3 -o tmp.nhdr

echo "space dimension: 2" >> tmp.nhdr
echo "space origin: (0.02,1.98)" >> tmp.nhdr
echo "space directions: none (0.04,0) (0,-0.04)" >> tmp.nhdr
echo "kinds: 2-vector space space" >> tmp.nhdr

unu project -i tmp.nhdr -a 0 -m l2 \
 | unu resample -s 300 200 \
 | unu pad -min -2 -2 -max M+2 M+2 -b pad -v 0.0 \
 | unu 1op rand -o vorttest-rand.nrrd

unu pad -i tmp.nhdr -min 0 -2 -2 -max M M+2 M+2 -o vorttest.nrrd

rm -f 0 1 x y vl vr za zb m tmp.nhdr tmp.raw