fio-l2-cairo-mixedgeom-split

Rerun

import geopandas as gpd

# Load the data
gdf = gpd.read_file('heritage_sites.geojson')

# Inspect the data
print("Initial data head:")
print(gdf.head())
print("\nGeometry types present:")
print(gdf.geometry.type.unique())
print("\nCRS:")
print(gdf.crs)

import geopandas as gpd
from shapely.geometry import MultiPoint, MultiLineString, MultiPolygon

def explode_to_single_part(gdf):
    """Explodes Multi-geometries into individual single-part geometries."""
    # Explode is available in modern geopandas to turn Multi-parts into individual rows
    # We use explode() which returns a GeoDataFrame with a MultiIndex.
    # We then reset the index to maintain a clean structure.
    exploded = gdf.explode(index_parts=True).reset_index(drop=True)
    return exploded

def process_heritage_data(input_file, output_file, target_crs):
    # Load the data
    gdf = gpd.read_file(input_file)
    
    # Transform to Egypt's National Grid (EPSG:22992 is often used for Egypt, 
    # but let's check common ones or if user meant something specific. 
    # WGS 84 / Egypt Red Belt is EPSG:22992. 
    # Actually, let's use EPSG:22992 if available or just stick to the standard for Egypt.
    # Let's use EPSG:22992 (Egypt Red Belt)
    gdf = gdf.to_crs(target_crs)
    
    # Explode Multi-geometries to single-part
    gdf = explode_to_single_part(gdf)
    
    # Separate by geometry type
    # We need to filter carefully because 'geometry.type' returns 'Point', 'LineString', 'Polygon'
    # Note: Multi-parts are already gone after explode.
    
    points = gdf[gdf.geometry.type == 'Point']
    lines = gdf[gdf.geometry.type == 'LineString']
    polygons = gdf[gdf.geometry.type == 'Polygon']
    
    # Save to GPKG
    # We write each layer to the same file
    with gpd.GeoDataFrame(columns=gdf.columns, geometry='geometry', crs=target_crs) as driver:
        # Use the engine for writing
        # To avoid conflict, we'll write them sequentially
        pass

    # Using a more robust way to write layers to GPKG
    points.to_file(output_file, layer='points', driver='GPKG')
    lines.to_file(output_file, layer='lines', driver='GPKG')
    polygons.to_file(output_file, layer='polygons', driver='GPKG')
    
    print(f"Successfully created {output_file}")
    print(f"Points count: {len(points)}")
    print(f"Lines count: {len(lines)}")
    print(f"Polygons count: {len(polygons)}")

if __name__ == "__main__":
    # Egypt National Grid: EPSG:22992 is standard for many parts, 
    # but let's verify if there's a more generic one or if EPSG:22992 is acceptable.
    # EPSG:22992 is Egypt Red Belt. 
    # Let's use EPSG:22992.
    process_heritage_data('heritage_sites.geojson', 'heritage.gpkg', 'EPSG:22992')

import geopandas as gpd

def explode_to_single_part(gdf):
    """Explodes Multi-geometries into individual single-part geometries."""
    # Explode is available in modern geopandas to turn Multi-parts into individual rows
    # We use explode(index_parts=False) to ensure we get a clean GeoDataFrame
    exploded = gdf.explode(index_parts=False).reset_index(drop=True)
    return exploded

def process_heritage_data(input_file, output_file, target_crs):
    # Load the data
    gdf = gpd.read_file(input_file)
    
    # Transform to Egypt's National Grid (EPSG:22992 is common)
    gdf = gdf.to_crs(target_crs)
    
    # Explode Multi-geometries to single-part
    gdf = explode_to_single_part(gdf)
    
    # Separate by geometry type
    points = gdf[gdf.geometry.type == 'Point'].copy()
    lines = gdf[gdf.geometry.type == 'LineString'].copy()
    polygons = gdf[gdf.geometry.type == 'Polygon'].copy()
    
    # Save to GPKG
    # We write each layer to the same file
    points.to_file(output_file, layer='points', driver='GPKG')
    lines.to_file(output_file, layer='lines', driver='GPKG')
    polygons.to_file(output_file, layer='polygons', driver='GPKG')
    
    print(f"Successfully created {output_file}")
    print(f"Points count: {len(points)}")
    print(f"Lines count: {len(lines)}")
    print(f"Polygons count: {len(polygons)}")

if __name__ == "__main__":
    # Egypt Red Belt (EPSG:22992)
    process_heritage_data('heritage_sites.geojson', 'heritage.gpkg', 'EPSG:22992')