Невозможно построить слой sf поверх ggmap

Тема некоторых предыдущих вопросов - например. это. Я не могу заставить слой sf отображаться поверх слоя ggmap.

require(ggmap)
require(raster)
require(sf)

target_bbox <- c(left = -0.65, bottom = 51.2, right = 0.45, top = 51.8)
map = get_stamenmap(target_bbox, zoom = 9, maptype = 'toner-2010')

# can test this with: ggmap(map)

e = extent(list(x = c(-0.65, 0.45), y = c(51.2, 51.8)))
r = raster(volcano)
crs(r) = CRS("+proj=longlat +datum=WGS84")
p = raster::rasterToPolygons(r)
sf = st_as_sf(p) %>% st_transform(3857)

# test with: ggplot() + geom_sf(data = sf, aes(fill = layer), col = NA)

ggmap(map) + 
   geom_sf(data = sf, aes(fill = layer), col = NA, alpha = 0.4, inherit.aes = FALSE)
Coordinate system already present. Adding new coordinate system, which will replace the existing one.

Конфликт системы координат, но ошибок нет, но слой не виден. Любые предложения высоко ценятся.

Это решение github — извлечение ограничивающей рамки из объекта ggmap с помощью пользовательской функции ggmap_bbox() — успешно добавляет еще одно предупреждение к приведенному выше:

Warning message:
In st_is_longlat(x) :
  bounding box has potentially an invalid value range for longlat data

r ggmap sf

24.12.2020 18:28

Стоит ли изучать PHP в 2023-2024 годах?

Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...

Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией

В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.

Приемы CSS-макетирования - floats и Flexbox

Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...

Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest

В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...

Концепция локализации и ее применение в приложениях React ⚡️

Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...

Пользовательский скаляр GraphQL

Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...

953

Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Похоже, вы пытаетесь нарисовать объект вулкана sf над Лондоном, но вы определили его долготу/широту 0,0. Кроме этого, ваш вызов ggmap, а затем geom_sf работает нормально. Проблема заключалась только в том, что объект sf не находился в ограничивающей рамке ggmap.

Мне пришлось изменить его на «пространственный» объект, чтобы переместить его с помощью maptools::elide, а затем обратно на объект «sf», так как я не мог найти простой способ сдвига полигонов sf. Сдвиг является приближением глазного яблока.

sf_moved <- sf %>% 
  st_transform(3857) %>% # <- may be unnecessary
  as('Spatial') %>% 
  maptools::elide(.,shift = c(-80000, 6650000)) %>% 
  as('sf') %>% 
  st_set_crs(3857) %>% 
  st_transform(4326) 

ggmap(map) + 
  geom_sf(data = sf_moved, inherit.aes = FALSE, 
          col = NA, aes(fill = layer), alpha = .4) + 
  scale_fill_viridis_c()

О, Боже. Ошибка новичка - я забыл присвоить экстент растровому объекту extent(r) = e

— 24.12.2020 22:33

Бывает со всеми. Назначение степени было бы намного проще, чем мой метод!

— 24.12.2020 22:41

24.12.2020 22:29