vtkMitkThickSlicesFilterTest.cpp

Go to the documentation of this file.

/*=========================================================================

Program:   Medical Imaging & Interaction Toolkit
Language:  C++
Date:      $Date: 2007-09-26 18:50:26 +0200 (Mi, 26 Sep 2007) $
Version:   $Revision: 9585 $

Copyright (c) German Cancer Research Center, Division of Medical and
Biological Informatics. All rights reserved.
See MITKCopyright.txt or https://www.mitk.org/copyright.html for details.

This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE.  See the above copyright notices for more information.

=========================================================================*/

#include <vtkImageData.h>

#include "mitkTestingMacros.h"

#include "vtkMitkThickSlicesFilter.h"


vtkImageData *GenerateTestImageForTSFilter()
{
  // a 2x2x2 image
  short myData[] = 
  {
   
    234,234,
    123,565,
   
    -213,800,
    1000,-20
  };
  
  vtkImageData *i = vtkImageData::New();
  
  i->SetExtent(0,1,0,1,0,1);
  
  i->SetScalarTypeToShort();

  i->AllocateScalars();
  
  short *p = (short*)i->GetScalarPointer();
  
  memcpy(p,myData,2*2*2*sizeof(short));
  
  return i;
}

void CheckResultImageForTSFilter(vtkImageData *i)
{
  int *e=i->GetExtent();
  
  MITK_TEST_CONDITION_REQUIRED( e[0] == 0 && e[1] == 1 && e[2] == 0 && e[3] == 1 && e[4] == 0 && e[5] == 0 , "output image has correct extent" )
  MITK_TEST_CONDITION_REQUIRED( i->GetScalarType() == VTK_SHORT , "output image has correct scalar type" )
  
  short expectedResult[] = 
  {
   
    234,800,
    1000,565
  };

  short *d = (short*)i->GetScalarPointer();

  MITK_TEST_CONDITION_REQUIRED( memcmp(d,expectedResult,sizeof(expectedResult)) == 0 , "output image has correct content" )
}


int vtkMitkThickSlicesFilterTest(int /*argc*/, char* /*argv*/[])
{
  MITK_TEST_BEGIN("ThickSlicesFilter")

  vtkImageData *i,*o;
  
  i = GenerateTestImageForTSFilter();

  vtkMitkThickSlicesFilter *f = vtkMitkThickSlicesFilter::New();
  f->SetThickSliceMode( 0 ); // MIP
  f->SetInput( i );
  f->Update();
  o = f->GetOutput();
  
  CheckResultImageForTSFilter(o);

  //Delete vtk variable correctly
  i->Delete();
  f->Delete();
  
  MITK_TEST_END()
}