[Insight-users] transformreader cannot read written transform of transformwriter

[Lodron, Gerald]

Hi

Following problem:




std::cout << "Matrix is orthogonal" << oRigid->MatrixIsOrthogonal(oRigid->GetMatrix(),1e-10) << std::endl;

itk::TransformFileWriter::Pointer oWriter = itk::TransformFileWriter::New();

oWriter->SetInput(oRigid);

oWriter->SetFileName(oPathName.str().c_str());

oWriter->Update();

itk::TransformFileReader::Pointer oR = itk::TransformFileReader::New();

oR->SetFileName(oPathName.str().c_str());

oR->Update();





Output:

Matrix is orthogonal1

itk::Errir: Rigid3DTransform(00000003F30A10): Attempt to set a non-orthogonal rotation matrix



The written file:

#Insight Transform File V1.0

# Transform 0

Transform: Rigid3DTransform_double_3_3

Parameters: 0.998636 0.0229049 -0.0469187 -0.0228043 0.999736 0.00267739 0.0469676 -0.00160379 0.998895 -50.8226 -21.0202 -3.69866

FixedParameters: 0 0 0





Any suggestions?






Gerald Lodron
Human-Centered Image Processing
Machine Vision Applications
DIGITAL - Institute of Information and Communication Technologies
JOANNEUM RESEARCH Forschungsgesellschaft mbH
Steyrergasse 17, 8010 Graz, AUSTRIA

phone:  +43-316-876-1751         fax: +43-316-876-1751

web:    http://www.joanneum.at/digital
e-mail: gerald.lodron at joanneum.at<mailto:gerald.lodron at joanneum.at>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20110505/f17bad71/attachment.htm>