vtkHighestDensityRegionsStatistics.h

Go to the documentation of this file.

/*=========================================================================

Program:   Visualization Toolkit
Module:    vtkHighestDensityRegionsStatistics.h

Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.

This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE.  See the above copyright notice for more information.

=========================================================================*/

#ifndef vtkHighestDensityRegionsStatistics_h
#define vtkHighestDensityRegionsStatistics_h

#include "vtkFiltersStatisticsModule.h" // For export macro
#include "vtkStatisticsAlgorithm.h"

class vtkMultiBlockDataSet;
class vtkVariant;

class VTKFILTERSSTATISTICS_EXPORT vtkHighestDensityRegionsStatistics :
  public vtkStatisticsAlgorithm
{
public:
  vtkTypeMacro(vtkHighestDensityRegionsStatistics, vtkStatisticsAlgorithm);
  virtual void PrintSelf( ostream& os, vtkIndent indent );
  static vtkHighestDensityRegionsStatistics* New();


  virtual void Aggregate(vtkDataObjectCollection*,
                         vtkMultiBlockDataSet*) { return; }

  void SetSigma(double sigma);


  vtkGetVectorMacro(SmoothHC1, double, 2);
  vtkSetVectorMacro(SmoothHC1, double, 2);
  vtkGetVectorMacro(SmoothHC2, double, 2);
  vtkSetVectorMacro(SmoothHC2, double, 2);

  double ComputeHDR(vtkDataArray *inObservations, vtkDataArray *outDensity);


  double ComputeHDR(vtkDataArray *inObs, vtkDataArray* inPOI,
                    vtkDataArray *outDensity);
protected:
  vtkHighestDensityRegionsStatistics();
  ~vtkHighestDensityRegionsStatistics();


  virtual void Learn(vtkTable*,
                     vtkTable*,
                     vtkMultiBlockDataSet*);

  virtual void Derive(vtkMultiBlockDataSet*);


  virtual void Assess(vtkTable*,
                      vtkMultiBlockDataSet*,
                      vtkTable*) { return; }


  virtual void Test(vtkTable*,
                    vtkMultiBlockDataSet*,
                    vtkTable*) { return; }

//BTX

  virtual void SelectAssessFunctor(vtkTable*,
                                   vtkDataObject*,
                                   vtkStringArray*,
                                   AssessFunctor*&) { return; }
//ETX


  double SmoothHC1[2];
  double SmoothHC2[2];

  vtkIdType NumberOfRequestedColumnsPair;

private :
  double ComputeSmoothGaussianKernel(int dimension, double khx, double khy);

  double ComputeStandardGaussianKernel(int dimension, double kx, double ky);

private:
  vtkHighestDensityRegionsStatistics(const vtkHighestDensityRegionsStatistics&); // Not implemented
  void operator = (const vtkHighestDensityRegionsStatistics&);  // Not implemented
};

#endif