Проєктування точки
У попередньому розділі ми працювали з місцезнаходженням Ейфелевої вежі. Знову надаємо координати в проєктованій системі координат, щоб ви, наприклад, могли обчислювати відстані. Повернімося до цієї знакової пам'ятки та задамо її положення в географічних координатах: 48°51′29,6″ пн. ш., 2°17′40,2″ сх. д. Або в десятковому вигляді: широта 48,8584 і довгота 2,2945.
Об'єкти геометрії Shapely не мають у собі відомостей про СКВ, тому їх не можна напряму перевести в іншу СКВ. Тож скористаємося GeoPandas, щоб трансформувати точку Ейфелевої вежі в альтернативну СКВ. Ми помістимо одну точку в GeoSeries, застосуємо метод to_crs(), а потім знову дістанемо точку.
GeoPandas уже імпортовано.
Ця вправа є частиною курсу
Робота з геопросторовими даними в Python
Інструкції до вправи
- Створіть об'єкт точки shapely з координатами Ейфелевої вежі та присвойте його змінній
eiffel_tower. - Створіть GeoSeries (назвіть її
s_eiffel_tower) з Ейфелевою вежею як єдиним елементом і вкажіть СКВ EPSG:4326. - Перетворіть
s_eiffel_towerдо EPSG:2154 і збережіть результат уs_eiffel_tower_projected.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Construct a Point object for the Eiffel Tower
from shapely.geometry import Point
eiffel_tower = ____
# Put the point in a GeoSeries with the correct CRS
s_eiffel_tower = geopandas.GeoSeries([____], crs={'init': '____'})
# Convert to other CRS
s_eiffel_tower_projected = s_eiffel_tower.____
# Print the projected point
print(s_eiffel_tower_projected)