<div dir="ltr"><div><div>Hey, <br>PostGIS is not magic and bound by numeric precision (double), as every other computer program.<br><br></div><div>For the moment we don't know if this is causing your issue. Trying a translation would help determine that.<br>
</div><div><br></div>When you do processing on geom , it is wise to start by translating them, do your stuff, then translate it back.<br><br></div><div>If you provide me a simple test case, I can try it if you want.<br></div>
<div>(simple test case : 2 geom as wkt, the intersection, not working)<br><br></div>Cheers,<br>Rémi-C<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-06-27 13:48 GMT+02:00 Komяpa <span dir="ltr"><<a href="mailto:me@komzpa.net" target="_blank">me@komzpa.net</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
2014-06-27 14:21 GMT+03:00 Rémi Cura <<a href="mailto:remi.cura@gmail.com">remi.cura@gmail.com</a>>:<br>
<div class="">> Seems like a precision issue.<br>
> Can you try to translate your  geometry and round theim (ST_SnpaToGrid) a<br>
> bit to reduce the number of digits?<br>
<br>
</div>I'd prefer not to, as<br>
a) this is intermediate step in more complex processing;<br>
b) I bumped into this issue multiple times, and now I have a testing<br>
case dataset.<br>
<br>
Can something be done inside PostGIS not to fail on such geometries?<br>
<div class="HOEnZb"><div class="h5"><br>
> Cheers,<br>
> Rémi-C<br>
><br>
><br>
> 2014-06-27 12:18 GMT+02:00 Komяpa <<a href="mailto:me@komzpa.net">me@komzpa.net</a>>:<br>
>><br>
>> Hi,<br>
>><br>
>> I managed to create a set of geometries that I'd like to intersect.<br>
>><br>
>> POSTGIS="2.1.3 r12547" GEOS="3.4.2-CAPI-1.8.2 r3921" PROJ="Rel. 4.8.0,<br>
>> 6 March 2012" GDAL="GDAL 1.11.0, released 2014/04/16" LIBXML="2.9.1"<br>
>> LIBJSON="UNKNOWN" (core procs from "2.1.2 r12389" need upgrade)<br>
>> TOPOLOGY (topology procs from "2.1.2 r12389" need upgrade) RASTER<br>
>> (raster procs from "2.1.2 r12389" need upgrade)<br>
>><br>
>> Test table:<br>
>> <a href="http://pastebin.com/qSYLGfvs" target="_blank">http://pastebin.com/qSYLGfvs</a><br>
>><br>
>> After that:<br>
>><br>
>> select ST_Intersection(p.geom, g.geom) from bad_intersections p,<br>
>> bad_intersections g;<br>
>><br>
>> ERROR:  XX000: Error performing intersection: TopologyException: found<br>
>> non-noded intersection between LINESTRING (1.45784e+07 -1.40333e+06,<br>
>> 1.45784e+07 -1.40333e+06) and LINESTRING (1.45784e+07 -1.40333e+06,<br>
>> 1.45784e+07 -1.40333e+06) at 14578380.693885569 -1403332.55464<br>
>> LOCATION:  pg_error, lwgeom_pg.c:162<br>
>> Time: 10,674 ms<br>
>><br>
>> All geometries are valid.<br>
>><br>
>> Can this be fixed?<br>
>><br>
>> --<br>
>> Darafei "Komяpa" Praliaskouski<br>
>> OSM BY Team - <a href="http://openstreetmap.by/" target="_blank">http://openstreetmap.by/</a><br>
>> <a href="mailto:xmpp%3Ame@komzpa.net">xmpp:me@komzpa.net</a> mailto:<a href="mailto:me@komzpa.net">me@komzpa.net</a><br>
>> _______________________________________________<br>
>> postgis-devel mailing list<br>
>> <a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a><br>
>> <a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel</a><br>
><br>
><br>
><br>
> _______________________________________________<br>
> postgis-devel mailing list<br>
> <a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel</a><br>
<br>
<br>
<br>
--<br>
Darafei "Komяpa" Praliaskouski<br>
OSM BY Team - <a href="http://openstreetmap.by/" target="_blank">http://openstreetmap.by/</a><br>
<a href="mailto:xmpp%3Ame@komzpa.net">xmpp:me@komzpa.net</a> mailto:<a href="mailto:me@komzpa.net">me@komzpa.net</a><br>
_______________________________________________<br>
postgis-devel mailing list<br>
<a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-devel</a></div></div></blockquote></div><br></div>