Mehr

Programmatische DEM-Downloads ähnlich wie landsat-util

Programmatische DEM-Downloads ähnlich wie landsat-util


Ich habe vor kurzem landsat-util entdeckt und mich gefragt, ob es eine gleichwertige Befehlszeilenanwendung für globale Höhendaten (für Daten mit einer Auflösung von etwa 1 Bogensekunde, vorzugsweise Geotiff) gibt oder ob es eine besonders saubere programmatische Möglichkeit gibt, diese Art von Daten herunterzuladen . Bisher habe ich nur Massen-Downloads gefunden (wobei ich nur relevante Kacheln herunterladen möchte, da ich sie mitten im Skript herunterladen werde) und die NED (die mich korrigieren, wenn ich falsch liege, aber ich glaube, sie deckt nur die USA ab? )

Ich hoffe, dies ist den verlinkten Abfragen nicht zu ähnlich.


Es gibt ein Tool aus den OTB-Anwendungen, um SRTM-Kacheln herunterzuladen, die sich auf eine Reihe von Bildern beziehen. Grundsätzlich liest es Bildinformationen und lädt dann die zugehörigen Kacheln von http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/ herunter. Dies ist eine Open-Source-Bibliothek in C++, sodass Sie sich den Code im Detail ansehen können, wenn Sie wissen möchten, wie es geht.

BEARBEITEN basierend auf @ Kerstens Kommentar: SRTM war früher überall bei 90 verfügbar, außer in den USA, wo es 30 m waren. In letzter Zeit wurden immer mehr SRTM-Kacheln mit einer Auflösung von ~30 Metern veröffentlicht, um die meisten Teile der Welt abzudecken, mit der größten Ausnahme im Nahen Osten und in Nordafrika.

ASTER DEM ist eine genauere Alternative, leidet jedoch an einigen Stellen unter großen Artefakten, also verwenden Sie es vorsichtig, wenn Sie dies tun.


Eine ähnliche Frage habe ich hier beantwortet

Sie können einen Teil des SRTM 30m DEM mit einem Befehl mit dem Python-Befehlszeilentool für die Höhe herunterladen und ausschneiden.

Installieren Sie es und führen Sie den Selbsttest durch mit:

$ pip Installationshöhe

Überprüfen Sie, ob alle Abhängigkeiten installiert sind (hauptsächlich GDAL-Tools):

$ eio-Selbsttest

Laden Sie einen Teil des SRTM 30m DEM herunter und schneiden Sie ihn aus, der die Begrenzungsbox als WGS84-Koordinaten angibt:

$ eio clip -o Rome-DEM.tif --bounds 12,35 41,8 12,65 42