OSM-Daten Gebiete ausschneiden

als Erweiterung zu meiner Anleitung OSM-Daten filtern geht es hier um das ausschneiden eigener Gebiete aus den OSM-Daten.

Vorbereitungen

Es werden wieder die Programme osmconvert und Wget benötigt.
Wer die Anleitung zum filtern der OSM-Daten gemacht hat, benötigt keine weiteren Dateien.

Rechteckige Ausschnitte

Die einfachste Möglichkeit ein Gebiet auszuschneiden ist die Angabe zweier Punkte.
Dazu gibt man die Koordinaten des südwestlichsten Punktes und des nordöstlichsten Punktes an.
Die zwei Punkte markieren dann die linke untere Ecke sowie die rechte obere Ecke eines Rechteckes das ausgeschnitten werden soll.

Um bei meinem Beispiel Andorra zu bleiben, muss die Befehlszeile um die Stadt Andorra (Andorra la Vella) aus der andorra-latest.osm.o5m auszuschneiden folgendermaßen lauten:

C:\Beispiel\osmconvert64 C:\Beispiel\andorra-latest.osm.o5m -b=1.48,42.48,1.56,42.53 --out-o5m > C:\Beispiel\Andorra_la_Vella_rechteck.o5m

Freie Asschnitte

Um eine 'freie' Fläche, also ein Vieleck oder Polygon auszuschneiden, benötigt man eine Polygon-Datei mit den entsprechenden Koordinaten.
Der einfachste Weg an eine Polygon-Datei zu kommen, ist der download einer Datei von der Geofabrik.de, die bieten für (fast) alle ihrer Dateien .poly zum download an.
Möchte man eigene Gebiete ausschneiden muss man sich natürlich selber eine Polygon-Datei erstellen, wie ich mir die Koordinaten besorgt habe, beschreibe ich unter 'Spielereien mit OSM-Daten'.
Dazu legen wir im Texteditor wieder eine neue Datei an und speichern diese, ich gebe ihr den Namen "Andorra_la_Vella.poly".

• In die erste Zeile kommt der Name der Datei.
• In die zweite Zeile kommt die Sektionsbezeichnung. Sie markiert den Beginn eines Gebietes / Polygon.
• Die jetzt folgenden Zeilen enthalten die Koordinaten des Gebietes / Polygons beginnend mit Leerzeichen, dann zuerst der Längengrad und dann der Breitengrad getrennt durch ein Leerzeichen. Normalerweise werden die Koordinaten im wissenschaftlichen Format angegeben, da osmconvert aber auch normale Gradzahlen akzeptiert werde ich der einfachheithalber Gradzahlen direkt aus Quovadis nehmen.
• Danach kann zum schließen des Polygons die erste Koordinatenzeile wiederholt werden. Macht man das nicht, zieht das Programm eine Linie vom letzten Punkt zum ersten und schließt damit das Polygon.
• In der nächsten Zeile muss mit dem Wort END das Gebiet abgeschlossen werden.
  • Man könnte jetzt noch ein weiteres Gebiet angeben das zusätzlich ausgeschnitten werden soll oder ein Gebiet angeben das aus dem ersten Gebiet ausgespart werden soll.
    Dazu müssen die Zeilen vom Schritt der Sektionsbezeichnung bis zum Abschluss mit END wiederholt werden.
    Im Falle das die Sektion eine andere Sektion aussparen soll, muss die Sektionsbezeichnung mit einem Ausrufezeichen beginnen.
• Nun muss die Datei noch mit einem weiteren END abgeschlossen werden.

Meine Polygon-Datei sieht nun so aus (die drei Punkte sind nur Platzhalter für weitere Koordinaten):

Andorra la Vella poly
1
   1.5224247 42.5024807
   1.5233568 42.5034325
   ...
   1.5218413 42.5023706
   1.5224247 42.5024807
END
END

Um ein Gebiet mittels dieser Polygon-Datei auszuschneiden lautet die Befehlszeile:

C:\Beispiel\osmconvert64 C:\Beispiel\andorra-latest.osm.o5m -B=C:\Beispiel\Andorra_la_Vella.poly --out-o5m > C:\Beispiel\Andorra_la_Vella_poly.o5m

Hier noch die Übersicht über die ausgeschnittenen Gebiete:

Integration in die Batchdatei zum Filtern

Natürlich lässt sich die Aufgabe in die Batchdatei integrieren und somit in einem Rutsch erledigen.
Dazu die Befehlszeile nach dem konvertieren der .pbf zu .o5m aber vor dem Filtern eintragen.

Arbeitsdateien

Zum Abschluss wieder die Ausschnitt.cmd, die Andorra_la_Vella.poly und die Parameter_POI_Camping.txt zum herunterladen:
Beispiel2.zip