cmtkTypedArrayNoiseEstimatorBrummer.h

Go to the documentation of this file.

/*
//
//  Copyright 2008-2010 SRI International
//
//  This file is part of the Computational Morphometry Toolkit.
//
//  http://www.nitrc.org/projects/cmtk/
//
//  The Computational Morphometry Toolkit is free software: you can
//  redistribute it and/or modify it under the terms of the GNU General Public
//  License as published by the Free Software Foundation, either version 3 of
//  the License, or (at your option) any later version.
//
//  The Computational Morphometry Toolkit is distributed in the hope that it
//  will be useful, but WITHOUT ANY WARRANTY; without even the implied
//  warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License along
//  with the Computational Morphometry Toolkit.  If not, see
//  <http://www.gnu.org/licenses/>.
//
//  $Revision: 2398 $
//
//  $LastChangedDate: 2010-10-05 14:54:37 -0700 (Tue, 05 Oct 2010) $
//
//  $LastChangedBy: torstenrohlfing $
//
*/

#include <cmtkconfig.h>

#include <Base/cmtkTypedArray.h>
#include <Base/cmtkTypedArrayNoiseEstimatorNaiveGaussian.h>
#include <Base/cmtkHistogram.h>

namespace 
cmtk
{

class
TypedArrayNoiseEstimatorBrummer :
  protected TypedArrayNoiseEstimatorNaiveGaussian
{
public:
  typedef TypedArrayNoiseEstimatorBrummer Self;

  typedef TypedArrayNoiseEstimatorNaiveGaussian Superclass;

  TypedArrayNoiseEstimatorBrummer( const TypedArray* data, const size_t histogramBins = 255 );
  
protected:
  TypedArrayNoiseEstimatorBrummer() {}

  static double SijbersBiasHat ( const Histogram<unsigned int>::SmartPtr histogram, const double sigmaHat, const int numBinsToUse );

  static double SijbersLogLikelihood ( const Histogram<unsigned int>::SmartPtr histogram, const double sigma, const int numBinsToUse );

  static double SijbersVarHat ( const Histogram<unsigned int>::SmartPtr histogram, const double sigmaHat, const int numBinsToUse );

  static double EstimateNumBinsToUse ( const TypedArray* data, const Histogram<unsigned int>::SmartPtr histogram, const double sigmaHat );
};

} // namespace cmtk