Tiger pugmark binary image analysis

Question

Suvidha 2013-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/69627-tiger-pugmark-binary-image-analysis

采纳的回答： Image Analyst

how could i find the centroid and area of each object in the image automatically.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Image Analyst 2013-4-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/69627-tiger-pugmark-binary-image-analysis#answer_80829

First run imclose() on it to fill in nooks and crannies. Then do imfill() to fill in holes. Then regionprops(). See im image segmentation tutorial for details if you need them: http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862

5 个评论
显示 3更早的评论隐藏 3更早的评论

Image Analyst 2013-4-4

编辑：Image Analyst 2013-4-4

在 MATLAB Online 中打开

Suvidha: OK here. This code pretty much does everything you want as far as I can tell, except the angles which I imagine you can get because you know how to use atand(). So just change the filename and foldername of your image and then run it. Let me know what you think.

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.
format longg;
format compact;
fontSize = 20;
% Read in a color demo image.
folder = 'C:\Users\Suvidha\Documents';
baseFileName = 'pawprint.png';
% 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
originalImage = imread(fullFileName);
% Get the dimensions of the image.  numberOfColorBands should be = 3.
[rows, columns, numberOfColorBands] = size(originalImage);
% Display the original color image.
subplot(2, 2, 1);
imshow(originalImage);
axis on;
title('Original Color Image', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
% Convert to grayscale.
if numberOfColorBands > 1
  originalImage = originalImage(:,:,2);
end
% Binarize the image
binaryImage = originalImage > 128;
% Get rid of stuff near border.
binaryImage = imclearborder(binaryImage);
% Get rid of small things
binaryImage = bwareaopen(binaryImage, 300);
% Do a morphological closing to smooth it some
se = strel('disk', 3, 0);
binaryImage = imclose(binaryImage, se);
% Fill holes
% binaryImage = imfill(binaryImage, 'holes');
% Display image.
subplot(2, 2, 2);
imshow(binaryImage);
axis on;
title('Binary Image', 'FontSize', fontSize);
% Get the convex hull of each blob.
binaryImage = bwconvhull(binaryImage, 'objects');
% Display image.
subplot(2, 2, 3);
imshow(binaryImage);
axis on;
title('Convex Hull Image with Centroids', 'FontSize', fontSize);
% Measure the image
measurements = regionprops(binaryImage, 'Area', 'Centroid');
allAreas = [measurements.Area] % Get all areas into one vector.
numberOfBlobs = length(allAreas)
% Plot centroids on to image
hold on;
for k = 1 : numberOfBlobs
  thisCentroid = measurements(k).Centroid
  plot(thisCentroid(1), thisCentroid(2), 'b+', 'MarkerSize', 10);
end
% Display image in lower right.
subplot(2, 2, 4);
imshow(originalImage);
axis on;
title('Original Image', 'FontSize', fontSize);
% Find the distances between every pair of blobs.
distance = zeros(numberOfBlobs, numberOfBlobs);
for b1 = 1 : numberOfBlobs
  thisCentroid1 = measurements(b1).Centroid
  x1 = thisCentroid1(1);
  y1 = thisCentroid1(2);
  for b2 = (b1 + 1) : numberOfBlobs
    thisCentroid2 = measurements(b2).Centroid
    x2 = thisCentroid2(1);
    y2 = thisCentroid2(2);
    distance(b1,b2) = sqrt((x1-x2)^2+(y1-y2)^2);
    distance(b2,b1) = distance(b1,b2);
    line([x1 x2], [y1 y2], 'Color', 'r', 'LineWidth', 3);
    message = sprintf('This distance is %.3f', distance(b1,b2));
    uiwait(helpdlg(message));
  end
end

Suvidha 2013-4-5

Thank you ... the code has done almost what i want .. i think i can do further from here.

Image Analyst 2013-4-5

Then please mark the answer as "Accepted".

请先登录，再进行评论。

Answer 2

Walter Roberson 2013-4-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/69627-tiger-pugmark-binary-image-analysis#answer_80820

regionprops() and ask for Centroid, and either Area or FilledArea.

Question: when you say "each object in the image", should each pugmark be considered one object, or should each toe be considered separate objects?

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2013-4-3

How do you want to define the distances between the toes? Closest point? Distance between centroids?

Suvidha 2013-4-3

distance between centroids

请先登录，再进行评论。

Tiger pugmark binary image analysis

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

5 个评论
显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

Tiger pugmark binary image analysis

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

5 个评论 显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

3 个评论 显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论