La Provincia magmatica Romana
Aprile 10, 2019
La carta di Casagrande: come costruirla con python
Maggio 15, 2019
Show all

Geologia e Python: uno script per le prove penetrometriche


Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 118

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 119

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 118

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 119

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 118

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 119

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 118

Notice: A non well formed numeric value encountered in /web/htdocs/www.intrageo.it/home/wp-content/plugins/crayon-syntax-highlighter/crayon_formatter.class.php on line 119

In un post precedente abbiamo parlato di prove penetrometriche con le relative modalità esecutive.

Inoltre abbiamo visto come sia molto utile l’impiego di Python nelle geoscienze. In questo post vediamo come elaborare e creare il grafico di una prova penetrometrica media – dopo la fase di lavoro in campagna – usando Python, nello specifico Pyplot e Pandas.

Prima di iniziare a scrivere il codice dobbiamo preparare un file csv da usare durante l’elaborazione. Dentro un foglio Excell – anche se io uso LibreOffice Calc, un foglio elettronico libero, componente del software di produttività personale LibreOffice – inseriamo le profondità raggiunte con i relativi numeri di colpi. Una volta completato il file lo esportiamo in formato csv.

python-geologia

Il codice Python per il grafico delle prove penetrometriche

Innanzitutto prepariamo l’ambiente di lavoro importando le librerie necessarie.

Adesso importiamo il file csv che contiene i dati, usando pandas, e creiamo due liste: una che contiene le profondità ed una che contiene il numero di colpi.

Inseriamo nello script le caratteristiche del penetrometro usato per questa prova penetrometrica (una DPM). Inoltre con un ciclo for compiliamo la lista Rd con i valori di resistenza dinamica alla punta, calcolata con la formula degli olandesi:

 \centering Rd=\frac{M^2*H}{(e*A)*(M*Q)}

Il codice Python per questa parte è:

Procediamo con il codice in modo da creare i grafici delle prove penetrometriche: posizioniamo a sinistra il grafico profondità-numero di colpi mentre a destra il grafico profondità-resistenza alla punta. Usiamo plt.subplot() per affiancare i due grafici.

Il grafico risultante é

grafico-prova-penetrometrica
Grafici della prova penetrometrica generati con Python

Approfondimenti

Per approfondire l’argomento trattatato in questo post consiglio i seguenti testi:

Antonio Nirta
Antonio Nirta
Ciao, mi chiamo Antonio Nirta e sono un geologo. Mi sono laureato all'Università di Pisa e dal 2017 svolgo la libera professione. Attraverso i post presenti in questo blog cercherò di trasmettere la mia passione per le Scienze Geologiche e di fornire informazioni che riguardano il bellissimo settore della Geologia.

Lascia un commento

Il tuo indirizzo email non sarà pubblicato. I campi obbligatori sono contrassegnati *