search for: sp2transformed

Displaying 1 result from an estimated 1 matches for "sp2transformed".

Did you mean: sp1transformed
2012 Oct 18
1
spTransform longlat to utm
....x1 coords.x2 [1,] 998395.0133 6319888.068 The result is an obvious nonsense. Well, after some deliberation I swapped the original coordinate columns: x <- 49.454735 y <- 8.968735 sp2 <- SpatialPoints(matrix(c(x,y), ncol=2), proj4string = CRS("+proj=longlat +datum=WGS84")) sp2Transformed <- spTransform(sp2, CRS("+proj=merc +zone=32u +datum=WGS84")) coordinates(sp2Transformed) coords.x1 coords.x2 [1,] 5505275.918 995840.692 Now the northing (which comes first after swapping) looks better, at least it falls within Germany, which is however not to say about the...