TEST: Adding test for ChangeNoDataFilter. Test is standalone and does not...

TEST: Adding test for ChangeNoDataFilter. Test is standalone and does not require ulages or baselines from OTB-Data

Modules/Filtering/ImageManipulation/test/CMakeLists.txt

@@ -70,6 +70,7 @@ otbSpectralAngleDistanceImageFilter.cxx
 otbFunctionWithNeighborhoodToImageFilterNew.cxx
 otbEuclideanDistanceMetricWithMissingValue.cxx
 otbEuclideanDistanceMetricWithMissingValueNew.cxx
+otbChangeNoDataValueFilter.cxx
 )
 
 add_executable(otbImageManipulationTestDriver ${OTBImageManipulationTests})
@@ -669,3 +670,6 @@ otb_add_test(NAME bfTvEuclideanDistanceMetricWithMissingValue COMMAND otbImageMa
   otbEuclideanDistanceMetricWithMissingValue)
 otb_add_test(NAME bfTuEuclideanDistanceMetricWithMissingValueNew COMMAND otbImageManipulationTestDriver
   otbEuclideanDistanceMetricWithMissingValueNew)
+
+otb_add_test(NAME filteringImageManipulationChangeNoDataValueFilter COMMAND otbImageManipulationTestDriver
+             otbChangeNoDataValueFilter)

Modules/Filtering/ImageManipulation/test/otbChangeNoDataValueFilter.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.
+
+
+     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 "otbImage.h"
+#include "otbChangeNoDataValueFilter.h"
+#include "itkImageRegionIterator.h"
+
+int otbChangeNoDataValueFilter(int itkNotUsed(argc),char * itkNotUsed(argv) [])
+{
+
+  // Build an image
+  typedef otb::Image<double> ImageType;
+  ImageType::Pointer img = ImageType::New();
+
+  ImageType::SizeType size;
+  size.Fill(20);
+
+  ImageType::RegionType region;
+  region.SetSize(size);
+
+  // Fill it with a default value
+  img->SetRegions(region);
+  img->Allocate();
+  img->FillBuffer(10);
+
+  // Write no-data flags to it
+  std::vector<bool> flags(1,true);
+  std::vector<double> values(1,-10.);
+  otb::WriteNoDataFlags(flags,values,img->GetMetaDataDictionary());
+
+  // Fill half of the pixels with no-data values
+  itk::ImageRegionIterator<ImageType> it(img,region);
+  unsigned int count = 0;
+  for(it.GoToBegin();!it.IsAtEnd();++it,++count)
+    {
+    if (count%2 == 0)
+      it.Set(-10.);
+    }
+
+  // Instanciate filter
+  typedef otb::ChangeNoDataValueFilter<ImageType,ImageType> FilterType;
+  FilterType::Pointer filter = FilterType::New();
+
+  std::vector<double> newValues(1,-20.);
+
+  filter->SetInput(img);
+  filter->SetNewNoDataValues(newValues);
+  filter->Update();
+
+  // Check output
+  it = itk::ImageRegionIterator<ImageType>(filter->GetOutput(),region);
+  count = 0;
+
+  bool failed = false;
+  
+  for(it.GoToBegin();!it.IsAtEnd();++it,++count)
+    {
+    if (count%2 == 0 and it.Get()!=-20.)
+      {
+      std::cerr<<"Pixel should have new no-data value"<<std::endl;
+      failed = true;
+      }
+    else if(count%2 == 1 and it.Get()!=10.)
+      {
+      std::cerr<<"Pixel value should be unchanged"<<std::endl;
+      failed = true;
+      }
+    }
+  
+  otb::ReadNoDataFlags(filter->GetOutput()->GetMetaDataDictionary(),flags,values);
+
+  if(flags.empty() || !flags[0])
+    {
+    std::cerr<<"Output no-data flag should be [1] "<<std::endl;
+    failed = true;
+    }
+
+  if(values.empty() || values[0]!=-20)
+    {
+    std::cerr<<"Output no-data value should be [-20.]"<<std::endl;
+    failed = true;
+    }
+
+  if(failed)
+    return EXIT_FAILURE;
+
+  return EXIT_SUCCESS;
+}

Modules/Filtering/ImageManipulation/test/otbImageManipulationTestDriver.cxx

@@ -76,4 +76,5 @@ void RegisterTests()
   REGISTER_TEST(otbFunctionWithNeighborhoodToImageFilterNew);
   REGISTER_TEST(otbEuclideanDistanceMetricWithMissingValue);
   REGISTER_TEST(otbEuclideanDistanceMetricWithMissingValueNew);
+  REGISTER_TEST(otbChangeNoDataValueFilter);
 }