Я пытаюсь построить шейп-файл как точки карты Folium, но выдергиваю волосы, пытаясь получить правильную проекцию. Данные в системе EPSG: 28992, новой системе координат RD. Это проекционная система координат, поэтому координаты указываются в метрах. Чтобы нанести это на карту листов, мне нужно преобразовать его в EPSG: 3857, который работает с градусами, если я прав. К сожалению, я получаю полную чушь после преобразования его с помощью функции to_crs().
Сначала загружаю данные:
import geopandas
WTG = geopandas.read_file('WTGs.shp')
print(WTG.head()) показывает это:
TEXTSTRING geometry
0 A27-09 POINT (151809.960 484552.520)
Что правильно. У него также есть правильный CRS, как показано в выводе print(WTG.crs):
epsg:28992
Теперь, чтобы преобразовать это в другой CRS, я использую следующий простой код:
WTG = WTG.to_crs(epsg=3857)
Но теперь проверка результата с помощью print(WTG.head()) дает полную чушь для столбца геометрии.
TEXTSTRING geometry
0 A27-09 POINT (594489.076 6863453.514)
Геометрические точки должны быть в градусах. Другими словами, координаты должны быть около 52 и 5. Вывод print(WTG.crs) верен:
epsg:3857
Кто-нибудь из вас знает, как это исправить? Или я что-то упускаю? Заранее спасибо!