Commit 6d0ba736 authored by MEULE Samuel's avatar MEULE Samuel
Browse files

Ajout d'un script TD solinst barologger

# Historique
* [X] Modification de tout les jupyters
* [X] ajouts des fichiers data NKE
* [X] retrait du fichier scriptNKE jupyter
* [X] ajout de la théorie linéaire dans spectral analysis
* [X] Modif TD7 vierge
* [X] Ajout TD lecture aquadopp
* [X] script NKE: Modification de la figure avec le temps pour que cela fontionne sur une ancienne version de matplotlib (sur les serveurs AMU)
* [X] Ajout package latex pour Mybinder dans requirements: ne marche pas
* [X] Ajout du fichier .ipynb pour TD NKE + fichier NKE
* [X] TD lecture de fichier NKE
* [X] Mise en place sous gitlab.osupytheas.fr
* [X] Ajout d'un debut de script TD7 (a remplir) + ajout des fichiers 2 à 10.txt (différentes series temporelles)
* [X] Ajout des notebooks TD2 à TD8 + ajout des fichiers
* [X] Ajout de requirements.txt pour gestion des dépendances
* [X] Ajout du fichier .py et du notebook jupyter .ipynb pour TD1
* [X] Création d'un git pour Stat
parent 118ced84
Serial_number:
2126743
Project ID:
Location:
LEVEL
UNIT: psi
TEMPERATURE
UNIT: C
Date;Time;ms;LEVEL;TEMPERATURE
28/02/2022;10:09:43;0;14,9029;18,229
28/02/2022;10:11:43;0;14,9097;18,858
28/02/2022;10:13:43;0;14,9012;21,490
28/02/2022;10:15:43;0;14,8947;23,736
28/02/2022;10:17:43;0;14,8938;25,067
28/02/2022;10:19:43;0;14,8936;26,348
28/02/2022;10:21:43;0;14,8926;26,540
28/02/2022;10:23:43;0;14,8944;26,753
28/02/2022;10:25:43;0;14,8965;26,532
28/02/2022;10:27:43;0;14,8950;27,290
28/02/2022;10:29:43;0;14,8940;27,263
28/02/2022;10:31:43;0;14,8943;27,384
28/02/2022;10:33:43;0;14,8946;27,585
28/02/2022;10:35:43;0;14,8939;27,337
28/02/2022;10:37:43;0;14,8952;27,573
28/02/2022;10:39:43;0;14,8958;27,727
28/02/2022;10:41:43;0;14,8947;28,641
28/02/2022;10:43:43;0;14,8914;29,818
28/02/2022;10:45:43;0;14,8899;30,595
28/02/2022;10:47:43;0;14,8897;31,250
28/02/2022;10:49:43;0;14,8890;31,851
28/02/2022;10:51:43;0;14,8882;32,400
28/02/2022;10:53:43;0;14,8878;32,665
28/02/2022;10:55:43;0;14,8881;33,151
28/02/2022;10:57:43;0;14,8870;33,717
28/02/2022;10:59:43;0;14,8872;34,136
28/02/2022;11:01:43;0;14,8851;34,234
28/02/2022;11:03:43;0;14,8866;34,216
28/02/2022;11:05:43;0;14,8854;34,459
28/02/2022;11:07:43;0;14,8851;34,599
28/02/2022;11:09:43;0;14,8867;34,469
28/02/2022;11:11:43;0;14,8865;34,605
28/02/2022;11:13:43;0;14,8859;34,707
28/02/2022;11:15:43;0;14,8862;34,688
28/02/2022;11:17:43;0;14,8868;34,629
28/02/2022;11:19:43;0;14,8863;34,610
28/02/2022;11:21:43;0;14,8869;34,668
28/02/2022;11:23:43;0;14,8867;34,508
28/02/2022;11:25:43;0;14,8862;34,471
28/02/2022;11:27:43;0;14,8855;34,531
28/02/2022;11:29:43;0;14,8852;34,616
28/02/2022;11:31:43;0;14,8859;34,442
28/02/2022;11:33:43;0;14,8852;34,458
28/02/2022;11:35:43;0;14,8853;34,069
28/02/2022;11:37:43;0;14,8861;33,968
28/02/2022;11:39:43;0;14,8863;34,009
28/02/2022;11:41:43;0;14,8842;34,051
28/02/2022;11:43:43;0;14,8847;34,063
28/02/2022;11:45:43;0;14,8836;33,984
28/02/2022;11:47:43;0;14,8859;33,624
28/02/2022;11:49:43;0;14,8840;33,483
28/02/2022;11:51:43;0;14,8854;33,202
28/02/2022;11:53:43;0;14,8862;32,798
28/02/2022;11:55:43;0;14,8860;32,470
28/02/2022;11:57:43;0;14,8851;32,215
28/02/2022;11:59:43;0;14,8862;32,075
28/02/2022;12:01:43;0;14,8852;32,061
28/02/2022;12:03:43;0;14,8848;31,672
28/02/2022;12:05:43;0;14,8843;31,174
28/02/2022;12:07:43;0;14,8845;31,183
28/02/2022;12:09:43;0;14,8852;31,034
28/02/2022;12:11:43;0;14,8845;30,842
28/02/2022;12:13:43;0;14,8831;30,716
28/02/2022;12:15:43;0;14,8823;30,470
28/02/2022;12:17:43;0;14,8828;30,067
28/02/2022;12:19:43;0;14,8815;29,863
28/02/2022;12:21:43;0;14,8816;29,793
28/02/2022;12:23:43;0;14,8821;29,733
28/02/2022;12:25:43;0;14,8812;29,614
28/02/2022;12:27:43;0;14,8812;29,624
28/02/2022;12:29:43;0;14,8799;29,436
28/02/2022;12:31:43;0;14,8801;28,978
28/02/2022;12:33:43;0;14,8811;28,799
28/02/2022;12:35:43;0;14,8790;28,680
28/02/2022;12:37:43;0;14,8795;28,310
28/02/2022;12:39:43;0;14,8823;27,376
28/02/2022;12:41:43;0;14,8834;26,568
28/02/2022;12:43:43;0;14,8860;25,418
28/02/2022;12:45:43;0;14,8880;24,046
28/02/2022;12:47:43;0;14,8894;22,664
28/02/2022;12:49:43;0;14,8899;21,541
28/02/2022;12:51:43;0;14,8894;20,634
28/02/2022;12:53:43;0;14,8907;19,854
28/02/2022;12:55:43;0;14,8901;19,371
28/02/2022;12:57:43;0;14,8892;19,066
28/02/2022;12:59:43;0;14,8879;18,831
28/02/2022;13:01:43;0;14,8883;18,556
28/02/2022;13:03:43;0;14,8875;18,361
28/02/2022;13:05:43;0;14,8876;18,223
28/02/2022;13:07:43;0;14,8860;18,073
28/02/2022;13:09:43;0;14,8855;17,943
28/02/2022;13:11:43;0;14,8851;17,963
28/02/2022;13:13:43;0;14,8850;17,861
28/02/2022;13:15:43;0;14,8842;17,757
28/02/2022;13:17:43;0;14,8842;17,654
28/02/2022;13:19:43;0;14,8835;17,626
28/02/2022;13:21:43;0;14,8835;17,541
28/02/2022;13:23:43;0;14,8825;17,514
28/02/2022;13:25:43;0;14,8832;17,441
28/02/2022;13:27:43;0;14,8831;17,369
28/02/2022;13:29:43;0;14,8826;17,366
28/02/2022;13:31:43;0;14,8816;17,366
28/02/2022;13:33:43;0;14,8805;17,324
28/02/2022;13:35:43;0;14,8798;17,275
28/02/2022;13:37:43;0;14,8795;17,210
28/02/2022;13:39:43;0;14,8797;17,013
28/02/2022;13:41:43;0;14,8794;16,945
28/02/2022;13:43:43;0;14,8791;16,914
28/02/2022;13:45:43;0;14,8795;16,890
28/02/2022;13:47:43;0;14,8780;16,862
28/02/2022;13:49:43;0;14,8785;16,645
28/02/2022;13:51:43;0;14,8780;16,545
28/02/2022;13:53:43;0;14,8780;16,507
28/02/2022;13:55:43;0;14,8778;16,491
28/02/2022;13:57:43;0;14,8771;16,537
28/02/2022;13:59:43;0;14,8767;16,515
28/02/2022;14:01:43;0;14,8767;16,478
28/02/2022;14:03:43;0;14,8761;16,509
28/02/2022;14:05:43;0;14,8757;16,501
28/02/2022;14:07:43;0;14,8750;16,519
28/02/2022;14:09:43;0;14,8752;16,380
28/02/2022;14:11:43;0;14,8745;16,305
28/02/2022;14:13:43;0;14,8747;16,306
28/02/2022;14:15:43;0;14,8738;16,276
28/02/2022;14:17:43;0;14,8732;16,245
28/02/2022;14:19:43;0;14,8736;16,232
28/02/2022;14:21:43;0;14,8734;16,216
28/02/2022;14:23:43;0;14,8729;16,164
28/02/2022;14:25:43;0;14,8722;16,156
28/02/2022;14:27:43;0;14,8719;16,120
28/02/2022;14:29:43;0;14,8717;16,049
28/02/2022;14:31:43;0;14,8715;16,028
28/02/2022;14:33:43;0;14,8719;15,980
28/02/2022;14:35:43;0;14,8711;15,936
28/02/2022;14:37:43;0;14,8711;15,765
28/02/2022;14:39:43;0;14,8711;15,691
28/02/2022;14:41:43;0;14,8700;15,575
28/02/2022;14:43:43;0;14,8704;15,459
28/02/2022;14:45:43;0;14,8694;15,415
28/02/2022;14:47:43;0;14,8695;15,371
28/02/2022;14:49:43;0;14,8695;15,322
28/02/2022;14:51:43;0;14,8689;15,281
28/02/2022;14:53:43;0;14,8691;15,222
28/02/2022;14:55:43;0;14,8684;15,160
28/02/2022;14:57:43;0;14,8675;15,094
28/02/2022;14:59:43;0;14,8676;15,048
28/02/2022;15:01:43;0;14,8675;15,002
28/02/2022;15:03:43;0;14,8679;14,929
28/02/2022;15:05:43;0;14,8668;14,899
28/02/2022;15:07:43;0;14,8674;14,899
28/02/2022;15:09:43;0;14,8666;14,832
28/02/2022;15:11:43;0;14,8669;14,764
28/02/2022;15:13:43;0;14,8654;14,692
28/02/2022;15:15:43;0;14,8664;14,650
28/02/2022;15:17:43;0;14,8653;14,647
28/02/2022;15:19:43;0;14,8646;14,668
28/02/2022;15:21:43;0;14,8640;14,706
28/02/2022;15:23:43;0;14,8645;14,755
28/02/2022;15:25:43;0;14,8633;14,785
28/02/2022;15:27:43;0;14,8639;14,823
28/02/2022;15:29:43;0;14,8631;14,788
28/02/2022;15:31:43;0;14,8628;14,777
28/02/2022;15:33:43;0;14,8632;14,738
28/02/2022;15:35:43;0;14,8625;14,743
28/02/2022;15:37:43;0;14,8622;14,724
28/02/2022;15:39:43;0;14,8618;14,668
28/02/2022;15:41:43;0;14,8613;14,595
28/02/2022;15:43:43;0;14,8607;14,506
28/02/2022;15:45:43;0;14,8606;14,465
28/02/2022;15:47:43;0;14,8601;14,439
28/02/2022;15:49:43;0;14,8596;14,396
28/02/2022;15:51:43;0;14,8591;14,399
28/02/2022;15:53:43;0;14,8592;14,377
28/02/2022;15:55:43;0;14,8579;14,310
28/02/2022;15:57:43;0;14,8585;14,199
28/02/2022;15:59:43;0;14,8576;14,091
28/02/2022;16:01:43;0;14,8581;14,020
28/02/2022;16:03:43;0;14,8580;13,900
28/02/2022;16:05:43;0;14,8572;13,755
28/02/2022;16:07:43;0;14,8570;13,680
28/02/2022;16:09:43;0;14,8565;14,607
28/02/2022;16:11:43;0;14,8550;15,089
28/02/2022;16:13:43;0;14,8483;15,405
28/02/2022;16:15:43;0;14,8448;15,778
28/02/2022;16:17:43;0;14,8431;16,299
28/02/2022;16:19:43;0;14,8454;16,511
28/02/2022;16:21:43;0;14,8474;15,851
28/02/2022;16:23:43;0;14,8506;15,123
28/02/2022;16:25:43;0;14,8478;15,887
28/02/2022;16:27:43;0;14,8512;15,425
Modification de tout les jupyters
Ajout d'un script TD solinst barologger
# Historique
* [X] Modification de tout les jupyters
* [X] ajouts des fichiers data NKE
* [X] retrait du fichier scriptNKE jupyter
* [X] ajout de la théorie linéaire dans spectral analysis
......
......@@ -506,8 +506,8 @@ aquadopp= AQD()
##########################################################
# Choose the working directory
os.chdir('../../../../workspace/marseille')
aquadopp.filename='18032102.hdr'
os.chdir('../../../../workspace/marseille/AQP')
aquadopp.filename='KIKI02.hdr'
#############################
### READ DATA FILES
......
......@@ -49,13 +49,13 @@ class TDNKE():
[adr, ext]=os.path.splitext(os.path.basename(self.filename))
dateparse = lambda X: datetime.strptime(X,'%d/%m/%Y\t%H:%M:%S:%f')
self.mydata = pd.read_table(self.filename,header=4,
decimal='.',
decimal=',',
usecols=[0,1,2, 3],
names=["Temperature", "level","thedate","thetime"],
parse_dates={"Date": ["thedate","thetime"]},
date_parser=dateparse,
delim_whitespace = True,
#delimiter='\s+',
#delim_whitespace = True,
delimiter='\s+|;',
#skiprows=100202,
#nrows=2537035,
skip_blank_lines=True)
......@@ -352,9 +352,9 @@ class TDNKE():
deg = u'\xb0'
data_NKE=TDNKE()
# Choose the working directory
os.chdir('../../../workspace/marseille')
data_NKE.filename='Sonde standard_30098_20190309_102905.txt'
data_NKE.Sample_Frequency=4
os.chdir('../../../../workspace/marseille/NKE')
data_NKE.filename='Sonde standard_32053_20220228_161959.txt'
data_NKE.Sample_Frequency=10
data_NKE.burst=2400 # Nombre de point pour l'analyse spectrale
data_NKE.altmes=0.01
##########################################################
......
#!/usr/bin/python
# coding: utf-8
#############################################"
## MODULE #####"
#############################################"
import matplotlib.pyplot as plt
from tkinter import filedialog
from datetime import datetime
import numpy.matlib as matlib
import matplotlib.dates as dt
import scipy.stats as stat
import tkinter as Tk
import pandas as pd
import numpy as np
import pylab as py
import scipy.io
import pickle
import math
import sys
import os
#
###
# # #
#
#
##########################################################################################################################################################
#
#
# # #
###
#
class meteo(object):
#############################################################
def _load_solinst(self, filename):
""" OPEN AND READ A FILE FROM SOLINST BAROLOGGER , VERSION PANDA """
dateparse = lambda X: datetime.strptime(X,'%d/%m/%Y \t%H:%M:%S')
self.mydata = pd.read_table(self.filename,header=11,
decimal=',',
usecols=[0,1,3,4],
names=["thedate","thetime","Pressure", 'Temperature'],
parse_dates={"Date": ["thedate","thetime"]},
date_parser=dateparse,
#delim_whitespace = True,
delimiter=';',
#skiprows=100202,
#nrows=2537035,
skip_blank_lines=True,
encoding='Latin-1')
self.mydata.set_index('Date', inplace=True)
#############################################################
def _figuremeteo(self, y, text):
# Figure creation
taille = (10,5) # 10*tdpi x 4*tdpi -> 1200x800 si tdpi =200
tdpi=200
fig1=plt.figure(figsize=taille)
ax = fig1.add_subplot(111)
plt.rc('text', usetex=True)
plt.rc('font', family='serif')
plt.plot(self.mydata.index,y, 'r.')
plt.xlabel(r'Time (m)', fontsize=16)
plt.ylabel(text, fontsize=16)
date_format = dt.DateFormatter('%y/%m/%d %H:%M')
plt.gca().xaxis.set_major_formatter(date_format)
plt.xticks(rotation=45)
plt.autoscale(enable=True, axis='both', tight=True)
plt.tight_layout()
#--------------------------------------------------
# Save the figure
#--------------------------------------------------
figname=self.adr+ '_'+text+ '.png'
fig1.savefig(figname,dpi=tdpi)
print( 'Figure :', figname, ' :', tdpi,' dpi......Ok')
#plt.show()
#############################
### DEF VARIABLES
#############################
deg = u'\xb0'
AtmoP=meteo()
##########################################################
### OPEN FILE ###
##########################################################
root=Tk.Tk()
root.update()
# Choose the working directory
os.chdir('../../../../workspace/marseille/Barologger')
AtmoP.filename = filedialog.askopenfilename(title="Ouvrir un fichier barolog solinst",
filetypes=[('solinst files','.csv'),('all files','.*')])
rep=os.path.dirname(AtmoP.filename)
[AtmoP.adr, ext]=os.path.splitext(os.path.basename(AtmoP.filename))
try:
os.mkdir(AtmoP.adr)
except:
pass
os.chdir(AtmoP.adr)
root.destroy()
# Read File
AtmoP._load_solinst(AtmoP.filename)
# Convert into HPa
AtmoP.mydata['Pressure']= 68.947572932 *AtmoP.mydata['Pressure']
# Figures
AtmoP._figuremeteo(AtmoP.mydata.Pressure,'Atmospheric Pressure (HPa)')
AtmoP._figuremeteo(AtmoP.mydata.Temperature,'Atmospheric Temperature (degree C)')
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment