def txtMidas(link, date_deb=-1):
"""
Fonction qui va analyser un fichier .xyz pour sortir des graphiques sur l'évolution de la vitesse calculer par
la méthode MIDAS en fonction du nombre de jour de mesure.
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
# on varie la longueur du pas pour faciliter le calcul
nb1 = np.arange(50, 400, 10)
nb2 = np.arange(400, 1000, 50)
nb3 = np.arange(1000, min(len(data), 1900) + 100, 100)
nb4 = np.arange(2000, min(len(data), 5000), 500)
nb_mesure = np.concatenate((nb1, nb2, nb3, nb4, np.array([len(data)])),
axis=0)
path = '../Resultatcsv/MIDAS/'
with open(path + link[-12:-8] + '.csv', 'w') as filename:
file = csv.writer(filename)
# pour chaque nombre de mesures, on effectue le calcul et on rempli chaques listes
for i in nb_mesure:
data = formatage(link, nb_jour=i, date_debut=date_deb)
vitesseMidas = globalMidas(data)
vitesseMidas_E = vitesseMidas[0][0]
ecartype_E = vitesseMidas[1][0]
vitesseMidas_N = vitesseMidas[0][1]
ecartype_N = vitesseMidas[1][1]
vitesseMidas_h = vitesseMidas[0][2]
ecartype_h = vitesseMidas[1][2]
file.writerow([
i, vitesseMidas_E, ecartype_E, vitesseMidas_N, ecartype_N,
vitesseMidas_h, ecartype_h
])
def boucle(file, link, periode, robust, nb_mesures, date_deb):
# pour chaque nombre de mesures, on effectue le calcul et on rempli chaques listes
for i in nb_mesures:
try:
data = formatage(link, nb_jour=i, date_debut=date_deb)
liste_MC = moindreCarres(data, periode, robust=robust)
vitesse_E = liste_MC[1][2][0]
vitesse_N = liste_MC[2][2][0]
vitesse_h = liste_MC[3][2][0]
ecart_E = liste_MC[1][2][1]
ecart_N = liste_MC[2][2][1]
ecart_h = liste_MC[3][2][1]
file.writerow([
i, vitesse_E, ecart_E, vitesse_N, ecart_N, vitesse_h, ecart_h
])
except np.linalg.LinAlgError:
pass
def graphiqueData(link):
"""
Fonction qui affiche la position d'un jeu de données sur les trois axes de coordonnées.
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
t = data[:, 1]
E = data[:, 2]
N = data[:, 3]
h = data[:, 4]
py.figure(5)
py.plot(t, E, 'g', label='position sur E')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe E en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\position\\" + link[-12:-8] + "_E")
py.close()
py.figure(6)
py.plot(t, N, 'g', label='position sur N')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe N en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\position\\" + link[-12:-8] + "_N")
py.close()
py.figure(7)
py.plot(t, h, 'g', label='position sur h')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe h en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\position\\" + link[-12:-8] + "_h")
py.close()
def txtMC(link, periode=[365.25, 365.25 / 2], robust=False, date_deb=-1):
data = formatage(link, date_deb)
# une liste avec un pas régulier pour afficher en fonction du nombre de mesure
nb1 = np.arange(50, 400, 10)
nb2 = np.arange(400, 1000, 50)
nb3 = np.arange(1000, min(len(data), 3000) + 100, 100)
nb4 = np.arange(3100, min(len(data), 7000), 300)
nb_mesures = np.concatenate((nb1, nb2, nb3, nb4, np.array([len(data)])),
axis=0)
path = '../Resultatcsv/MC_palier/'
if robust:
with open(path + link[-12:-8] + '_robuste.csv', 'w') as filename:
file = csv.writer(filename)
boucle(file, link, periode, robust, nb_mesures, date_deb)
else:
with open(path + link[-12:-8] + '.csv', 'w') as filename:
file = csv.writer(filename)
boucle(file, link, periode, robust, nb_mesures, date_deb)
from Moindres_Carres.MC_scipy import test_MC, test_ls, test_saut_ls
from Moindres_Carres.vitesse_MC import moindreCarres
from Traitement.formatage import formatage
import numpy as np
import pylab as py
path = "C:\\Users\\Hul\\Desktop\\Projet_recherche\\TOUL_igs.xyz"
data = formatage(path)
data[140][2] = -50000
a = test_MC(data)
b = moindreCarres(data, periode=[365.25, 365.25/2], robust=True, extend=True)
#Sur axe E
t = data[:,1]
E = data[:,2]
residu_sc = a[0].fun
residu_mc = b[2][len(b[2])-1]
py.figure(0)
#py.plot(t, E, 'g', label='position sur E')
py.plot(t, E + residu_sc, 'r', label='position sur E prédit par scipy')
py.plot(t, residu_mc[1], 'b', label='position sur E predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel("Position sur l'axe E en fonction du nombre de jour de mesure (m)")
py.legend()
py.show()
def graphiqueMidas(link):
"""
Fonction qui va analyser un fichier .xyz pour sortir des graphiques sur l'évolution de la vitesse calculer par
la méthode MIDAS en fonction du nombre de jour de mesure.
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
# les listes que l'on affichera ensuite
vitesseMidas_E = []
vitesseMidas_N = []
vitesseMidas_h = []
ecartype_E = []
ecartype_N = []
ecartype_h = []
# on varie la longueur du pas pour faciliter le calcul
nb1 = np.arange(50, 400, 10)
nb2 = np.arange(400, 1000, 50)
nb3 = np.arange(1000, min(len(data), 1900) + 100, 100)
nb4 = np.arange(2000, len(data), 500)
nb_mesure = np.concatenate((nb1, nb2, nb3, nb4), axis=0)
# pour chaque nombre de mesures, on effectue le calcul et on rempli chaques listes
for i in nb_mesure:
print(i)
data = formatage(link, nb_jour=i)
vitesseMidas = globalMidas(data)
vitesseMidas_E.append(vitesseMidas[0][0])
ecartype_E.append(vitesseMidas[1][0])
vitesseMidas_N.append(vitesseMidas[0][1])
ecartype_N.append(vitesseMidas[1][1])
vitesseMidas_h.append(vitesseMidas[0][2])
ecartype_h.append(vitesseMidas[1][2])
# figure avec les vitesses sur chaque axes
py.figure(3)
py.plot(nb_mesure, vitesseMidas_E, 'g', label="vitesse sur l'axe E")
py.plot(nb_mesure, vitesseMidas_N, 'r', label="vitesse sur l'axe N")
py.plot(nb_mesure, vitesseMidas_h, 'b', label="vitesse sur l'axe h")
py.xlabel("Temps en jour de mesure (j)")
py.ylabel("Vitesse en fonction du nombre de jour de mesure (m/j)")
py.legend()
py.savefig("..\\graph\\MIDAS\\" + link[-12:-8] + "_vitesse")
py.close()
# figure avec les écart-types de la mesures des vitesses
py.figure(4)
py.plot(nb_mesure, ecartype_E, 'g', label="ecart-type sur l'axe E")
py.plot(nb_mesure, ecartype_N, 'r', label="ecart-type sur l'axe N")
py.plot(nb_mesure, ecartype_h, 'b', label="ecart-type sur l'axe h")
py.xlabel("Temps en jour de mesure (j)")
py.ylabel("Ecart-type en fonction du nombre de jour de mesure (m/j)")
py.legend()
py.savefig("..\\graph\\MIDAS\\" + link[-12:-8] + "_ecart_type")
py.close()
def graphiqueCompMC(link, periode, robust=False):
"""
Fonction qui la position prédit par les moindres carrées et par la fonction least_squares de scipy
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
P = np.zeros((len(data), len(data)))
for i in range(len(data)):
P[i, i] = 1
t1 = time.time()
MC = moindreCarres(data, periode, robust=robust, extend=True)
t2 = time.time()
print("temps Simeon " + str(t2 - t1))
r = test_MC(data)
t3 = time.time()
print("temps ordi " + str(t3 - t2))
t = data[:, 1]
E = data[:, 2]
N = data[:, 3]
h = data[:, 4]
E_test = r[0].fun + E
N_test = r[1].fun + N
h_test = r[2].fun + h
t0 = np.mean(data[:, 1])
E_mc = MC[1][-1][1]
N_mc = MC[2][-1][1]
h_mc = MC[3][-1][1]
py.figure(12)
py.plot(t, E, 'g', label='position sur E')
py.plot(t, E_test, 'r', label='position sur E scipy')
py.plot(t,
E_mc,
'b',
label='position sur E predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe E en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\comparaison_simple\\" + link[-12:-8] + "_E")
py.close()
py.figure(13)
py.plot(t, N, 'g', label='position sur N')
py.plot(t, N_test, 'r', label='position sur N prédit par scipy')
py.plot(t,
N_mc,
'b',
label='position sur N predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe N en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\comparaison_simple\\" + link[-12:-8] + "_N")
py.close()
py.figure(14)
py.plot(t, h, 'g', label='position sur h')
py.plot(t, h_test, 'r', label='position sur h prédit par scipy')
py.plot(t,
h_mc,
'b',
label='position sur h predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe h en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\comparaison_simple\\" + link[-12:-8] + "_h")
py.close()
def graphiqueTot(link):
"""
Fonction qui affiche la position d'un jeu de données sur les trois axes de coordonnées.
Il affiche en plus la position prédit par les moindres carrées et par MIDAS
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
vitesseMidas = globalMidas(data)
MC = moindreCarres(data, [365.25, 365.25 / 2])
t = data[:, 1]
E = data[:, 2]
N = data[:, 3]
h = data[:, 4]
E_midas = []
N_midas = []
h_midas = []
E_mc = MC[1][-1][1]
N_mc = MC[2][-1][1]
h_mc = MC[3][-1][1]
py.figure(8)
py.plot(t, E, 'g', label='position sur E')
py.plot(t, E_midas, 'r', label='position sur E prédit par MIDAS')
py.plot(t,
E_mc,
'b',
label='position sur E predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe E en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\prediction\\" + link[-12:-8] + "_E")
py.close()
py.figure(9)
py.plot(t, N, 'g', label='position sur N')
py.plot(t, N_midas, 'r', label='position sur N prédit par MIDAS')
py.plot(t,
N_mc,
'b',
label='position sur N predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe N en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\prediction\\" + link[-12:-8] + "_N")
py.close()
py.figure(10)
py.plot(t, h, 'g', label='position sur h')
py.plot(t, h_midas, 'r', label='position sur h prédit par MIDAS')
py.plot(t,
h_mc,
'b',
label='position sur h predit par les paramètres moindres carrés')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Position sur l'axe h en fonction du nombre de jour de mesure (m)")
py.legend()
py.savefig("..\\graph\\prediction\\" + link[-12:-8] + "_h")
py.close()
def graphiqueUnique(link):
"""
Fonction qui va analyser un fichier .xyz pour sortir des graphiques sur l'évolution de la vitesse calculée par
moindres carrées en fonction du temps et comparer cela avec la vitesse obtenue avec la méthode MIDAS.
Elle les enregistre dans le fichier dédié.
:param link: chemin vers le fichier .xyz
:type link: str
"""
data = formatage(link)
vitesseMidas = globalMidas(data)
liste_MC_final = moindreCarres(data, [365.25, 365.25 / 2])
# une liste avec un pas régulier pour afficher en fonction du nombre de mesure
nb_mesures = np.arange(350, len(data), 50)
# les listes que l'on affichera ensuite
vitesse_E = []
vitesse_N = []
vitesse_h = []
# pour chaque nombre de mesures, on effectue le calcul et on rempli chaques listes
for i in nb_mesures:
data = formatage(link, nb_jour=i)
liste_MC = moindreCarres(data, [365.25, 365.25 / 2])
vitesse_E.append(liste_MC[1][2][0])
vitesse_N.append(liste_MC[2][2][0])
vitesse_h.append(liste_MC[3][2][0])
# figure de la vitesse sur l'axe E
py.figure(0)
py.plot(nb_mesures, vitesse_E, 'g', label='vitesse moindres carrées sur E')
py.plot(nb_mesures,
len(nb_mesures) * [vitesseMidas[0][0]],
'r',
label='vitesse MIDAS sur E')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Vitesse sur l'axe E en fonction du nombre de jour de mesure (m/j)")
py.legend()
py.savefig("..\\graph\\MC\\" + link[-12:-8] + "_E")
py.close()
# figure de la vitesse sur l'axe N
py.figure(1)
py.plot(nb_mesures, vitesse_N, 'g', label='vitesse moindres carrées sur N')
py.plot(nb_mesures,
len(nb_mesures) * [vitesseMidas[0][1]],
'r',
label='vitesse MIDAS sur N')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Vitesse sur l'axe N en fonction du nombre de jour de mesure (m/j)")
py.legend()
py.savefig("..\\graph\\MC\\" + link[-12:-8] + "_N")
py.close()
# figure de la vitesse sur l'axe h
py.figure(2)
py.plot(nb_mesures, vitesse_h, 'g', label='vitesse moindres carrées sur h')
py.plot(nb_mesures,
len(nb_mesures) * [vitesseMidas[0][2]],
'r',
label='vitesse MIDAS sur h')
py.xlabel("Temps en jour de mesure (j)")
py.ylabel(
"Vitesse sur l'axe h en fonction du nombre de jour de mesure (m/j)")
py.legend()
py.savefig("..\\graph\\MC\\" + link[-12:-8] + "_h")
py.close()