如何检测 MATLAB 图中的特定轮廓？

Question

我正在尝试从 MATLAB pcolor 图中提取特定轮廓（从一个名为 M 的矩阵生成，该矩阵也在此处上传），如下图所示： p>

^图。原图

但是我想只提取轮廓，如下所示：

^{图。所需的情节}

黄色线是我想提取的，这里为了说明目的而突出显示。

但我得到的结果如下：

^{图。获得的地块}

除了曲线之外，我不需要任何其他位/杂物。

我在matrix M 上应用以下代码以获得错误结果。

[~, threshold] = edge(M, 'sobel');
fudgeFactor = 0.6;
BWs = edge(M,'sobel', threshold * fudgeFactor);
se90 = strel('line',6,70);
se0 = strel('line',3,90);
BWsdil = imdilate(BWs, [se90 se0]);
BWnobord = imclearborder(BWsdil,18);
seD = strel('cube',3);
BWfinal = imerode(BWnobord,seD);
BWfinal = imerode(BWfinal,seD);
pcolor(BWfinal);colormap(jet);   shading interp;   colorbar;

如何只根据需要获取轮廓？我必须使用任何边缘检测过滤器还是必须在应用任何过滤器之前执行平滑处理？

更新：我应用了高斯滤波器并获得了以下图像。上部的噪音仍然存在。我只想提取曲线。 Image after applying a Gaussian filter

sigma = 30;
smoothImage = imgaussfilt(M,sigma);
smoothGradient = imgradient(smoothImage,'CentralDifference');
pcolor(smoothGradient);colormap(jet);   shading interp;   colorbar;

Answer 1

你可以做这样的事情，虽然这在数学上不是很合理，并且对于你提供的特定矩阵有点硬编码..

clear; close all; clc;
load('Matrix M.mat');

M(181:182,:)=0; % some noise (?) you can skip this line, it will still work..
M(M>0)=1; % binary

% smoothing
sigma = 10;
gauss_submatrix = imgaussfilt(M, sigma); 

% thresholding
threshold = 0.05;
gauss_submatrix(gauss_submatrix<threshold) = 0; 
gauss_submatrix(gauss_submatrix>=threshold) = 1; 

% getting the boundary
a = cumsum(gauss_submatrix);
[rows, cols] = find(a==1);

figure;
hold on;
imagesc(M(1:450,:));
line(cols, rows);
hold off;
ylim([0 450]);
xlim([0 size(M, 2)]);

您可以使用sigma 和threshold 让边界更接近您的实际需求。