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

 ZEVS ForumVremenske postajeDržavne vremenske postaje po Sloveniji (Moderatorji: ŠTROMAR, mat89)Nove ARSO samodejne postaje
Strani: 1 ... 10 11 [12]
Natisni
Avtor Tema: Nove ARSO samodejne postaje  (Prebrano 162124 krat)
0 Člani in 1 Gost gledajo to temo.
gregach
Član VD ZEVS
Cumulus
*
Odsoten Odsoten

RD: 1991-01-09
Sporočila: 9035


Lokacija: Rudnica (405 m.n.v.), bližina Podčetrtka

Email
« Odgovori #165 dne: 12. Maj 2021, 20:59:45 »

Ja, bo šlo tudi hitreje. Smiley Označi vse podatke iz obeh stolpcev, ki ti jih da Arsov arhiv in jih najprej prekopiraj v Beležnico (WordPad), verjetno bi šlo tudi v Word. Potem pa daš komando, da program avtomatsko zamenja vse pike z vejicami, potem pa takšne podatke prekopiraš v Excel (če bi uvozil podatke direktno s pikami namesto vejicami, ti jih namreč Excel spremeni v datumsko obliko, to si najbrž že sam opazil). Potem pa se lahko igraš naprej, samo pazi na morebitne luknje v podatkih, da jih izločiš, sicer se zna zgoditi, da ti bo pri izračunu manjkajočih polurnih povprečij upošteval 0,0°C, to ti pa najbrž ne ustreza.

Sem tudi sam pisal na ARSO in so mi celo poslali podatke za december 2020 na 1 decimalko natančno, za ostalo so mi pa rekli naj si sam izračunam.  icon_haha
Bom najprej poizkusil na roke izračunat za 1 mesec, pa da vidim, koliko časa mi vzame, ker tole na računalniku ne vem, če mi bo ratalo.
Prijavljen

"Baba je tako kot oblak: če gre stran bo morda še lep dan."
komar
Altostratus

Odsoten Odsoten

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


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

« Odgovori #166 dne: 12. Maj 2021, 23:16:27 »

December 2020 +2,558467742°C? A je dosti natančno? Ne vem, če je trajalo več kot eno minuto Grin. Je pa zanimivo, da cel mesec ni niti ene luknje v podatkih. Komaj čakam, da se lotim še Lesc. Pa marsikdo se bo zdaj lotil še česa drugega, Pandorina skrinjica se je odprla... Grin

Če pa zbrišeš prvi podatek (01. 12. 2020 ob 00.00, ki paše še pod november 2020) in pa dodaš še 1. januar 2021 ob 0.00, ki pa paše pod december 2020, pa pride +2,560147849°C.

Če se ti ne mudi, vsaj toliko počakaj, da pogruntamo, kako podatke prekopirati v Kalkulator. icon_confused
« Zadnje urejanje: 13. Maj 2021, 00:35:05 od komar » Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

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

Lokacija: Tolmin

WWW
« Odgovori #167 dne: 13. Maj 2021, 13:21:30 »

Ko ti izpiše tabelo s podatki, je nad to tabelo gumb "Shrani". Naredi ti eno datoteko, ki jo nato shraniš k sebi.
Tukaj ti pripenjam enostavno python kodo, ki ti fajl prežveči, napiše povprečno temperaturo in še nariše graf.
Koda:
#!/usr/bin/env python
# -*- coding: utf-8 -*-

# Ta skripta prebere log file ARSO (polurni podatki) in izpiše vrednosti za izbrani dan.

import sys
import time
import datetime
from datetime import datetime
import matplotlib.pyplot as plt

if len(sys.argv) < 1:
    print "\nUsage:   python parseARSO.py <filename>\n"
    exit()

datoteka=sys.argv[1]

print "Berem datoteko: %s" % datoteka

cas=[]
temp=[]

for line in open(datoteka, 'r'):
  fields = line.split(",")
  if fields[0] != "station id":
    c=datetime.strptime(fields[2], "%Y-%m-%d %H:%M")
    c1 = (c - datetime(1970, 1, 1)).total_seconds()
    cas.append(int(c1))
    temp.append(float(fields[3]))

cas = [x - cas[0] for x in cas]
cas = [x/3600 for x in cas]

print ("Povprečna temperatura: %.1f °C" % float(sum(temp)/len(temp)))

plt.figure()
plt.plot(cas,temp)
plt.grid(True)
plt.ylabel('temperatura [C]')
plt.xlabel('cas [h]')
plt.show()
Prijavljen
komar
Altostratus

Odsoten Odsoten

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


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

« Odgovori #168 dne: 13. Maj 2021, 15:21:50 »

December 2016: -0,929099462°C
December 2017: +2,253158602°C
December 2018: +1,644758065°C
December 2019: +2,896169355°C
December 2020: +2,560147849°C

