MATLAB Answers

Image Edges: Calculate complexity of a line segment?

22 views (last 30 days)
Philip
Philip on 6 Apr 2012
Can anyone suggest a good way to calculate the 'complexity' of a line segment (obtained from 'bwlabel' after 'canny' edge detection).
I want to define 'complexity' in 2 parts:
  1. How straight the line is
  2. How long the line is
and would like to weight each line segment such that a long and straight line is more favourable than a short and wavy line.
It is ok for the line to change direction so long as they do not change direction regularly (for example, a large square is fine because it is made up of 4 'simple' straight lines). However, edges corresponding to leaves on a tree for example are not good because they change direction too much.
As mentioned, I have used 'bwlabel' after 'canny' edge detection. I then want to iterate through each bwlabel number independently and rank them in order of how (predominately) straight and long they are.
Many thanks for your help!

  0 Comments

Sign in to comment.

Answers (3)

Ashish Uthama
Ashish Uthama on 6 Apr 2012
(I'll try to flesh this soon) How about:
  • Obtain the PixelList property of each line using regionprops
  • Compute the difference in location between pixels (gives you an idea of how much the line wiggles)
  • Maybe take the mean of this difference
  • And weight it with the number of pixels in this segment

  1 Comment

Walter Roberson
Walter Roberson on 6 Apr 2012
In a line that curves, you can get pixel lists where the list does not follow the "outside" (or "inside") of the curve. For example,
.x..
xxx.
x..x
then the ordering of the pixels in the section that touches itself is not necessarily going to "thread" the line.
It might make more sense to use http://www.mathworks.com/help/toolbox/images/ref/bwtraceboundary.html as at least then you know the order of the tracing.

Sign in to comment.


Image Analyst
Image Analyst on 6 Apr 2012
Usually the tortuosity ( http://en.wikipedia.org/wiki/Tortuosity) is what is used. I think it's sometimes called "curl" when dealing with fibers. It is the ratio of the length to the distance between the endpoints. A number of alternative definitions are also discussed on that web page.
You can also use fractal dimension ( http://en.wikipedia.org/wiki/Fractal_dimension).

  3 Comments

Philip
Philip on 10 Apr 2012
Thanks! I have used the simplest form of Tortuosity (length of line / distance between endpoints of the line), and this has helped to some extent... but when ordering the tortuosity of each line segment in ascending/descending order, it does not return results that are quite what I was hoping for.
I have attached an edge image (http://imageshack.us/photo/my-images/27/edgese.png/) with 6 line segments. I would hope that the line segment corresponding to the near vertical line (towards the right of the image) should have the highest tortuosity as it is closest to a straight line. Then I would expect to see the 3 line segments corresponding to the large wavy lines to appear next, followed by the 2 small horizontal blocks at the edge of the image... I think I will need to weight the tortuosity in favour of long lines that do not deviate much, but I am not sure how to...
Here is how I calculate tortuosity:
img = imread('edges.png');
canny_edges = edge(img,'canny');
[LS, qty] = bwlabel(canny_edges,8);
for i=1:qty
[xcoords ycoords] = find(LS==i);
route = bwtraceboundary(LS,[xcoords(1),ycoords(1)],'W',8,length(xcoords),'clockwise');
% Tortuoisity: ratio of length(route) / distance between start/endpoints
% distance = sqrt((x2-x1)^2 + (y2-y1)^2)
len = length(route);
dist = sqrt((route(end,2)-route(1,2))^2 + (route(end,1)-route(1,1))^2);
tort = len/dist;
profile(i,1) = i;
profile(i,2) = tort;
end
profile_sorted = sortrows(profile,-2);
Image Analyst
Image Analyst on 10 Apr 2012
Not how I'd do it. How did you get x1 and x2? I'd call bwmorph to get the endpoints, then calculate the distance for each labeled object. Then get the histogram - use as many bins as there are objects. This will get you the sum of all the lengths of all the lines in just one line of code: perimeters = hist(LS, qty). Then divide them. Sort them if desired.
Philip
Philip on 11 Apr 2012
x1 and x2 are obtained from the output of 'bwtraceboundary.' x1 = route(1,1); x2 = route(end,1); y1 = route(1,2); y2 = route(end,2);
Sorry, can I ask you to elaborate on what you mean when you say you would use 'bwmorph' to get the endpoints? And also, what do you mean by "then divide them"? Do you mean divide the lengths of each line with the corresponding distance between endpoints? or something else?

Sign in to comment.


Image Analyst
Image Analyst on 11 Apr 2012
Philip: Alright, try this code:
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
imtool close all; % Close all imtool figures.
clear; % Erase all existing variables.
workspace; % Make sure the workspace panel is showing.
fontSize = 20;
% Change the current folder to the folder of this m-file.
if(~isdeployed)
cd(fileparts(which(mfilename)));
end
% Check that user has the Image Processing Toolbox installed.
hasIPT = license('test', 'image_toolbox');
if ~hasIPT
% User does not have the toolbox installed.
message = sprintf('Sorry, but you do not seem to have the Image Processing Toolbox.\nDo you want to try to continue anyway?');
reply = questdlg(message, 'Toolbox missing', 'Yes', 'No', 'Yes');
if strcmpi(reply, 'No')
% User said No, so exit.
return;
end
end
% Read in a standard MATLAB gray scale demo image.
folder = 'C:\Users\Philip\Documents\Temporary';
baseFileName = 'edgese.png';
fullFileName = fullfile(folder, baseFileName);
% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);
if ~exist(fullFileName, 'file')
% Didn't find it there. Check the search path for it.
fullFileName = baseFileName; % No path this time.
if ~exist(fullFileName, 'file')
% Still didn't find it. Alert user.
errorMessage = sprintf('Error: %s does not exist.', fullFileName);
uiwait(warndlg(errorMessage));
return;
end
end
grayImage = imread(fullFileName);
% Get the dimensions of the image. numberOfColorBands should be = 1.
[rows columns numberOfColorBands] = size(grayImage);
% Display the original gray scale image.
subplot(2, 2, 1);
imshow(grayImage, []);
title('Original Image', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
set(gcf,'name','Demo by ImageAnalyst','numbertitle','off')
binaryImage = grayImage > 0;
[labeledImage numberOfLines] = bwlabel(binaryImage);
fprintf('numberOfLines = %d\n', numberOfLines);
endpoints = bwmorph(binaryImage, 'endpoints');
subplot(2, 2, 2);
imshow(endpoints, []);
title('Endpoints Image', 'FontSize', fontSize);
% Number of Endpoints should = 2 * numberOfLines because each line will have two endpoints.
numberOfEndpoints = sum(endpoints(:));
fprintf('numberOfEndpoints = %d\n', numberOfEndpoints);
% Find the length of each line
lineLengths = histc(labeledImage(:), 1:numberOfLines)
% Find the endpoints for each
% Then find the distance and the tortuosity.
for lineNumber = 1 : numberOfLines
theseEndpoints = (labeledImage .* endpoints) == lineNumber;
subplot(2, 2, 3);
imshow(theseEndpoints, []);
caption = sprintf('Endpoints for line #%d', lineNumber);
title(caption, 'FontSize', fontSize);
drawnow;
[rows columns] = find(theseEndpoints);
distanceBetweenEndpoints = sqrt((rows(1)-rows(2))^2 + (columns(1)-columns(2))^2);
tortuosity(lineNumber) = lineLengths(lineNumber) / distanceBetweenEndpoints;
fprintf('For line #%d, curve length = %.2f, distance = %.2f,tortuosity = %.3f\n',...
lineNumber, lineLengths(lineNumber), distanceBetweenEndpoints, tortuosity(lineNumber));
end

  9 Comments

Show 6 older comments
Image Analyst
Image Analyst on 12 Apr 2012
Philip: the perimeter computed by regionprops() seems to go center to center and so will work with slanted lines. Here, try this code
binaryImage = logical([...
0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 1 1 1 1 1 1 1 1 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 1 1 0 0 0 0 0 0
0 0 0 0 0 1 0 0 1 0 0 0 0 0
0 0 0 0 1 0 0 0 0 1 0 0 0 0
0 0 0 1 0 0 0 0 0 0 1 0 0 0
0 0 1 0 0 0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 ])
[labeledImage numberOfRegions] = bwlabel(binaryImage)
fprintf('numberOfRegions = %d\n', numberOfRegions);
measurements = regionprops(binaryImage, 'Perimeter', 'area');
allPerimeters = [measurements.Perimeter]/2
allAreas = [measurements.Area]
I trust you can apply this code to my earlier code to get the tortuosity with this perimeter array.
Image Analyst
Image Analyst on 12 Apr 2012
By the way, make sure your blobs have been skeletonized first (use bwmorph) so that they are only a single pixel wide.
Philip
Philip on 12 Apr 2012
Thanks, I have dropped that code into your previous demo, and it seems to be very fast and reliable. I'm trying to see if I can use the 'perimeter' data to help with measuring straight vs. diagonal lines (as per the example above). It seems that the perimeter data returns odd results for more complex lines, which may not be useful for some images that are processed.

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by