GeoTessPolygon3D.h

Go to the documentation of this file.

//- ****************************************************************************
//- 
//- 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.
//- 
//- BSD Open Source License.
//- All rights reserved.
//- 
//- Redistribution and use in source and binary forms, with or without
//- modification, are permitted provided that the following conditions are met:
//- 
//-    * Redistributions of source code must retain the above copyright notice,
//-      this list of conditions and the following disclaimer.
//-    * Redistributions in binary form must reproduce the above copyright
//-      notice, this list of conditions and the following disclaimer in the
//-      documentation and/or other materials provided with the distribution.
//-    * Neither the name of Sandia National Laboratories nor the names of its
//-      contributors may be used to endorse or promote products derived from
//-      this software without specific prior written permission.
//- 
//- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
//- AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
//- IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
//- ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
//- LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
//- CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
//- SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
//- INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
//- CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
//- ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//- POSSIBILITY OF SUCH DAMAGE.
//-
//- ****************************************************************************
 
#ifndef POLYGON3D_H_
#define POLYGON3D_H_
 
// **** _SYSTEM INCLUDES_ ******************************************************
 
#include <cstdio>
#include <iostream>
#include <iomanip>
#include <vector>
 
// use standard library objects
using namespace std;
 
// **** _LOCAL INCLUDES_ *******************************************************
 
#include "CPPUtils.h"
#include "GeoTessUtils.h"
#include "GeoTessGreatCircle.h"
#include "IFStreamAscii.h"
#include "GeoTessPolygon.h"
#include "GeoTessHorizon.h"
 
// **** _BEGIN GEOTESS NAMESPACE_ **********************************************
 
namespace geotess {
 
// **** _FORWARD REFERENCES_ ***************************************************
 
class GeoTessPosition;
 
// **** _CLASS DEFINITION_ *****************************************************
 
class GEOTESS_EXP_IMP GeoTessPolygon3D : public GeoTessPolygon
{
  private:
 
    GeoTessHorizon* bottom;
 
    GeoTessHorizon* top;
 
  public:
 
    ~GeoTessPolygon3D();
 
    GeoTessPolygon3D() : GeoTessPolygon(), bottom(NULL), top(NULL)
    {};
 
    GeoTessPolygon3D(vector<double*> points, GeoTessHorizon* h_bottom, GeoTessHorizon* h_top)
    : GeoTessPolygon(points), bottom(h_bottom), top(h_top)
    {}
 
    GeoTessPolygon3D(double* center, double radius, int nEdges, GeoTessHorizon* h_bottom, GeoTessHorizon* h_top)
    : GeoTessPolygon(center, radius, nEdges), bottom(h_bottom), top(h_top)
    {}
 
    GeoTessPolygon3D(string filename);
 
    virtual string class_name() { return "Polygon3D"; };
 
    GeoTessHorizon* getTop()
    {
        return top;
    }
 
    GeoTessHorizon* getBottom()
    {
        return bottom;
    }
 
    bool contains(const double* x, const double& radius, const int& layer, GeoTessProfile** profiles)
        {
        return (bottom->getLayerIndex() < 0 || layer >= bottom->getLayerIndex())
                && (top->getLayerIndex() < 0 || layer <= top->getLayerIndex())
                && radius > bottom->getRadius(x, profiles) - 1e-4
                && radius < top->getRadius(x, profiles) + 1e-4
                && GeoTessPolygon::contains(x);
    }
 
    bool contains(const double* x, const int& layer)
    {
        return ( bottom->getLayerIndex() < 0 || layer >= bottom->getLayerIndex())
                && ( top->getLayerIndex() < 0 || layer <= top->getLayerIndex())
                && GeoTessPolygon::contains(x);
    }
 
    bool contains(GeoTessPosition& position);
 
    bool containsAll(vector<double*>& points, vector<double>& radii, vector<int>& layers, vector<GeoTessProfile**>& profiles)
    {
        for (int i=0; i<(int)points.size(); ++i)
            if (!contains(points[i], radii[i], layers[i], profiles[i]))
                return false;
        return true;
    }
 
    bool containsAny(vector<double*>& points, vector<double>& radii, vector<int>& layers, vector<GeoTessProfile**>& profiles)
    {
        for (int i=0; i<(int)points.size(); ++i)
            if (contains(points[i], radii[i], layers[i], profiles[i]))
                return true;
        return false;
    }
 
    virtual void write(const string& outputFileName);
 
 
    virtual void loadAscii(vector<string>& records);
 
 
}; // end class Polygon
 
} // end namespace geotess
 
#endif /* POLYGON3D_H_ */