Question

我正在将来自 Bio-Oracle 的数据用于 SDM，并且正在努力更改坐标系以匹配我的观察结果。当前的 crs 是“ces”，我想要“longlat”。我是栅格数据的新手，因此尝试了一些解决方案，但收到以下错误消息：

# set up projection parameter for use throughout script
projection <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0" 

# Load equal area rasters
layercodes <- c("BO_ph", "MS_bathy_5m", 
                "BO_sstrange", "BO_sstmean", "BO_salinity")
env <- load_layers(layercodes, equalarea = TRUE,
                   datadir = "data/")
plot(env) # global
extent(env)
# class      : Extent 
# xmin       : -17402530 
# xmax       : 17401470 
# ymin       : -7378770 
# ymax       : 7377230 

crs(env)
# CRS arguments:
#   +proj=cea +lat_ts=30 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m
# +no_defs 

# Project to longlat
env.proj <- spTransform(env, 
                   crs(projection))

<块引用>

（函数（类，fdef，mtable）中的错误：无法为签名“RasterStack”、“CRS””找到函数“spTransform”的继承方法

env.proj <- projectRaster(env, 
                          crs = projection)

<块引用>

if (maxy == miny) { 中的错误：缺少需要 TRUE/FALSE 的值另外：警告信息：在 rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 108个投影点不是有限的

env.proj <- lapply(env, function(r) projectRaster(env, projection))
# ?projectRaster

<块引用>

h(simpleError(msg, call)) 中的错误：在为函数“范围”选择方法时评估参数“x”时出错：需要 TRUE/FALSE 时缺少值另外：警告信息：在 rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 108个投影点不是有限的

我错过了什么？

更改栅格堆栈上的投影

0 个答案: