Verschiebung durch mangelnde Transformation der Projektion

[KiwiKilian]

Bei unserem Anwendungsfall dem Zuordnen von Standort zu Gemeinden mit den 95% simplified Daten ist mir aufgefallen, dass viele der Datensätze einen leichten Versatz haben. Und hiermit ist nicht die Unschärfe durch die Simplifizierung gemeint sondern eindeutige Versätze. Die kann man z.B. bei den kleinere Dateien auch hier auf GitHub (Beispiel Simplified 95 Bezirke Topo) visualisiert sehen:

Gleiches Ergebnis mit den Daten in QGIS:

Betrifft nach stichprobenartiger Kontrolle Gemeinde, Bezirke und Länder sowie auch die unterschiedlichen Jahre, besonders deutlich bei 95 % simplified.

Gemäß der Doku in der Readme habe ich den Workflow nachvollzogen Ich denke das Problem ist die Transformation der Projektion mit mapshaper.org – dieser Schritt ist leider nicht dokumentiert. In den GeoJSON Files steht, dass EPSG:4326 bzw. WGS84 (zwei Namen für das gleiche Koordinatensystem) genutzt wird. Ich habe daher die Quelldaten von Statistik Austria zuerst mit folgendem Befehl in der Console transformiert:

-proj wgs84

Dies führte zu gleichen Versätzen wie aktuell in den Daten. Beheben lässt sich das Problem mit der Angabe der initialen Projektion wie sie von Statistik Austria (EPSG:31287) angegeben wird:

-proj from=EPSG:31287 crs=EPSG:4326

Sollte ich mich jetzt nicht vertan haben, wär hier wohl eine Korrektur sinnvoll? Ggf. konnte an dieser Stelle der Workflow noch besser dokumentiert werden? Besteht vielleicht sogar eine Möglichkeit die Generierung automatisiert mit einem Script durchzuführen oder steckt hier insbesondere bei der Benennung viel händische Arbeit drin?
Vielen Dank für deine Arbeit in diesem Repository, ich konnte dank dem dokumentierten Workflow jetzt meinen eigenen Datensatz für meine Zwecke erstellen.

[jfiala]

Hallo, wir hatten das gleiche Problem und haben es einfach direkt im TopoJson post-gefixed:

Originalwert: (Longitude, Latitude)
"translate":[9.527906274165764,46.37302230607365]}

das problem mit dem alignment ist gelöst - habe die koordinaten translate im topojson angepasst:

Geändert auf: (Longitude, Latitude)
"translate":[9.530906274165764,46.37232230607365]}

Geänderte Werte:
Longitude nach rechts verschoben:
Diff: +0,003

Latitude nach unten verschoben:
Diff: -0,0007

46,37302230607365 - 46,37232230607365 = -0,0007

[jfiala]

vorher:

nachher:

[jfiala]

Das Ergebnis sieht man hier, wenn man bei der Adressuche nach einem Bezirk/Gemeinde sucht:
https://www.zerowasteaustria.at/zero-waste-shops.html

[KiwiKilian]

Danke für den Workaround @jfiala. Leichter wäre es dennoch, wenn die Daten korrekt wären 😄.