Hi, I am trying to build a rtree index with a collection if waypoint (about 30.000) each represented by longitude and latitude. When I use Cartesian coordinates rtree.query works as expected. However the results are slightly wrong because of my use of longitude,latitude in the cartesian domain. I then tried to refine the datatypes to use spherical_equatorial and geographic CoordinateSystem. Now I encounter two problems: a) With geographic, it does not compile. b) With spherical_equatorial query(nearest) gives clearly wrong results. Do you have any clues on what I am doing wrong? Attached is two images illustrating the query results. Search target is dot with label on it. I am using boost 1.55.0 on Linux with g++ 4.4.3 (same result with 4.6.3). typedef bg::model::point<double, 2, bg::cs::cartesian> geo_point; typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > geo_point_with_wrong_results; typedef bg::model::point<double, 2, bg::cs::geographic<bg::degree> > geo_point_does_not_compile; typedef std::pair<geo_point, unsigned int> geo_vg_value; void test() { iterate_over_all_waypoints (wp) { geo_point pt(wp->on, wp->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)); } Kind regards, Ola Martin
