Bearbeitung von Multipolygon GeoJson

Hypothalamus · 2019-01-14T09:36:23.000Z

Moin Simon,
das Format ist egal.
Mir würde reichen wenn in Excel jeweils in einer Spalte untereinander die LatLon Werte aus den eckigen Klammer stehen.
Nehmen wir mal an wir haben folgende Werte

[5.8663877,51.0511197],[5.8668165,51.0512156],[5.8669173,51.0512381],[5.866737,51.0515704]
Und in Excel in
A1 5.8663877,51.0511197
B1 5.8668165,51.0512156
C1 5.8669173,51.0512381

usw …

Natürlich wäre es auch toll wenn ich dieselben Daten als shp oder json vorliegen habe mit den Attributtabellen in denen die Daten stehen.

Gruß
Dennis

SimonPoole · 2019-01-14T09:44:37.000Z

Hypothalamus:

Moin Simon,
das Format ist egal.
Mir würde reichen wenn in Excel jeweils in einer Spalte untereinander die LatLon Werte aus den eckigen Klammer stehen.
Nehmen wir mal an wir haben folgende Werte

[5.8663877,51.0511197],[5.8668165,51.0512156],[5.8669173,51.0512381],[5.866737,51.0515704]
Und in Excel in
A1 5.8663877,51.0511197
B1 5.8668165,51.0512156
C1 5.8669173,51.0512381

usw …

Natürlich wäre es auch toll wenn ich dieselben Daten als shp oder json vorliegen habe mit den Attributtabellen in denen die Daten stehen.

Irgendwo besteht da ein fundamentales Verständnisproblem: die Listen von Koordinatenpaare in dem GeoJOSN Format definieren Umrisse von Objekten, die willst du also nicht? Sondern die Flächenschwerpunkte (da gibt es dann noch eine gewisse zusätzliche Schwierigkeit) der einzelnen Umrisse?

seichter · 2019-01-14T10:10:58.000Z

SimonPoole:

Irgendwo besteht da ein fundamentales Verständnisproblem: die Listen von Koordinatenpaare in dem GeoJOSN Format definieren Umrisse von Objekten, die willst du also nicht?

Das ist eher ein Kommunikationsproblem: Er will nach meinem Verständnis die Umrisse als Koordinatenpaare, dazu noch “irgendwie” als Fläche möglichst mit Attributen gruppiert. Als Umriss will er nur Polygone, keine MP, damit zwingend ohne Exklaven/Enklaven (“Inseln”?), also “wambacher light”.

Hypothalamus · 2019-01-14T10:21:06.000Z

Moin,
doch ich will die Umrisse also die Ländergrenzen.
Das Polygon ist doch das Objekt, in dem Falle das Land oder ?
Da ich die Landclass Grenzen haben also die Küstenlinien bestehen hier die Grenzen an den Inseln auch.
Deswegen wahrscheinlich das Multipolygon ? Da ein normals Polygon ja clockwise geschlossen sein muss.
Multipolygon wahrscheinlich aufgrund mehrfacher geschlossener Polygone.

Ich brauch lediglich die Lat/Lon Daten der Ländergrenze.

wambacher · 2019-01-14T10:42:22.000Z

Hypothalamus:

Das Problem ist das z.B. die Daten von wambacher Multipolygone sind. Klar aufgrund der Inseln usw.

Das Problem ist, dass “meine” Daten Administrative Grenzen darstellen, die nun mal per se Polygone sind.
Und Polygone können aus mehreren Flächen bestehen (Outer) und zudem “Löcher” haben (Inner).

In anderen Worten: Du verwendest die falsche Quelle. So als ob du zum Metzger gehst und dort Brot kaufen willst.

Die Sache ist mMn nicht trivial und würde auch für mich mehrere Tage Arbeit bedeuten. Mein Ansatzt wäre:

Installation von PostgreSQL/Postgis
Download der OSM-Daten für Deutschland von Geofabrik.de
Import mit osm2pgsql oder osmconvert in die Datenbank
Abfragen mit PostGIS

Der letzte Punkt dürfte am schwierigsten sein, da man dafür Erfahrung braucht - das gilt aber für jede Software.

