Question

matlab中isosurface函数的原理

这几天一直在研究scilab的绘图并和MATLAB做了些比较，有个函数很让我头痛，就是MATLAB里面的isosurface函数不知道它是怎样实现的，求高手指导。
例如在MATLAB中画个单位圆可以用isosurface函数这样实现：
function solidsphere=Unitsphere(res)

if (res>50)
res=50;
elseif (res<1)
res=1;
end

res=1/res;
[x,y,z]=meshgrid(-1-res:res:1+res, ...
-1-res:res:1+res, -1-res:res:1+res);
solidsphere=sqrt( x.^2 + y.^2 + z.^2 ) - 1 ;
obj=isosurface(solidsphere);
obj.facecolor=[0 1 0];obj.edgecolor=[0 0 0];obj.facelighting='phong';
renderpatch(obj);

其中isosurface返回的是登高面上的点和面的信息，因为我想把它转到scilab中去，想知道这个过程是怎么计算的，有谁研究过没有啊，急求。
高手有写过这个函数的能否帮助下？？感激不尽！
有资料的麻烦发到邮箱吧！seven457@Yahoo.cn
我看过了帮助文档，讲的特简单，我要了解的是他的计算原理，不是他的用法简单点就是isosurface是怎么得到顶点和面的信息的，帮助文档里可没有这个吖。

Answer 1

　　函数定义：isosurface 是MATLAB中绘制三维隐函数图像的一个函数工具。
　　函数原理：
　　1、原型：isosurface（x，y，z，v，isovalue）；
　　2、参数：v是关于网格数据x，y，z的体数据，isovalue是对应于v的水平基下的关联数据；
　　3、举例：如果使用isosurface（x，y，z，v，isovalue，colors）就可以设置颜色，然后利用patch命令，对面，颜色值等进行显示。

Answer 2

同求

Answer 3

看isosurface的源代码

Answer 4

在一堆点(点云)中找等值面。
分步实现：
1。找等值点(找出“值==某某”的所有点)
2。把等值点连接成面。
已经完成。
第二步有点小难度，最简单是方式是三角化。
分步三角化（对等值点）：
1。取任意一点
2。找到最近的二点，连接成三角形。
3。用三角形的一条边，去外部找最近的一点，连接成一个新的三角形
4。重复3。直到所有的点都在三角形顶点中
这叫区域生长，或叫游行某某，行进某某的

Answer 5

自己查看帮助文档啊，里面有说明还有例子