mirror of https://github.com/CGAL/cgal
97 lines
3.8 KiB
C++
97 lines
3.8 KiB
C++
// Copyright (c) 2024 Max-Planck-Institute Saarbruecken (Germany), GeometryFactory (France)
|
|
// All rights reserved.
|
|
//
|
|
// This file is part of CGAL (www.cgal.org).
|
|
// You can redistribute it and/or modify it under the terms of the GNU
|
|
// General Public License as published by the Free Software Foundation,
|
|
// either version 3 of the License, or (at your option) any later version.
|
|
//
|
|
// Licensees holding a valid commercial license may use this file in
|
|
// accordance with the commercial license agreement provided with the software.
|
|
//
|
|
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
|
|
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
|
|
//
|
|
// $URL$
|
|
// $Id$
|
|
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
|
|
//
|
|
// Author(s) : André Nusser <anusser@mpi-inf.mpg.de>
|
|
// Marvin Künnemann <marvin@mpi-inf.mpg.de>
|
|
// Karl Bringmann <kbringma@mpi-inf.mpg.de>
|
|
// Andreas Fabri
|
|
// =============================================================================
|
|
|
|
#ifndef CGAL_FRECHET_DISTANCE_H
|
|
#define CGAL_FRECHET_DISTANCE_H
|
|
|
|
#include <CGAL/license/Frechet_distance.h>
|
|
#include <CGAL/basic.h>
|
|
#include <CGAL/Frechet_distance/internal/Frechet_distance.h>
|
|
|
|
#include <iterator>
|
|
|
|
namespace CGAL
|
|
{
|
|
|
|
/**
|
|
* \ingroup PkgFrechetDistanceFunctions
|
|
* decides if the Frechet distance between two polylines given as a range of
|
|
* points is less than a given distance.
|
|
* \param curve1 the first curve defined
|
|
* by the sequence of consecutive points along the polyline
|
|
* \param curve2 the
|
|
* second curve defined by the sequence of consecutive points along the polyline
|
|
* \param distance the decision distance
|
|
* \param traits the geometric traits object
|
|
* \tparam Traits a model of `FrechetDistanceTraits`
|
|
* \tparam PointRange a model of the concept `RandomAccessContainer`
|
|
* with `Traits::Point` as value type.
|
|
*/
|
|
template < class Traits, class PointRange>
|
|
bool continuous_Frechet_distance_less_than(const PointRange& curve1,
|
|
const PointRange& curve2,
|
|
const double distance,
|
|
const Traits& traits = Traits())
|
|
{
|
|
Protect_FPU_rounding<true> p;
|
|
auto icurve1 = Frechet_distance::internal::toCurve(curve1, traits);
|
|
auto icurve2 = Frechet_distance::internal::toCurve(curve2, traits);
|
|
auto idistance = Frechet_distance::internal::toDistance(distance);
|
|
|
|
return Frechet_distance::internal::lessThan(icurve1, icurve2, idistance, traits);
|
|
}
|
|
|
|
/**
|
|
* \ingroup PkgFrechetDistanceFunctions
|
|
* approximates the Fréchet distance between `curve1` and `curve2` up to an additive error
|
|
* of `precision` between two polylines given as a range of points.
|
|
*
|
|
* \param curve1 the first curve defined by the sequence of consecutive
|
|
* points along the polyline
|
|
* \param curve2 the second curve defined by the
|
|
* sequence of consecutive points along the polyline
|
|
* \param precision
|
|
* \param traits the geometric traits object
|
|
*
|
|
* \tparam Traits a model of `FrechetDistanceTraits`
|
|
* \tparam PointRange a model of the concept `RandomAccessContainer`
|
|
* with `Traits::Point` as value type.
|
|
*/
|
|
template <class Traits,class PointRange>
|
|
std::pair<double,double> continuous_Frechet_distance(const PointRange& curve1,
|
|
const PointRange& curve2,
|
|
const double precision,
|
|
const Traits& traits = Traits())
|
|
{
|
|
Protect_FPU_rounding<true> p;
|
|
auto icurve1 = Frechet_distance::internal::toCurve(curve1, traits);
|
|
auto icurve2 = Frechet_distance::internal::toCurve(curve2, traits);
|
|
|
|
return Frechet_distance::internal::calcDistance(icurve1, icurve2, traits, precision);
|
|
}
|
|
|
|
} // end of namespace CGAL
|
|
|
|
#endif // CGAL_FRECHET_DISTANCE_H
|