- wxString WxAnySimpleDlg::getInfoImage()
- {
- return infoimage;
- }
-
-
- // FCY : just to check but not needed, we hardly suppose that we load only same type and dim of images
- const std::type_info & WxAnySimpleDlg::getType(const std::string &i_name)
- {
- typedef itk::ImageIOBase::IOComponentType ScalarPixelType;
- itk::ImageIOBase::Pointer imageIO = itk::ImageIOFactory::CreateImageIO(i_name.c_str(), itk::ImageIOFactory::ReadMode);
- imageIO->SetFileName(i_name.c_str());
- imageIO->ReadImageInformation();
- return imageIO->GetComponentTypeInfo();//AsString( imageIO->GetComponentType());
- }
-
- const size_t WxAnySimpleDlg::getNumberOfDimensions(const std::string &i_name)
- {
- typedef itk::ImageIOBase::IOComponentType ScalarPixelType;
- itk::ImageIOBase::Pointer imageIO = itk::ImageIOFactory::CreateImageIO(i_name.c_str(), itk::ImageIOFactory::ReadMode);
- imageIO->SetFileName(i_name.c_str());
- imageIO->ReadImageInformation();
- return imageIO->GetNumberOfDimensions();
- }
-
- void WxAnySimpleDlg::readImg(const std::string &i_name)
- {
- size_t dims = getNumberOfDimensions(i_name);
- const std::type_info *type= &getType(i_name);
- switch(dims)
- {
- case 3:
- if (getType(i_name) == typeid(unsigned char))
- {
- typedef itk::Image<unsigned char, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed char))
- {
- typedef itk::Image<signed char, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(short))
- {
- typedef itk::Image<short, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned short))
- {
- typedef itk::Image<unsigned short, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned int))
- {
- typedef itk::Image<unsigned int, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed int))
- {
- typedef itk::Image<signed int, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned long))
- {
- typedef itk::Image<unsigned long, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed long))
- {
- typedef itk::Image<signed long, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(float))
- {
- typedef itk::Image<float, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(double))
- {
- typedef itk::Image<double, 3> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- }
- else
- {
- //????FCY, so what the type????
- }
- break;
- //assume that we have only one 4D file
- case 4:
- if (getType(i_name) == typeid(unsigned char))
- {
- typedef itk::Image<unsigned char, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed char))
- {
- typedef itk::Image<signed char, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned short))
- {
- typedef itk::Image<unsigned short, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if(getType(i_name) == typeid(short))
- {
- typedef itk::Image<short, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- try
- {
- reader->Update();
- }
- catch( itk::ExceptionObject & err )
- {
- std::cout << "Caught an exception reading" << i_name << ": " << std::endl;
- std::cout << err << " " << __FILE__ << " " << __LINE__ << std::endl;
- throw err;
- }
- catch(...)
- {
- std::cout << "Error while reading image " << i_name << std::endl;
- throw;
- }
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned short))
- {
- typedef itk::Image<unsigned short, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned int))
- {
- typedef itk::Image<unsigned int, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed int))
- {
- typedef itk::Image<signed int, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(unsigned long))
- {
- typedef itk::Image<unsigned long, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(signed long))
- {
- typedef itk::Image<signed long, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(float))
- {
- typedef itk::Image<float, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else if (getType(i_name) == typeid(double))
- {
- typedef itk::Image<double, 4> TImage;
- typedef itk::ImageFileReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetFileName( i_name );
- reader->Update();
- reader->GetOutput()->Register();
- m_AnyImages.push_back(reader->GetOutput());
- split4Din3Dvtk<TImage>(reader->GetOutput());
- }
- else
- {
- //????FCY, so what the type????
- }
- break;
- }
-
- // Test to know if we have to add new type
- bool btest = false;
- if( m_AnyDims.size() != 0)
- {
- btest = true;
- }
- else
- {
- m_AnyDims.push_back(dims);
- m_AnyType.push_back(&getType(i_name));
- }
- if (btest)
- {
- if(m_AnyDims.front() != dims || m_AnyType.front()->before(getType(i_name)))
- {
- m_AnyDims.push_back(dims);
- m_AnyType.push_back(&getType(i_name));
- }
- else{}
- }
- else {}
- }
-
-
-
-
- void WxAnySimpleDlg::readDicomImg(const std::vector<std::string> &i_names)
- {
- typedef itk::Image<short,3> TImage;
- typedef itk::GDCMImageIO GDCMType;
- typedef itk::DICOMSeriesFileNames dicnames;
- GDCMType::Pointer gdcmIO = GDCMType::New();
- dicnames::Pointer generator = dicnames::New();
-
-
- typedef itk::ImageSeriesReader<TImage> ReaderType;
- ReaderType::Pointer reader = ReaderType::New();
- reader->SetImageIO(gdcmIO);
- reader->SetFileNames(i_names);
- // reader->SetFileName( i_name );
- try
- {
- reader->Update();
- }
- catch( itk::ExceptionObject & err )
- {
- // std::cout << "Caught an exception reading" << i_name << ": " << std::endl;
- std::cout << err << " " << __FILE__ << " " << __LINE__ << std::endl;
- throw err;
- }
- catch(...)
- {
- //std::cout << "Error while reading image " << i_name << std::endl;
- throw;
- }
- //m_Iresults.push_back(reader->GetOutput());
+ wxString WxAnySimpleDlg::getInfoImage()
+ {
+ return infoimage;
+ }
+
+
+ // FCY : just to check but not needed, we hardly suppose that we load only same type and dim of images
+ const std::type_info & WxAnySimpleDlg::getType(const std::string &i_name)
+ {
+ typedef itk::ImageIOBase::IOComponentType ScalarPixelType;
+ itk::ImageIOBase::Pointer imageIO = itk::ImageIOFactory::CreateImageIO(i_name.c_str(), itk::ImageIOFactory::ReadMode);
+ imageIO->SetFileName(i_name.c_str());
+ imageIO->ReadImageInformation();
+ return imageIO->GetComponentTypeInfo();//AsString( imageIO->GetComponentType());
+ }
+
+ const size_t WxAnySimpleDlg::getNumberOfDimensions(const std::string &i_name)
+ {
+ typedef itk::ImageIOBase::IOComponentType ScalarPixelType;
+ itk::ImageIOBase::Pointer imageIO = itk::ImageIOFactory::CreateImageIO(i_name.c_str(), itk::ImageIOFactory::ReadMode);
+ imageIO->SetFileName(i_name.c_str());
+ imageIO->ReadImageInformation();
+ return imageIO->GetNumberOfDimensions();
+ }
+
+ void WxAnySimpleDlg::readImg(const std::string &i_name)
+ {
+ size_t dims = getNumberOfDimensions(i_name);
+ //const std::type_info *type= &getType(i_name);
+ switch(dims)
+ {
+ case 3:
+ if (getType(i_name) == typeid(unsigned char))
+ {
+ typedef itk::Image<unsigned char, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed char))
+ {
+ typedef itk::Image<signed char, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(short))
+ {
+ typedef itk::Image<short, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned short))
+ {
+ typedef itk::Image<unsigned short, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned int))
+ {
+ typedef itk::Image<unsigned int, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed int))
+ {
+ typedef itk::Image<signed int, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned long))
+ {
+ typedef itk::Image<unsigned long, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed long))
+ {
+ typedef itk::Image<signed long, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(float))
+ {
+ typedef itk::Image<float, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(double))
+ {
+ typedef itk::Image<double, 3> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ }
+ else
+ {
+ //????FCY, so what the type????
+ }
+ break;
+ //assume that we have only one 4D file
+ case 4:
+ if (getType(i_name) == typeid(unsigned char))
+ {
+ typedef itk::Image<unsigned char, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed char))
+ {
+ typedef itk::Image<signed char, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned short))
+ {
+ typedef itk::Image<unsigned short, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if(getType(i_name) == typeid(short))
+ {
+ typedef itk::Image<short, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ try
+ {
+ reader->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "Caught an exception reading" << i_name << ": " << std::endl;
+ std::cout << err << " " << __FILE__ << " " << __LINE__ << std::endl;
+ throw err;
+ }
+ catch(...)
+ {
+ std::cout << "Error while reading image " << i_name << std::endl;
+ throw;
+ }
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned short))
+ {
+ typedef itk::Image<unsigned short, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned int))
+ {
+ typedef itk::Image<unsigned int, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed int))
+ {
+ typedef itk::Image<signed int, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(unsigned long))
+ {
+ typedef itk::Image<unsigned long, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(signed long))
+ {
+ typedef itk::Image<signed long, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(float))
+ {
+ typedef itk::Image<float, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else if (getType(i_name) == typeid(double))
+ {
+ typedef itk::Image<double, 4> TImage;
+ typedef itk::ImageFileReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetFileName( i_name );
+ reader->Update();
+ reader->GetOutput()->Register();
+ m_AnyImages.push_back(reader->GetOutput());
+ split4Din3Dvtk<TImage>(reader->GetOutput());
+ }
+ else
+ {
+ //????FCY, so what the type????
+ }
+ break;
+ }
+
+ // Test to know if we have to add new type
+ bool btest = false;
+ if( m_AnyDims.size() != 0)
+ {
+ btest = true;
+ }
+ else
+ {
+ m_AnyDims.push_back(dims);
+ m_AnyType.push_back((std::type_info*)&getType(i_name)); //JPR cast
+ }
+ if (btest)
+ {
+ if(m_AnyDims.front() != dims || m_AnyType.front()->before(getType(i_name)))
+ {
+ m_AnyDims.push_back(dims);
+ m_AnyType.push_back((std::type_info*)&getType(i_name)); //JPR cast
+ }
+ else{}
+ }
+ else {}
+ }
+
+ void WxAnySimpleDlg::readDicomImg(const std::vector<std::string> &i_names)
+ {
+ typedef itk::Image<short,3> TImage;
+ typedef itk::GDCMImageIO GDCMType;
+ typedef itk::DICOMSeriesFileNames dicnames;
+ GDCMType::Pointer gdcmIO = GDCMType::New();
+ dicnames::Pointer generator = dicnames::New();
+
+
+ typedef itk::ImageSeriesReader<TImage> ReaderType;
+ ReaderType::Pointer reader = ReaderType::New();
+ reader->SetImageIO(gdcmIO);
+ reader->SetFileNames(i_names);
+ // reader->SetFileName( i_name );
+ try
+ {
+ reader->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ // std::cout << "Caught an exception reading" << i_name << ": " << std::endl;
+ std::cout << err << " " << __FILE__ << " " << __LINE__ << std::endl;
+ throw err;
+ }
+ catch(...)
+ {
+ //std::cout << "Error while reading image " << i_name << std::endl;
+ throw;
+ }
+ //m_Iresults.push_back(reader->GetOutput());