我正在尝试制作一个遵循“像素”边缘的轮廓pcolor在 Matlab 中绘图。这可能在图片中得到最好的解释。这是我的数据图。黄色数据（data==1）和蓝色数据（data==0）之间有明显的界限：
请注意，这是一个pcolor绘图，因此每个“正方形”本质上都是一个像素。我想返回一个轮廓faces黄色数据的pixels, not只是黄色数据的边缘。
因此，输出轮廓（绿线）穿过像素的面部（红点）的中点。
请注意，我不希望轮廓遵循数据的中心点（黑点），这会产生类似绿线的效果。这可以很容易地实现contour.
另外，如果有帮助的话，我有一些可能有用的网格。我在像素中间有点（显然，因为这就是我在这里绘制的），在角上也有点，并且在西/东面和北/南面上有点。如果你熟悉荒川网格，这是一个 Arakawa-C 网格，所以我有 rho-、u-、v- 和 psi- 点。
我尝试过插值、交织网格和其他一些方法，但我没有任何运气。任何帮助都会非常感激，并且会阻止我发疯。
干杯，戴夫
EDIT:
Sorry, I simplified the images to make what I was trying to explain more obvious, but here is a larger (zoomed out) image of the region I'm trying to separate:
正如您所看到的，它是一个复杂的轮廓，先向“西南”方向前进，然后环绕并移回“东北”方向。这是我想穿过黑点绘制的红线：
您可以通过一些修改来解决这个问题我发布的解决方案 to a 相关问题。我在问题中使用了示例图像蒙版的一部分data。首先，您需要填充面罩上的孔，您可以使用imfill从图像处理工具箱:
x = 1:15; % X coordinates for pixels
y = 1:17; % Y coordinates for pixels
mask = imfill(data, 'holes');
接下来，应用我的其他答案中的方法来计算一组有序的轮廓坐标（位于像素角上）：
% Create raw triangulation data:
[cx, cy] = meshgrid(x, y);
xTri = bsxfun(@plus, [0; 1; 1; 0], cx(mask).');
yTri = bsxfun(@plus, [0; 0; 1; 1], cy(mask).');
V = [xTri(:) yTri(:)];
F = reshape(bsxfun(@plus, [1; 2; 3; 1; 3; 4], 0:4:(4*nnz(mask)-4)), 3, []).';
% Trim triangulation data:
[V, ~, Vindex] = unique(V, 'rows');
V = V-0.5;
F = Vindex(F);
% Create triangulation and find free edge coordinates:
TR = triangulation(F, V);
freeEdges = freeBoundary(TR).';
xOutline = V(freeEdges(1, [1:end 1]), 1); % Ordered edge x coordinates
yOutline = V(freeEdges(1, [1:end 1]), 2); % Ordered edge y coordinates
最后，您可以在像素边缘的中心获得所需的坐标，如下所示：
ex = xOutline(1:(end-1))+diff(xOutline)./2;
ey = yOutline(1:(end-1))+diff(yOutline)./2;
这是显示结果的图：
imagesc(x, y, data);
axis equal
set(gca, 'XLim', [0.5 0.5+size(mask, 2)], 'YLim', [0.5 0.5+size(mask, 1)]);
hold on;
plot(ex([1:end 1]), ey([1:end 1]), 'r', 'LineWidth', 2);
plot(ex, ey, 'k.', 'LineWidth', 2);
