vtkPlanesIntersection.h

Go to the documentation of this file.

/*=========================================================================
 
  Program:   Visualization Toolkit
  Module:    vtkPlanesIntersection.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 (c) Sandia Corporation
 See Copyright.txt or http://www.paraview.org/HTML/Copyright.html for details.
----------------------------------------------------------------------------*/
 
#ifndef vtkPlanesIntersection_h
#define vtkPlanesIntersection_h
 
#include "vtkCommonDataModelModule.h" // For export macro
#include "vtkPlanes.h"
 
class vtkPoints;
class vtkPointsProjectedHull;
class vtkCell;
 
class VTKCOMMONDATAMODEL_EXPORT vtkPlanesIntersection : public vtkPlanes
{
  vtkTypeMacro(vtkPlanesIntersection, vtkPlanes);
 
public:
  void PrintSelf(ostream& os, vtkIndent indent) override;
 
  static vtkPlanesIntersection* New();
 
  void SetRegionVertices(vtkPoints* pts);
  void SetRegionVertices(double* v, int nvertices);
  int GetNumberOfRegionVertices();
  // Retained for backward compatibility
  int GetNumRegionVertices() { return this->GetNumberOfRegionVertices(); }
  int GetRegionVertices(double* v, int nvertices);
 
  int IntersectsRegion(vtkPoints* R);
 
  static int PolygonIntersectsBBox(double bounds[6], vtkPoints* pts);
 
  static vtkPlanesIntersection* Convert3DCell(vtkCell* cell);
 
protected:
  static void ComputeNormal(double* p1, double* p2, double* p3, double normal[3]);
  static double EvaluatePlaneEquation(double* x, double* p);
  static void PlaneEquation(double* n, double* x, double* p);
  static int GoodNormal(double* n);
  static int Invert3x3(double M[3][3]);
 
  vtkPlanesIntersection();
  ~vtkPlanesIntersection() override;
 
private:
  int IntersectsBoundingBox(vtkPoints* R);
  int EnclosesBoundingBox(vtkPoints* R);
  int EvaluateFacePlane(int plane, vtkPoints* R);
  int IntersectsProjection(vtkPoints* R, int direction);
 
  void SetPlaneEquations();
  void ComputeRegionVertices();
 
  void planesMatrix(int p1, int p2, int p3, double M[3][3]) const;
  int duplicate(double testv[3]) const;
  void planesRHS(int p1, int p2, int p3, double r[3]) const;
  int outsideRegion(double v[3]);
 
  // plane equations
  double* Planes;
 
  // vertices of convex regions enclosed by the planes, also
  //    the ccw hull of that region projected in 3 orthog. directions
  vtkPointsProjectedHull* RegionPts;
 
  vtkPlanesIntersection(const vtkPlanesIntersection&) = delete;
  void operator=(const vtkPlanesIntersection&) = delete;
};
#endif