I’m using GridBasedImageResampling using the same input image as in FimeRegistration and the displacement image calculated in that process.
I get “Requested region is (at least partially) outside the largest possible region of the displacement field”, nevertheless I’m using the size of the input image as out.sizex and out.sizey, so this should work.
BTW, why are not these values taken by default from the input image? These values are not requested by FineRegistration, which can also warp the input image (and that output also has the same size as I define for GridBasedImageResampling)
The reason of using GridBasedImageResampling is to apply the same displacement image (output of FineRegistration) to other images, thus saving the time of computing the displacement and making sure I apply always the same displacement.