[vtkusers] grid reslice

[Young, Ryan]

Hello,

I have a vector field that has the displacement vector for each voxel in a volume.  I am trying to apply this deformation and visualize the results.  The difficulty I am have is that the deformation only exist for the voxels inside the volume, not globally in the whole space. If I set the displacement to zeros outside the volume I get poor results.

My toy example is just do to a 10 voxel shift in the x direction The volume is padded with zeros so there is plenty of room.  The result is an image missing about ten voxels.  See attached picture.  If I make a similar deformation that exist at all voxels, the result is as expect.  Any advice. Also, a code snippet is below.

Thanks in advanced
Ryan

grid=vtk.vtkGridTransform()
grid.SetDisplacementGrid(defImporter.GetOutput())
grid.SetInterpolationModeToLinear()

reslice=vtk.vtkImageReslice()
 reslice.SetWrap(1)
 reslice.SetResliceTransform(grid)
 reslice.SetInterpolationModeToLinear()
 reslice.SetOutputSpacing(1,1,1)
 reslice.SetOutputOrigin(0,0,0)
 reslice.SetInputConnection(volImporter.GetOutputPort())
CONFIDENTIALITY NOTICE: This e-mail message, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information protected by law. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply e-mail and destroy all copies of the original message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20150127/c2540223/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: deformed_volume.png
Type: image/png
Size: 14726 bytes
Desc: deformed_volume.png
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20150127/c2540223/attachment.png>