Thank you @zivy, I will try that. maybe I don’t understand how to get the displacement field mapped correctly. I have a [2,2,slices] vectorfloat64 with 3 components per pixel. if the input ijk image has the regular direction [1 0 0 0 1 0 0 0 1], Only a displacement field that is [0 0 z] is required for each image? Does it need to generate a dense deformation field from that? is 4 corners all one needs?
when I resample, should the references image be representative in size and type of the correct spacing I want for the output?
I was able to build an image volume with the proper spacing but has blank slices in the gaps as opposed to interpolated images in between.
I will give it another whirl, thank you so much @lassoan as well.
The thing about the slicer example is that its using the Slicer api, and registers the transform so it updates things once the filter is called, I think.