Как выровнять ggmap CRS карты Google долготы широты - PullRequest
Новая
       35

Как выровнять ggmap CRS карты Google долготы широты

3 голосов
/ 05 февраля 2020

Несмотря на многочисленные сообщения о прогнозах CRS, et c Я не могу остановить затопление моего города.

Если я go до Google Maps и введите -35.016, 117.878 город Олбани находится на dry земле (как на следующем изображении): enter image description here

Если я введу широту / долготу в R и попытаюсь отобразить, используя пакет простых функций и ggmap Город в океане: 

library(tidyverse)
library(sf)
library(lwgeom) 
library(ggmap) 

lat <- c(-35.016)
lon <- c(117.878)
df <- tibble(lon,lat) %>% st_as_sf( coords = c("lon", "lat"), crs = 4326)

bbox_aus <- c(left = 113.338953078, bottom = -43.6345972634, right = 153.569469029, top = -10.6681857235)
ggmap_aus <- ggmap(get_stamenmap(bbox_aus, zoom = 5, maptype = "toner-background"))  

ggmap_aus + 
  geom_sf(data = df, colour = "red" , size = 3,  alpha = 0.5, inherit.aes = FALSE) +
 # coord_sf(datum = sf::st_crs(4326)) +
  labs(title = "Albany Sinking",
       x = NULL,
       y = NULL) +
   theme_bw()

enter image description here

1 Ответ

1 голос
/ 05 февраля 2020

Работает, если вы используете geom_point() с lon и lat в качестве x и y. 

df <- tibble(lon,lat) %>% st_as_sf( coords = c("lon", "lat"), crs = 4326,
                                    remove = FALSE)

ggmap_aus + 
  geom_point(data = df, colour = "red", size = 3,  alpha = 0.5,
             aes(x = lon, y = lat)) +
  # coord_sf(datum = sf::st_crs(4326)) +
  labs(title = "Albany is saved",
       x = NULL,
       y = NULL) +
  theme_bw()

enter image description here

На основании этого комментария , использование geom_point() с x и y aestheti c выравнивает больше тесно с тем, как ggmap создает ggplot.

К сожалению, я не уверен, как заставить его работать с geom_sf(), который строит график с использованием столбца geometry. В этом связанном комментарии есть некоторое обсуждение, но решение, похоже, заключается в использовании inherit.aes = FALSE, который вы уже попробовали.

Исходя из предупреждения Coordinate system already present. Adding new coordinate system, which will replace the existing one., я предполагаю, что у объекта ggmap есть некоторая система координат, которая не является 4326, но я не смог найти, как получить к ней доступ. Я попытался перепроектировать df в EPSG: 3857, но это не сработало.

...