R перепроецирует карту на широту/долготу

Я использовал пакет mapproj для преобразования набора широты/долготы в картографическую проекцию (в этом примере lagrange). Затем я использовал его, чтобы определить несколько других точек в пространстве. Теперь я хочу преобразовать эти новые точки обратно в широту/долготу, но mapproj не предлагает такой функции, и я не нашел ее в другом месте.

cord<-data.frame(x, y)

lagrange<-mapproject(cord$x, cord$y, projection="lagrange")
#Some  
#clustering
#here

#Here are cluster centers, but in the form of lagrange projection. I want a lat/long cord of them
cluster.centers

Как это сделать в Р?


r map map-projections
person Paweł W.    schedule 08.10.2014    source источник
comment
Проверьте пакет rgdal и функцию project (в частности, аргумент inv).   -  person nicola    schedule 08.10.2014

Ответы (1)


arrow_upward
2
arrow_downward

Я прочитал документацию по пакету и не верю, что автор намеревался использовать свой набор проекций помимо построения графиков.

«Стандартный» способ работы с геопространственными данными в R — через набор объектов sp. Преобразовать объекты sp можно с помощью spTransform. И объекты, и методы можно найти в пакете sp. Должно быть несколько доступных руководств, по общему признанию, это может быть крутой кривой обучения, если вы не знакомы с ГИС и R.

Пример преобразования данных «длинной широты» в другие проекции подробно описан в этом посте:

Ошибка при повторном проецировании пространственных точек с использованием spTransform в rgdal R

person Jonathan Lisic    schedule 08.10.2014