基于φ和θ范围用Open3D创建部分球体的技术咨询

方法一：Open3D原生手动生成部分球体 #

你可以通过球坐标转笛卡尔坐标的公式，手动生成指定φ、θ范围内的顶点和三角面，再构建TriangleMesh。具体实现如下：

import open3d as o3d
import numpy as np

def create_partial_sphere(radius=1.0,
                         phi_start=0.0, phi_end=np.pi/2,  # 极角范围，示例为0到π/2（上半球）
                         theta_start=0.0, theta_end=np.pi, # 方位角范围，示例为0到π（半圈）
                         phi_steps=30, theta_steps=60):
    # 生成φ和θ的网格采样点
    phi = np.linspace(phi_start, phi_end, phi_steps)
    theta = np.linspace(theta_start, theta_end, theta_steps)
    phi_grid, theta_grid = np.meshgrid(phi, theta)

    # 球坐标转笛卡尔坐标
    x = radius * np.sin(phi_grid) * np.cos(theta_grid)
    y = radius * np.sin(phi_grid) * np.sin(theta_grid)
    z = radius * np.cos(phi_grid)

    # 整理顶点数组
    vertices = np.stack([x.flatten(), y.flatten(), z.flatten()], axis=1)

    # 生成三角面索引
    faces = []
    for i in range(theta_steps - 1):
        for j in range(phi_steps - 1):
            # 第一个三角形
            v1 = i * phi_steps + j
            v2 = (i + 1) * phi_steps + j
            v3 = i * phi_steps + (j + 1)
            faces.append([v1, v2, v3])
            # 第二个三角形
            v1 = (i + 1) * phi_steps + j
            v2 = (i + 1) * phi_steps + (j + 1)
            v3 = i * phi_steps + (j + 1)
            faces.append([v1, v2, v3])

    # 创建TriangleMesh并计算法向量
    mesh = o3d.geometry.TriangleMesh()
    mesh.vertices = o3d.utility.Vector3dVector(vertices)
    mesh.triangles = o3d.utility.Vector3iVector(faces)
    mesh.compute_vertex_normals()
    return mesh

# 使用示例：创建1/4上半球（φ:0-π/2，θ:0-π/2）
partial_sphere = create_partial_sphere(phi_end=np.pi/2, theta_end=np.pi/2)
o3d.visualization.draw_plotly([partial_sphere])

方法二：使用trimesh库快速生成 #

如果觉得手动构建网格繁琐，可以用trimesh库实现更灵活的部分球体生成，步骤更简洁：

首先安装依赖：

pip install trimesh

示例代码：

import trimesh
import open3d as o3d
import numpy as np

# 方式1：通过裁剪完整球体得到指定范围的部分球
full_sphere = trimesh.creation.sphere(radius=1.0, subdivisions=4)
# 保留z>0且x>0的区域（对应φ:0-π/2，θ:-π/2-π/2）
partial_sphere_trimesh = full_sphere.slice_plane([1,0,0], 0).slice_plane([0,0,1], 0)

# 转换为Open3D格式（可选）
mesh_o3d = o3d.geometry.TriangleMesh()
mesh_o3d.vertices = o3d.utility.Vector3dVector(np.array(partial_sphere_trimesh.vertices))
mesh_o3d.triangles = o3d.utility.Vector3iVector(np.array(partial_sphere_trimesh.faces))
mesh_o3d.compute_vertex_normals()

o3d.visualization.draw_plotly([mesh_o3d])

补充说明 #

• 手动生成方式完全基于Open3D，无额外依赖，适合需要精确控制φ/θ范围的场景
• trimesh提供了更简洁的API和丰富的网格操作，适合快速开发需求

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