剧情从2D函数中的R

Question

一个3D线性表面我有以下数据：

x <- c(1000,2000,3000,4000,5000,6000,8000,12000) 

y_80 <- c(33276,33276,5913,2921,1052,411,219,146) 
y_60 <- c(14724,14724,3755,1958,852,372,211,140) 
y_40 <- c(9632,9632,2315, 1250,690,332,196,127) 
y_20 <- c(4672,4672,1051,562,387,213,129,81) 
y_5 <- c(825,825,210,118,88,44,27,17) 

从该数据创建5的样条函数：

f_80 <- splinefun(x, y_80, method=c("monoH.FC")) 
f_60 <- splinefun(x, y_60, method=c("monoH.FC")) 
f_40 <- splinefun(x, y_40, method=c("monoH.FC")) 
f_20 <- splinefun(x, y_20, method=c("monoH.FC")) 
f_5 <- splinefun(x, y_5, method=c("monoH.FC")) 

在Z轴是后的数f_即函数f_80，Z值为80，依此类推。

我需要做的是从这些函数绘制三维曲面，并在线条之间进行线性插值。 R中可能吗？ 我已经有过类似的问题，但我找不到答案。谢谢

Answer 1

有各种R伪3D绘图功能。基本图形有persp,wireframe是在格子中，最酷的是，surface3d在rgl中。它们在处理轴参数方面有所不同，但它们通常都接受z值的矩阵参数。 `persp是简单的：

z_80 <- f_80(x) 
z_60 <- f_60(x) 
z_40 <- f_40(x) 
z_20 <- f_20(x) 
z_5 <- f_5(x) 

png(); persp(matrix(c(z_80, z_60, z_40, z_20, z_5), nrow=length(x), 
       dimnames= list(X=x, Y=c(80,60,40,20,5)))) 
dev.off() 

你可能会发现，使用ticktype="detailed"是更多的信息：

png(); persp(matrix(c(z_80,z_60,z_40,z_20,z_5), nrow=length(x), 
         dimnames=list(X=x, Y=c(80,60,40,20,5))) , 
         ticktype="detailed", xlab="X axis", ylab="Y axis", 
         zlab=".    Z= f(y)", theta=45) ; dev.off() 

我通常需要观察角度theta得到的东西我喜欢玩。让z轴标签远离长z标签标签是一种破解。