细分生成球面三角网格（matlab）

最新推荐文章于 2023-04-15 15:10:04 发布

weixin_30247307

最新推荐文章于 2023-04-15 15:10:04 发布

阅读量1.9k

点赞数 1

文章标签： matlab

原文链接：http://www.cnblogs.com/yezhangxiang/archive/2011/03/22/1991323.html

版权

本文介绍如何利用Matlab生成球面三角网格，通过细分正多面体（尤其是正四面体和正二十面体）来逼近球面。文章提供了细分函数的实现，并展示了细分过程的图像，强调了正二十面体细分后的效果更佳。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

论文需要生成一个球面网格，最简单的想法是通过正多面体的细分来逼近球面。我们知道柏拉图多面体有五个，三个是由正三角形组成（正四，八，二十面体），一个由正方形组成（正六面体），一个由正五边形组成（正十二面体）。三角形的细分比较简单，而且保持拓扑可以连续细分。

由于只是为了得到图片，就直接用matlab来实现了，方便，matlab中画着三角网格可以用trimesh或者patch

Syntax

trimesh(Tri,X,Y,Z) trimesh(Tri,X,Y,Z,C) trimesh(...'PropertyName',PropertyValue...)

patch(X,Y,C) patch(X,Y,Z,C) patch(FV) patch(...'PropertyName',propertyvalue...) patch('PropertyName',propertyvalue,...)

但是都要用一个参数来指定每个三角面是由那些点组成（比如trimesh中的Tri）。我们用vertex储存点，face指定面由哪些点构成，细分的函数如下：

 1 function [vertex, face] = Subdivision(vertex, face)
 2 face_num = size(face, 1);
 3 vertex_num = size(vertex, 1);
 4 new_vertexs = zeros(face_num*3, 3);
 5 new_faces = zeros(face_num*3, 3);
 6 for i = 1 : face_num
 7     new_vertexs(3*(i-1)+1, :)=(vertex(face(i, 1), :)+vertex(face(i, 2), :))/2;
 8     new_vertexs(3*(i-1)+1, :)= new_vertexs(3*(i-1)+1, :)/norm(new_vertexs(3*(i-1)+1, :),2);
 9     new_vertexs(3*(i-1)+