[ITK] Write and read itkEuler3DTransform with m_ComputeZYX = 1

[Yaniv, Ziv Rafael (NIH/NLM/LHC) [C]]

Hello Silvain,

The patch addressing your issue is up on gerrit for review (http://review.source.kitware.com/#/c/21095/). Hopefully  it will be approved in the next couple of days and the problem will be resolved.

             Ziv

From: Silvain Bériault <silvain.beriault at zimmercas.com<mailto:silvain.beriault at zimmercas.com>>
Date: Wednesday, May 4, 2016 at 11:17 AM
To: "community at itk.org<mailto:community at itk.org>" <community at itk.org<mailto:community at itk.org>>
Subject: [ITK] Write and read itkEuler3DTransform with m_ComputeZYX = 1

Hello,

I and trying to write an itk::Euler3DTransform (with m_ComputeZYX = 1) to file (*.tfm file format). However, I noticed the flag m_ComputeZYX is not saved to the tfm file. Consequently, when I read-back the transform from file, I get a different transform (i.e. a transform with same angles but different matrix since m_ComputeZYX = 0 ). I am using ITK version 4.9.1 (latest release). Any help on this issue would be greatly appreciated. An example follows.

Regards,

Silvain Beriault, Ph.D.

code sample:

include<itkEuler3DTransform.h>
#include<itkTransformFileWriter.h>
#include<itkTransformFileReader.h>
#include<itkMath.h>

intmain()
{
   typedefitk::Euler3DTransform<double> TransformType;
   TransformType::PointertransformIn = TransformType::New();

   transformIn->SetComputeZYX(true);

   doubleangleX = 90;
   doubleangleY = 180;
   doubleangleZ = 0;

   doubledeg2rad = itk::Math::pi / 180.0;

   transformIn->SetRotation(angleX * deg2rad,
                          angleY * deg2rad,
                          angleZ * deg2rad);


   typedefitk::TransformFileWriterTemplate<double> WriterType;
   WriterType::Pointerwriter = WriterType::New();
   writer->SetInput(transformIn);
   writer->SetFileName("D:/tmp/transform.tfm");
   writer->Update();

   typedefitk::TransformFileReaderTemplate<double> ReaderType;
   ReaderType::Pointerreader = ReaderType::New();
   reader->SetFileName("D:/tmp/transform.tfm");
   reader->Update();
   autotransformOut = (reader->GetTransformList()->begin())->GetPointer();


   std::cout << "****TransformIn****" << std::endl;
   std::cout << transformIn << std::endl << std::endl << std::endl;

   std::cout << "****TransformOut****" << std::endl;
   std::cout << *transformOut << std::endl << std::endl << std::endl;
}

The screen output is:

****TransformIn****
Euler3DTransform (0000000000461BF0)
  RTTI typeinfo:   class itk::Euler3DTransform<double>
  Reference Count: 4
  Modified Time: 12
  Debug: Off
  Object Name:
  Observers:
    none
  Matrix:
    -1 1.22465e-016 7.4988e-033
    0 6.12323e-017 -1
    -1.22465e-016 -1 -6.12323e-017
  Offset: [0, 0, 0]
  Center: [0, 0, 0]
  Translation: [0, 0, 0]
  Inverse:
    -1 0 -1.22465e-016
    1.22465e-016 0 -1
    0 -1 0
  Singular: 0
  Euler's angles: AngleX=1.5708 AngleY=3.14159 AngleZ=0
  m_ComputeZYX = 1



****TransformOut****
Euler3DTransform (0000000000473180)
  RTTI typeinfo:   class itk::Euler3DTransform<double>
  Reference Count: 2
  Modified Time: 673
  Debug: Off
  Object Name:
  Observers:
    none
  Matrix:
    -1 0 1.22465e-016
    1.22465e-016 6.12323e-017 1
    -7.4988e-033 1 -6.12323e-017
  Offset: [0, 0, 0]
  Center: [0, 0, 0]
  Translation: [0, 0, 0]
  Inverse:
    -1 1.22465e-016 -7.4988e-033
    -7.4988e-033 4.59169e-049 1
    1.22465e-016 1 4.59169e-049
  Singular: 0
  Euler's angles: AngleX=1.5708 AngleY=3.14159 AngleZ=0
  m_ComputeZYX = 0


The saved transform is :

#Insight Transform File V1.0
#Transform 0
Transform: Euler3DTransform_double_3_3
Parameters: 1.5707963267948966 3.141592653589793 0 0 0 0
FixedParameters: 0 0 0

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/community/attachments/20160505/4b66c194/attachment-0001.html>