CRS(x) 中的 R 错误：PROJ4 参数-值对必须以 + 开头答案

【问题标题】：R Error in CRS(x): PROJ4 argument-value pairs must begin with +CRS(x) 中的 R 错误：PROJ4 参数-值对必须以 + 开头
【发布时间】：2020-09-27 00:56:02
【问题描述】：

将色彩空间图像读入变量

img <- stack('BlueMarbleJuly.tif')

将栅格投影到 nes 坐标系中。罗宾逊计划

img.merc <- projectRaster(img, crs = st_crs(3857)[[2]])

我得到一个错误

CRS(x) 中的错误： PROJ4 参数-值对必须以 + 开头：PROJCRS["WGS 84 / Pseudo-Mercator", BASEGEOGRCS["WGS 84", DATUM[“1984年世界大地测量系统”，椭圆体["WGS 84",6378137,298.257223563, LENGTHUNIT["米",1]]], PRIMEM["格林威治",0, ANGLEUNIT["度",0.0174532925199433]], ID["EPSG",4326]], CONVERSION["流行的可视化伪墨卡托", METHOD["流行的可视化伪墨卡托", ID["EPSG",1024]], PARAMETER["自然起源纬度",0, ANGLEUNIT["度",0.0174532925199433], ID["EPSG",8801]], PARAMETER["自然原点经度",0, ANGLEUNIT["度",0.0174532925199433], ID["EPSG",8802]], PARAMETER["假东移",0, LENGTHUNIT["米",1], ID["EPSG",8806]], PARAMETER["北偏",0, LENGTHUNIT["米",1], ID["EPSG",8807]]], CS[笛卡尔,2],

这是什么意思？

【问题讨论】：

标签： r r-raster

【解决方案1】：

您正在为projectRaster 提供坐标参考系的“wkt”表示，但它需要“proj4”表示。要解决此问题，您可以这样做

img.merc <- projectRaster(x, crs = "+init=epsg:3857")

如

library(raster)
f <- system.file("external/test.grd", package="raster")
x <- raster(f)
merc <- projectRaster(x, crs = "+init=epsg:3857")

【讨论】：

试一试，但将新 CRS 分配给 img 的操作给出了相同的错误：CRS(x) 中的错误：PROJ4 参数-值对必须以 + 开头：..."
我添加了一个示例，以便您可以看到它的工作原理。