Pri decembru 2020 sta prvi in zadnji podatek ročno popravljena, ostalo sem pustil, tako kot ima ARSO. Manjkajočih podatkov oz. lukenj ni videti. Shocked

Za december 2019 sem porabil minuto in 38 sekund, sem pa imel oba datuma že vnesena od prejšnjega leta, kar seveda prihrani kakšnih 10-15 sekund. Plus dodatno rabiš še enih 20 sekund za formulo, ki ti izračuna povprečje. Chuck Norris ti pa baje res izračuna v manj kot minuti, samo to so informacije iz neuradnih virov... Grin

Ti pa Excel datoteko ob priliki pošljem še po e-pošti.

@Tminc: ti pa potem zdaj lahko narediš kontrolo... A mogoče to kodo programer lahko vnese v program Vremenko (Ognjena veverica)? Tam so tudi to včasih imeli krasno rešeno, zdaj pa so od tabel in grafov ostale samo še ruševine...
« Zadnje urejanje: 13. Maj 2021, 15:27:14 od komar » Prijavljen
gregach
Član VD ZEVS
Cumulus
*
Odsoten Odsoten

RD: 1991-01-09
Sporočila: 9035


Lokacija: Rudnica (405 m.n.v.), bližina Podčetrtka

Email
« Odgovori #169 dne: 13. Maj 2021, 15:34:02 »

Uuu, hvala obema!  beerchug  praise
Ja me boš moral naučit, da te ne bom več obremenjeval ... ker bom enkrat v prihodnosti rabil še podatke za januar in februar ...
Prijavljen

"Baba je tako kot oblak: če gre stran bo morda še lep dan."
komar
Altostratus

Odsoten Odsoten

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


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

« Odgovori #170 dne: 29. Julij 2021, 19:36:52 »

No, to me zanima, kaj se zgodi z izračuni, ko se med podatki pojavi takšna luknja, kot je npr. nastala danes med 7.30 in 13.30 uro na letališču Lesce... Saj če računaš za cel mesec, ne vpliva toliko na končni izračun povprečja, če pa računaš za en sam dan...
« Zadnje urejanje: 30. Julij 2021, 14:20:56 od komar » Prijavljen
komar
Altostratus

Odsoten Odsoten

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


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

« Odgovori #171 dne: 01. Oktober 2021, 12:11:13 »

Ko ti izpiše tabelo s podatki, je nad to tabelo gumb "Shrani". Naredi ti eno datoteko, ki jo nato shraniš k sebi.
Tukaj ti pripenjam enostavno python kodo, ki ti fajl prežveči, napiše povprečno temperaturo in še nariše graf.
Koda:
#!/usr/bin/env python
# -*- coding: utf-8 -*-

# Ta skripta prebere log file ARSO (polurni podatki) in izpiše vrednosti za izbrani dan.

import sys
import time
import datetime
from datetime import datetime
import matplotlib.pyplot as plt

if len(sys.argv) < 1:
    print "\nUsage:   python parseARSO.py <filename>\n"
    exit()

datoteka=sys.argv[1]

print "Berem datoteko: %s" % datoteka

cas=[]
temp=[]

for line in open(datoteka, 'r'):
  fields = line.split(",")
  if fields[0] != "station id":
    c=datetime.strptime(fields[2], "%Y-%m-%d %H:%M")
    c1 = (c - datetime(1970, 1, 1)).total_seconds()
    cas.append(int(c1))
    temp.append(float(fields[3]))

cas = [x - cas[0] for x in cas]
cas = [x/3600 for x in cas]

print ("Povprečna temperatura: %.1f °C" % float(sum(temp)/len(temp)))

plt.figure()
plt.plot(cas,temp)
plt.grid(True)
plt.ylabel('temperatura [C]')
plt.xlabel('cas [h]')
plt.show()

Če ima kdo kaj viška volje in časa ter obvlada te kode, se priporočam za 10-minutne odčitke padavin tekom sredine povodnji v Ljubljani-Bežigrad. Graf bi tudi prav prišel. V Vidmu (Dobrepolje) je namreč npr. v 5 minutah te nevihte padlo 12,4 mm.icon_wink
Prijavljen
Tminc
Altokumulus

Odsoten Odsoten

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

Lokacija: Tolmin

WWW
« Odgovori #172 dne: 01. Oktober 2021, 21:46:18 »

V ARSO arhivu so na voljo samo polurni podatki s postaje Bežigrad. Verjetno je časovna resolucija pregroba in graf ne bo pomenil kaj dosti.
Prijavljen
komar
Altostratus

Odsoten Odsoten

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


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

« Odgovori #173 dne: 01. Oktober 2021, 22:13:45 »

Uf, hvala, ravno sem še ujel na spletu. 25,1 mm je padlo v 10 minutah.
Prijavljen
Strani: 1 ... 10 11 [12]
Natisni
Pojdi na: