itktwooutputexampleimagefiltertest.cxx

DTMK软件开发包,此为开源软件,是一款很好的医学图像开发资源.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkTwoOutputExampleImageFilterTest.cxx,v $
  Language:  C++
  Date:      $Date: 2007-08-20 07:50:39 $
  Version:   $Revision: 1.13 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm 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.

=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include "itkImage.h"
#include "itkRandomImageSource.h"
#include "itkTwoOutputExampleImageFilter.h"
#include "itkCommand.h"
#include "itkTextOutput.h"


int itkTwoOutputExampleImageFilterTest(int, char* [] )
{
  // Comment the following if you want to use the itk text output window
  itk::OutputWindow::SetInstance( itk::TextOutput::New() );

  // Uncomment the following if you want to see each message independently
  // itk::OutputWindow::GetInstance()->PromptUserOn();

  typedef itk::Image<float,2> FloatImage2DType;

  itk::RandomImageSource<FloatImage2DType>::Pointer random;
  random = itk::RandomImageSource<FloatImage2DType>::New();
  random->SetMin(0.0);
  random->SetMax(1000.0);

  float spacing[2] = {0.7, 2.1};
  random->SetSpacing( spacing );
  float origin[2] = {15, 400};
  random->SetOrigin( origin );
    
  itk::OStringStream *os;

  // Test #1, filter goes out of scope
  itk::OutputWindow::GetInstance()->DisplayText( "Test #1: Filter goes out of scope -----------------" );
  {
  itk::TwoOutputExampleImageFilter<FloatImage2DType>::Pointer threshold;
  threshold = itk::TwoOutputExampleImageFilter<FloatImage2DType>::New();
  threshold->SetInput(random->GetOutput());
  threshold->Update();

  std::cout << "Input spacing: " << random->GetOutput()->GetSpacing()[0]
            << ", "
            << random->GetOutput()->GetSpacing()[1] << std::endl;
  std::cout << "Output spacing: " << threshold->GetOutput()->GetSpacing()[0]
            << ", "
            << threshold->GetOutput()->GetSpacing()[1] << std::endl;

  os = new itk::OStringStream();
  *os << "Filter: " << threshold.GetPointer();
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #0: " << threshold->GetOutput(0);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #1: " << threshold->GetOutput(1);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;

  threshold->GetOutput(0)->DebugOn();
  threshold->GetOutput(1)->DebugOn();
  threshold->DebugOn();

  itk::OutputWindow::GetInstance()->DisplayText( "Ending Test #1: filter goes out of scope" );
  itk::OutputWindow::GetInstance()->DisplayText( "End of Test #1 -----------------------------------" );
  }
  
  // Test #2, user keeps an extra handle to an output
  itk::OutputWindow::GetInstance()->DisplayText( "Test #2: User keeps an extra hold on an output  -----------------" );
  {
  FloatImage2DType::Pointer keep;

  itk::TwoOutputExampleImageFilter<FloatImage2DType>::Pointer threshold;
  threshold = itk::TwoOutputExampleImageFilter<FloatImage2DType>::New();
  threshold->SetInput(random->GetOutput());
  threshold->Update();

  os = new itk::OStringStream();
  *os << "Filter: " << threshold.GetPointer();
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #0: " << threshold->GetOutput(0);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #1: " << threshold->GetOutput(1);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;

  threshold->GetOutput(0)->DebugOn();
  threshold->GetOutput(1)->DebugOn();
  threshold->DebugOn();

  keep = threshold->GetOutput(0);

  itk::OutputWindow::GetInstance()->DisplayText( "End of Test #2: last handle on output 0 should go out of scope");
  }
  itk::OutputWindow::GetInstance()->DisplayText( "End of Test #2 -----------------------------------");

  // Test #3, user disconnects a data object from the pipeline
  itk::OutputWindow::GetInstance()->DisplayText( "Test #3: application disconnects a data object from the pipeline  -----------------" );
  {
  FloatImage2DType::Pointer keep;

  itk::TwoOutputExampleImageFilter<FloatImage2DType>::Pointer threshold;
  threshold = itk::TwoOutputExampleImageFilter<FloatImage2DType>::New();
  threshold->SetInput(random->GetOutput());
  threshold->Update();

  os = new itk::OStringStream();
  *os << "Filter: " << threshold.GetPointer();
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #0: " << threshold->GetOutput(0);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;
  os = new itk::OStringStream();
  *os << "Output #1: " << threshold->GetOutput(1);
  itk::OutputWindow::GetInstance()->DisplayText( os->str().c_str() );
  delete os;

  threshold->GetOutput(0)->DebugOn();
  threshold->GetOutput(1)->DebugOn();
  threshold->DebugOn();

  keep = threshold->GetOutput(0);
  keep->DisconnectPipeline();

  itk::OutputWindow::GetInstance()->DisplayText( "End of Test #3: last handle on output 0 should go out of scope");
  }
  itk::OutputWindow::GetInstance()->DisplayText( "End of Test #3 -----------------------------------");
  

  // Test itkSetMacro and itkGetMacro
  itk::TwoOutputExampleImageFilter<FloatImage2DType>::Pointer threshold;
  threshold = itk::TwoOutputExampleImageFilter<FloatImage2DType>::New();
  const float outsideValue = 0;
  threshold->SetOutsideValue(outsideValue);
  float outsideValue2 = threshold->GetOutsideValue();
  std::cout << "threshold->GetOutsideValue(): " << outsideValue2 << std::endl;
  threshold->GetInverseOutput();

  return EXIT_SUCCESS;
}