[vtkusers] grid reslice

[David Gobbi]

Hi Ryan,

Can you attach your image both with and _without_ deformation?
I looked at the image but couldn't make any sense of it.

 - David


On Mon, Jan 26, 2015 at 6:29 PM, Young, Ryan <
Ryan.Young at seattlechildrens.org> wrote:

>  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/20150126/e9990112/attachment.html>