Find file
Fetching contributors…
Cannot retrieve contributors at this time
129 lines (97 sloc) 4.37 KB
Program: Visualization Toolkit
Module: vtkPlanesIntersection.h
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or for details.
This software is distributed WITHOUT ANY WARRANTY; without even
PURPOSE. See the above copyright notice for more information.
Copyright (c) Sandia Corporation
See Copyright.txt or for details.
// .NAME vtkPlanesIntersection - A vtkPlanesIntersection object is a
// vtkPlanes object that can compute whether the arbitrary convex region
// bounded by it's planes intersects an axis-aligned box.
// .SECTION Description
// A subclass of vtkPlanes, this class determines whether it
// intersects an axis aligned box. This is motivated by the
// need to intersect the axis aligned region of a spacial
// decomposition of volume data with various other regions.
// It uses the algorithm from Graphics Gems IV, page 81.
// .SECTION Caveat
// An instance of vtkPlanes can be redefined by changing the planes,
// but this subclass then will not know if the region vertices are
// up to date. (Region vertices can be specified in SetRegionVertices
// or computed by the subclass.) So Delete and recreate if you want
// to change the set of planes.
#ifndef __vtkPlanesIntersection_h
#define __vtkPlanesIntersection_h
#include "vtkPlanes.h"
class vtkPoints;
class vtkPointsProjectedHull;
class vtkCell;
class VTK_FILTERING_EXPORT vtkPlanesIntersection : public vtkPlanes
vtkTypeMacro(vtkPlanesIntersection, vtkPlanes);
void PrintSelf(ostream& os, vtkIndent indent);
static vtkPlanesIntersection *New();
// Description:
// It helps if you know the vertices of the convex region.
// If you don't, we will calculate them. Region vertices
// are 3-tuples.
void SetRegionVertices(vtkPoints *pts);
void SetRegionVertices(double *v, int nvertices);
int GetNumRegionVertices();
int GetRegionVertices(double *v, int nvertices);
// Description:
// Return 1 if the axis aligned box defined by R intersects
// the region defined by the planes, or 0 otherwise.
int IntersectsRegion(vtkPoints *R);
// Description:
// A convenience function provided by this class, returns
// 1 if the polygon defined in pts intersects the bounding
// box defined in bounds, 0 otherwise.
// The points must define a planar polygon.
static int PolygonIntersectsBBox(double bounds[6], vtkPoints *pts);
// Description:
// Another convenience function provided by this class, returns
// the vtkPlanesIntersection object representing a 3D
// cell. The point IDs for each face must be given in
// counter-clockwise order from the outside of the cell.
static vtkPlanesIntersection *Convert3DCell(vtkCell *cell);
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]);
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 *Plane;
// 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&); // Not implemented
void operator=(const vtkPlanesIntersection&); // Not implemented