Novice: Trenutno ni obvestil.
Dobrodošel, Gost. Prosim, prijavi se ali se registriraj.

 ZEVS ForumVremeSplošnoArhiv radarskih slik
Strani: 1 [2] 3 4
Natisni
Avtor Tema: Arhiv radarskih slik  (Prebrano 50588 krat)
0 Člani in 1 Gost gledajo to temo.
Smučar
Stratus

Odsoten Odsoten

RD: 1981-07-05
Sporočila: 5520


Lokacija: Notranje Gorice pri LJ (300 m.n.v.)

Email
« Odgovori #15 dne: 09. Junij 2022, 15:25:18 »

Sem probal tudi na telefonu in v domačem omrežju, ne da se prenesti. S formatom ni nič narobe, morda samo nimaš dobro nastavljenega share-anja.
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #16 dne: 09. Junij 2022, 16:02:23 »

Prosim poskusite še enkrat. Spremenil sem končnico v .tar . Meni je delovalo v Firefoxu, ampak verjamem, da so lahko problemi.
Prijavljen
komar
Nimbostratus

Odsoten Odsoten

RD: 1978-10-13
Sporočila: 2012


Lokacija: Radovljica (491 m.n.v.)

« Odgovori #17 dne: 09. Junij 2022, 16:15:49 »

Ja, v Firefoxu mi je uspelo. Točno tako kot sta oba zapisala.

Po prenosu je treba odpreti mapo s programom 8 ZIP in datoteko razširiti s programom 7 ali 8 ZIP.

Na telefonu tudi prenos deluje, na starem (Android 4.4.2) mi je uspelo odpreti, na novem pa moram prej še namestiti eno od številnih ZIP aplikacij.

EDIT: ne morem verjeti; zdaj pa deluje celo v Internet Explorerju 11; tu žal samo še do srede, ko sledi pokop z vojaškimi častmi... Cry

Iz slike spodaj (Microsoft Edge) je razvidno, da naj bi bil problem v končnici datoteke .gz in je zato včeraj ni odpiralo.

Povezava do slik od 8. junija 2022 naprej: https://drive.google.com/drive/folders/13CoYa5FHaJPeTox9sWQ8W3ufXw4i_idI?usp=sharing

« Zadnje urejanje: 10. Junij 2022, 16:06:02 od komar » Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #18 dne: 17. Junij 2022, 12:32:09 »

Sporočam, da sem danes naredil majhne spremembe v zajemih. Ime datoteke in pripadajoča dnevna mapa se zdaj skladata s timestampom v radarski sliki. To pomeni, da bi morale slike v mapi 2022-06-17 vse pripadati današnjemu datumu, njihova imena pa se so v formatu <HHMM>: HH ... ura iz radarske slike, MM ... minuta iz radarske slike
Prej je HHMMSS pripadalo OS času in je kasnilo za 10-15 minut se mi zdi.
Današnje datoteke bodo imensko malo pomešane, bi pa morale vseeno bit notri vse.
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #19 dne: 18. Junij 2022, 22:32:12 »

Vse skupaj je še enkrat posodobljeno, jutri bo "testni" dan, če bo vse delalo OK.
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #20 dne: 12. Julij 2022, 08:05:54 »

Testno sem dodal sliko akumulacije padavin za pretekli dan. Ne vem, kako to delajo/te drugi, moja se naredi takole: preberem vsako sliko posebej (.srd) in preračunam akumulacijo v zadnjih 5 minutah, nato prištejem naslednjo, itd. do zadnje. Na koncu zaokrožim vrednosti na cela števila in narišem png. Skala mi še manjka, to še sledi. Ker vse skupaj teče na nekem remote počasnem raspberryu, mi testiranje vzame čas in delam vsak dan po malo. Čas, potreben za izdelavo ene take slike, je okrog 12 minut.

http://tolmin.zevs.si/razno/
Prijavljen
komar
Nimbostratus

Odsoten Odsoten

RD: 1978-10-13
Sporočila: 2012


Lokacija: Radovljica (491 m.n.v.)

« Odgovori #21 dne: 14. Julij 2022, 10:00:54 »

Slika akumulacije je zdaj dodana tudi na prvo stran teme Trenutno vreme (1. 7. 2022 - 31. 7. 2022), nekje na sredino strani. Čakam soboto, a se bojim, da @Tminc-a čaka še kar nekaj testiranja... icon_wink
Prijavljen
stog
Skrbnik Portala VD ZEVS
Nimbostratus
*
Odsoten Odsoten

