35 lines
898 B
Python
35 lines
898 B
Python
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
|
|
)
|