Преобразовать шейп-файл в географические координаты - PullRequest
0 голосов
/ 04 февраля 2019

Очень новые данные для geospatial - в настоящее время я работаю над тем, чтобы получить набор почтовых индексов, помеченных в пределах их правильного географического региона, как определено NUTS1.Я подумал, что лучший способ сделать это будет:

  1. Получить координаты областей NUTS1;
  2. Преобразовать их в многоугольник;
  3. Получить координатыпочтового индекса и посмотрите, какой многоугольник находится внутри

Однако - я столкнулся с барьером с шагом (1) - данные от ONS относительно областей NUTS не в долготе /формат широты из того, что я собираю (см. источник и ошибку ниже).Из аналогичных вопросов я понимаю, что ogr2ogr - это способ конвертировать это - поэтому я скачал OSGeo4W, в котором есть ogr2ogr, однако я не совсем то, что я должен запускать в своей OSGeo4W Power Shell, чтобы сделатьсоответствующие преобразования.

Я также ознакомился с приведенной ниже шпаргалкой, но все же не повезло.

Я понимаю, что у геопанд есть очень простой ответ на мою проблему - однако после 2 дней попытки установить этот пакетМне все еще не повезло, поэтому я временно отступаю от этой авеню (на отдельной ноте, если кому-то повезло с установкой геопанд с Python 3.7, укажите мне правильное направление. Я попыталсяКонда установила и загрузила все whl-файлы по отдельности, ни один из них, похоже, не работал).

Я подумал, может ли кто-нибудь написать для новичка точную команду, которую мне нужно выполнить, чтобы решить мою проблему.

Заранее спасибо!

Источник данных: http://geoportal.statistics.gov.uk/datasets/nuts-level-1-january-2018-full-extent-boundaries-in-the-united-kingdom/data

Ошибка: ValueError: шейп-файл должен иметь широту / долготу

Похоже, у этого есть вершины в координатах проекции карты.Вы можете преобразовать шейп-файл в географические координаты, используя утилиту shpproj из инструментов shapelib (http://shapelib.maptools.org/shapelib-tools.html)

Шпаргалка: http://www.bostongis.com/PrinterFriendly.aspx?content_name=ogr_cheatsheet

1 Ответ

0 голосов
/ 07 февраля 2019

Ваш первый шаг - преобразование исходных данных во вездесущую систему координат WGS84 , используя опцию -t_srs для программы ogr2ogr cli.Подробнее обо всех текущих параметрах ogr2ogr здесь .

ogr2ogr boundaries.shp NUTS_Level_1_January_2018_Full_Extent_Boundaries_in_the_United_Kingdom.shp -t_srs epsg:4326

Набор данных уже находится в форме 'ST_Polygon', поэтому шаг 2 является спорным.

Я оставлю это в качестве упражнения для вас, чтобы выяснить шаг 3самостоятельно;подсказка: вы можете делать пространственные объединения с помощью ogr2ogr, используя -sql -dialect sqlite.

...