RD: 1986-11-19
Sporočila: 2669


Lokacija: Ljubljana

Email
« Odgovori #22 dne: 17. Julij 2022, 16:35:11 »

Čestitam Tminc.  love praise clap
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #23 dne: 24. Julij 2022, 09:19:01 »

@stog: ti obvladaš te akumulacije. Jaz sem prišel do nekega rezultata, kot je v pripet. Gre za akumulacijo danes od polnoči do 0650 CEST. Prvi preračun mi da maksimume do 50 mm, vendar so bili realno okrog 20 (rečmo 25), kot si preračunal ti. Še enkrat sem preveril razlago, ki so mi jo podali na ARSO o preračunu iz dBZ v mm/h in mi vse matematično štima. Vrednosti iz posamezne slike pa pomnožim s faktorjem 5/60, ki je moj dt (5 min). V priloženi akumulaciji sem nato vse skupaj zdelil z 2. Kaj delam narobe? Postopek ali kaj drugega?
Kako akumulacije preračunaš ti?
Hvala.
Prijavljen
Mendos
Cirus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 35



« Odgovori #24 dne: 24. Julij 2022, 15:44:21 »

Zdravo,

predlagam, da za akumulacijo namesto maksimalne odbojnosti v dBZ uporabiš talno jakost padavin (RRG), kjer vrednosti pretvoriš iz dBR v mm/h, potem pa samo narediš vsoto vseh slik in jo deliš z 12.
Tole dodajam iz ARSO dokumentacije: RRG = 2.0*(n-64) - 8.0, n = 64..79 (https://meteo.arso.gov.si/uploads/meteo/help/sl/SRD3Format.html).

Primer pretvorbe med razredi:
Koda:
dBR = 10 log (R[1 mm]/1 mm)

quant=  rr
unit=   dBR
start=  -8.0
slope=  2.0
offset= 64
nlevels=16

level  +offset ASCII       middle   lower    upper  
                             dBR      dBR      dBR    

00 64 @    undef     -inf    -7.00   # clear-sky
01 65 A    -6.00    -7.00    -5.00
02 66 B    -4.00    -5.00    -3.00
03 67 C    -2.00    -3.00    -1.00
04 68 D     0.00    -1.00     1.00
05 69 E     2.00     1.00     3.00
06 70 F     4.00     3.00     5.00
07 71 G     6.00     5.00     7.00
08 72 H     8.00     7.00     9.00
09 73 I    10.00     9.00    11.00
10 74 J    12.00    11.00    13.00
11 75 K    14.00    13.00    15.00
12 76 L    16.00    15.00    17.00
13 77 M    18.00    17.00    19.00
14 78 N    20.00    19.00    21.00
15 79 O    undef    21.00     +inf

level  +offset ASCII       middle   lower    upper  
                             mm       mm       mm    

00 64 @    undef     -inf     0.20   # clear-sky
01 65 A     0.25     0.20     0.32
02 66 B     0.40     0.32     0.50
03 67 C     0.63     0.50     0.79
04 68 D     1.00     0.79     1.26
05 69 E     1.58     1.26     2.00
06 70 F     2.51     2.00     3.16
07 71 G     3.98     3.16     5.01
08 72 H     6.31     5.01     7.94
09 73 I    10.00     7.94    12.59
10 74 J    15.85    12.59    19.95
11 75 K    25.12    19.95    31.62
12 76 L    39.81    31.62    50.12
13 77 M    63.10    50.12    79.43
14 78 N   100.00    79.43   125.89
15 79 O    undef   125.89     +inf

Intenziteto padavin RR (mm/h) se izračuna po formuli: RR = 10^(x/10). Prvi razred oz. level 00 (kjer velja RRG = -8 dBR) mora biti 0 mm/h.

Akumulacija iz talne jakosti v primerjavi s tisto na podlagi maksimalne odbojnosti sicer pogosto podcenjuje padavine, sploh na hribovitih območjih, vendar je po drugi strani pri močnejših padavinah nekoliko bolj realna. V poročilih ARSO so akumulacije narejene iz talne jakosti (primer): https://meteo.arso.gov.si/uploads/probase/www/climate/text/sl/weather_events/neurje_15jul2022.pdf.

Pri akumulaciji iz maksimalne odbojnosti dobimo boljšo razporeditev padavin, saj je slika sestavljena iz največjih zaznanih odbojev po vertikali nad določeno točko ne glede na njihovo višino nad tlemi (boljše izpade predvsem nad območjem Kamniško-Savinjskih Alp, Karavank in Julijskih Alp), so pa količine padavin zelo pogosto precenjene.

Pri maksimalni odbojnosti dBZ v mm/h izračunaš po formuli: RR = (10^((x-1.5)/10)/250)^(2/3). Pri tem je dobro anulirat nekaj najnižjih razredov odbojnosti oz. odvisno od namena (jaz recimo izločim prve tri, do vključno 24 dBZ oz. intenziteto pod 1,2 mm/h). Šibke padavine vsaj v primeru poletnih konvekcijskih padavin do tal večinoma izhlapijo, zato jih je najbolje kar zanemariti. Nadaljevanje je enako (vsota vseh slik/12).

Dodajam primerjavo današnje akumulacije do 8. ure. Prva karta je narejena na podlagi talne jakosti (RRG), druga na podlagi prilagojene maksimalne odbojnosti (ZM).

Upam, da je kaj v pomoč.  icon_smile
« Zadnje urejanje: 24. Julij 2022, 15:48:37 od Mendos » Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #25 dne: 24. Julij 2022, 16:33:26 »

@Mendos: hvala za opis in nasvet. Res uporabljam maksimalno odbojnost. Preklopil bom na talno odbojnost.
Zakaj omenjaš, da sešteješ vse slike in deliš z 12? Si mislil 12 slik iz zadnje ure? Jaz naredim načeloma enako, ker vsako "sliko" (podatke) delim z 12 in seštevam. Je to v redu?
Imam pa še (najmanj) dva problema tehnične narave... kje dobiti umerjen zemljevid za točno območje radarja? :-) In še - kako lepše narisati podatke / barvne razrede. Tu je moje znanje čisto na dnu, ker uporabljam enostaven pythonov matplotlib in plt.imshow. Če ne želiš ali imaš časa, se ne sekiraj s to razlago. Bom skušal sam najti rešitev.
Prijavljen
Mendos
Cirus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 35



