<div dir="ltr"><div><div><div><br><br>On Mon, Feb 12, 2018 at 6:22 PM, Paolo Cavallini <<a href="mailto:cavallini@faunalia.it">cavallini@faunalia.it</a>> wrote:<br>><br>> Il 11/02/2018 18:34, Markus Metz ha scritto:<br>><br>> > This implies that the data have been reprojected at some stage somewhere<br>> > from one SRS to the other. The most likely explanation for the different<br>> > results, particularly the topological errors reported by v.in.ogr, is<br>> > that non-topological polygons were reprojected. Can you make these geoms<br>> > in EPSG:3857 and EPSG:3003 available for testing?<br>><br>> Hi Markus,<br>> thanks for your thoughts. Here the data, original and simplified, in two<br>> CRS (3003 is simplified correctly, 3857 not):<br>> <a href="http://www.faunalia.eu/~paolo/test_generalize.tar.xz">http://www.faunalia.eu/~paolo/test_generalize.tar.xz</a><br><br></div>The original data in EPSG:3003 and EPSG:3857 are not identical. The original data in EPSG:3003 are topologically correct, while the original data in EPSG:3857 are not. Further on, the original data in EPSG:3857 contain features that are not present in the original data in EPSG:3003 which is quite strange.<br><br></div><div>Reprojecting the original data in EPSG:3003 to EPSG:3857 within GRASS works fine, also with subsequent v.generalize.<br></div><div><br></div>That means that the original data in EPSG:3857 are some reprojected version of the original original data (which are these?). The reprojection step, apparently performed on polygons, not GRASS areas (how did you reproject?), introduced topological errors. Please use native GRASS v.in.ogr + v.proj to reproject polygons.<br><br></div>Markus M<br><div><div><div><br>> All the best.<br>> --<br>> Paolo Cavallini - <a href="http://www.faunalia.eu">www.faunalia.eu</a><br>> QGIS & PostGIS courses: <a href="http://www.faunalia.eu/training.html">http://www.faunalia.eu/training.html</a><br>> <a href="https://www.google.com/trends/explore?date=all&geo=IT&q=qgis,arcgis">https://www.google.com/trends/explore?date=all&geo=IT&q=qgis,arcgis</a><br><br></div></div></div></div>