Your are correct:<br><br>  typedef unsigned char                                                 PixelType;<br>  typedef itk::Point&lt;PixelType, ImageDimension&gt;            PointType;<br>  typedef itk::Vector&lt;PixelType, ImageDimension&gt;          VectorType;<br>
<br>solved the problem. I looked at the examples in the itksoftwareguide and they used float for the point and vectortype so assumed incorrectly that this was unrelated.<br><br>Another thing. The deformation field looks like this:<br>
<br><a href="http://img18.imageshack.us/img18/5054/deform.jpg">http://img18.imageshack.us/img18/5054/deform.jpg</a><br><br>and like this when I open it with paraview:<br><br><a href="http://img23.imageshack.us/img23/5054/deform.jpg">http://img23.imageshack.us/img23/5054/deform.jpg</a><br>
<br>If I use the calculator to change to vector nothing happens.<br><br>In the itkSoftwareGuide it says:<br><br>&quot;It may be also desirable to write the deformation field as an image of vectors. This can be done<br>with the following code.&quot;<br>
<br>        typedef itk::ImageFileWriter&lt; DeformationFieldType &gt; FieldWriterType;<br>        FieldWriterType::Pointer fieldWriter = FieldWriterType::New();<br>        fieldWriter-&gt;SetFileName( argv[4] );<br>        fieldWriter-&gt;SetInput( filter-&gt;GetOutput() );<br>
        fieldWriter-&gt;Update();<br><br>&quot;Note that the file format used for writing the deformation field must be capable of representing<br>multiple components per pixel. This is the case for the MetaImage and VTK file formats for<br>
example.&quot;<br><br><br>As I understand I need to store the deformation field in another format than PNG (unless this format supports &quot;multiple components per pixel&quot;) to make the conversion to vectors in paraview. Is that correct?<br>
<br><br><br><br><div class="gmail_quote">2009/7/22 Bill Lorensen <span dir="ltr">&lt;<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>&gt;</span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
The source code he posted shows that he is trying to write the type:<br>
<div class="im">typedef itk::Image&lt;VectorType, ImageDimension&gt;<br>
                 DeformationFieldType;<br>
</div>where VectorType is float.<br>
<div class="im"><br>
2009/7/21 Ramón Casero Cañas &lt;<a href="mailto:ramon.casero@comlab.ox.ac.uk">ramon.casero@comlab.ox.ac.uk</a>&gt;:<br>
</div><div><div></div><div class="h5">&gt; Bill Lorensen wrote:<br>
&gt;&gt;<br>
&gt;&gt; The thrown exception is very clear...<br>
&gt;&gt; Exception thrown<br>
&gt;&gt;<br>
&gt;&gt; itk::ExceptionObject (00CAE438)<br>
&gt;&gt; Location: &quot;void __thiscall itk::PNGImageIO::WriteSlice(const class<br>
&gt;&gt; std::basic_string&lt;char,struct std::char_traits&lt;c<br>
&gt;&gt; har&gt;,class std::allocator&lt;char&gt; &gt; &amp;,const void *)&quot;<br>
&gt;&gt; File: ..\..\..\Code\IO\itkPNGImageIO.cxx<br>
&gt;&gt; Line: 473<br>
&gt;&gt; Description: PNG supports unsigned char and unsigned short<br>
&gt;&gt;<br>
&gt;&gt; You must choose an output file type that supports the type you are<br>
&gt;&gt; trying to write. I would suggest .mhd. If you can think of a better<br>
&gt;&gt; exception message, please let us know.<br>
&gt;<br>
&gt;<br>
&gt; Hi Bill,<br>
&gt;<br>
&gt; Thaks for your reply. I don&#39;t think we have expressed any issues about the<br>
&gt; exception itself, but motes said that he&#39;s defining the PixelType as<br>
&gt; unsigned char/short, yet writing the data fails.<br>
&gt;<br>
&gt; Oh! Is this a problem with any of the intermediate steps changing the pixel<br>
&gt; type?<br>
&gt;<br>
&gt; On a related note: Where is the method GetComponentType() called by<br>
&gt; PNGImageIO::WriteSlice declared? I&#39;m sure it&#39;s very obvious where it is, but<br>
&gt; I have been searching for it in the ITK code (grep and [1]), and I couldn&#39;t<br>
&gt; find it.<br>
&gt;<br>
&gt; [1] <a href="http://www.itk.org/Doxygen/html/functions.html" target="_blank">http://www.itk.org/Doxygen/html/functions.html</a><br>
&gt;<br>
&gt; Cheers,<br>
&gt;<br>
&gt; R.<br>
&gt;<br>
&gt; --<br>
&gt; Ramón Casero Cañas, DPhil<br>
&gt;<br>
&gt; Computational Biology, Computing Laboratory<br>
&gt; University of Oxford<br>
&gt; Wolfson Building, Parks Rd<br>
&gt; Oxford OX1 3QD<br>
&gt;<br>
&gt; tlf     +44 (0) 1865 610807<br>
&gt; web     <a href="http://web.comlab.ox.ac.uk/people/Ramon.CaseroCanas" target="_blank">http://web.comlab.ox.ac.uk/people/Ramon.CaseroCanas</a><br>
&gt; photos  <a href="http://www.flickr.com/photos/rcasero/" target="_blank">http://www.flickr.com/photos/rcasero/</a><br>
&gt;<br>
</div></div></blockquote></div><br>