パッケージを使用してplot3D、4 番目の変数に基づく色で xyz データをプロットしようとしています。
私のデータは次のようになります。
data.frame(xx,yy,zz,Ammpl)
xx yy zz Ammpl
1 63021.71 403105.0 1.181028516 1170
2 63021.71 403105.0 0.977028516 1381
3 63021.71 403105.0 0.861028516 807
4 63021.71 403105.0 0.784028516 668
5 63021.71 403105.0 0.620028516 19919
6 63021.71 403105.0 0.455028516 32500
7 63021.71 403105.0 0.446028516 32500
8 63021.71 403105.0 0.436028516 32500
9 63021.71 403105.0 0.426028516 32500
10 63021.71 403105.0 0.281028516 17464
私が最初にしたことは、パッケージakimaから関数を使用することでしたinterp。通常のグリッドを作成するには
s100<-interp(xx,yy,zz, xo=seq(min(xx):max(xx), length=30, yo=seq(min(yy):max(yy), length=30), duplicate="mean")
パッケージからplot3D関数を使用しましたslice3D
slice3D(x = s100$x, y = s100$y, z = s100$z, colvar = ampl,
ys = 1:10, zs = NULL, NAcol = "black",
expand = 0.4, theta = 45, phi = 45)
ampl私の問題は、データを 作成する方法がわからないcolvarため、このプロットで色付けを提供することです。
いろいろ試しましたが、ここで行き詰まると思います。
from をexample(slice3D)使用して を作成します。oxsat$valdata(oxsat)colvar
> dim(Oxsat$val)
[1] 180 90 33