Получение координат зданий в населенном пункте из источника OpenStreetMap

Получение координат зданий в населенном пункте из источника OpenStreetMap

В открытой базе данных https://www.openstreetmap.org/ содержатся географические координаты множества зданий. Используя эту информацию можно получить базу соответствия:
Адрес Координаты.
Для стран СНГ данные в виде файлов PBF размещены на сайте https://needgeo.com/
Файлы разбиты по областям

Скачиваем себе на ПК файл, который содержит искомый данные нужного города. Я для примера возьму файл по Самарской области для получения координат домов для города Самара.
Создав на жестком диске папку «Координаты», положим туда файл RU-SAM.pbf.
Для извлечения из этого файла координат адресов понадобятся утилиты:
osmconvert (https://wiki.openstreetmap.org/wiki/RU:Osmconvert)
osmfilter (https://wiki.openstreetmap.org/wiki/Osmfilter)
Скачиваем их и помещаем в ту же созданную ранее папку «Координаты». Для ОС Windows 64x файлы будут такие: osmconvert64.exe и osmfilter.exe
Для начала вырезания из файла RU-SAM.pbf данных необходимо указать область координат, которую требуется вырезать. Область координат указывается в виде набор координат полигона. Область координат можно сформировать на сайте http://geojson.io/ с помощью инструмента «Draw a polygon»

Чтобы понять, что выделять можно воспользоваться картами от Яндекс. Здесь при поиске города (населенного пункта на карте выделяются его границы.
Для Самары:

Повторяем что-то подобное на сайте http://geojson.io/. Главное здесь – не захватить данные других городов, где могут повторятся адреса. Улица Ленина, например, есть везде, в каждом населенном пункте.

Набор получившихся координат справа необходимо преобразовать к виду, который поймет программа osmconvert. Для этого создается обычный текстовый файл city.poly в который вписываются координаты полигона города. Пример содержимого файла для Самары:
city
1
              50.1416015625               53.373497900586855
              50.17181396484375           53.30297979468481
              50.15533447265625           53.238920640924974
              50.05096435546875           53.19616119954287
              50.020751953125             53.16324027106289
              50.22125244140625           53.1335898292448
              50.32012939453125           53.19122467094173
              50.34759521484375           53.29641404214567
              50.34210205078125           53.39479437775073
              50.2130126953125            53.46516070239169
              50.07843017578125           53.50438429010248
              50.009765625                53.46025520362287
              50.06744384765625           53.420990805007634
              50.1416015625               53.373497900586855
END
END
Основное здесь – это координаты долготы и широты, взятые с сайта http://geojson.io/ после выделения полигона города.
Файл с координатами полигона города city.poly нужно положить в ту же папку «Координаты».
Переходим в режим командной строки и выполняем последовательно три команды:
osmconvert64.exe -B=city.poly --all-to-nodes RU-SAM.pbf -o=1samara.o5m
osmfilter.exe 1samara.o5m --keep="building AND addr*" --drop-author --drop-version -o=2samara.o5m
osmconvert64.exe 2samara.o5m -o=3samara.csv --csv-headline --csv-separator=; --csv="@id addr:street addr:housenumber @lat @lon"

После этого в папке «Координаты» будут созданы файлы: 1samara.o5m, 2samara.o5m и 3samara.csv
Первые два файла можно удалить. Они промежуточные. Третий файл 3samara.csv содержит нужные нам данные: Его можно открыть в EXCEL. 
Для EXCEL 2016 переходим во вкладку «Данные»  «Получение внешних данных»  «Из текста». Выбираем наш итоговый файл (3samara.csv), Проходим мастер импорта текста. На втором этапе установить галочку «точка с запятой».

Получаем представление адресов и координат, которое можно использовать. В файле есть информация о ID адреса в базе OpenStreetMap, название улицы, номер дома и координаты.