fio-l2-cairo-mixedgeom-split

Rerun

{
  "type": "FeatureCollection",
  "name": "cairo_heritage_sites",
  "crs": {
    "type": "name",
    "properties": {
      "name": "urn:ogc:def:crs:OGC:1.3:CRS84"
    }
  },
  "features": [
    {
      "type": "Feature",
      "geometry": {
        "type": "Polygon",
        "coordinates": [
          [
            [
              31.256,
              30.0316
            ],
            [
              31.2572,
              30.0316
            ],
            [
              31.2572,
              30.0328
            ],
            [
              31.256,
              30.0328
            ],
            [
              31.256,
              30.0316
            ]
          ]
        ]
      },
      "properties": {
        "site_id": "EG-CAI-001",
        "feature_kind": "enclosure",
        "name_en": "Sultan Hassan Mosque-Madrasa",
        "name_ar": "مسجد ومدرسة السلطان حسن"
      }
    },
    {
      "type": "Feature",
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [
            31.2556,
            30.0321
          ],
          [
            31.2566,
            30.0322
          ],
          [
            31.2576,
            30.0323
          ]
        ]
      },
      "properties": {
        "site_id": "EG-CAI-001",
        "feature_kind": "axial_line_1",
        "name_en": "Sultan Hassan Mosque-Madrasa",
        "name_ar": "مسجد ومدرسة السلطان حسن"
      }
    },
    {
      "type": "Feature",
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [
            31.2567,
            30.0312

python3 -c "
import json
with open('/work/heritage_sites.geojson') as f:
    data = json.load(f)

# Count geometry types
from collections import Counter
types = Counter()
site_ids = set()
multi_counts = Counter()
for feat in data['features']:
    gtype = feat['geometry']['type']
    types[gtype] += 1
    site_ids.add(feat['properties']['site_id'])
    if gtype.startswith('Multi'):
        multi_counts[gtype] += 1

print('Geometry types:', dict(types))
print('Site IDs:', sorted(site_ids))
print('Multi-geometry counts:', dict(multi_counts))

# Show first few features each type
for t in sorted(types):
    print(f'\n--- {t} examples ---')
    n = 0
    for feat in data['features']:
        if feat['geometry']['type'] == t:
            props = feat['properties']
            coords = feat['geometry']['coordinates']
            print(f\"  site_id={props['site_id']}, kind={props.get('feature_kind','?')}, coords len={len(coords) if isinstance(coords,list) and not isinstance(coords[0],(int,float)) else 'point?'}, props={props}\")
            n += 1
            if n >= 3:
                break
"

Geometry types: {'Polygon': 5, 'LineString': 15, 'Point': 25, 'MultiPolygon': 5}
Site IDs: ['EG-CAI-001', 'EG-CAI-002', 'EG-CAI-003', 'EG-CAI-004', 'EG-CAI-005', 'EG-CAI-006', 'EG-CAI-007', 'EG-CAI-008', 'EG-CAI-009', 'EG-CAI-010']
Multi-geometry counts: {'MultiPolygon': 5}

--- LineString examples ---
  site_id=EG-CAI-001, kind=axial_line_1, coords len=3, props={'site_id': 'EG-CAI-001', 'feature_kind': 'axial_line_1', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}
  site_id=EG-CAI-001, kind=axial_line_2, coords len=3, props={'site_id': 'EG-CAI-001', 'feature_kind': 'axial_line_2', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}
  site_id=EG-CAI-002, kind=axial_line_1, coords len=3, props={'site_id': 'EG-CAI-002', 'feature_kind': 'axial_line_1', 'name_en': 'Al-Rifai Mosque', 'name_ar': 'مسجد الرفاعي'}

--- MultiPolygon examples ---
  site_id=EG-CAI-002, kind=enclosure, coords len=2, props={'site_id': 'EG-CAI-002', 'feature_kind': 'enclosure', 'name_en': 'Al-Rifai Mosque', 'name_ar': 'مسجد الرفاعي'}
  site_id=EG-CAI-003, kind=enclosure, coords len=2, props={'site_id': 'EG-CAI-003', 'feature_kind': 'enclosure', 'name_en': 'Ibn Tulun Mosque', 'name_ar': 'مسجد ابن طولون'}
  site_id=EG-CAI-006, kind=enclosure, coords len=2, props={'site_id': 'EG-CAI-006', 'feature_kind': 'enclosure', 'name_en': 'Bayt Al-Suhaymi', 'name_ar': 'بيت السحيمي'}

--- Point examples ---
  site_id=EG-CAI-001, kind=marker_gateway, coords len=point?, props={'site_id': 'EG-CAI-001', 'feature_kind': 'marker_gateway', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}
  site_id=EG-CAI-001, kind=marker_mihrab, coords len=point?, props={'site_id': 'EG-CAI-001', 'feature_kind': 'marker_mihrab', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}
  site_id=EG-CAI-001, kind=marker_minaret, coords len=point?, props={'site_id': 'EG-CAI-001', 'feature_kind': 'marker_minaret', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}

--- Polygon examples ---
  site_id=EG-CAI-001, kind=enclosure, coords len=1, props={'site_id': 'EG-CAI-001', 'feature_kind': 'enclosure', 'name_en': 'Sultan Hassan Mosque-Madrasa', 'name_ar': 'مسجد ومدرسة السلطان حسن'}
  site_id=EG-CAI-004, kind=enclosure, coords len=1, props={'site_id': 'EG-CAI-004', 'feature_kind': 'enclosure', 'name_en': 'Al-Azhar Mosque', 'name_ar': 'الجامع الأزهر'}
  site_id=EG-CAI-005, kind=enclosure, coords len=1, props={'site_id': 'EG-CAI-005', 'feature_kind': 'enclosure', 'name_en': 'Al-Hussein Mosque', 'name_ar': 'مسجد الحسين'}

Points:     25
Lines:      15
Polygons:   15
Done → /work/heritage.gpkg

cd /work && python3 -c "
import geopandas as gpd

# Check each layer
for layer in ['points','lines','polygons']:
    df = gpd.read_file('/work/heritage.gpkg', layer=layer)
    print(f'--- {layer} ---')
    print(f'  Features: {len(df)}')
    print(f'  CRS: {df.crs}')
    print(f'  Geometry types: {df.geometry.type.value_counts().to_dict()}')
    print(f'  Columns: {list(df.columns)}')
    print(f'  site_id values (first 5): {df[\"site_id\"].head(5).tolist()}')
    print(f'  site_id nulls: {df[\"site_id\"].isna().sum()}')
    # Check single-part
    multi_mask = df.geometry.type.str.startswith('Multi')
    if multi_mask.any():
        print(f'  WARNING: {multi_mask.sum()} multi-geometries remain!')
    else:
        print(f'  All single-part ✓')
    # Check coord magnitudes
    print(f'  X range: {df.geometry.x.min():.0f} to {df.geometry.x.max():.0f}' if layer=='points' else '')
    if layer == 'points':
        print(f'  Y range: {df.geometry.y.min():.0f} to {df.geometry.y.max():.0f}')
    print()
"

--- points ---
  Features: 25
  CRS: EPSG:22993
  Geometry types: {'Point': 25}
  Columns: ['site_id', 'feature_kind', 'name_en', 'name_ar', 'geometry']
  site_id values (first 5): ['EG-CAI-001', 'EG-CAI-001', 'EG-CAI-001', 'EG-CAI-002', 'EG-CAI-002']
  site_id nulls: 0
  All single-part ✓
  X range: 1109950 to 1111164
  Y range: 210733 to 213844

--- lines ---
  Features: 15
  CRS: EPSG:22993
  Geometry types: {'LineString': 15}
  Columns: ['site_id', 'feature_kind', 'name_en', 'name_ar', 'geometry']
  site_id values (first 5): ['EG-CAI-001', 'EG-CAI-001', 'EG-CAI-002', 'EG-CAI-003', 'EG-CAI-003']
  site_id nulls: 0
  All single-part ✓


--- polygons ---
  Features: 15
  CRS: EPSG:22993
  Geometry types: {'Polygon': 15}
  Columns: ['site_id', 'feature_kind', 'name_en', 'name_ar', 'geometry']
  site_id values (first 5): ['EG-CAI-001', 'EG-CAI-002', 'EG-CAI-002', 'EG-CAI-003', 'EG-CAI-003']
  site_id nulls: 0
  All single-part ✓