Пакетное геокодирование в QGIS
Геокодирование в QGIS с помощью сервисов Google, Yandex, OSM(Nominatim).
Введение
Одним из способов получения геоданных является операция геокодирования. Она позволяет получить координаты объектов, у которых заранее известны адреса. В связи с большим объемом адресных БД, эту операцию обычно производят сервисы геокодирования. Наиболее известные сервисы на данный момент - Google Geocoding, Яндекс.Карты Геокодирование и OSM Nominatim.
QGIS не имеет встроенной возможности пользоваться сервисами геокодирования. Но гибкая система расширений позволяет добавлять практически любые функции. Для выполнения пакетного геокодирования можно использовать расширение RuGeocoder.
Подготовка к работе
Для выполнения геокодирования понадобятся:
- QGIS версии 1.6 и выше
- Расширение RuGeocoder (Присутствует в официальном репозитории)
- Соединение с Интернетом (желательно быстрое)
При подготовке данных так же может понадобиться расширенный калькулятор атрибутивных полей. В примерах будет использоваться расширение FieldPyculator.
Процесс установки QGIS подробно описан в статье Установка QGIS/GRASS с помощью OSGeo4W.
Более подробно об установке расширений можно прочитать тут.
Подготовка данных
Для работы геокодер требует точечный слой, с адресами, находящимися в одном или нескольких атрибутивных полях. Компонентная форма записи адреса (название города, название улицы, номер дома и т.д. находятся в отдельных полях) более предпочтительна, и обычно позволяет добиться значительно лучшего результата при использовании любого геокодера.
Преобразование таблицы в SHP файл
Достаточно часто, входные данные представлены в виде обычных таблиц. Наиболее удобный формат для обработки - CSV. С данным форматом могут работать практически все современные табличные редакторы: MS Excel, OpenOffice Calc и другие.
Для геокодирования таблицы её необходимо конвертировать в любой формат геоданных, которые можно редактировать в QGIS - SHP, sqlite, mapinfo и тд. Это преобразование можно сделать с помощью ogr2ogr, подробное описание можно прочитать в статье Конвертация данных из CSV в SHP и обратно с OGR. Но данный метод требует создания vrt файла и использование командной строки.
Для удобства RuGeocoder имеет встроенную функцию, с подобным преобразованием. На панели "Модули" необходимо нажать кнопку "Конвертировать CSV в SHP" . Тоже действие можно выполнить через меню "Модули" -> "RuGeocoder" -> "Конвертировать CSV в SHP". В открывшейся форме нужно выбрать входной CSV файл и указать путь, по которому будет сохранен результирующий SHP файл.
В результате конвертации, будет получен точечный shp файл с той же самой атрибутикой, и нулевыми координатами.
Если включена соответствующая опция, то к атрибутике будут добавлены пустые поля типа string: settlement, street, building_n и geocoded. Данные поля можно будет использовать в целях подготовки адреса для геокодера.
Кодировка данных выходного слоя будет такой же, как у входного. Если CSV файл был сохранен в UTF-8, то SHP будет иметь ту же кодировку. Возможно, после добавление слоя на карту, у него будет выставлена неверная кодировка (зависит от настроек по умолчанию). Изменить кодировку можно в свойствах слоя, на вкладке Общие, Параметры источника.
Рассмотрим для примера CSV файл (Файл:RuGeocoder Ufa-leninsky utf8.zip), содержащий адреса УИК Ленинского района г. Уфы. После конвертирования получаем следующий SHP файл: Файл:Ufa-leninsky utf8 shp.zip.
Подготовка адреса
Для повышения качества геокодирования необходимо разбить адрес на составные части: населенный пункт, улица, строение. Каждая из частей должна находиться в отдельном поле (поле с названием населенного пункта можно не помещать в отдельное поле, а указать вручную, если все адреса в слое принадлежат этому населенному пункту). Если исходные данные изначально были в такой форме, то этот пункт обработки можно пропустить.
Как примером воспользуемся ранее подготовленным SHP файлом. Добавим его в проект и с помощью FieldPyculator подготовим необходимые поля. Как видно по полю addr_vote, адреса указаны единообразно, но содержат лишнюю информацию - номера школ, название организаций и тд. Выделим название улиц, использовав следующее выражение для Pyculator:
value = <addr_vote>.split(',')[1].replace(u'ул.','').strip()
Для получения номеров строений используем следующий код:
value = <addr_vote>.split(',')[2].strip().split(' ')[0]
В результате, получаем SHP файл (Файл:Ufa-leninsky utf8 shp attr.zip) с заполненными полями street и building_n. Исключением является одна запись (УИК с номером 187), в поле street которой попало название поселка. Данный недостаток легко исправить вручную.
Больше примеров обработки адресов можно найти на странице в вики: Проект_УИК_Республики_Башкортостан.
Геокодирование
После предварительной подготовки адресов можно приступать к самому геокодированию. Геокодировать можно любой точечный слой, редактирование которого возможно в QGIS (PostGIS, sqlite, MSSQL и т.д.). Слой для обработки должен быть добавлен в проект и открыт на редактирование.
Форма обработки вызывается нажатием кнопки "Пакетное геокодирование" на панели "Модули". Тоже действие можно выполнить через меню "Модули".
В открывшейся форме необходимо выбрать слой для геокодирования.
Если в слое присутствуют адреса из нескольких населенных пунктов, то необходимо выбрать поле, содержащее название поселения. Если все адреса относятся к одному населенному пункту, то его название можно указать вручную. Аналогичная ситуация с районам, с той лишь разницей, что это поле необязательное.
Указание адреса несколько отличается. Если адрес записан в одном поле, то достаточно выбрать это поле.
Если же адрес компонентный, то необходимо включить переключатель "Номер здания", и соответственно указать поле содержащее название улицы, и поле содержащее номер строения.
Выбрав геокодер, можно запускать процесс обработки. Время выполнения геокодирования зависит от количества записей в слое, скорости соединения и загруженности выбранного сервиса.
После геокодирования
Результатом геокодирования является изменение геометрии. Координаты точки в таблицу не дописываются. Если вам необходимо табличное представление координат, выберите в меню Вектор->Обработка геометрии->Экспортировать\добавить поле геометрии