<div dir="ltr">Hello Marco,<div><br></div><div>Subdivision should be done in target projection as a final step of processing. If you expect to keep topology after reprojection, you should be either merging the polygon back before the reprojection, or investing in using the topology framework of postgis. That's the base of all overlay operations.<br><br></div><div>Insertion of an extra vertex breaks the logic of "every time you do a split on a node, you have less nodes in each half" and will lead to undesired behavior in sawtooth polygons. Thus behavior you're requesting has to be under a flag.</div><div><br></div><div>Please send in the pull request for keeping the nodes.</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Feb 24, 2021 at 1:08 PM Marco Boeringa <<a href="mailto:marco@boeringa.demon.nl">marco@boeringa.demon.nl</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">Hi all,<br>
<br>
I don't know how much of an effort it would be, but I think it would be <br>
nice if PostGIS's "Subdivide" was enhanced to insert a vertex where <br>
subdivision lines "intersect". Currently, when a large shape is <br>
subdivided, if a vertical sudivide line "hits" a horizontal one, the <br>
sub-part of the large geometry being subdivided that is above or below <br>
the horizontal subdivision line, won't receive a vertex where the <br>
vertical subdivision line "hits" it.<br>
<br>
This has the negative consequence that as soon as you start reprojecting <br>
the geometries resulting of this subdivide operation to another <br>
projection and depending on that projection's properties, e.g. with <br>
on-the-fly projection in QGIS, that the three geometries around this <br>
"virtual" intersection of subdivision lines, won't match, and show a <br>
clearly visible gap.<br>
<br>
This could be prevented if "Subdivide" always inserted an extra vertex <br>
at the intersection of the subdivision lines, even if the geometry part <br>
above or below the horizontal subdivision line, won't be cut vertically <br>
at that point.<br>
<br>
Marco<br>
<br>
_______________________________________________<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><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Darafei "Komяpa" Praliaskouski<br>OSM BY Team - <a href="http://openstreetmap.by/" target="_blank">http://openstreetmap.by/</a><br></div></div>