id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
9759,boost::geometry::index::rtree does not work properly spherical_equatorial,ola.lykkja@…,Barend Gehrels,"When spherical_equatorial coordinate system is used with index::rtree, the query(nearest) gives wrong results.

My example insert 30.000 waypoint and query(nearest) mostly returns point east of the query point. When using cartesian, it works ok, except for the expected domain mismatch. Tested with boost 1.55.0 on Linux with g++ 4.4.3 and g++ 4.6.3 (identical results). The attached map plots shows the results from the two tests (cartesian and spherical_equatorial). Full source code, inclduing all waypoint data available upon request.

{{{
typedef bg::model::point<double, 2, bg::cs::cartesian> geo_point_ok;
typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<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(lon, 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));
}
}}}",Bugs,closed,To Be Determined,geometry,Boost 1.55.0,Problem,fixed,,