Как вырезать сабсет города (любого отношения) из OSM данных

Однажды мне потребовалось получить из сырых ОСМ данных чистый сабсет города (потому что так удобно, компактно и просто красиво). К моему удивлению я не нашел готового рецепта, из-за чего для выполнения этой задачи потребовалось немного попотеть.

Ввиду вышеизложенного считаю своим долгом поделиться накопленным опытом.



Допустим я хочу получить сабсет Кишинёва (почему? … там тепло, там мой дом, там моя мама)

Planet.osm нас, конечно, не интересует, зато есть geofabrik, но там планета нарезана на достаточно крупные шматки — http://download.geofabrik.de/europe/moldova.html

Значит, граница города есть в сырых данных.

Из отношения несложно получить полигон, делаем это либо руками, либо при помощи любой подходящей утилиты, например, этой.

Иии… теперь в игру вступает osmosis.

У нас есть осм всея Молдовы и ограничивающий полигон Кишинёва, этого вполне достаточно для получения сабсета города из сырого XML (OSM) файла:

osmosis --read-xml file="moldova.osm" --bounding-polygon file="kishinev.poly" --write-xml file="kishinev.osm"

Аналогично можно выковырять данные из PBF (Protocolbuffer Binary Format) файла:

osmosis --read-pbf file="moldova.osm.pbf --bounding-polygon file="kishinev.poly" --write-xml file="kishinev.osm"

Немногим больше про osmosis можно почитать здесь.

Ну а готовый сабсет нужно срочно открывать в JOSM или импортировать в личный проект!