<div dir="ltr"><div><div><div><div><div><div>Hi,<br><br></div>I had a strange behavior regarding otsu filter. My code is as simple as this:<br><br><br>    typedef itk::Statistics::ScalarImageToHistogramGenerator<ImageType> ScalarImageToHistogramGeneratorType;<br>    typedef ScalarImageToHistogramGeneratorType::HistogramType HistogramType;<br>    typedef itk::OtsuMultipleThresholdsCalculator<HistogramType> OtsuCalculatorType;<br><br>    ScalarImageToHistogramGeneratorType::Pointer scalarImageToHistogramGenerator = ScalarImageToHistogramGeneratorType::New();<br>    scalarImageToHistogramGenerator->SetNumberOfBins(128);<br>    scalarImageToHistogramGenerator->SetInput(itkImage);<br>    scalarImageToHistogramGenerator->Compute();<br><br><br>    OtsuCalculatorType::Pointer otsuCalculator = OtsuCalculatorType::New();<br>    otsuCalculator->SetNumberOfThresholds(numberOfThresholds);<br>    otsuCalculator->SetInputHistogram(scalarImageToHistogramGenerator->GetOutput());<br>    otsuCalculator->Update();<br><br>    const OtsuCalculatorType::OutputType& thresholdVector = otsuCalculator->GetOutput();<br>    HistogramType::MeasurementType threshold = thresholdVector[usedThreshold];<br><br><br></div>I tested this with itk configured with ITK_USE_REVIEW_STATISTICS set to ON and once to OFF.<br></div>With some images, the results is not the same. <br><br></div>Is this a known issue?<br><br></div>thanks,<br></div>Zein<br></div>