unu resample -i ddro-80.nrrd -s x1 =  -k cubicd:1,0 -o dx
unu resample -i ddro-80.nrrd -s =  x1 -k cubicd:1,0 -o dy

unu join -i dx dy -a 0 -incr -o tmp

unu axinsert -i ddro-80.nrrd -a 0 \
 | unu pad -min 0 0 0 -max 1 M M \
 | unu axinfo -a 0 -k 2-vector \
 | unu inset -i - -s tmp -min 0 0 0 \
 | ../bin/dnorm -i - -o ddro-80-grad.nrrd




unu resample -i ddro-80.nrrd -s x1 =  -k cubicdd:1,0 -o dxx
unu resample -i ddro-80.nrrd -s x1 x1 -k cubicd:1,0 -o dxy
unu resample -i ddro-80.nrrd -s =  x1 -k cubicdd:1,0 -o dyy

unu join -i dxx dxy dxy dyy -a 0 -incr -o tmp

unu axinsert -i ddro-80.nrrd -a 0 \
 | unu pad -min 0 0 0 -max 3 M M \
 | unu axinfo -a 0 -k 2D-matrix \
 | unu inset -i - -s tmp -min 0 0 0 \
 | ../bin/dnorm -i - -o ddro-80-hess.nrrd





rm -f dx dy dxx dxy dyy tmp