如何用Python生成指定同心圆环数的平滑线框球体并标记交点？

要实现指定数量的线框球体并标记交点，核心思路是给每个圆环足够的采样点保证平滑，而非增加圆环总数。以下是基于Matplotlib的实现代码，完全匹配你3个水平圆环、22个垂直圆环的需求：

import matplotlib.pyplot as plt
import numpy as np

# 配置参数
sphere_radius = 1
num_horizontal_rings = 3  # 水平圆环数
num_vertical_rings = 22   # 垂直圆环数
ring_samples = 100        # 每个圆环的采样点数量（控制平滑度）
point_size = 10           # 交点标记大小

# 生成水平圆环（纬度圈）
# 这里选择3个典型纬度：赤道、北纬45°、南纬45°，可根据需求调整theta值
theta_list = np.array([0, np.pi/4, -np.pi/4])  # 纬度角度（弧度）

# 生成垂直圆环（经度圈）
phi_list = np.linspace(0, 2*np.pi, num_vertical_rings, endpoint=False)

# 创建3D绘图对象
fig = plt.figure(figsize=(8,8))
ax = fig.add_subplot(111, projection='3d')
ax.set_box_aspect([1,1,1])  # 保证球体比例正确

# 绘制水平圆环
for theta in theta_list:
    phi = np.linspace(0, 2*np.pi, ring_samples)
    x = sphere_radius * np.cos(phi) * np.cos(theta)
    y = sphere_radius * np.sin(phi) * np.cos(theta)
    z = np.full_like(phi, sphere_radius * np.sin(theta))
    ax.plot(x, y, z, color='gray', linewidth=1)

# 绘制垂直圆环
for phi in phi_list:
    theta = np.linspace(-np.pi/2, np.pi/2, ring_samples)
    x = sphere_radius * np.cos(phi) * np.cos(theta)
    y = sphere_radius * np.sin(phi) * np.cos(theta)
    z = sphere_radius * np.sin(theta)
    ax.plot(x, y, z, color='gray', linewidth=1)

# 计算并标记所有交点
intersection_x = []
intersection_y = []
intersection_z = []
for theta in theta_list:
    for phi in phi_list:
        x = sphere_radius * np.cos(phi) * np.cos(theta)
        y = sphere_radius * np.sin(phi) * np.cos(theta)
        z = sphere_radius * np.sin(theta)
        intersection_x.append(x)
        intersection_y.append(y)
        intersection_z.append(z)

ax.scatter(intersection_x, intersection_y, intersection_z, color='red', s=point_size)

# 隐藏坐标轴刻度，优化视觉效果
ax.set_xticks([])
ax.set_yticks([])
ax.set_zticks([])

plt.show()

关键说明： #

• 平滑度控制：通过ring_samples参数给每个圆环设置足够多的采样点（比如100个），单个圆环就能呈现平滑曲线，无需依赖大量圆环。
• 圆环位置调整：theta_list可以自定义水平圆环的纬度位置，比如如果需要对应六元环糖的结构，可替换为符合需求的角度值。
• 交点标记：通过遍历所有水平/垂直圆环的角度组合，计算球面坐标交点并统一标记，确保每个交叉点都被精准定位。
• 视觉优化：设置box_aspect保证球体比例正确，隐藏坐标轴刻度让焦点集中在线框和交点上。

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