Sta per nascere la pagina web per che permetterà di stimare l'epicentro dei terremoti. Occorre selezionare la stazione e poi mettere il momento dell'arrivo dell'onda P (si vede nel Waterfall ne avevamo già parlato).
Si inviano i dati al server che risponderà con una visualizzazione di una cartografia georeferenziata con delle iperboli...poi passeremo ala visualizzazione di un cerchio...il più piccolo possibile con l'elenco delle città e la distanza dal punto.

Non tutti gli utenti potranno fare questa operazione, anche perchè gli allarmi arrivano mediante telegram, e solo chi è registrato riceve gli allarmi in tempo reale e può vedere le pagine web mentre il Waterfall viene costruito.

Questo è il metodo manuale che ci permetterà di testare il metodo automatico...cioè se il metodo automatico commette qualche errore...mediante questa pagina potremmo effettuare delle correzioni.

Penso di essere stato chiaro...ma se non lo sono stato ...fate un fischio!

Questo è un esempio:

Questo è un esempio di un messaggio sul calcolo dell'epicentro.Arriva al cellulare sulla chat TERREMOTIForumEA di Telegram :

Passiamo dalle iperboli a localizzare la zona dell'epicentro in modo automatico.Vengono calcolate le intersezioni delle iperboli e poi il baricentro delle intersezioni...e a quel punto viene fatto un cerchio della zona di interesse.

Questi sono i dati...tutto si poggia sui tempi di arrivo delle onde P. Con 3 orari vengono generate 3 iperboli:

Come è andato il primo calcolo automatico dell'epicentro??
Due iperboli sono state identificate per bene...la terza ha preso il ramo sbagliato
L'epicentro è stato sbagliato di 17 km...non molto considerando l'errore che ha combinato

Ho inserito anche i tempi di Fabriano...esce fuori un casino...ma nonostante i rami sbagliati e i calcoli che non dovrebbero riguardare iperboli che non si incrociano...l'epicentro calcolato non è molto lontano...perchè spostare il baricentro dei punti non è facile!

CITAZIONE (Bolle, 07/06/2018 15:32:51 )

Non ci avete capito nulla...state sereni...è tutto sotto controllo!

Questo il sisma di Amatrice all'ora di pranzo con il Waterfall di alcune postazioni:

Ci siamoooooooooooo...ho cambiato l'implementazione...ho stravolto totalmente approccio semplificandolo concettualmente!

La spiego al volo...ricordate come viene calcolata l'iperbole? Si hanno due tempi di arrivo alle postazioni e si calcolano tutte le possibili soluzioni ...da questa elaborazione esce fuori la curva dove unendo tutti i segmenti appare l'iperbole. Per calcolare l'intersezione delle due iperboli si passa a calcolare l'intersezione di ogni segmento di una iperbole per tutti gli altri dell'altra iperbole. Il processo di elaborazione verifica tutte le intersezioni delle iperboli a due a due. In questo modo otteniamo tutti i punti in questione e nell'immagine viene visualizzato un cerchio giallo. Aumenta sicuramente il tempo di elaborazione, che ovviamente dura qualche istante, ma questo metodo rispetto il precedente, semplifica la gestione che era diventata troppo complessa...in sintesi si è scomposto un problema di difficile soluzione in sotto-problemi semplici da risolvere! Non la vedere così? Vabbè...ci ho provato!


Questi dei tempi messi a casaccio per testare l'algoritmo!

Questa notte...un fulmine deve aver cotto la scheda amplificatrice...la postazione della sabina è OFF!!!

Piano piano riesco ad inquadrare il problema...chiano chiano.
Dopo un po di consulti con il prof la mia zucca ha partorito questo pseudo algoritmo per calcolare l'epicentro...suddiviso in almeno DIECI FASI.

1. calcolo iperboli a velocità media dell'onda P delle stazioni
2. Calcolo intersezione di tutte le iperboli
3. Calcolo baricentro delle intersezioni
4. Sfrondamento intersezioni dal baricentro ottenuto
5. Calcolo baricentro intersezioni dell'epicentro provvisorio
6. Rielaborazione iperboli ipotizzando l'epicentro provvisorio con velocità dell'onda P progressiva
7. Ricalcolo intersezioni delle iperboli
8. Ricalcolo baricentro delle intersezioni
9. Sfrondamento intersezioni
10. Calcolo Epicentro definitivo

Si potrebbero ripetere le fasi da 6 a 10 se c'è uno spostamento importante dell'epicentro!
Tutti i calcoli , anche se possono sembrare molti e complessi vengono effettuati in meno di un secondo

Facciamo un esempio relativo al 28 maggio 2018 al sisma avvenuto a Pizzoli (AQ) di cui potete trovare i dati del WF nei messaggi precendenti. Cmq c'è qualcosa che non mi convince... si può ancora migliorare.

Ripartito la postazione SABINA...sono saltati entrambi gli ADC per il fulmine delle 2:44 di notte !

Le sorprese non finiscono mai, ho inserito i dati del terremoto 1.9 avvenuto poco fa ad Arquata del Tronto...il sistema con le tre onde p prese dalla stazione di Fabriano,Sabina e Viterbo ha una indeterminazione!!!!

NUOVO ALGORITMO

Ipotesi terremoto vicino

a)si tassella la superficie utile ad es con quadrati di 20km di lato

b)si simula il terremoto con v standard si calcola lo scarto quadratico medio tra i tempi simulati e i tempi effettivamente rilevati

c) il tassello che minimizza lo scarto quadratico medio contiene l'epicentro.

d) si restrige l'area al tassello (e a quelli vicini) con un nuovo tassello ad es di 500m di lato

e) si ripete il calcolo con velocita prese da un data base e si trova la zona che minimizza lo scato quadratico medio
si può raffinare ulteriormente la ricerca per trovare il punto più probabile