fixed_im = sitk.ReadImage(fixed_fn)
print(fixed_im.GetPixelIDTypeAsString())
fixed_pixelID = fixed_im.GetPixelID()
# set same origins
f_org = fixed_im.GetOrigin()
m_org = moving_im.GetOrigin()
f_sp = fixed_im.GetSpacing()
m_sp = moving_im.GetSpacing()
# fixed_im.SetOrigin (f_sp)
moving_im.SetOrigin(f_org)
moving_im.SetSpacing(f_sp)
diff_fm = sitk.AbsoluteValueDifference(fixed_im, moving_im)
sitk.WriteImage(diff_fm, "diff_before.nii.gz")
# undo the modifications.
moving_im.SetOrigin(m_org)
moving_im.SetSpacing(m_sp)
# elastixImageFilter = sitk.ElastixImageFilter()
# elastixImageFilter.SetFixedImage(fixed_im)
# elastixImageFilter.SetMovingImage(moving_im)
# parameterMapVector = sitk.VectorOfParameterMap()
# parameterMapVector.append(sitk.GetDefaultParameterMap("affine"))
# parameterMapVector.append(sitk.GetDefaultParameterMap("bspline"))
# p=sitk.ParameterMap()
