- typedef itk::MinimumMaximumImageCalculator< _TScalarImage > _TMinMax;
- typename _TMinMax::Pointer minMax = _TMinMax::New( );
- minMax->SetImage( this->GetInputVesselness( ) );
- minMax->Compute( );
- this->m_MinVesselness =
- ( this->m_VesselnessThreshold / double( 100 ) ) *
- double( minMax->GetMaximum( ) );
+ typedef itk::StatisticsImageFilter< _TScalarImage > _TStats;
+ typename _TStats::Pointer stats = _TStats::New( );
+ stats->SetInput( this->GetInputVesselness( ) );
+ stats->Update( );
+ double vAvg = double( stats->GetMean( ) );
+ double vStd = double( stats->GetSigma( ) );
+ double vMin = double( stats->GetMinimum( ) );
+ double vMax = double( stats->GetMaximum( ) );
+ std::exit( 1 );
+ /* TODO
+ this->m_MinVesselness =
+ ( this->m_VesselnessThreshold / double( 100 ) ) *
+ double( stats->GetMaximum( ) );
+ */