fio-l2-cairo-mixedgeom-split

Rerun

import geopandas as gpd

def solve():
    # Load the data
    gdf = gpd.read_file('heritage_sites.geojson')
    
    # Print info for debugging
    print("Initial Data Info:")
    print(gdf.info())
    print("\nGeometry types present:")
    print(gdf.geometry.type.value_counts())
    print("\nCRS:", gdf.crs)
    
    # 1. Convert to single-part geometries
    # explode() is the standard way to turn Multi-part into Single-part
    gdf = gdf.explode(index_parts=False)
    
    # 2. Project to Egypt's National Grid
    # Egypt National Grid is EPSG:22992 (UTM zone 36N) or similar? 
    # Actually, Egypt uses several. Let's check common ones or use a robust one like EPSG:32636 (UTM 36N).
    # However, "Egypt National Grid" often refers to a specific local projection.
    # Let's check if EPSG:22992 is correct or if I should use a standard UTM.
    # Looking up "Egypt National Grid EPSG": EPSG:22992 is often cited.
    # I'll use EPSG:32636 as a safe bet for UTM 36N if I'm unsure, but I'll try to find 22992.
    # Wait, EPSG:22992 is for Egypt (Egypt Grid). Let's try that.
    target_crs = "EPSG:22992"
    gdf = gdf.to_crs(target_crs)
    print(f"\nProjected to {target_crs}")
    
    # 3. Split into layers
    points = gdf[gdf.geometry.type == 'Point']
    lines = gdf[gdf.geometry.type.isin(['LineString', 'MultiLineString'])] # explode handled MultiLineString
    polygons = gdf[gdf.geometry.type.isin(['Polygon', 'MultiPolygon'])] # explode handled MultiPolygon
    
    # Re-filtering after explode to ensure we only have single parts
    # Since we exploded, they should already be single parts.
    # But let's be explicit about types.
    
    # Note: explode() on a MultiLineString gives LineString.
    # Let's refine the split logic.
    
    points = gdf[gdf.geometry.type == 'Point']
    lines = gdf[gdf.geometry.type == 'LineString']
    polygons = gdf[gdf.geometry.type == 'Polygon']
    
    # 4. Save to GPKG
    # Use pyogrio or fiona (geopandas uses them)
    output_file = 'heritage.gpkg'
    
    if not points.empty:
        points.to_file(output_file, layer='points', driver='GPKG')
    else:
        print("Warning: No points found.")
        
    if not lines.empty:
        lines.to_file(output_file, layer='lines', driver='GPKG')
    else:
        print("Warning: No lines found.")
        
    if not polygons.empty:
        polygons.to_file(output_file, layer='polygons', driver='GPKG')
    else:
        print("Warning: No polygons found.")

    print(f"\nSuccessfully saved to {output_file}")
    print("Layers created:")
    print(f"Points: {len(points)}")
    print(f"Lines: {len(lines)}")
    print(f"Polygons: {len(polygons)}")

if __name__ == "__main__":
    solve()