ENH: add extract before the process, test timeout on some plateform

Testing/Fa/00134-feature_extraction.cxx

@@ -27,14 +27,17 @@
 #include "otbImageListToVectorImageFilter.h"
 #include "itkMeanImageFilter.h"
 #include "otbMultiToMonoChannelExtractROI.h"
+#include "otbMultiChannelExtractROI.h"
 #include "otbMultiChannelRAndNIRIndexImageFilter.h"
 #include "otbVegetationIndicesFunctor.h"
 
 int main(int argc, char* argv[])
 {
-  const char *       inputName         = argv[1];
-  const char *       outputName        = argv[2];
+  const char *       inputName   = argv[1];
+  const char *       outputName  = argv[2];
   const unsigned int radius      = atoi(argv[3]);
+  const unsigned int id          = atoi(argv[4]);
+  const unsigned int size        = atoi(argv[5]);
 
   typedef double                                PixelType;
   typedef otb::Image<PixelType, 2>              ImageType;
@@ -45,22 +48,35 @@ int main(int argc, char* argv[])
   typedef otb::ImageList<ImageType>                                                             ImageListType;
   typedef otb::ImageListToVectorImageFilter<ImageListType, VectorImageType>                     ListToImageFilterType;
   typedef otb::MultiToMonoChannelExtractROI<PixelType, PixelType>                               ExtractorFilterType;
+  typedef otb::MultiChannelExtractROI<PixelType, PixelType>                                     MultiChannelExtractorFilterType;
   typedef itk::MeanImageFilter<ImageType, ImageType>                                            MeanFilterType;
   typedef otb::Functor::NDVI<PixelType, PixelType, PixelType>                                   NDVIFunctorType;
   typedef otb::MultiChannelRAndNIRIndexImageFilter<VectorImageType, ImageType, NDVIFunctorType> NDVIFilterType;
 
-  ReaderType::Pointer            reader    = ReaderType::New();
-  WriterType::Pointer            writer    = WriterType::New();
-  ImageListType::Pointer         imListOut = ImageListType::New();
-  ExtractorFilterType::Pointer   extract   = ExtractorFilterType::New();
-  MeanFilterType::Pointer        meanner   = MeanFilterType::New();
-  NDVIFilterType::Pointer        ndvi      = NDVIFilterType::New();
-  ListToImageFilterType::Pointer caster    = ListToImageFilterType::New();
+  ReaderType::Pointer                      reader    = ReaderType::New();
+  WriterType::Pointer                      writer    = WriterType::New();
+  ImageListType::Pointer                   imListOut = ImageListType::New();
+  ExtractorFilterType::Pointer             extract   = ExtractorFilterType::New();
+  MultiChannelExtractorFilterType::Pointer MCExtract = MultiChannelExtractorFilterType::New();
+  MeanFilterType::Pointer                  meanner   = MeanFilterType::New();
+  NDVIFilterType::Pointer                  ndvi      = NDVIFilterType::New();
+  ListToImageFilterType::Pointer           caster    = ListToImageFilterType::New();
 
   reader->SetFileName(inputName);
   reader->GenerateOutputInformation();
 
-  extract->SetInput(reader->GetOutput());
+  ImageType::RegionType region;
+  ImageType::SizeType   extSize;
+  ImageType::IndexType  extId;
+  extSize.Fill(size);
+  extId.Fill(id);
+  region.SetSize(extSize);
+  region.SetIndex(extId);
+  MCExtract->SetInput(reader->GetOutput());
+  MCExtract->SetExtractionRegion(region);
+  MCExtract->UpdateOutputInformation();
+
+  extract->SetInput(MCExtract->GetOutput());
   extract->SetChannel(1);
   extract->UpdateOutputInformation();
 
@@ -70,7 +86,7 @@ int main(int argc, char* argv[])
   meanner->SetInput(extract->GetOutput());
   meanner->SetRadius(rad);
 
-  ndvi->SetInput(reader->GetOutput());
+  ndvi->SetInput(MCExtract->GetOutput());
 
   imListOut->PushBack(meanner->GetOutput());
   imListOut->PushBack(ndvi->GetOutput());

Testing/Fa/CMakeLists.txt

@@ -232,6 +232,8 @@ ADD_TEST(FA-00134-feature_extraction ${CXX_TEST_PATH}/00134-feature_extraction
 	${IMAGEDATA}/QUICKBIRD/TOULOUSE/000000128955_01_P001_MUL/02APR01105228-M1BS-000000128955_01_P001.TIF
         ${TEMP}/extracted_filter_134.tif
         5
+        500 # extract start index
+        300 # extract size
 	)
 ENDIF(OTB_DATA_USE_LARGEINPUT)