python - <<'PY'
import tifffile,cv2,numpy as np
s='HH-1206';a=tifffile.imread('/workspace/input/photos/'+s+'.tif').astype('float32')/65535;b=tifffile.imread('/workspace/input/charts/'+s+'_chart.tif').astype('float32')/65535
for rect in [(0,1100,400,1550),(1900,1000,2300,1550),(0,0,400,1000)]:
 x0,y0,x1,y1=rect;aa=a[y0:y1,x0:x1].mean(2);bb=b[y0:y1,x0:x1].mean(2);print(rect,cv2.phaseCorrelate(aa,bb),np.std(aa-bb),np.mean(aa-bb))
PY