盖特，这就是我的做法。请注意我使用
isocontour
为了one step. Just a simple MATLAB "contour" call may also do the job, but that requires plotting to a figure so I went with an FEX function.
％获得二进制图像
i = phantom（'Modified Shepp-Logan'，nx）> 0;
用零的％垫图像以在图像边界处启用边框
temp = zeros(size(I)+2);
温度（2：end-1,2：end-1）= i;
i = temp;
% Get an isocontour
轮廓线= 0.5;
[Lines,Vertices,Objects] = isocontour(I,contourThreshold);
Vertices = fliplr(Vertices);% Get it back in XY from IJ
％三角剖分Isocontour中的所有PT，并检查哪些三重分输入/输出
DT = delaunayTriangulation(Vertices);
fc = dt.incenter;
in = Interp2（i，fc（：：，1），fc（：：，2））> =轮廓tresthreshold;
% Show the result
图，imagsc（i），保持on，，，，
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(in,:),“ faceColor”，，，，'g'）
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(~in,:),“ faceColor”，，，，'c'）
plot(fc(in,1),fc(in,2),'b。'，，，，fc(~in,1),fc(~in,2),'y.'）
为了i=1:length(Objects)
点=对象{i};
绘图（顶点（点，1），顶点（点，2），，'颜色'，，，，'M'）；
结尾
请注意，您也可以通过
bwperim
而不是Isocontour ...看起来像：
% Get an isocontour
[a，b] = find（bwperim（i））;
Vertices = [b,a];
％三角剖分Isocontour中的所有PT，并检查哪些三重分输入/输出
DT = delaunayTriangulation(Vertices);
fc = dt.incenter;
in = Interp2（i，fc（：：，1），fc（：：，2））== 1;
% Show the result
图，imagsc（i），保持on，，，，
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(in,:),“ faceColor”，，，，'g'）
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(~in,:),“ faceColor”，，，，'c'）
plot(fc(in,1),fc(in,2),'b。'，，，，fc(~in,1),fc(~in,2),'y.'）
And if you were going for
minimal
traingulation, you could try something like this:
% Get a reduced set of boundary vertices
bb = bwboundaries(I);
为了k = 1：长度（bb）
dp = diff（bb {k}，[]，1）;
pdiff = bsxfun(@rdivide, dP, sum(abs(dP),2));
idx = find（任何（pdiff -pdiff -circshift（pdiff，1），2））;
bb {k} = bb {k}（idx，:);
结尾
顶点= fliplr（cat（1，bb {：}））;
％三角剖分Isocontour中的所有PT，并检查哪些三重分输入/输出
DT = delaunayTriangulation(Vertices);
fc = dt.incenter;
in = Interp2（i，fc（：：，1），fc（：：，2））> 0;
图，imagsc（i），保持on，，，，
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(in,:),“ faceColor”，，，，'g'）
修补（'vertices'，，，，DT.Points,'faces'，，，，DT.ConnectivityList(~in,:),“ faceColor”，，，，'c'）
plot(fc(in,1),fc(in,2),'b。'，，，，fc(~in,1),fc(~in,2),'y.'）
