evFeatureDetectorSwirl.h

Go to the documentation of this file.

/*
 *
 * EVITA: Efficient Visualization of Terascale Datasets
 * Copyright (C) 2000-2016  Team Evita
 * 
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Library General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 * 
 * This library 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
 * Library General Public License for more details.
 * 
 * You should have received a copy of the GNU Library General Public
 * License along with this library; if not, write to the
 * Free Software Foundation, Inc., 675 Mass Ave, Cambridge,
 * MA 02139, USA.
 *
 */


#ifndef __evFeatureDetectorSwirl_h
#define __evFeatureDetectorSwirl_h


#include "evFeatureDetector.h"


class evFeatureDetectorSwirl : public evFeatureDetector
{
 public:
  evFeatureDetectorSwirl(void);
  virtual ~evFeatureDetectorSwirl(void);
  
  virtual const evString GetClassName(void) const
    {
      return(evString("evFeatureDetectorSwirl"));
    }
  
  int DetectFeatures(const evDataset &dataset,
                     evScalarField &significance_map) const;

 protected:

 private:

  void Initialize(void);
  int CalculateSwirlParameter(const evVectorField &velocity,
                              const evVectorField &velocity_gradient_u,
                              const evVectorField &velocity_gradient_w,
                              evScalarField &significance_map) const;
  double ComputeCompressibleTauBar(const double velocity_u,
                                   const double velocity_v,
                                   const double velocity_w,
                                   const evVector
                                   &velocity_gradient_u,
                                   const evVector
                                   &velocity_gradient_v,
                                   const evVector
                                   &velocity_gradient_w) const;
  double Determinant(const evVector &velocity_gradient_u,
                     const evVector &velocity_gradient_v,
                     const evVector &velocity_gradient_w) const;

  int VectorCoreDenoising(const evVectorField &velocity,
                          evScalarField &significance_map) const;
  bool Sperner(const evVectorField &velocity, const evAddress &address) const;
  unsigned int GetQuadrant(const evVector &vector) const;

};


#endif