<div dir="ltr"><div>Italy? I am not aware of any Geoid grid file for Italy in PROJ-Data: <a href="https://cdn.proj.org/">https://cdn.proj.org/</a></div><div><br></div><div>EPSG:5214 is the vertical coordinate reference system, not a geoid model. <a href="https://epsg.org/crs_5214/Genoa-1942-height.html">https://epsg.org/crs_5214/Genoa-1942-height.html</a></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 17 Oct 2022 at 10:59, Diogo <<a href="mailto:diogo@photopills.com">diogo@photopills.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"><div dir="ltr">Hi Andrew. Thanks for taking the time to read my email and reply it. <br></div><div dir="ltr"><div>I'm
 aware of Copernicus land and other datasets like SMRT3 and ASTER 
however I need the best possible elevation resolution available. To do 
so, I'm getting the elevation dataset directly from countries geographic
 institutes and converting them into the coordinate system 4326+3855. 
But as described in my previous email, I'm having some trouble with some
 geoids. Do you have an idea why I'm having these transformation issues?
 To be honest, I'm a bit lost because I can't understand it.<br></div><div><br></div><div>Thanks again for your help,</div><div>Diogo</div><div><br></div></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 14 Oct 2022 at 15:57, Andrew C Aitchison <<a href="mailto:andrew@aitchison.me.uk" target="_blank">andrew@aitchison.me.uk</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">On Fri, 14 Oct 2022, Diogo wrote:<br>
<br>
> Hello everyone!<br>
><br>
> I'm merging multiple elevation datasets from different countries and I<br>
> would like to re-project all into the same coordinates system, 4326+3855<br>
> using gdalwarp. I'm using `gdaltransform` to confirm that there is in fact<br>
> a conversion from the original projection to the aspected one. However, for<br>
> some datasets, there is only conversion in the XY axis, while the Z output<br>
> is the same as the input.<br>
<br>
Is there a good reason to start with these particular datasets ?<br>
<br>
There are international bodies who have produced elevation datasets covering the whole of <br>
Europe, eg<br>
<a href="https://land.copernicus.eu/user-corner/how-to-access-our-data" rel="noreferrer" target="_blank">https://land.copernicus.eu/user-corner/how-to-access-our-data</a><br>
<a href="https://www.eea.europa.eu/data-and-maps/data#c0=5&c11=&c5=all&b_start=0&c12=elevation" rel="noreferrer" target="_blank">https://www.eea.europa.eu/data-and-maps/data#c0=5&c11=&c5=all&b_start=0&c12=elevation</a><br>
<a href="https://www.eea.europa.eu/data-and-maps/data/copernicus-land-monitoring-service-eu-dem" rel="noreferrer" target="_blank">https://www.eea.europa.eu/data-and-maps/data/copernicus-land-monitoring-service-eu-dem</a><br>
<br>
which would save you trying to stitch together data in different CRSes.<br>
<br>
> Below there are some examples (the geoid code is extracted from the dataset<br>
> documentation)<br>
><br>
> Genoa 1942 height (Italy)<br>
><br>
> $ echo "14.995 37.755 3357" | PROJ_NETWORK=ON gdaltransform --debug on<br>
> -s_srs "EPSG:4326+3855" -t_srs "EPSG:32632+5214">> 1028333.73796622<br>
> 4195591.29690422 3357<br>
><br>
> N2000 (Finland)<br>
><br>
> $ echo "25.15296 62.50195 269" | PROJ_NETWORK=ON gdaltransform --debug<br>
> on -s_srs "EPSG:4326+3855" -t_srs "EPSG:3067+3900">> 404861.868388335<br>
> 6931460.32557422 269<br>
><br>
> EVRF2007 (Estonia)<br>
><br>
> $ echo "27.059167 57.714444 318" | PROJ_NETWORK=ON gdaltransform<br>
> --debug on -s_srs "EPSG:4326+3855" -t_srs "EPSG:3301+5621">><br>
> 682284.606829111 6401087.4159131 318<br>
><br>
><br>
> I don't have much experience with the projections and I'm a bit lost<br>
> to understand the possible issues, or interpreted the Z output. What<br>
> can be the possible issues?<br>
><br>
> Any help would be appreciated.<br>
><br>
> Best regards,<br>
><br>
> Diogo<br>
><br>
<br>
-- <br>
Andrew C. Aitchison                      Kendal, UK<br>
                    <a href="mailto:andrew@aitchison.me.uk" target="_blank">andrew@aitchison.me.uk</a><br>
</blockquote></div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div>