cmtkTransformChangeFromSpaceAffine.cxx

Go to the documentation of this file.
00001 /*
00002 //
00003 //  Copyright 1997-2009 Torsten Rohlfing
00004 //
00005 //  Copyright 2004-2010 SRI International
00006 //
00007 //  This file is part of the Computational Morphometry Toolkit.
00008 //
00009 //  http://www.nitrc.org/projects/cmtk/
00010 //
00011 //  The Computational Morphometry Toolkit is free software: you can
00012 //  redistribute it and/or modify it under the terms of the GNU General Public
00013 //  License as published by the Free Software Foundation, either version 3 of
00014 //  the License, or (at your option) any later version.
00015 //
00016 //  The Computational Morphometry Toolkit is distributed in the hope that it
00017 //  will be useful, but WITHOUT ANY WARRANTY; without even the implied
00018 //  warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019 //  GNU General Public License for more details.
00020 //
00021 //  You should have received a copy of the GNU General Public License along
00022 //  with the Computational Morphometry Toolkit.  If not, see
00023 //  <http://www.gnu.org/licenses/>.
00024 //
00025 //  $Revision: 2479 $
00026 //
00027 //  $LastChangedDate: 2010-10-20 15:35:37 -0700 (Wed, 20 Oct 2010) $
00028 //
00029 //  $LastChangedBy: torstenrohlfing $
00030 //
00031 */
00032 
00033 #include "cmtkTransformChangeFromSpaceAffine.h"
00034 
00035 cmtk::TransformChangeFromSpaceAffine
00036 ::TransformChangeFromSpaceAffine
00037 ( const AffineXform& xform, const UniformVolume& reference, const UniformVolume& floating, const char* forceSpace )
00038 {
00039   // adapt transformation to Slicer's image coordinate systems as defined in the Nrrd files we probably read
00040   UniformVolume::SmartPtr refVolumeOriginalSpace( reference.CloneGrid() );
00041   UniformVolume::SmartPtr fltVolumeOriginalSpace( floating.CloneGrid() );
00042   
00043   // first bring volumes back into their native coordinate space, or forced space if one was provided.
00044   if ( forceSpace )
00045     {
00046     refVolumeOriginalSpace->ChangeCoordinateSpace( forceSpace );
00047     fltVolumeOriginalSpace->ChangeCoordinateSpace( forceSpace );
00048     }
00049   else
00050     {
00051     refVolumeOriginalSpace->ChangeCoordinateSpace( reference.m_MetaInformation[META_SPACE_ORIGINAL] );
00052     fltVolumeOriginalSpace->ChangeCoordinateSpace( floating.m_MetaInformation[META_SPACE_ORIGINAL] );
00053     }
00054   
00055   // now determine image-to-physical transformations and concatenate these.
00056   AffineXform::MatrixType concatMatrix = refVolumeOriginalSpace->GetImageToPhysicalMatrix ();
00057   AffineXform::MatrixType fltMatrix = fltVolumeOriginalSpace->GetImageToPhysicalMatrix ();
00058   
00059   (concatMatrix *= xform.Matrix) *= fltMatrix.GetInverse();
00060   
00061   // create output transformation and write
00062   this->m_NewXform.SetMatrix( concatMatrix );
00063 }
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines