draw horizontal lines on edges of a image

Question

Tommy Schumacher 2021-1-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/724573-draw-horizontal-lines-on-edges-of-a-image

评论： Image Analyst 2021-1-29

采纳的回答： Image Analyst

Hello,

i have a picture ("problem.png) and i want to draw horizontal lines on the edges in the image, like "result.png"

I have tried it with edge detection but unfortunately no useful results. Is there another way to do this?

I am grateful for any help. THX

2 个评论
显示无隐藏无

KALYAN ACHARJYA 2021-1-23

Please try with other segmentation approaches, may take time, but it is absolutely possible to get those bars band to draw lines.

J. Alex Lee 2021-1-23

I assume one of the greatest challenges with an edge detection is all those vertical polishing lines...maybe you can go a bit beyond edge detection to determine edge angle. In principle, this is taking the arctangent of the gradient vector components, but you'll still be bit by noise...so you can do some more analysis of the eigenvalues of the image gradient (see e.g., structure tensor), and you can get something like below. Combined with some other ways to filter based on other properties of the image, maybe you can achieve your goal

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2021-1-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/724573-draw-horizontal-lines-on-edges-of-a-image#answer_604612

在 MATLAB Online 中打开

Try this:

clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
%--------------------------------------------------------------------------------------------------------
%    READ IN IMAGE
folder = pwd;
baseFileName = 'problem.png';
grayImage = imread(baseFileName);
% Get the dimensions of the image.
% numberOfColorChannels should be = 1 for a gray scale image, and 3 for an RGB color image.
[rows, columns, numberOfColorChannels] = size(grayImage)
if numberOfColorChannels > 1
	% It's not really gray scale like we expected - it's color.
	% Use weighted sum of ALL channels to create a gray scale image.
	grayImage = min(grayImage, [], 3);
end
%--------------------------------------------------------------------------------------------------------
% Display the image.
subplot(2, 3, 1);
imshow(grayImage, []);
axis('on', 'image');
title('Gray Scale Image', 'FontSize', fontSize, 'Interpreter', 'None');
impixelinfo;
hFig = gcf;
hFig.WindowState = 'maximized'; % May not work in earlier versions of MATLAB.
drawnow;
subplot(2, 3, 2);
imhist(grayImage);
grid on;
% Threshold the image.
grayImage = adapthisteq(grayImage);
subplot(2, 3, 3);
imshow(grayImage);
axis('on', 'image');
title('Flattened Image', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
verticalProfile = mean(grayImage, 2);
% Display the image.
subplot(2, 3, 4);
plot(verticalProfile);
grid on;
threshold = 110;
yline(threshold, 'Color', 'r', 'LineWidth', 2);
title('Vertical Profile', 'FontSize', fontSize);
% Threshold the image
binaryImage = imfill(grayImage < threshold, 'holes');
% Take the 3 largest blobs
binaryImage = bwareafilt(binaryImage, 3);
% Snip off small tendrils using imopen()
binaryImage = imopen(binaryImage, true(1, 3));
% Take the 3 largest blobs
binaryImage = bwareafilt(binaryImage, 3);
subplot(2, 3, 5);
imshow(binaryImage);
axis('on', 'image');
title('Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
% If we assume the bands are perfectly horizontal, we cna use strfind() to find the tops and bottoms.
binaryProfile = (verticalProfile > threshold)';
bandStarts = strfind(binaryProfile, [0, 1])
bandStops = strfind(binaryProfile, [1, 0])
for k = 1 : length(bandStarts)
	yline(bandStarts(k), 'Color', 'r', 'LineWidth', 2);
	yline(bandStops(k), 'Color', 'r', 'LineWidth', 2);
end
% If they are tilted, we can use fitPolynomialRANSAC() in the Computer Vision Toolbox.
% Code for this is not shown but is straightforward.

There are lines going across the binary image in red. I assumed they were perfectly horizontal. If they are tilted, you can find the tops and bottoms of each band and use polyfit or fitPolynomialRANSAC() to get the equation of the tilted lines. I didn't give code for that.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Image Analyst 2021-1-24

在 MATLAB Online 中打开

radon_rotate.m

If you don't have the Computer Vision Toolbox you can use the radon transform to rotate an image. See attached demo.

Otherwise scan across the image column by column getting the starting and stopping points of the 3 bands. I'm sure you've already done this by now Tommy since it's so easy, but anyway (for others) here is how we did it:

for col = 1 : length(columns)
    thisColumn = binaryImage(:, col)
    bandStarts = strfind(thisColumn, [0, 1])
    bandStops = strfind(thisColumn, [1, 0])
    topLine1(col) = bandStarts(1);
    topLine2(col) = bandStarts(2);
    topLine3(col) = bandStarts(3);
    bottomLine1(col) = bandStarts(1);
    bottomLine2(col) = bandStops(2);
    bottomLine3(col) = bandStops(3);
end
x = 1 : columns;
topCoeffs1 = polyfit(x, topLine1, 2); % Fit boundary to a quadratic.
topCoeffs2 = polyfit(x, topLine2, 2); % Fit boundary to a quadratic.
topCoeffs3 = polyfit(x, topLine3, 2); % Fit boundary to a quadratic.
bottomCoeffs1 = polyfit(x, bottomLine1, 2); % Fit boundary to a quadratic.
bottomCoeffs2 = polyfit(x, bottomLine2, 2); % Fit boundary to a quadratic.
bottomCoeffs3 = polyfit(x, bottomLine3, 2); % Fit boundary to a quadratic.
% Get smoothed/fitted values.
topY1 = polyfit(topCoeffs1, x);
topY2 = polyfit(topCoeffs2, x);
topY3 = polyfit(topCoeffs3, x);
bottomY1 = polyfit(bottomCoeffs1, x);
bottomY2 = polyfit(bottomCoeffs2, x);
bottomY3 = polyfit(bottomCoeffs3, x);
hold on;
plot(x, topY1, 'r-', 'LineWidth', 2);
plot(x, topY2, 'r-', 'LineWidth', 2);
plot(x, topY3, 'r-', 'LineWidth', 2);
plot(x, bottomY1, 'r-', 'LineWidth', 2);
plot(x, bottomY2, 'r-', 'LineWidth', 2);
plot(x, bottomY3, 'r-', 'LineWidth', 2);

Of course in the loop you need to check to make sure that you actually have 3 bands, not 0 or 1 or 2.

Image Analyst 2021-1-26

Tommy - what's the status on this? Did my Answer solve it? If so, can you "Accept this answer". If my Answer did not help, what's the remaining problem?

Image Analyst 2021-1-29

bandStarts and bandStops are the coordinates of the actual lines. Call polyfit if you want them to be smoother.

请先登录，再进行评论。

draw horizontal lines on edges of a image

2 个评论
显示无隐藏无

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

draw horizontal lines on edges of a image

2 个评论 显示 无隐藏 无

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

2 个评论
显示无隐藏无

4 个评论
显示 2更早的评论隐藏 2更早的评论