Noch ein Tip: mit Shapefiles (shp) bist du besser bedient, wenn du es dennoch mit GIS-Software versuchen willst. Ganz am Ende kannst du das Ergebnis ja nach CSV konvertieren, was ich allerdings bei Geodaten nicht für sinnvoll halte. Von Excel und Windows ganz abgesehen

Gruss
walter

maxbe · 2019-01-14T10:58:31.000Z

Ich hab mir mal bei Naturalearth die Ländergrenzen (mittlere Auflösung, 1:50m) geholt und in QGis geladen. Da gibt es ein Werkzeug “Vektor” → “Geometriewerkzeuge” → “Mehr- zu Einteilig”. Das zerteilt Deutschland in 6 Stücke (viele Inseln sind in der Auflösdung nicht drauf). Das Ergebnis als JSON abgespeichert (das lag am nächsten zu einem Text…) sieht dann so aus (gezipptes geojson 1MB).

Sowas kann man sicher auch mit Wambachers Multipolygonen machen, man müsste noch sehen, wie diese “Mehr- zu Einteilig”-Funktion mit Enklaven umgeht…

Grüße
Max

wambacher · 2019-01-14T11:26:02.000Z

Hypothalamus:

Moin,
doch ich will die Umrisse also die Ländergrenzen.
Das Polygon ist doch das Objekt, in dem Falle das Land oder ?
Da ich die Landclass Grenzen haben also die Küstenlinien bestehen hier die Grenzen an den Inseln auch.
Deswegen wahrscheinlich das Multipolygon ? Da ein normals Polygon ja clockwise geschlossen sein muss.
Multipolygon wahrscheinlich aufgrund mehrfacher geschlossener Polygone.

Ich brauch lediglich die Lat/Lon Daten der Ländergrenze.
https://i.imgur.com/v9Zh4FP.jpg

sorry, hab ich erst nach meiner 1. Antwort gelesen.

Ja, das sind die Ländergrenzen. Geklippt an der Küstenlinie. Und da wirst du auch keine anderen Grenzen bekommen, da die nun mal so aussehen. Und diese sind IN JEDEM FORMAT nun mal Polygone mit Outern und ggf Innern.

Es kommt durchaus vor, dass Länder Exklaven (Outer) bzw Enklaven (Inner) haben können, ohne dass da Inseln involviert sind.

siehe Angola

oder Südafrika

Ansonsten hab ich das Gefühl, dass du selber noch nicht weisst, was du eigentlich willst. So ist z.B. die Frage, wie man Inseln benennen kann, völlig verschwunden. Und ohne diese Frage ist die Lösung wesentlich einfacher: Die Grenzen sind so, wie sie sind und du solltest die genau so nehmen.

Gruss
walter

SimonPoole · 2019-01-14T11:42:56.000Z

seichter:

SimonPoole:

Irgendwo besteht da ein fundamentales Verständnisproblem: die Listen von Koordinatenpaare in dem GeoJOSN Format definieren Umrisse von Objekten, die willst du also nicht?

Das ist eher ein Kommunikationsproblem: Er will nach meinem Verständnis die Umrisse als Koordinatenpaare, dazu noch “irgendwie” als Fläche möglichst mit Attributen gruppiert. Als Umriss will er nur Polygone, keine MP, damit zwingend ohne Exklaven/Enklaven (“Inseln”?), also “wambacher light”.

Siehe das Posting das ich gequotet habe:

[5.8663877,51.0511197],[5.8668165,51.0512156],[5.8669173,51.0512381],[5.866737,51.0515704]

ist nun mal -nicht-

A1 5.8663877,51.0511197
B1 5.8668165,51.0512156
C1 5.8669173,51.0512381

Drum habe ich den Verdacht, dass das immer noch eine erhebliche Verständnislücke ist.

Hypothalamus · 2019-01-14T12:27:31.000Z

Moin Walter,
ok das habe ich soweit verstanden. Gäbe es denn die Möglichkeit aus dem angezeigten Multipolygon mit einem Programm wie ArcGis oder so die Inseln usw einzelnd herauszulösen und diese als extra Layer zu speicher ?
Gruß
Dennis

Hypothalamus · 2019-01-14T12:30:30.000Z

Siehe das Posting das ich gequotet habe:

[5.8663877,51.0511197],[5.8668165,51.0512156],[5.8669173,51.0512381],[5.866737,51.0515704]

ist nun mal -nicht-

