Re: [Boost-users] Enquiry on boost.geometry and intersection()

Subject: Re: [Boost-users] Enquiry on boost.geometry and intersection()
From: Zhang Qun (zhangq.rhy_at_[hidden])
Date: 2018-11-19 09:27:44

Hi Adam
I think I have figured it out. Because of the precision, the two polygons
are not the same.

That also explains that why sometimes the intersection() gives me two
intersection points looking "exactly" the same but if you look at the
double precision, they are actually two different points. In reality, these
two points can be safely considered the same point. I can solve this by
comparing the two points distance. However, sometimes also because of the
precision, the intersection() does not even give me any intersection points
especially when the line is right across the vertex of the polygon. How do
I solve this kind of problem?

Ethan

On Mon, Nov 19, 2018 at 12:01 PM Zhang Qun <zhangq.rhy_at_[hidden]> wrote:

> Thanks Adam. I did use bg::reverse to reverse when necessary.
>
> Sorry that I still encounter a problem with the intesection() function. It
> got me weird results when the line just right across the vertices of the
> polygon.
>
> Here is what I did, trying to make it clear.
>
> 1. In my program, the intersection between line and polygon is zero. So I
> print out the polygon and the linestring. (Actually from my visualization,
> it should has one intersection at (9.29497, 2.29204))
>
> terminal printout:
> Linestring: ((10, 3.43403), (8.4973, 1))
> Polygon: (((7, 1.52705), (2.06343, 3.17257), (5.16667, 5.5), (7.69098,
> 5.5), (9.29497, 2.29204), (7, 1.52705)))
>
> 2. I reconstruct polygon and linestring using the above printout value,
> and redo the intersection(), weird thing happened, it got me two same
> intersection points.
>
> terminal printout:
> Intersections: ((9.29497, 2.29204), (9.29497, 2.29204))
>
> 3. I compare the polygon/linestring in my program and the reconstructed
> polygon/linestring, it says these two are not even equal.
>
> code:
> std::cout << "Polygons are spatially " << (boost::geometry::equals(poly1,
> poly2) ? "equal" : "not equal") << std::endl;
> std::cout << "Lines are spatially " << (boost::geometry::equals(line1,
> line2) ? "equal" : "not equal") << std::endl;
>
> terminal printout:
> Polygons are spatially not equal
> Lines are spatially not equal
>
> It really got me here. How come the geometry I re-constructed from the
> printout is not equal to the original? The code I used to re-construct is
> as follows:
>
> bgPolygon poly2{{{7, 1.52705},
>                         {2.06343, 3.17257},
>                         {5.16667, 5.5},
>                         {7.69098, 5.5},
>                         {9.29497, 2.29204},
>                         {7, 1.52705}}};  // C++11 initialization list
> support
>
>   bg::correct(poly2); // make it clockwise and closed polygon
>   bgLinestring line2{{10, 3.43403}, {8.4973, 1.0}};
>
> Could you help me on this issue?
>
> Thank you.
>
> On Fri, Nov 16, 2018, 9:25 PM Adam Wulkiewicz <adam.wulkiewicz_at_[hidden]
> wrote:
>
>> Please do not top-post.
>>
>> Zhang Qun wrote:
>>
>> Hi, Adam thanks for the reply and detailed code and explanation. Yes, it
>> worked.
>> I implemented as per your instructions but encountered a small problem.
>> That is sometimes I realized the intersections points between the
>> linestrings and the polygon are not in a "consistent" order, e.g. as shown
>> in the below image, lingstring 1 has P1->P2; However, instead of P3->P4,
>> linestring 2 is P4->P3. I actually expect linestring 2 is also along P3->P4
>> direction. Is this how it is supposed to be? Any workaround to solve this
>> problem? Thanks.
>>
>> [image: Presentation1.png]
>>
>> I don't know what happens in this case exactly because I don't know the
>> data but I would rather assume that Boost.Geometry doesn't care about the
>> order of linestring points (both directions represent the same geometry).
>> The simplest solution would be to check the direction of the output and if
>> needed reverse the linestrings manually or with bg::reverse().
>>
>> Adam
>>
>