import geopandas as gpd
from shapely.geometry import Point

## input shapefile path
in_shp = r"path\to\input\shapefile\input.shp"

## the output shapefile path for the mean center point
out_shp = r"path\to\output\shapefile\output.shp"

## read in the shapefile to a GeoDataFrame
gdf = gpd.read_file(in_shp)

## get the geometry type from the first record
geom_type = gdf.geom_type[0]

## get the EPSG code
crs = gdf.crs

## for Point geometry
if geom_type == "Point":
    ## get all x and y values
    gdf["x"] = gdf.geometry.x
    gdf["y"] = gdf.geometry.y

## for LineString geometry
elif geom_type == "LineString":
    ## get all x and y values for the midpoints
    gdf["midpoint"] = gdf.geometry.interpolate(0.5, normalized=True)
    gdf["x"] = gdf["midpoint"].x
    gdf["y"] = gdf["midpoint"].y

## for Polygon geometry
elif geom_type == "Polygon":
    ## get all x and y values for the centroids
    gdf["centroid"] = gdf.geometry.centroid
    gdf["x"] = gdf["centroid"].x
    gdf["y"] = gdf["centroid"].y

## get the mean of the x and y values
mean_x = gdf['x'].mean()
mean_y = gdf['y'].mean()

## create a point geometry representing the mean center
mean_center = Point(mean_x, mean_y)

## create a GeoDataFrame with the point geometry
gdf_mean_center = gpd.GeoDataFrame(geometry=[mean_center], crs=crs)

## add fields for the XCoord and YCoord
gdf_mean_center["XCoord"] = mean_x
gdf_mean_center["YCoord"] = mean_y

## write the mean center point to the output shapefile
gdf_mean_center.to_file(out_shp, driver="ESRI Shapefile")