From 3320112ebb474394a5ec4c9bcea5c57438298582 Mon Sep 17 00:00:00 2001
From: Christophe Palmann <christophe.palmann@c-s.fr>
Date: Wed, 1 Jul 2015 10:22:09 +0200
Subject: [PATCH] WRG: despeckling SAR filters (SAR)

---
 .../Filtering/ImageNoise/include/otbGammaMAPImageFilter.txx   | 2 +-
 Modules/Filtering/ImageNoise/include/otbKuanImageFilter.txx   | 4 +---
 2 files changed, 2 insertions(+), 4 deletions(-)

diff --git a/Modules/Filtering/ImageNoise/include/otbGammaMAPImageFilter.txx b/Modules/Filtering/ImageNoise/include/otbGammaMAPImageFilter.txx
index 125f533e8e..10883ddc81 100644
--- a/Modules/Filtering/ImageNoise/include/otbGammaMAPImageFilter.txx
+++ b/Modules/Filtering/ImageNoise/include/otbGammaMAPImageFilter.txx
@@ -119,7 +119,7 @@ void GammaMAPImageFilter<TInputImage, TOutputImage>::ThreadedGenerateData(
   InputRealType sum;
   InputRealType sum2;
 
-  double Ci, Ci2, Cu, Cu2, w, E_I, I, Var_I, dPixel, alpha, b, d, Cmax;
+  double Ci, Ci2, Cu, Cu2, E_I, I, Var_I, dPixel, alpha, b, d, Cmax;
 
   //Compute the ratio using the number of looks
   Cu2 = 1.0/m_NbLooks;
diff --git a/Modules/Filtering/ImageNoise/include/otbKuanImageFilter.txx b/Modules/Filtering/ImageNoise/include/otbKuanImageFilter.txx
index b1190fcbbe..351aed9dcd 100644
--- a/Modules/Filtering/ImageNoise/include/otbKuanImageFilter.txx
+++ b/Modules/Filtering/ImageNoise/include/otbKuanImageFilter.txx
@@ -119,11 +119,10 @@ void KuanImageFilter<TInputImage, TOutputImage>::ThreadedGenerateData(
   InputRealType sum;
   InputRealType sum2;
 
-  double Ci, Ci2, Cu, Cu2, w, E_I, I, Var_I, dPixel, alpha, b, d, Cmax;
+  double  Ci2, Cu2, w, E_I, I, Var_I, dPixel;
 
   //Compute the ratio using the number of looks
   Cu2 = 1.0/m_NbLooks;
-  Cu = vcl_sqrt(Cu2);
 
   // Process each of the boundary faces.  These are N-d regions which border
   // the edge of the buffer.
@@ -161,7 +160,6 @@ void KuanImageFilter<TInputImage, TOutputImage>::ThreadedGenerateData(
       I = static_cast<double>(bit.GetCenterPixel());
       
       Ci2 = Var_I / (E_I * E_I);
-      Ci  = vcl_sqrt(Ci2);
 
       const double epsilon = 0.0000000001;
       if (vcl_abs(E_I) < epsilon)
-- 
GitLab