CEN38/Python_scripts
2
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

clean file

Browse Source
This commit is contained in:
Colas Geier 2024-02-26 15:36:37 +01:00
parent 244399633c
commit 4849f1c6ed
1 changed files with 14 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
3_AZALEE/update_psInf1200.py
View File

@ -9,7 +9,6 @@ import geopandas as gpd
def update_psInf1200(df):
from rasterstats import zonal_stats
from pycen import update_to_sql
path0 = '/home/colas/Documents/9_PROJETS/2_PS/'
path2_mnt = path0+'MNT/'
file_mnt5 = 'MNT_5m.tif'
@ -17,9 +16,10 @@ def update_psInf1200(df):
zonal_tmp = zonal_stats(df,path2_mnt+file_mnt5,stats="max",nodata=0,all_touched=True)
zonal = gpd.pd.DataFrame(zonal_tmp)
zonal.columns = ['max_alti']
tmp = gpd.pd.concat([df,zonal], axis=1)
tmp = gpd.pd.concat([df.reset_index(drop=True),zonal], axis=1)
tmp = tmp.set_geometry('geom', crs=2154)
tmp.infeq_1200 = tmp.max_alti <= 1200
PS = tmp[['site_code','infeq_1200']]\
.rename(columns={'site_code':'id_site'})
@ -27,10 +27,21 @@ def update_psInf1200(df):
# update_to_sql(PS,con,'r_infeq_1200m','ps','id_site')
ps_all = gpd.read_postgis('''
SELECT vp.*,s7.infeq_1200
FROM ps."v_pelouseseches_all" vp
LEFT JOIN ps.r_infeq_1200m s7 ON vp.site_code::text = s7.id_site::text
ORDER BY vp.site_code DESC NULLS LAST;
''',con)
isnot_1200 = ps_all.infeq_1200.isna()
update_psInf1200(ps_all[isnot_1200])
metro = gpd.read_file("/home/colas/Documents/tmp/PS_ACTU_METRO/Perimètre_GAM.geojson")
metro.rename_geometry('geom', inplace=True)
ps_agg = gpd.read_postgis('SELECT * FROM "pelouse_seche"."c_ps_inventaire_agreg"',con_bdcen)
ps_all = gpd.read_postgis('SELECT * FROM ps."v_pelouseseches_all"',con)
ps_vps = gpd.read_postgis('SELECT * FROM ps."v_pelouseseches"',con)
ps_noalti = gpd.read_postgis('SELECT * FROM ps."v_pelouseseches_noalti"',con)