コード例 #1
0
# =============================================================================
# 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)):
コード例 #2
0
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")