2013-03-30 3 views
6

私はパラメトリック3Dカーブしているデータを持っている:私はscatterplot3dplot3dは本当にクールたRで平滑化したパラメトリック3Dカーブを描くには?

t  x  y  z 
0.000 3.734 2.518 -0.134  
0.507 2.604 9.059 0.919 
0.861 1.532 11.584 -0.248 
1.314 1.015 1.886 -0.325 
1.684 2.815 4.596 3.275 
1.938 1.359 8.015 2.873 
2.391 1.359 8.015 2.873 
.............................. 

を。しかし、私は滑らかな3D曲線が必要です。

どうすればRでプロットできますか?

+0

あなたは 'xyz'値を生成 'T'のパラメトリック機能を持っていますか?しかし、そうでない場合でも、あなたはすでに持っているポイントの間で3Dで補間したいのですが、簡単な線形補間よりも優れていますか? – Spacedman

+1

'scatterplot3d'には' type'引数があります。これはあなたのポイントをリンクして曲線を描くために '' l ''に設定することができます。 –

+0

@Spacedmanはい、単純な線形補間よりも優れたものが欲しいです。 – midas

答えて

10

splineを使用すると、ポイントを補間して曲線を滑らかにすることができます。 Spacedmanさんのコメント@を1として

d <- read.delim(textConnection(
"t x y z 
0.000 3.734 2.518 -0.134 
0.507 2.604 9.059 0.919 
0.861 1.532 11.584 -0.248 
1.314 1.015 1.886 -0.325 
1.684 2.815 4.596 3.275 
1.938 1.359 8.015 2.873 
2.391 1.359 8.015 2.873" 
), sep=" ") 
ts <- seq(from = min(d$t), max(d$t), length=100) 
d2 <- apply(d[,-1], 2, function(u) spline(d$t, u, xout = ts)$y) 
library(scatterplot3d) 
p <- scatterplot3d(d2, type="l", lwd=3) 
p$points3d(d[,-1], type="h") 

Smoothed curve

、あなたもrglを使用することができます。 これはあなたが対話的にシーンを回転させることができます。

library(rgl) 
plot3d(d2, type="l", lwd=5, col="navy") 
points3d(d[,-1]) 
spheres3d(d[,-1], radius=.1, col="orange") 
segments3d(matrix(t(cbind(d[,-1], d[,2:3], 0)), nc=3, byrow=TRUE)) 
planes3d(0,0,1,0, col="yellow", alpha=.5)  # Plane z=0 

enter image description here

+0

ブラボー!もし私ができるなら、私はこのより多くの信用を与えるだろう。 –

+1

3次元プロットのrglパッケージを試しましたか? – Spacedman

+0

これは、彼らがローラーコースターを作成する方法ですね! – PascalVKooten

関連する問題