MCIndexSegmentSetMutualIntersector.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2006 Refractions Research Inc.
 *
 * 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.
 *
 *
 **********************************************************************
 *
 * Last port: noding/MCIndexSegmentSetMutualIntersector.java r388 (JTS-1.12)
 *
 **********************************************************************/
 
#ifndef GEOS_NODING_MCINDEXSEGMENTSETMUTUALINTERSECTOR_H
#define GEOS_NODING_MCINDEXSEGMENTSETMUTUALINTERSECTOR_H
 
#include <geos/noding/SegmentSetMutualIntersector.h> // inherited
#include <geos/index/chain/MonotoneChainOverlapAction.h> // inherited
 
namespace geos {
        namespace index {
                class SpatialIndex;
 
                namespace chain {
                        class MonotoneChain;
                }
                namespace strtree {
                        //class STRtree;
                }
        }
        namespace noding {
                class SegmentString;
                class SegmentIntersector;
        }
}
 
//using namespace geos::index::strtree;
 
namespace geos {
namespace noding { // geos::noding

class MCIndexSegmentSetMutualIntersector : public SegmentSetMutualIntersector 
{
public:
 
        MCIndexSegmentSetMutualIntersector();
 
        ~MCIndexSegmentSetMutualIntersector();
 
        /* Returns a reference to a vector of MonotoneChain objects owned
         * by this class and destroyed on next call to ::process.
         * Copy them if you need them alive for longer.
         */
        std::vector<index::chain::MonotoneChain *>& getMonotoneChains() 
        { 
                return monoChains; 
        }
 
        index::SpatialIndex* getIndex() 
        { 
                return index; 
        }
 
        void setBaseSegments(SegmentString::ConstVect* segStrings);
  
        // NOTE: re-populates the MonotoneChain vector with newly created chains
        void process(SegmentString::ConstVect* segStrings);
 
    class SegmentOverlapAction : public index::chain::MonotoneChainOverlapAction
    {
    private:
        SegmentIntersector & si;
 
        // Declare type as noncopyable
        SegmentOverlapAction(const SegmentOverlapAction& other);
        SegmentOverlapAction& operator=(const SegmentOverlapAction& rhs);
 
    public:
        SegmentOverlapAction(SegmentIntersector & si) :
          index::chain::MonotoneChainOverlapAction(), si(si)
          {}
 
          void overlap(index::chain::MonotoneChain& mc1, std::size_t start1,
              index::chain::MonotoneChain& mc2, std::size_t start2);
    };
 
private:
 
        typedef std::vector<index::chain::MonotoneChain *> MonoChains;
        MonoChains monoChains;
 
        /*
         * The {@link SpatialIndex} used should be something that supports
         * envelope (range) queries efficiently (such as a {@link Quadtree}
         * or {@link STRtree}.
         */
        index::SpatialIndex * index;
        int indexCounter;
        int processCounter;
        // statistics
        int nOverlaps;
        
        /* memory management helper, holds MonotoneChain objects used
         * in the SpatialIndex. It's cleared when the SpatialIndex is
         */
        MonoChains chainStore;
      
        void addToIndex( SegmentString * segStr);
 
        void intersectChains();
 
        void addToMonoChains( SegmentString * segStr);
 
};
 
} // namespace geos::noding
} // namespace geos
 
#endif // GEOS_NODING_MCINDEXSEGMENTSETMUTUALINTERSECTOR_H