DiscreteHausdorffDistance.h

/**********************************************************************
 *
 * GEOS - Geometry Engine Open Source
 * http://geos.osgeo.org
 *
 * Copyright (C) 2009  Sandro Santilli <strk@keybit.net>
 *
 * 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: algorithm/distance/DiscreteHausdorffDistance.java 1.5 (JTS-1.10)
 *
 **********************************************************************/
 
#ifndef GEOS_ALGORITHM_DISTANCE_DISCRETEHAUSDORFFDISTANCE_H
#define GEOS_ALGORITHM_DISTANCE_DISCRETEHAUSDORFFDISTANCE_H
 
#include <geos/export.h>
#include <geos/algorithm/distance/PointPairDistance.h> // for composition
#include <geos/algorithm/distance/DistanceToPoint.h> // for composition
#include <geos/util/IllegalArgumentException.h> // for inlines
#include <geos/geom/Geometry.h> // for inlines
#include <geos/util/math.h> // for inlines
#include <geos/geom/CoordinateFilter.h> // for inheritance
#include <geos/geom/CoordinateSequenceFilter.h> // for inheritance
 
#include <cstddef>
#include <vector>
 
#ifdef _MSC_VER
#pragma warning(push)
#pragma warning(disable: 4251) // warning C4251: needs to have dll-interface to be used by clients of class
#endif
 
namespace geos {
        namespace algorithm {
                //class RayCrossingCounter;
        }
        namespace geom {
                class Geometry;
                class Coordinate; 
                //class CoordinateSequence; 
        }
        namespace index {
                namespace intervalrtree {
                        //class SortedPackedIntervalRTree;
                }
        }
}
 
namespace geos {
namespace algorithm { // geos::algorithm
namespace distance { // geos::algorithm::distance

class GEOS_DLL DiscreteHausdorffDistance
{
public:
 
        static double distance(const geom::Geometry& g0,
                               const geom::Geometry& g1);
 
        static double distance(const geom::Geometry& g0,
                               const geom::Geometry& g1, double densifyFrac);
 
        DiscreteHausdorffDistance(const geom::Geometry& g0,
                                  const geom::Geometry& g1)
                :
                g0(g0),
                g1(g1),
                ptDist(),
                densifyFrac(0.0)
        {}

        void setDensifyFraction(double dFrac)
        {
                if ( dFrac > 1.0 || dFrac <= 0.0 )
                {
                        throw util::IllegalArgumentException(
                                "Fraction is not in range (0.0 - 1.0]");
                }
 
                densifyFrac = dFrac;
        }
 
        double distance()
        {
                compute(g0, g1);
                return ptDist.getDistance();
        }
 
        double orientedDistance()
        {
                computeOrientedDistance(g0, g1, ptDist);
                return ptDist.getDistance();
        }
 
        const std::vector<geom::Coordinate> getCoordinates() const
        {
                return ptDist.getCoordinates();
        }
 
        class MaxPointDistanceFilter : public geom::CoordinateFilter
        {
        public:
                MaxPointDistanceFilter(const geom::Geometry& geom)
                        :
                        geom(geom)
                {}
 
                void filter_ro(const geom::Coordinate* pt)
                {
                        minPtDist.initialize();
                        DistanceToPoint::computeDistance(geom, *pt,
                                                                 minPtDist);
                        maxPtDist.setMaximum(minPtDist);
                }
 
                const PointPairDistance& getMaxPointDistance() const
                {
                        return maxPtDist;
                }
 
        private:
                PointPairDistance maxPtDist;
                PointPairDistance minPtDist;
                DistanceToPoint euclideanDist;
                const geom::Geometry& geom;
 
        // Declare type as noncopyable
        MaxPointDistanceFilter(const MaxPointDistanceFilter& other);
        MaxPointDistanceFilter& operator=(const MaxPointDistanceFilter& rhs);
        };
 
        class MaxDensifiedByFractionDistanceFilter
                        : public geom::CoordinateSequenceFilter
        {
        public:
 
                MaxDensifiedByFractionDistanceFilter(
                                const geom::Geometry& geom, double fraction)
                        :
                        geom(geom),
                numSubSegs( std::size_t(util::round(1.0/fraction)) )
                {
                }
 
                void filter_ro(const geom::CoordinateSequence& seq,
                               std::size_t index);
 
                bool isGeometryChanged() const { return false; }
 
                bool isDone() const { return false; }
 
                const PointPairDistance& getMaxPointDistance() const {
                        return maxPtDist;
                }
 
        private:
                PointPairDistance maxPtDist;
                PointPairDistance minPtDist;
                const geom::Geometry& geom;
                std::size_t numSubSegs; // = 0;
 
        // Declare type as noncopyable
        MaxDensifiedByFractionDistanceFilter(const MaxDensifiedByFractionDistanceFilter& other);
        MaxDensifiedByFractionDistanceFilter& operator=(const MaxDensifiedByFractionDistanceFilter& rhs);
        };
 
private:
 
        void compute(const geom::Geometry& g0,
                     const geom::Geometry& g1)
        {
                computeOrientedDistance(g0, g1, ptDist);
                computeOrientedDistance(g1, g0, ptDist);
        }
 
        void computeOrientedDistance(const geom::Geometry& discreteGeom,
                                     const geom::Geometry& geom,
                                     PointPairDistance& ptDist);
 
        const geom::Geometry& g0;
 
        const geom::Geometry& g1;
 
        PointPairDistance ptDist;

    double densifyFrac; // = 0.0;
 
    // Declare type as noncopyable
    DiscreteHausdorffDistance(const DiscreteHausdorffDistance& other);
    DiscreteHausdorffDistance& operator=(const DiscreteHausdorffDistance& rhs);
};
 
} // geos::algorithm::distance
} // geos::algorithm
} // geos
 
#ifdef _MSC_VER
#pragma warning(pop)
#endif
 
#endif // GEOS_ALGORITHM_DISTANCE_DISCRETEHAUSDORFFDISTANCE_H