vtkHighestDensityRegionsStatistics.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause

#ifndef vtkHighestDensityRegionsStatistics_h
#define vtkHighestDensityRegionsStatistics_h

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

VTK_ABI_NAMESPACE_BEGIN
class vtkStatisticalModel;
class vtkVariant;

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

  int GetMaximumNumberOfColumnsPerRequest() const override { return 2; }

  bool Aggregate(vtkDataObjectCollection*, vtkStatisticalModel*) override { return false; }

  void AppendAlgorithmParameters(std::string& algorithmParameters) const override;

  std::size_t ConsumeNextAlgorithmParameter(
    vtkStringToken parameterName, const std::string& algorithmParameters) override;

  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*, vtkStatisticalModel*) override;

  void Derive(vtkStatisticalModel*) override;

  void Assess(vtkTable*, vtkStatisticalModel*, vtkTable*) override {}

  void Test(vtkTable*, vtkStatisticalModel*, vtkTable*) override {}

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


  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);

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

VTK_ABI_NAMESPACE_END
#endif