[ITK] [ITK-users] Transforming a point
Yaniv, Ziv Rafael (NIH/NLM/LHC) [C]
zivrafael.yaniv at nih.gov
Wed Apr 12 09:41:13 EDT 2017
Hello Swetha,
In general, the transformations computed using the ITK registration framework are from the fixed coordinate system to the moving coordinate system. Thus, the transformation should be applied as T(point_fixed) = point_moving.
I did not read your code so I may be missing something specific to your work, but this mapping is something that many newcomers to ITK seem to get confused with.
hope this helps
Ziv
From: Swetha Sharma <swetha.bsharma at gmail.com>
Date: Wednesday, April 12, 2017 at 4:11 AM
To: Insight-users <insight-users at itk.org>
Subject: [ITK] [ITK-users] Transforming a point
Hi,
I have a affine transformation matrix with a set of points as fixed landmarks and moving landmarks. I was trying to verify the transformation by using the transform point method in the transformation class,ideally i should get back the fixed points when i give the moving points as input , but I am not getting it as the output. The code which i used is as below:
typedef itk::AffineTransform<double,3> AffineTransformType;
AffineTransformType::Pointer transf;
transf = dynamic_cast<AffineTransformType*>( transfm.GetPointer() );
typedef itk::Similarity3DTransform<double> SimilarityTransformType;
SimilarityTransformType::Pointer transform = SimilarityTransformType::New();
typedef itk::LandmarkBasedTransformInitializer< SimilarityTransformType, InputImageType, InputImageType >
LandmarkBasedTransformInitializerType;
LandmarkBasedTransformInitializerType::Pointer landmarkBasedTransformInitializer =
LandmarkBasedTransformInitializerType::New();
typedef LandmarkBasedTransformInitializerType::LandmarkPointContainer LandmarkContainerType;
typedef LandmarkBasedTransformInitializerType::LandmarkPointType LandmarkPointType;
LandmarkContainerType movingLandmarks;
LandmarkPointType movingPoint;
double *ptrf_moving_point = ptf_moving_landmark_pointslist;
for(int iPointIndex=0;iPointIndex<iNumberofLandMarkPoints_moving;iPointIndex++){
movingPoint[0] = (*( ptrf_moving_point )) ;
movingPoint[1] = (*( ptrf_moving_point + 1));
movingPoint[2] = (*( ptrf_moving_point + 2));
movingLandmarks.push_back( movingPoint );
ptrf_moving_point = ptrf_moving_point + 3;
}
typedef itk::PointSet< double, 3 > PointSetType;
PointSetType::Pointer registeredPointSet = PointSetType::New();
typedef PointSetType::PointType PointType;
PointType registeredPoint;
typedef PointSetType::PointsContainer PointsContainer;
PointsContainer::Pointer registeredPointContainer = PointsContainer::New();
for(int i=0;i<movingLandmarks.size();i++){
registeredPoint = transf->TransformPoint( movingLandmarks[i] );
registeredPointContainer->InsertElement( i, registeredPoint );
}
registeredPointSet->SetPoints(registeredPointContainer);
//Write the registered pointset into file
std::ofstream registeredFile;
registeredFile.open( "registered_points.txt");
if( registeredFile.fail() )
{
return -1;
}
for( int i = 0; i<registeredPointSet->GetNumberOfPoints(); i++)
{
registeredPointSet->GetPoint( i, & registeredPoint );
for(int j = 0; j<3; j++)
{
registeredFile<<registeredPoint[j]<<"\t";
}
registeredFile<<"\n";
}
Is there a mistake in the code ?How do I transform the points?
-swetha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/community/attachments/20170412/c1e7dd11/attachment.html>
-------------- next part --------------
_____________________________________
Powered by www.kitware.com
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php
Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ
Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/insight-users
More information about the Community
mailing list