Postgis + Boost :: геометрия + C ++ - PullRequest
Новая
       19

Postgis + Boost :: геометрия + C ++

0 голосов
/ 04 мая 2018

У меня есть следующие wtk: 

POLYGON((0 0, 10 0, 10 11, 11 10, 0 10))

В boost :: geometry мы можем получить это представление многоугольника, используя 

boost::geometry::dsv(polygon," "," "," ");

Как я могу применить функцию postgis st_makevalid к многоугольнику boost :: geometry? Я полагаю, что-то вроде: 

#include "libpq/libpq-fs.h"
#include "../libpq-fe.h"
.
.
.
std::string geom = "POLYGON" + boost::geometry::dsv(multipoly," "," "," ");

res = PQexecParams(conn, "SELECT st_makevalid(geom) FROM ....

Я не хочу подключаться к базе данных, я просто хочу восстановить многоугольник или многоугольник с помощью функции postgis st_makevalid.

1 Ответ

0 голосов
/ 04 мая 2018

Я сомневаюсь, что вам нужен PostGIS для этой операции.

Я также сомневаюсь, что есть способ "сделать его действительным". Поскольку многоугольник имеет четкое самопересечение:

enter image description here

Вот как вы делаете проверку и исправление в самой Boost Geometry:

Live On Coliru 

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/io/io.hpp>
#include <boost/geometry/algorithms/equals.hpp>
#include <iostream>

namespace bg = boost::geometry;
namespace bgm = boost::geometry::model;

template<typename G>
bool check(G const& g) {
    std::string reason;
    bool valid = bg::is_valid(g, reason);

    if (valid) std::cout << "Valid (dsv): " << bg::dsv(g) << "\n";
    else       std::cout << "Invalid: " << reason << "\n";

    return valid;
}

int main() {
    using pt = bgm::d2::point_xy<double>;
    using poly = bgm::polygon<pt>;

    poly p;
    bg::read_wkt("POLYGON((0 0, 10 0, 10 11, 11 10, 0 10))", p);

    while (!check(p)) {
        auto same = p;
        bg::correct(p);

        if (bg::equals(p, same)) {
            std::cout << "Out of ideas\n";
            break;
        }
    }
}

И обратите внимание на вывод: 

Invalid: Geometry is defined as closed but is open
Invalid: Geometry has invalid self-intersections. A self-intersection point was found at (10, 10); method: i; operations: u/i; segment IDs {source, multi, ring, segment}: {0, -1, -1, 1}/{0, -1, -1, 3}
Out of ideas

Если ваш источник на самом деле содержит подобные самопересечения, трудно сказать, что вы хотите. Возможно, вы хотите посмотреть на

...