[Insight-users] Display deformation field

Karthik Krishnan karthik.krishnan at kitware.com
Mon Apr 14 15:40:40 EDT 2008


Why don't you just change the meta image header that you mailed
earlier. The datapointer is going to be the same if the extents in the
third dimension are 1. So your new header is.

ObjectType = Image
NDims = 3
BinaryData = True
BinaryDataByteOrderMSB = False
CompressedData = False
Offset = 0 0 0
CenterOfRotation = 0 0 0
ElementSpacing = 1 1 1
DimSize = 768 576 1
AnatomicalOrientation = ??
ElementNumberOfChannels = 3
ElementType = MET_FLOAT
ElementDataFile = kalip_grey_Deformtionfield_ReciprocalSpline_new.raw


On Mon, Apr 14, 2008 at 3:17 PM, Boettcher, Dr. Peter
<Boettcher at kleintierklinik.uni-leipzig.de> wrote:
>
>
>
>
>
> When computing deformation fields on a 2D-image it seems to me, that by
> declaring the a two-dimensional vector ParaView can not properly display the
> computed deformation field. Therefore I declared a 3D-vector  instead of the
> 2D-vector in the registation examples Computation of the deformation field
> is not affected by that change and the program runs fine (I get the same
> deformation as with a 2D-vector). However the third component of the vector
> (which is not used) is filled with enormous values (e+8). Is there an easy
> method to fill the third component of the vector with a zero value (sfter
> the deformationfieldsource has been updated)?
>
>
>
> const     unsigned int   Dimension = 2;
>
> typedef   float          VectorComponentType;
>
> typedef   itk::Vector< VectorComponentType, 3 >    VectorType;
>
> typedef   itk::Image< VectorType,  Dimension >   DeformationFieldType;
>
>
>
>
>
> Regards, Peter.
>
>
> ------------------------------------
> Peter Böttcher, Dr med vet, DipECVS
> European Veterinary Specialist in Surgery
> Fachtierarzt für Kleintierchirurgie
> Klinik für Kleintiere
> Universität Leipzig
> An den Tierkliniken 23
> D-04103 Leipzig (Germany)
> Tel: +49-341-9738700
> Fax: +49-341-9738799
> email: boettcher at kleintierklinik.uni-leipzig.de
>
> _______________________________________________
>  Insight-users mailing list
>  Insight-users at itk.org
>  http://www.itk.org/mailman/listinfo/insight-users
>
>


More information about the Insight-users mailing list