r - 将坐标引用系统设置为在 R 中创建的多边形

Question

我正在尝试在使用 ggplot2::geom_sf() 绘制的 map 中创建放大正方形(多边形)因此，我需要将多边形作为一个简单的特征，以便我可以将正方形和 map 一起绘制。我可以使用以下代码创建多边形。

zoom_square_1<- rbind(c(-10.821079,-68.403542),
                  c(-10.821079,-68.367060),
                  c(-10.862918,-68.367060),
                  c(-10.862918,-68.403542),
                  c(-10.821079,-68.403542)) # add the first point again to it wrap up.
zoom_square_1_pol <- st_polygon(list(zoom_square_1))
plot(zoom_square_1_pol) # it returns the polygon.

但是，当我尝试为其分配坐标引用系统 (CRS) 时，它不起作用。这些是我到目前为止尝试过的代码。

sp::proj4string(zoom_square_1_pol) <- CRS("+proj=longlat +datum=WGS84 +no_defs")
zoom_square_1_pol_ex1 <- st_as_sf(x = zoom_square_1_pol[[1]],
                              crs = "+proj=longlat +datum=WGS84 +no_defs")
zoom_square_1_pol_ex2 <- st_transform(x = zoom_square_1_pol,
                                  crs = "+proj=longlat +datum=WGS84 +no_defs")

欢迎任何帮助。

在 Chris 发表评论后，我实际上可以使用 ggplot2::geom_sf() 绘制多边形，并意识到我的纬度和经度坐标自开始以来就处于错误的位置。所以，这是代码的最终版本:

 zoom_square_1<- rbind(c(-68.403542,-10.821079),
                       c(-68.367060,-10.821079),
                       c(-68.367060,-10.862918),
                       c(-68.403542,-10.862918),
                       c(-68.403542,-10.821079)) #add the first point again to it wrap up.
zoom_square_1_pol <- st_polygon(list(zoom_square_1))
zoom_square_1_pol_CRS<- st_sfc(x = zoom_square_1_pol,
                               crs = "+proj=longlat +datum=WGS84 +no_defs")

Answer 1

我能够通过将正方形作为列表提供给 st_as_sfc 来设置 crs:

res <- st_as_sfc(list(zoom_square_1_pol), 
                 crs = "+proj=longlat +datum=WGS84 +no_defs")