前言

本篇学习使用numpy生成球面斐波那契网格。

斐波那契网格公式

对于一个球体，要在球体表面上均匀选取N个点，就可以使用斐波那契网格公式：

其中，表示球半径，即黄金分割比例。

代码

通过numpy形成斐波那契网格点云，然后通过open3d来可视化点云图：

import open3d as o3d
import cv2
import numpy as np


if __name__ == '__main__':
    vis = o3d.visualization.Visualizer()
    N = 1000
    r = 1
    phi = (np.sqrt(5) - 1) / 2
    n = np.arange(0, N)
    z = ((2*n + 1) / N - 1)
    x = (np.sqrt(1 - z**2)) * np.cos(2 * np.pi * (n + 1) * phi)
    y = (np.sqrt(1 - z**2)) * np.sin(2 * np.pi * (n + 1) * phi)
    print(z.shape, x.shape, y.shape)
    pcd = r * np.stack([x, y, z], axis=-1)
    print(pcd.shape)
    point_cloud = o3d.geometry.PointCloud()
    point_cloud.points = o3d.utility.Vector3dVector(pcd)
    vis.add_geometry(point_cloud)
    o3d.visualization.draw_geometries([point_cloud])

取100个点的效果：