A1 5.8663877,51.0511197
B1 5.8668165,51.0512156
C1 5.8669173,51.0512381

Drum habe ich den Verdacht, dass das immer noch eine erhebliche Verständnislücke ist.

Deswegen mein usw… . also Ziel wäre es schon aus der Json alle Lat/Lon Koordinate der Grenzen diese untereinander stehen.
Habe das Beispiel nur nicht fertig geschrieben.

Stelle dir einfach vor alle Lat Lon Daten aus deiner Json Datei stehen jeweils in einer zeile untereinander anstatt aneinandergereiht in einer Zeile.

Ob es als text Datei oder Word oder sonst ist eigentlich egal.

Im Endeffekt brauch ich die die Ländergrenze von Deutschland mit den Koordinaten das es einen geschlossenen Kreis gibt, und dies möglichst genau. Deswegen beste Quelle die ich gefunden habe sind die von Walter.

Störend deshalb das Multipolygon aufgrund der Inseln usw. Wenn ich natürlich die Inseln rausziehen könnte und als einzelnes Polygon speichern und die anderen Löcher schliesse.

Gruß
Dennis

wambacher · 2019-01-14T13:17:09.000Z

Hypothalamus:

Gäbe es denn die Möglichkeit aus dem angezeigten Multipolygon mit einem Programm wie ArcGis oder so die Inseln usw einzelnd herauszulösen und diese als extra Layer zu speicher ?

Hi Dennis,

wie schon geschrieben solltest du “meine” Daten als Shapes nehmen. Diese kannst du - wenn du es kannst - in QGIS bearbeiten.

Nur wie genau das geht, musst du schon selber herausfinden. Ich bin mir aber sicher, dass das an sich geht, da QGIS einen Editor beinhaltet.

QGIS funzt auch unter Windows und kostet nix.

Gruss
walter

maxbe · 2019-01-14T14:20:38.000Z

Ich hab das, was ich in #14 beschreiben hab auch mal mit Deutschland aus Wambachers Webseite ausprobiert: runtergeladen als Shapefile, Enklaven gelöscht, Geometriewerkzeug angewendet und als json gespeichert.

Wenn man sich vorher um die einzige “Enklave” in Jungholz kümmert und die Grenze dort ummalt, funktionierts. Man bekommt 189 Polygone ohne Löcher: Ein grosses Stück Festland und 188 Inseln. Die meisten davon liegen im Meer, manche in Belgien, ein paar in der Schweiz (eine echte Exklave in Büsingen und ein paar Stellen, wo ein Stück Deutschland auf dem Schweizer Bodenseeufer liegt).

Mit eingemeindetem Jungholz (ich hab die Grenze dort gelöscht), sieht die Json für Deutschland so aus (gezipptes json 1.7MB)

Ich denke das wäre mit vertretbarem Aufwand auch mit den anderen Ländern hinzubekommen, sooo viele Enklaven gibts ja nicht. Und dass man sich um die Gedanken machen muss, liegt auf der Hand. Wenn man San Marino und den Vatikan behalten will, bekommt man eben kein ungelochtes Italien und wer Büsingen und Campione berücksichtigt hat Löcher in der Schweiz.

Grüße
Max

gormo · 2019-01-14T14:22:58.000Z

Hypothalamus:

Im Endeffekt brauch ich die die Ländergrenze von Deutschland mit den Koordinaten das es einen geschlossenen Kreis gibt, und dies möglichst genau. Deswegen beste Quelle die ich gefunden habe sind die von Walter.