« Odgovori #26 dne: 24. Julij 2022, 18:58:32 »

Vsoto za katerokoli časovno obdobje (število slik) deliš z 12. Saj v bistvu je vseeno, lahko tudi vsako sliko posebej in potem sešteješ. Oboje je prav, rezultat je enak.

Pri zemljevidu je pomembno, da so podatki in podlaga v istem koordinatnem sistemu. Najbolje je naredit transformacijo radarske slike v kakšno splošno uporabljeno projekcijo, recimo WGS 84 / Pseudo-Mercator (EPSG:3857) – ta se uporablja za spletne zemljevide Google Maps in OpenStreetMap. S tem lahko kombiniraš sliko z različnimi podlagami. Ujemanje na približno preveriš na tem zemljevidu https://meteo.arso.gov.si/uploads/meteo/app/inca/?par=si0zm. Transformacija sliko malo popači (zgoraj je bolj raztegnjeno, robovi so ukrivljeni).

Glede izrisa in grafike je v bistvu čisto poljubna izbira. Jaz vse skupaj sestavim kar v programskem jeziku R z nekaj osnovnimi knjižnicami za izrisovanje. Verjamem, da ima matplotlib podobne možnosti. Uporaben bi znal bit tudi paket cartopy, čeprav nimam izkušenj z risanjem zemljevidov v Pythonu. Razredi barvne skale so čisto poljubni, lahko se tudi preliva in še malo zgladi piksle. Najbolj enostavno je kar pustit originalne vrednosti in naredit dovolj široke razpone, da ni preveč odtenkov. Pa še glede na možna velika odstopanja od realnega stanje je bolje prikazovat malo manj podrobnosti. icon_mrgreen
Prijavljen
stog
Skrbnik Portala VD ZEVS
Nimbostratus
*
Odsoten Odsoten

RD: 1986-11-19
Sporočila: 2669


Lokacija: Ljubljana

Email
« Odgovori #27 dne: 25. Julij 2022, 22:03:18 »

