BUG: Region merging: set number of components for clustered output image

Code/BasicFilters/otbLabelImageRegionMergingFilter.h

@@ -111,6 +111,8 @@ public:
 
 protected:
 
+   virtual void GenerateOutputInformation(void);
+
    virtual void GenerateData();
 
   /** Constructor */

Code/BasicFilters/otbLabelImageRegionMergingFilter.txx

@@ -131,6 +131,21 @@ LabelImageRegionMergingFilter<TInputLabelImage, TInputSpectralImage, TOutputLabe
   return static_cast<OutputClusteredImageType *>(this->itk::ProcessObject::GetOutput(1));
 }
 
+template <class TInputLabelImage, class TInputSpectralImage, class TOutputLabelImage, class TOutputClusteredImage>
+void
+LabelImageRegionMergingFilter<TInputLabelImage, TInputSpectralImage, TOutputLabelImage, TOutputClusteredImage>
+::GenerateOutputInformation()
+{
+  Superclass::GenerateOutputInformation();
+
+  unsigned int numberOfComponentsPerPixel = this->GetInputSpectralImage()->GetNumberOfComponentsPerPixel();
+
+  if(this->GetClusteredOutput())
+    {
+    this->GetClusteredOutput()->SetNumberOfComponentsPerPixel(numberOfComponentsPerPixel);
+    }
+
+}
 
 template <class TInputLabelImage, class TInputSpectralImage, class TOutputLabelImage, class TOutputClusteredImage>
 void