delaunay

近期论坛上有不少讨论delaunay函数的帖子。似乎主要有以下问题：
1、delaunay函数各参数的意义
2、知道几何边界时，用delaunay函数划分三角形网格由于区域 内部没有点，质量很差，怎么改进
3、怎样避免产生过于狭长的delaunay 三角形
4、 凹多边形的情况怎么处理
第1个问题，看看帮助应该能解决。第2个问题， delaunay本来是用来对离散点进行三角剖分，内部没有点时并不合适。除非特别处理。第3个问题，估计是利用 delaunay和meshgrid来划网格，边界附近会产生 狭长的delaunay 三角形 ，这个也可以做特别处理。第4个问题，可以用在划分好网格后删掉域外的三角形即可。

由于我也经常使用 delaunay来处理背景积分问题，因此仔细琢磨了一下用 delaunay来划分已知边界的几何区域的可行方案，在此和大家分享一下，也是抛砖引玉，希望大家有更好的方法。

方案一： 先对区域 delaunay剖分，删掉域外的三角形，然后将剩下的三角形的边细分，得到新的离散点，然后再次 delaunay 剖分，然后再次细分边，这样循环下去，直到达到一定的尺寸为止

方案二： 利用 delaunay和meshgrid函数。将边界细分得到相比原区域边界更加密集边界点，用 meshgrid得到包含整个区域的点，将域内的点和边界点一起 delaunay 剖分。

讨论：
方案一对于一开始就有很小边界段的情况情况较差，容易出现狭长单元（比如边界有圆弧的话属于这种情况）。还有就是前一步的边界轮廓很清楚，看着别扭。方案二中间的网格能搞保证形状较好。对于边界附近的内部点，容易导致边界单元畸变，可以将离边界太近的点进行删除，这样得到的形状比较好

综合来说，方案二较好，尤其是当删掉离边界太近的内部点。贴出程序，望大家多多指点，共同进步。

P.S. 当然，matlab自身也有很好的网格划分函数，在pdetool中有用到，不过关于几何描述那块比较难以理解（我不是很理解）。另外matlab语言写的划分网格的程序很多，网上可以找到不少很优秀的。这里仅限于简单的使用delaunay来划分。

也就是约束Delaunay三角形
X = [0 0; 16 0; 16 2; 2 2; 2 3; 8 3; 8 5; 0 5];
C = [1 2; 2 3; 3 4; 4 5; 5 6; 6 7; 7 8; 8 1];
dt = DelaunayTri(X, C);
subplot(3,1,1);
io=dt.inOutStatus();
triplot(dt(io,:),X(:,1),X(:,2));
axis([-1 17 -1 6]);
xlabel('Constrained Delaunay triangulation(inside)', 'fontweight','b');
hold on;
plot(X(C'),X(C'+size(X,1)),'-r', 'LineWidth', 2);
hold off;
subplot(3,1,2)
triplot(dt(~io,:),X(:,1),X(:,2));
axis([-1 17 -1 6]);
xlabel('Constrained Delaunay triangulation(outside)', 'fontweight','b');
hold on;
plot(X(C'),X(C'+size(X,1)),'-r', 'LineWidth', 2);
hold off;
dt.Constraints = [];
subplot(3,1,3);
triplot(dt);
axis([-1 17 -1 6]);
xlabel('Unconstrained Delaunay triangulation', 'fontweight','b');