Форумы GIS-Lab.info

Добрый день,

подскажите, пожалуйста, как средствами только R перепроектировать два векторых слоя (оба WGS'84, dd) так, чтобы центральный меридиан был сдвинут на отметку -200 Long, потому что хочется иметь картинку (см вложение) - центрированную на Тихом океане, а не на Гринвиче, чтобы ареал вида не разбивался на два куска (оранжевый полигон)

Спасибо большое

Посмотрите в документации elide-methods: пакет maptools, функция elide, аргумент shift.

к сожалению, не получилось таким образом.

Зато получилось вот так (нужны пакеты rgdal/maptools):
cnt <- readShapePoly("myshapefile", proj4string=CRS("+proj=longlat +ellps=WGS84"))
cnt2 <- spTransform(cnt, CRS("+proj=merc +lon_0=240w"))
plot(cnt2)