Pagina 13 di 22 PrimaPrima ... 31112131415 ... UltimaUltima
Risultati da 121 a 130 di 220
  1. #121
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da leccese_a_pavia Visualizza Messaggio
    Hai perfettamente ragione, ma dubito che andrà a finire così, ahimè...

    Non in questo Paese. Valido o non valido che sia il metodo, sarà organizzato un coro di derisione e beffeggio. La cosa sarà presa sul serio e controllata come si deve all'estero. Ma naturalmente sui massmedia (e sulla letteratura scientifica) in Italia uscirano solo le notizie sull'eventuale fallimento.

  2. #122
    Brezza leggera
    Data Registrazione
    12/12/05
    Località
    Pavia ovest (PV)
    Età
    42
    Messaggi
    368
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da Borat Visualizza Messaggio
    La cosa sarà presa sul serio e controllata come si deve all'estero.
    E' questa la cosa che mi fa sperare...perchè, se il metodo fosse esatto, sarebbe una cosa troppo grossa!


  3. #123
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da leccese_a_pavia Visualizza Messaggio
    E' questa la cosa che mi fa sperare...perchè, se il metodo fosse esatto, sarebbe una cosa troppo grossa!

    Chiaramente se il metodo funziona, sarà sviluppato e implementato altrove, e poi noi ce lo dovremo comprare da imprese e dipartimenti stranieri, pagandolo a peso d'oro.
    Come successe a suo tempo con la radio e una sua importante applicazione: i radar. Durante la II Guerra Mondiale, siccome il nostro paese in precedenza aveva sostanzialmente snobbato Marconi, ci trovammo ad affrontare navi e aerei nemici dotati di radar, e noi senza (con risultati terribilmente noti dei vari scontri specie navali). Se ricordo bene, pare che Marconi avesse offerto alle nostre forze armate i suoi strumenti.
    Devo dire che non so se all'epoca Marconi fosse stato cacato poco perché non-accademico (e quindi la sua invenzione era patagonia, pregiudizi etnico-culturali, cabala, magia nera, non-scientifica, non controllata) o per altre ragioni. Ma ci sono pochi dubbi che dire di no al radar costò un'ecatombe.
    Mi viene in mente una interessante analogia. Durante i bombardamenti, noi sentivamo i bombardieri solo quando ce li avevamo ormai sopra la testa. Restava il tempo solo di tirare fuori il rosario. Se qualcuno ha parenti anziani a Bolzano si faccia raccontare un po' com'era vivere in una città bombardata con tanta insistenza.
    Col radar si sarebbero potuti avvistare con quel pochino di tempo in più per mettere in salvo qualcuno. Chiaramente se si sfolla bisogna sfollare bene, altrimenti è meglio lasciare la gente crepare sotto le bombe perché sennò si creerebbe disordine eccetera (), ma forse in certi contesti si sarebbe salvata un sacco di gente.

  4. #124
    Uragano L'avatar di Conte
    Data Registrazione
    06/03/04
    Località
    Grottammare (AP) 0-170 mt slm
    Età
    63
    Messaggi
    17,250
    Menzionato
    10 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da belli83 Visualizza Messaggio
    Dai, non prendiamoci in giro..
    O l'evacuazione la fai bene (e mandi via tutti), comprese case di riposo, ospedali, carcere, anziani, ecc.. o non la fai.
    Te la immagini la PC che dice: andate fuori casa e dormite in auto nei parcheggi?
    E' quello che sta facendo la popolazione di Ascoli e Teramo ultimamente tanto per dirne una visto che si sono resi conto che i tanto conclamati esperti non c'hanno capito purtroppo una mazza affermando che la situazione era tranquilla (ma non era vero che i terremoti non si possono prevedere? Questa non è una previsone? Una brutta previsione direi). Quindi per te è meglio morire nel letto di casa che dormire in auto per qualche notte.

    Mi sembra che anche tu ti stia aggrappando ai vetri pur di screditare Giuliani senza nemmeno aver letto o visto qualcosa detta e relazionata di suo pugno.
    Se sai che esiste un metodo che ha buone percentuali di successo previsionale nel breve (2-3 ore), anche se non hai il 100% di affidabilità, predisponi opportunamente un piano di evacuazione. Giuliani avrebbe voluto essere sostenuto nelle sue osservazioni e ricerche per collaborare (esistono le prove) con la Protezione Civile e quant'altro. Non sarebbe spettato a lui proporre un piano di emergenza ma solo mettere a disposizione le sue consocenze per giungere ad un miglioramento della salvaguardia della popolazione a rischio.

    Quindi se vuoi capire bene altrimenti rischi di crearti un bel muro mentale come tanti eminenti scienziati che abbiamo visto in TV.
    Sono veramente stanco di sparate a zero su Giuliani da parte di chi non sa nulla se non per sentito dire qua e là.
    Giuliani non mi è parente e manco lo conosco ma basta indagare un pò (in rete si trovano alcuni interessanti documenti a parte le sue interviste) per capire che non è l'ultimo dei cretini, che i suoi intenti sono degni di attenzione e che alcuni esperti sismologi "sparano" in malafede per timore di sentirsi defraudati del lustro goduto fino al 5 aprile.
    Un bagno di umiltà e buon senso non guasterebbe.

    La cosa curiosa è che noi partecipiamo ad un forum di ricercatori appassionati dove in maniera evidente eccellono dei bravissimi personaggi ormai noti, persone che magari nella vita non sono meteorologi di professione ma che hanno doti inequivocabilmente riconosciute. Se Giuliani dovesse essere zittito, e con lui anche i suoi studi e rilevamenti, dovremmo anche noi per logica spegnere le ns stazioni e chiudere il forum. Non credi?
    Ultima modifica di Conte; 10/04/2009 alle 16:07
    Davis Vantage Pro2 wireless. LaCrosse WS 2300 cablata con schermo artigianale autoventilabile.
    "L'uomo non ha avuto il mondo in regalo dai suoi genitori, lo ha avuto in prestito dai propri figli"


  5. #125
    Brezza tesa
    Data Registrazione
    19/01/07
    Località
    Fonte Nuova (RM NordEst 150 s.l.m)
    Messaggi
    862
    Menzionato
    1 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Dopo giorni di esperimenti, ricerche e studi, sono riuscito a scrivere un programmino python che grafica i dati dell'INGV ( http://cnt.rm.ingv.it/~earthquake/index_web_cnt.php ) mantenendo la proporzione dei tempi (impossibile farlo in Excel per esempio).

    Quando ho visto per la prima volta il grafico nudo e crudo, mi è preso un colpo.
    Quando ci ho aggiunto il grafico del radon di Giuliani, mi è cascata la mascella.

    Riporto qui un'immagine preliminare con indicazioni approssimative della scala temporale, fatta misurando i pixel... Adesso vedo se riesco a spiegare allo script in python come disegnare qualche linea verticale di riferimento.



    Questo è lo script (richiede WxPython):
    Codice:
    import wx
    
    FACTORX=20 # Fattore di scala orizzontale per il grafico dell'intensita
    FACTORY=20 # Fattore di scala verticale per il grafico dell'intensita
    FACTORY2=10 # Fattore di scala verticale per il grafico della profondita
    
     
    class MyFrame(wx.Frame):
       def __init__(self):
     
          self.frame1 = wx.Frame(None, title="Studio preliminare terremoto L'Aquila 2009", id=-1, size=(820, 620))
          self.panel1 = wx.Panel(self.frame1)
      
          self.statbmp = wx.StaticBitmap(self.panel1, id=-1, pos=(0,0), size=(640,150))
     
          self.image = wx.Bitmap("photo.bmp") # E' necessaria nella cartella del programma un'immagine bianca di 800x600 pixel
          canvas_dc = wx.MemoryDC(self.image)
          canvas_dc.DrawBitmap(self.image, 0, 0)
    
          file = open('c://temp//terr.txt', 'r') # Dati dei terremoti presi da http://cnt.rm.ingv.it/~earthquake/index_web_cnt.php
          # Inzializza alcune variabili
          number=0
          MAX=200
          x=0
          y=0
          x2=0
          y2=0
          Depth=0
          
          # Inizia ciclo lettura/disegno
          Line=file.readline()
          while Line<>"" and number<MAX:
            number=number+1
            Date=Line[13:32]
            Year=Date[0:4]
            Month=Date[5:7]
            Day=Date[8:10]
            Hour=Date[11:13]
            Min=Date[11:12]
            Total=int(Day)*1440+int(Hour)*60+int(Min) # Converte la data in numero (non considera anno e mese)
            error=0
            try: 
              Mag = float(AfterTab(Line,6,4,7))   # Magnitudo
              Depth = float(AfterTab(Line,5,1,4)) # Profondita
              Lat = float(AfterTab(Line,3,1,7)) # Latitudine
              Lon = float(AfterTab(Line,4,1,7)) # Longitudine
            except Exception,e:
              Mag=0
              Lat=0
              Lon=0
              error=1
              print Line, "ERRORE: ", str(e)
            #print Day,Hour,Min,Total,Mag,Lat,Lon
            
            # Memorizza coordinate precedenti per permettere tracciamento linea tra due punti consecutivi:
            prevx=x 
            prevy=y
            prevx2=x2
            prevy2=y2
            
            # Prepara coordinate da disegnare:
            x=Total/FACTORX
            y=200-Mag*FACTORY
            x2=x
            y2=600-Depth*FACTORY2
            print x,y, "   ",
            
            if (prevx<>0 and prevy<>0):
              if error==0:
                canvas_dc.DrawLine(x1=prevx, y1=prevy, x2=x, y2=y)
                canvas_dc.DrawLine(x1=prevx2, y1=prevy2, x2=x2, y2=y2)
            else: # First line
              if error==0:
                canvas_dc.DrawLine(x1=x, y1=y, x2=x+1, y2=y+1)
                canvas_dc.DrawLine(x1=x2, y1=y2, x2=x2+1, y2=y2+1)
                
            # Disegna rettangolo su ogni vertice            
            rect = wx.Rect(x, y, 4, 4) 
            canvas_dc.DrawRoundedRectangleRect(rect, 8)        
            rect = wx.Rect(x2, y2, 4, 4) 
            canvas_dc.DrawRoundedRectangleRect(rect, 8)        
            Line = file.readline() # Legge nuova riga
          file.close()
    
          self.statbmp.SetBitmap(self.image) 
          self.frame1.Show(True)
     
    
    def AfterTab(s,n,inizio,fine):
    # Legge dalla riga S il valore presente dopo il tab numero N (da carattere INIZIO a carattere FINE)
      count = 0
      index = s.find(chr(9))  # use default start
      while index != -1 and count<n:
          count += 1
          start = index + 1
          index = s.find(chr(9), start)  # set new start
      return s[index+inizio:index+fine]
    
      
    app = wx.PySimpleApp()
    f = MyFrame()
    app.MainLoop()
    NOTA:
    E' necessaria nella cartella del programma un'immagine bianca di 800x600 pixel chiamata "photo.bmp" (è il mio primo programma grafico in python, non sono molto pratico...)

    Questi sono i dati (non mi pare di averli riformattati dopo averli copiati dal sito, ma nn si sa mai...)
    Codice:
    2206563060	 	2009/04/10	15:46:17	42.349	13.376	9.9	Ml:3.3	 	Aquilano
    2206562820	 	2009/04/10	15:22:43	42.248	13.484	10.2	Ml:3	 	Velino-Sirente
    2206562670	 	2009/04/10	15:07:20	42.359	13.311	8.4	Ml:2.9	 	Aquilano
    2206560730	 	2009/04/10	11:53:09	42.243	13.484	9.6	Ml:3.1	 	Velino-Sirente
    2206557610	 	2009/04/10	06:41:31	42.52	13.339	9.3	Ml:3	 	Monti_della_Laga
    2206556330	 	2009/04/10	04:33:04	42.458	13.35	10.5	Ml:3.3	 	Gran_Sasso
    2206555610	 	2009/04/10	03:22:22	42.47	13.417	9.4	Ml:3.7	 	Gran_Sasso
    2206552790	 	2009/04/09	22:40:06	42.481	13.298	10.9	Ml:3.6	 	Gran_Sasso
    2206551670	 	2009/04/09	20:47:01	42.493	13.305	10	Ml:3.1	 	Gran_Sasso
    2206550980	 	2009/04/09	19:38:16	42.501	13.356	17.2	Ml:4.9	 	Gran_Sasso
    2206550090	 	2009/04/09	18:09:52	44.258	11.946	29.6	Ml:2.6	 	Zona_Forli
    2206548380	 	2009/04/09	15:18:14	42.308	13.495	10	Ml:3.2	 	Valle_dell'Aterno
    2206547790	 	2009/04/09	14:19:20	42.45	13.454	9.3	Ml:2.5	 	Gran_Sasso
    2206547180	 	2009/04/09	13:19:33	42.338	13.259	10.0	Ml:3.6	 	Aquilano
    2206543700	 	2009/04/09	07:31:44	44.247	11.95	27.7	Ml:2.7	 	Zona_Forli
    2206542030	 	2009/04/09	04:43:09	42.506	13.366	9.2	Ml:3.7	 	Gran_Sasso
    2206541910	 	2009/04/09	04:32:44	42.445	13.42	8.1	Ml:4	 	Gran_Sasso
    2206541410	 	2009/04/09	03:41:55	42.51	13.332	9.8	Ml:3.1	 	Monti_della_Laga
    2206541141	 	2009/04/09	03:14:52	42.338	13.437	18.0	Ml:4.2	 	Aquilano
    2206539720	 	2009/04/09	00:52:59	42.484	13.343	15.4	Ml:5.1	 	Gran_Sasso
    2206538780	 	2009/04/08	23:18:06	42.391	13.325	10.8	Ml:3.3	 	Aquilano
    2206538560	 	2009/04/08	22:56:50	42.507	13.364	10.2	Ml:4.3	 	Gran_Sasso
    2206538190	 	2009/04/08	22:19:56	42.364	13.395	9.4	Ml:2.7	 	Aquilano
    2206535580	 	2009/04/08	17:58:35	42.364	13.396	8.8	Ml:3.2	 	Aquilano
    1206531750	 	2009/04/08	11:35:57	42.355	13.328	10.6	Ml:3.5	 	Aquilano
    1206531130	 	2009/04/08	10:34:09	42.352	13.381	9.2	Ml:3.1	 	Aquilano
    2206528390	 	2009/04/08	06:00:24	42.38	13.36	9.9	Ml:2.1	 	Aquilano
    2206527470	 	2009/04/08	04:27:41	42.305	13.467	10.7	Ml:3.7	 	Valle_dell'Aterno
    2206526590	 	2009/04/08	03:00:34	42.299	13.459	10	Ml:3.5	 	Aquilano
    1206524030	 	2009/04/07	22:43:42	42.328	13.486	9.3	Ml:3	 	Valle_dell'Aterno
    1206523891	 	2009/04/07	22:29:50	42.284	13.489	9.6	Ml:3.1	 	Valle_dell'Aterno
    2206523390	 	2009/04/07	21:39:06	42.361	13.363	10.4	Ml:3.7	 	Aquilano
    1206523330	 	2009/04/07	21:34:29	42.38	13.376	7.4	Ml:4.2	 	Aquilano
    2206523150	 	2009/04/07	21:15:13	44.159	11.974	10	Ml:2.5	 	Appennino_forlivese
    2206522630	 	2009/04/07	20:24:54	39.177	16.795	13.8	Ml:3.2	 	La_Sila
    1206521070	 	2009/04/07	17:47:37	42.275	13.464	15.1	Ml:5.3	 	Valle_dell'Aterno
    1206517890	 	2009/04/07	12:29:00	42.461	13.412	8.3	Ml:3.3	 	Gran_Sasso
    2206516690	 	2009/04/07	10:29:10	42.321	13.4	9.5	Ml:3.2	 	Aquilano
    8206516100	 	2009/04/07	09:30:56	42.334	13.355	10.0	Ml:3.5	 	Aquilano
    2206516040	 	2009/04/07	09:26:28	42.342	13.388	10.2	Ml:4.7	 	Aquilano
    2206515190	 	2009/04/07	08:00:52	42.383	13.401	15.1	Ml:3.2	 	Aquilano
    2206514140	 	2009/04/07	06:15:07	42.378	13.32	11.2	Ml:3	 	Aquilano
    2206509940	 	2009/04/06	23:15:37	42.451	13.364	8.6	Ml:4.8	 	Gran_Sasso
    1206509670	 	2009/04/06	22:47:13	42.349	13.293	11.6	Ml:3.6	 	Aquilano
    1206509160	 	2009/04/06	21:56:53	42.396	13.323	9.7	Ml:3.8	 	Aquilano
    2206505980	 	2009/04/06	16:38:09	42.362	13.333	10.2	Ml:4	 	Aquilano
    2206502360	 	2009/04/06	10:36:18	42.343	13.402	10.1	Ml:3.5	 	Aquilano
    2206501990	 	2009/04/06	09:59:29	42.322	13.381	10	Ml:3.1	 	Aquilano
    1206500360	 	2009/04/06	07:17:10	42.355	13.367	9.2	Ml:3.9	 	Aquilano
    1206499150	 	2009/04/06	05:15:28	42.482	13.358	4.2	Ml:3.1	 	Gran_Sasso
    1206498870	 	2009/04/06	04:47:53	42.352	13.347	9.4	Ml:3.7	 	Aquilano
    2206498420	 	2009/04/06	04:02:27	42.322	13.399	10.8	Ml:2.7	 	Aquilano
    1206498100	 	2009/04/06	03:30:59	42.335	13.388	9	Ml:2.8	 	Aquilano
    1206498020	 	2009/04/06	03:23:10	42.333	13.331	9	Ml:2.8	 	Aquilano
    1206497960	 	2009/04/06	03:16:04	42.315	13.377	9.8	Ml:2.8	 	Aquilano
    1206497750	 	2009/04/06	02:55:33	42.366	13.328	9.2	Ml:3.3	 	Aquilano
    2206497640	 	2009/04/06	02:44:32	42.302	13.354	10.0	Ml:3.2	 	Aquilano
    2206497570	 	2009/04/06	02:37:04	42.366	13.34	10.1	Ml:4.6	 	Aquilano
    2206497510	 	2009/04/06	02:31:33	42.405	13.394	8.2	Ml:3.4	 	Aquilano
    2206497470	 	2009/04/06	02:27:46	42.374	13.342	10.0	Ml:3.9	 	Aquilano
    2206497450	 	2009/04/06	02:25:45	42.427	13.357	13.3	Ml:3.3	 	Aquilano
    2206496920	 	2009/04/06	01:32:39	42.334	13.334	8.8	Ml:5.8	 	Aquilano
    2206495180	 	2009/04/05	22:39:41	42.341	13.38	8.5	Ml:3.5	 	Aquilano
    1206494510	 	2009/04/05	21:31:23	44.161	11.932	10	Ml:2.5	 	Appennino_forlivese
    2206494080	 	2009/04/05	20:48:54	42.332	13.372	8.4	Ml:3.9	 	Aquilano
    2206493790	 	2009/04/05	20:20:53	44.236	11.999	28.2	Ml:4.6	 	Zona_Forli
    2206475490	 	2009/04/04	13:50:41	38.561	14.774	19.3	Ml:2.6	 	Isole_Lipari
    2206456830	 	2009/04/03	06:43:24	42.321	13.361	10.3	Ml:2.6	 	Aquilano
    2206455640	 	2009/04/03	04:44:42	42.327	13.358	9.4	Ml:3	 	Aquilano
    2206445110	 	2009/04/02	11:11:54	42.307	13.382	10.1	Ml:2.4	 	Aquilano
    2206435450	 	2009/04/01	19:05:11	43.589	10.433	9.7	Ml:1.8	 	Costa_toscana_settentrionale
    2206435130	 	2009/04/01	18:33:58	42.333	13.368	10.4	Ml:2.0	 	Aquilano
    2206433150	 	2009/04/01	15:15:38	42.355	13.378	10.6	Ml:2.2	 	Aquilano
    2206430230	 	2009/04/01	10:23:28	42.324	13.361	9.9	Ml:2.5	 	Aquilano
    2206425040	 	2009/04/01	01:44:46	42.336	13.38	10.1	Ml:1.8	 	Aquilano

  6. #126
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Cacchio non hai messo l'intestazione delle colonne.

    ???, data, ora, e le ultime tre? E le successive? Vorrei provare a importare quei dati in un programma statistico e correlarli in qualche modo. Oimé la regressione logistica per dati seriali non la so fare.

    Non so Giuliani, ma penso che l'INGV dia l'orario UTC.

  7. #127
    Brezza tesa
    Data Registrazione
    19/01/07
    Località
    Fonte Nuova (RM NordEst 150 s.l.m)
    Messaggi
    862
    Menzionato
    1 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da Borat Visualizza Messaggio
    Cacchio non hai messo l'intestazione delle colonne.

    ???, data, ora, e le ultime tre? E le successive? Vorrei provare a importare quei dati in un programma statistico e correlarli in qualche modo. Oimé la regressione logistica per dati seriali non la so fare.

    Non so Giuliani, ma penso che l'INGV dia l'orario UTC.
    Event-Id Data Ora (UTC) Lat Lon Prof(Km) Mag Distretto Sismico

    Noto pero' ora che le scosse "incriminate" hanno epicentro nelle isole Lipari e e a Forli'!

    Quindi mi sa che è solo una coincidenza....
    Pero' è una coincidenza VERAMENTE bislacca!!!!

    Ora sto provando a tracciare un grafico delle coordinate, per vedere come so è spostato l'epicentro nel tempo, ma non so se mi riesce.

    Sarebbe anche interessante conoscere le DURATE di ogni scossa (non capisco perche' non le inseriscono nelle varie liste, come se un terremoto di un secondo o un minuto fossero la stessa cosa!!!)

  8. #128
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Citazione Originariamente Scritto da jumpjack Visualizza Messaggio
    Event-Id Data Ora (UTC) Lat Lon Prof(Km) Mag Distretto Sismico

    Noto pero' ora che le scosse "incriminate" hanno epicentro nelle isole Lipari e e a Forli'!

    Quindi mi sa che è solo una coincidenza....
    Pero' è una coincidenza VERAMENTE bislacca!!!!

    Ora sto provando a tracciare un grafico delle coordinate, per vedere come so è spostato l'epicentro nel tempo, ma non so se mi riesce.

    Sarebbe anche interessante conoscere le DURATE di ogni scossa (non capisco perche' non le inseriscono nelle varie liste, come se un terremoto di un secondo o un minuto fossero la stessa cosa!!!)
    Fra un minuto ti posto i dati dei soli terremonti nell'Aquilano (cioè Aquilano, Aterno, Velino-Sirente e Gran Sasso, togliendo Forlì, Sila, Toscana eccetera).

    Dammi 5 minuti per elaborarli.

  9. #129
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Ficata, sono riuscito a importare in un programma per l'analisi dei dati quantitativi i dati INGV di aprile.
    Qui vediamo la serie dei terremoti del mese, etichettati ciascuno col luogo dove si è verificato. Per migliorare il look bisognerebbe aggiungere un record con magnitudo 0 un minuto dopo ogni record con magnitudo non-zero (così il grafico scende e non fa quelle planate o salite assurde).




    Qui vediamo il numero di eventi per luogo:
    Codice:
    -----------------------------------
    area                              n
    -----------------------------------
    Aquilano                         46
    Gran_Sasso                       13
    Valle_dell'Aterno                 5
    Zona_Forli                        3
    Monti_della_Laga                  2
    Velino-Sirente                    2
    Appennino_forlivese               2
    Costa_toscana_settentrionale      1
    Isole_Lipari                      1
    La_Sila                           1
    -----------------------------------
    somma                            76
    -----------------------------------
    e qui la sommatoria delle magnitudo e le loro medie per zona:
    Codice:
    ------------------------------------------
    area                          somma  media
    ------------------------------------------
    Aquilano                      152.1    3.3
    Gran_Sasso                     49.4    3.8
    Valle_dell'Aterno              18.3    3.7
    Zona_Forli                      9.9    3.3
    Monti_della_Laga                6.1    3.0
    Velino-Sirente                  6.1    3.0
    Appennino_forlivese             5.0    2.5
    Isole_Lipari                    2.6    2.6
    La_Sila                         3.2    3.2
    Costa_toscana_settentrionale    1.8    1.8
    ------------------------------------------
    totale                        254.5    3.3
    ------------------------------------------
    Abbiamo visto un bell'esempio del fatto che generalmente non bisogna farsi fregare dalle medie. Se badiamo alla media, l'aquilano pare che abbia avuto una sismicità nella media, quanto La Sila o la zona di Forlì. Solo che Forlì e Sila sono OK, e L'Aquila sembra che sia stata bombardata.
    A ogni modo, mi hanno detto che la scala Richter è logaritmica e non so molto bene quanto corretto sia calcolare banalmente la media dei logaritmi e che cosa la media dei logarimti significhi rispetto alle grandezze da cui i logaritmi sono tratti. Dovrei rinfrescarmi un po' in matematica.

    Adesso però vorrei sovrapporre alla line dei terremoti quella del radon di Giuliani. Dove trovo i dati in formato tabulare/listato?

  10. #130
    Burrasca L'avatar di Borat
    Data Registrazione
    08/06/08
    Località
    -999
    Messaggi
    5,669
    Menzionato
    0 Post(s)

    Predefinito Re: Previsione terremoti non scientificamente provata...(qui le notizie, non le polemiche)

    Ecco i dati con solo siti aquilani, cioè Aquilano, Gran Sasso, Laga, Aterno e Velino-Sirente.
    Basta che sostituisci in bloc notes le virgole con un tab per ricavare un tab-delimited (azz, quando parlo di queste cose anni Ottanta mi sento veramente un matusa) e puoi incollare direttamente in excel
    Codice:
    data,ora,latit,longit,profond,magnit,dove
    2009/04/01,01:44:46,42.336,13.38,10.1,1.8,Aquilano
    2009/04/01,10:23:28,42.324,13.361,9.9,2.5,Aquilano
    2009/04/01,15:15:38,42.355,13.378,10.6,2.2,Aquilano
    2009/04/01,18:33:58,42.333,13.368,10.4,2,Aquilano
    2009/04/02,11:11:54,42.307,13.382,10.1,2.4,Aquilano
    2009/04/03,04:44:42,42.327,13.358,9.4,3,Aquilano
    2009/04/03,06:43:24,42.321,13.361,10.3,2.6,Aquilano
    2009/04/05,20:48:54,42.332,13.372,8.4,3.9,Aquilano
    2009/04/05,22:39:41,42.341,13.38,8.5,3.5,Aquilano
    2009/04/06,01:32:39,42.334,13.334,8.8,5.8,Aquilano
    2009/04/06,02:25:45,42.427,13.357,13.3,3.3,Aquilano
    2009/04/06,02:27:46,42.374,13.342,10,3.9,Aquilano
    2009/04/06,02:31:33,42.405,13.394,8.2,3.4,Aquilano
    2009/04/06,02:37:04,42.366,13.34,10.1,4.6,Aquilano
    2009/04/06,02:44:32,42.302,13.354,10,3.2,Aquilano
    2009/04/06,02:55:33,42.366,13.328,9.2,3.3,Aquilano
    2009/04/06,03:16:04,42.315,13.377,9.8,2.8,Aquilano
    2009/04/06,03:23:10,42.333,13.331,9,2.8,Aquilano
    2009/04/06,03:30:59,42.335,13.388,9,2.8,Aquilano
    2009/04/06,04:02:27,42.322,13.399,10.8,2.7,Aquilano
    2009/04/06,04:47:53,42.352,13.347,9.4,3.7,Aquilano
    2009/04/06,05:15:28,42.482,13.358,4.2,3.1,Gran_Sasso
    2009/04/06,07:17:10,42.355,13.367,9.2,3.9,Aquilano
    2009/04/06,09:59:29,42.322,13.381,10,3.1,Aquilano
    2009/04/06,10:36:18,42.343,13.402,10.1,3.5,Aquilano
    2009/04/06,16:38:09,42.362,13.333,10.2,4,Aquilano
    2009/04/06,21:56:53,42.396,13.323,9.7,3.8,Aquilano
    2009/04/06,22:47:13,42.349,13.293,11.6,3.6,Aquilano
    2009/04/06,23:15:37,42.451,13.364,8.6,4.8,Gran_Sasso
    2009/04/07,06:15:07,42.378,13.32,11.2,3,Aquilano
    2009/04/07,08:00:52,42.383,13.401,15.1,3.2,Aquilano
    2009/04/07,09:26:28,42.342,13.388,10.2,4.7,Aquilano
    2009/04/07,09:30:56,42.334,13.355,10,3.5,Aquilano
    2009/04/07,10:29:10,42.321,13.4,9.5,3.2,Aquilano
    2009/04/07,12:29:00,42.461,13.412,8.3,3.3,Gran_Sasso
    2009/04/07,17:47:37,42.275,13.464,15.1,5.3,Valle_dell'Aterno
    2009/04/07,21:34:29,42.38,13.376,7.4,4.2,Aquilano
    2009/04/07,21:39:06,42.361,13.363,10.4,3.7,Aquilano
    2009/04/07,22:29:50,42.284,13.489,9.6,3.1,Valle_dell'Aterno
    2009/04/07,22:43:42,42.328,13.486,9.3,3,Valle_dell'Aterno
    2009/04/08,03:00:34,42.299,13.459,10,3.5,Aquilano
    2009/04/08,04:27:41,42.305,13.467,10.7,3.7,Valle_dell'Aterno
    2009/04/08,06:00:24,42.38,13.36,9.9,2.1,Aquilano
    2009/04/08,10:34:09,42.352,13.381,9.2,3.1,Aquilano
    2009/04/08,11:35:57,42.355,13.328,10.6,3.5,Aquilano
    2009/04/08,17:58:35,42.364,13.396,8.8,3.2,Aquilano
    2009/04/08,22:19:56,42.364,13.395,9.4,2.7,Aquilano
    2009/04/08,22:56:50,42.507,13.364,10.2,4.3,Gran_Sasso
    2009/04/08,23:18:06,42.391,13.325,10.8,3.3,Aquilano
    2009/04/09,00:52:59,42.484,13.343,15.4,5.1,Gran_Sasso
    2009/04/09,03:14:52,42.338,13.437,18,4.2,Aquilano
    2009/04/09,03:41:55,42.51,13.332,9.8,3.1,Monti_della_Laga
    2009/04/09,04:32:44,42.445,13.42,8.1,4,Gran_Sasso
    2009/04/09,04:43:09,42.506,13.366,9.2,3.7,Gran_Sasso
    2009/04/09,13:19:33,42.338,13.259,10,3.6,Aquilano
    2009/04/09,14:19:20,42.45,13.454,9.3,2.5,Gran_Sasso
    2009/04/09,15:18:14,42.308,13.495,10,3.2,Valle_dell'Aterno
    2009/04/09,19:38:16,42.501,13.356,17.2,4.9,Gran_Sasso
    2009/04/09,20:47:01,42.493,13.305,10,3.1,Gran_Sasso
    2009/04/09,22:40:06,42.481,13.298,10.9,3.6,Gran_Sasso
    2009/04/10,03:22:22,42.47,13.417,9.4,3.7,Gran_Sasso
    2009/04/10,04:33:04,42.458,13.35,10.5,3.3,Gran_Sasso
    2009/04/10,06:41:31,42.52,13.339,9.3,3,Monti_della_Laga
    2009/04/10,11:53:09,42.243,13.484,9.6,3.1,Velino-Sirente
    2009/04/10,15:07:20,42.359,13.311,8.4,2.9,Aquilano
    2009/04/10,15:22:43,42.248,13.484,10.2,3,Velino-Sirente
    2009/04/10,15:46:17,42.349,13.376,9.9,3.3,Aquilano
    2009/04/10,19:07:21,42.375,13.393,9.5,3.1,Aquilano
    68 terremonti misurabili in 10 giorni.

    Adesso, se avessimo nello stesso formato i dati di Giuliani sarebbe molto semplice sovrapporre le linee e vedere che cosa ne esce. Ma suppongo che l'INGV abbia già mandato una squadra di ninja a richiederli.

Segnalibri

Permessi di Scrittura

  • Tu non puoi inviare nuove discussioni
  • Tu non puoi inviare risposte
  • Tu non puoi inviare allegati
  • Tu non puoi modificare i tuoi messaggi
  •