From: Leonardo Flórez-Valencia <florez-l@javeriana.edu.co>
Date: Tue, 7 Mar 2017 16:50:07 +0000 (-0500)
Subject: ...
X-Git-Tag: v0.1^0
X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=cpPlugins.git;a=commitdiff_plain;h=61b3659afe961ed248f30e26f9ca8f28fcfafddc

...
---

diff --git a/lib/cpExtensions/Algorithms/DiscontinuityMapImageFilter.hxx b/lib/cpExtensions/Algorithms/DiscontinuityMapImageFilter.hxx
index 877edbc..f1786b4 100644
--- a/lib/cpExtensions/Algorithms/DiscontinuityMapImageFilter.hxx
+++ b/lib/cpExtensions/Algorithms/DiscontinuityMapImageFilter.hxx
@@ -57,7 +57,7 @@ ThreadedGenerateData( const TRegion& region, itk::ThreadIdType threadId )
   for( fit = faceList.begin( ); fit != faceList.end( ); ++fit )
   {
     bit = itk::ConstNeighborhoodIterator< TInputImage >( this->GetRadius( ), input, *fit );
-    neighborhoodSize = bit.Size( );
+    neighborhoodSize = TOutputPixel( bit.Size( ) );
     it = itk::ImageRegionIterator< TOutputImage >( output, *fit );
     bit.OverrideBoundaryCondition( &nbc );
     bit.GoToBegin( );
@@ -71,10 +71,9 @@ ThreadedGenerateData( const TRegion& region, itk::ThreadIdType threadId )
         sum1 += v;
         sum2 += v * v;
       }
-      TOutputPixel mean = ( sum1 / TOutputPixel( neighborhoodSize ) );
       TOutputPixel var =
-        ( sum2 - ( ( sum1 * sum1 ) / TOutputPixel( neighborhoodSize ) ) ) /
-        ( TOutputPixel( neighborhoodSize ) - TOutputPixel( 1 ) );
+        ( sum2 - ( ( sum1 * sum1 ) / neighborhoodSize ) ) /
+        ( neighborhoodSize - TOutputPixel( 1 ) );
       if( this->m_UseSquareRoot )
         it.Set( std::sqrt( double( var ) ) );
       else