ScanConversionPolyhedron.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_ScanConversionPolyhedron_h__)
#define __geom_ScanConversionPolyhedron_h__

#include "IndexedEdgePolyhedron.h"
#include "ScanConversionPolygon.h"

#include "../kernel/SegmentMath.h"
#include "../grid/RegularGrid.h"

#include <iostream>
#include <iomanip>
#include <vector>

#include <cassert>
#include <cmath>

BEGIN_NAMESPACE_GEOM


template<typename T = double>
class ScanConversionPolyhedron {
  //
  // Public types.
  //

public:

  typedef T Number;
  typedef ads::FixedArray<3,Number> Point;
  typedef SegmentMath<3,Number> Segment;

  //
  // Private types.
  //

private:

  typedef ScanConversionPolygon<Number> Polygon;

  // 
  // Data
  //

private:

  std::vector<Segment> _edges;

public:

  //--------------------------------------------------------------------------
  // @{

  ScanConversionPolyhedron() :
    _edges()
  {}

  ScanConversionPolyhedron(const ScanConversionPolyhedron& other) : 
    _edges(other._edges) 
  {}

  ScanConversionPolyhedron& 
  operator=(const ScanConversionPolyhedron& other);
  
  ScanConversionPolyhedron& 
  operator=(const IndexedEdgePolyhedron<Number>& x);
  
  ~ScanConversionPolyhedron()
  {}

  // @}
  //--------------------------------------------------------------------------
  // @{

  void 
  computeBBox(BBox<3,Number>* bb) const;

  void
  convertLocationsToIndices(const RegularGrid<3,Number>& grid);


  template<typename IndexOutputIterator>
  void
  scanConvert(IndexOutputIterator coordinates, 
              const RegularGrid<3,Number>& grid) const;

  // @}
  //--------------------------------------------------------------------------
  // @{

  const std::vector<Segment>& 
  getEdges() const { 
    return _edges; 
  }

  // @}
  //--------------------------------------------------------------------------
  // @{

  void 
  clear() {
    _edges.clear();
  }

  void 
  insertEdge(const Point& p, const Point& q);

  // @}

private:

  // Intersect with the plane of constant z value to obtain a polygon.
  void 
  computeZIntersection(Polygon* polygon, Number z) const;

};


//-----------------------------------------------------------------------------
// File I/O Operators
//-----------------------------------------------------------------------------



template<typename T>
std::ostream& 
operator<<(std::ostream& out, const ScanConversionPolyhedron<T>& polyhedron);



template<typename T>
void 
mathematicaPrint(std::ostream& out, 
                 const ScanConversionPolyhedron<T>& polyhedron);



template<typename T>
std::istream& 
operator>>(std::istream& in, ScanConversionPolyhedron<T>& polyhedron);



//-----------------------------------------------------------------------------
// Equality Operators
//-----------------------------------------------------------------------------


template<typename T>
inline
bool 
operator==(const ScanConversionPolyhedron<T>& a, 
           const ScanConversionPolyhedron<T>& b) { 
  return a.getEdges() == b.getEdges(); 
}



template<typename T>
inline
bool 
operator!=(const ScanConversionPolyhedron<T>& a, 
           const ScanConversionPolyhedron<T>& b) { 
  return !(a == b);
}


END_NAMESPACE_GEOM

#define __geom_ScanConversionPolyhedron_ipp__
#include "ScanConversionPolyhedron.ipp"
#undef __geom_ScanConversionPolyhedron_ipp__

#endif

---