r - 添加投影时，向 geom_sf 形状添加点不起作用

Question

当我尝试将点图层添加到geom_sf()图层和投影时，这些点似乎最终位于德克萨斯州南部的一个位置。以下是重现此问题的最小示例。

library(sf)
library(ggplot2)
# devtools::install_github("hrbrmstr/albersusa")
library(albersusa)
crs_use = "+proj=laea +lat_0=30 +lon_0=-95"

d_points = data.frame(long = c(-110, -103, -84), 
                      lat  = c(45, 40, 41))

A = ggplot(data = usa_sf()) +
  geom_sf() + 
  geom_point(data = d_points, 
             aes(x = long, y = lat), 
             color = "red", size = 5) + 
  theme_minimal() + 
  ggtitle("(A) right point position, wrong projection")

B = ggplot(data = usa_sf()) +
  geom_sf() + 
  geom_point(data = d_points, 
             aes(x = long, y = lat), 
             color = "red", size = 5) + 
  coord_sf(crs = crs_use) + 
  theme_minimal() + 
  ggtitle("(B) right projection, wrong points using geom_point()")

C = ggplot() + 
  geom_sf(data = usa_sf()) + 
  geom_sf(data = st_as_sf(d_points,
                          coords = c("long", "lat"), crs = crs_use), 
          color = "red", size = 5) + 
  coord_sf(crs = crs_use) + 
  theme_minimal() + 
  ggtitle("(C) right projection, wrong points using geom_sf() points")

cowplot::plot_grid(A, B, C, nrow = 3)

我想使用自定义投影向美国地图添加一个点图层。但是，每当我使用投影时，我指定的点就会变成德克萨斯州东南部的一个奇怪位置，它不在我指定的位置。

任何有关解决此问题的建议表示赞赏。谢谢！

标签： rggplot2data-visualizationgeospatialsf

Answer 1

您还需要重新规划您的观点；sf::st_transform()应该做的工作。

我无法访问 {albersusa} 包，所以我使用 {USABoundaries} 来获取较低 48 的地图，但这不是重点；st_transform(crs = crs_use)在绘图之前，该点适用于您的两个空间对象。

library(sf)
library(ggplot2)
library(USAboundaries)

crs_use <- "+proj=laea +lat_0=30 +lon_0=-95"

usa_sf <- us_boundaries(type="state", resolution = "low") %>% 
  dplyr::filter(!state_abbr %in% c("PR", "AK", "HI")) %>% 
  st_transform(crs = crs_use)

d_points <- data.frame(long = c(-110, -103, -84), 
                      lat  = c(45, 40, 41)) %>% 
  st_as_sf(coords = c("long", "lat"), crs = 4326) %>% 
  st_transform(crs = crs_use)

ggplot(data = usa_sf) +
  geom_sf() + 
  geom_sf(data = d_points, 
             color = "red", size = 5) + 
  theme_minimal()