boost_mirror/geometry
1
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2025-05-10 15:44:07 +00:00
Code Issues Packages Projects Releases Wiki Activity
geometry/test/algorithms/test_difference.hpp
Barend Gehrels 8d202fd319 [geometry] Verified compatibility with/without robust rescaling 
We now have the define BOOST_GEOMETRY_NO_ROBUSTNESS to avoid rescaling
(probably should be renamed)
2014-02-12 11:45:02 +01:00

326 lines
10 KiB
C++
Raw Blame History

// Boost.Geometry (aka GGL, Generic Geometry Library)
// Unit Test
// Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
// 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_DIFFERENCE_HPP
#define BOOST_GEOMETRY_TEST_DIFFERENCE_HPP
#include <fstream>
#include <iomanip>
#include <boost/foreach.hpp>
#include <geometry_test_common.hpp>
#include <boost/range/algorithm/copy.hpp>
#include <boost/geometry/algorithms/correct.hpp>
#include <boost/geometry/algorithms/difference.hpp>
#include <boost/geometry/algorithms/sym_difference.hpp>
#include <boost/geometry/algorithms/area.hpp>
#include <boost/geometry/algorithms/length.hpp>
#include <boost/geometry/algorithms/num_points.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/multi/geometries/multi_point.hpp>
#include <boost/geometry/multi/geometries/multi_linestring.hpp>
#include <boost/geometry/multi/geometries/multi_polygon.hpp>
#include <boost/geometry/strategies/strategies.hpp>
#include <boost/geometry/io/wkt/wkt.hpp>
#if defined(TEST_WITH_SVG)
# define BOOST_GEOMETRY_DEBUG_SEGMENT_IDENTIFIER
# define BOOST_GEOMETRY_DEBUG_IDENTIFIER
# include <boost/geometry/io/svg/svg_mapper.hpp>
# include <boost/geometry/algorithms/detail/overlay/debug_turn_info.hpp>
#endif
template <typename Output, typename G1, typename G2>
void difference_output(std::string const& caseid, G1 const& g1, G2 const& g2, Output const& output)
{
#if defined(TEST_WITH_SVG)
{
typedef typename bg::coordinate_type<G1>::type coordinate_type;
typedef typename bg::point_type<G1>::type point_type;
std::ostringstream filename;
filename << "difference_"
<< caseid << "_"
<< string_from_type<coordinate_type>::name()
#if defined(BOOST_GEOMETRY_NO_ROBUSTNESS)
<< "_no_rob"
#endif
<< ".svg";
std::ofstream svg(filename.str().c_str());
bg::svg_mapper<point_type> mapper(svg, 500, 500);
mapper.add(g1);
mapper.add(g2);
mapper.map(g1, "fill-opacity:0.3;fill:rgb(51,51,153);stroke:rgb(51,51,153);stroke-width:3");
mapper.map(g2, "fill-opacity:0.5;fill:rgb(153,204,0);stroke:rgb(153,204,0);stroke-width:3");
for (typename Output::const_iterator it = output.begin(); it != output.end(); ++it)
{
mapper.map(*it,
//sym ? "fill-opacity:0.2;stroke-opacity:0.4;fill:rgb(255,255,0);stroke:rgb(255,0,255);stroke-width:8" :
"fill-opacity:0.2;stroke-opacity:0.4;fill:rgb(255,0,0);stroke:rgb(255,0,255);stroke-width:8");
}
}
#endif
}
template <typename OutputType, typename G1, typename G2>
void test_difference(std::string const& caseid, G1 const& g1, G2 const& g2,
int expected_count, int expected_point_count,
double expected_area,
double percentage = 0.0001,
bool sym = false)
{
std::vector<OutputType> clip;
typedef typename bg::coordinate_type<G1>::type coordinate_type;
if (sym)
{
bg::sym_difference(g1, g2, clip);
}
else
{
bg::difference(g1, g2, clip);
}
typename bg::default_area_result<G1>::type area = 0;
int n = 0;
for (typename std::vector<OutputType>::iterator it = clip.begin();
it != clip.end();
++it)
{
if (expected_point_count >= 0)
{
n += bg::num_points(*it);
}
area += bg::area(*it);
}
difference_output(caseid, g1, g2, clip);
#ifndef BOOST_GEOMETRY_DEBUG_ASSEMBLE
{
// Test inserter functionality
// Test if inserter returns output-iterator (using Boost.Range copy)
typedef typename bg::point_type<G1>::type point_type;
typedef typename bg::rescale_policy_type<point_type>::type
rescale_policy_type;
rescale_policy_type rescale_policy
= bg::get_rescale_policy<rescale_policy_type>(g1, g2);
std::vector<OutputType> inserted, array_with_one_empty_geometry;
array_with_one_empty_geometry.push_back(OutputType());
if (sym)
{
boost::copy(array_with_one_empty_geometry,
bg::detail::sym_difference::sym_difference_insert<OutputType>
(g1, g2, rescale_policy, std::back_inserter(inserted)));
}
else
{
boost::copy(array_with_one_empty_geometry,
bg::detail::difference::difference_insert<OutputType>(
g1, g2, rescale_policy, std::back_inserter(inserted)));
}
BOOST_CHECK_EQUAL(boost::size(clip), boost::size(inserted) - 1);
}
#endif
#if ! defined(BOOST_GEOMETRY_NO_BOOST_TEST)
if (expected_point_count >= 0)
{
BOOST_CHECK_MESSAGE(bg::math::abs(n - expected_point_count) < 3,
"difference: " << caseid
<< " #points expected: " << expected_point_count
<< " detected: " << n
<< " type: " << (type_for_assert_message<G1, G2>())
);
}
if (expected_count >= 0)
{
BOOST_CHECK_MESSAGE(int(clip.size()) == expected_count,
"difference: " << caseid
<< " #outputs expected: " << expected_count
<< " detected: " << clip.size()
<< " type: " << (type_for_assert_message<G1, G2>())
);
}
BOOST_CHECK_CLOSE(area, expected_area, percentage);
#endif
}
#ifdef BOOST_GEOMETRY_CHECK_WITH_POSTGIS
static int counter = 0;
#endif
template <typename OutputType, typename G1, typename G2>
void test_one(std::string const& caseid,
std::string const& wkt1, std::string const& wkt2,
int expected_count1,
int expected_point_count1,
double expected_area1,
int expected_count2,
int expected_point_count2,
double expected_area2,
double percentage = 0.0001)
{
#ifdef BOOST_GEOMETRY_CHECK_WITH_SQLSERVER
std::cout
<< "-- " << caseid << std::endl
<< "with qu as (" << std::endl
<< "select geometry::STGeomFromText('" << wkt1 << "',0) as p," << std::endl
<< "geometry::STGeomFromText('" << wkt2 << "',0) as q)" << std::endl
<< "select " << std::endl
<< " p.STDifference(q).STNumGeometries() as cnt1,p.STDifference(q).STNumPoints() as pcnt1,p.STDifference(q).STArea() as area1," << std::endl
<< " q.STDifference(p).STNumGeometries() as cnt2,q.STDifference(p).STNumPoints() as pcnt2,q.STDifference(p).STArea() as area2," << std::endl
<< " p.STDifference(q) as d1,q.STDifference(p) as d2 from qu" << std::endl << std::endl;
#endif
G1 g1;
bg::read_wkt(wkt1, g1);
G2 g2;
bg::read_wkt(wkt2, g2);
bg::correct(g1);
bg::correct(g2);
test_difference<OutputType>(caseid + "_a", g1, g2,
expected_count1, expected_point_count1,
expected_area1, percentage);
#ifdef BOOST_GEOMETRY_DEBUG_ASSEMBLE
return;
#endif
test_difference<OutputType>(caseid + "_b", g2, g1,
expected_count2, expected_point_count2,
expected_area2, percentage);
test_difference<OutputType>(caseid + "_s", g1, g2,
expected_count1 + expected_count2,
expected_point_count1 >= 0 && expected_point_count2 >= 0
? (expected_point_count1 + expected_point_count2) : -1,
expected_area1 + expected_area2,
percentage, true);
#ifdef BOOST_GEOMETRY_CHECK_WITH_POSTGIS
std::cout
<< (counter > 0 ? "union " : "")
<< "select " << counter++
<< ", '" << caseid << "' as caseid"
<< ", ST_NumPoints(ST_Difference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
<< ", ST_NumGeometries(ST_Difference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
<< ", ST_Area(ST_Difference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
//<< ", " << expected_area1 << " as expected_area_a"
//<< ", " << expected_count1 << " as expected_count_a"
<< ", ST_NumPoints(ST_Difference(ST_GeomFromText('" << wkt2 << "'), "
<< " ST_GeomFromText('" << wkt1 << "'))) "
<< ", ST_NumGeometries(ST_Difference(ST_GeomFromText('" << wkt2 << "'), "
<< " ST_GeomFromText('" << wkt1 << "'))) "
<< ", ST_Area(ST_Difference(ST_GeomFromText('" << wkt2 << "'), "
<< " ST_GeomFromText('" << wkt1 << "'))) "
//<< ", " << expected_area2 << " as expected_area_b"
//<< ", " << expected_count2 << " as expected_count_b"
<< ", ST_NumPoints(ST_SymDifference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
<< ", ST_NumGeometries(ST_SymDifference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
<< ", ST_Area(ST_SymDifference(ST_GeomFromText('" << wkt1 << "'), "
<< " ST_GeomFromText('" << wkt2 << "'))) "
//<< ", " << expected_area1 + expected_area2 << " as expected_area_s"
//<< ", " << expected_count1 + expected_count2 << " as expected_count_s"
<< std::endl;
#endif
}
template <typename OutputType, typename G1, typename G2>
void test_one_lp(std::string const& caseid,
std::string const& wkt1, std::string const& wkt2,
std::size_t expected_count,
int expected_point_count,
double expected_length)
{
G1 g1;
bg::read_wkt(wkt1, g1);
G2 g2;
bg::read_wkt(wkt2, g2);
bg::correct(g1);
std::vector<OutputType> pieces;
bg::difference(g1, g2, pieces);
typename bg::default_length_result<G1>::type length = 0;
std::size_t n = 0;
std::size_t piece_count = 0;
for (typename std::vector<OutputType>::iterator it = pieces.begin();
it != pieces.end();
++it)
{
if (expected_point_count >= 0)
{
n += bg::num_points(*it);
}
piece_count++;
length += bg::length(*it);
}
BOOST_CHECK_MESSAGE(piece_count == expected_count,
"difference: " << caseid
<< " #outputs expected: " << expected_count
<< " detected: " << pieces.size()
);
if (expected_point_count >= 0)
{
BOOST_CHECK_EQUAL(n, std::size_t(expected_point_count));
}
BOOST_CHECK_CLOSE(length, expected_length, 0.001);
std::string lp = "lp_";
difference_output(lp + caseid, g1, g2, pieces);
}
#endif
Reference in New Issue View Git Blame Copy Permalink