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);
  void PrintSelf( ostream& os, vtkIndent indent ) override;
  static vtkHighestDensityRegionsStatistics* New();

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

  void SetSigma(double sigma);

  void SetSigmaMatrix(double s11, double s12, double s21, double s22);

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

  double ComputeHDR(vtkDataArray *inObs, vtkDataArray* inPOI,
                    vtkDataArray *outDensity);

protected:
  vtkHighestDensityRegionsStatistics();
  ~vtkHighestDensityRegionsStatistics() override;

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

  void Derive(vtkMultiBlockDataSet*) override;

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

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

  void SelectAssessFunctor(vtkTable*,
                           vtkDataObject*,
                           vtkStringArray*,
                           AssessFunctor*&) override { return; }


  double SmoothHC1[2];
  double SmoothHC2[2];
  double InvSigmaC1[2];
  double InvSigmaC2[2];
  double Determinant;

  vtkIdType NumberOfRequestedColumnsPair;

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

private:
  vtkHighestDensityRegionsStatistics(const vtkHighestDensityRegionsStatistics&) = delete;
  void operator = (const vtkHighestDensityRegionsStatistics&) = delete;
};

#endif