<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Jochen,</p>
<p>I don't have a ready-made theory to explain that. I've checked that recent and older PROJ versions (the library that powers QGIS for coordinate transformation) have identical output for EPSG:31287 both in proj.4 strings, WKT2 and WKT:ESRI (the later is the variant use in the .prj file of shapefiles). Looking at the EPSG database for <a class="moz-txt-link-freetext" href="https://epsg.org/crs_31287/MGI-Austria-Lambert.html">https://epsg.org/crs_31287/MGI-Austria-Lambert.html</a> also shows no changes in ages in its definition. The definition in the open ESRI CRS database hasn't changed either in years.
</p>
<p>The fact you are observing the same issue with ArcGIS Pro is really intriguing. Whereas ESRI software uses part of the open source GIS stack in their software, to the best of my knowledge, they don't use PROJ but their own Projection Engine. Sorry for asking a question that may sound like the one a typical operator of a help center would ask you :-), but are you really sure you are testing the exact same files between QGIS 3.28 and 3.36/ArcGIS Pro ...?</p>
<p>Does the issue if on-the-fly reprojection occurs (and if so to which CRS are you reprojecting to?), or also if you keep your project in EPSG:31287
</p>
<p>Otherwise please share all the files consisting of his shapefile.</p>
<p>For reference, the expected content for a .prj file using EPSG:31287 should be:</p>
<p>PROJCS["MGI_Austria_Lambert",GEOGCS["GCS_MGI",DATUM["D_MGI",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",400000.0],PARAMETER["False_Northing",400000.0],PARAMETER["Central_Meridian",13.3333333333333],PARAMETER["Standard_Parallel_1",49.0],PARAMETER["Standard_Parallel_2",46.0],PARAMETER["Latitude_Of_Origin",47.5],UNIT["Meter",1.0]]</p>
<p>In QGIS, both in the "wrong" and "good" version, please go to the Project Properties dialog, CRS tab, and filter on 31287 and paste the output of the bottom left pane with the WKT and proj.4 strings. It should be:
</p>
<p>WKT
PROJCRS["MGI / Austria Lambert",
BASEGEOGCRS["MGI",
DATUM["Militar-Geographische Institut",
ELLIPSOID["Bessel 1841",6377397.155,299.1528128,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4312]],
CONVERSION["Austria Lambert",
METHOD["Lambert Conic Conformal (2SP)",
ID["EPSG",9802]],
PARAMETER["Latitude of false origin",47.5,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8821]],
PARAMETER["Longitude of false origin",13.3333333333333,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8822]],
PARAMETER["Latitude of 1st standard parallel",49,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8823]],
PARAMETER["Latitude of 2nd standard parallel",46,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8824]],
PARAMETER["Easting at false origin",400000,
LENGTHUNIT["metre",1],
ID["EPSG",8826]],
PARAMETER["Northing at false origin",400000,
LENGTHUNIT["metre",1],
ID["EPSG",8827]]],
CS[Cartesian,2],
AXIS["northing (X)",north,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["easting (Y)",east,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["Topographic mapping (medium and small scale)."],
AREA["Austria."],
BBOX[46.4,9.53,49.02,17.17]],
ID["EPSG",31287]]</p>
<p>Proj4
+proj=lcc +lat_0=47.5 +lon_0=13.3333333333333 +lat_1=49 +lat_2=46 +x_0=400000 +y_0=400000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
</p>
<p>
</p>
<p>If you are doing a reprojection in QGIS, please also go to the Transformations tab in the Project Properties dialog, and show the details of the Datum transformation which is used, in particular the string that starts with "+proj=pipeline"
<style type="text/css">p, li { white-space: pre-wrap; }</style></p>
<p>(not sure what you referred to with the "bottom of this email". an attachment? if so, it didn't get through)
</p>
<p></p>
<p>Even
</p>
<div class="moz-cite-prefix">Le 24/09/2024 à 20:13, Jochen Albrecht
via osgeolive a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CAKRQQ2TbP+dfnapp3XWbRMf2v5Kf6tAAjZoWWfnEUHhO_NH9ow@mail.gmail.com">
<div dir="ltr">
<div>Hi</div>
<div><br>
</div>
<div>I am teaching with the latest version of OSgeo-Live 16,
which runs QGIS 3.28.</div>
<div>I am using for my class a small demo shapefile, which
happens to be in EPSG:31287 (Austria, Styria). QGIS 3.28
places this file into the Mediterranean Sea, off the coast of
Italy/France.</div>
<div>Earlier versions of QGIS (in previous years of teaching),
as well as my current stand-alone version 3.36 are placing the
data correctly into the province of Styria in Austria. The
problem occurs only in version 3.28.</div>
<div><br>
</div>
<div>The problem is aggravated by the fact that any attempt at
transforming the data is rendered futile in this particular
version of QGIS because none of the reprojection dialogues
allow me to specify a transformation. I even tried to define a
new coordinate system using the transformation parameters
given by QGIS 3.38 (see the bottom of this email). Applying
these, places the data again into the Mediterranean.</div>
<div><br>
</div>
<div>I have been teaching with QGIS for over 20 years, and this
is a new one for me. The easy solution would be to just not
use version QGIS 3.28. But it is the one that has been burned
into the current OSgeo-Live distribution. The whole purpose of
OSgeo-Live is to work with versions of the individual software
packages that have proven to play nice with each other and I
do not want to experiment and have students upgrade to a
newer version of GIS.</div>
<div><br>
</div>
<div>At this point, I am at a loss to advise my students what to
do, or to even explain what is happening here. It is as if the
particular coordinate system library for this version has a
bug that nobody has come across before. Am I missing
something?</div>
<div>Last minute update: the same error occurs in the latest
version of ArcGIS Pro 3.3. I have no idea what those two have
in common, but it suggests to me that it has something to do
with a flawed coordinate system library rather than the
application software itself.</div>
<div>Again, any suggestion for how to address this (or even
explain it to my students) would be appreciated.</div>
<div><br>
</div>
<div>Cheers,</div>
<div> Jochen</div>
<div><br>
</div>
<div>
<div dir="ltr" class="gmail_signature"
data-smartmail="gmail_signature">
<div dir="ltr"><span>
<div dir="ltr">
<table>
<colgroup><col><col><col></colgroup><tbody>
<tr>
<td colspan="2">
<p dir="ltr"><span>Dr. Jochen Albrecht, GISP (he/him/his)</span></p>
<p dir="ltr"><span>Professor for Computational and Theoretical Geography</span></p>
<p dir="ltr"><a
href="http://www.geo.hunter.cuny.edu/"
target="_blank" moz-do-not-send="true"><span>Department of Geography and Environmental Science</span></a>
</p>
<p dir="ltr"><span>Hunter College CUNY</span></p>
<p dir="ltr"><span>695 Park Avenue</span></p>
<p dir="ltr"><span>New York, NY 10065</span></p>
<p dir="ltr"><span>Immediate Past-President, </span><a
href="https://www.gisci.org/" target="_blank"
moz-do-not-send="true"><span>GIS Certification Institute</span></a></p>
<img moz-do-not-send="true" width="200"
height="80"><br>
</td>
<td>
<p dir="ltr"><span><span></span></span></p>
<br>
</td>
</tr>
</tbody>
</table>
<span>Recently published: </span><a
href="https://sites.google.com/rojasap.com/gisandhousing/"
target="_blank" moz-do-not-send="true">GIS and
Housing: Principles and Practices</a><b><br>
</b><br>
</div>
</span></div>
</div>
</div>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
osgeolive mailing list
<a class="moz-txt-link-abbreviated" href="mailto:osgeolive@lists.osgeo.org">osgeolive@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/osgeolive">https://lists.osgeo.org/mailman/listinfo/osgeolive</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>