{"id":686,"date":"2022-01-27T15:48:57","date_gmt":"2022-01-27T14:48:57","guid":{"rendered":"https:\/\/sigeo.cerege.fr\/?p=686"},"modified":"2026-03-03T14:44:59","modified_gmt":"2026-03-03T13:44:59","slug":"convert-raster-from-ellispoidal-height-to-orthometric-altitude","status":"publish","type":"post","link":"https:\/\/sigeo.cerege.fr\/?p=686","title":{"rendered":"Convert raster from ellispoidal height to orthometric altitude"},"content":{"rendered":"\n<p>Differences between ellipsoidal height and orthometric altitude (or geoid height) are explained here https:\/\/spatialthoughts.com\/2019\/10\/26\/convert-between-orthometric-and-ellipsoidal-elevations-using-gdal\/<\/p>\n\n\n\n<p>The definition of the projection is the PROJ4 line that can be found on <a href=\"https:\/\/spatialreference.org\/\" target=\"_blank\" rel=\"noreferrer noopener\">https:\/\/spatialreference.org\/<\/a><\/p>\n\n\n\n<h3 class=\"wp-block-heading\">examples<\/h3>\n\n\n\n<ul class=\"wp-block-list\">\n<li>using gdalwarp to convert from UTM33N-WGS84 ellipsoidal height to UTM33N-EGM96 geoid height<\/li>\n<\/ul>\n\n\n\n<p><code>gdalwarp -s_srs \"+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs\" -t_srs \"+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs +geoidgrids=egm96_15.gtx\" input.tif output.tif<\/code><\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>from Lambert93 planimetric and NGF-IGN69 (with RAF20) to WGS84 ellipsoidal heights<\/li>\n<\/ul>\n\n\n\n<p><code>gdalwarp -s_srs \"+proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=49 +lat_2=44 +x_0=700000 +y_0=6600000 +ellps=GRS80 +units=m +no_defs +type=crs +geoidgrids=RAF20.gtx\" -t_srs \"+proj=longlat +datum=WGS84 +no_defs +type=crs\" input.tif output.tif<\/code><\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>from altitude to ellipsoidal height use such command (change srs and geoid grid if needed): <\/li>\n<\/ul>\n\n\n\n<p><code>gdalwarp -s_srs \"+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs <code>+geoidgrids=egm96_15.gtx\"<\/code> -t_srs \"+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs\" input.tif output.tif<\/code><\/p>\n\n\n\n<p>you can choose your input and output spatial reference using the proj4 syntax https:\/\/epsg.io\/32633 , and note the add of the geoid file with <em>+geoidgrids=egm96_15.gtx<\/em> in the command.<\/p>\n\n\n\n<p>if you want to use EGM2008 geoid, get the file from https:\/\/github.com\/OSGeo\/proj-datumgrid\/blob\/master\/world\/egm08_25.gtx and copy it to the proj directory. In Linux it&#8217;s <\/p>\n\n\n\n<p>\/<code>usr\/share\/proj<\/code><\/p>\n\n\n\n<p>then you can simply change the geoid file in the command <em>+geoidgrids=egm08_25.gtx<\/em><\/p>\n\n\n\n<h3 class=\"wp-block-heading\">examples of proj definitions<\/h3>\n\n\n\n<p>WGS84 ellipsoid with EGM96 elevation<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>+proj=longlat +datum=WGS84 +no_defs +geoidgrids=egm96_15.gtx<\/code><\/pre>\n\n\n\n<p>WGS84 ellipsoid with WGS84 elevation<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>+proj=longlat +datum=WGS84 +no_def<\/code><\/pre>\n\n\n\n<p>RGF93\/Lambert 93 with NGF-IGN69 elevation<\/p>\n\n\n\n<p><code>\"+proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=49 +lat_2=44 +x_0=700000 +y_0=6600000 +ellps=GRS80 +units=m +no_defs +type=crs +geoidgrids=RAF20.gtx\"<\/code><\/p>\n\n\n\n<p><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Differences between ellipsoidal height and orthometric altitude (or geoid height) are explained here https:\/\/spatialthoughts.com\/2019\/10\/26\/convert-between-orthometric-and-ellipsoidal-elevations-using-gdal\/ The definition of the projection is the PROJ4 line that can be found on https:\/\/spatialreference.org\/ examples gdalwarp -s_srs &#8220;+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs&#8221; -t_srs &#8220;+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs +geoidgrids=egm96_15.gtx&#8221; input.tif output.tif gdalwarp -s_srs &#8220;+proj=lcc +lat_0=46.5 +lon_0=3 +lat_1=49 +lat_2=44 +x_0=700000 +y_0=6600000 &hellip; <a href=\"https:\/\/sigeo.cerege.fr\/?p=686\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;Convert raster from ellispoidal height to orthometric altitude&#8221;<\/span><\/a><\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[12,7,1],"tags":[],"class_list":["post-686","post","type-post","status-publish","format-standard","hentry","category-gdal-ogr","category-geo-traitements","category-sigeo"],"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/686","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=686"}],"version-history":[{"count":8,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/686\/revisions"}],"predecessor-version":[{"id":983,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/686\/revisions\/983"}],"wp:attachment":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=686"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=686"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=686"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}