Início / coding / Perguntas / 78946558
Accepted
bill999
Asked: 2024-09-04 09:56:08 +0800 CST

Como obter coordenadas de latitude e longitude no mesmo sistema de coordenadas do mapa de limites estaduais

  • 772

Tenho algumas coordenadas, digamos:

library(tidyverse)
library(haven)
library(tidycensus)
library(tigris)

coords <- data.frame(lat = c(38.09720, 36.85298, 31.31517, 21.48344), long = c(-121.38785, -75.97742, -85.85522, -158.03648))

Então recebo um mapa dos EUA:

geo <- get_acs(geography = "state",
               variables = c(x = "B04006_036"),
               year = 2021, 
               geometry = TRUE, 
               keep_geo_vars=TRUE) %>%
    filter(STATEFP!="72")

#to get alaska and hawaii in the picture
geo <- shift_geometry(geo)

Em seguida, tento traçar o mapa do estado, com as coordenadas sobrepostas:

ggplot(data = coords) +
    geom_point(aes(x=lat, y=long)) +
    geom_sf(fill = "transparent", color = "gray50", size = 1, data = geo %>% group_by(STATEFP) %>% summarise()) +   
    theme(panel.background = element_rect(fill = 'white')) +
    theme(panel.grid = element_blank(),axis.title = element_blank(),
          axis.text = element_blank(),axis.ticks = element_blank(),
          panel.border = element_blank())

Isso produz:

insira a descrição da imagem aqui

No entanto, isso não funciona, pois produz um mapa com todas as coordenadas aparentemente no mesmo lugar. Como posso modificar as coisas para que o mapa e as coordenadas fiquem na mesma escala?

1 respostas

  1. Best Answer

    Isto é simplesmente um erro de projeção! Você precisa dizer ao R (mais especificamente ao sf) em qual sistema de referência de coordenadas suas coordenadas estão. Eu assumirei que elas estão no típico WGS 1984 (EPSG:4326). Então podemos criar um objeto sf a partir do seu coordsdf e plotar com geom_sf.

    library(tidyverse)
library(haven)
library(tidycensus)
library(tigris)
library(sf)

coords <- data.frame(lat = c(38.09720, 36.85298, 31.31517, 21.48344), long = c(-121.38785, -75.97742, -85.85522, -158.03648)) |>
  st_as_sf(coords = c("long", "lat"), crs = 4326) %>% # specify coordinate columns and coordinate reference system
  shift_geometry() # remember to shift the coords too!

geo <- get_acs(geography = "state",
               variables = c(x = "B04006_036"),
               year = 2021, 
               geometry = TRUE, 
               keep_geo_vars=TRUE) %>%
    filter(STATEFP!="72")

geo <- shift_geometry(geo)

ggplot() +
    geom_sf(data = coords) + # now you can use geom_sf
    geom_sf(fill = "transparent", color = "gray50", size = 1, data = geo %>% group_by(STATEFP) %>% summarise()) +   
    theme(panel.background = element_rect(fill = 'white')) +
    theme(panel.grid = element_blank(),axis.title = element_blank(),
          axis.text = element_blank(),axis.ticks = element_blank(),
          panel.border = element_blank())

    enredo correto

    • 1

relate perguntas