TEST: add validation test for wavelet

Testing/Code/MultiScale/CMakeLists.txt

@@ -360,6 +360,22 @@ ADD_TEST(msTvGeodesicMorphologyIterativeDecompositionImageFilter ${MULTISCALE_TE
 ADD_TEST(msTuWaveletTransformNew ${MULTISCALE_TESTS4}
          otbWaveletTransformNew)
 
+ADD_TEST(msTvWaveletTransform ${MULTISCALE_TESTS4}
+        --compare-n-images ${TOL} 4
+              ${BASELINE}/msTvWaveletTransform-0.tif
+              ${TEMP}/msTvWaveletTransform-0.tif
+              ${BASELINE}/msTvWaveletTransform-1.tif
+              ${TEMP}/msTvWaveletTransform-1.tif
+              ${BASELINE}/msTvWaveletTransform-2.tif
+              ${TEMP}/msTvWaveletTransform-2.tif
+              ${BASELINE}/msTvWaveletTransform-3.tif
+              ${TEMP}/msTvWaveletTransform-3.tif
+         otbWaveletTransform
+          -in ${INPUTDATA}/ROI_IKO_PAN_LesHalles.tif
+          -l 1
+          -out ${TEMP}/msTvWaveletTransform
+)
+
 # -------       Fichiers sources CXX -----------------------------------
 SET(BasicMultiScale_SRCS1
 otbMorphologicalPyramidResamplerNew.cxx

Testing/Code/MultiScale/otbWaveletTransform.cxx

+/*=========================================================================
+
+  Program:   ORFEO Toolbox
+  Language:  C++
+  Date:      $Date$
+  Version:   $Revision$
+
+
+  Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
+  See OTBCopyright.txt for details.
+
+  Copyright (c) Institut Telecom / Telecom Bretagne. All rights reserved.
+  See ITCopyright.txt 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 "otbImage.h"
+#include "otbImageFileReader.h"
+#include "otbImageFileWriter.h"
+#include "otbImageViewer.h"
+
+#include "otbStationaryFilterBank.h"
+#include "otbWaveletForwardTransform.h"
+
+#include "otbHaarOperator.h"
+#include "otb_9_7_Operator.h"
+
+#include "otbCommandLineArgumentParser.h"
+#include "otbCommandProgressUpdate.h"
+
+int otbWaveletTransform( int argc, char ** argv )
+{
+  typedef otb::CommandLineArgumentParser ParserType;
+  ParserType::Pointer parser = ParserType::New();
+  parser->AddInputImage() ;
+  parser->AddOption( "--Level", "Decomposition level through the low pass branch (def. 1)", "-l", 1, false );
+  parser->AddOption( "--Output", "Output filename prefix", "-out", 1, true );
+  typedef otb::CommandLineArgumentParseResult ParserResultType;
+  ParserResultType::Pointer  parseResult = ParserResultType::New();
+
+  try
+  {
+    parser->ParseCommandLine( argc, argv, parseResult );
+  }
+  catch( itk::ExceptionObject & err )
+  {
+    std::cerr << "Test foward wavelet transform \n";
+    std::string descriptionException = err.GetDescription();
+    if ( descriptionException.find("ParseCommandLine(): Help Parser") != std::string::npos )
+      return EXIT_SUCCESS;
+    if ( descriptionException.find("ParseCommandLine(): Version Parser") != std::string::npos )
+      return EXIT_SUCCESS;
+    return EXIT_FAILURE;
+  }
+
+  const char * inputFileName = parseResult->GetInputImage().c_str();
+  unsigned int level = parseResult->IsOptionPresent("--Level" ) ?
+      parseResult->GetParameterUInt("--Level") : 1;
+
+  const int Dimension = 2;
+  typedef double PixelType;
+  typedef otb::Image< PixelType, Dimension >  ImageType;
+  typedef otb::ImageFileReader< ImageType > ReaderType;
+
+  /* Reading */
+  ReaderType::Pointer reader = ReaderType::New();
+  reader->SetFileName( inputFileName );
+
+  /* Transformation */
+  //typedef otb::LowPassHaarOperator< PixelType, Dimension > LowPassOperator;
+  //typedef otb::HighPassHaarOperator< PixelType, Dimension > HighPassOperator;
+  typedef otb::LowPass_9_7_Operator< PixelType, Dimension > LowPassOperator;
+  typedef otb::HighPass_9_7_Operator< PixelType, Dimension > HighPassOperator;
+
+  typedef otb::StationaryFilterBank< ImageType, ImageType, LowPassOperator, HighPassOperator > WaveletFilterType;
+  typedef otb::WaveletForwardTransform< ImageType, ImageType, WaveletFilterType > FilterType;
+  FilterType::Pointer filter = FilterType::New();
+  filter->SetInput( reader->GetOutput() );
+  filter->SetNumberOfDecompositions( level );
+
+  /* Observer */
+  typedef otb::CommandProgressUpdate< FilterType > CommandType;
+  CommandType::Pointer observer = CommandType::New();
+  filter->AddObserver( itk::ProgressEvent(), observer );
+
+  filter->Update();
+
+
+  typedef otb::ImageViewer< PixelType > ViewerType;
+
+  std::string prefix = parseResult->GetParameterString("--Output");
+  for ( unsigned int i = 0; i < filter->GetNumberOfOutputs(); i++ )
+  {
+    std::stringstream filename;
+    filename << prefix;
+    filename << "-" << i << ".tif";
+
+    typedef otb::ImageFileWriter<ImageType> WriterType;
+    WriterType::Pointer writer = WriterType::New();
+    writer->SetFileName(filename.str());
+    writer->SetInput(filter->GetOutput(i));
+    writer->Update();
+
+  }
+
+  return EXIT_SUCCESS;
+}
+
+
+