Странно себя ведет использование boost :: geometry :: inside () для многоугольника и кольца

Я разрабатываю приложение на Debian 9 с подталкивание-1.62.0 установлен sudo apt-get install команда. Я нашел странное поведение, когда звонить boost::geometry::with_in(polygon, ring), Полигон пересекается с кольцом, но не с in. Вот код:

#include <iostream>
#include <boost/geometry/algorithms/intersection.hpp>
#include <boost/version.hpp>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/ring.hpp>
#include <boost/geometry/geometries/polygon.hpp>

using namespace std;
using namespace boost;
using namespace boost::geometry;

int main()
{
using Point = boost::geometry::model::d2::point_xy<double>;
using Polygon = boost::geometry::model::polygon<Point>;
using Ring = boost::geometry::model::ring<Point>;

Ring outer;
append(outer, Point(-6.46720129179205, 15.61591992211971));
append(outer, Point(-3.617204647384145, 15.61591992211971));
append(outer, Point(-3.617204647384145, 7.615919922119708));
append(outer, Point(-6.46720129179205, 7.615919922119708));
append(outer, Point(-6.46720129179205, 15.61591992211971));

Polygon polygon;
polygon.outer() = outer;
correct(polygon);

Ring ring;
append(ring, Point(-3.542202969588097, 7.615919922119709));
append(ring, Point(-3.542202969588097, 6.81591992211971));
append(ring, Point(-3.542202969588097, 6.015919922119708));
append(ring, Point(-3.542202969588097, 5.215919922119708));
append(ring, Point(-3.542202969588097, 4.415919922119708));
append(ring, Point(-3.542202969588097, 3.615919922119708));
append(ring, Point(-3.542202969588097, 2.815919922119709));
append(ring, Point(-3.542202969588097, 2.015919922119709));
append(ring, Point(-3.542202969588097, 1.215919922119709));
append(ring, Point(-3.542202969588097, 0.4159199221197081));
append(ring, Point(-3.542202969588097, -0.3840800778802915));
append(ring, Point(-6.542202969588097, -0.3840800778802916));
append(ring, Point(-6.542202969588097, 0.415919922119708));
append(ring, Point(-6.542202969588097, 1.215919922119709));
append(ring, Point(-6.542202969588097, 2.015919922119709));
append(ring, Point(-6.542202969588097, 2.815919922119709));
append(ring, Point(-6.542202969588097, 3.615919922119708));
append(ring, Point(-6.542202969588096, 4.41591992211971));
append(ring, Point(-6.542202969588097, 5.215919922119709));
append(ring, Point(-6.542202969588097, 6.015919922119709));
append(ring, Point(-6.542202969588097, 6.815919922119708));
append(ring, Point(-6.542202969588097, 7.615919922119708));
append(ring, Point(-3.542202969588097, 7.615919922119709));
correct(ring);

cout << boolalpha
<< within(polygon, ring) << endl /* should return false */
<< covered_by(polygon, ring) << endl; /* should retrurn false */

return 0;
}

Сборка и запуск:

 g++ -g -o wi wi.cpp
./wi

Получите результат:

 true
true

Многоугольник не является внутри кольца, но boost::geometry::within(polygon, ring) вернуть true,

boost_within

Так это ошибка Boost-1.62? И как это исправить.

1

Решение

Разница 7,615919922119709 — 7,615919922119708 составляет 1 × 10 ^ -15.

Это означает, что геометрия перекрывается по вертикали. Согласно спецификации стандарта Simple Feature это означает, что они удовлетворяют within:

Создать сплошной многоугольник в геометрии наддува

1

Другие решения

Других решений пока нет …

По вопросам рекламы [email protected]