geometry/test/algorithms/similarity/test_hausdorff_distance.hpp

139 lines
4.6 KiB
C++

// Boost.Geometry
// Copyright (c) 2018 Yaghyavardhan Singh Khangarot, Hyderabad, India.
// Contributed and/or modified by Yaghyavardhan Singh Khangarot,
// as part of Google Summer of Code 2018 program.
// This file was modified by Oracle on 2018.
// Modifications copyright (c) 2018, Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_GEOMETRY_TEST_HAUSDORFF_DISTANCE_HPP
#define BOOST_GEOMETRY_TEST_HAUSDORFF_DISTANCE_HPP
#include <geometry_test_common.hpp>
#include <boost/geometry/algorithms/discrete_hausdorff_distance.hpp>
#include <boost/geometry/io/wkt/wkt.hpp>
#include <boost/geometry/strategies/strategies.hpp>
#include <boost/variant/variant.hpp>
template <typename Geometry1, typename Geometry2, typename Expected>
void test_hausdorff_distance(Geometry1 const& geometry1, Geometry2 const& geometry2,
Expected const& expected_hausdorff_distance)
{
typedef typename bg::distance_result
<
typename bg::point_type<Geometry1>::type,
typename bg::point_type<Geometry2>::type
>::type result_type;
result_type h_distance = bg::discrete_hausdorff_distance(geometry1, geometry2);
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::ostringstream out;
out << typeid(typename bg::coordinate_type<Geometry1>::type).name()
<< std::endl
<< typeid(typename bg::coordinate_type<Geometry2>::type).name()
<< std::endl
<< typeid(h_distance).name()
<< std::endl
<< "hausdorff_distance : " << h_distance
<< std::endl;
std::cout << out.str();
#endif
BOOST_CHECK_CLOSE(h_distance, result_type(expected_hausdorff_distance), 0.01);
}
template <typename Geometry1, typename Geometry2, typename Expected>
void test_geometry(std::string const& wkt1, std::string const& wkt2,
Expected const& expected_hausdorff_distance)
{
Geometry1 geometry1;
bg::read_wkt(wkt1, geometry1);
Geometry2 geometry2;
bg::read_wkt(wkt2, geometry2);
test_hausdorff_distance(geometry1, geometry2, expected_hausdorff_distance);
#if defined(BOOST_GEOMETRY_TEST_DEBUG)
test_hausdorff_distance(boost::variant<Geometry1>(geometry1),
boost::variant<Geometry2>(geometry2),
expected_hausdorff_distance);
#endif
}
template <typename Geometry1, typename Geometry2, typename Strategy, typename Expected>
void test_hausdorff_distance(Geometry1 const& geometry1, Geometry2 const& geometry2,
Strategy strategy, Expected const& expected_hausdorff_distance)
{
typedef typename bg::distance_result
<
typename bg::point_type<Geometry1>::type,
typename bg::point_type<Geometry2>::type,
Strategy
>::type result_type;
result_type h_distance = bg::discrete_hausdorff_distance(geometry1, geometry2, strategy);
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::ostringstream out;
out << typeid(typename bg::coordinate_type<Geometry1>::type).name()
<< std::endl
<< typeid(typename bg::coordinate_type<Geometry2>::type).name()
<< std::endl
<< typeid(h_distance).name()
<< std::endl
<< "hausdorff_distance : " << h_distance
<< std::endl;
std::cout << out.str();
#endif
BOOST_CHECK_CLOSE(h_distance, result_type(expected_hausdorff_distance), 0.01);
}
template <typename Geometry1, typename Geometry2, typename Strategy, typename Expected>
void test_geometry(std::string const& wkt1, std::string const& wkt2,
Strategy strategy, Expected const& expected_hausdorff_distance)
{
Geometry1 geometry1;
bg::read_wkt(wkt1, geometry1);
Geometry2 geometry2;
bg::read_wkt(wkt2, geometry2);
test_hausdorff_distance(geometry1, geometry2, strategy, expected_hausdorff_distance);
#if defined(BOOST_GEOMETRY_TEST_DEBUG)
test_hausdorff_distance(boost::variant<Geometry1>(geometry1),
boost::variant<Geometry2>(geometry2),
strategy, expected_hausdorff_distance);
#endif
}
template <typename Geometry1,typename Geometry2>
void test_empty_input(Geometry1 const& geometry1, Geometry2 const& geometry2)
{
try
{
bg::discrete_hausdorff_distance(geometry1, geometry2);
}
catch(bg::empty_input_exception const& )
{
return;
}
BOOST_CHECK_MESSAGE(false, "A empty_input_exception should have been thrown");
}
#endif