X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FcpExtensions%2FAlgorithms%2FDiscontinuityMapImageFilter.hxx;fp=lib%2FcpExtensions%2FAlgorithms%2FDiscontinuityMapImageFilter.hxx;h=f1786b4d6f43e43fa5bed4d68d424a85e00168b9;hb=61b3659afe961ed248f30e26f9ca8f28fcfafddc;hp=877edbc4c42754e88939e0a08ea5e8a711653999;hpb=c467ae02a4568c1c245d06ecd5b9ebfd2d0a3b99;p=cpPlugins.git

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