You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

149 lines
5.1 KiB
C

/*=========================================================================
Program: Visualization Toolkit
Module: vtkComputeHistogram2DOutliers.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.
=========================================================================*/
/*-------------------------------------------------------------------------
Copyright 2009 Sandia Corporation.
Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
the U.S. Government retains certain rights in this software.
-------------------------------------------------------------------------*/
/**
* @class vtkComputeHistogram2DOutliers
* @brief compute the outliers in a set
* of 2D histograms and extract the corresponding row data.
*
*
* This class takes a table and one or more vtkImageData histograms as input
* and computes the outliers in that data. In general it does so by
* identifying histogram bins that are removed by a median (salt and pepper)
* filter and below a threshold. This threshold is automatically identified
* to retrieve a number of outliers close to a user-determined value. This
* value is set by calling SetPreferredNumberOfOutliers(int).
*
* The image data input can come either as a multiple vtkImageData via the
* repeatable INPUT_HISTOGRAM_IMAGE_DATA port, or as a single
* vtkMultiBlockDataSet containing vtkImageData objects as blocks. One
* or the other must be set, not both (or neither).
*
* The output can be retrieved as a set of row ids in a vtkSelection or
* as a vtkTable containing the actual outlier row data.
*
* @sa
* vtkExtractHistogram2D vtkPComputeHistogram2DOutliers
*
* @par Thanks:
* Developed by David Feng at Sandia National Laboratories
*------------------------------------------------------------------------------
*/
#ifndef vtkComputeHistogram2DOutliers_h
#define vtkComputeHistogram2DOutliers_h
//------------------------------------------------------------------------------
#include "vtkFiltersImagingModule.h" // For export macro
#include "vtkSelectionAlgorithm.h"
//------------------------------------------------------------------------------
class vtkCollection;
class vtkDoubleArray;
class vtkIdTypeArray;
class vtkImageData;
class vtkTable;
//------------------------------------------------------------------------------
class VTKFILTERSIMAGING_EXPORT vtkComputeHistogram2DOutliers : public vtkSelectionAlgorithm
{
public:
static vtkComputeHistogram2DOutliers* New();
vtkTypeMacro(vtkComputeHistogram2DOutliers, vtkSelectionAlgorithm);
void PrintSelf(ostream& os, vtkIndent indent) override;
vtkSetMacro(PreferredNumberOfOutliers, int);
vtkGetMacro(PreferredNumberOfOutliers, int);
//
vtkTable* GetOutputTable();
enum InputPorts
{
INPUT_TABLE_DATA = 0,
INPUT_HISTOGRAMS_IMAGE_DATA,
INPUT_HISTOGRAMS_MULTIBLOCK
};
enum OutputPorts
{
OUTPUT_SELECTED_ROWS = 0,
OUTPUT_SELECTED_TABLE_DATA
};
/**
* Set the source table data, from which data will be filtered.
*/
void SetInputTableConnection(vtkAlgorithmOutput* cxn)
{
this->SetInputConnection(INPUT_TABLE_DATA, cxn);
}
/**
* Set the input histogram data as a (repeatable) vtkImageData
*/
void SetInputHistogramImageDataConnection(vtkAlgorithmOutput* cxn)
{
this->SetInputConnection(INPUT_HISTOGRAMS_IMAGE_DATA, cxn);
}
/**
* Set the input histogram data as a vtkMultiBlockData set
* containing multiple vtkImageData objects.
*/
void SetInputHistogramMultiBlockConnection(vtkAlgorithmOutput* cxn)
{
this->SetInputConnection(INPUT_HISTOGRAMS_MULTIBLOCK, cxn);
}
protected:
vtkComputeHistogram2DOutliers();
~vtkComputeHistogram2DOutliers() override;
int PreferredNumberOfOutliers;
vtkTimeStamp BuildTime;
int RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*) override;
int FillInputPortInformation(int port, vtkInformation* info) override;
int FillOutputPortInformation(int port, vtkInformation* info) override;
/**
* Compute the thresholds (essentially bin extents) that contain outliers for
* a collection of vtkImageData histograms.
*/
virtual int ComputeOutlierThresholds(vtkCollection* histograms, vtkCollection* thresholds);
/**
* Compute the thresholds (bin extents) that contain outliers for a single vtkImageData histogram
*/
virtual int ComputeOutlierThresholds(
vtkImageData* histogram, vtkDoubleArray* thresholds, double threshold);
/**
* Take a set of range thresholds (bin extents) and filter out rows from the input table data that
* fits inside those thresholds.
*/
virtual int FillOutlierIds(
vtkTable* data, vtkCollection* thresholds, vtkIdTypeArray* rowIds, vtkTable* outTable);
private:
vtkComputeHistogram2DOutliers(const vtkComputeHistogram2DOutliers&) = delete;
void operator=(const vtkComputeHistogram2DOutliers&) = delete;
};
#endif