Weihnachtsangebot 2025! 🎄

[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.

[Chris 1]

@chris-1

export "BROOK=10.2076 54.5257"

gdallocationinfo dwd-source-bw.tif -wgs84 % BROOK
-> value 231

Grau muss nicht grau sein, Kathegorisiert Blau-Stufen

Für mein "Schnorchel-Script" könnte ich aktuell also "aktueller Wellengang" auf fast 50cm setzen, eher panne

Nun müssen weitere Berechnungen erfolgen - ich will ja auch die Welle-Richtung und -Geschwindigkeit wissen

Obiges basiert übrigens auf "Peak Dünung" gemäß numerischen Modell CWAM +00:00 des DWD. Man kann da mit "mittlerer Windsee ....", sowie "Peak" und "Signifikant" beliebig ergänzen, wenn es um Extremwerte geht.
Aber ich will ja nur ruhige See :-)

@Armilar - wie würde ich den Return-Wert von gdallocationinfo in ibroker speichern ?
Denke da (vorläufig) an Trigger <25cm-> Schnorcheln=true

[Chris 1]

Wie bekomme ich den RETURN Wert von

In 0_userdata.Ostsee.Brook.Wert ?
Sprich, Return von letzter EXEC ANweisung

[DJMarc75]

@chris-1 sagte in Schnorchel Bedingungen Ostsee:

RETURN Wert

Variable "result" erstellen

[Chris 1]

@djmarc75 Danke

Habe mittlerweile beschlossen, das ganze einfach extern in einem Script zu machen,
speziell nachdem ich über "iobroker state set stateId newValue ack" stolperte. Das
gefällt mir, macht das Leben leichter

[Chris 1]

@chris-1 orrrrrr - lernen kann auch frustrierend sein

Jetzt habe ich nicht nur über re-klassifizierung von Karten vom Prompt aus (bzw aus'm Script) gelernt, NEIN, ich habe auch mehr über WMS (WebMappingService) und hier speziell "GetFeatureInfo" gelernt.

Alles, was ich wollte, bekomme ich ohne 1 Zeile Script zu schreiben frei Haus!, da bei der Feature-Abfrage nicht auf die Kartendarstellung zurück gegriffen wird, sondern auf die zugrunde liegenden Daten. Es ist also weit genauer!

Die URL für die Abfrage ist zwar monströs (wenn ich mehr gelernt hab, wirds hoffendlich kürzer), aber sieht derweil so aus (also, für meinen Lieblings-Strand:

https://maps.dwd.de/geoserver/dwd/ows?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureInfo&LAYERS=QUERY_LAYERS=dwd:Cwam_reg025_fd_sl_SWH&QUERY_LAYERS=dwd:Cwam_reg025_fd_sl_SWH&WIDTH=10&HEIGHT=10&X=1&Y=1&SRS=EPSG:4326&BBOX=11.0407,54.0141,11.0417,54.0151&INFO_FORMAT=application/json&FORMAT=text/text&TILESIZE=256&OVERVIEWCOUNT=17&MINRESOLUTION=0.0000053644180298&TILED=true

Bevor ich es vergesse: Wichtig ist, dass der Abfragewert (X und Y) innerhalb des "dargestellten Fensters" (WIDTH HEIGHT, hier also 1-10) liegt. X und Y sind keine geografischen Koordinaten!

Die BoundingBox (BBOX, Abgefragter Karten-Ausschnitt in geographischen Koordinaten (default EPSG:4326 bzw WGS84) sollte klein gehalten werden, das reduziert das Datenvolumen. Ich habe, wie man sieht, auf der 3ten Nachkomma-Stelle um 1 erhöht.

Jetzt muss ich nur noch mal schauen/lernen, wie ich hier im ioBroker aus dem JSON-Object (bei mir in 0_userdate.0.Ostsee.Brook.HWS) den Wert "Grayscale" auslese und freundlich runde

[Chris 1]

@nordischerjung Für's Angeln kannst Du auf diesem Wege auch viele Informationen abfragen

Windgescheindigkeit (über Wasser), Richtung, Wassertemperatur, Wellenrichtung,... und und und

Bei mir sollen im Sommerhalbjahr bei passenden Bedingungen Flossen im Dashboard erscheinen, ein Fisch wär bei Dir bestimmt ein guter hinweis :grin: