- //-------------------------------------------------------------------
- // Update with the number of dimensions and the pixeltype
- //-------------------------------------------------------------------
- template<class args_info_type>
- template <unsigned int Dimension, class PixelType>
- void
- AverageTemporalDimensionGenericFilter<args_info_type>::UpdateWithDimAndPixelType()
- {
-
- // ImageTypes
- typedef itk::Image<PixelType, Dimension> InputImageType;
- typedef itk::Image<PixelType, Dimension-1> OutputImageType;
- typename InputImageType::Pointer input;
-
- if (m_ArgsInfo.input_given ==1 )
- {
- // Read the input
- typedef itk::ImageFileReader<InputImageType> InputReaderType;
- typename InputReaderType::Pointer reader = InputReaderType::New();
- reader->SetFileName( m_InputFileName);
- reader->Update();
- input= reader->GetOutput();
- }
-
- else
- {
- // Read and join multiple inputs
- if (m_Verbose) std::cout<<m_ArgsInfo.input_given<<" inputs given..."<<std::endl;
- std::vector<std::string> filenames;
- for(unsigned int i=0; i<m_ArgsInfo.input_given;i++)
- {
- if (m_Verbose) std::cout<<m_ArgsInfo.input_arg[i]<<std::endl;
- filenames.push_back(m_ArgsInfo.input_arg[i]);
- }
-
- typedef itk::ImageSeriesReader<InputImageType> ImageReaderType;
- typename ImageReaderType::Pointer reader= ImageReaderType::New();
- reader->SetFileNames(filenames);
- reader->Update();
- input =reader->GetOutput();
- }
-
-
- // Output properties
- typename OutputImageType::RegionType region;
- typename OutputImageType::RegionType::SizeType size;
- typename OutputImageType::IndexType index;
- typename OutputImageType::SpacingType spacing;
- typename OutputImageType::PointType origin;
- typename InputImageType::RegionType region4D=input->GetLargestPossibleRegion();
- typename InputImageType::RegionType::SizeType size4D=region4D.GetSize();
- typename InputImageType::IndexType index4D=region4D.GetIndex();
- typename InputImageType::SpacingType spacing4D=input->GetSpacing();
- typename InputImageType::PointType origin4D=input->GetOrigin();
-
- for (unsigned int i=0; i< Dimension-1; i++)
- {
- size[i]=size4D[i];
- index[i]=index4D[i];
- spacing[i]=spacing4D[i];
- origin[i]=origin4D[i];
- }
- region.SetSize(size);
- region.SetIndex(index);
- typename OutputImageType::Pointer output= OutputImageType::New();
- output->SetRegions(region);
- output->SetSpacing(spacing);
- output->SetOrigin(origin);
- output->Allocate();
-
-
- // Region iterators
- typedef itk::ImageRegionIterator<InputImageType> IteratorType;
- std::vector<IteratorType> iterators(size4D[Dimension-1]);
- for (unsigned int i=0; i< size4D[Dimension-1]; i++)
- {
- typename InputImageType::RegionType regionIt=region4D;
- typename InputImageType::RegionType::SizeType sizeIt=regionIt.GetSize();
- sizeIt[Dimension-1]=1;
- regionIt.SetSize(sizeIt);
- typename InputImageType::IndexType indexIt=regionIt.GetIndex();
- indexIt[Dimension-1]=i;
- regionIt.SetIndex(indexIt);
- iterators[i]=IteratorType(input, regionIt);
- }
-
- typedef itk::ImageRegionIterator<OutputImageType> OutputIteratorType;
- OutputIteratorType avIt(output, output->GetLargestPossibleRegion());
-
- // Average
- PixelType vector;
- PixelType zeroVector=itk::NumericTraits<PixelType>::Zero;
- //zeroVector.Fill(0.0);
- while (!(iterators[0]).IsAtEnd())
- {
- vector=zeroVector;
- for (unsigned int i=0; i<size4D[Dimension-1]; i++)
- {
- vector+=iterators[i].Get();
- ++(iterators[i]);
- }
- vector/=size4D[Dimension-1];
- avIt.Set(vector);
- ++avIt;
- }
-
- // Output
- typedef itk::ImageFileWriter<OutputImageType> WriterType;
- typename WriterType::Pointer writer = WriterType::New();
- writer->SetFileName(m_ArgsInfo.output_arg);
- writer->SetInput(output);
- writer->Update();
+ // Output properties
+ typename OutputImageType::RegionType region;
+ typename OutputImageType::RegionType::SizeType size;
+ typename OutputImageType::IndexType index;
+ typename OutputImageType::SpacingType spacing;
+ typename OutputImageType::PointType origin;
+ typename InputImageType::RegionType region4D=input->GetLargestPossibleRegion();
+ typename InputImageType::RegionType::SizeType size4D=region4D.GetSize();
+ typename InputImageType::IndexType index4D=region4D.GetIndex();
+ typename InputImageType::SpacingType spacing4D=input->GetSpacing();
+ typename InputImageType::PointType origin4D=input->GetOrigin();
+
+ for (unsigned int i=0; i< Dimension-1; i++) {
+ size[i]=size4D[i];
+ index[i]=index4D[i];
+ spacing[i]=spacing4D[i];
+ origin[i]=origin4D[i];
+ }
+ region.SetSize(size);
+ region.SetIndex(index);
+ typename OutputImageType::Pointer output= OutputImageType::New();
+ output->SetRegions(region);
+ output->SetSpacing(spacing);
+ output->SetOrigin(origin);
+ output->Allocate();
+
+
+ // Region iterators
+ typedef itk::ImageRegionIterator<InputImageType> IteratorType;
+ std::vector<IteratorType> iterators(size4D[Dimension-1]);
+ for (unsigned int i=0; i< size4D[Dimension-1]; i++) {
+ typename InputImageType::RegionType regionIt=region4D;
+ typename InputImageType::RegionType::SizeType sizeIt=regionIt.GetSize();
+ sizeIt[Dimension-1]=1;
+ regionIt.SetSize(sizeIt);
+ typename InputImageType::IndexType indexIt=regionIt.GetIndex();
+ indexIt[Dimension-1]=i;
+ regionIt.SetIndex(indexIt);
+ iterators[i]=IteratorType(input, regionIt);
+ }