投影rasterLayer以使用R中的terrain（）函数计算斜率

Question

我想要投射以下栅格图层，因此我可以使用terrain()函数计算斜率。这不是真正的空间物体，只是理论表面，所以我需要一些简单的投影。 我尝试在代码的第二部分执行此操作但由于某种原因我收到错误消息，我不知道为什么会这样。

location <- 0       # location/distance parameter
s <- 1          # scale
inc <- 0.1  # increment for sequencing y axis
increment <- rev(seq(-6,6,inc))   # Create a sequence of x values (reverse)

y <- 1/(1+ exp(-(increment/s)))

# Create matrix:
rep <- length(increment)
valuematrix <- replicate(rep,y)

# Create raster
raster <- raster(valuematrix, xmn=0,xmx=13,ymn=0,ymx=13)

分配投影和计算斜率：

newproj <-  CRS("+init=epsg:4326 +proj=longlat +ellps=WGS84")

proj4string(raster) <- newproj

terrain(raster,alt='slope')

错误讯息：

Error in isLonLat(out) : 
  4 arguments passed to .Internal(nchar) which requires 3