Ich habe ein paar Scripts geschrieben, die die oben genannten Befehle für den Kanton Luzern ausführt. Die Scripts findet man hier [1], falls es jemand interessiert und darauf aufbauen möchte.
Wer das Resultat ausprobieren möchte, kann folgendes URL Template verwenden: https://marzili.duckdns.org/tmsluzern/%7Bzoom%7D/%7Bx%7D/%7B-y%7D.jpg Momentan sind nur Zoom Level 16-19 berechnet. D.h. man muss schon recht Zoomen bis etwas angezeigt wird. Die Tiles sind momentan auf einem Rasperry Pi, und meine Internetleitung ist nicht sehr schnell. D.h. rechnet nicht damit, dass sie erreichbar sind.
Die Qualität ist nicht vergleichbar mit dem SwissAlti3D Hillshade. Aber bis dieses für den Kanton Luzern verfügbar ist, wird man es durchaus nutzen können, um z.B. Strassen/Wege/Flüsse im Wald zu mappen.
lg rene
[1] https://github.com/rbuffat/swisssurface_dtm
On Thu, 13 Jan 2022 at 22:23, René Buffat buffat@gmail.com wrote:
Eine Möglichkeit wäre, aus der Punktwolke selber ein DTM herzustellen. Ich habe das noch nie gemacht, aber nach ein bisschen Googlen scheint es durchaus machbar. Nur benötigt man wohl recht viel CPU / Speicherplatz und noch einiges an Parameter Tuning. Von dem her ist es vielleicht besser zu warten, bis SwissTopo es macht.
Folgender aus Tutorials zusammengesetzter Workflow produziert dieses Hillshade: https://i.imgur.com/vhF4KPb.png
- Filtern des Bodens:
https://pdal.io/workshop/exercises/analysis/ground/ground.html
pdal translate 2663_1206.las -o 2663_1206_ground.laz outlier smrf range --filters.outlier.method="statistical" --filters.outlier.mean_k=8 --filters.outlier.multiplier=3.0 --filters.smrf.ignore="Classification[7:7]" --filters.range.limits="Classification[2:2]" --writers.las.compression=true --verbose 4
- Erstellen des DTM:
https://pdal.io/workshop/exercises/analysis/dtm/dtm.html
pdal pipeline pipeline.json
pipeline.json:
{ "pipeline": [ "2663_1206_ground.laz", { "filename":"2663_1206_dtm.tif", "gdaldriver":"GTiff", "output_type":"all", "resolution":"0.25", "type": "writers.gdal" } ] }
- Füllen der Lücken: https://gdal.org/programs/gdal_fillnodata.html
gdal_fillnodata.py -md 50 2663_1206_dtm.tif 2663_1206_dtm_fill.tif
- Erstellen eines Reliefs/Hillshade:
https://gdal.org/programs/gdaldem.html
gdaldem hillshade 2663_1206_dtm_fill.tif 2663_1206_dtm_hill.tif -z 0.5 -s 1.0 -multidirectional -of GTiff
On Thu, 13 Jan 2022 at 21:20, René Buffat buffat@gmail.com wrote:
Swisstopo hat den Kanton Luzern [1] in der Periode 2019-2020 mit Lidar abgeflogen. Die Punktwolke kann man bereits herunterladen[2] und im swissSurface3d hillshade sind die Daten auch schon enthalten[3]. Ich vermute, das swissAlti3D DTM wird in einem der nächsten Releases mit diesen Daten aktualisiert.
lg rene
[1] https://www.swisstopo.admin.ch/content/swisstopo-internet/de/geodata/height/... [2] https://www.swisstopo.admin.ch/de/geodata/height/surface3d.html [3] https://map.geo.admin.ch/?lang=en&topic=ech&bgLayer=ch.swisstopo.pix...
On Thu, 13 Jan 2022 at 18:46, Simon Poole simon@poole.ch wrote:
Das DTM scheint ein Bezahlprodukt zu sein https://geoportal.lu.ch/-/media/RAWI/Dokumente/Downloads/geoinformation/mb_d... Am 13.01.2022 um 14:49 schrieb Sentalize:
Das ist genial ... Luzern fiel mir immer als ziemlich rückständig auf, was Versatz usw angeht .. da half schon Swissimage enorm. Wie ich sehe, hat es auch ein Terrainmodell aus 2018. Das wär äusserst hilfreich. Mal sehen, ob man das irgendwo als WMS findet.
Am Donnerstag, 13. Januar 2022, 14:22:56 MEZ hat Simon Poole simon@poole.ch simon@poole.ch Folgendes geschrieben:
Eine gute Nachricht aus der Innerschweiz: der Kanton Luzern hat Anfangs Jahr auch angefangen seine Geodaten zu offen Bedingungen freizugeben.
Siehe https://geoportal.lu.ch/geodaten
Dies ist umso mehr bemerkenswert als der Kanton über Jahre hinweg massiv gegen Open Data opponiert hat (von den Gerüchten die jeweils gestreut wurden bezüglich OSM will ich gar nicht anfangen).
Natürlich gibt es auch hier eine Schattenseite: es sieht aktuell so aus als ob die Nutzungsbedingungen eine direkte Nutzung in OSM nicht erlauben, das übliche Quellenhinweisproblem. Falls aber jemand einen Datensatz speziell interessant findet, kläre ich das gerne ab.
Simon
talk-ch mailing list talk-ch@openstreetmap.ch http://lists.openstreetmap.ch/mailman/listinfo/talk-ch
talk-ch mailing listtalk-ch@openstreetmap.chhttp://lists.openstreetmap.ch/mailman/listinfo/talk-ch
talk-ch mailing list talk-ch@openstreetmap.ch http://lists.openstreetmap.ch/mailman/listinfo/talk-ch