Merge pull request #120 from jjimenezshaw/pt-geoid

Add Portuguese geoid model GeodPT08

Author: rouault

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,,

files.geojson

@@ -271,9 +271,10 @@
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/pl_gugik_geoid2011-PL-EVRF2007-NH.tif", "name": "pl_gugik_geoid2011-PL-EVRF2007-NH.tif", "area_of_use": "Poland - onshore", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:9701", "source_crs_name": "ETRF2000-PL", "target_crs_code": "EPSG:9651", "target_crs_name": "EVRF2007-PL height", "source": "Główny Urząd Geodezji i Kartografii (GUGiK)", "source_country": "Poland", "source_id": "pl_gugik", "source_url": "https://www.gov.pl/web/gugik/", "description": "ETRF2000-PL (EPSG:9701) to EVRF2007-PL height (EPSG:9651)", "file_size": 699407, "sha256sum": "345a7fb62adffab6c8e864b18239d64e519d19b815418c5af522a254630e3182", "version_added": "1.9" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 14.05, 49.0 ], [ 24.2, 49.0 ], [ 24.2, 54.88 ], [ 14.05, 54.88 ], [ 14.05, 49.0 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/pl_gugik_geoid2011-PL-KRON86-NH.tif", "name": "pl_gugik_geoid2011-PL-KRON86-NH.tif", "area_of_use": "Poland - onshore", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:9701", "source_crs_name": "ETRF2000-PL", "target_crs_code": "EPSG:9650", "target_crs_name": "Baltic 1986 height", "source": "Główny Urząd Geodezji i Kartografii (GUGiK)", "source_country": "Poland", "source_id": "pl_gugik", "source_url": "https://www.gov.pl/web/gugik/", "description": "ETRF2000-PL (EPSG:9701) to Baltic 1986 height (EPSG:9650)", "file_size": 661511, "sha256sum": "16bb69272fdce5e812f7d1c81facd3fb0bf1489e298cbdf82abfe7e73690daa8", "version_added": "1.9" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 14.05, 49.0 ], [ 24.2, 49.0 ], [ 24.2, 54.88 ], [ 14.05, 54.88 ], [ 14.05, 49.0 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/pl_gugik_geoid2021-PL-EVRF2007-NH.tif", "name": "pl_gugik_geoid2021-PL-EVRF2007-NH.tif", "area_of_use": "Poland - onshore", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:9701", "source_crs_name": "ETRF2000-PL", "target_crs_code": "EPSG:9651", "target_crs_name": "EVRF2007-PL height", "source": "Główny Urząd Geodezji i Kartografii (GUGiK)", "source_country": "Poland", "source_id": "pl_gugik", "source_url": "https://www.gov.pl/web/gugik/", "description": "ETRF2000-PL (EPSG:9701) to EVRF2007-PL height (EPSG:9651)", "file_size": 1033566, "sha256sum": "8b806d084296c0b8070a9cc9ad6d28218141e089e3ff0573254870447f361788", "version_added": "1.15" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 13.0, 48.0 ], [ 25.0, 48.0 ], [ 25.0, 56.0 ], [ 13.0, 56.0 ], [ 13.0, 48.0 ] ] ] } },
-{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/pt_dgt_README.txt", "name": "pt_dgt_README.txt", "source": "DG Territorio", "source_country": "Portugal", "source_id": "pt_dgt", "source_url": "http://www.dgterritorio.pt/", "file_size": 1571, "sha256sum": "d6bd95a83240f94e1655c6dd5cb1e7ab2eb2c53d59294986a60c38ee8221b27d" }, "geometry": null },
+{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/pt_dgt_README.txt", "name": "pt_dgt_README.txt", "source": "DG Territorio", "source_country": "Portugal", "source_id": "pt_dgt", "source_url": "http://www.dgterritorio.pt/", "file_size": 2403, "sha256sum": "b727d2a6b170286042fa015f5978b7627cb47ba331dace8ca631b99af78bf1ac" }, "geometry": null },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/pt_dgt_D73_ETRS89_geo.tif", "name": "pt_dgt_D73_ETRS89_geo.tif", "area_of_use": "Portugal", "type": "HORIZONTAL_OFFSET", "source_crs_code": "EPSG:4274", "source_crs_name": "Datum 73", "target_crs_code": "EPSG:4258", "target_crs_name": "ETRS89", "source": "DG Territorio", "source_country": "Portugal", "source_id": "pt_dgt", "source_url": "http://www.dgterritorio.pt/", "description": "Datum 73 (EPSG:4274) to ETRS89 (EPSG:4258)", "file_size": 222863, "sha256sum": "802c7f4702e2fcf68b8720028182513efdd2af3d5154487ac2e225c188e02e07" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ -9.9305556, 36.7638889 ], [ -5.7505556, 36.7638889 ], [ -5.7505556, 42.3638889 ], [ -9.9305556, 42.3638889 ], [ -9.9305556, 36.7638889 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/pt_dgt_DLx_ETRS89_geo.tif", "name": "pt_dgt_DLx_ETRS89_geo.tif", "area_of_use": "Portugal", "type": "HORIZONTAL_OFFSET", "source_crs_code": "EPSG:4207", "source_crs_name": "Lisbon", "target_crs_code": "EPSG:4258", "target_crs_name": "ETRS89", "source": "DG Territorio", "source_country": "Portugal", "source_id": "pt_dgt", "source_url": "http://www.dgterritorio.pt/", "description": "Lisbon (EPSG:4207) to ETRS89 (EPSG:4258)", "file_size": 217126, "sha256sum": "a2581158f257b86712d884d003a3d55f013f5fed07a2afc9865a35b08e462d9a" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ -9.9305556, 36.7638889 ], [ -5.7505556, 36.7638889 ], [ -5.7505556, 42.3638889 ], [ -9.9305556, 42.3638889 ], [ -9.9305556, 36.7638889 ] ] ] } },
+{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/pt_dgt_GeodPT08.tif", "name": "pt_dgt_GeodPT08.tif", "area_of_use": "Portugal - mainland - onshore.", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4937", "source_crs_name": "ETRS89", "target_crs_code": "EPSG:5780", "target_crs_name": "Cascais height", "source": "DG Territorio", "source_country": "Portugal", "source_id": "pt_dgt", "source_url": "http://www.dgterritorio.pt/", "description": "ETRS89 (EPSG:4937) to Cascais height (EPSG:5780)", "file_size": 81638, "sha256sum": "8c006ae3afa037c214d30e7ec97736452705773ccb663274e003ced1d08ccaa0", "version_added": "1.17" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ -10.0125, 36.5125 ], [ -5.9875, 36.5125 ], [ -5.9875, 42.2375 ], [ -10.0125, 42.2375 ], [ -10.0125, 36.5125 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/se_lantmateriet_README.txt", "name": "se_lantmateriet_README.txt", "source": "Lantmäteriet", "source_country": "Sweden", "source_id": "se_lantmateriet", "source_url": "https://www.lantmateriet.se/", "file_size": 840, "sha256sum": "3cc363d80243ed9a464e243706df93f6ef0fbb39e1bb372a70851267a974164f" }, "geometry": null },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/se_lantmateriet_SWEN17_RH2000.tif", "name": "se_lantmateriet_SWEN17_RH2000.tif", "area_of_use": "Sweden", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4977", "source_crs_name": "SWEREF99", "target_crs_code": "EPSG:5613", "target_crs_name": "RH2000 height", "source": "Lantmäteriet", "source_country": "Sweden", "source_id": "se_lantmateriet", "source_url": "https://www.lantmateriet.se/", "description": "SWEREF99 (EPSG:4977) to RH2000 height (EPSG:5613)", "file_size": 2317290, "sha256sum": "c5160019c9d2f64e0379e7725be795952b5005759e6f51e5322f6607af92a5a3" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 10.0, 54.0 ], [ 25.0, 54.0 ], [ 25.0, 70.0 ], [ 10.0, 70.0 ], [ 10.0, 54.0 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/si_gurs_README.txt", "name": "si_gurs_README.txt", "source": "Geodetska uprava Republika Slovenija (GURS)", "source_country": "Slovenia", "source_id": "si_gurs", "source_url": "http://www.gu.gov.si", "file_size": 802, "sha256sum": "0cfa89cafe012ed043ac1054b999b98070eb4313c300b28d10e8294d00dc80a7", "version_added": "1.13" }, "geometry": null },

index.html

@@ -822,9 +822,10 @@ <h2>Content</h2>
 <li><a style="word-break: break-word" href="pl_gugik_geoid2011-PL-KRON86-NH.tif">pl_gugik_geoid2011-PL-KRON86-NH.tif</a> - Poland - onshore - ETRF2000-PL (EPSG:9701) to Baltic 1986 height (EPSG:9650). Last modified: 2021-11-20</li>
 <li><a style="word-break: break-word" href="pl_gugik_geoid2021-PL-EVRF2007-NH.tif">pl_gugik_geoid2021-PL-EVRF2007-NH.tif</a> - Poland - onshore - ETRF2000-PL (EPSG:9701) to EVRF2007-PL height (EPSG:9651). Last modified: 2023-07-10</li>
 </ul><hr><h3><a href="http://www.dgterritorio.pt/">DG Territorio</a></h3><ul>
-<li><a style="word-break: break-word" href="pt_dgt_README.txt">pt_dgt_README.txt</a>. Last modified: 2020-01-24</li>
+<li><a style="word-break: break-word" href="pt_dgt_README.txt">pt_dgt_README.txt</a>. Last modified: 2024-02-25</li>
 <li><a style="word-break: break-word" href="pt_dgt_D73_ETRS89_geo.tif">pt_dgt_D73_ETRS89_geo.tif</a> - Portugal - Datum 73 (EPSG:4274) to ETRS89 (EPSG:4258). Last modified: 2020-01-24</li>
 <li><a style="word-break: break-word" href="pt_dgt_DLx_ETRS89_geo.tif">pt_dgt_DLx_ETRS89_geo.tif</a> - Portugal - Lisbon (EPSG:4207) to ETRS89 (EPSG:4258). Last modified: 2020-01-24</li>
+<li><a style="word-break: break-word" href="pt_dgt_GeodPT08.tif">pt_dgt_GeodPT08.tif</a> - Portugal - mainland - onshore. - ETRS89 (EPSG:4937) to Cascais height (EPSG:5780). Last modified: 2024-02-25</li>
 </ul><hr><h3><a href="https://www.lantmateriet.se/">Lantmäteriet</a></h3><ul>
 <li><a style="word-break: break-word" href="se_lantmateriet_README.txt">se_lantmateriet_README.txt</a>. Last modified: 2020-01-24</li>
 <li><a style="word-break: break-word" href="se_lantmateriet_SWEN17_RH2000.tif">se_lantmateriet_SWEN17_RH2000.tif</a> - Sweden - SWEREF99 (EPSG:4977) to RH2000 height (EPSG:5613). Size: 2.2 MB. Last modified: 2020-01-24</li>

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