fio-l1-vienna-shapefile-recovery
Map (reference ← swipe → agent)
0:00
I'm migrating an old parcel shapefile (`parcels.shp`) into our web viewer, and a reference file `column_map.csv` is provided alongside it. Can you produce a clean GeoJSON `parcels.geojson` with the correct full column names and proper character encoding? The features should come through as Polygon or MultiPolygon.
0:00
System info
{
"adapter": "openrouter",
"model": "google/gemma-4-26b-a4b-it",
"base_url": "https://openrouter.ai/api/v1"
}
0:03
Bash
ls -F
column_map.csv parcels.cpg parcels.dbf parcels.prj parcels.shp parcels.shx
0:04
Read
column_map.csv
1 truncated,original
2 KATASTRALG,KATASTRALGEMEINDE_NAME
3 GRUNDSTUEC,GRUNDSTUECKSNUMMER
4 EIGENTUEME,EIGENTUEMER_NAME
5 WIDMUNG_BE,WIDMUNG_BEZEICHNUNG
6 STRASSE_NA,STRASSE_NAME
7 FLAECHE_M2,FLAECHE_M2
0:25
Write
/work/solve.py
content · 98 lines · py
import geopandas as gpd
import pandas as pd
import json
def main():
# Load the column mapping
# The CSV seems to have a leading space and some potential issues with header/content
# Based on the Read output:
# 1 truncated,original
# 2 KATASTRALG,KATASTRALGEMEINDE_NAME
# ...
# It looks like it might be a CSV where the first line is a header.
# Let's try reading it with pandas.
mapping_df = pd.read_csv('column_map.csv')
# Rename columns if necessary, but based on 'Read', it looks like:
# index, truncated, original
# Create a dictionary for mapping: {truncated_name: original_name}
# We need to make sure we skip the first row if it's just a header.
# Looking at the 'Read' output:
# 1 truncated,original
# This means the first line is "truncated,original".
# Wait, the line numbers in 'Read' are 1-indexed.
# line 1: truncated,original
# line 2: KATASTRALG,KATASTRALGEMEINDE_NAME
# It seems the file has no index column in the CSV itself,
# but the 'Read' tool added line numbers.
# Let's inspect the first few lines properly to be sure about the structure.
mapping_dict = {}
for _, row in mapping_df.iterrows():
# Based on 'Read' output:
# line 1: truncated,original (This is the header)
# line 2: KATASTRALG,KATASTRALGEMEINDE_NAME
# So column 0 is 'truncated', column 1 is 'original'
# But wait, the first line of Read says "1 truncated,original".
# If 'truncated,original' is the first line, then column 0 is 'truncated' and column 1 is 'original'.
# However, if the file actually starts with 1 truncated,original
# it might mean there's an extra column or it's just how 'Read' presents it.
# Let's assume the CSV has columns 'truncated' and 'original'.
if pd.notna(row['truncated']) and pd.notna(row['original']):
# The first row might be the header if it contains 'truncated'
if row['truncated'] == 'truncated':
continue
mapping_dict[str(row['truncated'])] = str(row['original'])
print(f"Mapping dictionary: {mapping_dict}")
# Load the shapefile
gdf = gpd.read_file('parcels.shp')
print("Shapefile loaded.")
print(f"Original columns: {gdf.columns.tolist()}")
# Rename columns based on the mapping
# We want the 'original' names (the values in the dictionary)
# The mapping is {truncated_name: original_name}
# Create a new dataframe with renamed columns
# We only keep columns that are in the mapping or are the geometry column
new_columns = {}
for col in gdf.columns:
if col in mapping_dict:
new_columns[col] = mapping_dict[col]
elif col == 'geometry':
new_columns[col] = 'geometry'
else:
# If it's not in the mapping, should we keep it?
# The prompt says "with the correct full column names".
# Usually this implies using the mapping for what's provided.
# Let's keep other columns as they are, or drop them if "clean" implies only mapped ones.
# "produce a clean GeoJSON ... with the correct full column names"
# suggests we should use the mapping.
new_columns[col] = col
gdf = gdf.rename(columns=new_columns)
# Ensure features are Polygon or MultiPolygon
# Geopandas usually handles this, but let's check.
# If there are any geometries that are not Polygons/MultiPolygons, we might want to filter them.
# But the prompt says "The features should come through as Polygon or MultiPolygon".
# This might imply filtering or ensuring they are valid.
# Filtering for only Polygon and MultiPolygon
gdf = gdf[gdf.geometry.type.isin(['Polygon', 'MultiPolygon'])]
print(f"Columns after renaming: {gdf.columns.tolist()}")
print(f"Number of features: {len(gdf)}")
# Save to GeoJSON
# GeoJSON supports UTF-8 by default.
gdf.to_file('parcels.geojson', driver='GeoJSON')
print("Successfully saved to parcels.geojson")
if __name__ == "__main__":
main()
Wrote /work/solve.py (4174 chars).
0:26
Bash
python solve.py
Mapping dictionary: {'KATASTRALG': 'KATASTRALGEMEINDE_NAME', 'GRUNDSTUEC': 'GRUNDSTUECKSNUMMER', 'EIGENTUEME': 'EIGENTUEMER_NAME', 'WIDMUNG_BE': 'WIDMUNG_BEZEICHNUNG', 'STRASSE_NA': 'STRASSE_NAME', 'FLAECHE_M2': 'FLAECHE_M2'}
Shapefile loaded.
Original columns: ['KATASTRALG', 'GRUNDSTUEC', 'EIGENTUEME', 'WIDMUNG_BE', 'STRASSE_NA', 'FLAECHE_M2', 'geometry']
Columns after renaming: ['KATASTRALGEMEINDE_NAME', 'GRUNDSTUECKSNUMMER', 'EIGENTUEMER_NAME', 'WIDMUNG_BEZEICHNUNG', 'STRASSE_NAME', 'FLAECHE_M2', 'geometry']
Number of features: 60
Successfully saved to parcels.geojson
0:36
Bash
head -n 20 parcels.geojson
{
"type": "FeatureCollection",
"name": "parcels",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31287" } },
"features": [
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Innere Stadt", "GRUNDSTUECKSNUMMER": "0013/1", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Mariahilfer Straße", "FLAECHE_M2": 750.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625700.0, 483400.0 ], [ 625700.0, 483425.0 ], [ 625730.0, 483425.0 ], [ 625730.0, 483400.0 ], [ 625700.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Mariahilf", "GRUNDSTUECKSNUMMER": "0020/2", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Döblinger Hauptstraße", "FLAECHE_M2": 751.25 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625732.0, 483400.0 ], [ 625732.0, 483425.0 ], [ 625762.0, 483425.0 ], [ 625762.0, 483400.0 ], [ 625732.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Währing", "GRUNDSTUECKSNUMMER": "0027/3", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Währinger Gürtel", "FLAECHE_M2": 752.5 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625764.0, 483400.0 ], [ 625764.0, 483425.0 ], [ 625794.0, 483425.0 ], [ 625794.0, 483400.0 ], [ 625764.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Döbling", "GRUNDSTUECKSNUMMER": "0034/4", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Höfergasse", "FLAECHE_M2": 753.75 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625796.0, 483400.0 ], [ 625796.0, 483425.0 ], [ 625826.0, 483425.0 ], [ 625826.0, 483400.0 ], [ 625796.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Hütteldorf", "GRUNDSTUECKSNUMMER": "0041/5", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Schönbrunner Allee", "FLAECHE_M2": 755.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625828.0, 483400.0 ], [ 625828.0, 483425.0 ], [ 625858.0, 483425.0 ], [ 625858.0, 483400.0 ], [ 625828.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Floridsdorf", "GRUNDSTUECKSNUMMER": "0048/6", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Bäckerstraße", "FLAECHE_M2": 750.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625860.0, 483400.0 ], [ 625860.0, 483425.0 ], [ 625890.0, 483425.0 ], [ 625890.0, 483400.0 ], [ 625860.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Brigittenau", "GRUNDSTUECKSNUMMER": "0055/7", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Naschmarkt", "FLAECHE_M2": 751.25 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625892.0, 483400.0 ], [ 625892.0, 483425.0 ], [ 625922.0, 483425.0 ], [ 625922.0, 483400.0 ], [ 625892.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Schönbrunn", "GRUNDSTUECKSNUMMER": "0062/8", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Mariahilfer Straße", "FLAECHE_M2": 752.5 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625924.0, 483400.0 ], [ 625924.0, 483425.0 ], [ 625954.0, 483425.0 ], [ 625954.0, 483400.0 ], [ 625924.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Innere Stadt", "GRUNDSTUECKSNUMMER": "0069/9", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Döblinger Hauptstraße", "FLAECHE_M2": 753.75 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625956.0, 483400.0 ], [ 625956.0, 483425.0 ], [ 625986.0, 483425.0 ], [ 625986.0, 483400.0 ], [ 625956.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Mariahilf", "GRUNDSTUECKSNUMMER": "0076/1", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Währinger Gürtel", "FLAECHE_M2": 755.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625988.0, 483400.0 ], [ 625988.0, 483425.0 ], [ 626018.0, 483425.0 ], [ 626018.0, 483400.0 ], [ 625988.0, 483400.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Währing", "GRUNDSTUECKSNUMMER": "0083/2", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Höfergasse", "FLAECHE_M2": 750.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625700.5, 483427.0 ], [ 625700.5, 483452.0 ], [ 625730.5, 483452.0 ], [ 625730.5, 483427.0 ], [ 625700.5, 483427.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Döbling", "GRUNDSTUECKSNUMMER": "0090/3", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Schönbrunner Allee", "FLAECHE_M2": 751.25 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625732.5, 483427.0 ], [ 625732.5, 483452.0 ], [ 625762.5, 483452.0 ], [ 625762.5, 483427.0 ], [ 625732.5, 483427.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Hütteldorf", "GRUNDSTUECKSNUMMER": "0097/4", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Bäckerstraße", "FLAECHE_M2": 752.5 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625764.5, 483427.0 ], [ 625764.5, 483452.0 ], [ 625794.5, 483452.0 ], [ 625794.5, 483427.0 ], [ 625764.5, 483427.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Floridsdorf", "GRUNDSTUECKSNUMMER": "0104/5", "EIGENTUEMER_NAME": "Schönbrunner Bauges. m.b.H.", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Naschmarkt", "FLAECHE_M2": 753.75 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625796.5, 483427.0 ], [ 625796.5, 483452.0 ], [ 625826.5, 483452.0 ], [ 625826.5, 483427.0 ], [ 625796.5, 483427.0 ] ] ] } },
{ "type": "Feature", "properties": { "KATASTRALGEMEINDE_NAME": "Brigittenau", "GRUNDSTUECKSNUMMER": "0111/6", "EIGENTUEMER_NAME": "Müller GmbH", "WIDMUNG_BEZEICHNUNG": "Wohngebiet", "STRASSE_NAME": "Mariahilfer Straße", "FLAECHE_M2": 755.0 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 625828.5, 483427.0 ], [ 625828.5, 483452.0 ], [ 625858.5, 483452.0 ], [ 625858.5, 483427.0 ], [ 625828.5, 483427.0 ] ] ] } },
0:55