У меня есть отверстие для размыва почвы, которое развивается со временем, то есть глубина размытия со временем становится глубже. Похоже, что отверстие для размытия имеет параболическую форму, как показано здесь:
Запрошенная поверхность показана синей линией на рисунке ниже:
Поскольку отверстие для размытия симметрично относительно вертикальной оси, я обрезал изображение, оставив только левую половину. Я использовал следующий код:
frm1 = imread('ScourHoleImage.jpg');
imshow(frm1);
CrpWndo = zeros(1, 4);
CrpWndo = [1017, 1127.6, 1606.2, 905.1];
frm1Crp = imcrop(frm1, CrpWndo);
BW1 = im2bw(frm1Crp);
boundaries = bwboundaries(BW1);
boundaries = cell2mat(boundaries);
imshow(BW1);
hold on
plot(boundaries(:,2), boundaries(:,1), '.g');
Приведенный выше код дал мне следующий результат:
Мои вопросы: как удалить ненужные пиксели, отмеченные выбросами? А как проследить поверхность почвы только по полилинии?