Python numpy使用记录5.生成球面斐波那契网格

最新推荐文章于 2023-05-11 19:25:08 发布

RuiH.AI

最新推荐文章于 2023-05-11 19:25:08 发布

阅读量1k

点赞数 1

分类专栏： python相关学习文章标签： python 开发语言机器学习

本文链接：https://blog.csdn.net/qq_41035283/article/details/124689331

版权

python相关学习专栏收录该内容

41 篇文章 10 订阅

订阅专栏

Python numpy使用记录5.生成球面斐波那契网格

前言
斐波那契网格公式
代码

前言

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

斐波那契网格公式

对于一个球体，要在球体表面上均匀选取N个点，就可以使用斐波那契网格公式：
$z_n = r(\frac{2n-1}{N} - 1) \\ \quad \\ x_n = r\sqrt{1-z_n^2}\cdot \cos(2\pi n\phi) \\ \quad \\ y_n = r\sqrt{1-z_n^2}\cdot \sin(2\pi n\phi) \\ \quad \\ \phi = \frac{\sqrt{5} - 1}{2}$
其中， $r$ 表示球半径， $\phi$ 即黄金分割比例。

代码

通过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个点的效果：
在这里插入图片描述
取1000个点的效果：

RuiH.AI

关注

1
点赞
踩
6

收藏

觉得还不错? 一键收藏
0
评论
Python numpy使用记录5.生成球面斐波那契网格

Python numpy使用记录5.生成球面斐波那契网格前言斐波那契网格公式代码前言本篇学习使用numpy生成球面斐波那契网格。斐波那契网格公式对于一个球体，要在球体表面上均匀选取N个点，就可以使用斐波那契网格公式：zn=r(2n−1N−1)xn=r1−zn2⋅cos⁡(2πnϕ)yn=r1−zn2⋅sin⁡(2πnϕ)ϕ=5−12z_n = r(\frac{2n-1}{N} - 1) \\ \quad \\x_n = r\sqrt{1-z_n^2}\cdot \cos(2\pi n\phi)
复制链接

扫一扫