diff --git a/Examples/ChangeDetection/CMakeLists.txt b/Examples/ChangeDetection/CMakeLists.txt
index 351757ce0b6e0a9dab4ec07b0d89ab7192b80520..9abd8c7b2c563510fd16fe12a0675df38cdb6bfa 100644
--- a/Examples/ChangeDetection/CMakeLists.txt
+++ b/Examples/ChangeDetection/CMakeLists.txt
@@ -4,4 +4,20 @@ INCLUDE_REGULAR_EXPRESSION("^.*$")
ADD_EXECUTABLE(ChangeDetectionFrameworkExample ChangeDetectionFrameworkExample.cxx )
TARGET_LINK_LIBRARIES(ChangeDetectionFrameworkExample OTBCommon OTBIO ITKCommon ITKIO)
+ADD_EXECUTABLE(DiffChDet DiffChDet.cxx)
+TARGET_LINK_LIBRARIES(DiffChDet OTBIO OTBCommon ${CAI_LIBRARY} ${GDAL_LIBRARY} ITKCommon ITKIO)
+
+ADD_EXECUTABLE(RatioChDet RatioChDet.cxx)
+TARGET_LINK_LIBRARIES(RatioChDet OTBIO OTBCommon ${CAI_LIBRARY} ${GDAL_LIBRARY} ITKCommon ITKIO)
+
+ADD_EXECUTABLE(CorrelChDet CorrelChDet.cxx)
+TARGET_LINK_LIBRARIES(CorrelChDet OTBIO OTBCommon ${CAI_LIBRARY} ${GDAL_LIBRARY} ITKCommon ITKIO)
+
+ADD_EXECUTABLE(LHMIChDet LHMIChDet.cxx)
+TARGET_LINK_LIBRARIES(LHMIChDet OTBIO OTBCommon ${CAI_LIBRARY}
+${GDAL_LIBRARY} ITKCommon ITKIO ITKStatistics)
+
+ADD_EXECUTABLE(JHMIChDet JHMIChDet.cxx)
+TARGET_LINK_LIBRARIES(JHMIChDet OTBIO OTBCommon ${CAI_LIBRARY}
+${GDAL_LIBRARY} ITKCommon ITKIO ITKStatistics)
diff --git a/Examples/ChangeDetection/CorrelChDet.cxx b/Examples/ChangeDetection/CorrelChDet.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..69bb05c9d6cbaef5c3119704eea89be9a4afe0eb
--- /dev/null
+++ b/Examples/ChangeDetection/CorrelChDet.cxx
@@ -0,0 +1,106 @@
+/*=========================================================================
+
+ Programme : RatioChDet
+ Auteurs : CNES - J. INGLADA
+ Language : C++
+ Date : 8 juin 2006
+
+ Role : Détection de changements par corrélation locale
+
+
+=========================================================================*/
+
+#include "otbImageFileReader.h"
+#include "otbStreamingImageFileWriter.h"
+#include "itkImage.h"
+#include "itkShiftScaleImageFilter.h"
+#include "otbCorrelationChangeDetector.h"
+#include "otbCommandProgressUpdate.h"
+
+int main(int argc, char* argv[] )
+{
+
+ if( argc < 5 )
+ {
+ std::cerr << "Usage: " << std::endl;
+ std::cerr << argv[0] << " inputImageFile1 inputImageFile2 radius outputImageFile " << std::endl;
+ return -1;
+ }
+
+ // Define the dimension of the images
+ const unsigned int Dimension = 2;
+
+ // Declare the types of the images
+ typedef float InternalPixelType;
+ typedef unsigned char OutputPixelType;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType1;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType2;
+ typedef itk::Image<InternalPixelType, Dimension> ChangeImageType;
+ typedef itk::Image<OutputPixelType, Dimension> OutputImageType;
+
+ typedef otb::ImageFileReader< InputImageType1 > ReaderType1;
+ typedef otb::ImageFileReader< InputImageType2 > ReaderType2;
+ typedef otb::StreamingImageFileWriter< OutputImageType > WriterType;
+ typedef itk::ShiftScaleImageFilter< ChangeImageType,
+ OutputImageType > RescalerType;
+
+
+
+ // Declare the type for the filter
+ typedef otb::CorrelationChangeDetector<
+ InputImageType1,
+ InputImageType2,
+ ChangeImageType > FilterType;
+
+
+ ReaderType1::Pointer reader1 = ReaderType1::New();
+ ReaderType2::Pointer reader2 = ReaderType2::New();
+ WriterType::Pointer writer = WriterType::New();
+ FilterType::Pointer filter = FilterType::New();
+ RescalerType::Pointer rescaler = RescalerType::New();
+
+ const char * inputFilename1 = argv[1];
+ const char * inputFilename2 = argv[2];
+ const char * outputFilename = argv[4];
+
+ reader1->SetFileName( inputFilename1 );
+ reader2->SetFileName( inputFilename2 );
+ writer->SetFileName( outputFilename );
+
+ float scale = itk::NumericTraits< OutputPixelType >::max();
+ rescaler->SetScale( scale );
+
+ filter->SetInput1( reader1->GetOutput() );
+ filter->SetInput2( reader2->GetOutput() );
+ filter->SetRadius( atoi(argv[3]) );
+
+ rescaler->SetInput( filter->GetOutput() );
+ writer->SetInput( rescaler->GetOutput() );
+
+
+ typedef otb::CommandProgressUpdate<FilterType> CommandType;
+
+ CommandType::Pointer observer = CommandType::New();
+ filter->AddObserver(itk::ProgressEvent(), observer);
+
+
+
+ try
+ {
+ writer->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "ExceptionObject caught !" << std::endl;
+ std::cout << err << std::endl;
+ return -1;
+ }
+
+
+ return 0;
+
+}
+
+
+
+
diff --git a/Examples/ChangeDetection/DiffChDet.cxx b/Examples/ChangeDetection/DiffChDet.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..63837604defea65b96e5df074b55e8ed81b6b76b
--- /dev/null
+++ b/Examples/ChangeDetection/DiffChDet.cxx
@@ -0,0 +1,238 @@
+/*=========================================================================
+
+ Programme : DiffChDet
+ Auteurs : CNES - J. INGLADA
+ Language : C++
+ Date : 8 juin 2006
+
+ Role : Détection de changements par différence de moyennes locales
+
+
+=========================================================================*/
+
+#if defined(_MSC_VER)
+#pragma warning ( disable : 4786 )
+#endif
+
+#ifdef __BORLANDC__
+#define ITK_LEAN_AND_MEAN
+#endif
+
+
+// Software Guide : BeginCommandLineArgs
+// INPUTS: {SpotBefore.png}, {SpotAfter.png}
+// OUTPUTS: {DiffChDet.png}
+// 5
+// Software Guide : EndCommandLineArgs
+
+
+
+// Software Guide : BeginLatex
+// This example illustrates the class
+// \doxygen{otb::MeanDifferenceImageFilter} for detecting changes
+// between pairs of images. This filter computes the mean intensity in
+// the neighborhood of each pixel of the pair of images to be compared
+// and uses the difference of means as a change indicator. This
+// example will use the images shown in
+// figure ~\ref{fig:DIFFCHDETINIM}. These correspond to the near
+// infrared band of two Spot acquisitions before and during a flood.
+// \begin{figure}
+// \center
+// \includegraphics[width=0.35\textwidth]{SpotBefore.eps}
+// \includegraphics[width=0.35\textwidth]{SpotAfter.eps}
+// \itkcaption[Spot Images for Change Detection]{Images used for the
+// change detection. Left: Before the flood. Right: during the flood.}
+// \label{fig:DIFFCHDETINIM}
+// \end{figure}
+//
+// We start by including the corresponding header file.
+//
+// Software Guide : EndLatex
+
+// Software Guide : BeginCodeSnippet
+#include "otbMeanDifferenceImageFilter.h"
+// Software Guide : EndCodeSnippet
+
+#include "otbImageFileReader.h"
+#include "otbStreamingImageFileWriter.h"
+#include "itkImage.h"
+#include "itkShiftScaleImageFilter.h"
+
+#include "otbCommandProgressUpdate.h"
+
+
+
+int main(int argc, char* argv[] )
+{
+
+ if( argc < 5 )
+ {
+ std::cerr << "Usage: " << std::endl;
+ std::cerr << argv[0] << " inputImageFile1 inputImageFile2 outputImageFile radius" << std::endl;
+ return -1;
+ }
+
+ // Define the dimension of the images
+ const unsigned int Dimension = 2;
+
+ // Software Guide : BeginLatex
+ // We start by declaring the types for the two input images, the
+ // change image and the image to be stored in a file for visualization.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ typedef float InternalPixelType;
+ typedef unsigned char OutputPixelType;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType1;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType2;
+ typedef itk::Image<InternalPixelType, Dimension> ChangeImageType;
+ typedef itk::Image<OutputPixelType, Dimension> OutputImageType;
+
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // We can now declare the types for the readers. Since the images
+ // can be vey large, we will force the pipeline to use
+ // streaming. For this purpose, the file writer will be
+ // streamed. This is achieved by using the
+ // \doxygen{otb::StreamingImageFileWriter} class. Also, before
+ // saving the image to a file in, for instance, PNG format, we will
+ // rescale the results of the change detection in order to use all
+ // the output pixel type range of values.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ typedef otb::ImageFileReader< InputImageType1 > ReaderType1;
+ typedef otb::ImageFileReader< InputImageType2 > ReaderType2;
+ typedef otb::StreamingImageFileWriter< OutputImageType > WriterType;
+ typedef itk::ShiftScaleImageFilter< ChangeImageType,
+ OutputImageType > RescalerType;
+
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // The \doxygen{otb::MeanDifferenceImageFilter} is templated over
+ // the types of the two input images and the type of the generated change
+ // image.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ typedef otb::MeanDifferenceImageFilter<
+ InputImageType1,
+ InputImageType2,
+ ChangeImageType > FilterType;
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // The different elements of the pipeline can now be instantiated.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ ReaderType1::Pointer reader1 = ReaderType1::New();
+ ReaderType2::Pointer reader2 = ReaderType2::New();
+ WriterType::Pointer writer = WriterType::New();
+ FilterType::Pointer filter = FilterType::New();
+ RescalerType::Pointer rescaler = RescalerType::New();
+ // Software Guide : EndCodeSnippet
+
+ const char * inputFilename1 = argv[1];
+ const char * inputFilename2 = argv[2];
+ const char * outputFilename = argv[3];
+
+ // Software Guide : BeginLatex
+ //
+ // We set the parametters of the different elements of the pipeline.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ reader1->SetFileName( inputFilename1 );
+ reader2->SetFileName( inputFilename2 );
+ writer->SetFileName( outputFilename );
+ rescaler->SetShift( itk::NumericTraits< InternalPixelType >::min());
+
+ float scale = itk::NumericTraits< OutputPixelType >::max()/
+ float(itk::NumericTraits< OutputPixelType >::max()+
+ itk::NumericTraits< InternalPixelType >::min());
+ rescaler->SetScale( scale );
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // The only parametter for this change detector is the radius of
+ // the window used for computing the mean of the intensities.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ filter->SetRadius( atoi(argv[4]) );
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // We build the pipeline by plugging all the elements together.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ filter->SetInput1( reader1->GetOutput() );
+ filter->SetInput2( reader2->GetOutput() );
+ rescaler->SetInput( filter->GetOutput() );
+ writer->SetInput( rescaler->GetOutput() );
+ // Software Guide : EndCodeSnippet
+
+ // Software Guide : BeginLatex
+ //
+ // Since the processing time of large images can be long, it is
+ // interesting to monitor the evolution of the computation. In
+ // order to do so, the change detectors can use the
+ // command/observer design pattern. This is easily done by
+ // attaching an observer to the filter.
+ //
+ // Software Guide : EndLatex
+
+ // Software Guide : BeginCodeSnippet
+ typedef otb::CommandProgressUpdate<FilterType> CommandType;
+
+ CommandType::Pointer observer = CommandType::New();
+ filter->AddObserver(itk::ProgressEvent(), observer);
+ // Software Guide : EndCodeSnippet
+
+
+ try
+ {
+ writer->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "ExceptionObject caught !" << std::endl;
+ std::cout << err << std::endl;
+ return -1;
+ }
+// Software Guide : BeginLatex
+// Figure \ref{fig:RESDIFFCHDET} shows the result of the change
+// detection by difference of local means.
+// \begin{figure}
+// \center
+// \includegraphics[width=0.35\textwidth]{DiffChDet.eps}
+// \itkcaption[Difference Change Detection Results]{Result of the mean
+// difference change detector}
+// \label{fig:RESDIFFCHDET}
+// \end{figure}
+// Software Guide : EndLatex
+
+
+ return 0;
+
+}
+
+
+
+
diff --git a/Examples/ChangeDetection/JHMIChDet.cxx b/Examples/ChangeDetection/JHMIChDet.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..04c5101e4d4ac174986b02778a190c15e8655ed0
--- /dev/null
+++ b/Examples/ChangeDetection/JHMIChDet.cxx
@@ -0,0 +1,115 @@
+/*=========================================================================
+
+ Programme : RatioChDet
+ Auteurs : CNES - J. INGLADA
+ Language : C++
+ Date : 8 juin 2006
+
+ Role : Détection de changements par information mutuelle semi-locale
+
+
+=========================================================================*/
+
+#include "otbImageFileReader.h"
+#include "otbStreamingImageFileWriter.h"
+#include "itkImage.h"
+#include "itkShiftScaleImageFilter.h"
+#include "otbJoinHistogramMIImageFilter.h"
+#include "otbCommandProgressUpdate.h"
+
+int main(int argc, char* argv[] )
+{
+
+ if( argc < 5 )
+ {
+ std::cerr << "Usage: " << std::endl;
+ std::cerr << argv[0] << " inputImageFile1 inputImageFile2 radius outputImageFile " << std::endl;
+ return -1;
+ }
+
+ // Define the dimension of the images
+ const unsigned int Dimension = 2;
+
+ // Declare the types of the images
+ typedef float InternalPixelType;
+ typedef unsigned char OutputPixelType;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType1;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType2;
+ typedef itk::Image<InternalPixelType, Dimension> ChangeImageType;
+ typedef itk::Image<OutputPixelType, Dimension> OutputImageType;
+
+ typedef otb::ImageFileReader< InputImageType1 > ReaderType1;
+ typedef otb::ImageFileReader< InputImageType2 > ReaderType2;
+ typedef otb::StreamingImageFileWriter< OutputImageType > WriterType;
+
+ typedef itk::ShiftScaleImageFilter< ChangeImageType,
+ OutputImageType > RescalerType;
+
+
+
+ // Declare the type for the filter
+ typedef otb::JoinHistogramMIImageFilter<
+ InputImageType1,
+ InputImageType2,
+ ChangeImageType > FilterType;
+
+
+ ReaderType1::Pointer reader1 = ReaderType1::New();
+ ReaderType2::Pointer reader2 = ReaderType2::New();
+ WriterType::Pointer writer = WriterType::New();
+ FilterType::Pointer filter = FilterType::New();
+ RescalerType::Pointer rescaler = RescalerType::New();
+
+ const char * inputFilename1 = argv[1];
+ const char * inputFilename2 = argv[2];
+ const char * outputFilename = argv[4];
+
+ reader1->SetFileName( inputFilename1 );
+ reader2->SetFileName( inputFilename2 );
+ writer->SetFileName( outputFilename );
+
+
+ rescaler->SetShift( itk::NumericTraits< InternalPixelType >::min());
+
+ float scale = itk::NumericTraits< OutputPixelType >::max()/ 100.0;
+ /*float(itk::NumericTraits< OutputPixelType >::max()+
+ itk::NumericTraits< InternalPixelType >::min());*/
+ rescaler->SetScale( scale );
+
+
+ filter->SetInput1( reader1->GetOutput() );
+ filter->SetInput2( reader2->GetOutput() );
+ filter->SetRadius( atoi(argv[3]) );
+
+
+
+ rescaler->SetInput( filter->GetOutput() );
+ writer->SetInput( rescaler->GetOutput() );
+
+
+ typedef otb::CommandProgressUpdate<FilterType> CommandType;
+
+ CommandType::Pointer observer = CommandType::New();
+ filter->AddObserver(itk::ProgressEvent(), observer);
+
+
+
+ try
+ {
+ writer->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "ExceptionObject caught !" << std::endl;
+ std::cout << err << std::endl;
+ return -1;
+ }
+
+
+ return 0;
+
+}
+
+
+
+
diff --git a/Examples/ChangeDetection/LHMIChDet.cxx b/Examples/ChangeDetection/LHMIChDet.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..e21c28ac4c0b9a4f0f2120af564ef5a9568f9f64
--- /dev/null
+++ b/Examples/ChangeDetection/LHMIChDet.cxx
@@ -0,0 +1,106 @@
+/*=========================================================================
+
+ Programme : RatioChDet
+ Auteurs : CNES - J. INGLADA
+ Language : C++
+ Date : 8 juin 2006
+
+ Role : Détection de changements par information mutuelle locale
+
+
+=========================================================================*/
+
+#include "otbImageFileReader.h"
+#include "otbStreamingImageFileWriter.h"
+#include "itkImage.h"
+#include "itkShiftScaleImageFilter.h"
+#include "otbLHMIChangeDetector.h"
+#include "otbCommandProgressUpdate.h"
+
+int main(int argc, char* argv[] )
+{
+
+ if( argc < 5 )
+ {
+ std::cerr << "Usage: " << std::endl;
+ std::cerr << argv[0] << " inputImageFile1 inputImageFile2 radius outputImageFile " << std::endl;
+ return -1;
+ }
+
+ // Define the dimension of the images
+ const unsigned int Dimension = 2;
+
+ // Declare the types of the images
+ typedef float InternalPixelType;
+ typedef unsigned char OutputPixelType;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType1;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType2;
+ typedef itk::Image<InternalPixelType, Dimension> ChangeImageType;
+ typedef itk::Image<OutputPixelType, Dimension> OutputImageType;
+
+ typedef otb::ImageFileReader< InputImageType1 > ReaderType1;
+ typedef otb::ImageFileReader< InputImageType2 > ReaderType2;
+ typedef otb::StreamingImageFileWriter< OutputImageType > WriterType;
+ typedef itk::ShiftScaleImageFilter< ChangeImageType,
+ OutputImageType > RescalerType;
+
+
+
+ // Declare the type for the filter
+ typedef otb::LHMIChangeDetector<
+ InputImageType1,
+ InputImageType2,
+ ChangeImageType > FilterType;
+
+
+ ReaderType1::Pointer reader1 = ReaderType1::New();
+ ReaderType2::Pointer reader2 = ReaderType2::New();
+ WriterType::Pointer writer = WriterType::New();
+ FilterType::Pointer filter = FilterType::New();
+ RescalerType::Pointer rescaler = RescalerType::New();
+
+ const char * inputFilename1 = argv[1];
+ const char * inputFilename2 = argv[2];
+ const char * outputFilename = argv[4];
+
+ reader1->SetFileName( inputFilename1 );
+ reader2->SetFileName( inputFilename2 );
+ writer->SetFileName( outputFilename );
+
+ float scale = itk::NumericTraits< OutputPixelType >::max();
+ rescaler->SetScale( scale );
+
+ filter->SetInput1( reader1->GetOutput() );
+ filter->SetInput2( reader2->GetOutput() );
+ filter->SetRadius( atoi(argv[3]) );
+
+ rescaler->SetInput( filter->GetOutput() );
+ writer->SetInput( rescaler->GetOutput() );
+
+
+ typedef otb::CommandProgressUpdate<FilterType> CommandType;
+
+ CommandType::Pointer observer = CommandType::New();
+ filter->AddObserver(itk::ProgressEvent(), observer);
+
+
+
+ try
+ {
+ writer->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "ExceptionObject caught !" << std::endl;
+ std::cout << err << std::endl;
+ return -1;
+ }
+
+
+ return 0;
+
+}
+
+
+
+
diff --git a/Examples/ChangeDetection/RatioChDet.cxx b/Examples/ChangeDetection/RatioChDet.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..cda8d84d5d673229c5fc08e81b79cfb254185509
--- /dev/null
+++ b/Examples/ChangeDetection/RatioChDet.cxx
@@ -0,0 +1,106 @@
+/*=========================================================================
+
+ Programme : RatioChDet
+ Auteurs : CNES - J. INGLADA
+ Language : C++
+ Date : 8 juin 2006
+
+ Role : Détection de changements par rapport de moyennes locales
+
+
+=========================================================================*/
+
+#include "otbImageFileReader.h"
+#include "otbStreamingImageFileWriter.h"
+#include "itkImage.h"
+#include "itkShiftScaleImageFilter.h"
+#include "otbMeanRatioImageFilter.h"
+#include "otbCommandProgressUpdate.h"
+
+int main(int argc, char* argv[] )
+{
+
+ if( argc < 5 )
+ {
+ std::cerr << "Usage: " << std::endl;
+ std::cerr << argv[0] << " inputImageFile1 inputImageFile2 radius outputImageFile " << std::endl;
+ return -1;
+ }
+
+ // Define the dimension of the images
+ const unsigned int Dimension = 2;
+
+ // Declare the types of the images
+ typedef float InternalPixelType;
+ typedef unsigned char OutputPixelType;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType1;
+ typedef itk::Image<InternalPixelType, Dimension> InputImageType2;
+ typedef itk::Image<InternalPixelType, Dimension> ChangeImageType;
+ typedef itk::Image<OutputPixelType, Dimension> OutputImageType;
+
+ typedef otb::ImageFileReader< InputImageType1 > ReaderType1;
+ typedef otb::ImageFileReader< InputImageType2 > ReaderType2;
+ typedef otb::StreamingImageFileWriter< OutputImageType > WriterType;
+ typedef itk::ShiftScaleImageFilter< ChangeImageType,
+ OutputImageType > RescalerType;
+
+
+
+ // Declare the type for the filter
+ typedef otb::MeanRatioImageFilter<
+ InputImageType1,
+ InputImageType2,
+ ChangeImageType > FilterType;
+
+
+ ReaderType1::Pointer reader1 = ReaderType1::New();
+ ReaderType2::Pointer reader2 = ReaderType2::New();
+ WriterType::Pointer writer = WriterType::New();
+ FilterType::Pointer filter = FilterType::New();
+ RescalerType::Pointer rescaler = RescalerType::New();
+
+ const char * inputFilename1 = argv[1];
+ const char * inputFilename2 = argv[2];
+ const char * outputFilename = argv[4];
+
+ reader1->SetFileName( inputFilename1 );
+ reader2->SetFileName( inputFilename2 );
+ writer->SetFileName( outputFilename );
+
+ float scale = itk::NumericTraits< OutputPixelType >::max();
+ rescaler->SetScale( scale );
+
+ filter->SetInput1( reader1->GetOutput() );
+ filter->SetInput2( reader2->GetOutput() );
+ filter->SetRadius( atoi(argv[3]) );
+
+ rescaler->SetInput( filter->GetOutput() );
+ writer->SetInput( rescaler->GetOutput() );
+
+
+ typedef otb::CommandProgressUpdate<FilterType> CommandType;
+
+ CommandType::Pointer observer = CommandType::New();
+ filter->AddObserver(itk::ProgressEvent(), observer);
+
+
+
+ try
+ {
+ writer->Update();
+ }
+ catch( itk::ExceptionObject & err )
+ {
+ std::cout << "ExceptionObject caught !" << std::endl;
+ std::cout << err << std::endl;
+ return -1;
+ }
+
+
+ return 0;
+
+}
+
+
+
+
diff --git a/Examples/Data/SpotAfter.png b/Examples/Data/SpotAfter.png
new file mode 100644
index 0000000000000000000000000000000000000000..ac5fd8408effa7bb236a00b0cd8e2ceffbb01d98
Binary files /dev/null and b/Examples/Data/SpotAfter.png differ
diff --git a/Examples/Data/SpotBefore.png b/Examples/Data/SpotBefore.png
new file mode 100644
index 0000000000000000000000000000000000000000..e5427961404b7cdb9bb6eb3367c9649a091de2ba
Binary files /dev/null and b/Examples/Data/SpotBefore.png differ