Los datos en total hacen una matriz de 468x4. La idea es sacar el promedio del caudal de cada año, por lo cual se deben tomar los 12 valores (mes 1 al 12) de cada año y sacarle el promedio, así se debe obtener para el año 1981, 1982...
Estación | Año | Mes | Caudal |
---|---|---|---|
23087210 | 1981 | 1 | 287.6 |
23087210 | 1981 | 2 | 246.7 |
23087210 | 1981 | 3 | 231.7 |
23087210 | 1981 | 4 | 257 |
23087210 | 1981 | 5 | 376.3 |
23087210 | 1981 | 6 | 312 |
23087210 | 1981 | 7 | 284.6 |
23087210 | 1981 | 8 | 294.5 |
23087210 | 1981 | 9 | 348.7 |
23087210 | 1981 | 10 | 528.6 |
23087210 | 1981 | 11 | 531.1 |
23087210 | 1981 | 12 | 448.7 |
23087210 | 1982 | 1 | 286.9 |
23087210 | 1982 | 2 | 380.5 |
23087210 | 1982 | 3 | 314 |
23087210 | 1982 | 4 | 385.3 |
23087210 | 1982 | 5 | 520.7 |
23087210 | 1982 | 6 | 493.9 |
23087210 | 1982 | 7 | 366.1 |
23087210 | 1982 | 8 | 331.8 |
23087210 | 1982 | 9 | 410 |
23087210 | 1982 | 10 | 475.5 |
23087210 | 1982 | 11 | 484.8 |
23087210 | 1982 | 12 | 434.7 |
He intentado plantear el código pero sólo obtengo el valor de los primeros 12 elementos correspondientes al año 1981.
import numpy as np
import matplotlib.pyplot as plt
mat = np.loadtxt('23087211.txt', skiprows = 1, unpack = False)
matc1 = mat[:, 3:468]
# Matriz para los promedios de los caudales
prom_anual = np.zeros(468)
for i in range(468):
vector_anual = matc1[0:12,i]
panual = np.average(vector_anual)
prom_anual[i] = panual