Fitting lines and polynomialon on edge after edge extraction and finding a certain distance
3 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Ron Herman
il 25 Lug 2020
Commentato: Image Analyst
il 26 Lug 2020
I got the following edge after extracting it from a image which has been provided by me in attachmnent. I need to fit a red line along the lower edge and green horizontal line as shown in 2nd figure and find the distance.
What approach should I use???
The desired result I seek??
1 Commento
Matt J
il 25 Lug 2020
I have attached the matlab code for extracting the edge.
It would be more convenient for us if you simply attach a .mat file containing h.
Risposta accettata
Image Analyst
il 25 Lug 2020
Modificato: Image Analyst
il 25 Lug 2020
Really easy. Call bwconvhull(bw, 'union'). Then scan column-by-column using find() to get the last/lowest white point. Then fit those bottom rows to a quadratic with polyfit(), if you want - it's optional because the curve is not really needed to find the height of the last column. Then you simple need to find the last column where there is some white and subtract it from the max value. Something like (untested)
mask = bwconvhull(mask, 'union');
[rows, columns] = size(mask);
bottomRows = zeros(1, columns);
for col = 1 : columns
thisCol = mask(:, col);
r = find(thisCol, 1, 'last');
if ~isempty(r)
bottomRows(col) = r;
end
end
lastColumn = find(bottomRows > 0, 1, 'last');
bottomLine = max(bottomRows);
yline(bottomLine, 'Color', 'g', 'LineWidth', 2)
distance = bottomLine - bottomRows(lastColumn)
% Optional fitting
goodIndexes = bottomRows > 0;
c = 1:columns;
x = c(goodIndexes);
y = bottomRows(goodIndexes);
coefficients = polyfit(x, y, 2);
xFit = c;
yFit = polyval(coefficients, xFit);
hold on;
plot(xFit, yFit, 'r-', 'LineWidth', 2)
yLast = yFit(lastColumn);
plot(lastColumn, yLast, 'r+', 'MarkerSize', 30)
fittedDistance = bottomLine - yLast
% Will be different than other distance because the other one was actual
% and this one is from the regression.
Più risposte (0)
Vedere anche
Prodotti
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!