{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"collapsed": true
},
"source": [
"# Lectura de archivo, calculo de percentiles y graficado de histograma\n",
"\n",
"## Integrantes:\n",
"\n",
"* Vanesa Pantano\n",
"* Alejandro Godoy\n",
"\n",
"## DescripciĆ³n:\n",
"\n",
"Vamos a leer los datos de precipitacion diarios de estaciones de Argentina, y nos vamos a quedar con el mes seleccionado y realizaremos un histograma de la estacion Azul\n",
"\n",
"[Este](https://github.com/pewen/wpc/tree/master/Trabajos_finales/Pantano_Godoy/Precipitacion.txt) es el archivo que analizamos"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"#Leo Archivo\n",
"import numpy as np\n",
"import pandas as pd\n",
"fname='../Precipitacion.txt'\n",
"#f=open('../Precipitacion_2.txt','r')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#leo con np fromtxt\n",
"#datos=np.genfromtxt(fname,dtype=None,skip_header=1,usecols = (4))"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#Leo con Pandas\n",
"#data = pd.read_csv(fname, skiprows=1, header=None, sep='\\t')\n",
"data = pd.read_csv(fname, header=0, index_col=0, sep='\\t', dtype=float)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" anio mes dia Azul BahiaBlanca 'Bella Vista' 'Buenos Aires' \\\n",
"nro \n",
"1.0 1970.0 1.0 1.0 0.0 10.0 0.0 0.0 \n",
"2.0 1970.0 1.0 2.0 9.4 -999.0 0.9 13.5 \n",
"3.0 1970.0 1.0 3.0 0.0 0.0 87.4 1.4 \n",
"4.0 1970.0 1.0 4.0 3.2 49.0 0.0 0.2 \n",
"5.0 1970.0 1.0 5.0 13.6 3.7 0.0 27.4 \n",
"6.0 1970.0 1.0 6.0 0.0 0.0 0.0 0.0 \n",
"7.0 1970.0 1.0 7.0 0.5 9.4 38.3 11.6 \n",
"8.0 1970.0 1.0 8.0 0.0 0.0 5.5 0.0 \n",
"9.0 1970.0 1.0 9.0 0.0 0.0 0.0 0.0 \n",
"10.0 1970.0 1.0 10.0 0.0 0.0 0.0 0.0 \n",
"11.0 1970.0 1.0 11.0 0.0 0.0 0.0 0.0 \n",
"12.0 1970.0 1.0 12.0 0.0 0.0 0.0 0.0 \n",
"13.0 1970.0 1.0 13.0 0.0 -999.0 13.0 0.4 \n",
"14.0 1970.0 1.0 14.0 0.0 0.0 0.0 0.0 \n",
"15.0 1970.0 1.0 15.0 0.0 4.4 0.0 0.0 \n",
"16.0 1970.0 1.0 16.0 0.0 0.0 0.0 0.0 \n",
"17.0 1970.0 1.0 17.0 6.5 0.0 0.0 0.0 \n",
"18.0 1970.0 1.0 18.0 0.6 0.0 0.0 0.0 \n",
"19.0 1970.0 1.0 19.0 0.0 0.0 0.0 0.0 \n",
"20.0 1970.0 1.0 20.0 0.0 0.0 0.0 0.0 \n",
"21.0 1970.0 1.0 21.0 0.0 8.0 0.0 1.6 \n",
"22.0 1970.0 1.0 22.0 0.0 0.0 0.0 0.0 \n",
"23.0 1970.0 1.0 23.0 55.7 0.0 19.5 13.8 \n",
"24.0 1970.0 1.0 24.0 0.0 20.0 0.0 23.4 \n",
"25.0 1970.0 1.0 25.0 0.0 0.0 0.0 0.0 \n",
"26.0 1970.0 1.0 26.0 0.0 0.0 0.0 29.5 \n",
"27.0 1970.0 1.0 27.0 0.0 0.0 0.0 0.0 \n",
"28.0 1970.0 1.0 28.0 0.2 0.0 60.3 0.6 \n",
"29.0 1970.0 1.0 29.0 0.0 0.0 0.0 0.0 \n",
"30.0 1970.0 1.0 30.0 0.0 0.0 0.0 0.0 \n",
"... ... ... ... ... ... ... ... \n",
"14946.0 2010.0 12.0 2.0 0.0 0.0 -999.0 1.4 \n",
"14947.0 2010.0 12.0 3.0 0.0 0.0 -999.0 0.0 \n",
"14948.0 2010.0 12.0 4.0 0.5 0.0 -999.0 0.0 \n",
"14949.0 2010.0 12.0 5.0 0.0 0.0 -999.0 0.0 \n",
"14950.0 2010.0 12.0 6.0 2.0 0.0 -999.0 0.0 \n",
"14951.0 2010.0 12.0 7.0 0.5 0.0 -999.0 0.0 \n",
"14952.0 2010.0 12.0 8.0 0.0 0.0 -999.0 0.0 \n",
"14953.0 2010.0 12.0 9.0 2.0 0.0 -999.0 0.0 \n",
"14954.0 2010.0 12.0 10.0 2.0 0.0 -999.0 11.3 \n",
"14955.0 2010.0 12.0 11.0 0.0 0.0 -999.0 28.0 \n",
"14956.0 2010.0 12.0 12.0 0.0 0.0 -999.0 0.0 \n",
"14957.0 2010.0 12.0 13.0 0.0 0.0 -999.0 0.0 \n",
"14958.0 2010.0 12.0 14.0 0.0 0.0 -999.0 0.0 \n",
"14959.0 2010.0 12.0 15.0 0.0 0.0 -999.0 0.0 \n",
"14960.0 2010.0 12.0 16.0 0.0 3.0 -999.0 0.0 \n",
"14961.0 2010.0 12.0 17.0 0.1 0.0 -999.0 1.4 \n",
"14962.0 2010.0 12.0 18.0 0.0 0.0 -999.0 0.0 \n",
"14963.0 2010.0 12.0 19.0 0.0 0.0 -999.0 0.0 \n",
"14964.0 2010.0 12.0 20.0 0.0 0.0 -999.0 0.0 \n",
"14965.0 2010.0 12.0 21.0 0.0 0.0 -999.0 0.0 \n",
"14966.0 2010.0 12.0 22.0 0.0 0.0 -999.0 0.0 \n",
"14967.0 2010.0 12.0 23.0 3.0 5.0 -999.0 0.0 \n",
"14968.0 2010.0 12.0 24.0 0.0 0.0 -999.0 0.0 \n",
"14969.0 2010.0 12.0 25.0 0.0 0.0 -999.0 0.0 \n",
"14970.0 2010.0 12.0 26.0 0.0 0.0 -999.0 0.0 \n",
"14971.0 2010.0 12.0 27.0 0.0 0.0 -999.0 0.0 \n",
"14972.0 2010.0 12.0 28.0 0.0 0.0 -999.0 0.0 \n",
"14973.0 2010.0 12.0 29.0 0.0 -999.0 -999.0 0.0 \n",
"14974.0 2010.0 12.0 30.0 0.0 0.0 -999.0 0.0 \n",
"14975.0 2010.0 12.0 31.0 0.0 0.0 -999.0 0.0 \n",
"\n",
" 'Caseros' 'Ceres' 'Concordia' ... 'Pergamino' \\\n",
"nro ... \n",
"1.0 0.0 0.0 0.0 ... 0.0 \n",
"2.0 46.8 29.2 2.0 ... 0.0 \n",
"3.0 4.0 0.9 0.0 ... 0.0 \n",
"4.0 0.0 0.0 0.0 ... 2.8 \n",
"5.0 65.0 10.0 155.0 ... 0.9 \n",
"6.0 46.1 0.9 11.0 ... 0.0 \n",
"7.0 49.6 0.0 11.5 ... 2.3 \n",
"8.0 14.9 3.9 6.2 ... 0.0 \n",
"9.0 0.0 0.0 0.0 ... 0.0 \n",
"10.0 0.0 0.0 0.0 ... 0.0 \n",
"11.0 0.0 0.0 0.0 ... 0.0 \n",
"12.0 0.0 0.0 0.0 ... 1.7 \n",
"13.0 3.8 0.0 8.0 ... 0.0 \n",
"14.0 0.0 0.0 0.0 ... 0.0 \n",
"15.0 0.0 0.0 0.0 ... 0.0 \n",
"16.0 0.0 0.0 0.0 ... 0.0 \n",
"17.0 0.0 0.6 0.0 ... 0.0 \n",
"18.0 0.0 1.7 0.0 ... 0.0 \n",
"19.0 0.0 0.0 0.0 ... 0.0 \n",
"20.0 0.0 0.0 0.0 ... 0.0 \n",
"21.0 0.0 8.0 1.0 ... 32.0 \n",
"22.0 0.0 0.0 0.0 ... 0.0 \n",
"23.0 15.5 28.4 21.5 ... 26.0 \n",
"24.0 0.0 0.0 0.0 ... 0.0 \n",
"25.0 0.0 0.0 0.0 ... 0.0 \n",
"26.0 0.0 0.0 0.0 ... 0.2 \n",
"27.0 0.0 0.0 0.0 ... 0.0 \n",
"28.0 22.7 0.6 9.0 ... 3.3 \n",
"29.0 0.0 0.0 0.0 ... 0.0 \n",
"30.0 0.0 0.0 0.0 ... 0.0 \n",
"... ... ... ... ... ... \n",
"14946.0 -999.0 0.0 0.0 ... 0.0 \n",
"14947.0 -999.0 0.0 0.0 ... 0.0 \n",
"14948.0 -999.0 0.0 0.0 ... 0.0 \n",
"14949.0 -999.0 0.0 0.0 ... 0.0 \n",
"14950.0 -999.0 0.0 0.0 ... 0.0 \n",
"14951.0 -999.0 0.0 0.0 ... 0.0 \n",
"14952.0 -999.0 0.0 0.0 ... 0.0 \n",
"14953.0 -999.0 0.0 0.0 ... 2.5 \n",
"14954.0 -999.0 33.0 0.0 ... 11.0 \n",
"14955.0 -999.0 15.0 3.0 ... 0.0 \n",
"14956.0 -999.0 0.0 -999.0 ... 0.0 \n",
"14957.0 -999.0 0.0 0.0 ... 0.0 \n",
"14958.0 -999.0 0.0 0.0 ... 0.0 \n",
"14959.0 -999.0 0.0 0.0 ... 0.0 \n",
"14960.0 -999.0 0.0 0.0 ... 0.0 \n",
"14961.0 -999.0 0.0 0.0 ... 2.0 \n",
"14962.0 -999.0 74.0 10.0 ... 1.0 \n",
"14963.0 -999.0 4.0 17.0 ... 0.0 \n",
"14964.0 -999.0 0.0 0.0 ... 0.0 \n",
"14965.0 -999.0 0.0 0.0 ... 0.0 \n",
"14966.0 -999.0 0.0 0.0 ... 0.0 \n",
"14967.0 -999.0 4.0 32.0 ... 16.0 \n",
"14968.0 -999.0 0.0 47.0 ... 0.0 \n",
"14969.0 -999.0 0.0 0.0 ... 0.0 \n",
"14970.0 -999.0 0.0 0.0 ... 0.0 \n",
"14971.0 -999.0 0.0 0.0 ... 0.0 \n",
"14972.0 -999.0 0.0 0.0 ... 0.0 \n",
"14973.0 -999.0 0.0 0.0 ... 0.0 \n",
"14974.0 -999.0 0.0 0.0 ... 0.0 \n",
"14975.0 -999.0 0.0 3.0 ... 0.0 \n",
"\n",
" 'Posadas' 'Reconquista' 'Rio Cuarto' 'Rosario' 'San Luis' \\\n",
"nro \n",
"1.0 12.4 0.0 0.1 0.0 12.6 \n",
"2.0 1.2 92.0 0.3 2.3 0.0 \n",
"3.0 0.9 43.9 87.9 0.0 69.0 \n",
"4.0 0.0 0.0 16.0 42.5 11.8 \n",
"5.0 0.0 0.0 3.0 1.0 0.0 \n",
"6.0 5.7 0.5 0.0 0.0 0.0 \n",
"7.0 0.0 107.1 0.0 4.0 0.0 \n",
"8.0 7.6 79.3 0.0 0.0 0.0 \n",
"9.0 0.0 0.0 0.0 0.0 1.4 \n",
"10.0 0.0 0.0 0.0 0.0 0.0 \n",
"11.0 0.0 0.0 0.0 0.0 0.0 \n",
"12.0 0.0 0.0 8.3 2.3 17.0 \n",
"13.0 0.0 0.0 0.0 3.2 0.0 \n",
"14.0 0.0 0.0 0.0 0.0 0.0 \n",
"15.0 0.0 0.0 0.0 0.0 0.0 \n",
"16.0 0.0 0.0 0.0 0.0 0.0 \n",
"17.0 0.0 0.0 23.0 0.0 13.2 \n",
"18.0 0.0 0.0 0.0 0.0 0.0 \n",
"19.0 -999.0 0.0 37.0 0.0 0.0 \n",
"20.0 0.0 0.0 0.0 0.0 0.0 \n",
"21.0 0.0 10.2 0.0 28.5 0.0 \n",
"22.0 0.0 0.0 0.4 0.0 10.5 \n",
"23.0 0.0 24.5 17.0 63.0 30.0 \n",
"24.0 -999.0 0.0 0.0 0.0 0.0 \n",
"25.0 0.0 0.0 13.6 0.5 15.5 \n",
"26.0 0.0 0.6 0.0 0.1 0.0 \n",
"27.0 0.0 0.0 0.0 0.0 0.0 \n",
"28.0 1.0 26.3 0.8 21.0 21.0 \n",
"29.0 39.1 0.0 0.0 0.0 0.0 \n",
"30.0 0.0 0.0 0.0 0.0 0.0 \n",
"... ... ... ... ... ... \n",
"14946.0 27.0 0.2 -999.0 0.5 0.0 \n",
"14947.0 31.0 1.0 18.0 17.0 0.0 \n",
"14948.0 61.0 0.0 9.0 0.6 16.0 \n",
"14949.0 0.0 0.0 0.0 0.0 0.0 \n",
"14950.0 0.0 0.0 23.0 0.0 18.0 \n",
"14951.0 0.0 0.0 1.0 0.0 0.0 \n",
"14952.0 0.0 0.0 0.0 0.0 0.0 \n",
"14953.0 0.0 0.0 -999.0 0.6 0.0 \n",
"14954.0 34.0 49.0 0.0 23.7 0.0 \n",
"14955.0 46.0 77.0 0.0 23.2 0.0 \n",
"14956.0 33.0 0.0 0.0 0.0 0.0 \n",
"14957.0 0.0 0.0 0.0 0.0 0.0 \n",
"14958.0 0.0 0.0 0.0 0.0 0.0 \n",
"14959.0 0.0 0.0 0.0 0.0 0.0 \n",
"14960.0 0.0 0.0 0.0 0.0 0.0 \n",
"14961.0 0.0 0.0 9.0 16.5 0.0 \n",
"14962.0 0.0 3.0 -999.0 12.6 0.0 \n",
"14963.0 16.0 65.0 0.0 0.5 0.0 \n",
"14964.0 -999.0 4.0 0.0 0.0 0.0 \n",
"14965.0 -999.0 0.0 0.0 0.0 0.0 \n",
"14966.0 19.0 0.2 0.0 0.0 0.0 \n",
"14967.0 6.0 22.0 1.0 29.0 20.0 \n",
"14968.0 25.0 7.0 0.0 0.0 0.0 \n",
"14969.0 0.0 0.0 0.0 0.0 0.0 \n",
"14970.0 0.0 0.0 0.0 0.0 0.0 \n",
"14971.0 0.0 0.0 0.0 0.0 16.6 \n",
"14972.0 0.0 0.0 13.0 0.0 1.6 \n",
"14973.0 0.0 0.0 0.0 0.0 22.1 \n",
"14974.0 20.0 0.0 0.0 0.0 4.0 \n",
"14975.0 2.0 0.0 0.0 0.0 -999.0 \n",
"\n",
" 'Sauce Viejo' 'Santa Rosa' 'Tandil' 'Tres Arroyos' \n",
"nro \n",
"1.0 0.0 28.8 0.0 3.8 \n",
"2.0 4.0 4.6 27.5 0.0 \n",
"3.0 0.6 0.0 0.0 0.0 \n",
"4.0 0.0 40.0 3.0 33.2 \n",
"5.0 1.0 0.0 47.0 29.7 \n",
"6.0 0.0 0.0 0.0 0.0 \n",
"7.0 23.1 0.0 3.8 3.5 \n",
"8.0 0.0 0.0 0.0 0.0 \n",
"9.0 0.0 0.0 0.0 0.0 \n",
"10.0 0.0 0.0 0.0 0.0 \n",
"11.0 0.0 0.0 0.0 0.0 \n",
"12.0 0.3 0.0 5.3 5.9 \n",
"13.0 0.6 0.0 0.0 0.2 \n",
"14.0 0.0 0.0 0.0 0.0 \n",
"15.0 0.0 0.0 0.0 0.0 \n",
"16.0 0.0 0.0 1.0 0.0 \n",
"17.0 0.0 0.0 2.0 7.8 \n",
"18.0 0.3 0.0 2.0 0.0 \n",
"19.0 0.0 16.9 0.0 0.0 \n",
"20.0 0.0 23.7 0.0 0.0 \n",
"21.0 15.0 0.0 0.0 0.2 \n",
"22.0 0.0 9.0 0.0 0.0 \n",
"23.0 31.0 15.6 34.0 1.9 \n",
"24.0 0.0 0.0 0.0 9.8 \n",
"25.0 0.0 0.0 13.5 0.0 \n",
"26.0 0.0 0.0 0.0 0.0 \n",
"27.0 2.0 0.0 0.0 0.0 \n",
"28.0 37.7 0.0 0.0 0.0 \n",
"29.0 0.0 0.0 0.0 0.0 \n",
"30.0 0.0 0.0 0.0 0.0 \n",
"... ... ... ... ... \n",
"14946.0 0.0 0.0 -999.0 -999.0 \n",
"14947.0 0.0 0.0 -999.0 -999.0 \n",
"14948.0 0.0 0.0 -999.0 -999.0 \n",
"14949.0 0.0 0.0 -999.0 -999.0 \n",
"14950.0 0.0 0.4 -999.0 -999.0 \n",
"14951.0 0.0 0.0 -999.0 -999.0 \n",
"14952.0 0.0 0.0 -999.0 -999.0 \n",
"14953.0 0.0 0.0 -999.0 -999.0 \n",
"14954.0 10.5 0.0 -999.0 -999.0 \n",
"14955.0 11.0 0.0 -999.0 -999.0 \n",
"14956.0 0.0 0.0 -999.0 -999.0 \n",
"14957.0 0.0 0.0 -999.0 -999.0 \n",
"14958.0 0.0 0.0 -999.0 -999.0 \n",
"14959.0 0.0 0.0 -999.0 -999.0 \n",
"14960.0 0.0 2.1 -999.0 -999.0 \n",
"14961.0 0.0 0.0 -999.0 -999.0 \n",
"14962.0 22.0 0.0 -999.0 -999.0 \n",
"14963.0 13.0 2.0 -999.0 -999.0 \n",
"14964.0 0.0 0.0 -999.0 -999.0 \n",
"14965.0 0.0 0.0 -999.0 -999.0 \n",
"14966.0 0.0 0.0 -999.0 -999.0 \n",
"14967.0 4.0 2.0 -999.0 -999.0 \n",
"14968.0 0.0 0.0 -999.0 -999.0 \n",
"14969.0 0.0 0.0 -999.0 -999.0 \n",
"14970.0 0.0 0.0 -999.0 -999.0 \n",
"14971.0 0.0 0.0 -999.0 -999.0 \n",
"14972.0 0.0 4.0 -999.0 -999.0 \n",
"14973.0 0.0 0.0 -999.0 -999.0 \n",
"14974.0 0.0 3.2 -999.0 -999.0 \n",
"14975.0 0.0 0.0 -999.0 -999.0 \n",
"\n",
"[14975 rows x 33 columns]\n"
]
}
],
"source": [
"#Imprimo los datos por pantalla:\n",
"print(data)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"data=data.replace(-999, np.nan)\n",
"len(data)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"nro\n",
"93.0 8.1\n",
"101.0 7.4\n",
"102.0 3.5\n",
"110.0 3.6\n",
"113.0 18.4\n",
"114.0 1.6\n",
"457.0 0.2\n",
"460.0 2.6\n",
"461.0 1.9\n",
"466.0 14.0\n",
"469.0 22.5\n",
"470.0 65.7\n",
"471.0 19.5\n",
"475.0 0.2\n",
"478.0 0.2\n",
"482.0 0.6\n",
"823.0 15.0\n",
"829.0 5.3\n",
"830.0 2.6\n",
"834.0 32.6\n",
"844.0 2.2\n",
"1188.0 5.5\n",
"1191.0 0.2\n",
"1200.0 3.0\n",
"1209.0 2.6\n",
"1210.0 45.0\n",
"1211.0 5.0\n",
"1212.0 0.5\n",
"1559.0 7.0\n",
"1567.0 37.0\n",
" ... \n",
"11796.0 3.0\n",
"11799.0 0.2\n",
"11802.0 5.0\n",
"12880.0 5.0\n",
"13242.0 10.0\n",
"13252.0 32.0\n",
"13253.0 0.4\n",
"13258.0 24.0\n",
"13259.0 0.4\n",
"13261.0 0.4\n",
"13267.0 6.0\n",
"13269.0 1.0\n",
"13605.0 8.0\n",
"13613.0 6.0\n",
"13614.0 15.0\n",
"13615.0 9.0\n",
"13620.0 5.0\n",
"13621.0 0.6\n",
"13623.0 53.0\n",
"13625.0 3.0\n",
"13626.0 19.0\n",
"13634.0 1.0\n",
"13980.0 0.2\n",
"13981.0 2.6\n",
"13998.0 0.1\n",
"14336.0 3.0\n",
"14338.0 40.0\n",
"14352.0 4.0\n",
"14712.0 16.0\n",
"14714.0 17.0\n",
"Name: Azul, dtype: float64"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Datos del mes utilizado en el histograma para la estacion Azul\n",
"mesh=4\n",
"enero=data[data.mes==mesh]\n",
"enero=enero.Azul[enero.Azul>0]\n",
"len(enero)\n",
"enero"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#Vemos el mes elegido\n",
"if mesh == 1:\n",
" tmes=\"Mes de Enero\"\n",
" print(\"Mes de Enero\") \n",
"elif mesh ==2: \n",
" tmes=\"Mes de Febrero\"\n",
" print(\"Mes de Febrero\")\n",
"elif mesh ==3: \n",
" tmes=\"Mes de Marzo\"\n",
" print(\"Mes de Marzo\")\n",
"elif mesh ==4: \n",
" tmes=\"Mes de Abril\"\n",
" print(\"Mes de Abril\")\n",
"elif mesh ==5: \n",
" tmes=\"Mes de Mayo\"\n",
" print(\"Mes de Mayo\")\n",
"elif mesh ==6: \n",
" tmes=\"Mes de Junio\"\n",
" print(\"Mes de Junio\")\n",
"elif mesh ==7: \n",
" tmes=\"Mes de Julio\"\n",
" print(\"Mes de Julio\")\n",
"elif mesh ==8: \n",
" tmes=\"Mes de Agosto\"\n",
" print(\"Mes de Agosto\")\n",
"elif mesh ==9: \n",
" tmes=\"Mes de Septiembre\"\n",
" print(\"Mes de Septiembre\")\n",
"elif mesh ==10: \n",
" tmes=\"Mes de Octubre\"\n",
" print(\"Mes de Octubre\")\n",
"elif mesh ==11: \n",
" tmes=\"Mes de Noviembre\"\n",
" print(\"Mes de Noviembre\")\n",
"elif mesh ==12: \n",
" tmes=\"Mes de Diciembre\"\n",
" print(\"Mes de Diciembre\")\n",
"else:\n",
" print(\"Ese numero no corresponde a un mes\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#Calculo de Percentiles\n",
"p10=np.percentile(enero,10)\n",
"p20=np.percentile(enero,20)\n",
"p30=np.percentile(enero,30)\n",
"p40=np.percentile(enero,40)\n",
"p50=np.percentile(enero,50)\n",
"p60=np.percentile(enero,60)\n",
"p70=np.percentile(enero,70)\n",
"p80=np.percentile(enero,80)\n",
"p90=np.percentile(enero,90)\n",
"p100=np.percentile(enero,100)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#Ejemplo\n",
"print(p100)\n",
"#Record historico de Azul del mes de enero:\n",
"#anio\n",
"var1=data[data.mes==mesh]\n",
"a=var1.anio[var1.Azul==p100]\n",
"aa=a.values[0]\n",
"#mes\n",
"mesh2=str(mesh)\n",
"#dia\n",
"var2=data[data.mes==mesh]\n",
"d=var2.dia[var2.Azul==p100]\n",
"dd=d.values[0]\n",
"\n",
"#fecha\n",
"aa=np.array(aa,dtype=np.int)\n",
"dd=np.array(dd,dtype=np.int)\n",
"aa=str(aa)\n",
"dd=str(dd)\n",
"fecha=dd+\"/\"+mesh2+\"/\"+aa"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#Plotear el histograma de la estacion 1\n",
"import matplotlib.pyplot as plt\n",
"\n",
"plt.hist(enero, bins=12,facecolor='green', alpha=0.5)\n",
"plt.title(\"Histograma de precipitacion de Azul - \" + tmes)\n",
"plt.xlabel(\"Precipitacion\")\n",
"plt.ylabel(\"Frecuencia\")\n",
"\n",
"#\n",
"plt.text(p100-20, 25, \"Record\", fontsize = 14, weight='bold')\n",
"plt.text(p100-25, 12, fecha, fontsize = 14)\n",
"#plt.show()\n",
"#Guardo\n",
"#plt.savefig(\"hist_azul.eps\")\n",
"plt.savefig(\"../output/hist_azul.png\",bbox_inches='tight')\n",
"plt.show()\n",
"#plt.close()\n",
"\n",
"# VIEJO\n",
"#import potly.potly as py\n",
"#ax.annotate(\"129\",xy=(125,10),arrowprops=dict(arrowstyle=\"->\"))\n",
"#plt.hist(enero[:,1],bins=[0,20,40,60,80,100])\n",
"#EJEMPLO FRANCO\n",
"#a = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])\n",
"#a\n",
"#a[a > 4]"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
"\n",
""
],
"text/plain": [
""
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"#Este css es trabajo de @LorenaABarba y su grupo\n",
"from IPython.core.display import HTML\n",
"css_file = '../../css/personal.css'\n",
"HTML(open(css_file, \"r\").read())"
]
}
],
"metadata": {
"anaconda-cloud": {},
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 0
}