在 R 中从 Behrmann CRS 映射转换为 +proj=longlat WGS84

Question

我有一个 200 公里 x 200 公里的网格单元形状文件，覆盖贝尔曼等面积圆柱投影中的世界陆地区域。我的目标是将 shapefile 转换为 +proj=longlat WGS84 格式，以便我可以将其与常用投影中的地图相匹配，例如 maptools 中的 wrld_simpl。但是，我没有成功，希望能得到一些帮助。

rm(list = ls())
library(RCurl)
library(raster)
library(maptools)
library(rgdal)
data("wrld_simpl")

tmp <- tempfile() download.file("https://github.com/darunabas/extras/blob/master/temp_shapefile.zip?raw=true", destfile = tmp)
unzip(tmp, exdir = ".")

s <- rgdal::readOGR("temp_shapefile")
proj4string(s) = CRS("+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +datum=WGS84 +ellps=WGS84 +units=m +no_defs")

p <- spTransform(s, CRS("+proj=longlat +datum=WGS84"))

我收到以下错误： 检测到非有限变换： [,1] [,2] [,3] [,4] .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, 中的错误： 多边形失败 1106 多边形 1 分 另外：警告信息： 在 .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, : 2 个投影点不是有限的

Answer 1

我在地理世界上不是很有经验，但这可能会有所帮助：

library( sf )
sf <- read_sf( "./temp_shapefile.shp")
st_crs( sf ) <- "+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +datum=WGS84 +ellps=WGS84 +units=m +no_defs"
sf2 <- st_transform( sf, crs = "+proj=longlat +datum=WGS84" )