【发布时间】:2017-10-06 15:49:42
【问题描述】:
我有下面的网格,其中一个函数只在点进行评估:
我无法设置数据 x、y 和 z,因此我可以使用
interp(x,y,z)
有谁知道这是否可能? Akima 说它支持不规则间隔的网格。
数据:
categoryvector=21
ninterp=4
points = seq(0,categoryvector,by=ceil(((categoryvector - 0)/(ninterp - 1))))
if(!categoryvector%in%points){points=c(points,categoryvector)}
if(!1%in%points){points=c(points,1); points=sort(points)}
keepers = data.table(expand.grid(prod1vint1=points,prod1vint2=points))
points1 = c(keepers$prod1vint1,0:categoryvector,rep(0,categoryvector+1),rep(categoryvector,categoryvector+1),0:categoryvector)
points2 = c(keepers$prod1vint2,rep(0,categoryvector+1),0:categoryvector,0:categoryvector,rep(categoryvector,categoryvector+1))
z =data.table(points1=points1,points2=points2,z=rnorm(length(points2),10,3))
【问题讨论】:
-
请分享一个可重现的例子。
-
@Axeman ok 添加了一个可重现的示例,尝试在 z 上插值或使用 interp
标签: r matrix grid interpolation