I'm merging a list of shapefiles with this code:
from pathlib import Path import pandas import geopandas from tqdm import tqdm folder = Path(r"read path") shapefiles = folder.glob("PARCELA(*).SHP") gdf = pandas.concat([ geopandas.read_file(shp) for shp in tqdm(shapefiles) ]).pipe(geopandas.GeoDataFrame) gdf.to_file(folder / r'write path')
The shapefile is created correctly, the problem is that I have diferents projections on some of the shapefiles i'm merging and I want to normalize them. What I thought is to add a column with the SRID of each geometry for later reproject each geometry to a unique SRID.
I know how to extract the EPSG of each geometry:
geom_srid_num = gdf.crs.to_epsg()
But I don't know how to add a new column for each row of the concatenation shown before.