dd-l1-vienna-gpkg-manifest

Rerun

# /// script
# requires-python = ">=3.12"
# dependencies = ["geopandas", "fiona", "pyproj"]
# ///

import geopandas as gpd
import fiona
import json

def get_epsg_string(crs):
    if crs is None:
        return "Unknown"
    # Handle if it's a pyproj object or a string
    try:
        # Try to get EPSG code from CRS object
        return f"EPSG:{crs.to_epsg()}"
    except Exception:
        # If it's already a string like 'EPSG:31287'
        if isinstance(crs, str) and crs.startswith("EPSG:"):
            return crs
        # Fallback for other formats if possible
        return str(crs)

def generate_manifest(gpkg_path, output_path):
    layers_info = []
    
    # Use fiona to list layers and metadata without loading all data
    layers = fiona.listlayers(gpkg_path)
    
    for layer_name in layers:
        with fiona.open(gpkg_path, layer=layer_name) as src:
            # Get CRS
            crs_obj = src.crs
            epsg_str = get_epsg_string(crs_obj)
            
            # Get geometry type
            # fiona geometry types are usually 'Polygon', 'LineString', etc.
            # We'll check the first feature or the schema
            geom_type = src.schema['geometry']
            # Note: fiona/geopandas might return 'MultiPolygon' or 'Polygon'
            
            # Get feature count
            feature_count = len(src)
            
            # Get bbox [xmin, ymin, xmax, ymax]
            # fiona provides bounds: (minx, miny, maxx, maxy)
            bbox = list(src.bounds)
            
            layers_info.append({
                "layer_name": layer_name,
                "crs": epsg_str,
                "geometry_type": geom_type,
                "feature_count": feature_count,
                "bbox": bbox
            })

    with open(output_path, 'w') as f:
        json.dump(layers_info, f, indent=4)

if __name__ == "__main__":
    generate_manifest("vienna_planning.gpkg", "manifest.json")

# /// script
# requires-python = ">=3.12"
# dependencies = ["geopandas", "fiona", "pyproj"]
# ///

import geopandas as gpd
import fiona
import json

def get_epsg_string(crs):
    if crs is None:
        return "Unknown"
    # Fiona CRS is often a dict or a tuple. 
    # If it is a dict, we might need to extract it.
    # If it's a pyproj object or a CRS object, we use to_epsg.
    
    # Let's try a more robust way to get the EPSG code
    try:
        # If it's a dictionary (common in Fiona)
        if isinstance(crs, dict):
            # Check if 'init' or 'name' or something similar is there
            # Or if it's a proj string
            import pyproj
            p = pyproj.CRS(crs)
            epsg = p.to_epsg()
            if epsg:
                return f"EPSG:{epsg}"
            else:
                return str(crs)
        
        # If it's a string
        if isinstance(crs, str):
            return crs
            
        # Try pyproj directly
        import pyproj
        p = pyproj.CRS(crs)
        epsg = p.to_epsg()
        if epsg:
            return f"EPSG:{epsg}"
        else:
            return str(crs)
    except Exception:
        return str(crs)

def generate_manifest(gpkg_path, output_path):
    layers_info = []
    
    layers = fiona.listlayers(gpkg_path)
    
    for layer_name in layers:
        with fiona.open(gpkg_path, layer=layer_name) as src:
            # Get CRS
            crs_obj = src.crs
            epsg_str = get_epsg_string(crs_obj)
            
            # Get geometry type
            # fiona schema geometry type
            geom_type = src.schema['geometry']
            
            # Get feature count
            feature_count = len(src)
            
            # Get bbox [xmin, ymin, xmax, ymax]
            bbox = list(src.bounds)
            
            layers_info.append({
                "layer_name": layer_name,
                "crs": epsg_str,
                "geometry_type": geom_type,
                "feature_count": feature_count,
                "bbox": bbox
            })

    with open(output_path, 'w') as f:
        json.dump(layers_info, f, indent=4)

if __name__ == "__main__":
    generate_manifest("vienna_planning.gpkg", "manifest.json")