add GeodPT08 for Portugal

Author: jjimenezshaw

copyright_and_licenses.csv

@@ -273,6 +273,7 @@ pl_gugik_geoid2021-PL-EVRF2007-NH.tif,Główny Urząd Geodezji i Kartografii,CC-
 pl_gugik_README.txt,Disclaimed,Public domain,1.9,
 pt_dgt_D73_ETRS89_geo.tif,Direcao Geral de Territorio,CC-BY-4.0,,
 pt_dgt_DLx_ETRS89_geo.tif,Direcao Geral de Territorio,CC-BY-4.0,,
+pt_dgt_GeodPT08.tif,Direcao Geral de Territorio,CC-BY-4.0,1.17,
 pt_dgt_README.txt,Disclaimed,Public domain,,
 README.DATA,Disclaimed,Public domain,,
 se_lantmateriet_README.txt,Disclaimed,Public domain,,

pt_dgt/build_GeodPT08.sh

@@ -0,0 +1,42 @@
+#!/bin/bash -e
+
+# https://www.dgterritorio.gov.pt/geodesia/modelo-geoide
+# Download https://www.dgterritorio.gov.pt/sites/default/files/ficheiros-geodesia/GeodPT08.dat
+
+# Usage:
+# PROJ_DATA_DIR=/path/to/PROJ-data ./build_GeodPT08.sh
+
+def_data_dir=$(dirname $0)/..
+
+PROJ_DATA_DIR="${PROJ_DATA_DIR:-$def_data_dir}"
+
+orig=GeodPT08.dat
+
+# Setup build directory
+mkdir -p build
+
+tmp=./build/GeodPT08.xyz
+sort -s -n -k1,1 -k2,2r "$orig" > $tmp
+input=$tmp
+output=pt_dgt_GeodPT08.tif
+
+docker run --user $(id -u):$(id -g) --workdir $PWD \
+            --rm -v /home:/home ghcr.io/osgeo/gdal:alpine-normal-latest \
+            sh -c " \
+            # Call vertoffset_grid_to_gtiff-script
+            python3 ${PROJ_DATA_DIR}/grid_tools/vertoffset_grid_to_gtiff.py \
+            --type GEOGRAPHIC_TO_VERTICAL \
+            --parameter-name geoid_undulation \
+            --source-crs \"EPSG:4937\" \
+            --target-crs \"EPSG:5780\" \
+            --description \"ETRS89 (EPSG:4937) to Cascais height (EPSG:5780). Converted from $orig\" \
+            --area-of-use \"Portugal - mainland - onshore.\" \
+            --copyright \"Derived from work DGT+FCUL. CC-BY-4.0 https://creativecommons.org/licenses/by/4.0/\" \
+            ./$input ./$output && \
+            # Run check_gtiff_grid-script
+            python3 ${PROJ_DATA_DIR}/grid_tools/check_gtiff_grid.py ./$output \
+            # Show info
+            gdalinfo ./$output "
+
+# Remove build directory
+rm -rf build

pt_dgt/pt_dgt_GeodPT08.tif

@@ -34,3 +34,21 @@ Residuals at 130 further test points average 0.06m, maximum 0.16m.
 Recommended to transform from Portuguese Datum 73 (EPSG:27493) to ETRS89 in Portugal (EPSG:3763).
 
 * pt_dgt_D73_ETRS89_geo.tif
+
+### Portugal: Vertical grid GeodPT08
+
+*Source*: [DGT](https://www.dgterritorio.gov.pt/geodesia/modelo-geoide)  
+*Format*: GeoTIFF converted from XYZ (adjusted)  
+*License*: [Creative Commons Attribution 4.0](http://creativecommons.org/licenses/by/4.0/)  
+*Credit*: (c)2024 Geoid for Portugal Mainland by DGT+FCUL  
+*Horizontal CRS*: EPSG:4258 (ETRS89)  
+*Accuracy*: 0.04 m
+
+Built based on the ICAGM07 gravimetric geoid model using 137 leveling marks and 1020 geodetic vertices.
+Estimated global accuracy of 4 cm, determined with reference to the continent's geodetic vertex and geometric leveling networks.
+
+Vertical transformation (EPSG:10544) for Geoid model GeodPT08.
+Used to make the transition from heights in vertical CRS Cascais height (EPSG:5780)
+to heights abobe the ellipsoid in ETRS89 (EPSG:4937)
+
+* pt_dgt_GeodPT08.tif

pt_dgt/pt_dgt_README.txt

@@ -270,6 +270,7 @@ pl_gugik_geoid2021-PL-EVRF2007-NH.tif
 pl_gugik_README.txt
 pt_dgt_D73_ETRS89_geo.tif
 pt_dgt_DLx_ETRS89_geo.tif
+pt_dgt_GeodPT08.tif
 pt_dgt_README.txt
 README.DATA
 se_lantmateriet_README.txt