X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=vv%2FvvImageReader.cxx;h=1f33fc72b2abe997ff119adf7f90b156edc7b0e7;hb=34a775104b11d2cab1fa48db1c6e254591b0fabe;hp=c28b2e0a4ca4761e4bb1d0acba3e86f69c08bbb8;hpb=934a64f5b51382ad6bbc78d9c8734691e599cf09;p=clitk.git
diff --git a/vv/vvImageReader.cxx b/vv/vvImageReader.cxx
index c28b2e0..1f33fc7 100644
--- a/vv/vvImageReader.cxx
+++ b/vv/vvImageReader.cxx
@@ -1,116 +1,178 @@
/*=========================================================================
+ Program: vv http://www.creatis.insa-lyon.fr/rio/vv
- Program: vv
- Language: C++
- Author : Pierre Seroul (pierre.seroul@gmail.com)
+ Authors belong to:
+ - University of LYON http://www.universite-lyon.fr/
+ - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
+ - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
-Copyright (C) 2008
-Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
-CREATIS-LRMN http://www.creatis.insa-lyon.fr
+ This software is distributed WITHOUT ANY WARRANTY; without even
+ the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ PURPOSE. See the copyright notices for more information.
-This program is free software: you can redistribute it and/or modify
-it under the terms of the GNU General Public License as published by
-the Free Software Foundation, version 3 of the License.
+ It is distributed under dual licence
-This program is distributed in the hope that it will be useful,
-but WITHOUT ANY WARRANTY; without even the implied warranty of
-MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-GNU General Public License for more details.
-
-You should have received a copy of the GNU General Public License
-along with this program. If not, see .
-
-=========================================================================*/
-#ifndef vvImageReader_CXX
-#define vvImageReader_CXX
+ - BSD See included LICENSE.txt file
+ - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+ ======================================================================-====*/
+#ifndef VVIMAGEREADER_CXX
+#define VVIMAGEREADER_CXX
#include
#include
#include "vvImageReader.h"
#include "vvImageReader.txx"
+#include "clitkTransformUtilities.h"
-
-//====================================================================
+//------------------------------------------------------------------------------
vvImageReader::vvImageReader()
{
- mImage = NULL;
- mInputFilenames.resize(0);
- mLastError = "";
- mType = UNDEFINEDIMAGETYPE;
+ mImage = NULL;
+ mInputFilenames.resize(0);
+ mLastError = "";
+ mType = UNDEFINEDIMAGETYPE;
+ mSlice = 0;
}
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
vvImageReader::~vvImageReader() { }
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
void vvImageReader::Update()
{
- Update(mType);
+ Update(mType);
}
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
void vvImageReader::Update(LoadedImageType type)
{
- itk::ImageIOBase::Pointer reader = itk::ImageIOFactory::CreateImageIO(mInputFilenames[0].c_str(), itk::ImageIOFactory::ReadMode);
- if (!reader) {
- mLastError="Unable to read file.";
- }
- else {
- reader->SetFileName(mInputFilenames[0]);
- reader->ReadImageInformation();
- if (mInputFilenames.size() > 1)
- Update(reader->GetNumberOfDimensions()+1,reader->GetComponentTypeAsString(reader->GetComponentType()),type);
- else
- Update(reader->GetNumberOfDimensions(),reader->GetComponentTypeAsString(reader->GetComponentType()),type);
- }
+ itk::ImageIOBase::Pointer reader = itk::ImageIOFactory::CreateImageIO(mInputFilenames[0].c_str(), itk::ImageIOFactory::ReadMode);
+ if (!reader) {
+ mLastError="Unable to read file.";
+ } else {
+ reader->SetFileName(mInputFilenames[0]);
+ reader->ReadImageInformation();
+ if (mInputFilenames.size() > 1)
+ Update(reader->GetNumberOfDimensions()+1,reader->GetComponentTypeAsString(reader->GetComponentType()),type);
+ else
+ Update(reader->GetNumberOfDimensions(),reader->GetComponentTypeAsString(reader->GetComponentType()),type);
+ }
}
+//------------------------------------------------------------------------------
-//====================================================================
-void vvImageReader::Update(int dim,std::string inputPixelType, LoadedImageType type) {
- //CALL_FOR_ALL_DIMS(dim,UpdateWithDim,inputPixelType);
- mType = type;
- mDim = dim;
- mInputPixelType=inputPixelType;
- this->start(); //Start heavy read operation in a separate thread
- while (this->isRunning())
- {
- qApp->processEvents();
- this->wait(50);
- }
+
+//------------------------------------------------------------------------------
+void vvImageReader::Update(int dim,std::string inputPixelType, LoadedImageType type)
+{
+ //CALL_FOR_ALL_DIMS(dim,UpdateWithDim,inputPixelType);
+ mType = type;
+ mDim = dim;
+ mInputPixelType=inputPixelType;
+ this->start(); //Start heavy read operation in a separate thread
+ while (this->isRunning()) {
+ qApp->processEvents();
+ this->wait(50);
+ }
}
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
void vvImageReader::run()
{
- switch(mDim)
- {
- case 2:
- UpdateWithDim<2>(mInputPixelType);
- break;;
- case 3:
- UpdateWithDim<3>(mInputPixelType);
- break;;
- case 4:
- UpdateWithDim<4>(mInputPixelType);
- break;;
- default:
- std::cerr << "dimension unknown in Update ! " << std::endl;
- }
+ switch(mDim) {
+ case 2:
+ UpdateWithDim<2>(mInputPixelType);
+ break;;
+ case 3:
+ UpdateWithDim<3>(mInputPixelType);
+ break;;
+ case 4:
+ UpdateWithDim<4>(mInputPixelType);
+ break;;
+ default:
+ std::cerr << "dimension unknown in Update ! " << std::endl;
+ }
}
-//====================================================================
+//------------------------------------------------------------------------------
+
-//====================================================================
-/*void vvImageReader::Extract(int dim, std::string inputPixelType, int slice) {
- CALL_FOR_ALL_DIMS(dim, ExtractWithDim, inputPixelType, slice);
-}*/
-//====================================================================
+//------------------------------------------------------------------------------
void vvImageReader::SetInputFilename(const std::string & filename)
{
- mInputFilenames.resize(0);
- mInputFilenames.push_back(filename);
+ mInputFilenames.resize(0);
+ mInputFilenames.push_back(filename);
}
+//------------------------------------------------------------------------------
-//====================================================================
-void vvImageReader::SetInputFilenames(const std::vector & filenames) {
- mInputFilenames = filenames;
+
+//------------------------------------------------------------------------------
+void vvImageReader::SetInputFilenames(const std::vector & filenames)
+{
+ mInputFilenames = filenames;
+}
+//------------------------------------------------------------------------------
+
+
+//------------------------------------------------------------------------------
+//Read transformation in NKI format (Xdr, transposed, cm)
+void vvImageReader::ReadNkiImageTransform()
+{
+ bool bRead=false;
+ std::string filename = mInputFilenames[0]+".MACHINEORIENTATION";
+ if(itksys::SystemTools::FileExists(filename.c_str())){
+ typedef itk::ImageFileReader< itk::Image< double, 2 > > MatrixReaderType;
+ MatrixReaderType::Pointer readerTransfo = MatrixReaderType::New();
+ readerTransfo->SetFileName(filename);
+ try {
+ bRead = true;
+ readerTransfo->Update();
+ } catch( itk::ExceptionObject & err ) {
+ bRead=false;
+ std::cerr << "Cannot read " << filename << std::endl
+ << "The error is: " << err << std::endl;
+ }
+
+ if (bRead) {
+ //Transpose matrix (NKI format)
+ for(int j=0; j<4; j++)
+ for(int i=0; i<4; i++)
+ mImage->GetTransform()->GetMatrix()->SetElement(j,i,readerTransfo->GetOutput()->GetBufferPointer()[4*i+j]);
+
+ //From cm to mm
+ for(int i=0; i<3; i++)
+ mImage->GetTransform()->GetMatrix()->SetElement(i,3,10*mImage->GetTransform()->GetMatrix()->GetElement(i,3));
+
+ mImage->GetTransform()->Inverse();
+ mImage->UpdateReslice();
+ }
+ }
}
-//====================================================================
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
+//Read transformation in ASCII format
+void vvImageReader::ReadMatImageTransform()
+{
+ std::string filename(mInputFilenames[0]+".mat");
+ std::ifstream f(filename.c_str());
+ if(f.is_open()) {
+ f.close();
+
+ itk::Matrix itkMat = clitk::ReadMatrix3D(filename);
+ for(int j=0; j<4; j++)
+ for(int i=0; i<4; i++)
+ mImage->GetTransform()->GetMatrix()->SetElement(j,i,itkMat[j][i]);
+ mImage->UpdateReslice();
+ }
+}
+//------------------------------------------------------------------------------
#endif