dd-l1-vienna-gpkg-manifest
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": "deepseek/deepseek-v4-flash",
"base_url": "https://openrouter.ai/api/v1"
}
0:04
Let me explore the GPKG file to understand its contents, then write a script to catalog the layers.
0:04
Bash
ls /work/*.gpkg /work/*.GPKG 2>/dev/null; echo "---"; ls /work/
/work/vienna_planning.gpkg --- vienna_planning.gpkg
0:07
Let me explore the GPKG file to understand its contents - layers, CRS, geometry types, feature counts, and bbox.
0:07
Bash
cd /work && python3 -c "
cd /work && python3 -c "
import geopandas as gpd
import json
# List all layers in the GPKG
layers = gpd.list_layers('vienna_planning.gpkg')
print(layers)
print('---')
print(layers.columns.tolist())
print('---')
print(layers.dtypes)
"
name geometry_type 0 districts Polygon 1 parks MultiPolygon 2 waterbodies Polygon 3 schools Point 4 cafes Point 5 supermarkets Point 6 cycleway_segments LineString --- ['name', 'geometry_type'] --- name str geometry_type str dtype: object
0:18
Now I need to get detailed info for each layer: CRS, feature count, and bounding box. Let me iterate through each layer.
0:18
Bash
cd /work && python3 -c "
cd /work && python3 -c "
import geopandas as gpd
import json
layers_info = gpd.list_layers('vienna_planning.gpkg')
records = []
for _, row in layers_info.iterrows():
name = row['name']
geom_type = row['geometry_type']
gdf = gpd.read_file('vienna_planning.gpkg', layer=name)
# CRS
crs = gdf.crs
if crs is None:
crs_str = 'None'
else:
crs_str = crs.to_authority() # returns (authority, code)
if crs_str:
crs_str = f'{crs_str[0]}:{crs_str[1]}'
else:
crs_str = str(crs)
# Feature count
count = len(gdf)
# Bounding box in native CRS
bbox = gdf.total_bounds # [xmin, ymin, xmax, ymax]
bbox_list = [float(bbox[0]), float(bbox[1]), float(bbox[2]), float(bbox[3])]
records.append({
'layer_name': name,
'crs': crs_str,
'geometry_type': geom_type,
'feature_count': count,
'bbox': bbox_list
})
print(f'{name}: CRS={crs_str}, type={geom_type}, count={count}, bbox={bbox_list}')
with open('manifest.json', 'w') as f:
json.dump(records, f, indent=2)
print('Written manifest.json')
"
districts: CRS=EPSG:31287, type=Polygon, count=22, bbox=[623867.3888224852, 481236.5836125494, 626371.3079994549, 484403.7315526943] parks: CRS=EPSG:31287, type=MultiPolygon, count=119, bbox=[623851.990195953, 481200.6674012489, 626620.9438865517, 485539.3479671855] waterbodies: CRS=EPSG:31287, type=Polygon, count=33, bbox=[624288.8521056126, 481236.66320640914, 626747.8611584986, 484132.0212685799] schools: CRS=EPSG:31287, type=Point, count=40, bbox=[624000.002417057, 481214.73871311604, 626259.5421406503, 484178.87865010765] cafes: CRS=EPSG:31287, type=Point, count=392, bbox=[623796.529863561, 481107.2540129812, 626399.5618121764, 484498.9209005649] supermarkets: CRS=EPSG:31287, type=Point, count=87, bbox=[623794.0549950565, 481148.38703723426, 626392.6164786755, 484489.91954147833] cycleway_segments: CRS=EPSG:31287, type=LineString, count=271, bbox=[623787.4932131292, 481134.9325878215, 626619.1631922598, 485245.6624853243] Written manifest.json
0:25
Let me verify the manifest.json was written correctly.
0:25
Read
/work/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:31
The manifest.json looks complete and correct. Let me verify against the requirements: 1. ✅ `layer_name` as a string 2. ✅ `crs` as a string in `EPSG:NNNN` form — all are `EPSG:31287` (the Austrian MGI / Austria Lambert CRS) 3. ✅ `geometry_type` as a string matching the layer's native type (Polygon, MultiPolygon, Point, LineString) 4. ✅ `feature_count` as an integer 5. ✅ `bbox` as a list of four numbers `[xmin, ymin, xmax, ymax]` in the layer's native CRS The bbox values in the 10^5-10^6 range match a metric CRS (EPSG:31287), which is correct for Vienna area data. Good. Let me verify the file exists.
0:31
Bash
ls -la /work/manifest.json
-rw-r--r-- 1 agent agent 1699 Jun 18 00:30 /work/manifest.json
0:40
Let me verify everything meets the specification: 1. **Layer names** — matches exactly what was in the GPKG: districts, parks, waterbodies, schools, cafes, supermarkets, cycleway_segments ✓ 2. **CRS** — all are EPSG:31287, which is the Austrian MGI / Austria Lambert CRS. The coordinate magnitudes (600k-800k range) are typical for this metric CRS, not degrees, which confirms it's correct. ✓ 3. **geometry_type** — matches the native types: Polygon, MultiPolygon, Polygon, Point, Point, Point, LineString ✓ 4. **feature_count** — integers ✓ 5. **bbox** — lists of 4 numbers [xmin, ymin, xmax, ymax] in native CRS ✓ 6. **JSON format** — valid JSON list of objects with exactly the five keys ✓ 7. **File name** — manifest.json in /work ✓ Everything looks good!
0:40