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.
Any ideas?