Что является альтернативой функции npts базовой карты в картопии?

Question

Я хочу эмулировать функцию ntps, которая находится в классе basemap.Geod.Его функция описывается следующим образом:

Given a single initial point and terminus point (specified by
        python floats lon1,lat1 and lon2,lat2), returns a list of
        longitude/latitude pairs describing npts equally spaced
        intermediate points along the geodesic between the initial and
        terminus points.

Я хочу сделать то же самое в картопии, очевидно, один из способов - это вычислить расстояние, необходимое между каждой парой последовательных координат, и использовать это расстояние для вычисления координат с начала.до конечной точки.Есть ли другой способ сделать это?

Answer 1

Исходя из источника cartopy, https://github.com/SciTools/cartopy/tree/master/lib/cartopy/geodesic, очевидно, что из cartopy.

не существует эквивалента функции ntps().вычислить расположение точек вдоль геодезической, если вы знаете направление от одной из конечных точек.Вот шаги, которые вы можете выполнить:

1. Используйте Geodesic.inverse () с вводом (long, lat) точек 1 и 2 и получите в качестве результата (forward_bearing, backward_bearing, geodesic_dist).
2. Предположим, что вы хотите получить (long, lat) точки (скажем, 1/4 от всего расстояния) вдоль этой геодезической (в 1), используйте Geodesic.direct () с long1, lat1, полученные forward_bearing и geodesic_dist/4.

Лучшая альтернатива - использовать pyproj и забыть все вышеперечисленное.