SwissMetNet Contourmap Applikation/Karten

[Bernhard Oker]

Erstellt mit Testdaten von MeteoSchweiz vom 20.06.2021:
https://www.meteoswiss.admin.ch/weather ... twork.html
https://www.meteoswiss.admin.ch/dam/jcr ... g-2023.zip





Gruss
Bernhard

[helios]

Erstellt mit Testdaten von MeteoSchweiz vom 20.06.2021:
https://www.meteoswiss.admin.ch/weather ... twork.html
https://www.meteoswiss.admin.ch/dam/jcr ... g-2023.zip

Sieht gut aus, hast du pyrad benutzt?

[Bernhard Oker]

Sieht gut aus, hast du pyrad benutzt?

Danke. Ich benutze aktuell nur h5dump um die Daten zu dekodieren. Den Rest programmiere ich selber in Delphi (Pascal). Einer der nächsten Schritte wird die Darstellung in globalen WGS84 Koordinaten sein.

Code: Alles auswählen

  // https://www.swisstopo.admin.ch/dam/en/sd-web/KLRCX9XIdXDu/ch1903wgs84-EN.pdf
  procedure ConvertGradToCHKoord(var chLat : double;
                                 var chLon : double;
                                 const gradLat : double;
                                 const gradLon : double);
  var
    p, l : double;
  begin
    p := ((gradLat * 3600.0) - 169028.66) / 10000.0;  // Breitenkoordinaten
    l := ((gradLon * 3600.0) - 26782.5) / 10000.0;    // Längskoordinaten

    chLat := 200147.07
           + (308807.95 * p)
           + (3745.25 * Power(l, 2.0))
           + (76.63 * Power(p, 2.0))
           - (194.56 * Power(l, 2.0) * p)
           + (119.79 * Power(l, 3.0));
           
    chLon := 600072.37
           + (211455.93 * l)
           - (10938.51 * l * p)
           - (0.36 * l * Power(p, 2.0))
           - (44.54 * Power(l, 3.0));
  end;

  procedure ConvertCHKoordToGrad(var gradLat : double;
                                 var gradLon : double;
                                 const chLat : double;
                                 const chLon : double);
  var
    p, l : double;
  begin
    p := ((chLon * 1000.0) - 600000.0) / 1000000.0;  // Breitenkoordinaten
    l := ((chLat * 1000.0) - 200000.0) / 1000000.0;  // Längskoordinaten

    gradLon := 2.6779094
            + 4.728982 * p
            + 0.791484 * p * l
            + 0.1306 * p * Power(l, 2.0)
            - 0.0436 * Power(p, 3.0);

    gradLat := 16.9023892
            + 3.238272 * l
            - 0.270978 * Power(p, 2.0)
            - 0.002528 * Power(l, 2.0)
            - 0.0447 * Power(p, 2.0) * l
            - 0.0140 * Power(l, 3.0);

    gradLat := gradLat * 100.0 / 36.0;
    gradLon := gradLon * 100.0 / 36.0;
  end;

Gruss
Benhard

[Bernhard Oker]

Auf der Hauptseite der Contourmap Karten wird neu der Gewitter Typ der Soundings farbig angezeigt.
Schwarz = Keine Gewitter
Grün = Einzelzellen Gewitter
Türkis = Multizellen Gewitter
Gelb = Low-Precipitation Supezellen
Rot = Classic Superzellen
Pink = High-Precipitation Superzellen

Beispiel: 1600Z 03.02.2026


Disclaimer:
Das Ermittlung des Gewitter Typ ist sehr experimentell und mit Vorsicht zu betrachten.

Gruss
Bernhard

[Bernhard Oker]

Wenn man die vor Jahren (2017) erstellte Relief Karte endlich brauchen kann:






Gruss
Bernhard

[Bernhard Oker]

Eine zusätzliche Version der Contourmap Karte mit optimierter Radar Niederschlag Darstellung (Enhanced Version).







Gruss
Bernhard

[Rontaler]

Hi Bernhard und helios

Ich danke euch für eure Programmierarbeiten und hoffe, dass die Experimente bis zur nächsten Gewittersaison "marktreif" sein werden!

[Bernhard Oker]

Ich habe die Enhanced Version optimiert. in der 1. Version wurden nur die Maximalwerte speziell berücksichtigt. Neu gibt es eine Berechnung von Min-Max Niederschlags Intensität.


Gruss
Bernhard

[Bernhard Oker]

Ich bin am experimentieren mit QGIS und Swisstopo/EU Daten für Grenzen/Gewässer um ein Hintergrundbild für den ganzen Radarbereich zu erstellen.


Gruss
Bernhard

[Bernhard Oker]

Coming soon:

Dateninterval = 1%


Daten ab 2cm Korngrösse; Dateninterval = 0.1cm


Gruss
Bernhard