IndexedPointInAreaLocator.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2006 Refractions Research Inc.
 * Copyright (C) 2018 Daniel Baston <dbaston@gmail.com>
 *
 * This is free software; you can redistribute and/or modify it under
 * the terms of the GNU Lesser General Public Licence as published
 * by the Free Software Foundation.
 * See the COPYING file for more information.
 *
 *
 **********************************************************************/
 
#ifndef GEOS_ALGORITHM_LOCATE_INDEXEDPOINTINAREALOCATOR_H
#define GEOS_ALGORITHM_LOCATE_INDEXEDPOINTINAREALOCATOR_H
 
#include <geos/geom/LineSegment.h>
#include <geos/algorithm/locate/PointOnGeometryLocator.h> // inherited
#include <geos/index/ItemVisitor.h> // inherited
#include <geos/index/intervalrtree/SortedPackedIntervalRTree.h> // inherited
 
#include <memory>
#include <vector> // composition
 
namespace geos {
namespace algorithm {
class RayCrossingCounter;
}
namespace geom {
class Geometry;
class Coordinate;
class CoordinateSequence;
}
}
 
namespace geos {
namespace algorithm { // geos::algorithm
namespace locate { // geos::algorithm::locate
 
class IndexedPointInAreaLocator : public PointOnGeometryLocator {
private:
    class IntervalIndexedGeometry {
    private:
        index::intervalrtree::SortedPackedIntervalRTree index;
        bool isEmpty;
 
        void init(const geom::Geometry& g);
        void addLine(const geom::CoordinateSequence* pts);
 
        // To keep track of LineSegments
        std::vector< geom::LineSegment > segments;
 
    public:
        IntervalIndexedGeometry(const geom::Geometry& g);
 
        void query(double min, double max, index::ItemVisitor* visitor);
    };
 
 
    class SegmentVisitor : public index::ItemVisitor {
    private:
        algorithm::RayCrossingCounter* counter;
 
    public:
        SegmentVisitor(algorithm::RayCrossingCounter* p_counter)
            :   counter(p_counter)
        { }
 
        ~SegmentVisitor() override
        { }
 
        void visitItem(void* item) override;
    };
 
 
    const geom::Geometry& areaGeom;
    std::unique_ptr<IntervalIndexedGeometry> index;
 
    void buildIndex(const geom::Geometry& g);
 
    // Declare type as noncopyable
    IndexedPointInAreaLocator(const IndexedPointInAreaLocator& other) = delete;
    IndexedPointInAreaLocator& operator=(const IndexedPointInAreaLocator& rhs) = delete;
 
public:
    IndexedPointInAreaLocator(const geom::Geometry& g);
 
    geom::Location locate(const geom::Coordinate* /*const*/ p) override;
 
};
 
} // geos::algorithm::locate
} // geos::algorithm
} // geos
 
#endif // GEOS_ALGORITHM_LOCATE_INDEXEDPOINTINAREALOCATOR_H