CRS（x）中的R错误：PROJ4参数-值对必须以+

Question

将色彩空间图像读入变量

img <- stack('BlueMarbleJuly.tif')

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

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

我得到一个错误

CRS（x）中的错误： PROJ4参数值对必须以+开头：PROJCRS [“ WGS 84 / Pseudo-Mercator”， BASEGEOGCRS [“ WGS 84”， DATUM [“世界大地测量系统1984”， ELLIPSOID [“ WGS 84”，6378137,298.257223563， LENGTHUNIT [“ metre”，1]]]， PRIMEM [“格林威治”，0， ANGLEUNIT [“ degree”，0.0174532925199433]]， ID [“ EPSG”，4326]]， CONVERSION [“大众可视化伪墨卡托”， 方法[“大众可视化伪墨卡托”， ID [“ EPSG”，1024]]， PARAMETER [“自然起源的纬度”，0， ANGLEUNIT [“ degree”，0.0174532925199433]， ID [“ EPSG”，8801]]， PARAMETER [“自然起源的经度”，0， ANGLEUNIT [“ degree”，0.0174532925199433]， ID [“ EPSG”，8802]]， PARAMETER [“虚假东移”，0， LENGTHUNIT [“ metre”，1]， ID [“ EPSG”，8806]]， PARAMETER [“ False northing”，0， LENGTHUNIT [“ metre”，1]， ID [“ EPSG”，8807]]]， CS [Cartesian，2]，

这是什么意思？

Answer 1

您正在为projectRaster提供坐标参考系统的“ wkt”表示形式，但是它期望使用“ proj4”表示形式。要解决此问题，您可以执行

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