dd-l1-vienna-gpkg-manifest

Status: done
Score: 100
Duration: 0:26
Cost: 0.20¢
Model: google/gemma-4-26b-a4b-it

Map (reference ← swipe → agent)

0:00
Inherited a GPKG (vienna_planning) from a colleague who retired last week. Before I commit to scripting against it I need a one-page manifest so I can see what's in there and decide which layers actually feed next month's councillor briefing. Catalog the contents and give me one record per layer. Output manifest.json — a JSON list of these records. Each record in the list must be a JSON object with exactly five keys: `layer_name` as a string, `crs` as a string in `EPSG:NNNN` form, `geometry_type` as a string matching the layer's native type (e.g. `Point`, `LineString`, `Polygon`, `MultiPolygon`), `feature_count` as an integer, and `bbox` as a list of four numbers `[xmin, ymin, xmax, ymax]` in the layer's native CRS.
0:00
System info 
{
  "adapter": "openrouter",
  "model": "google/gemma-4-26b-a4b-it",
  "base_url": "https://openrouter.ai/api/v1"
}
0:00
Bash ls -F 
vienna_planning.gpkg
0:06
Write /work/solve.py
content · 59 lines · py 
import geopandas as gpd
import json

def catalog_gpkg(gpkg_path):
    import fiona
    
    layers = fiona.listlayers(gpkg_path)
    manifest = []
    
    for layer_name in layers:
        # Read only header/metadata to be efficient if possible, 
        # but geopandas is easy for getting bbox and count.
        # Using pyogrio or fiona directly is often faster for just metadata.
        
        gdf = gpd.read_file(gpkg_path, layer=layer_name)
        
        # Get CRS in EPSG:NNNN form
        epsg_code = gdf.crs.to_epsg()
        if epsg_code is None:
            # Try to get it from the wkt if to_epsg() fails
            # This is a fallback, though instructions ask for EPSG:NNNN
            crs_str = str(gdf.crs)
        else:
            crs_str = f"EPSG:{epsg_code}"
            
        # Get geometry type (native type of the first feature or from gdf.geometry.type)
        # gdf.geometry.type returns a series, we want the common type.
        # If it's mixed, we might have to be careful, but usually layers have one type.
        geom_type = gdf.geometry.type.iloc[0] if not gdf.empty else "None"
        
        # feature_count
        feature_count = len(gdf)
        
        # bbox [xmin, ymin, xmax, ymax]
        if not gdf.empty:
            bbox = list(gdf.total_bounds)
        else:
            bbox = [0.0, 0.0, 0.0, 0.0]
            
        manifest.append({
            "layer_name": layer_name,
            "crs": crs_str,
            "geometry_type": geom_type,
            "feature_count": feature_count,
            "bbox": bbox
        })
        
    return manifest

if __name__ == "__main__":
    gpkg_file = "vienna_planning.gpkg"
    try:
        manifest_data = catalog_gpkg(gpkg_file)
        with open("manifest.json", "w") as f:
            json.dump(manifest_data, f, indent=4)
        print("Successfully created manifest.json")
    except Exception as e:
        print(f"Error: {e}")

Wrote /work/solve.py (1942 chars).
0:07
Bash python solve.py 
Successfully created manifest.json
0:09
Read manifest.json 
     1	[
     2	    {
     3	        "layer_name": "districts",
     4	        "crs": "EPSG:31287",
     5	        "geometry_type": "Polygon",
     6	        "feature_count": 22,
     7	        "bbox": [
     8	            623867.3888224852,
     9	            481236.5836125494,
    10	            626371.3079994549,
    11	            484403.7315526943
    12	        ]
    13	    },
    14	    {
    15	        "layer_name": "parks",
    16	        "crs": "EPSG:31287",
    17	        "geometry_type": "MultiPolygon",
    18	        "feature_count": 119,
    19	        "bbox": [
    20	            623851.990195953,
    21	            481200.6674012489,
    22	            626620.9438865517,
    23	            485539.3479671855
    24	        ]
    25	    },
    26	    {
    27	        "layer_name": "waterbodies",
    28	        "crs": "EPSG:31287",
    29	        "geometry_type": "Polygon",
    30	        "feature_count": 33,
    31	        "bbox": [
    32	            624288.8521056126,
    33	            481236.66320640914,
    34	            626747.8611584986,
    35	            484132.0212685799
    36	        ]
    37	    },
    38	    {
    39	        "layer_name": "schools",
    40	        "crs": "EPSG:31287",
    41	        "geometry_type": "Point",
    42	        "feature_count": 40,
    43	        "bbox": [
    44	            624000.002417057,
    45	            481214.73871311604,
    46	            626259.5421406503,
    47	            484178.87865010765
    48	        ]
    49	    },
    50	    {
    51	        "layer_name": "cafes",
    52	        "crs": "EPSG:31287",
    53	        "geometry_type": "Point",
    54	        "feature_count": 392,
    55	        "bbox": [
    56	            623796.529863561,
    57	            481107.2540129812,
    58	            626399.5618121764,
    59	            484498.9209005649
    60	        ]
    61	    },
    62	    {
    63	        "layer_name": "supermarkets",
    64	        "crs": "EPSG:31287",
    65	        "geometry_type": "Point",
    66	        "feature_count": 87,
    67	        "bbox": [
    68	            623794.0549950565,
    69	            481148.38703723426,
    70	            626392.6164786755,
    71	            484489.91954147833
    72	        ]
    73	    },
    74	    {
    75	        "layer_name": "cycleway_segments",
    76	        "crs": "EPSG:31287",
    77	        "geometry_type": "LineString",
    78	        "feature_count": 271,
    79	        "bbox": [
    80	            623787.4932131292,
    81	            481134.9325878215,
    82	            626619.1631922598,
    83	            485245.6624853243
    84	        ]
    85	    }
    86	]
0:10
Bash ls -F 
manifest.json
solve.py
vienna_planning.gpkg
0:25