; SHIFTVOL.BAT
; Compares two volumes. If cross correlation peak is not at center
; of volume, X,Y,Z shift values are written to a doc file.
;
; Inputs:
;  reference_volume : reference volume
;  test_volume : this volume will be shifted if they are not aligned.
;

FR G
<ref>../reference_volume

FR G
<vol>test_volume  

FR G
<out>aligned_volume  ; shifted version of test_volume

CC N
<vol>
<ref>
_1

PK 3D X11,X12,X13,x14,x15,x16,x17
_1
+
(1,0)
N
N
shiftdoc   ; doc file with shift values

; The first 3 columns in the docfile (shiftdoc) are the X,Y,Z shifts
; (integer values) needed to bring test_volume into alignment with
; the reference_volume.
; X14,x15,x16 also contain these values in floating point format.

SH
<vol>   ; volume to be shifted
<out>   ; aligned_volume
-x14,-x15,-x16

EN D