Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add 4 height correction / geoid models for Austria #13

Merged
merged 2 commits into from
Feb 26, 2020
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
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