Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
triple mapa
import geopandas as gp
import time
hdf=gp.read_file('Humedales', encoding='utf-8')
hdf.region.value_counts()
hdf.cut_reg.value_counts()
#hdf.set_crs(epsg=4326, inplace=True)
#hdf=hdf.to_crs(epsg=4326) # doesn not work for every region!
t0=time.time()
for reg, rdf in hdf.groupby('cut_reg'):
rrdf = rdf.to_crs(epsg=4326)
rrdf.to_file('Humedales_Region_%s.json' %reg, driver='GeoJSON')
print(reg, len(rdf), round(time.time()-t0,2))
@sergiolucero
Copy link
Author

sergiolucero commented Jun 8, 2022

Debieran estar limites comunales, rios y humedales. Más aeropuertos, etc:
https://www.bcn.cl/siit/mapas_vectoriales/index_html

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment