#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on September 2022

@author: Phil Vernant philippe.vernant@umontpellier.fr
"""

# import libraries
import numpy as np
import pandas as pd
import matplotlib
import matplotlib.pyplot as plt

#-----------------------------------------------------------------------------------------

gps = pd.read_csv('gcl26092022.csv',names=['point','E','N','Z'],index_col=False) # lecture du fichier
gps = gps[gps['point'].str.contains('base')==False] # on elimine la ligne contenant la base
gps['sta'] =  gps.point.str[0] # on creer une nouvelle colonne contenant juste la premiere lettre
stations = gps.sta.unique() # on liste les stations
print ('points mesurées :', stations)

#creation d'un tableau de donnees pour les valeurs moyennes  
gps_m = pd.DataFrame(columns=['sta','n_obs','E','N','Z','sE','sN','sZ'])
gps_m.sta = stations

for i, sta in enumerate(stations):   # calcul des moyennes et écarts types
  print ('processing :',sta)
  gps_m.n_obs[i] = len(gps.E[gps.sta==sta])
  gps_m.E[i] = np.mean(gps.E[gps.sta==sta])
  gps_m.N[i] = np.mean(gps.N[gps.sta==sta])
  gps_m.Z[i] = np.mean(gps.Z[gps.sta==sta])
  gps_m.sE[i] = np.std(gps.E[gps.sta==sta])
  gps_m.sN[i] = np.std(gps.N[gps.sta==sta])
  gps_m.sZ[i] = np.std(gps.Z[gps.sta==sta])

print(gps_m)