r - 获得R上的花键表面

Question

我如何生成b样条曲面，比如说:

x=attitude$rating
y=attitude$complaints
z=attitude$privileges

对于样条曲线将为x和y。 z是控制点集。

Answer 1

据我了解，您有x，y和z数据，并且想在x和y上使用双变量样条插值，并使用z作为控制点。您可以使用interp(...)包中的akima来执行此操作。

library(akima)
spline <- interp(x,y,z,linear=FALSE)
# rotatable 3D plot of points and spline surface
library(rgl)
open3d(scale=c(1/diff(range(x)),1/diff(range(y)),1/diff(range(z))))
with(spline,surface3d(x,y,z,alpha=.2))
points3d(x,y,z)
title3d(xlab="rating",ylab="complaints",zlab="privileges")
axes3d()

由于x，y和x是高度相关的，因此图本身对您的数据集就没有兴趣。

编辑OP的评论。

如果要使用b样条曲面，请尝试不幸地命名为 mba.surf(...)包中的 MBA。

library(MBA)
spline <- mba.surf(data.frame(x,y,z),100,100)

library(rgl)
open3d(scale=c(1/diff(range(x)),1/diff(range(y)),1/diff(range(z))))
with(spline$xyz,surface3d(x,y,z,alpha=.2))
points3d(x,y,z)
title3d(xlab="rating",ylab="complaints",zlab="privileges")
axes3d()