求Matlab中sphere()函数的R语言等效实现方法

Matlab的sphere(n)函数会生成单位球面的三个坐标矩阵X、Y、Z，每个矩阵维度为(n+1)×(n+1)，对应将球面分成n段的网格点坐标。要在R中实现[X,Y,Z]=sphere(200)的等效效果，有两种常用方式：

一、手动实现（完全匹配Matlab输出格式） #

直接通过三角函数生成球面坐标网格，代码如下：

n <- 200
# 生成极角（0到π）和方位角（0到2π）的网格点
phi <- seq(0, pi, length.out = n + 1)
theta <- seq(0, 2 * pi, length.out = n + 1)

# 计算对应坐标矩阵
X <- outer(cos(theta), sin(phi))
Y <- outer(sin(theta), sin(phi))
Z <- outer(rep(1, n + 1), cos(phi))

说明： #

• phi和theta分别对应球面坐标系的极角和方位角，分成n+1个点以匹配Matlab的分段逻辑
• outer()函数用于生成二维网格矩阵，等价于Matlab的meshgrid操作
• 最终生成的X、Y、Z均为201×201的矩阵，和Matlab输出的维度、坐标值完全一致

二、使用rgl包（适合可视化场景） #

如果需要结合3D可视化，可以用rgl包的sphere3d()函数，若要提取坐标数据：

library(rgl)
n <- 200
# 生成球面对象，subdivision参数对应分段数
sphere_obj <- sphere3d(subdivision = n - 1)
# 提取坐标数据（注意：这里的坐标是顶点数组，结构和Matlab的网格矩阵不同）
X <- sphere_obj$vb[1,,]
Y <- sphere_obj$vb[2,,]
Z <- sphere_obj$vb[3,,]

注意： #

rgl返回的坐标是三角化的顶点数据，并非Matlab那样的规则网格矩阵，若需要和原Matlab代码的矩阵格式完全兼容，优先选择手动实现的方法。

内容的提问来源于stack exchange，提问作者Bensstats