<div dir="ltr">Due to finite precision of coordinates (numeric rounding), generally the polygons created by ST_Subdivide (or any other function which adds vertices) are never going to match exactly with the source polygons.<div><br></div><div>How you handle this depends on what you are trying to compute.  One thing you might be able to do is to filter out polygons which have a very small area of overlap.</div><div><br></div><div>Hopefully a future version of PostGIS will provide the ability to specify a snapping tolerance for spatial predicates.  In that case using relate with an "interiorIntersects" pattern [1] and a tolerance should work for this.</div><div><br></div><div>[1] <a href="https://postgis.net/docs/ST_Relate.html">https://postgis.net/docs/ST_Relate.html</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 26, 2019 at 8:25 AM Eduardo Ramos <<a href="mailto:duda.ramos@gmail.com">duda.ramos@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hello all.</div><div><br></div><div>I'm exploring the function ST_Subdivide(with type geometry) and noticed that the nodes that this function creates to divide the original polygon are missing the polygon edges by ~6cm.</div><div><br></div><div>Here is an example:<br><br><div><img src="cid:ii_k3g26mti0" alt="image.png" width="543" height="314"></div><div>This is originating some issues when using the ST_Intersects because the new polygons generated by ST_Subdivide are intersecting polygons that were not being intersected before.</div><div><br></div><div>I guess this must be a matter of numeric precision but I we can't control the precision or the tolerance in either of the two functions (ST_Subdivide and ST_Intersects).</div><div><br></div><div>Is there a way to configure the precision of ST_Subdivide or to play with the ST_Intersects tolerance?</div><div><br></div><div>Any ideas or suggestions?</div><div><br></div><div>Thanks,</div><div>Eduardo<br></div><div><br></div></div></div>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/postgis-users</a></blockquote></div>