Início / coding / Perguntas / 77244364
Accepted
geotheory
Asked: 2023-10-06 20:16:34 +0800 CST

O ggmap/sf ainda está plotando o ponto no lugar errado?

  • 772

Esse bug antigo ainda não foi resolvido ou já existe uma solução alternativa razoável?

Algumas cidades:

require(ggplot2)
require(ggmap)
require(sf)

cities = sf::read_sf('{
  "type": "FeatureCollection",
  "name": "cities",
  "crs": {
    "type": "name",
    "properties": {
      "name": "urn:ogc:def:crs:OGC:1.3:CRS84"
    }
  },
  "features": [{
      "type": "Feature",
      "properties": {
        "city": "London"
      },
      "geometry": {
        "type": "Point",
        "coordinates": [-0.13, 51.51]
      }
    },
    {
      "type": "Feature",
      "properties": {
        "city": "Rome"
      },
      "geometry": {
        "type": "Point",
        "coordinates": [12.48, 41.89]
      }
    },
    {
      "type": "Feature",
      "properties": {
        "city": "Stockholm"
      },
      "geometry": {
        "type": "Point",
        "coordinates": [18.07, 59.33]
      }
    },
    {
      "type": "Feature",
      "properties": {
        "city": "Istanbul"
      },
      "geometry": {
        "type": "Point",
        "coordinates": [28.96, 41.01]
      }
    }
  ]
}
')

E um mapa base ggmap:

bb = c(left = -11, bottom = 35, right = 42, top = 65)
europe = ggmap::get_stamenmap(bbox = bb, maptype = 'toner-lite', zoom = 3)

ggmap(europe) + 
  geom_sf(data = cities, inherit.aes = FALSE, col = 'red', size = 3)

insira a descrição da imagem aqui

1 respostas

  1. Best Answer

    Acho que o problema é que ggmapusa o objeto coordenado CoordMappara decidir onde colocar os pixels de um raster. Ao mudar para coord_sf(o que acontece quando você adiciona uma geom_sfcamada), você está desalinhando o raster. Se você extrair as coordenadas dos seus pontos, poderá simplesmente sobrepô-los como uma geom_pointcamada:

    ggmap(europe) + 
  geom_point(data = as.data.frame(st_coordinates(cities)), aes(X, Y), 
             inherit.aes = FALSE, col = 'red', size = 3)

    insira a descrição da imagem aqui

    Se você quiser desenhar sfobjetos mais complexos, eu ficaria longe ggmape usaria maptiles:

    library(maptiles)
library(tidyterra)

bb <- c(left = -11, bottom = 35, right = 42, top = 65)

europe <- matrix(bb, 2, byrow = TRUE) |>
  st_multipoint()       |> 
  st_sfc(crs = 4326)    |>
  st_transform(3857) |>
  get_tiles(provider = "Stadia.Stamen.TonerLite", zoom = 4, crop = TRUE)

ggplot() +
  geom_spatraster_rgb(data = europe) +
  geom_sf(data = cities, size = 3, color = "red") +
  coord_sf(crs = 3857, expand = FALSE, ylim = st_bbox(europe)[c(2, 4)])

    insira a descrição da imagem aqui

    • 1

relate perguntas