fio-l2-cairo-mixedgeom-split

Rerun

import geopandas as gpd

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

# Print info for verification
print("Initial Data Info:")
print(gdf.info())
print("\nGeometry types present:")
print(gdf.geometry.type.value_counts())
print("\nFirst few rows:")
print(gdf.head())

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

def explode_to_single_parts(gdf):
    """
    Explodes multi-part geometries into single-part geometries.
    """
    # Explode multi-geometries into single ones
    # Note: explode() in recent geopandas works on multi-geometries
    exploded_gdf = gdf.explode(index_parts=True).reset_index(drop=True)
    return exploded_gdf

def main():
    # 1. Load the data
    gdf = gpd.read_file('heritage_sites.geojson')
    
    # 2. Reproject to Egypt National Grid (EPSG:22992 - WGS 84 / UTM zone 36N is common, 
    # but let's check if there is a specific 'Egypt National Grid'.
    # Egypt uses various UTM zones. However, EPSG:22992 is a common UTM zone for Egypt.
    # Actually, common Egyptian grids often refer to WGS 84 / UTM zone 36N (EPSG:32636).
    # Since the user asked for 'Egypt's national grid', I'll check for a standard one.
    # Many users mean UTM Zone 36N for Egypt. Let's use EPSG:32636.
    # Wait, a more precise 'Egypt National Grid' might be something else, 
    # but EPSG:32636 is standard for the Cairo area (where the sample seems to be).
    # I'll use EPSG:32636 as a reliable choice for Egypt.
    egypt_crs = "EPSG:32636"
    gdf = gdf.to_crs(egypt_crs)
    print(f"Reprojected to {egypt_crs}")

    # 3. Explode Multi-geometries to Single-part
    # The requirement is "every feature has to be single-part"
    gdf_single = explode_to_single_parts(gdf)
    
    # 4. Separate into layers: points, lines, polygons
    # We need to be careful about what counts as a Point/Line/Polygon
    # after exploding.
    
    points_gdf = gdf_single[gdf_single.geometry.type == 'Point'].copy()
    lines_gdf = gdf_single[gdf_single.geometry.type == 'LineString'].copy()
    polygons_gdf = gdf_single[gdf_single.geometry.type == 'Polygon'].copy()

    # 5. Write to GeoPackage
    output_file = 'heritage.gpkg'
    
    # Remove existing file if it exists to start clean
    import os
    if os.path.exists(output_file):
        os.remove(output_file)

    # Write each layer
    # Note: If a layer is empty, we might still want to write it or skip it.
    # The user asked for three layers named points, lines, and polygons.
    
    if not points_gdf.empty:
        points_gdf.to_file(output_file, layer='points', driver="GPKG")
        print(f"Wrote {len(points_gdf)} points to layer 'points'")
    else:
        print("No points found.")

    if not lines_gdf.empty:
        lines_gdf.to_file(output_file, layer='lines', driver="GPKG")
        print(f"Wrote {len(lines_gdf)} lines to layer 'lines'")
    else:
        print("No lines found.")

    if not polygons_gdf.empty:
        polygons_gdf.to_file(output_file, layer='polygons', driver="GPKG")
        print(f"Wrote {len(polygons_gdf)} polygons to layer 'polygons'")
    else:
        print("No polygons found.")

    print(f"Successfully created {output_file}")

if __name__ == "__main__":
    main()

import geopandas as gpd

def main():
    # Load the data
    gdf = gpd.read_file('heritage_sites.geojson')
    
    # Egypt National Grid / UTM zone 36N
    egypt_crs = "EPSG:32636"
    gdf = gdf.to_crs(egypt_crs)
    print(f"Reprojected to {egypt_crs}")

    # Explode Multi-geometries to Single-part
    # gdf.explode(index_parts=True) is standard in modern geopandas
    gdf_exploded = gdf.explode(index_parts=True).reset_index(drop=True)
    
    # Separate into layers based on geometry type
    # We filter based on the individual part type
    points_gdf = gdf_exploded[gdf_exploded.geometry.type == 'Point'].copy()
    lines_gdf = gdf_exploded[gdf_exploded.geometry.type == 'LineString'].copy()
    polygons_gdf = gdf_exploded[gdf_exploded.geometry.type == 'Polygon'].copy()

    # Output filename
    output_file = 'heritage.gpkg'

    # Write each layer to the same GPKG file
    # We use the GPKG driver to ensure multiple layers can be saved in one file
    points_gdf.to_file(output_file, layer='points', driver="GPKG")
    lines_gdf.to_file(output_file, layer='lines', driver="GPKG")
    polygons_gdf.to_file(output_file, layer='polygons', driver="GPKG")

    print(f"Created {output_file} with layers: points, lines, polygons.")
    print(f"Counts - Points: {len(points_gdf)}, Lines: {len(lines_gdf)}, Polygons: {len(polygons_gdf)}")

if __name__ == "__main__":
    main()

import geopandas as gpd

def main():
    # Load the data
    gdf = gpd.read_file('heritage_sites.geojson')
    
    # Egypt National Grid (UTM Zone 36N)
    egypt_crs = "EPSG:32636"
    gdf = gdf.to_crs(egypt_crs)
    print(f"Reprojected to {egypt_crs}")

    # Explode Multi-geometries to Single-part
    # index_parts=False ensures we don't get a multi-index that complicates the output
    gdf_exploded = gdf.explode(index_parts=False).reset_index(drop=True)
    
    # Separate into layers
    points_gdf = gdf_exploded[gdf_exploded.geometry.type == 'Point'].copy()
    lines_gdf = gdf_exploded[gdf_exploded.geometry.type == 'LineString'].copy()
    polygons_gdf = gdf_exploded[gdf_exploded.geometry.type == 'Polygon'].copy()

    output_file = 'heritage.gpkg'

    # Ensure we write to the same file. 
    # Note: to_file with driver='GPKG' overwrites by default if not careful, 
    # but layer-by-layer addition works in standard geopandas/fiona for GPKG.
    # However, the cleanest way to ensure a fresh file is to delete it if it exists.
    import os
    if os.path.exists(output_file):
        os.remove(output_file)

    points_gdf.to_file(output_file, layer='points', driver="GPKG")
    lines_gdf.to_file(output_file, layer='lines', driver="GPKG")
    polygons_gdf.to_file(output_file, layer='polygons', driver="GPKG")

    print(f"Created {output_file} with layers: points, lines, polygons.")
    print(f"Counts - Points: {len(points_gdf)}, Lines: {len(lines_gdf)}, Polygons: {len(polygons_gdf)}")
    
    # Verify site_id is still present
    for name, layer in [('points', points_gdf), ('lines', lines_gdf), ('polygons', polygons_gdf)]:
        if not layer.empty:
            print(f"Layer '{name}' site_id sample: {layer['site_id'].iloc[0]}")

if __name__ == "__main__":
    main()