import geopandas as gpd
from os import path
from pycen import con

PATH = '/home/colas/Documents/9_PROJETS/1_ZH/MAJ/GERECO - Schéma_strat_ZH_Huez/Inventaire_2016_GERECO'
F1 = 'Inf1000m2_NonSurfaciques.shp'
F2 = 'Inf1000m2_Surfaciques.shp'

df1 = gpd.read_file(path.join(PATH,F1))
del df1['id']
df1.rename_geometry('geom',inplace=True)
df2 = gpd.read_file(path.join(PATH,F2))
del df2['id']
del df2['SurfHa']
df2['geom'] = df2.centroid
df2.set_geometry('geom', crs=2154, inplace=True)
del df2['geometry']

df = (gpd.pd.concat([df1,df2])
      .rename(columns={'NOM_ZH':'comment'}))

df['qui'] = 'GERECKO'
df['source'] = "Bureau d'étude"
df['ponctuelle'] = 'Oui'
df['type'] = 'Inconnu'
df['classe'] = 'Inconnu'
df['x'] = df.geom.x
df['y'] = df.geom.y
df['datemodif'] = '2016-01-01'
df['observer'] = 'GERECKO'

df.to_postgis(
    'zh_ponctuelle',con,'zones_humides',if_exists='append', index=False
)