Chat-Fenster: Wetter und Off-Topic

[helios]

Die Koordinaten der Radarstandorte sind mir natürlich bekannt.
Aber anhand der Artefakte im Radarbild kann man erkennen dass das Radarbild nicht exakt mit den tatsächlichen Koordinaten übereinstimmt und zumindest bei Lema ein paar wenige km abweicht.

Wenn ich das jetzt richtig verstehe, hat es vermutlich damit zu tun, dass LV95 (radardaten) eine Zylinderprojektion ist, während WGS84 (Karte) sphärisch ist.

Wobei der Fehler da nicht so gross sein dürfte. Möglicherweise bin ich von LV03-Koordinaten ausgegangen.

[Federwolke]

Link zum Datenblatt geht nicht mehr.

Danke für den Hinweis. Betrifft alle Links weltweit, dabei habe ich die gerade erst im letzten Jahr alle aktualisieren müssen

Schweiz habe ich geflickt. Wann ich Zeit für den Rest habe, weiss ich nicht - wird wieder eine wochenlange Fleissarbeit
PS: Dank Lukas' Testversion weiss ich jetzt, dass im ganzen Atlas 856 Radarstationen erfasst sind. Na dann viel Vergnügen beim Links flicken...

[helios]

Aktuelle Radardaten von MeteoSchweiz OpenData mit tiles von geo.admin.ch
https://rrdstats.ddns.net/openrad/
(Beta)

[Bernhard Oker]

mit tiles von geo.admin.ch

Das ging ja schnell.

Wie berechnest Du die PNG welche die Radardaten enthalten? Aktuell komme ich da nicht weiter.

CH-Koordinaten:


WGS84 mit Mercator-Projektion:


Vergleich:


Grus
Bernhard

[helios]

Ich habe alles in python gemacht. Für das importieren von Wetterradardaten gibt es wradlib.
Das einlesen des h5 files geschieht mit wradlib.io.hdf.read_generic_hdf5
https://docs.wradlib.org/en/2.2.0/noteb ... _hdf5.html
und das erzeugen der Grafik danach mit matplotlib.


wradlib kann übrigens auch die Metadaten auslesen. Darin enthalten sind der Status der Radare, die Boundaries für die Karte (falls sich das in anderen Radarprodukten ändert) und die Projektion.

Code: Alles auswählen

[WHERE] Attributes:
  LL_lat: 43.62900161743164
  LL_lon: 3.1687800884246826
  LR_lat: 43.61899948120117
  LR_lon: 11.955599784851074
  UL_lat: 49.3744010925293
  UL_lon: 2.689419984817505
  UR_lat: 49.36330032348633
  UR_lon: 12.462300300598145
  projdef: b'+proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +k_0=1 +x_0=2600000 +y_0=1200000 +ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs'
  xscale: 1000.0
  xsize: 710
  yscale: 1000.0
  ysize: 640

[WHAT] Attributes:
  date: b'20260206'
  object: b'COMP'
  source: b'ORG:215, CTY:644, CMT:MeteoSwiss (Switzerland)'
  time: b'174000'
  version: b'H5rad 2.4'

[HOW] Attributes:
  nodes: b'WMO:06661,WMO:06699,WMO:06768,WMO:06726,WMO:06776'

Ein Hindernis für die Nützlichkeit ist derzeit die Aktualität der Daten bei Opendata. Es gibt manchmal Verzögerungen von über 20 Minuten, dafür kommen mehrere Bilder in einem Rutsch.

[Bernhard Oker]

Habe die Projektion endlich hinbekommen:

Update: Nach einer Korrektur stimmt es jetzt perfekt.




Wenn man die Berechnung mit einem Intervall von 500m statt 1000m macht, verschwinden die leeren Stellen die durch die Mercator Umrechnung entstehen.


Testapplikation:


Gruss
Bernhard

[helios]

Update von meiner Seite:

Die grafana seite für die Messwerte von Meteoschweiz ist jetzt erreichbar unter https://grafana.openrad.ch

Der Wetterradar mit den Radardaten von Meteoschweiz Opendata ist erreichbar unter https://openrad.ch

Das Interface ist aufgeteilt in ein Realtime Tab und einen Archiv Tab (bis 8 Stunden zurück).
Im Archiv Tab kann bei den ab jetzt neu auffüllenden Daten mit einem Mouse-Over die Niederschlagsrate angezeigt werden:




Auf Mobilgeräten muss man auf die Stelle tippen, die man "messen" will.

Zwischenfrage: Wie hoch kann die Regenrate in einer starken Gewitterzelle steigen? Was sind da realistische Werte? Kann ich leider unter den aktuellen Umständen nicht testen

[Bernhard Oker]

Zwischenfrage: Wie hoch kann die Regenrate in einer starken Gewitterzelle steigen?

Das frage ich mich auch noch, was die maximale Regenrate ist, welche von MeteoSchweiz geliefert wird. Ich habe bis jetzt dazu keine Definition gefunden.

Bei Meteoradar wird 380 mm/h als Maximalwert auf der Skala angezeigt, daher bin ich auf meiner Grafik bis 400 m/h gegangen.
https://www.metradar.ch/2009/tgn/aktuell.php

Bei Kachelmannwetter wohl ähnlich hoch.
https://kachelmannwetter.com/ch/regenradar/schweiz

Wenn es diese Gewitterzelle aber nur auf 118.4 mm/h schafft, nehme ich an dass dort Schluss ist. Habe kurz die Stunde mit den Testdaten geprüft und überall war dort Schluss. Könnte sich aber inzwischen geändert haben...


R Max = Radarbild Maximalwert
S Max = Side-View Maximalwert

Update: ich habe eine Anfrage an MeteoSchweiz gestellt.

Gruss
Bernhard

[helios]

Ich arbeitete gerade an der Vorhersage der Wolkenbewegungen mit https://pysteps.github.io/
Aber 8 Stunden Radarausfall ist schon etwas viel zum überbrücken
Leider ist die prognose ziemlich rechenaufwändig, wohl etwas zu viel für den kleinen vserver.

[Bernhard Oker]

Wenn ich die Projektion von "LV95" auf "WGS84 Mercator" ändere, gehen durch das Komprimieren Datenpunkte verloren. Daher versuche ich gerade herauszufinden, wie stark man die Grösse der Daten (= des Bildes) ausdehnen muss, damit bei der Umwandlung zumindest jeder Datenpunkt des LV95 einmal im WGS84 Mercator vorkommt.







Gruss
Bernhard