<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.Shkpostityyli18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 2.0cm 70.85pt 2.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="FI" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Maybe there is no transformation for converting heights from N2000 into WGS84 heights registered into the EPSG database. I found only two N2000 related transformations EPSG:9573 and
EPSG:9574 and they seem convert to EVRF2019 height.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">Proj can convert between our national height systems N2000, N60, and N43 but I fear that we at the National Land Survey of Finland have not managed to register them yet.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">This publication is too heavy to read for me but please have a look and contact Mirjam if you need more advice
<a href="https://www.researchgate.net/publication/256481988_Development_of_the_Finnish_Height_Conversion_Surface_FIN2005N00#read">
https://www.researchgate.net/publication/256481988_Development_of_the_Finnish_Height_Conversion_Surface_FIN2005N00#read</a><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US">-Jukka Rahkonen-<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> gdal-dev <gdal-dev-bounces@lists.osgeo.org>
<b>Puolesta </b>Diogo<br>
<b>Lähetetty:</b> maanantai 17. lokakuuta 2022 11.59<br>
<b>Vastaanottaja:</b> gdal-dev@lists.osgeo.org<br>
<b>Aihe:</b> Re: [gdal-dev] Transformations with same horizontal datum but different geoid outputs the same Z value<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Hi Andrew. Thanks for taking the time to read my email and reply it.
<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">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.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks again for your help,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Diogo<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">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:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<p class="MsoNormal">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://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fland.copernicus.eu%2Fuser-corner%2Fhow-to-access-our-data&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Caa18fc3bad0649824f9708dab01ed49b%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638015943753338265%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000%7C%7C%7C&sdata=mfo6NuPABxZJkU3MRQeYfLRAsV0TeX%2FlNkAlinnKN8w%3D&reserved=0" target="_blank">https://land.copernicus.eu/user-corner/how-to-access-our-data</a><br>
<a href="https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.eea.europa.eu%2Fdata-and-maps%2Fdata%23c0%3D5%26c11%3D%26c5%3Dall%26b_start%3D0%26c12%3Delevation&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Caa18fc3bad0649824f9708dab01ed49b%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638015943753338265%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000%7C%7C%7C&sdata=1b7YhCgROWYUTq8TBXIivDLKmzEHuj5NpNRCW%2BeVQ%2Fg%3D&reserved=0" 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://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.eea.europa.eu%2Fdata-and-maps%2Fdata%2Fcopernicus-land-monitoring-service-eu-dem&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Caa18fc3bad0649824f9708dab01ed49b%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638015943753338265%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000%7C%7C%7C&sdata=PKBB6EfcZ5qniFw6gKkcENOzoPvnfTiQr53ErOZV4DA%3D&reserved=0" 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><o:p></o:p></p>
</blockquote>
</div>
</div>
</div>
</body>
</html>