# ============================================================================= # Lecture du répertoire ; construction d'une liste avec nom des stations # ============================================================================= directory = os.path.join('/Users/antoineleblevec/Desktop/hodor_wgtn') rep = os.path.abspath(os.path.expanduser(directory)) files = os.listdir(rep) files.sort() os.chdir(directory) sat = [] df = [] rlon = [] rlat = [] for i in range(len(files)): name = files[i].split('_') sat.append(name[1]) df.append(f1.read(files[i], epoque1, epoque2)) df[i] = df[i].set_index("tsn") df[i].columns = [ 'el_{0}'.format(sat[i]), 'az_{0}'.format(sat[i]), 'tec_{0}'.format(sat[i]) ] rlon.append(f1.lecture_lon(files[i])) rlat.append(f1.lecture_lat(files[i])) dftot = reduce(lambda x, y: pd.merge(x, y, on="tsn"), df) df_rlat = pd.DataFrame([rlat]) df_rlon = pd.DataFrame([rlon]) df_el = pd.DataFrame() df_az = pd.DataFrame() df_tec = pd.DataFrame() for j in range(len(files)):
tda = [] # ============================================================================= # Accès, lecture des données et création des DataFrames # ============================================================================= # chemin vers données directory = os.path.join('/Users/antoineleblevec/Desktop/2004_article') rep = os.path.abspath(os.path.expanduser(directory)) files = os.listdir(rep) files.sort() os.chdir(directory) # lecture des données et stockage des latitudes et longitudes des stations présentes for file in files: name = file.split('_') da_fr = f1.read(file, e1, e2) if len(da_fr) < 10: # données incomplètes pour cette station continue else: df.append(da_fr) station.append(name[0]) lon_station.append(f1.lecture_lat_lon(file)[1]) lat_station.append(f1.lecture_lat_lon(file)[0]) df_lat_station = pd.DataFrame([lat_station]) df_lon_station = pd.DataFrame([lon_station]) df_lat_station.columns = station df_lon_station.columns = station if len(df) == 0: # satellite inutilisable sys.exit("aucunes données correspondent")