Не могу правильно повернуть точку, используя GeoPandas - PullRequest
Новая
       44

Не могу правильно повернуть точку, используя GeoPandas

0 голосов
/ 06 мая 2019

Мой код Python: 

from geopandas import GeoSeries
from shapely.geometry import Point

WGS84 = {'init': 'epsg:4326'}
p1 = (41.8121843, 45.2178516,)
p2 = (41.8124878, 45.2177536,)

point2 = Point(p2)

s2 = GeoSeries(point2, crs=WGS84)
rotated = s2.rotate(90, p1)  # Try to rotate p2 around p1

print(rotated)  # 41.8122823 45.2181551

Мой код JS (Turf.js): 

var p1 = [41.8121843, 45.2178516];
var p2 = [41.8124878, 45.2177536];

var l = turf.lineString( [p1, p2] );
L.geoJson(l, {style: {color: '#0000ff', fillOpacity: 1, fillColor: '#0000ff', opacity: 1}}).addTo(map);

// Rotate line around p1
var trufRotated = turf.transformRotate(l, -90, {pivot: p1});
L.geoJson(trufRotated, {style: {color: '#ff0000'}}).addTo(map);

Иллюстрация. Синий - это оригинальная линия. Красная линия повернута с использованием Turf.js (ок). Зеленая линия повернута с помощью GeoPandas (не в порядке):

enter image description here

Зеленая линия, повернутая с помощью геопанд, имеет неправильное смещение. Возможно, это проблема проецирования, но я не могу понять, как это исправить.

У меня вопрос, как правильно вращаться с помощью GeoPandas?

1 Ответ

0 голосов
/ 06 мая 2019

решаемая.Благодаря @PaulH.Мне нужно сделать вращение в Mercator и перевести его снова в lon / lat: 

WGS84 = {'init': 'epsg:4326'}
MERC = {'init': 'epsg:3857'}

p1 = (41.8121843, 45.2178516,)
p2 = (41.8124878, 45.2177536,)

point1 = Point(p1)
point2 = Point(p2)

s1 = GeoSeries(point1, crs=WGS84)
s2 = GeoSeries(point2, crs=WGS84)

s11 = s1.to_crs(MERC)
s22 = s2.to_crs(MERC)

rotated = s22.rotate(90, s11[0])
rotated = rotated.to_crs(WGS84)
print(rotated)
...