<div dir="ltr"><div>Hi Vedran</div><div><br></div><div>No, I do not want to change the SRS of the geometry. The geometries are defined in EPSG:3395. I want to do the reprojection on the fly with QGIS, and it is not doing what I would expect.</div><div>If I do the reprojection with "ogr2ogr -t_srs EPSG:3857" for both files and open them in QGIS, then yes, I see the difference. But I thought that QGIS was able to reproject on the fly.</div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Wed, 23 Jul 2025 at 08:40, Vedran Stojnović <<a href="mailto:phidrho@gmail.com">phidrho@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">Hi Javier,<div><br></div><div>if I understood properly what you are trying to achieve - you want overlap drawings (coordinates) from two coordinate systems as one image without transformation - then you need to override "Assigned coordinate system" (Right click on Layer -> Properties -> Source, and change Assigned Coordinate Reference System.</div><div>That way QGIS will think that these two datasets are in the same coordinate system and won't transform them "On the fly", otherwise it will reproject all data to a coordinate system currently set in project properties - see attached image. </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">uto, 22. srp 2025. u 17:03 Javier Jimenez Shaw via QGIS-User <<a href="mailto:qgis-user@lists.osgeo.org" target="_blank">qgis-user@lists.osgeo.org</a>> napisao je:<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>Hi</div><div><br></div><div>I am trying to compare web-mercator (EPSG:3857) with Mercator (EPSG:3395) projections.</div><div>For that purpose I think that a straight line in EPSG:3395 (loxodromic) shouldn't be exactly straight in EPSG:3857.</div><div><br></div><div>To show that, I created two linestrings in EPSG:3395, one with one segment, and the second with two segments co-linear (just adding a point in the middle). If I then select EPSG:3857 in QGIS, there should be a small difference... but I do not see it!</div><div><br></div><div>Only if I use the "Identify features", click in the area (at the proper zoom level), and hover over "identify all" or "2p". Then a second line appears! But it disappears when I select anything. (see attached screenshot)</div><div><br></div><div>The attached gpkg files are just the geometries as one segment (2p.gpkg) and 2 segments (3p.gpkg), both defined in EPSG:3395</div><div><br></div><div>The coordinates of the line are (0,0) and (8e6, 10e6) Yes, a long line to make it more visible.</div><div>the midpoint in EPSG:3395 is (4e6, 5e6). That is easy.</div><div><br></div><div>Reprojecting the two points into EPSG:3857 we have (0,0) and</div><div><br></div><div><span style="font-family:monospace">echo 8e6 10e6 | cs2cs EPSG:3395 EPSG:3857<br>8000000.00 10039255.88 0.00</span></div><div><br></div><div>The midpoint is</div><div></div><div><span style="font-family:monospace">echo "scale=2; 8000000.00 / 2; 10039255.88 / 2" | bc <br>4000000.00<br>5019627.94<br></span></div><div><br></div><div>The midpoint in EPSG:3395 and then projected to EPSG:3857 is this:</div><div><span style="font-family:monospace">echo 4e6 5e6 | cs2cs EPSG:3395 EPSG:3857<br>4000000.00 5028099.31 0.00</span></div><div><br></div><div>The difference is about 8 km. But I cannot see it in QGIS.</div><div><br></div><div>Jakub suggested that it could be the simplification. I disabled it in the settings. No difference.</div><div><br></div><div>Is QGIS doing something strange? Can I disable it?</div><div><br></div><div>If I select EPSG:4326 instead of EPSG:3857, then you can see clearly the two linestrings.</div><div><br></div><div>Thanks,</div><div>Javier.</div><div><br></div></div>
_______________________________________________<br>
QGIS-User mailing list<br>
<a href="mailto:QGIS-User@lists.osgeo.org" target="_blank">QGIS-User@lists.osgeo.org</a><br>
List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><br>
</blockquote></div><div><br clear="all"></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr">Srdačan pozdrav / Kind regards,<div>Vedran Stojnović.</div></div></div>
</blockquote></div>