Получение координат точек из многоугольника Boost Geometry
У меня есть простая DLL, делающая некоторые вычисления с полигонами Boost Geometry. (В основном пересечения и различия.) Поскольку DLL, скорее всего, будет вызываться из кода C# и из Delphi, и кто знает откуда, я должен преобразовать результат в массивы, с которыми может справиться все.
ОБНОВЛЕНИЕ: я упростил и несколько исправил свой код. Новый код выглядит совершенно по-другому, использует совершенно другой подход (for_each_point) и как-то до сих пор не компилируется.
Мой новый код:
#include <vector>
#include <boost/range.hpp>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/polygon.hpp>
using namespace boost::geometry;
typedef boost::geometry::model::point
<
double, 2, boost::geometry::cs::spherical_equatorial<boost::geometry::degree>
> spherical_point;
class PointAggregator {
private :
double *x, *y;
int count;
public :
PointAggregator(int size) {
x = (double*) malloc(sizeof(double) * size);
y = (double*) malloc(sizeof(double) * size);
count = 0;
}
~PointAggregator() {
free(x);
free(y);
}
inline void operator()(spherical_point& p) {
x[count] = get<0>(p);
y[count] = get<1>(p);
count++;
}
void GetResult(double *resultX, double *resultY) {
resultX = x;
resultY = y;
}
};
void VectorToArray(std::vector<model::polygon<spherical_point>> resultVector, double x[], double y[], int *count) {
int i = 0;
for (std::vector<model::polygon<spherical_point>>::iterator it = resultVector.begin(); it != resultVector.end(); ++it) {
if (boost::size(*it) >= 2) {
*count = boost::size(*it);
PointAggregator* pa = new PointAggregator(*count);
boost::geometry::for_each_point(*it, *pa);
pa->GetResult(x, y);
delete(pa);
break;
}
}
}
Текущие ошибки компиляции:
- ошибка C2039: "тип": не является членом "boost::mpl::eval_if_c" iterator.hpp 63
- ошибка C3203: "тип": шаблон специализированного класса нельзя использовать в качестве аргумента шаблона для параметра шаблона "Итератор", ожидается реальный тип diff_type.hpp 25
- ошибка C2955: 'boost::type': использование шаблона класса требует списка аргументов шаблона diff_type.hpp 25
- ошибка C2955: 'boost::iterator_difference': использование шаблона класса требует списка аргументов шаблона diff_type.hpp 26
Какие из них не выглядят так, как будто они имеют какое-либо отношение к этой части кода (мое имя файла - geometry.cpp), но все остальное, использующее Boost Geometry, закомментировано, и я все еще получаю эти ошибки, так что...
Вот мой плохой код, который я имел ранее (отредактированный sehe)
(Я новичок в C++ и Boost, так что, возможно, я пропустил некоторую базовую концепцию при объединении кода из Интернета.) Я предполагаю, что я просто не могу перебрать многоугольник, который легко, и я пропустил нетривиальную часть, или что многоугольник нельзя использовать в качестве кольца, или итерации просто не такие, как я думал, или я понятия не имею, что еще может быть не так. Что я сделал не так?
2 ответа
Я нашел несколько вещей, которые нужно исправить:
- Я вижу одну проблему в ваших шаблонах. Обязательно ставьте пробелы!
- диапазон увеличения работает на контейнерах или диапазонах, которые содержат начало, конец пары
- Итераторы представляют собой что-то вроде указателя на объект. Получение размера итератора не будет делать то, что вы хотите. Вам нужно либо использовать boost:: size целого контейнера, либо std::distance(begin_iterator,end_iterator).
Вот версия, которая компилируется:
#include <vector>
#include <boost/range.hpp>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/polygon.hpp>
using namespace boost::geometry;
typedef boost::geometry::model::point
<
double, 2, boost::geometry::cs::spherical_equatorial<boost::geometry::degree>
> spherical_point;
class PointAggregator {
private :
double *x, *y;
int count;
public :
PointAggregator(int size) {
x = (double*) malloc(sizeof(double) * size);
y = (double*) malloc(sizeof(double) * size);
count = 0;
}
~PointAggregator() {
free(x);
free(y);
}
inline void operator()(spherical_point& p) {
x[count] = get<0>(p);
y[count] = get<1>(p);
count++;
}
void GetResult(double *resultX, double *resultY) {
resultX = x;
resultY = y;
}
};
// added spaces to the close brackets >> becomes > >
void VectorToArray(std::vector<model::polygon<spherical_point> > resultVector, double x[], double y[], int *count) {
for (std::vector<model::polygon<spherical_point> >::iterator it = resultVector.begin(); it != resultVector.end(); ++it) {
if (boost::size(resultVector) >= 2) {
// getting the size of the whole container
*count = boost::size(resultVector);
PointAggregator* pa = new PointAggregator(*count);
boost::geometry::for_each_point(*it, *pa);
pa->GetResult(x, y);
delete(pa);
break;
}
}
}
Хорошо, я думаю, у меня есть то, что вы ищете здесь. Я до сих пор не совсем понимаю, почему вы ищете этот диапазон, который, как я предполагаю, равен точкам, большим или равным 2, но я выяснил, как заставить его скомпилироваться, по крайней мере, с использованием boost::size().
Прежде всего, поймите, что первый параметр функции
void VectorToArray(std::vector<model::polygon<spherical_point> > resultVector, double x[], double y[], int *count)
{
...
}
является std:: vector, содержащим экземпляры типа model:: polygon.
Это означает, что когда вы разыменовываете свой итератор... определяется как
std::vector<model::polygon<spherical_point> >::iterator it
значение - это модель:: полигон.
boost:: model:: polygon НЕ является сам по себе Boost.Range-совместимым. boost:: model:: polygon - это тип, содержащий 5 функций-членов....
inline ring_type const& outer() const { return m_outer; }
inline inner_container_type const& inners() const { return m_inners; }
inline ring_type& outer() { return m_outer; }
inline inner_container_type & inners() { return m_inners; }
inline void clear()
{
m_outer.clear();
m_inners.clear();
}
Это означает, что ваш * it (т.е. модель:: полигон) ограничен вызовом только этих функций.
Похоже, что вы хотите сделать, это захватить либо внешнее кольцо, либо одно из внутренних колец каждого многоугольника в векторе (не знаю, какой, внутренний или внешний), и посмотреть, больше ли диапазон чего-либо в этом кольце или равно 2.
Чтобы сделать это, мы должны сделать еще несколько mpl и typedef.
typedef boost::geometry::model::point<double, 2, boost::geometry::cs::spherical_equatorial<boost::geometry::degree> > spherical_point; // your definition of a spherical_point
typedef boost::geometry::model::polygon<spherical_point> polygon; //consolidation of template args for a polygon
typedef boost::geometry::ring_type<polygon>::type ring_type; // define a ring_type that can handle your spherical_point by way of the polygon typedef.
typedef boost::geometry::interior_type<polygon>::type int_type; //define a interior_type that can handle your spherical_point
Чтобы завершить это и просто приступить к работе, я решил предположить, что вам нужно "внешнее" кольцо для вашего ограничения диапазона.
Вот для меня компиляция кода на gcc 4.1.1 с boost 1.48. Я оставляю правильной ли логику до кого-то еще.
using namespace boost::geometry;
typedef boost::geometry::model::point<double, 2, boost::geometry::cs::spherical_equatorial<boost::geometry::degree> > spherical_point;
typedef boost::geometry::model::polygon<spherical_point> polygon;
typedef boost::geometry::ring_type<polygon>::type ring_type;
typedef boost::geometry::interior_type<polygon>::type int_type;
class PointAggregator
{
private :
double *x, *y;
int count;
public :
PointAggregator(int size)
{
x = (double*) malloc(sizeof(double) * size);
y = (double*) malloc(sizeof(double) * size);
count = 0;
}
~PointAggregator()
{
free(x);
free(y);
}
inline void operator()(spherical_point& p)
{
x[count] = get<0>(p);
y[count] = get<1>(p);
count++;
}
void GetResult(double *resultX, double *resultY)
{
resultX = x;
resultY = y;
}
};
void VectorToArray(std::vector<model::polygon<spherical_point> > resultVector, double x[], double y[], int *count)
{
for (std::vector<model::polygon<spherical_point> >::iterator it = resultVector.begin(); it != resultVector.end(); ++it)
{
model::polygon<spherical_point> tmpPoly;
tmpPoly = (*it);
boost::geometry::ring_type<polygon>::type somering = tmpPoly.outer(); //typed it all out again instead of using ring_type since the complier was complaining and i didn't wanna get into it.
int ringsize = boost::size(somering);
if(ringsize >= 2)
{
*count = ringsize;
PointAggregator* pa = new PointAggregator(*count);
boost::geometry::for_each_point(*it, *pa);
pa->GetResult(x, y);
delete(pa);
break;
}
}
}