Geometry Boost List: Re: [geometry] Querying index using intersection with custom geometry

Subject: Re: [geometry] Querying index using intersection with custom geometry
From: Adam Wulkiewicz (adam.wulkiewicz_at_[hidden])
Date: 2014-03-28 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/boost-geometry-spatial-query-shapes).
>
> 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