zonal_stats.org 19.7 KB

Construction d’un data frame avec les stats zonales

Init

Petit bout de python vide pour initialiser la session. Ca permet d’avoir la version de python utilisée


Import des fichiers dans des DF

Les fichiers shape avec les stats sont

zs_dir='/work/OT/theia/oso/jordi/maassp_data'
zs_rf1=zs_dir+'/zonal_stats_refl_1.shp'
zs_rf2=zs_dir+'/zonal_stats_refl_2.shp'
zs_dem=zs_dir+'/zonal_stats_dem.shp'
zs_mask=zs_dir+'/zonal_stats_mask.shp'

Lire les fichiers dans un DF pandas

import geopandas as gp
df_rf1 = gp.read_file(zs_rf1)
df_rf2 = gp.read_file(zs_rf2)
df_dem = gp.read_file(zs_dem)
df_mask = gp.read_file(zs_mask)
print(len(df_rf1))
print(len(df_rf2))
print(len(df_dem))
print(len(df_mask))

Il paraît que le fichier de DEM a perdu quelques polygones?

Fusion dans un seul DF

On utilise overlay avec intersection

res_intersection = gp.overlay(df_rf1, df_mask, how='identity')
print(len(df_rf1.columns))
print(len(df_mask.columns))
print(len(res_intersection.columns))

Pandas renomme les colonnes qui ont le mm nom dans les 2 DF

for c in res_intersection.columns:
    print(c)

print(res_intersection[['CODE_CULTU_1', 'CODE_CULTU_2']])

Templates


import matplotlib.pyplot as plt
# plt.figure(figsize=(10,5))
# df31.hist(column='RENDNORME')
# plt.savefig(matplot_lib_filename)
matplot_lib_filename