我有以下格,僅在點評價函數:使用阿克瑪包線性內插:interp的一個非常不規則網格
我無法設置數據的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
@Axeman ok添加了一個可重現的例子,試圖在z上插入或使用interp – robertevansanders