Schnorchel Bedingungen Ostsee

[Chris 1]

@chris-1 Kurzes Rekapitulieren für morgen,
bevor ich gleich einschlafe...

curl https://www.dwd.de/DWD/wetter/wv_spez/seewetter/ostsa_00.png | convert - -remap dwd-ostsee-legende-quer.png dwd-ostsee.png

gdal_translate -of GTiff -a_ullr 3225261.8474067864008248 5742205.9424629490822554 4689520.2210977496579289 6869384.3093071654438972 -a_srs EPSG:31467 dwd-ostsee.png dwd-ostsee-gk3.tif

gdalwarp -s_srs EPSG:31467 -t_srs EPSG:4326 dwd-ostsee-gk3.tif dwd-ostsee-wgs84.tif

Zwischenstand: Georeferrenzierung & ColorMap & WGS84

Ziel:

• kleiner Kartenausschnitt
• Farben in Text-Matrix (Tabelle) überführen
• Blaualgen >= 18/25°C Wassertempoeratur

Im letzten Schritt könnte es noch spannend werden, die vorhandenen (etwas groben) Daten an die doch recht akkorat vorhandene Küstenlinie an zu passen, aber das müsste mit masking funktionieren - hoffendlich

Und wenn das alles klappt wie geplant..... wär natürlich ein Traum, wenn ein Java-Guru daraus einen Adapter bauen würde

FEHLER: 2ter cmd als "-of GTiff" definiert, aber File Typ png // korrigiert

[Armilar]

@chris-1 sagte in Schnorchel Bedingungen Ostsee:

wär natürlich ein Traum, wenn eijn Java-Guru daraus einen Adapter bauen würde

Den kannst du lange suchen - aber evtl. kann ja jemand JavaScript/TypeScript

[Chris 1]

@armilar Ich muss jetzt mal eben Schnitzel machen.

Junior war letztes mal begeistertm, wie ich gemäß "Wiener Schnitzel" Blasen in die Pannade machte lach

[Chris 1]

@chris-1 Kartoffel (für Stampf) kochen, Spinat in der Mache, Schnitzel sind geklopft

[Chris 1]

@chris-1 Das ganze endet ja nicht mit dem aktuellen Bericht.... für "Schnorcheln" muss ich ja die Daten Aufzeichnen, um dann über einen Zeitraum X (3-5 Tage?) die Werte zu vergleichen

Es muss so etwas wie eine "Unterwasser-Sicht-Prognose" erstellt werden

Prognose != Vorhersage, nur vermuteter Zustand

[Chris 1]

@chris-1 Seufz.... ob es die Euphorie war, Unachtsamkeit oder der Rotwein ..... habe gestern Abend leider Fehler gemacht. Musste ich feststellen, als ich mit den Postings von gestern weiter machen wollte.
Ich habe wohl mit einer Datei aus dem GIS weiter gearbeitet und NICHT mit einer Datei aus dem Scripting

Mit etwas Ruhe und mehr Aufmerksamkeit hab ich das ganze nochmal durch gearbeitet, so dass es per Script und ohne QGIS Hilfe durchlaufen kann

Ich ergänze hier gleich das aktuallisierte Script

