Merge pull request #105 from AlexBass05/add-polish-geoid2021

Add Polish geoid model PL-geoid2021

copyright_and_licenses.csv

@@ -256,6 +256,7 @@ nz_linz_taraht1970-nzvd2016.tif,Land Information New Zealand,CC-BY-4.0,,
 nz_linz_wellht1953-nzvd2016.tif,Land Information New Zealand,CC-BY-4.0,,
 pl_gugik_geoid2011-PL-EVRF2007-NH.tif,Główny Urząd Geodezji i Kartografii,CC-BY-4.0,1.9,
 pl_gugik_geoid2011-PL-KRON86-NH.tif,Główny Urząd Geodezji i Kartografii,CC-BY-4.0,1.9,
+pl_gugik_geoid2021-PL-EVRF2007-NH.tif,Główny Urząd Geodezji i Kartografii,CC-BY-4.0,1.15,
 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,,

files.geojson

@@ -255,9 +255,10 @@
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/nz_linz_stisht1977-nzvd2016.tif", "name": "nz_linz_stisht1977-nzvd2016.tif", "area_of_use": "New Zealand - Stewart Island", "type": "VERTICAL_OFFSET_VERTICAL_TO_VERTICAL", "source_crs_code": "EPSG:7839", "source_crs_name": "NZVD2016 height", "target_crs_code": "EPSG:5772", "target_crs_name": "Stewart Island 1977 height", "source": "Land Information New Zealand", "source_country": "New Zealand", "source_id": "nz_linz", "source_url": "https://www.linz.govt.nz/", "description": "NZVD2016 height (EPSG:7839) to Stewart Island 1977 height (EPSG:5772)", "file_size": 1206, "sha256sum": "f1018025ee18707b9330f9f73efaba72028b4c7dd98341cce0f37312f225efe3" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 167.2, -47.5 ], [ 168.8, -47.5 ], [ 168.8, -46.5 ], [ 167.2, -46.5 ], [ 167.2, -47.5 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/nz_linz_taraht1970-nzvd2016.tif", "name": "nz_linz_taraht1970-nzvd2016.tif", "area_of_use": "New Zealand - Taranaki", "type": "VERTICAL_OFFSET_VERTICAL_TO_VERTICAL", "source_crs_code": "EPSG:7839", "source_crs_name": "NZVD2016 height", "target_crs_code": "EPSG:5769", "target_crs_name": "Taranaki 1970 height", "source": "Land Information New Zealand", "source_country": "New Zealand", "source_id": "nz_linz", "source_url": "https://www.linz.govt.nz/", "description": "NZVD2016 height (EPSG:7839) to Taranaki 1970 height (EPSG:5769)", "file_size": 7040, "sha256sum": "6b0018eb2e0e1c0f413707c6a9116b4377f8fe0b7bcbfb1f94239e3372d627cf" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 173.6, -41.1 ], [ 176.4, -41.1 ], [ 176.4, -38.3 ], [ 173.6, -38.3 ], [ 173.6, -41.1 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/nz_linz_wellht1953-nzvd2016.tif", "name": "nz_linz_wellht1953-nzvd2016.tif", "area_of_use": "New Zealand - Wellington", "type": "VERTICAL_OFFSET_VERTICAL_TO_VERTICAL", "source_crs_code": "EPSG:7839", "source_crs_name": "NZVD2016 height", "target_crs_code": "EPSG:5770", "target_crs_name": "Wellington 1953 height", "source": "Land Information New Zealand", "source_country": "New Zealand", "source_id": "nz_linz", "source_url": "https://www.linz.govt.nz/", "description": "NZVD2016 height (EPSG:7839) to Wellington 1953 height (EPSG:5770)", "file_size": 5845, "sha256sum": "7dfd7f545692301956338387982f5ba4e65fb41217883486347d0965cea04824" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 174.4, -41.6 ], [ 176.4, -41.6 ], [ 176.4, -39.1 ], [ 174.4, -39.1 ], [ 174.4, -41.6 ] ] ] } },
-{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/pl_gugik_README.txt", "name": "pl_gugik_README.txt", "source": "Główny Urząd Geodezji i Kartografii (GUGiK)", "source_country": "Poland", "source_id": "pl_gugik", "source_url": "https://www.gov.pl/web/gugik/", "file_size": 823, "sha256sum": "132f1fe2334d17eeec031246e6259ab017e251fbc12defd205bbf45524ebfae8", "version_added": "1.9" }, "geometry": null },
+{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/pl_gugik_README.txt", "name": "pl_gugik_README.txt", "source": "Główny Urząd Geodezji i Kartografii (GUGiK)", "source_country": "Poland", "source_id": "pl_gugik", "source_url": "https://www.gov.pl/web/gugik/", "file_size": 1070, "sha256sum": "8b03a28591398872146786f1dfef0938c885d2d26d31625dab6569f3bcaf58b5", "version_added": "1.9" }, "geometry": null },
 { "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_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 ] ] ] } },

index.html

@@ -465,8 +465,8 @@ <h2>Content</h2>
 <ul>
 <li><a href="README.DATA">README.DATA</a></li>
 </ul><hr><h3><a href="https://www.ign.gob.ar">Instituto Geográfico Nacional (IGN)</a></h3><ul>
-<li><a style="word-break: break-word" href="ar_ign_README.txt">ar_ign_README.txt</a>. Last modified: 2023-06-20</li>
-<li><a style="word-break: break-word" href="ar_ign_GEOIDE-Ar16.tif">ar_ign_GEOIDE-Ar16.tif</a> - Argentina - onshore - POSGAR 2007 (EPSG:5342) to SRVN16 height (EPSG:9255). Size: 6.1 MB. Last modified: 2023-06-20</li>
+<li><a style="word-break: break-word" href="ar_ign_README.txt">ar_ign_README.txt</a>. Last modified: 2023-06-23</li>
+<li><a style="word-break: break-word" href="ar_ign_GEOIDE-Ar16.tif">ar_ign_GEOIDE-Ar16.tif</a> - Argentina - onshore - POSGAR 2007 (EPSG:5342) to SRVN16 height (EPSG:9255). Size: 6.1 MB. Last modified: 2023-06-23</li>
 </ul><hr><h3><a href="http://www.bev.gv.at/portal/page?_pageid=713,2157075&amp;_dad=portal&amp;_schema=PORTAL">Austria Bundesamt für Eich- und Vermessungswessen</a></h3><ul>
 <li><a style="word-break: break-word" href="at_bev_README.txt">at_bev_README.txt</a>. Last modified: 2022-11-08</li>
 <li><a style="word-break: break-word" href="at_bev_AT_GIS_GRID.tif">at_bev_AT_GIS_GRID.tif</a> - Austria - MGI (EPSG:4312) to ETRS89 (EPSG:4258). Last modified: 2020-01-24</li>
@@ -767,9 +767,10 @@ <h2>Content</h2>
 <li><a style="word-break: break-word" href="nz_linz_taraht1970-nzvd2016.tif">nz_linz_taraht1970-nzvd2016.tif</a> - New Zealand - Taranaki - NZVD2016 height (EPSG:7839) to Taranaki 1970 height (EPSG:5769). Last modified: 2020-02-02</li>
 <li><a style="word-break: break-word" href="nz_linz_wellht1953-nzvd2016.tif">nz_linz_wellht1953-nzvd2016.tif</a> - New Zealand - Wellington - NZVD2016 height (EPSG:7839) to Wellington 1953 height (EPSG:5770). Last modified: 2020-02-02</li>
 </ul><hr><h3><a href="https://www.gov.pl/web/gugik/">Główny Urząd Geodezji i Kartografii (GUGiK)</a></h3><ul>
-<li><a style="word-break: break-word" href="pl_gugik_README.txt">pl_gugik_README.txt</a>. Last modified: 2021-11-25</li>
+<li><a style="word-break: break-word" href="pl_gugik_README.txt">pl_gugik_README.txt</a>. Last modified: 2023-07-10</li>
 <li><a style="word-break: break-word" href="pl_gugik_geoid2011-PL-EVRF2007-NH.tif">pl_gugik_geoid2011-PL-EVRF2007-NH.tif</a> - Poland - onshore - ETRF2000-PL (EPSG:9701) to EVRF2007-PL height (EPSG:9651). Last modified: 2021-11-20</li>
 <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_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>
@@ -931,7 +932,7 @@ <h2>Content</h2>
 <li><a style="word-break: break-word" href="za_cdngi_sageoid2010.tif">za_cdngi_sageoid2010.tif</a> - South Africa - ITRF2005 (EPSG:7910) to SA LLD height (EPSG:9279). Last modified: 2021-10-18</li>
 </ul>
 <p>
-Total size of content: 694 MB
+Total size of content: 695 MB
 </p>
 
 <h2>Logs</h2>

pl_gugik/build.sh → pl_gugik/build_geoid2011.sh

@@ -1,8 +1,8 @@
 #!/bin/bash
 
 # Usage:
-# PROJ_DATA_DIR=/path/to/PROJ-data ./build.sh gugik-geoid2011-PL-EVRF2007-NH.txt EPSG:9651 "EVRF2007-PL height" pl_gugik_geoid2011-PL-EVRF2007-NH.tif
-# PROJ_DATA_DIR=/path/to/PROJ-data ./build.sh gugik-geoid2011-PL-KRON86-NH.txt EPSG:9650 "Baltic 1986 height" pl_gugik_geoid2011-PL-KRON86-NH.tif
+# PROJ_DATA_DIR=/path/to/PROJ-data ./build_geoid2011.sh gugik-geoid2011-PL-EVRF2007-NH.txt EPSG:9651 "EVRF2007-PL height" pl_gugik_geoid2011-PL-EVRF2007-NH.tif
+# PROJ_DATA_DIR=/path/to/PROJ-data ./build_geoid2011.sh gugik-geoid2011-PL-KRON86-NH.txt EPSG:9650 "Baltic 1986 height" pl_gugik_geoid2011-PL-KRON86-NH.tif
 
 # Setup build directory
 mkdir build

pl_gugik/build_geoid2021.sh

@@ -0,0 +1,41 @@
+#!/bin/bash -e
+
+# Usage:
+# PROJ_DATA_DIR=/path/to/PROJ-data ./build_geoid2021.sh
+
+input=Model_quasi-geoidy-PL-geoid2021-PL-EVRF2007-NH.txt
+output=pl_gugik_geoid2021-PL-EVRF2007-NH.tif
+
+# Setup build directory
+mkdir build
+
+# Copy input file into build directory
+cp $input ./build/$input
+
+# Replace first line with a header to swap lat and long
+sed -i '1cy x z' ./build/$input
+
+docker run --user $(id -u):$(id -g) --workdir $PWD \
+            --rm -v /home:/home ghcr.io/osgeo/gdal:alpine-normal-latest \
+            sh -c " \
+            # Convert to GeoTIFF
+            gdal_translate -a_nodata 0 ./build/$input ./build/temp.tif && \
+            # Set nodata to -32768 (this also makes sure that the output is "north-up" instead of "south-up")
+            gdalwarp -srcnodata 0 -dstnodata -32768 ./build/temp.tif ./build/temp_nodata_corrected.tif && \
+            # 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:9701\" \
+            --target-crs \"EPSG:9651\" \
+            --description \"ETRF2000-PL (EPSG:9701) to EVRF2007-PL height (EPSG:9651). Converted from $input\" \
+            --area-of-use \"Poland - onshore\" \
+            --copyright \"Derived from work by GUGiK. CC-BY-4.0 https://creativecommons.org/licenses/by/4.0/\" \
+            ./build/temp_nodata_corrected.tif ./$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

pl_gugik/pl_gugik_README.txt

@@ -19,3 +19,9 @@ to heights above the ellipsoid in ETRF2000-PL (EPSG:9701).
 
 * pl_gugik_geoid2011-PL-EVRF2007-NH.tif
 * pl_gugik_geoid2011-PL-KRON86-NH.tif
+
+Vertical transformation for Geoid model PL-geoid-2021. Used to make the transitions
+from heights in vertical CRS (EPSG:9651 - EVRF2007-PL height) to heights above the ellipsoid
+in ETRF2000-PL (EPSG:9701).
+
+* pl_gugik_geoid2021-PL-EVRF2007-NH.tif

travis/expected_main.lst

@@ -254,6 +254,7 @@ nz_linz_taraht1970-nzvd2016.tif
 nz_linz_wellht1953-nzvd2016.tif
 pl_gugik_geoid2011-PL-EVRF2007-NH.tif
 pl_gugik_geoid2011-PL-KRON86-NH.tif
+pl_gugik_geoid2021-PL-EVRF2007-NH.tif
 pl_gugik_README.txt
 pt_dgt_D73_ETRS89_geo.tif
 pt_dgt_DLx_ETRS89_geo.tif