BUG: Correction of the binarisation

the binarisation is now working with m_Min > 0

Author: jbvimort

include/itkCoocurrenceTextureFeaturesImageFilter.h

@@ -278,17 +278,18 @@ class ITK_TEMPLATE_EXPORT CoocurrenceTextureFeaturesImageFilter:public ImageToIm
     inline typename DigitalisedImageType::PixelType operator()(const MaskPixelType & maskPixel,
                                                                const PixelType & inputPixel) const
       {
+
         if(maskPixel != m_MaskValue)
           {
-          return this->m_Min - 10;
+          return - 10;
           }
         else if (inputPixel < this->m_Min || inputPixel >= m_Max)
           {
-          return m_Min - 1;
+          return - 1;
           }
         else
           {
-          return (inputPixel - m_Min)/((m_Max-m_Min)/ (float)m_NumberOfBinsPerAxis);
+          return ((inputPixel - m_Min)/((m_Max-m_Min)/ (float)m_NumberOfBinsPerAxis));
           }
       }
 

include/itkCoocurrenceTextureFeaturesImageFilter.hxx

@@ -178,7 +178,7 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage>
     while( !inputNIt.IsAtEnd() )
       {
       // If the voxel is outside of the mask, don't treat it
-      if( inputNIt.GetCenterPixel() < (this->m_Min - 5) ) //the pixel is outside of the mask
+      if( inputNIt.GetCenterPixel() < ( - 5) ) //the pixel is outside of the mask
         {
         progress.CompletedPixel();
         ++inputNIt;
@@ -198,7 +198,7 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage>
           {
           // Test if the curent voxel is in the mask and is the range of the image intensity sepcified
           curentInNeighborhoodPixelIntensity =  inputNIt.GetPixel(nb);
-          if( curentInNeighborhoodPixelIntensity < this->m_Min )
+          if( curentInNeighborhoodPixelIntensity < 0 )
             {
             continue;
             }
@@ -222,7 +222,7 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage>
 
           // Test if the pointed voxel is in the mask and is the range of the image intensity sepcified
           pixelIntensity = inputNIt.GetPixel(tempOffset);
-          if(pixelIntensity< this->m_Min )
+          if(pixelIntensity< 0 )
             {
             continue;
             }

include/itkRunLengthTextureFeaturesImageFilter.hxx

@@ -91,11 +91,11 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage>
     {
     if( maskPointer && maskPointer->GetPixel( inputIt.GetIndex() ) != this->m_InsidePixelValue )
       {
-      digitIt.Set(this->m_Min - 10);
+      digitIt.Set(-10);
       }
     else if(inputIt.Get() < this->m_Min || inputIt.Get() >= this->m_Max)
       {
-      digitIt.Set(this->m_Min - 1);
+      digitIt.Set(-1);
       }
     else
       {
@@ -196,7 +196,7 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage>
     while( !inputNIt.IsAtEnd() )
       {
       // If the voxel is outside of the mask, don't treat it
-      if( inputNIt.GetCenterPixel() < (this->m_Min - 5) ) //the pixel is outside of the mask
+      if( inputNIt.GetCenterPixel() < ( - 5) ) //the pixel is outside of the mask
         {
         progress.CompletedPixel();
         ++inputNIt;
@@ -224,7 +224,7 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage>
           curentInNeighborhoodPixelIntensity =  inputNIt.GetPixel(nb);
           tempOffset = inputNIt.GetOffset(nb);
           // Cecking if the value is out-of-bounds or is outside the mask.
-          if( curentInNeighborhoodPixelIntensity < this->m_Min || //the pixel is outside of the mask or outside of bounds
+          if( curentInNeighborhoodPixelIntensity < 0 || //the pixel is outside of the mask or outside of bounds
             alreadyVisitedImage->GetPixel( boolCurentInNeighborhoodIndex + tempOffset) )
             {
             continue;