проецировании SpatialLines объекта от широты/долготы к UTM

Question

Я создал большого круга траекторию широты/долготы с использованием geosphere пакета:

flightTraj = greatCircle(c(originAptLon, originAptLat), c(destinAptLon, destinAptLat), n = nPts, sp = TRUE) 

его свойства:

class  : SpatialLines 
features : 1 
extent  : -180, 180, -52.74719, 52.74719 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +ellps=WGS84 

Я хочу проект по картографии UTM 14 с использованием эллипсоида WGS84. Я пробовал:

projectedTrajectories <- CRS("+proj=utm +zone=14 +datum=WGS84 +units=km +no_defs") %>% 
    spTransform(flightTraj, .) 

Но консоль отображается:

non finite transformation detected: 
    V1 lat 
    Error in .spTransform_Line(input[[i]], to_args = to_args, from_args = from_args, : 
    failure in Lines 1 Line 1 points 
    In addition: Warning message: 
    In .spTransform_Line(input[[i]], to_args = to_args, from_args =  from_args, : 
    6 projected point(s) not finite 

Любая помощь, пожалуйста.

Answer 1

"Все хорошо, что хорошо кончается ..."

distOrigDest = distGeo(c(originAptLon, originAptLat), c(destinAptLon, destinAptLat), a = earthRadius, f = earthFlattening)/km2m 
nPts = 20 * floor(distOrigDest/cellCentroidDist) 
flightTraj = gcIntermediate(c(originAptLon, originAptLat), c(destinAptLon, destinAptLat), n = nPts, sp = TRUE) 
flightTraj = spTransform(flightTraj, "+proj=utm +zone=14 +datum=WGS84 +units=km") 
plot(flightTraj, col= "blue", add = TRUE, lwd = 2)