Diskussion:Höhenmodell: Unterschied zwischen den Versionen
Aus Geoinformation HSR
Stefan (Diskussion | Beiträge) (Die Seite wurde neu angelegt: „DHM-Geowebserivces: * Den OSM-Objekten z.B. fehlen ja fast durchwegs die Höhen. Um diese hinzuzufügen gibt es Webservices, z.B. von Geonames http://www.geoname…“) |
Stefan (Diskussion | Beiträge) K |
||
| Zeile 13: | Zeile 13: | ||
# die GeoTIFF-Tiles/Files wohl zu einem File astergdem.tif zusammengesetzt werden: Vgl. http://www.gis.hsr.ch/wiki/GDAL_-_Tipps_und_Tricks | # die GeoTIFF-Tiles/Files wohl zu einem File astergdem.tif zusammengesetzt werden: Vgl. http://www.gis.hsr.ch/wiki/GDAL_-_Tipps_und_Tricks | ||
# Die TIFF-Datei astergdem in PostGIS 2.0 in eine Datenbank astergdem importiert werden mittels raster2pgsql.py (dauert). Raster-Parameter 50x50 Pixels (oder 32). Vgl. http://trac.osgeo.org/postgis/wiki/WKTRaster | # Die TIFF-Datei astergdem in PostGIS 2.0 in eine Datenbank astergdem importiert werden mittels raster2pgsql.py (dauert). Raster-Parameter 50x50 Pixels (oder 32). Vgl. http://trac.osgeo.org/postgis/wiki/WKTRaster | ||
| − | # Die Datenbank astergdem kann dann wie folgt abgefragt werden, gegeben eine Koordinate long, lat: <tt><nowiki>% SELECT ST_Value(rast, ST_SetSRID(ST_MakePoint(long, lat), 4326)) FROM astergdem;</nowiki><tt> | + | # Die Datenbank astergdem kann dann wie folgt abgefragt werden, gegeben eine Koordinate long, lat: <tt><nowiki>% SELECT ST_Value(rast, ST_SetSRID(ST_MakePoint(long, lat), 4326)) FROM astergdem;</nowiki></tt> |
ASTER: | ASTER: | ||
Aktuelle Version vom 18. Dezember 2011, 22:10 Uhr
DHM-Geowebserivces:
- Den OSM-Objekten z.B. fehlen ja fast durchwegs die Höhen. Um diese hinzuzufügen gibt es Webservices, z.B. von Geonames http://www.geonames.org/export/web-services.html#astergdem .
- Das wird aber langsam und ist limitiert (man müsste zu Geonames Premium wechseln).
Hillshading:
Variante Java:
- http://www.coderanch.com/t/551724/Streams/java/Reading-tif-files-system
- Abfrage gegeben eine Koordinate long, lat: mit einem Random Access read?
Variante DB:
- die GeoTIFF-Tiles/Files wohl zu einem File astergdem.tif zusammengesetzt werden: Vgl. http://www.gis.hsr.ch/wiki/GDAL_-_Tipps_und_Tricks
- Die TIFF-Datei astergdem in PostGIS 2.0 in eine Datenbank astergdem importiert werden mittels raster2pgsql.py (dauert). Raster-Parameter 50x50 Pixels (oder 32). Vgl. http://trac.osgeo.org/postgis/wiki/WKTRaster
- Die Datenbank astergdem kann dann wie folgt abgefragt werden, gegeben eine Koordinate long, lat: % SELECT ST_Value(rast, ST_SetSRID(ST_MakePoint(long, lat), 4326)) FROM astergdem;
ASTER:
- ASTER-Daten im PostGIS Raster (Datentyp Raster). Kritische Betrachtung [1]
- Vorgehen:
- Höhenmodell-Daten „ASTER GDEM Version 2 (October 17, 2011)“ herunterladen: Vgl. http://www.gdem.aster.ersdac.or.jp/ (geht nicht?), bzw. https://wist.echo.nasa.gov/wist-bin/api/ims.cgi?mode=MAINSRCH&JS=1 :
- links „Search“ klicken,
- dann „By coordinates...“ folgendes eingeben,
- 5.959163 45.816873 , 10.560794 47.768146 (das sind die CH-Eck-Koordinaten der Schweiz in WGS84-Koordinaten/EPSG:4326)
SRTM3 v2:
- Download direkt mit Python-Skript und Import. Legt die Punkte als Koordinaten einzeln in eine PG-Tabelle mit dem normalen PostGIS-Datentyp „POINT“ ab:
- https://masterbranch.com/srtm2postgis-project/395291
- https://github.com/Sjors/srtm2postgis
- http://wiki.openstreetmap.org/wiki/Setup_Route_altitude_profiles_SRTM
