# make sure it looks right post-reshape
unu reshape -i synth-u.txt -s 2 9 9 | unu slice -a 0 -p 0 | unu save -f text
unu reshape -i synth-u.txt -s 2 9 9 | unu slice -a 0 -p 1 | unu save -f text

unu reshape -i synth-u.txt -s 2 9 9 \
 | unu axinfo -a 0 -k 2-vector \
 | ~/diderot/diderot/trunk/bin/dnorm -i - -o - \
 | unu pad -min 0 -3 -3 -max M M+3 M+3 -o imp-u.nrrd

unu shuffle -i imp-u.nrrd -a 0 -p 1 0 -o imp-v.nrrd

# upsampling with bspln3 == cubic:1,0
unu crop -i imp-u.nrrd -min 0 3 3 -max M M-3 M-3 \
 | unu resample -s = 500 500 -k cubic:1,0 -c node -co \
 | unu pad -min 0 0 0 -max 2 M M -b wrap -o teem.nrrd

# upsampling with ctmr == cubic:0.0.5
unu crop -i imp-u.nrrd -min 0 3 3 -max M M-3 M-3 \
 | unu resample -s = 500 500 -k cubic:0,0.5 -c node -co \
 | unu pad -min 0 0 0 -max 2 M M -b wrap -o teem.nrrd

# 1st derivative along X  d/dx ( cubic:B,C) == cubicd:B,C
unu crop -i imp-u.nrrd -min 0 3 3 -max M M-3 M-3 \
 | unu resample -s = 500 = -k cubicd:0,0.5 -c node -co \
 | unu resample -s = = 500 -k cubic:0,0.5 -c node -co \
 | unu slice -a 0 -p 0 -o tmp-dx.nrrd

unu crop -i imp-u.nrrd -min 0 3 3 -max M M-3 M-3 \
 | unu resample -s = 500 = -k cubic:0,0.5 -c node -co \
 | unu resample -s = = 500 -k cubicd:0,0.5 -c node -co \
 | unu slice -a 0 -p 0 -o tmp-dy.nrrd

unu join -i tmp-d{x,y,x}.nrrd -a 0 -incr -o teem.nrrd