FIX: Encapsulate output metadata in GenericRSResampleImageFilter

Modules/Core/Metadata/src/otbImageMetadata.cxx

@@ -481,7 +481,7 @@ void ImageMetadata::Merge(const ImageMetadata& imd)
 {
   ImageMetadataBase::Fuse(imd);
   
-  for (unsigned int i = 0; i < std::min(Bands.size(), imd.Bands.size()); i++)
+  for (unsigned int i = 0; i < std::min(Bands.size(), imd.Bands.size()); ++i)
   {
     Bands[i].Fuse(imd.Bands[i]);
   }

Modules/Filtering/Projection/include/otbGenericRSResampleImageFilter.hxx

@@ -99,8 +99,10 @@ void GenericRSResampleImageFilter<TInputImage, TOutputImage>::GenerateOutputInfo
   m_Resampler->UpdateOutputInformation();
   this->GraftOutput(m_Resampler->GetOutput());
 
-  // Encapsulate output projRef
-  this->GetOutput()->m_Imd.Add(MDGeom::ProjectionProj, this->GetOutputProjectionRef());
+  // Encapsulate output projRef and metadata
+  if (this->GetOutputImageMetadata() != nullptr)
+    this->GetOutput()->m_Imd.Merge(*(this->GetOutputImageMetadata()));
+  this->GetOutput()->m_Imd.Add(MDGeom::ProjectionWKT, this->GetOutputProjectionRef());
 }
 
 /**
@@ -119,8 +121,8 @@ void GenericRSResampleImageFilter<TInputImage, TOutputImage>::EstimateOutputRpcM
   tempPtr->SetRegions(region);
 
   // Encapsulate the output metadata in the temp image
+  tempPtr->m_Imd.Add(MDGeom::ProjectionWKT, this->GetOutputProjectionRef());
   tempPtr->SetImageMetadata(*(this->GetOutputImageMetadata()));
-  tempPtr->m_Imd.Add(MDGeom::ProjectionProj, this->GetOutputProjectionRef());
 
   // Estimate the rpc model from the temp image
   m_OutputRpcEstimator->SetInput(tempPtr);