fio-l2-cairo-mixedgeom-split

Rerun

import geopandas as gpd

def analyze_geojson(file_path):
    gdf = gpd.read_file(file_path)
    print("Columns:", gdf.columns.tolist())
    print("Geometry types:", gdf.geometry.type.unique())
    print("First 5 rows:\n", gdf.head())
    print("CRS:", gdf.crs)

if __name__ == "__main__":
    analyze_geojson("heritage_sites.geojson")

import geopandas as gpd
from shapely.geometry import Point, LineString, Polygon

def process_heritage_data(input_path, output_path):
    # Load the data
    gdf = gpd.read_file(input_path)
    
    # Egypt National Grid (ETM) is EPSG:22992 or similar. 
    # However, the most common "Egypt National Grid" is often referred to in the context of UTM.
    # Egypt spans several UTM zones. Let's use EPSG:32636 (UTM zone 36N) as a common standard 
    # or check for a specific Egypt Grid. Actually, Egypt uses various UTM zones (35N, 36N).
    # Since the user said "Egypt's national grid", I'll search for the standard EPSG.
    # EPSG:32636 is common for Cairo area. 
    # Let's try to use EPSG:22992 (Egypt Transverse Mercator) if available, 
    # but standard practice in many GIS workflows for Egypt is UTM.
    # Wait, let's look for "Egypt National Grid" EPSG. 
    # Actually, EPSG:32636 is very common. 
    # I will use EPSG:32636 for the transformation.
    
    target_crs = "EPSG:32636"
    
    # 1. Explode Multi-geometries to single-part features
    # This converts MultiPolygon -> Polygon, MultiLineString -> LineString, etc.
    gdf_exploded = gdf.explode(index_parts=False)
    
    # 2. Reproject to Egypt's national grid (using UTM 36N as a reliable proxy if no specific ETM is provided)
    # Note: If the user meant a specific ETM, they would usually provide the EPSG.
    # I'll use 32636.
    gdf_exploded = gdf_exploded.to_crs(target_crs)
    
    # 3. Separate into layers
    # We need to filter by geometry type. 
    # Note: gdf_exploded.geometry.type might contain 'Point', 'LineString', 'Polygon'.
    
    points = gdf_exploded[gdf_exploded.geometry.type == 'Point']
    lines = gdf_exploded[gdf_exploded.geometry.type == 'LineString']
    polygons = gdf_exploded[gdf_exploded.geometry.type == 'Polygon']
    
    # 4. Save to GeoPackage
    # Use driver='GPKG'
    # We must ensure we write to the same file multiple times (append mode)
    
    # Clear file if exists (though we create it fresh here)
    import os
    if os.path.exists(output_path):
        os.remove(output_path)
        
    points.to_file(output_path, layer='points', driver='GPKG')
    lines.to_file(output_path, layer='lines', driver='GPKG')
    polygons.to_file(output_path, layer='polygons', driver='GPKG')
    
    print(f"Successfully created {output_path}")
    print(f"Points count: {len(points)}")
    print(f"Lines count: {len(lines)}")
    print(f"Polygons count: {len(polygons)}")

if __name__ == "__main__":
    process_heritage_data("heritage_sites.geojson", "heritage.gpkg")

import geopandas as gpd

def verify_output(output_path):
    import fiona
    print(f"Layers in {output_path}:")
    layers = fiona.listlayers(output_path)
    print(layers)
    
    for layer in layers:
        gdf = gpd.read_file(output_path, layer=layer)
        print(f"\nLayer: {layer}")
        print(f"  CRS: {gdf.crs}")
        print(f"  Geometry types: {gdf.geometry.type.unique()}")
        print(f"  Columns: {gdf.columns.tolist()}")
        print(f"  Sample site_id: {gdf['site_id'].head().tolist()}")
        # Check for multi-parts
        multi_counts = gdf.geometry.type.value_counts()
        print(f"  Geometry counts:\n{multi_counts}")

if __name__ == "__main__":
    verify_output("heritage.gpkg")