Merge pull request #109 from jjimenezshaw/german-geoid

Add GCG2016 German geoid model file

CMakeLists.txt

@@ -27,6 +27,7 @@ set(CPACK_SOURCE_INSTALLED_DIRECTORIES
     "${CMAKE_CURRENT_SOURCE_DIR}/ca_que_mern/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/ch_swisstopo/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/de_adv/" "."
+    "${CMAKE_CURRENT_SOURCE_DIR}/de_bkg/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/de_geosn/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/de_hvbg/" "."
     "${CMAKE_CURRENT_SOURCE_DIR}/de_lgl_bw/" "."

agency.json

@@ -59,6 +59,12 @@
         "agency":"Arbeitsgemeinschaft der Vermessungsverwaltungender der Länder der Bundesrepublik Deutschland (AdV)",
         "url": "http://www.adv-online.de/icc/extdeu/nav/9ae/9ae594bb-a094-311a-3b21-718a438ad1b2&sel_uCon=3c10e056-7955-311a-3b21-718a438ad1b2&uTem=73d607d6-b048-65f1-80fa-29f08a07b51a.htm"
     },
+    {
+        "id": "de_bkg",
+        "country": "Germany",
+        "agency":"Bundesamt für Kartographie und Geodäsie (BKG)",
+        "url": "https://gdz.bkg.bund.de/index.php/default/digitale-geodaten.html"
+    },
     {
         "id": "de_geosn",
         "country": "Germany",

copyright_and_licenses.csv

@@ -75,6 +75,8 @@ ch_swisstopo_chgeo2004_ETRS89_LN02.tif,Swisstopo,CC0-1.0,,
 ch_swisstopo_README.txt,Disclaimed,Public domain,,
 de_adv_BETA2007.tif,ADV,Free redistribution is allowed and welcome,,
 de_adv_README.txt,Disclaimed,Public domain,,
+de_bkg_gcg2016.tif,BKG,CC-BY-4.0,1.16,
+de_bkg_README.txt,Disclaimed,Public domain,1.16,
 de_geosn_NTv2_SN.tif,Staatsbetrieb Geobasisinformation und Vermessung Sachsen (GEOSN),Data licence Germany - attribution - version 2.0,,
 de_geosn_README.txt,Disclaimed,Public domain,,
 de_hvbg_hessen_HeTA2010.tif,Hessisches Landesamt für Geoinformation und Bodenmanagement (HVBG),Data licence Germany - attribution - version 2.0,1.14,

de_bkg/.github/README.md

@@ -0,0 +1 @@
+../de_bkg_README.txt

de_bkg/build_GCG2016.sh

@@ -0,0 +1,40 @@
+#!/bin/bash -e
+
+# Usage:
+# PROJ_DATA_DIR=/path/to/PROJ-data ./build_GCG2016.sh [path/to/GCG2016.txt]
+
+# Setup build directory
+mkdir -p build
+
+src="${1:-GCG2016.txt}"
+dst="de_bkg_gcg2016.tif"
+# Copy input file into build directory
+cp $src ./build/GCG2016.txt
+
+# Add header
+sed -i '1i\y x z' ./build/GCG2016.txt
+
+docker run  --pull=always --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 999999 ./build/GCG2016.txt ./build/temp.tif && \
+            # Set nodata to -32768
+            gdalwarp -srcnodata 999999 -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:10283\" \
+            --target-crs \"EPSG:7837\" \
+            --description \"ETRS89/DREF91/2016 (EPSG:10283) to DHHN2016 height (EPSG:7837). Converted from GCG2016.txt\" \
+            --area-of-use \"Germany\" \
+            --copyright \"Derived from work by BKG. CC-BY-4.0 https://creativecommons.org/licenses/by/4.0/\" \
+            ./build/temp_nodata_corrected.tif ./$dst && \
+            # Run check_gtiff_grid-script
+            python3 ${PROJ_DATA_DIR}/grid_tools/check_gtiff_grid.py ./$dst \
+            # Show info
+            gdalinfo ./$dst "
+
+# Remove build directory
+rm -rf build

de_bkg/de_bkg_README.txt

@@ -0,0 +1,21 @@
+# de_bkg_README.txt
+
+The files in this section result from the conversion of datasets originating
+from [BKG](https://www.bkg.bund.de/)
+
+## Included grids
+
+### German vertical grid:
+
+*Source*: [BKG](https://gdz.bkg.bund.de/index.php/default/digitale-geodaten/geodaetische-basisdaten/quasigeoid-der-bundesrepublik-deutschland-quasigeoid.html)
+*Format*: GeoTIFF converted from 'XYZ/ASCII' (adjusted)
+*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)
+*Credit*: (c) Bundesamt für Kartographie und Geodäsie - BKG - Deutschland
+*Horizontal CRS*: EPSG:10283 (ETRS89/DREF91/2016)
+
+Vertical transformation for Geoid model GCG2016. Used to make the transition
+from heights in vertical CRS DHHN2016 height (EPSG:7837)
+to heights above the ellipsoid in ETRS89/DREF91/2016 (EPSG:10283).
+
+* de_bkg_gcg2016.tif
+

files.geojson

@@ -77,6 +77,8 @@
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/ch_swisstopo_chgeo2004_ETRS89_LN02.tif", "name": "ch_swisstopo_chgeo2004_ETRS89_LN02.tif", "area_of_use": "Switzerland", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4937", "source_crs_name": "ETRS89", "target_crs_code": "EPSG:5728", "target_crs_name": "LN02 height", "source": "Swisstopo Federal Office of Topography", "source_country": "Switzerland", "source_id": "ch_swisstopo", "source_url": "https://www.swisstopo.admin.ch/en/knowledge-facts/surveying-geodesy/reference-frames/transformations-position.html", "description": "ETRS89 (EPSG:4937) to LN02 height (EPSG:5728)", "file_size": 295192, "sha256sum": "ec9638c152f4519278b3e7a781d1a7c2769bd54de39accec099bd5036f4fd214" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 5.85, 45.75 ], [ 10.5, 45.75 ], [ 10.5, 47.85 ], [ 5.85, 47.85 ], [ 5.85, 45.75 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_adv_README.txt", "name": "de_adv_README.txt", "source": "Arbeitsgemeinschaft der Vermessungsverwaltungender der Länder der Bundesrepublik Deutschland (AdV)", "source_country": "Germany", "source_id": "de_adv", "source_url": "http://www.adv-online.de/icc/extdeu/nav/9ae/9ae594bb-a094-311a-3b21-718a438ad1b2&sel_uCon=3c10e056-7955-311a-3b21-718a438ad1b2&uTem=73d607d6-b048-65f1-80fa-29f08a07b51a.htm", "file_size": 802, "sha256sum": "69c387720cba6f596671de999d8a9683194fa5787294d2e10f81ed90e04ba11c" }, "geometry": null },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_adv_BETA2007.tif", "name": "de_adv_BETA2007.tif", "area_of_use": "Germany", "type": "HORIZONTAL_OFFSET", "source_crs_code": "EPSG:4314", "source_crs_name": "DHDN", "target_crs_code": "EPSG:4258", "target_crs_name": "ETRS89", "source": "Arbeitsgemeinschaft der Vermessungsverwaltungender der Länder der Bundesrepublik Deutschland (AdV)", "source_country": "Germany", "source_id": "de_adv", "source_url": "http://www.adv-online.de/icc/extdeu/nav/9ae/9ae594bb-a094-311a-3b21-718a438ad1b2&sel_uCon=3c10e056-7955-311a-3b21-718a438ad1b2&uTem=73d607d6-b048-65f1-80fa-29f08a07b51a.htm", "description": "DHDN (EPSG:4314) to ETRS89 (EPSG:4258)", "file_size": 24379, "sha256sum": "46e681fcc7d022dde1db1f9d0a3426a9bfb1d4a151af69a81b3c30104c9388e2" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 5.5, 47.0 ], [ 15.6666667, 47.0 ], [ 15.6666667, 55.3 ], [ 5.5, 55.3 ], [ 5.5, 47.0 ] ] ] } },
+{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_bkg_README.txt", "name": "de_bkg_README.txt", "source": "Bundesamt für Kartographie und Geodäsie (BKG)", "source_country": "Germany", "source_id": "de_bkg", "source_url": "https://gdz.bkg.bund.de/index.php/default/digitale-geodaten.html", "file_size": 840, "sha256sum": "032f7f735dd47f381488cc9476b4ff491caea4e094fbe8d4688da3e88515be51", "version_added": "1.16" }, "geometry": null },
+{ "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_bkg_gcg2016.tif", "name": "de_bkg_gcg2016.tif", "area_of_use": "Germany", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:10283", "source_crs_name": "ETRS89/DREF91/2016", "target_crs_code": "EPSG:7837", "target_crs_name": "DHHN2016 height", "source": "Bundesamt für Kartographie und Geodäsie (BKG)", "source_country": "Germany", "source_id": "de_bkg", "source_url": "https://gdz.bkg.bund.de/index.php/default/digitale-geodaten.html", "description": "ETRS89/DREF91/2016 (EPSG:10283) to DHHN2016 height (EPSG:7837)", "file_size": 995647, "sha256sum": "598f18324dea7f8e72421d18add7ac6228259adf91eeb335cc9c27d98484f7ac", "version_added": "1.16" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 3.25625, 47.2208333 ], [ 15.11875, 47.2208333 ], [ 15.11875, 55.9791667 ], [ 3.25625, 55.9791667 ], [ 3.25625, 47.2208333 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_geosn_README.txt", "name": "de_geosn_README.txt", "source": "Staatsbetrieb Geobasisinformation und Vermessung Sachsen GeoSN", "source_country": "Germany", "source_id": "de_geosn", "source_url": "https://www.landesvermessung.sachsen.de/download-5608.html?_cp=%7B%22accordion-content-4389%22%3A%7B%222%22%3Atrue%7D%2C%22previousOpen%22%3A%7B%22group%22%3A%22accordion-content-4389%22%2C%22idx%22%3A2%7D%7D", "file_size": 610, "sha256sum": "0c3da042e030577eba00a41575706d3bb40449b7a92dc62e86e104e94772f51e" }, "geometry": null },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_geosn_NTv2_SN.tif", "name": "de_geosn_NTv2_SN.tif", "area_of_use": "Germany - Saxony", "type": "HORIZONTAL_OFFSET", "source_crs_code": "EPSG:4745", "source_crs_name": "RD/83", "target_crs_code": "EPSG:4258", "target_crs_name": "ETRS89", "source": "Staatsbetrieb Geobasisinformation und Vermessung Sachsen GeoSN", "source_country": "Germany", "source_id": "de_geosn", "source_url": "https://www.landesvermessung.sachsen.de/download-5608.html?_cp=%7B%22accordion-content-4389%22%3A%7B%222%22%3Atrue%7D%2C%22previousOpen%22%3A%7B%22group%22%3A%22accordion-content-4389%22%2C%22idx%22%3A2%7D%7D", "description": "RD/83 (EPSG:4745) to ETRS89 (EPSG:4258)", "file_size": 4250734, "sha256sum": "211caa07476e351146117beb9281184b97caeae58ff9724517b4630ac03f7028" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 15.6666667, 49.8 ], [ 15.6666667, 54.7 ], [ 9.8333333, 54.7 ], [ 9.8333333, 49.8 ], [ 15.6666667, 49.8 ] ] ] } },
 { "type": "Feature", "properties": { "url": "https://cdn.proj.org/de_hvbg_README.txt", "name": "de_hvbg_README.txt", "source": "Hessisches Landesamt für Bodenmanagement und Geoinformation", "source_country": "Germany", "source_id": "de_hvbg", "source_url": "https://www.hvbg.hessen.de/", "file_size": 924, "sha256sum": "b7815d9872c089f82027120d63202faf159d55582a92268c0b742fb83c3e4e95" }, "geometry": null },

index.html

@@ -549,6 +549,9 @@ <h2>Content</h2>
 </ul><hr><h3><a href="http://www.adv-online.de/icc/extdeu/nav/9ae/9ae594bb-a094-311a-3b21-718a438ad1b2&amp;sel_uCon=3c10e056-7955-311a-3b21-718a438ad1b2&amp;uTem=73d607d6-b048-65f1-80fa-29f08a07b51a.htm">Arbeitsgemeinschaft der Vermessungsverwaltungender der Länder der Bundesrepublik Deutschland (AdV)</a></h3><ul>
 <li><a style="word-break: break-word" href="de_adv_README.txt">de_adv_README.txt</a>. Last modified: 2020-01-24</li>
 <li><a style="word-break: break-word" href="de_adv_BETA2007.tif">de_adv_BETA2007.tif</a> - Germany - DHDN (EPSG:4314) to ETRS89 (EPSG:4258). Last modified: 2020-01-24</li>
+</ul><hr><h3><a href="https://gdz.bkg.bund.de/index.php/default/digitale-geodaten.html">Bundesamt für Kartographie und Geodäsie (BKG)</a></h3><ul>
+<li><a style="word-break: break-word" href="de_bkg_README.txt">de_bkg_README.txt</a>. Last modified: 2023-10-14</li>
+<li><a style="word-break: break-word" href="de_bkg_gcg2016.tif">de_bkg_gcg2016.tif</a> - Germany - ETRS89/DREF91/2016 (EPSG:10283) to DHHN2016 height (EPSG:7837). Last modified: 2023-10-14</li>
 </ul><hr><h3><a href="https://www.landesvermessung.sachsen.de/download-5608.html?_cp=%7B%22accordion-content-4389%22%3A%7B%222%22%3Atrue%7D%2C%22previousOpen%22%3A%7B%22group%22%3A%22accordion-content-4389%22%2C%22idx%22%3A2%7D%7D">Staatsbetrieb Geobasisinformation und Vermessung Sachsen GeoSN</a></h3><ul>
 <li><a style="word-break: break-word" href="de_geosn_README.txt">de_geosn_README.txt</a>. Last modified: 2020-01-24</li>
 <li><a style="word-break: break-word" href="de_geosn_NTv2_SN.tif">de_geosn_NTv2_SN.tif</a> - Germany - Saxony - RD/83 (EPSG:4745) to ETRS89 (EPSG:4258). Size: 4.1 MB. Last modified: 2020-01-24</li>
@@ -936,7 +939,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: 700 MB
+Total size of content: 701 MB
 </p>
 
 <h2>Logs</h2>

travis/expected_main.lst

@@ -74,6 +74,8 @@ ch_swisstopo_README.txt
 copyright_and_licenses.csv
 de_adv_BETA2007.tif
 de_adv_README.txt
+de_bkg_gcg2016.tif
+de_bkg_README.txt
 de_geosn_NTv2_SN.tif
 de_geosn_README.txt
 de_hvbg_hessen_HeTA2010.tif