19 lines
671 B
Python
19 lines
671 B
Python
import pycen
|
|
import geopandas as gpd
|
|
|
|
par_id = '38412000ZD0136'
|
|
PATH_OUT = '/home/colas/Documents/tmp/'+par_id
|
|
|
|
sql = "SELECT par_id, geom FROM cadastre.parcelles_38 WHERE par_id = '%s'"%par_id
|
|
parc = gpd.read_postgis(sql,pycen.con_fon)
|
|
parc.to_crs(4326).to_file(PATH_OUT+'.geojson', driver='GeoJSON')
|
|
parc.to_crs(4326).to_file(PATH_OUT+'.gpkg', driver='GPKG')
|
|
|
|
sql = "SELECT * FROM zones_humides.v_zoneshumides WHERE ST_INTERSECTS(geom,'SRID=2154;%s')"%parc.geometry[0]
|
|
zh =gpd.read_postgis(sql,pycen.con)
|
|
|
|
|
|
sql = "SELECT * FROM saisie.vm_synthese_observations WHERE ST_INTERSECTS(geom,'SRID=2154;%s')"%parc.geometry[0]
|
|
obs_sicen =gpd.read_postgis(sql,pycen.con_sicen)
|
|
|