2017-05-08 50 views
0

ためINTERP私は点のみで評価機能付き以下グリッド有する:アキマパッケージ線形補間を使用して:非常に不規則グリッド

をIトラブル私ができるように、データX、Y、およびZの設定を有していuse

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)) 

enter image description here

+0

再現可能な例をご紹介ください。 – Axeman

+0

@Axeman okは再現可能な例を追加しました。補間やzのinterpを使用しようとしました。 – robertevansanders

答えて

1

私はあなたがinterppないinterpを探していると思います。

+0

これは私が行列を設定する必要はなく、補間を行うための実際の凸集合を持つことができますか? – robertevansanders

+0

正しく覚えていれば、それに座標と値のセットを投げれば、戻り値として素晴らしいグリッド値が返されます(もちろん、グリッドを指定する必要があります)。 – Lyngbakr

+0

これは廃止予定で、interppを使用しています。それは私が使っていたものなので、奇妙なwarninigのように思える。 – robertevansanders

関連する問題