在R

Question

我想使用arcos变换在R中绘制一个球体，其表面上的网格线对应于球体的等面积网格化。

我一直在尝试R packakge rgl，并从以下方面获得了一些帮助： Plot points on a sphere in R

以相等的纬度长距绘制网格线。

我有以下函数，该函数返回点的数据框，这些点是我想要的网格线的交叉点，但不确定如何进行。

plot_sphere <- function(theta_num,phi_num){

  theta <- seq(0,2*pi,(2*pi)/(theta_num))

  phi <-  seq(0,pi,pi/(phi_num))

  tmp <- seq(0,2*phi_num,2)/phi_num

  phi <- acos(1-tmp)

  tmp <- cbind(rep(seq(1,theta_num),each = phi_num),rep(seq(1,phi_num),times = theta_num))

  results <- as.data.frame(cbind(theta[tmp[,1]],phi[tmp[,2]]))
  names(results) <- c("theta","phi")

  results$x <- cos(results$theta)*sin(results$phi)
  results$y <- sin(results$theta)*sin(results$phi)
  results$z <- cos(results$phi)
  return(results)
}

sphere <- plot_sphere(10,10)

任何人都可以帮忙，总的来说，我发现很难使用rgl函数。

Answer 1

如果使用lines3d或plot3d(..., type="l")，则会得到一个将数据框中的点连接起来的图。要休息一下（您不需要长行），请添加包含NA值的行。

您的plot_sphere函数中的代码似乎真的搞砸了（您两次计算phi，没有生成所请求长度的向量，等等），但是基于此函数的工作原理是：

function(theta_num,phi_num){

  theta0 <- seq(0,2*pi, len = theta_num)

  tmp <- seq(0, 2, len = phi_num)

  phi0 <- acos(1-tmp)

  i <- seq(1, (phi_num + 1)*theta_num) - 1

  theta <- theta0[i %/% (phi_num + 1) + 1]
  phi <- phi0[i %% (phi_num + 1) + 1]

  i <- seq(1, phi_num*(theta_num + 1)) - 1
  theta <- c(theta, theta0[i %% (theta_num + 1) + 1])
  phi <- c(phi, phi0[i %/% (theta_num + 1) + 1])

  results <- data.frame( x = cos(theta)*sin(phi),
                         y = sin(theta)*sin(phi),
                         z = cos(phi))
  lines3d(results)
}