Skip to content

Instantly share code, notes, and snippets.

@mickaelandrieu
Created February 3, 2023 20:46
Show Gist options
  • Save mickaelandrieu/06b9dc5b91c4e82adbe1cbc2e6cba9f9 to your computer and use it in GitHub Desktop.
Save mickaelandrieu/06b9dc5b91c4e82adbe1cbc2e6cba9f9 to your computer and use it in GitHub Desktop.
import geopandas as gdp
dep971 = gdp.read_file('./ZFU/ZFU_Dep971_Scan25_W84.shp')
dep972 = gdp.read_file('./ZFU/ZFU_Dep972_Scan25_W84.shp')
dep973 = gdp.read_file('./ZFU/ZFU_Dep973_Scan25_W84.shp')
dep974 = gdp.read_file('./ZFU/ZFU_Dep974_Scan25_W84.shp')
metropole = gdp.read_file('./ZFU/ZFU_FRM_Scan25_L93.shp')
# convert metropole to lon/lat system
dep971 = dep971.to_crs('EPSG:4326')
dep972 = dep972.to_crs('EPSG:4326')
dep973 = dep973.to_crs('EPSG:4326')
dep974 = dep974.to_crs('EPSG:4326')
metropole = metropole.to_crs('EPSG:4326')
total_zfu = gdp.pd.concat([dep971, dep972, dep973, dep974, metropole])
total_zfu.to_file('total_zfu.shp')
total_zfu.to_file('zfus.geojson', driver='GeoJSON')
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment