From 7582fd5a2a89db58d04684e80ff0a223d1ffd9c0 Mon Sep 17 00:00:00 2001 From: Colas Geier Date: Mon, 26 Feb 2024 15:33:29 +0100 Subject: [PATCH] =?UTF-8?q?int=C3=A9gration=20des=20donn=C3=A9es=20zh=20PL?= =?UTF-8?q?U=20metro?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- 3_AZALEE/tmp/zh_plu_metro.py | 165 +++++++++++++++++++++++++++++++++++ 1 file changed, 165 insertions(+) diff --git a/3_AZALEE/tmp/zh_plu_metro.py b/3_AZALEE/tmp/zh_plu_metro.py index 25a6526..c28b119 100644 --- a/3_AZALEE/tmp/zh_plu_metro.py +++ b/3_AZALEE/tmp/zh_plu_metro.py @@ -635,6 +635,24 @@ if __name__ == "__main__": ############################ ######### NEW maj ########## + ### 38GAM0030 + GAM030 = maj_dt_old[maj_dt_old.id_origine=='38GAM0030'].copy() + GAM030.date = '2017-01-01' + GAM030geom = (GAM030[['site_code','id_origine','date','auteur','geom']] + .copy() + .merge(info,left_on='id_origine',right_on='CODE_SITE',how='left') + .set_geometry('geom',crs=2154) + .drop(columns=['CODE_SITE','nom'])) + insertNewGeom(GAM030geom) + ### 38GAM0106 + GAM106 = maj_dt_old[maj_dt_old.id_origine=='38GAM0106'].copy() + GAM106.date = '2017-01-01' + GAM106geom = (GAM106[['site_code','id_origine','date','auteur','geom']] + .copy() + .merge(info,left_on='id_origine',right_on='CODE_SITE',how='left') + .set_geometry('geom',crs=2154) + .drop(columns=['CODE_SITE','nom'])) + insertNewGeom(GAM106geom) ### 38GAM0108 GAM108 = cor_dt_new[cor_dt_new.id_origine=='38GAM0108'].copy() GAM108geom = (GAM108[['site_code','id_origine','date','auteur','geom']] @@ -729,6 +747,64 @@ if __name__ == "__main__": ############################# ######## New maj ############ + ### 38GAM0107 + replace_geom = maj_dt_old[maj_dt_old.id_origine=='38GAM0107'].copy() + GAM0107_geom = replace_geom[['site_code','id_origine','date','geom','auteur']] + GAM0107_geom.date = '2017-01-01' + insertNewGeom((GAM0107_geom + .merge(info,left_on='id_origine',right_on='CODE_SITE',how='left') + .drop(columns=['CODE_SITE','nom'])) + ) + + rpl_delim = (GAM0107_geom + .merge(delim,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + hactivhum = gpd.pd.read_sql('SELECT * FROM zones_humides.param_activ_hum',con=pycen.con) + rpl_activHum = (GAM0107_geom + .merge(usg,left_on='id_origine',right_on='CODE_SITE',how='inner') + .rename(columns={ + 'crit_usg':'activite_hum', + 'remarques':'remarks'}) + .astype({'activite_hum':int}) + .replace({'activite_hum':dict(zip(hactivhum.id,hactivhum.nom))}) + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_activHum['impact'] = 'Inconnu' + rpl_activHum['localisation'] = 'Inconnu' + rpl_RHyd = (GAM0107_geom + .merge(reghydro,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_subm = (GAM0107_geom + .merge(subm,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_cnx = (GAM0107_geom + .merge(cnx,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_fct = (GAM0107_geom + .merge(fct,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom']) + .rename(columns={ 'description':'memo'})) + insertAttrsDelim(rpl_delim) + insertAttrsUsgPrss(rpl_activHum) + insert_regHydro(rpl_RHyd) + insert_subm(rpl_subm) + insert_cnx(rpl_cnx) + insertAttrsFct(rpl_fct,True) + ## ADD GEOM 2023 + PATH = '/media/colas/SRV/FICHIERS/OUTILS/CARTOGRAPHIE/ESPACE DE TRAVAIL/ETUDES/PLUI METRO/INV_ZH_PLUI_METRO/' + FILE_2023 = 'zh_gam_cc45_modifByMJ&GC.gpkg' + df2023 = gpd.read_file(path.join(PATH,FILE_2023)) + df2023.rename_geometry('geom',inplace=True) + geom2023 = (df2023 + .loc[df2023.CODE_SITE=='38GAM0107',['CODE_SITE','DATEMODIF','geom']] + .rename(columns={ + 'CODE_SITE':'id_site', 'DATEMODIF':'date', + })) + geom2023.date = '2023-11-28' + geom2023.id_site = '38VE0213' + geom2023['auteur'] = 'JUTON Mathieu' + insertNewGeom(geom2023) + + ### 38GAM0115 replace_geom = maj_dt_old[maj_dt_old.id_origine=='38GAM0115'].copy() replace_geom = replace_geom[['site_code','id_origine','date','geom','auteur']] @@ -773,6 +849,65 @@ if __name__ == "__main__": insert_cnx(rpl_cnx) insertAttrsFct(rpl_fct,True) + ### 38GAM0128 + replace_geom = cor_dt_new[cor_dt_new.id_origine=='38GAM0128'].copy() + GAM128_geom = (replace_geom[['site_code','id_origine','date','geom','auteur']] + .set_geometry('geom',crs=2154)) + + insertNewGeom((GAM128_geom + .merge(info,left_on='id_origine',right_on='CODE_SITE',how='left') + .drop(columns=['CODE_SITE','nom'])) + ) + + rpl_delim = (GAM128_geom + .merge(delim,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + hactivhum = gpd.pd.read_sql('SELECT * FROM zones_humides.param_activ_hum',con=pycen.con) + rpl_activHum = (GAM128_geom + .merge(usg,left_on='id_origine',right_on='CODE_SITE',how='inner') + .rename(columns={ + 'crit_usg':'activite_hum', + 'remarques':'remarks'}) + .astype({'activite_hum':int}) + .replace({'activite_hum':dict(zip(hactivhum.id,hactivhum.nom))}) + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_activHum['impact'] = 'Inconnu' + rpl_activHum['localisation'] = 'Inconnu' + rpl_RHyd = (GAM128_geom + .merge(reghydro,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_subm = (GAM128_geom + .merge(subm,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_cnx = (GAM128_geom + .merge(cnx,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom'])) + rpl_fct = (GAM128_geom + .merge(fct,left_on='id_origine',right_on='CODE_SITE',how='inner') + .drop(columns=['CODE_SITE','id_origine','geom']) + .rename(columns={ 'description':'memo'})) + insertAttrsDelim(rpl_delim) + insertAttrsUsgPrss(rpl_activHum) + insert_regHydro(rpl_RHyd) + insert_subm(rpl_subm) + insert_cnx(rpl_cnx) + insertAttrsFct(rpl_fct,True) + + ## ADD GEOM 2023 + PATH = '/media/colas/SRV/FICHIERS/OUTILS/CARTOGRAPHIE/ESPACE DE TRAVAIL/ETUDES/PLUI METRO/INV_ZH_PLUI_METRO/' + FILE_2023 = 'zh_gam_cc45_modifByMJ&GC.gpkg' + df2023 = gpd.read_file(path.join(PATH,FILE_2023)) + df2023.rename_geometry('geom',inplace=True) + geom2023 = (df2023 + .loc[df2023.CODE_SITE=='38GAM0128',['CODE_SITE','DATEMODIF','geom']] + .rename(columns={ + 'CODE_SITE':'id_site', 'DATEMODIF':'date', + })) + geom2023.date = '2023-11-28' + geom2023.id_site = '38GL0011' + geom2023['auteur'] = 'JUTON Mathieu' + insertNewGeom(geom2023) + ### 38GAM0091 replace_geom = cor_dt_old[cor_dt_old.id_origine=='38GAM0091'].copy() replace_geom = replace_geom[['site_code','id_origine','date','geom','auteur']] @@ -1023,3 +1158,33 @@ if __name__ == "__main__": insert_cnx(COR_maj_newCnx) insertAttrsFct(COR_maj_newFct,True) + + ########################## + ##### ADD 2010 GEOM ###### + ########################## + PATH = '/media/colas/SRV/FICHIERS/OUTILS/CARTOGRAPHIE/VECTEURS/ETUDES/ZONES HUMIDES/INVENTAIRE_ZH/' + FILE_2010 = 'ZH_2010 (re-save).shp' + df2010 = gpd.read_file(PATH+FILE_2010) + df2010.to_crs(2154,inplace=True) + df2010.rename_geometry('geom',inplace=True) + GL11 = (df2010.loc[df2010.SITE_CODE=='38GL0011',['SITE_CODE','geom','DATEMODIF']] + .copy() + .rename(columns={ + 'DATEMODIF':'date','SITE_CODE':'id_site' + })) + GL11['auteur'] = 'CEN Isère' + insertNewGeom(pycen.tools.Polygons_to_MultiPolygon(GL11)) + GL12 = (df2010.loc[df2010.SITE_CODE=='38GL0012',['SITE_CODE','geom','DATEMODIF']] + .copy() + .rename(columns={ + 'DATEMODIF':'date','SITE_CODE':'id_site' + })) + GL12['auteur'] = 'CEN Isère' + insertNewGeom(pycen.tools.Polygons_to_MultiPolygon(GL12)) + VE213 = (df2010.loc[df2010.SITE_CODE=='38VE0213',['SITE_CODE','geom','DATEMODIF']] + .copy() + .rename(columns={ + 'DATEMODIF':'date','SITE_CODE':'id_site' + })) + VE213['auteur'] = 'CEN Isère' + insertNewGeom(pycen.tools.Polygons_to_MultiPolygon(VE213)) \ No newline at end of file