У меня есть точки в географической системе координат, и я хотел преобразовать их в швейцарскую сетку (CH1903 +).
Образец данных:
id lon lat
2 7.173500 45.86880
3 7.172540 45.86887
4 7.171636 45.86924
5 7.180180 45.87158
6 7.178070 45.87014
7 7.177229 45.86923
8 7.175240 45.86808
9 7.181409 45.87177
10 7.179299 45.87020
Уважаемые результаты:
id E N
2 2579408.2431 1079721.1499
3 2579333.7158 1079729.1852
4 2579263.6502 1079770.1125
5 2579928.0358 1080028.4605
6 2579763.6471 1079868.9218
7 2579698.0756 1079767.9762
8 2579543.1019 1079640.6512
9 2580023.6226 1080049.2672
10 2579859.1889 1079875.2740
convert
coordinate-system
r
Topdombili
источник
источник
proj4string(d)
соответственно. Мое внимание было в первую очередь обращено на ложные параметры восточного и северного направления,x0
иy0
, поскольку некоторые популярные ссылки в Интернете (например, в первом комментарии в коде) опустили свои наиболее значимые цифры, сместив все точки на несколько тысяч километров. :-).Ответы:
Используйте пакет RGDAL . Есть проблема, какую CRS использовать; RGDAL не распознает код EPSG. Вы должны предоставить параметры явно, как показано здесь. (Очевидно, это приблизительные значения, но они должны быть довольно хорошими. Кажется, они находятся в пределах 0,1 м от предполагаемых значений.)
источник