<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 16, 2020 at 5:15 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.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"><u></u>
<div style="font-family:"monospace";font-size:9pt;font-weight:400;font-style:normal">
<p style="margin:0px;text-indent:0px">> Well, when providing that coordinate, I got them from a sample Dataset in</p>
<p style="margin:0px;text-indent:0px">> EPSG:31255 so they should be valid.</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> Weird thing:</p>
<p style="margin:0px;text-indent:0px">> cs2cs "EPSG:31255" "EPSG:4326"</p>
<p style="margin:0px;text-indent:0px">> 88000 273000</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> 45d52'39.123"N 16d50'57.28"E 0.000 (as you said)</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> However:</p>
<p style="margin:0px;text-indent:0px">> gdaltransform -s_srs EPSG:31255 -t_srs EPSG:4326</p>
<p style="margin:0px;text-indent:0px">> 88000 273000</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> 14.5025684938377 47.5898013589131 46.4711159374565</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> And same thing if I draw a point in QGIS: SRID=31255;POINT(88000 273000)</p>
<p style="margin:0px;text-indent:0px">> and I switch the map projection to 4326, that rendered point on the map has</p>
<p style="margin:0px;text-indent:0px">> the same coordinates reported by gdaltransform.</p>
<p style="margin:0px;text-indent:0px">> </p>
<p style="margin:0px;text-indent:0px">> What do you think about it?</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Ah, the reason is axis order. EPSG:31255 uses Northing,Easting, and cs2cs honours it. Whereas gdaltransform uses the "traditional GIS axis order", ie easting, northing.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">So with PROJ master:</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">$ echo 273000 88000 600 | PROJ_NETWORK=ON src/cs2cs "EPSG:31255+5778" "EPSG:3035+9274"</p>
<p style="margin:0px;text-indent:0px">2729881.62        4659643.84 599.75</p></div></blockquote><div><br></div><div>Ah. That makes sense. Thanks.</div><div>Final question (sorry about that):</div><div>Are you aware of any involvement of this kind of reprojection in gdal_warp processing?</div><div>I mean, I just did a test by building PROJ master and GDAL 3.1 against that PROJ version and tried a <br></div><div><br></div><div>gdalwarp -s_srs "EPSG:31255+5778" -t_srs "EPSG:3035+9274" dem.tif warpeddem.tif</div><div><br></div><div>Is it expected to reproject the DEM values? :-)</div><div>Asking that since I'm trying to figure out a way to see that vertical transformation used in some practical use-case with DEM raster but I didn't get the pixel values modified so I probably went too far and that's not how it should be used.</div><div><br></div><div>Thanks again for your feedbacks so far.</div><div>Regards,</div><div>Daniele</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="font-family:"monospace";font-size:9pt;font-weight:400;font-style:normal">
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">With EPSG:5621 "EVRF2007 height" instead of EPSG:5778 "EVRF2000 height", I don't get any change in Z, because there's no registered transformation between EPSG:5778 "GHA height" and it.</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">Even</p>
<p style="margin:0px;text-indent:0px"> </p>
<p style="margin:0px;text-indent:0px">-- </p>
<p style="margin:0px;text-indent:0px">Spatialys - Geospatial professional services</p>
<p style="margin:0px;text-indent:0px"><a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a></p></div></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Regards,</div><div dir="ltr">Daniele Romagnoli</div><div dir="ltr">==</div><div dir="ltr">GeoServer Professional Services from the experts! Visit <a href="http://goo.gl/it488V" target="_blank">http://goo.gl/it488V</a> for more information.</div><div dir="ltr">==</div><div dir="ltr"><br></div><div dir="ltr">Ing. Daniele Romagnoli</div><div dir="ltr">Senior Software Engineer</div><div dir="ltr"><br></div><div dir="ltr">GeoSolutions S.A.S.</div><div dir="ltr">Via di Montramito 3/A</div><div dir="ltr">55054  Massarosa (LU)</div><div dir="ltr">Italy</div><div dir="ltr">phone: +39 0584 962313</div><div dir="ltr">fax:      +39 0584 1660272</div><div dir="ltr"><br></div><div dir="ltr"><a href="http://www.geo-solutions.it" target="_blank">http://www.geo-solutions.it</a></div><div dir="ltr"><a href="http://twitter.com/geosolutions_it" target="_blank">http://twitter.com/geosolutions_it</a></div><div dir="ltr"><br></div><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">-------------------------------------------------------</span></div><div dir="ltr"><span style="font-size:12.8px"><br></span></div><div dir="ltr"><span style="font-size:12.8px">Con riferimento alla normativa sul trattamento dei dati personali (Reg. UE 2016/679 - Regolamento generale sulla protezione dei dati “GDPR”), si precisa che ogni circostanza inerente alla presente email (il suo contenuto, gli eventuali allegati, etc.) è un dato la cui conoscenza è riservata al/i solo/i destinatario/i indicati dallo scrivente. Se il messaggio Le è giunto per errore, è tenuta/o a cancellarlo, ogni altra operazione è illecita. Le sarei comunque grato se potesse darmene notizia.</span></div><div dir="ltr"><span style="font-size:12.8px"><br></span></div><div dir="ltr"><span style="font-size:12.8px">This email is intended only for the person or entity to which it is addressed and may contain information that is privileged, confidential or otherwise protected from disclosure. We remind that - as provided by European Regulation 2016/679 “GDPR” - copying, dissemination or use of this e-mail or the information herein by anyone other than the intended recipient is prohibited. If you have received this email by mistake, please notify us immediately by telephone or e-mail.</span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>