Wie füge ich hier eine Datei hinzu? (will GeoRef-Datei (ASC, "georef.gcp") hinzufügen

Krücke.... folgenden ASC in Textdatei kopieren.
Bitte daruf achten, keine Zeilenumbrüche o.ä. zu ergänzen.
Es ist "ein" String (TXT-Variable)

#!/bin/bash

echo \
-gcp 30.3072 390.49 3.36699e+6 5.91262e+6 -gcp 68.4839 391.037 3.44856e+6 5.91899e+6 -gcp 72.221 380.293 3.45411e+6 5.94344e+6 -gcp 170.085 372.703 3.66034e+6 5.9779e+6 -gcp 161.793 339.419 3.63486e+6 6.04714e+6 -gcp 160.742 358.105 3.6365e+6 6.00955e+6 -gcp 239.279 340.938 3.802e+6 6.05844e+6 -gcp 254.811 375.505 3.84268e+6 5.99003e+6 -gcp 424.198 344.523 4.1991e+6 6.0889e+6 -gcp 516.375 340.782 4.39305e+6 6.11356e+6 -gcp 466.396 296.489 4.27881e+6 6.19369e+6 -gcp 485.325 95.4487 4.27779e+6 6.63218e+6 -gcp 471.059 105.918 4.24861e+6 6.60753e+6 -gcp 506.377 98.4398 4.32258e+6 6.63054e+6 -gcp 625.791 54.9539 4.56989e+6 6.75027e+6 -gcp 624.64 37.6976 4.56455e+6 6.78602e+6 -gcp 407.844 32.4057 4.09683e+6 6.75037e+6 -gcp 368.729 29.8748 4.01465e+6 6.74996e+6 -gcp 358.376 26.4235 3.98989e+6 6.75469e+6 -gcp 89.422 314.906 3.47678e+6 6.08705e+6 -gcp 77.8524 321.75 3.45418e+6 6.06856e+6 -gcp 87.7925 298.611 3.46261e+6 6.10329e+6 -gcp 150.855 283.456 3.60196e+6 6.1661e+6 -gcp 149.552 273.679 3.59723e+6 6.18377e+6 -gcp 159.296 171.293 3.59739e+6 6.40366e+6 -gcp 111.605 198.165 3.50062e+6 6.33771e+6 -gcp 98.0099 200.608 3.47617e+6 6.3336e+6 -gcp 87.2823 220.257 3.45151e+6 6.28594e+6 -gcp 28.214 90.0733 3.304e+6 6.54882e+6 -gcp 144.788 106.376 3.55712e+6 6.53888e+6 -gcp 50.3229 144.564 3.35947e+6 6.44221e+6 -gcp 179.403 102.244 3.62009e+6 6.54268e+6 -gcp 186.103 115.42 3.6429e+6 6.53302e+6 -gcp 339.554 195.706 3.98539e+6 6.38972e+6 -gcp 321.467 255.836 3.9591e+6 6.25412e+6 -gcp 371.004 218.357 4.05566e+6 6.34698e+6 -gcp 367.419 186.499 4.04477e+6 6.41314e+6

[Chris 1]

@chris-1 Die aktuallisiert Version des SCript siegt nun so aus.....

#!/bin/bash

curl https://www.dwd.de/DWD/wetter/wv_spez/seewetter/ostsa_00.png | convert - -remap dwd-ostsee-legende-quer.png dwd-source.png

GCP=source georef.gcp

gdal_translate -of GTiff $GCP "dwd-source.png" "dwd-conv.tif"

gdalwarp -r near -order 1 -co COMPRESS=NONE -t_srs EPSG:31467 "dwd-conv.tif" "wellengang.tif"

Das Bild sollte nun inetwa so aussehen:

fooo... scheint, ich muss einen SCreenshot machen, , weil das System hier mit den Modifikationen über fordert ist...

[Chris 1]

@chris-1 Wir schauen hier übrigens gerade auf eine Karte,
welche im Prinzip "verbogen" ist

[Chris 1]

@chris-1 Als nöchstes werden wir unser Untersuchungs-Gebiet festlegen und dies aufarbeiten

[Chris 1]

@chris-1 @Armilar - ist dies bis hierhin nachvollziehbar oder bin ich über manches zu schnell hinweg gegangen?

[Armilar]

@chris-1

Sehe ich mir an wenn es fertig ist, oder wenn du nicht weiter kommst. Ich denke dann sollte noch eine Zusammenfassung kommen.

[Chris 1]

@armilar Eine Zusammenfassung wird gewiss folgen.
Würde mich nur interessieren, ob Du es bis hierhin technisch nachvollziehen kannst

[Armilar]

@chris-1 ist doch super beschrieben

[Chris 1]

@armilar Was ich gewiss morgen oder übermorgen benötigen werde, ist eine Überführung der Werte in.... denke JSON Tabelle macht hier am meisten Sinn.

Also ich lese die pixel-Werte aus und DU sagst mir, wie ich die in eine JSON Tabellé packe

[Armilar]

@chris-1 bekommen wir hin

[Chris 1]

@armilar said in Schnorchel Bedingungen Ostsee:

@chris-1 bekommen wir hin

Etwas schwierig werden noch, dass "Land" den selben Farbcode (weiss) verwendet wie "3m Welle"und das Masking zu Schwarz-Werten (Wind-Pfeile)

[Chris 1]

@chris-1 Zomm auf süd-westliche Ostsee, um einen Hint zu geben, worum ich mich kümmern muss. Vielleicht hat iwer gute Ideen....
Hier Lübecker Bucht/Fehmarn bis Darss.

[Chris 1]

@chris-1 Wer schonmal in die Zukunft schauen will....

nach Wellhöhe (graphische Auswertung) erfolgt die Auswertung von Wellengeschwindigkeit & Richtung auf Basis von Text-Daten im GRIP-Format (https://filebox.bsh.de/index.php/s/Z8k9NMB3dKVpOXA/download?path=%2Fgrib1%2FOstsee&files=Current_ba_2022111100_00.grb.bz2)

Dann habwen wir "Wellenhöhe", "Wellengeschwindigkeit" und "Wellen-Richtung"

[Chris 1]

Oh fein, ich habe einen WMS Service des BSH gefunden.
Bedingt zwar immernoch umwandlung/Kategorisierung von Farbwerten in zugeordnete Wellenhöhen, aber ist schon deutlich feiner als diese kleine Karte. Wellenhöhe <=25cm wird klassifiziert und die Auflösung ist weit feiner

Basis:
https://marineforecast.bsh.de/Seewetter/#
(Erst "Themen" alles deaktivieren, Dann "Aktuell" "Wellenhöhe aktivieren)

Ausschnitt & Download (Beispiel-Code süd-westliche Ostsee, Lübeck/Wismar) für Download:
gdal_translate -of GTiff -outsize 1024 800 "WMS:https://maps.dwd.de/geoserver/dwd/ows?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=dwd:Cwam_reg025_fd_sl_SWH&SRS=EPSG:4326&BBOX=10.5,53.8,12,54.5&FORMAT=image/jpeg&TILESIZE=256&OVERVIEWCOUNT=17&MINRESOLUTION=0.0000053644180298&TILED=true" dwd-source.tif
Da müsste jetzt, hoffendlich, nur noch eine Kathegorisierung erfolgen

[Chris 1]

@chris-1 Viele Wege führen zum Ziel, nicht alle sind geradninig

---

gdal_translate -of GTiff -outsize 1024 800 "WMS:https://maps.dwd.de/geoserver/dwd/ows?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&LAYERS=dwd:Cwam_reg025_fd_sl_SWH&SRS=EPSG:4326&BBOX=6,53,12,55&FORMAT=image/jpeg&TILESIZE=256&OVERVIEWCOUNT=17&MINRESOLUTION=0.0000053644180298&TILED=true" dwd-source.tif

gdal_calc.py -A dwd-source.tif --A_band=1 -B dwd-source.tif --B_band=2 -C dwd-source.tif --C_band=3 --outfile dwd-source-bw.tif --calc=A0.2989+B0.5870+C*0.1140

Die GrauWert Zuordnung werde ich ergänzen, aber wenn Ihr Euch "dwd-source-bw".tif an schaut, dann sehr Ihr, dass es nicht ganz schlecht ist.