Skip to content

Commit

Permalink
Merge pull request #13 from rouault/at_bev_height_grids
Browse files Browse the repository at this point in the history
Add 4 height correction / geoid models for Austria
  • Loading branch information
rouault authored Feb 26, 2020
2 parents b256af6 + 13a7624 commit af78e35
Show file tree
Hide file tree
Showing 11 changed files with 97 additions and 3 deletions.
Binary file added at_bev/at_bev_GEOID_BESSEL_Oesterreich.tif
Binary file not shown.
Binary file added at_bev/at_bev_GEOID_GRS80_Oesterreich.tif
Binary file not shown.
Binary file added at_bev/at_bev_GV_Hoehengrid_V1.tif
Binary file not shown.
Binary file added at_bev/at_bev_GV_Hoehengrid_plus_Geoid_V2.tif
Binary file not shown.
45 changes: 45 additions & 0 deletions at_bev/at_bev_README.txt
Original file line number Diff line number Diff line change
Expand Up @@ -12,3 +12,48 @@ from [Austria Bundesamt für Eich- und Vermessungswessen](www.bev.gv.at)
*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)

* at_bev_AT_GIS_GRID.tif

### Austria :GHA height -> EVRF2000 Austria height

*Source*: [BEV](http://www.bev.gv.at/portal/page?_pageid=713,2601281&_dad=portal&_schema=PORTAL)
*Format*: GeoTIFF converted from original GeoTIFF
*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)
*Source CRS*: EPSG:5778 (GHA height)
*Target CRS*: EPSG:9274 (EVRF2000 Austria height)
*Interpolation CRS*: EPSG:4312 (MGI)
*Used by*: [EPSG:9275 Transformation](https://www.epsg-registry.org/export.htm?gml=urn:ogc:def:coordinateOperation:EPSG::9275)

* at_bev_GV_Hoehengrid_V1.tif

### Austria :ETRS89 -> EVRF2000 Austria height

*Source*: [BEV](http://www.bev.gv.at/portal/page?_pageid=713,2601285&_dad=portal&_schema=PORTAL)
*Format*: GeoTIFF converted from original CSV file GEOID_GRS80_Oesterreich.csv from http://www.bev.gv.at/pls/portal/docs/PAGE/BEV_PORTAL_CONTENT_ALLGEMEIN/0200_PRODUKTE/UNENTGELTLICHE_PRODUKTE_DES_BEV/GV_GEOID_Oesterreich.zip
*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)
*Source CRS*: EPSG:4937 (ETRS89)
*Target CRS*: EPSG:9274 (EVRF2000 Austria height)
*Used by*: [EPSG:9276 Transformation](https://www.epsg-registry.org/export.htm?gml=urn:ogc:def:coordinateOperation:EPSG::9276)

* at_bev_GEOID_GRS80_Oesterreich.tif

### Austria :MGI -> EVRF2000 Austria height

*Source*: [BEV](http://www.bev.gv.at/portal/page?_pageid=713,2601285&_dad=portal&_schema=PORTAL)
*Format*: GeoTIFF converted from original CSV file GEOID_BESSEL_Oesterreich.csv from http://www.bev.gv.at/pls/portal/docs/PAGE/BEV_PORTAL_CONTENT_ALLGEMEIN/0200_PRODUKTE/UNENTGELTLICHE_PRODUKTE_DES_BEV/GV_GEOID_Oesterreich.zip
*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)
*Source CRS*: EPSG:9267 (MGI)
*Target CRS*: EPSG:9274 (EVRF2000 Austria height)
*Used by*: [EPSG:9277 Transformation](https://www.epsg-registry.org/export.htm?gml=urn:ogc:def:coordinateOperation:EPSG::9277)

* at_bev_GEOID_BESSEL_Oesterreich.tif

### Austria :ETRS89 -> GHA Austria height

*Source*: [BEV](http://www.bev.gv.at/portal/page?_pageid=713,2823796&_dad=portal&_schema=PORTAL)
*Format*: GeoTIFF converted from original GeoTIFF
*License*: [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/)
*Source CRS*: EPSG:4937 (ETRS89)
*Target CRS*: EPSG:5778 (GHA Austria height)
*Used by*: [EPSG:9278 Transformation](https://www.epsg-registry.org/export.htm?gml=urn:ogc:def:coordinateOperation:EPSG::9278)

* at_bev_GV_Hoehengrid_plus_Geoid_V2.tif
28 changes: 28 additions & 0 deletions at_bev/build.sh
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
#!/bin/sh

rm -rf build_tmp
mkdir build_tmp
cd build_tmp

wget http://www.bev.gv.at/pls/portal/docs/PAGE/BEV_PORTAL_CONTENT_ALLGEMEIN/0200_PRODUKTE/UNENTGELTLICHE_PRODUKTE_DES_BEV/Hoehen-Grid_GeoTIFF.zip
unzip -o Hoehen-Grid_GeoTIFF.zip
python3 ../../grid_tools/vertoffset_grid_to_gtiff.py --type VERTICAL_TO_VERTICAL --source-crs EPSG:5778 --target-crs EPSG:9274 --interpolation-crs EPSG:4312 --copyright "Derived from work by BEV. Creative Commons Attribution 4.0 https://creativecommons.org/licenses/by/4.0/" --description "GHA height (EPSG:5778) to EVRF2000 Austria height (EPSG:9274). Converted from GV_Hoehengrid_V1.tif" --area-of-use Austria GV_Hoehengrid_V1.tif ../at_bev_GV_Hoehengrid_V1.tif

wget "http://www.bev.gv.at/pls/portal/docs/PAGE/BEV_PORTAL_CONTENT_ALLGEMEIN/0200_PRODUKTE/UNENTGELTLICHE_PRODUKTE_DES_BEV/GV_Hoehengrid_plus_Geoid_V2(GeoTIFF).zip"
unzip -o "GV_Hoehengrid_plus_Geoid_V2(GeoTIFF).zip"
python3 ../../grid_tools/vertoffset_grid_to_gtiff.py --type GEOGRAPHIC_TO_VERTICAL --source-crs EPSG:4937 --target-crs EPSG:5778 --copyright "Derived from work by BEV. Creative Commons Attribution 4.0 https://creativecommons.org/licenses/by/4.0/" --description "ETRS89 (EPSG:4937) to GHA Austria height (EPSG:5778). Converted from GV_Hoehengrid_plus_Geoid_V2.tif" --area-of-use Austria GV_Hoehengrid_plus_Geoid_V2.tif ../at_bev_GV_Hoehengrid_plus_Geoid_V2.tif

wget http://www.bev.gv.at/pls/portal/docs/PAGE/BEV_PORTAL_CONTENT_ALLGEMEIN/0200_PRODUKTE/UNENTGELTLICHE_PRODUKTE_DES_BEV/GV_GEOID_Oesterreich.zip
unzip -o GV_GEOID_Oesterreich.zip

ogr2ogr GEOID_GRS80_Oesterreich.gpkg GEOID_GRS80_Oesterreich.csv -oo Y_POSSIBLE_NAMES=BREITE -oo X_POSSIBLE_NAMES=LAENGE
gdal_grid GEOID_GRS80_Oesterreich.gpkg GEOID_GRS80_Oesterreich.tif -zfield UNDULATION -ot Float32 -txe 9.4791666666666667 17.270833333333333 -tye 49.0875 46.3125 -outsize 187 111 -a nearest:radius1=0.001:radius2=0.001:nodata=-32768
python3 ../../grid_tools/vertoffset_grid_to_gtiff.py --type GEOGRAPHIC_TO_VERTICAL --source-crs EPSG:4937 --target-crs EPSG:9274 --copyright "Derived from work by BEV. Creative Commons Attribution 4.0 https://creativecommons.org/licenses/by/4.0/" --description "ETRS89 (EPSG:4937) to EVRF2000 Austria height (EPSG:9274). Converted from GEOID_GRS80_Oesterreich.csv" --area-of-use Austria GEOID_GRS80_Oesterreich.tif ../at_bev_GEOID_GRS80_Oesterreich.tif

ogr2ogr GEOID_BESSEL_Oesterreich.gpkg GEOID_BESSEL_Oesterreich.csv -oo Y_POSSIBLE_NAMES=BREITE -oo X_POSSIBLE_NAMES=LAENGE
gdal_grid GEOID_BESSEL_Oesterreich.gpkg GEOID_BESSEL_Oesterreich.tif -zfield UNDULATION -ot Float32 -txe 9.4791666666666667 17.270833333333333 -tye 49.0875 46.3125 -outsize 187 111 -a nearest:radius1=0.001:radius2=0.001:nodata=-32768
python3 ../../grid_tools/vertoffset_grid_to_gtiff.py --type GEOGRAPHIC_TO_VERTICAL --source-crs EPSG:9267 --target-crs EPSG:9274 --copyright "Derived from work by BEV. Creative Commons Attribution 4.0 https://creativecommons.org/licenses/by/4.0/" --description "MGI (EPSG:9267) to EVRF2000 Austria height (EPSG:9274). Converted from GEOID_BESSEL_Oesterreich.csv" --area-of-use Austria GEOID_BESSEL_Oesterreich.tif ../at_bev_GEOID_BESSEL_Oesterreich.tif

cd ..
rm -rf build_tmp

4 changes: 4 additions & 0 deletions copyright_and_licenses.csv
Original file line number Diff line number Diff line change
@@ -1,5 +1,9 @@
filename,copyright,license
at_bev_AT_GIS_GRID.tif,"Austria Bundesamt für Eich- und Vermessungswessen (BEV)",CC-BY-4.0
at_bev_GEOID_BESSEL_Oesterreich.tif,"Austria Bundesamt für Eich- und Vermessungswessen (BEV)",CC-BY-4.0
at_bev_GEOID_GRS80_Oesterreich.tif,"Austria Bundesamt für Eich- und Vermessungswessen (BEV)",CC-BY-4.0
at_bev_GV_Hoehengrid_plus_Geoid_V2.tif,"Austria Bundesamt für Eich- und Vermessungswessen (BEV)",CC-BY-4.0
at_bev_GV_Hoehengrid_V1.tif,"Austria Bundesamt für Eich- und Vermessungswessen (BEV)",CC-BY-4.0
at_bev_README.txt,Disclaimed,Public domain
au_ga_AUSGeoid09_V1.01.tif,Geoscience Australia,CC-BY-4.0
au_ga_AUSGeoid2020_20180201.tif,Geoscience Australia,CC-BY-4.0
Expand Down
6 changes: 5 additions & 1 deletion files.geojson
Original file line number Diff line number Diff line change
Expand Up @@ -2,8 +2,12 @@
"type": "FeatureCollection",
"name": "files",
"features": [
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_README.txt", "name": "at_bev_README.txt", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "file_size": 476 }, "geometry": null },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_README.txt", "name": "at_bev_README.txt", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "file_size": 3042 }, "geometry": null },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_AT_GIS_GRID.tif", "name": "at_bev_AT_GIS_GRID.tif", "area_of_use": "Austria", "type": "HORIZONTAL_OFFSET", "source_crs_code": "EPSG:4312", "source_crs_name": "MGI", "target_crs_code": "EPSG:4258", "target_crs_name": "ETRS89", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "description": "MGI (EPSG:4312) to ETRS89 (EPSG:4258)", "file_size": 216878 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 9.5, 46.35 ], [ 17.1625, 46.35 ], [ 17.1625, 49.05 ], [ 9.5, 49.05 ], [ 9.5, 46.35 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_GEOID_BESSEL_Oesterreich.tif", "name": "at_bev_GEOID_BESSEL_Oesterreich.tif", "area_of_use": "Austria", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:9267", "source_crs_name": "MGI", "target_crs_code": "EPSG:9274", "target_crs_name": "EVRF2000 Austria height", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "description": "MGI (EPSG:9267) to EVRF2000 Austria height (EPSG:9274)", "file_size": 38547 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 9.5, 46.325 ], [ 17.25, 46.325 ], [ 17.25, 49.075 ], [ 9.5, 49.075 ], [ 9.5, 46.325 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_GEOID_GRS80_Oesterreich.tif", "name": "at_bev_GEOID_GRS80_Oesterreich.tif", "area_of_use": "Austria", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4937", "source_crs_name": "ETRS89", "target_crs_code": "EPSG:9274", "target_crs_name": "EVRF2000 Austria height", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "description": "ETRS89 (EPSG:4937) to EVRF2000 Austria height (EPSG:9274)", "file_size": 29948 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 9.5, 46.325 ], [ 17.25, 46.325 ], [ 17.25, 49.075 ], [ 9.5, 49.075 ], [ 9.5, 46.325 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_GV_Hoehengrid_V1.tif", "name": "at_bev_GV_Hoehengrid_V1.tif", "area_of_use": "Austria", "type": "VERTICAL_OFFSET_VERTICAL_TO_VERTICAL", "source_crs_code": "EPSG:5778", "source_crs_name": "GHA height", "target_crs_code": "EPSG:9274", "target_crs_name": "EVRF2000 Austria height", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "description": "GHA height (EPSG:5778) to EVRF2000 Austria height (EPSG:9274)", "file_size": 1108072 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 9.5, 46.35 ], [ 17.1875, 46.35 ], [ 17.1875, 49.05 ], [ 9.5, 49.05 ], [ 9.5, 46.35 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/at_bev_GV_Hoehengrid_plus_Geoid_V2.tif", "name": "at_bev_GV_Hoehengrid_plus_Geoid_V2.tif", "area_of_use": "Austria", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4937", "source_crs_name": "ETRS89", "target_crs_code": "EPSG:5778", "target_crs_name": "GHA height", "source": "Austria Bundesamt für Eich- und Vermessungswessen", "source_country": "Austria", "source_id": "at_bev", "source_url": "http:\/\/www.bev.gv.at\/portal\/page?_pageid=713,2157075&_dad=portal&_schema=PORTAL", "description": "ETRS89 (EPSG:4937) to GHA Austria height (EPSG:5778)", "file_size": 770542 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 9.5, 46.35 ], [ 17.1875, 46.35 ], [ 17.1875, 49.05 ], [ 9.5, 49.05 ], [ 9.5, 46.35 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/au_ga_README.txt", "name": "au_ga_README.txt", "source": "Geoscience Australia", "source_country": "Australia", "source_id": "au_ga", "source_url": "http:\/\/www.ga.gov.au", "file_size": 1468 }, "geometry": null },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/au_ga_AUSGeoid09_V1.01.tif", "name": "au_ga_AUSGeoid09_V1.01.tif", "area_of_use": "Australia", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:4939", "source_crs_name": "GDA94", "target_crs_code": "EPSG:5711", "target_crs_name": "AHD height", "source": "Geoscience Australia", "source_country": "Australia", "source_id": "au_ga", "source_url": "http:\/\/www.ga.gov.au", "description": "GDA94 (EPSG:4939) to AHD height (EPSG:5711)", "file_size": 12781759 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 108.0, -45.9833333 ], [ 159.9833333, -45.9833333 ], [ 159.9833333, -8.0 ], [ 108.0, -8.0 ], [ 108.0, -45.9833333 ] ] ] } },
{ "type": "Feature", "properties": { "url": "https:\/\/cdn.proj.org\/au_ga_AUSGeoid2020_20180201.tif", "name": "au_ga_AUSGeoid2020_20180201.tif", "area_of_use": "Australia", "type": "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL", "source_crs_code": "EPSG:7843", "source_crs_name": "GDA2020", "target_crs_code": "EPSG:5711", "target_crs_name": "AHD height", "source": "Geoscience Australia", "source_country": "Australia", "source_id": "au_ga", "source_url": "http:\/\/www.ga.gov.au", "description": "GDA2020 (EPSG:7843) to AHD height (EPSG:5711)", "file_size": 5839686 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 93.0, -60.9833333 ], [ 173.9833333, -60.9833333 ], [ 173.9833333, -8.0 ], [ 93.0, -8.0 ], [ 93.0, -60.9833333 ] ] ] } },
Expand Down
8 changes: 6 additions & 2 deletions index.html
Original file line number Diff line number Diff line change
Expand Up @@ -432,8 +432,12 @@ <h2>Content</h2>
<ul>
<li><a href="README.DATA">README.DATA</a></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 href="at_bev_README.txt">at_bev_README.txt</a>. Last modified: 2020-01-24</li>
<li><a href="at_bev_README.txt">at_bev_README.txt</a>. Last modified: 2020-02-26</li>
<li><a 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>
<li><a href="at_bev_GEOID_BESSEL_Oesterreich.tif">at_bev_GEOID_BESSEL_Oesterreich.tif</a> - Austria - MGI (EPSG:9267) to EVRF2000 Austria height (EPSG:9274). Last modified: 2020-02-26</li>
<li><a href="at_bev_GEOID_GRS80_Oesterreich.tif">at_bev_GEOID_GRS80_Oesterreich.tif</a> - Austria - ETRS89 (EPSG:4937) to EVRF2000 Austria height (EPSG:9274). Last modified: 2020-02-26</li>
<li><a href="at_bev_GV_Hoehengrid_V1.tif">at_bev_GV_Hoehengrid_V1.tif</a> - Austria - GHA height (EPSG:5778) to EVRF2000 Austria height (EPSG:9274). Size: 1.1 MB. Last modified: 2020-02-26</li>
<li><a href="at_bev_GV_Hoehengrid_plus_Geoid_V2.tif">at_bev_GV_Hoehengrid_plus_Geoid_V2.tif</a> - Austria - ETRS89 (EPSG:4937) to GHA Austria height (EPSG:5778). Last modified: 2020-02-26</li>
</ul><hr><h3><a href="http://www.ga.gov.au">Geoscience Australia</a></h3><ul>
<li><a href="au_ga_README.txt">au_ga_README.txt</a>. Last modified: 2020-01-24</li>
<li><a href="au_ga_AUSGeoid09_V1.01.tif">au_ga_AUSGeoid09_V1.01.tif</a> - Australia - GDA94 (EPSG:4939) to AHD height (EPSG:5711). Size: 12.2 MB. Last modified: 2020-01-24</li>
Expand Down Expand Up @@ -727,7 +731,7 @@ <h2>Content</h2>
<li><a href="us_noaa_wyhpgn.tif">us_noaa_wyhpgn.tif</a> - USA - Wyoming - NAD83 (EPSG:4269) to NAD83(HARN) (EPSG:4152). Last modified: 2020-01-24</li>
</ul>
<p>
Total size of content: 487 MB
Total size of content: 489 MB
</p>

<h2>Logs</h2>
Expand Down
5 changes: 5 additions & 0 deletions regenerate_index_html.py
100644 → 100755
Original file line number Diff line number Diff line change
@@ -1,3 +1,5 @@
#!/usr/bin/env python3

from osgeo import gdal, ogr, osr
import glob
import os
Expand Down Expand Up @@ -58,6 +60,9 @@
readme_filename = None
for f in glob.glob(dirname + '/*'):
f = os.path.basename(f)
if f.endswith('.sh'):
continue

if f.endswith('README.txt'):
assert not readme_filename
readme_filename = f
Expand Down
4 changes: 4 additions & 0 deletions travis/expected_main.lst
Original file line number Diff line number Diff line change
@@ -1,4 +1,8 @@
at_bev_AT_GIS_GRID.tif
at_bev_GEOID_BESSEL_Oesterreich.tif
at_bev_GEOID_GRS80_Oesterreich.tif
at_bev_GV_Hoehengrid_plus_Geoid_V2.tif
at_bev_GV_Hoehengrid_V1.tif
at_bev_README.txt
au_ga_AUSGeoid09_V1.01.tif
au_ga_AUSGeoid2020_20180201.tif
Expand Down

0 comments on commit af78e35

Please sign in to comment.