如何在R?
中绘制球面坐标(r,theta和phi)是否可以使用persp()
功能?
我有一个积分和一个网格。
答案 0 :(得分:3)
只要你知道your math,这是一个相当微不足道的转变:
spher_to_cart <- function(r, theta, phi) list(x=r*cos(phi)*sin(theta),
y=r*sin(theta)*sin(phi),
z=r*cos(theta))
#An example dataset
data <- data.frame(r=1:10,
theta = seq(0,2*pi,length=10),
phi = seq(2*pi, 0,length=10))
spher_to_cart(data$r, data$theta, data$phi)
$x
[1] 0.000000e+00 9.848078e-01 5.130302e-01 -1.732051e+00 -1.606969e+00 1.928363e+00 3.031089e+00 -1.368081e+00 -4.431635e+00 -2.449294e-15
$y
[1] 0.0000000 -0.8263518 -2.9095389 -3.0000000 -0.5848889 -0.7018667 -5.2500000 -7.7587705 -3.7185832 0.0000000
$z
[1] 1.0000000 1.5320889 0.5209445 -2.0000000 -4.6984631 -5.6381557 -3.5000000 1.3891854 6.8944000 10.0000000
注意以弧度而不是度数使用theta和phi的值
然后,您可以使用plot3d
绘制包rgl
,例如:
s <- spher_to_cart(data$r, data$theta, data$phi)
library(rgl)
plot3d(s$x,s$y,s$z)