vv now builds even with all options turned off

[clitk.git] / itk / clitkExtractSliceFilter.txx

diff --git a/itk/clitkExtractSliceFilter.txx b/itk/clitkExtractSliceFilter.txx
index 40f7982b71e62445477f9b102f21f76a848f3f44..594fa106f4ea2397ddab19b1cad25db2cc301642 100644 (file)
--- a/itk/clitkExtractSliceFilter.txx
+++ b/itk/clitkExtractSliceFilter.txx
@@ -19,13 +19,15 @@
 // clitk
 #include "clitkCommon.h"
 
+// itk
+#include <itkExtractImageFilter.h>
 
 //--------------------------------------------------------------------
 template <class ImageType>
 clitk::ExtractSliceFilter<ImageType>::
 ExtractSliceFilter():
   clitk::FilterBase(),
-  itk::ImageToImageFilter<ImageType, ImageType>()
+  Superclass()
 {
   this->SetNumberOfRequiredInputs(1);
   SetDirection(2);
@@ -37,35 +39,45 @@ ExtractSliceFilter():
 template <class ImageType>
 void 
 clitk::ExtractSliceFilter<ImageType>::
-SetInput(const ImageType * image) {
+SetInput(const ImageType * image) 
+{
   // Process object is not const-correct so the const casting is required.
   this->SetNthInput(0, const_cast<ImageType *>(image));
 }
 //--------------------------------------------------------------------
   
 
+//--------------------------------------------------------------------
+template <class ImageType>
+void
+clitk::ExtractSliceFilter<ImageType>::
+GetOutputSlices(std::vector<typename SliceType::Pointer> & o)
+{
+  DD("GetOutputSlices");
+  o.clear();
+  for(unsigned int i=0; i<this->GetNumberOfOutputs(); i++) {
+    o.push_back(this->GetOutput(i));
+    // writeImage<SliceType>(this->GetOutput(i), "extractB"+clitk::toString(i)+".mhd");
+  }
+}
+//--------------------------------------------------------------------
+  
+
 //--------------------------------------------------------------------
 template <class ImageType>
 void 
 clitk::ExtractSliceFilter<ImageType>::
-GenerateOutputInformation() { 
+GenerateOutputInformation() 
+{ 
   DD("GenerateOutputInformation");
   ImagePointer input = dynamic_cast<ImageType*>(itk::ProcessObject::GetInput(0));
-  //   ImagePointer outputImage = this->GetOutput(0);
-  //   outputImage->SetRegions(input->GetLargestPossibleRegion());
-  
-  output = this->GetOutput(0);
   
-  // create vector
-  typename SliceType::RegionType SliceRegionType;
-  typename SliceType::SizeType SliceSizeType;
-  typename SliceType::IndexType SliceIndexType;
-  // SliceRegionType region;
-
-  // create region
-  // loop ExtractImageFilter with region updated, push_back
-
-
+  // Create region to extract
+  m_region = input->GetLargestPossibleRegion();
+  m_size   = m_region.GetSize();
+  m_index  = m_region.GetIndex();
+  m_NumberOfSlices = m_region.GetSize()[GetDirection()];
+  DD(m_NumberOfSlices);
 }
 //--------------------------------------------------------------------
 
@@ -77,7 +89,7 @@ clitk::ExtractSliceFilter<ImageType>::
 GenerateInputRequestedRegion() {
   DD("GenerateInputRequestedRegion");
   // Call default
-  itk::ImageToImageFilter<ImageType, ImageType>::GenerateInputRequestedRegion();
+  Superclass::GenerateInputRequestedRegion();
   // Get input pointers and set requested region to common region
   ImagePointer input = dynamic_cast<ImageType*>(itk::ProcessObject::GetInput(0));
   input->SetRequestedRegion(input->GetLargestPossibleRegion());
@@ -91,15 +103,53 @@ clitk::ExtractSliceFilter<ImageType>::
 GenerateData() {
   DD("GenerateData");
 
+  //--------------------------------------------------------------------
   // Get input pointer
   input = dynamic_cast<ImageType*>(itk::ProcessObject::GetInput(0));
 
-
+  //--------------------------------------------------------------------
+  // Create region to extract in 3D (and 2D = not used)
+  m_size[GetDirection()] = 0;
+  m_region.SetSize(m_size);
+  int start = m_index[GetDirection()];
+  DD(start);
+  this->SetNumberOfOutputs(m_NumberOfSlices);
+  /*
+    typename SliceType::RegionType regionSlice;
+    typename SliceType::SizeType   sizeSlice;
+    typename SliceType::IndexType  indexSlice;
+    int j=0;
+    for(int i=0; i<ImageDimension; i++) {
+    if (i != GetDirection()) {
+    sizeSlice[j] = m_size[i];
+    indexSlice[j] = m_index[j];
+    j++;
+    }
+    }
+    regionSlice.SetSize(sizeSlice);
+    regionSlice.SetIndex(indexSlice);
+    DD(regionSlice);
+  */
 
   //--------------------------------------------------------------------
-  //--------------------------------------------------------------------  
-  // Final Step -> set output
-  //this->SetNthOutput(0, working_image);
+  // loop ExtractImageFilter with region updated, push_back
+  typedef itk::ExtractImageFilter<ImageType, SliceType> ExtractImageFilterType;
+  typename ExtractImageFilterType::Pointer extract = ExtractImageFilterType::New();
+  extract->SetInput(input);
+  for(int i=0; i<m_NumberOfSlices; i++) {
+    extract = ExtractImageFilterType::New();
+    extract->SetInput(input);
+    // DD(i);
+    m_index[GetDirection()] = start + i;
+    m_region.SetIndex(m_index);
+    // DD(m_region);
+    extract->SetExtractionRegion(m_region);
+    extract->Update();
+    // DD("set output");
+    SetNthOutput(i, extract->GetOutput());
+    // writeImage<SliceType>(extract->GetOutput(), "extractA"+clitk::toString(i)+".mhd");
+  }
+
   return;
 }
 //--------------------------------------------------------------------