Soren Christensen

Hi, If this is RGB images and you are using bspline interpolation for the final resampling then you can get values outside of the valid range (do you could have...

Hi, You can switch on output of the iteration log file for this: SimpleElastix.SetParameter("WriteIterationInfo",["true"]) and then td = tempfile.mkdtemp() SimpleElastix.SetOutputDirectory(td) .. execute registration .... then xfms = sorted(glob.glob(os.path.join(td, 'TransformParameters.*.R*'))) dbg_xfms...

I seem to have this issue too. Here's a self contained minimal example below. I don't see an obvious way to incorporate the workaround above. Cheers, Soren `myimg=sitk.GetImageFromArray(np.zeros((64,64,64))) sitk.WriteImage(myimg,'myimg.nii') #...

Which transform are you using? Does the error message apply to your input? **The number of pixels along direction 2 is less than 4. This filter requires a minimum of...

Then you are using the default transform(s) which you can inspect here sitk.PrintParameterMap(se.GetParameterMap()) #(where se=sitk.ElastixImageFilter()) So it seems as if that applies a gaussian blur which requires 3D data whereas...