Geometry Boost List: [ggl] Polygon DP simplification

Subject: [ggl] Polygon DP simplification
From: V D (zedxz2)
Date: 2011-10-04 20:25:07

Hi,

Anyone had any luck specifying a distance strategy to DP ? 

Thanks,
-v

On Sep 30, 2011, at 11:06 AM, V D wrote:

> *Typo:
> boost::geometry(p1,p2, my_dist_strat);
> 
> Should read
> boost::geometry::distance(p1,p2, my_dist_strat);
> 
> 
> Sorry about that
> -v
> 
> On Sep 30, 2011, at 10:56 AM, V D wrote:
> 
>> Hi,
>> 
>> I'm trying to supply my own strategy to work around the problem for now.
>> Using the same code as before, except now I'm trying:
>> 
>> 
>> 27  typedef boost::geometry::strategy::distance::pythagoras<point_xy, point_xy, double> dist_strat;
>> 28  typedef boost::geometry::strategy::simplify::douglas_peucker<point_xy, dist_strat > simp_strat;
>> 
>> 30  simp_strat my_simplify_strat_test;  
>> 31  boost::geometry::simplify(original_1, simplified_1, 1.0 / 2048.0, my_simplify_strat_test);
>> 
>> 
>> But I cannot get the strategy to compile, I get a pretty nasty compile error :
>> 
>> 
>> instantiated from ?void boost::geometry::simplify(const Geometry&, Geometry&, const Distance&, const Strategy&) [with Geometry = main()::multi_polygon, Distance = double, Strategy = main()::simp_strat]?
>> simp_test.cpp:31:   instantiated from here
>> /opt/local/include/boost/geometry/strategies/concepts/distance_concept.hpp:232: error: no type named ?type? in ?struct boost::geometry::strategy::distance::services::strategy_point_point<boost::geometry::strategy::distance::comparable::pythagoras<boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, double> >?
>> /opt/local/include/boost/geometry/strategies/concepts/distance_concept.hpp:233: error: no type named ?type? in ?struct boost::geometry::strategy::distance::services::strategy_point_point<boost::geometry::strategy::distance::comparable::pythagoras<boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, double> >?
>> /opt/local/include/boost/geometry/strategies/concepts/distance_concept.hpp:244: error: no matching function for call to ?boost::geometry::strategy::distance::comparable::pythagoras<boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, double>::apply(const boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>&, const boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>&, const boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>&)?
>> /opt/local/include/boost/geometry/strategies/cartesian/distance_pythagoras.hpp:92: note: candidates are: static typename boost::geometry::select_calculation_type<S1, S2, CalculationType>::type boost::geometry::strategy::distance::comparable::pythagoras<Point1, Point2, CalculationType>::apply(const Point1&, const Point2&) [with Point1 = boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, Point2 = boost::geometry::model::d2::point_xy<double, boost::geometry::cs::cartesian>, CalculationType = double]
>> 
>> 
>> The dist_strat type works fine, for example:
>> 
>>   point_xy p1(0.0, 0.0)
>>   point_xy p2(5.0, 0.0)
>> 
>>   dist_strat my_dist_strat;
>>   boost::geometry(p1,p2, my_dist_strat);
>> 
>> But I can't seem to supply it to my deuglas peucker strategy ?
>> 
>> -v
>> 
>> 
>> 
>> 
>> On Sep 30, 2011, at 6:04 AM, Mateusz Loskot wrote:
>> 
>>> On 30/09/11 03:47, V D wrote:
>>>> Hum, if the problem is with the pythagoras distance, why is the
>>>> following snippet returning the expected value, 5.0.
>>>> Given the current error in the source, shouldn't it return 25 ?
>>>> 
>>>> #include <boost/geometry.hpp>
>>>> #include <boost/geometry/strategies/cartesian/distance_pythagoras.hpp>
>>>> #include <boost/geometry/geometries/point_xy.hpp>
>>>> #include <iostream>
>>>> 
>>>> int main()
>>>> {
>>>> typedef boost::geometry::model::d2::point_xy<double> point_xy;
>>>> 
>>>> point_xy p1(0.0, 0.0);
>>>> point_xy p2(5.0, 0.0);
>>>> 
>>>> boost::geometry::strategy::distance::pythagoras<point_xy, point_xy,
>>>> double> my_strategy;
>>>> std::cout << boost::geometry::distance(p1,p2, my_strategy) << std::endl;
>>>> 
>>>> return 0;
>>>> }
>>>> 
>>>> 
>>>> Is it only in the use of DP simplification that the pythagoras strategy
>>>> is problematic ?
>>> 
>>> The difference is in the actual strategies that kick in.
>>> I explained it in my last comment:
>>> 
>>> https://svn.boost.org/trac/boost/ticket/5954#comment:3
>>> 
>>> So, yes, my patch does not fix the problem, but sheds some light where
>>> is the bug. Though, I'm having problems with finding out why the
>>> different strategy is generated for simplify, not having enough time to
>>> dig it now.
>>> 
>>> p.s. Please, take care to respond to the mailing list too. I added ggl address to my response here.
>>> 
>>> Best regards,
>>> -
>>> Mateusz Loskot, http://mateusz.loskot.net
>>> Charter Member of OSGeo, http://osgeo.org
>>> Member of ACCU, http://accu.org
>> 
>