boost::geometry::index::rtree does not compile with CoordinateSystem geographic<

[ola.lykkja@…]

boost::geometry::index::rtree does not compile with CoordinateSystem geographic<

CoordinateSystem cartesian is ok

Example program:

typedef bg::model::point<double, 2, bg::cs::cartesian> geo_point_ok;
typedef bg::model::point<double, 2, bg::cs::geographic<bg::degree> > geo_point;
typedef std::pair<geo_point, unsigned int> geo_vg_value;

void test() {
      double lon=10.0, lat=60.0;
      geo_point pt(on, lat);
      rtree.insert(std::make_pair(pt, wp->id));

      std::vector<geo_vg_value> result;
      geo_point here(targetLat, targetLat);
      rtree.query(bgi::nearest(here, 20), std::back_inserter(result));
}

[awulkiew]

In 1.56 the geometry::comparable_distance() is used by the nearest query. After adding the support for geographic CS in this algorithm the rtree should work automatically.

[djh]

Is there a workaround in the meantime?

That is how do I tell Boost.Geometry to derive the distance between two points in the geographic coordinate system? From an implementation standpoint, I think this should be easily doable using e.g. the Haversine formula.

Is there a quick fix that lets the RTree work with (longitude, latitude)-pairs or similar in the geographical coordinate system?

[awulkiew]

If haversine distance is good for you, instead of geographic CS you could use spherical_equatorial.

[anonymous]

You need a point to MBR distance: ​https://link.springer.com/chapter/10.1007/978-3-642-40235-7_9