tools/clitkMatrixInverse.cxx

   1 /*=========================================================================
   2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
   3
   4   Authors belong to:
   5   - University of LYON              http://www.universite-lyon.fr/
   6   - Léon Bérard cancer center       http://www.centreleonberard.fr
   7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
   8
   9   This software is distributed WITHOUT ANY WARRANTY; without even
  10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  11   PURPOSE.  See the copyright notices for more information.
  12
  13   It is distributed under dual licence
  14
  15   - BSD        See included LICENSE.txt file
  16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  17 ===========================================================================**/
  18
  19 // clitk
  20 #include "clitkMatrixInverse_ggo.h"
  21 #include "clitkTransformUtilities.h"
  22 #include "clitkIO.h"
  23 #include "clitkMatrix.h"
  24
  25 //--------------------------------------------------------------------
  26 int main(int argc, char * argv[])
  27 {
  28   // Init command line
  29   GGO(clitkMatrixInverse, args_info);
  30   CLITK_INIT;
  31
  32   // Read matrix
  33   itk::Matrix<double, 4, 4> matrix;
  34   try {
  35     matrix = clitk::ReadMatrix3D(args_info.input_arg);
  36   }
  37   catch (itk::ExceptionObject & err) {
  38     std::cerr << "Error reading " << args_info.input_arg << std::endl;
  39     std::cerr << err.GetDescription() << std::endl;
  40     exit(-1);
  41   }
  42
  43   matrix = matrix.GetInverse();
  44
  45   // Print matrix
  46   std::ofstream os;
  47   clitk::openFileForWriting(os, args_info.output_arg);
  48   os << clitk::Get4x4MatrixDoubleAsString(matrix, 16);
  49   os.close();
  50
  51   return EXIT_SUCCESS;
  52 }// end main
  53
  54 //--------------------------------------------------------------------