函数“矢量图”显示矩阵中的梯度,但不清楚如何在单个像素处得到梯度值。例如:
library(raster)
library(rasterVis)
ex <- matrix( c(1,2,3,4,2,3,4,5,3,4,5,6,4,5,6,7), nrow=4)
r=raster(t(ex[,ncol(ex):1]), xmn=0.5,xmx=nrow(ex)+.5, ymn=0.5,ymx=ncol(ex)+0.5)
projection(r)=CRS("+init=epsg:27700")
persp(r,theta=-50,phi=20, shade=0.23,col="red")
vectorplot(r, scaleSlope=.1)如何得到梯度的值?函数“矢量图”返回一个对象,但所有字段似乎都不包含与所显示的梯度相关的信息。文档中说显示是副作用,所以如果没有返回梯度信息,是否还有其他包计算梯度并将该信息返回给用户?
谢谢!
发布于 2018-03-22 06:40:17
rasterVis在raster包中使用terrain方法
您可以通过下面的代码来找出这一点:
showMethods("vectorplot")
getMethod("vectorplot", "Raster")
rasterVis:::fooSlopeAspect
rasterVis:::extractSAhttps://stackoverflow.com/questions/49420066
复制相似问题