import geopandas as gpd

# 1. Charger les IRIS
iris = gpd.read_file("iris.gpkg", layer="contours_iris")

# 2. Projeter en Lambert-93 (EPSG:2154) pour calculer les centroïdes
iris_projected = iris.to_crs(epsg=2154)

# 3. Calculer les centroïdes en Lambert-93
iris_projected['x_centroid'] = iris_projected.geometry.centroid.x
iris_projected['y_centroid'] = iris_projected.geometry.centroid.y

# 4. Reprojeter en WGS84 (EPSG:4326) pour obtenir latitude/longitude
iris_wgs84 = iris_projected.to_crs(epsg=4326)
iris_wgs84['longitude'] = iris_wgs84.geometry.centroid.x
iris_wgs84['latitude'] = iris_wgs84.geometry.centroid.y

# 5. Exporter les colonnes utiles
iris_wgs84[['code_insee', 'nom_iris', 'latitude', 'longitude']].to_csv(
    "/var/www/html/iris_france.csv",
    index=False
)