
Geometry : 
Subject: Re: [geometry] Querying index using intersection with custom geometry
From: Adam Wulkiewicz (adam.wulkiewicz_at_[hidden])
Date: 20140328 09:00:56
Hi Michael,
Michael Nett wrote:
> Hello,
>
> We're considering using Boost.Geometry in one of our projects, but we
> require the ability to query an index using a super set of the
> geometric shapes provided by the library (frustum, conical frustums
> and complex geometry).
>
> We have tried implementing template specializations as discussed in
> this thread
> (http://stackoverflow.com/questions/19490435/boostgeometryspatialqueryshapes).
>
> namespace boost { namespace geometry {
> template <typename Box> inline
> bool intersects(Box const& b, MyFrustum const& f)
> {
> // your implementation
> }
> }}
>
> However, compiling this provides the following error.
>
> /usr/include/boost/geometry/core/point_type.hpp:45:5: error: no
> matching function for call to
> 'assertion_failed(mpl_::failed************
> (boost::geometry::traits::point_type<MyGeometry>::NOT_IMPLEMENTED_FOR_THIS_POIN?
> T_TYPE::************)(mpl_::assert_::types<MyGeometry, mpl_::na,
> mpl_::na, mpl_::na>))'
>
> I tried digging into the source files, but my understanding of the
> geometry library is too limited to figure out what's going on.
>
> Can anyone advise on how to approach the problem?
>
This assertion fails if you try to get a point_type for the type which
is not adapted to any of the Boost.Geometry concepts. So you must call
something like geometry::point_type<MyGeometry>::type or pass an object
of MyGeometry to one of the existing functions of Boost.Geometry which
tries to use is as a properly adapted Geometry.
Did you check if your implementation of intersects() work without the rtree?
Btw, what's the Value type stored in the rtree?
Below is the fully working (for me) example, maybe it will help you to
find the cause:
#include <boost/geometry/index/rtree.hpp>
struct MyFrustum
{
MyFrustum(int d) : dummy(d) {}
int dummy;
};
namespace boost { namespace geometry {
// This will be called for Nodes Bounds and Indexables!
template <typename Box> inline
bool intersects(Box const& b, MyFrustum const& f)
{
std::cout << "checking the intersection with " << f.dummy << std::endl;
return true;// always intersects
}
}}
int main()
{
namespace bg = boost::geometry;
namespace bgi = bg::index;
namespace bgm = bg::model;
typedef bgm::point<float, 2, bg::cs::cartesian> pt;
typedef bgm::box<pt> box;
// check your implementation
bool ok = bg::intersects(box(pt(0, 0), pt(1, 1)), MyFrustum(5));
std::cout << "" << std::endl;
// now use the rtree
bgi::rtree<box, bgi::rstar<8> > rt;
// insert some values
rt.insert(box(pt(0, 0), pt(1, 1)));
rt.insert(box(pt(2, 2), pt(3, 3)));
// performa a query
std::vector<box> vec;
rt.query(bgi::intersects(MyFrustum(10)), std::back_inserter(vec));
std::cout << vec.size() << std::endl;
}
The above should print:
checking the intersection with 5

checking the intersection with 10
checking the intersection with 10
2
Let me know if it helps.
Regards,
Adam
Geometry list run by mateusz at loskot.net