@stog: ti obvladaš te akumulacije. Jaz sem prišel do nekega rezultata, kot je v pripet. Gre za akumulacijo danes od polnoči do 0650 CEST. Prvi preračun mi da maksimume do 50 mm, vendar so bili realno okrog 20 (rečmo 25), kot si preračunal ti. Še enkrat sem preveril razlago, ki so mi jo podali na ARSO o preračunu iz dBZ v mm/h in mi vse matematično štima. Vrednosti iz posamezne slike pa pomnožim s faktorjem 5/60, ki je moj dt (5 min). V priloženi akumulaciji sem nato vse skupaj zdelil z 2. Kaj delam narobe? Postopek ali kaj drugega?
Kako akumulacije preračunaš ti?
Hvala.
Vsak padavinski dogodek je svoja zgodba. Včasih imamo celo znotraj padavinskega dogodka dve zgodbi, ki bi potrebovali vsaka svoje vrednotenje barv.  Wink
Res je. Niti približno ni tako enostavno. Včasih ti bo ista akumulacija dala 50mm pa je padlo v resnici 25mm, drugič bo kazala 50mm pa bo padlo 60mm. Tako se dogaja na arso in osmer akumulacijah. Ter vseh ostalih. So pač avtomatično generirane. Moja edina ni.  icon_cool Moj način je edini ki omogoča kolikor toliko točno sezonsko seštevanje poslabšanj, ker če bi za 2 meseca združil vse radarske slike bi bila natančnost vsaj 50% slabša.

Če še malo zakompliciram, pozimi bi ti isto vrednotenje barv, kot je primerno za poletje dalo situacijo, da bi v resnici padlo 100mm, karta bi ti pa kazala 10mm.  Grin

Jaz imam 5 različnih vrednotenj barv in jih uporabljam po občutku. Za večino poletnih dogodkov se zadnje čase dobro izkažeta tisti dve, kjer akumulacija ignorira barve do svetlo modre ali do temno zelene. Prav tako je problem z vijolčno barvo saj je velika razlika, če padajo teniške žogice ali pa neumorno zliva. Ne vemo.

Jaz torej izberem uteži, primerne za tisto poslabšanje, ter končni rezultat skupaj s celo karto priredim kolikor je mogoče da se sklada s padlimi padavinami na arso in drugih postajah. Na nekaterih bo karta kazala preveč, drugih premalo. Če nimam na kaj vpeti pač vsekam povprečje.
Oddaljena območja radarja so večna težava ampak tu nimamo kaj. Teoretično je možno to ročno popravljati kar izjemoma kdaj naredim ampak če nisi nežen vse pokvariš. In zelo je zamudno.

Težko ti povem točno pravilno razporeditev uteži za določene barve, to boš moral preračunati sam, ker ne obstaja absolutno pravilna, ampak iz napisanega boš morda bolje razmišljal. Vsekakor pa moje uteži upoštevajo razmerja, ki jih navaja arso med posameznimi barvami. Posegam le v vrednosti prvih par barv in zadnje.

Močno odsvetujem akumulacije iz talnih jakosti ker boste imeli luknje v kartah.  icon_arrow icon_arrow  icon_surprised
Poskusite pa boste videli. Morda za kakšno območje res pridobite na natančnosti a celotno karto spremenite v sir.  icon_lol Še večja past je, ker se za kakšna območja ne bo videlo da je rezultat "sirovih lukenj"  Grin shrug
« Zadnje urejanje: 25. Julij 2022, 22:07:18 od stog » Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #28 dne: 26. Julij 2022, 21:38:03 »

@stog: hvala tudi za tvojo razlago in priporočila. Zdaj sem res preklopil na talno odbojnost. Bom videl, kako bo. Meni je to bolj v zabavo, vseeno pa mora biti vsaj malo uporabno.
Sem pa zdaj verjetno dosegel svoj maksimum v risanju kart. Moral bi prešaltat na projekcijo, ampak tega ne bom osvojil.
Me pa fascinira, kako dobro imaš ti in tudi Mendos to predstavljeno in verjetno arhivirano.

Lp.
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

RD: 0000-00-00
Sporočila: 629

Lokacija: Tolmin

WWW
« Odgovori #29 dne: 02. Avgust 2022, 18:18:07 »

Spremenil sem čas za generacijo slike akumulacije. Zdaj je termin skladen z ARSO meritvami: 8:00 - 20:00 ter 20:00 - 8:00. Slike se naložijo okrog 16 minut čez termin.
Prijavljen
Strani: 1 [2] 3 4
Natisni
Pojdi na: