PLZ von Gemeinden - Datenbankabfrage?

Hallo zusammen,
ich bin auf der Suche nach einem Weg, die PLZ von Gemeinden aus der Datenbank auszulesen. Ich möchte (falls möglich automatisiert) eine Liste Gemeinden anhand ihres Schlüssels in eine Abfrage stecken und eine Liste aller PLZ erhalten, deren Gebiet die jeweilige Gemeinde überschneidet. Leider kenne ich mich bisher kaum mit Overpass API aus.
Kann mir jemand sagen, ob solch eine Abfrage machbar ist?

Ja.

Sie ist. Mit der CSV-Ausgabe kann man auch nur die Postleizahl-Liste erhalten.

/edit: Link für CSV

Geht bei mir mit den meisten Gemeinden nicht - Output Menge ist leer.

Wie ist denn bei Overpass “(area.searchArea)” nach DE-9IM definiert?

F*… Kann man irgendwie einen ganz ganz kleinen Buffer um die Grenze legen?

Prinzipiell geht das ja schonmal in die richtige Richtung.

Blöd, dass der area Filter die relations nicht erfasst, wenn sie nur Elemente auf dem Rand der area haben. Da muss es doch irgendeine Lösung für geben?

Kann ich die Area auch über den Gemeindeschlüssel definieren? geocodeArea akzeptiert ja nur den Namen der Gemeinde oder? Und der ist ja leider nicht immer eindeutig…

Irgendwie mit around, aber das verstehe ich nicht

http://wiki.openstreetmap.org/wiki/Overpass_API/Overpass_QL#Relative_to_other_elements_.28around.29

Das ist doch alles Gefrickel. Gibt es in Overpass also kein vernünftiges Intersection-Konzept? Dann doch besser PostGIS. Ist aber natürlich auch nicht so einfach…

http://overpass-turbo.eu/s/6gI

[out:json][timeout:50];
{{geocodeArea:Landkreis Goslar}}->.searchArea;
(
  relation["boundary"="administrative"]["de:amtlicher_gemeindeschluessel"="03153017"](area.searchArea)->.b;
  relation["boundary"="postal_code"](around.b:1)
);
out body;
>;
out skel qt;

Man muss anpassen

  • die Nominatim-Area, in der man sucht (im Zweifelsfall Deutschland nehmen)
  • den AGS, nachdem man sucht

…Aber damit kriegt man auch die PLZ, die direkt an die Gemeinde mit dem AGS angrenzen.

Als Basis zum Weiterarbeiten ist das aber vielleicht schon nützlich

Eben. Einiges geht schnell, einiges nicht. So wie immer :wink:

Hier mal ein ganz anderer Ansatz mit is_in: http://overpass-turbo.eu/s/6iR


[out:json][timeout:300];

// Area für Merchweiler in .a1 merken
area["boundary"="administrative"]["de:amtlicher_gemeindeschluessel"="10043113"]->.a1;
// Area in .a1 umwandeln in eine Relation und dafür die Wege/ Knoten ermitteln
// Ergebnis in .n1 merken. In .n1 sind nun alle Wege/Knoten, die den
// Rand der boundary=administrative Relation beschreiben
rel(pivot.a1); > -> .n1;
// Für die Area .a1 (Merchweiler) alle Knoten ermitteln, die einen Tag haben
// Annahme: es gibt mindestens einen Knoten in der Area
// Ergebnis in .n2 merken. Knoten, die den Rand der Area beschreiben
// sind dort noch enthalten
node[~"."~"."](area.a1) -> .n2;
//Randknoten rauswerfen
(.n2; - .n1;);
//Alle Areas ermitteln, in denen innere Knoten in Merchweiler vorkommen
is_in;
//Filtern auf Areas mit boundary=postal_code
area._["boundary"="postal_code"];
//Area wieder in eine Relation umwandeln
rel(pivot);
//und ausgeben
out geom;

Ich bin gerade auf die openGeoDB nodes gestoßen, die sowohl die Gemeindeschlüssel als auch die der Gemeinde zugeordneten PLZ enthalten.

k="openGeoDB:community_identification_number"
k="openGeoDB:postal_codes"

Wenn man davon ausgeht, dass diese Daten in der OSM Datenbank noch aktuell sind, könnte man doch über diese nodes leicht an die PLZ kommen. Ich habe mal stichprobenartig 15 Gemeinden quer durch Deutschland überprüft und die Daten sind sowohl vorhanden als auch korrekt.

[out:csv("openGeoDB:postal_codes";false)][timeout:25];
(
  node["openGeoDB:community_identification_number"="16069001"];
);
out;

Wie würde ich denn nun vorgehen, wenn ich sagen wir die PLZ aller ~11000 Gemeinden abfragen möchte? Gibt es eine Möglichkeit, mit der Overpass API eine Liste aller Gemeindeschlüssel zu importieren und dann beispielsweise über eine for-Schleife alle Gemeinden abzuarbeiten?

Vorsicht: die openGeoDB wurden irgendwann mal importiert, werden aber kaum gepflegt.

Nahmd,

<mode=“Advocatus diaboli”>
Die — rechnet man seine Arbeitszeit mit — kostengünstigste Lösung ist, bei Postdirekt die CD “DATAFACTORY BASIC” zu bestellen, die Daten in eine MySQL/PostgreSQL/YouNameItSQL-DB einzuspielen und darauf Selects loszulassen.

Gruß Wolf

Also für ganz DE wird dir Roland sicher früher oder später auf die Finger hauen (->lese: eigene Instanz aufbauen!!), aber für’s Saarland kann man das schon noch laufen lassen:

http://overpass-turbo.eu/s/6iS

Prinzipiell würde ich empfehlen, einen Extrakt von Geofabrik zu nehmen und damit eine eigene DB zu füttern.

Oder anders gesagt: Wenn du die Daten aus diesen Tags nutzen willst nutze lieber direkt die openGeoDB.

Eine Lösung für die Grenzsuchanfrage (aka “zu welchem Landkreis gehört Kleinstadt”) würde mich jedenfalls dennoch interessieren.

OpenGeoDB ist da zu veraltet, soweit ich informiert bin.

Das geht auch mit is_in - einfach mal selbst ausprobieren… :sunglasses:

Ich beginne zu verstehen. Vielen Dank!

is_in wird aber doch nicht funktionieren, wenn sich z.B. 2 Gemeinden eine PLZ teilen, oder?

nö, und ausserdem ist is_in am Aussterben und auf keinen Fall komplett.

Gruss
walter

Wie bitte?

Ich rede hiervon: http://wiki.openstreetmap.org/wiki/Overpass_API/Overpass_QL#Query_for_areas_.28is_in.29

Beispiel: http://overpass-turbo.eu/s/6iu