Weil “die Grenze” nicht so genau ist: du willst alle Landflächen, die zu Deutschland gehören? Also eine Fläche für den “Kontinent-Teil” (wir ignorieren https://de.wikipedia.org/wiki/B%C3%BCsingen_am_Hochrhein)), sowie eine Fläche pro Insel (für eine noch zu bestimmende Definition von Insel)?

Seegrenzen interessieren Dich nicht? Grenze im Bodensee ist egal?

Du willst jeden Stützpunkt dieser Grenzlinien dann als ein Koordinatenpaar extrahieren, und die Koordinatenpaare geordnet in einer Datei haben?

Es sollte das Plugin MMQGIS für QGis geben, das das kann: https://gis.stackexchange.com/questions/109553/how-to-export-polygon-coordinates-to-csv-not-wkt

Wenn Du willst könntest Du noch sagen, warum Du das so brauchst - vielleicht gibts einfachere Wege.

Die_immer_lacht · 2019-01-14T15:48:27.000Z

Ich finde es wirklich gefährlich dass eine ganze Riege von erfahrenen OSMlern nicht begreift, dass es sich bei Hypothalamus um einen Hilfevampir handelt, dem OSM völlig egal und der sich lediglich von euch KOSTENLOS SEINE ARBEIT MACHEN LÄSST.

Viel Spaß dabei noch!

GerdP · 2019-01-14T16:15:41.000Z

Was ist denn daran schlecht? Es gibt Menschen, denen es Spaß macht, Probleme zu lösen und es gibt andere Menschen, die dafür bezahlt werden, Lösungen zu finden. Wenn auf diese Art eine Lösung herauskommt, dann ist der Arbeitgeben von Hypothalamus zufrieden, er selbst auch und die Problemlöser hatten ihren Spaß.

maxbe · 2019-01-14T17:52:59.000Z

Die immer lacht:

… um einen Hilfevampir handelt …

Vampire sind ja nicht mehr die Nosferatus von früher, inzwischen sind das richtig nette Leute geworden.

Erst gestern hab ich eine Kochbuchseite ausgesaugt. War lecker und ich hab mich nicht mal bedankt… Solange mich irgendwelche unbekannten Typen mit know how versorgen, das mir fehlt, habe ich nix dagegen, Daten oder Wissen oder Code kostenlos abzugeben.

Hypothalamus · 2019-01-14T21:13:15.000Z

Die immer lacht:

Ich finde es wirklich gefährlich dass eine ganze Riege von erfahrenen OSMlern nicht begreift, dass es sich bei Hypothalamus um einen Hilfevampir handelt, dem OSM völlig egal und der sich lediglich von euch KOSTENLOS SEINE ARBEIT MACHEN LÄSST.

Viel Spaß dabei noch!

Woher hast du diese Informationen ?
Ja es handelt sich um ein Projekt. Eine Arbeit wie im meinem Post beschrieben ist nicht gleich eine Arbeit meines Arbeitsgebers ! Es handelt sich um ein privates Projekt im Rahmen meines Hobbies.
Wie kommst du darauf das es sich um meinen Arbeitgeber handelt ?

LG
Dennis

Hypothalamus · 2019-01-14T21:14:22.000Z

Vielen Dank für diejenigen die geholfen haben. Ich unterlass dann mal lieber das illegale Nachfragen bevor es noch weiter Probleme gibt und wurschtel mich mal weiter durch.

Hypothalamus · 2019-01-16T12:35:49.000Z

So ich wollte Euch mit euren Hinweisen und meinen selbsterarbeiteden gewurschtel ein Ergebnis präsentieren.
Ich habe mit Hilfe von Qgis (die Version 2.18) und diesem Python Script https://gis.stackexchange.com/questions/166934/python-library-for-converting-geojson-multi-polygon-to-polygon es geschafft die Multipolygone aufzulösen.

Leider habe ich des nicht geschafft ein Shape File wie von Walter empfohlen zu bearbeiten. Leider fehlen hier die Attribute und Layer. Wenn es denn damit zusammenhängt.

Ich werden jetzt mit den aufgeschlüsselten einzelnen Polygonen ein neues Shapefile bzw. Geojson mit mehreren Layern bauen wenn es denn so einfach klappt. Das ich dann die Inseln, Enklaven usw. dementsprechen bennenen kann. Wenn das denn so ohne weiteres klappt. Aber das wisst ihr wahrscheinlich besser.

Gruß
Dennis

maxbe · 2019-01-16T13:10:07.000Z

Hypothalamus:

Leider habe ich des nicht geschafft ein Shape File wie von Walter empfohlen zu bearbeiten. Leider fehlen hier die Attribute und Layer. Wenn es denn damit zusammenhängt.

Du hast aber schon wirklich ein shp-File runtergeladen oder dein JSON als SHP gespeichert? Mein Qgis (2.8.1) weigert sich nämlich, JSONs zu bearbeiten. Wenn ich ihm Shapedateien gebe oder mein JSON als SHP speicher und das dann lade, habe ich auch “rechte Maustaste” → “Bearbeitungsstatus umschalten” auf dem Layer und kann das editieren…