Перепроецирование GeoPandas с 4326 на 3857 не работает должным образом - PullRequest
0 голосов
/ 08 октября 2019

У меня небольшие проблемы с проекцией DataFrame в GeoPandas из EPSG: 4326 в EPSG: 3857 в блокноте. Мой оригинальный набор данных выглядит следующим образом:

0       POLYGON ((-97.44573128938707 25.889635, -97.35...
1       POLYGON ((-97.61263173798922 25.995165, -97.57...
2       POLYGON ((-97.583158 25.96192090708935, -97.58...
3       POLYGON ((-97.461286 25.89270423248976, -97.46...
4       POLYGON ((-97.33941 25.926101719766, -97.33941...
                              ...                        
4971    POLYGON ((-102.906172 36.475082, -102.906172 3...
4972    POLYGON ((-102.772786 36.481234, -102.772786 3...
4973    POLYGON ((-102.639385 36.48728000000001, -102....
4974    POLYGON ((-102.50597 36.49322100000001, -102.5...
4975    POLYGON ((-102.37254 36.499056, -102.37254 36....
Name: geometry, Length: 4976, dtype: object

визуализируется:

Map version of data

Когда я пытаюсь перепроектировать это в EPSG: 3857:

reprojected = texas.to_crs(epsg=3857)

В качестве вывода я получаю несколько огромных чисел (которые, очевидно, не являются правильными координатами):

0       POLYGON ((-10847609.18711273 2985418.174237921...
1       POLYGON ((-10866188.46006429 2998482.12034216,...
2       POLYGON ((-10862907.45855956 2994365.44865411,...
3       POLYGON ((-10849340.7295776 2985797.960665139,...
4       POLYGON ((-10835773.55531768 2989931.200110516...
                              ...                        
4971    POLYGON ((-11455462.66652503 4366190.37784271,...
4972    POLYGON ((-11440614.20492607 4367042.077273317...
4973    POLYGON ((-11425764.07353476 4367879.167666097...
4974    POLYGON ((-11410912.38367058 4368701.784074328...
4975    POLYGON ((-11396059.02401403 4369509.784691155...
Name: geometry, Length: 4976, dtype: object

Любые идеи о том, почему такое простое перепроецирование не будетработать правильно? Я попытался использовать различные формы функции to_crs, включая texas.to_crs({'init': 'epsg:3857'}), и трижды проверил, что texas.crs установлен на {'init': 'epsg:4326'}.

Для справки я использую GeoPandas 0.5.1 и PyProj2.3.1. Редактировать: Я также пытался перейти на GeoPandas 0.6 только ради забавы, но не для игры в кости: (.

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

1 Ответ

3 голосов
/ 08 октября 2019

Я думаю, что эти цифры верны. Если вы посмотрите на этот сайт (https://epsg.io/3857), то увидите, что ваши значения находятся в прогнозируемых пределах.

...