Opened 11 years ago
Closed 11 years ago
#2423 closed defect (fixed)
[REGRESSION] ST_LineToCurve gets semicircle wrong
Reported by: | strk | Owned by: | pramsey |
---|---|---|---|
Priority: | medium | Milestone: | PostGIS 2.0.4 |
Component: | postgis | Version: | 2.0.x |
Keywords: | curves | Cc: |
Description
Passing a semicircle polygon to ST_LineToCurve in 2.0.x results in a full circle with about half the area outside of the original lienear shape and not covering half of the original shape.
Result from 1.5.x is perfect.
Attachments (1)
Change History (8)
by , 11 years ago
Attachment: | curve_from_semicircle_line.png added |
---|
comment:1 by , 11 years ago
comment:2 by , 11 years ago
Keywords: | curves added |
---|
comment:3 by , 11 years ago
Beside breaking the arc at the appropriate place (when angle between edges changes drastically) I think we should also avoid to define circles with a single curve, because such a curve would effectively be only defined by 2 points (as first and last would be the same) thus allowing for infinite arcs. When a circle is detected (which would be wrong, in this case) the code should produce an output composed by two semicircle curves (but this is probably the topic for another ticket).
comment:4 by , 11 years ago
Or, by convention, we could assume that the control point of a full circle must always be on the other side of the diameter starting from the first/last point of the curve. In this case the control point should then be 10,0 rather than 0,10.
The actual output:
CURVEPOLYGON(CIRCULARSTRING(-10 0,-3.73640463187386e-14 10,-10 0))
The input:
st_intersection(st_buffer(st_point(0,0),10),'POLYGON((-10 0,10 0,10 10,-10 10,-10 0))')
comment:5 by , 11 years ago
The current code takes the num_points/2 vertex as the control point, which is just a guess. A better guess would be the vertex which the least difference in distance from last and first points. The ideal control point is equidistant from endpoints.
comment:6 by , 11 years ago
I've a tentative fix for this that checks angle between consecutive arc edges. Looks promising. The only failing test is for ticket #183 but the obtained result with the patch is actually what the reporter of #183 was hoping for (a circularstring followed by a linestring, in a compound curve) so I'd say it's a more correct output.
comment:7 by , 11 years ago
Resolution: | → fixed |
---|---|
Status: | new → closed |
I think the problem is with 2.0.x only considering distance from circumcenter rather than also angle between segments composing it. If angles were considered, like in the 1.5.x times, it will be clear that the last sement (the horizontal one at the bottom) would form a much smaller angle than the precedent ones.