Как преобразовать EPSG: 4326 координат в EPSG: 3857 в r, используя rgdal? - PullRequest
Новая
       56

Как преобразовать EPSG: 4326 координат в EPSG: 3857 в r, используя rgdal?

0 голосов
/ 09 января 2019

Мне удалось преобразовать мой набор координат, спроецированный в EPSG: 4326, в EPSG: 3857, однако результат не совпадает с теми, которые я получаю при

https://epsg.io/transform#s_srs=4326&t_srs=3857&x=14.5172200&y=46.0658300 или непосредственно в ArcGIS, хотя особенности систем EPSG, представленные в r (см. Скрипт ниже), соответствуют описанию в программном обеспечении ГИС.

Я пробовал это со следующим кодом: 

library("rgdal", lib.loc="~/Library/R/3.4/library")


orig_coords <- data.frame(lat=c(46.065830, 46.042211, 46.094612), lon=c(14.517220, 14.487756, 14.597046))
coordinates(orig_coords) <- c('lat', 'lon')

#Determine the projection of the lat-long coordinates, by default it is EPSG:4326
proj4string(orig_coords) <- CRS("+init=epsg:4326")
print(summary(orig_coords))

#Convert the coordinates to the used metric system (EPSG:3857)
Metric_coords<-spTransform(orig_coords,CRS("+init=epsg:3857"))
print(summary(Metric_coords))

Правильные координаты будут 

5790904.807 1616049.538;
5787116.145 1612769.621;
5795523.844 1624935.728;

Вместо этого я получаю: 

lat - lon;
[1,] 5128025 1633624;
[2,] 5125395 1630236;
[3,] 5131229 1642804;

Я не могу найти, где проблема ...

1 Ответ

0 голосов
/ 25 июня 2019

Я считаю, что CRS("+init=epsg:3857") должен быть прописным CRS("+init=EPSG:3857")

...