Как изменить проекцию proj = laea на longlat, используя spTransform, не деформируя карту? - PullRequest
Новая
       79

Как изменить проекцию proj = laea на longlat, используя spTransform, не деформируя карту?

0 голосов
/ 20 апреля 2020

У меня есть набор данных с широтой / долготой, который я собираюсь нанести на карту рек АН. У меня проблема с проекцией карты. Шейп-файл, который я использую, является proj = laea, и при построении с помощью ggplot оси находятся в UTM (?), А не в широте / долготе. Я не смог понять, как изменить проекцию на что-то, что будет отображаться с широтой / долготой. Прямо сейчас, если я попытаюсь изменить проекцию (скажем, "+ init = epsg: 4269 + proj = longlat + ellps = GRS80 + datum = NAD83 + no_defs + towgs84 = 0,0,0"), то он закручивает карту. 

rivers <- readOGR(dsn=path.expand("shapefile/data/hydrography_l_rivers_v2.shp"))

proj4string(rivers)
[1] "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs"

rivers<-spTransform(rivers, " +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997+units=m +no_defs")

ggplot()+
  geom_line(data=rivers, aes(x=long, y=lat, group=group))
...