Calculate the distance between two objects automatically using matlab.
조회 수: 26 (최근 30일)
이전 댓글 표시
hi there!
i am currently doing project on image processing in matlab, in which I'm assigned to calculate the Euclidean distance between two object automatically.here is the image link https://imageshack.com/i/pcUJ9tQgj i first converted the image to binary using 'im2bw', then i applied 'canny' edge detection method and obtained the edges, and by using 'find' function i got the corresponding pixel values. And i am not able to processed further, problem is i have to get array of pixel coordinates for object A and object B separately. and the code has to find which pixel coordinates is nearer between both objects and calculate the nearer distance. I'd very much appreciate your input.
Sincrely, lakshya
댓글 수: 0
채택된 답변
Image Analyst
2014년 8월 25일
You forgot to attach your image, which people usually do when they want image processing advice. So I don't know if edge detection is what you want to do or not. Perhaps it is if you have something like phase contrast or DIC microscopy images, but if you have objects of one intensity on a background of a different intensity then you can probably use thresholding. That method is explained in my Image Segmentation Tutorial http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862 I find the blobs' centroids. From that you can find the distances between the centroids.
You didn't really explain your definition of distance clearly. I think you mean Euclidean /Pythagorean distance, but there's a chance you mean the Hausdorf distance. See this page to understand the Hausdorf distance: http://cgm.cs.mcgill.ca/~godfried/teaching/cg-projects/98/normand/main.html, though chances are you've never heard of it.
In my File Exchange there are also color segmentation methods. How you do your task really depends on what your images look like.
댓글 수: 24
Image Analyst
2020년 11월 13일
No of course not. It's too basic and simple for anyone to have written a paper just on that.
추가 답변 (6개)
Ahmed Mehar
2015년 7월 1일
You Can calculate the Distance of multiple Objects of Center point by using this code:
k=2; for i=1:1:length(g)-1 x(i) = g(i).Centroid(1); y(i) = g(i).Centroid(2);
x(k)=g(k).Centroid(1); y(k)=g(k).Centroid(2);
distance=sqrt((x(i)-x(k))^2+(y(i)-y(k))^2);
댓글 수: 0
Vihan P
2021년 4월 23일
@Image Analyst I have been following the code you attached above 'test3' for finding the distance between two objects. I am unable to find the distance between the objects, instead I am able to find the distance between the boudary of the image to one of the objects. I am attaching the input as well as the output image that I am getting after applying your code, please help me find the minimum distance between these two bones. ![Output Image](https://www.mathworks.com/matlabcentral/answers/uploaded_files/594935/Output%20Image%20.png)
![Output Image](https://www.mathworks.com/matlabcentral/answers/uploaded_files/594935/Output%20Image%20.png)
![Input Image](https://www.mathworks.com/matlabcentral/answers/uploaded_files/594940/Input%20Image.png)
댓글 수: 5
Vihan P
2021년 4월 25일
Hello @Image Analyst Yes, he did help me find a solution. Although, I would love to understand how things can be changed in the code that you have provided.
Image Analyst
2021년 4월 26일
@Vihan P, see code adapted for your image. I also did it both ways, using sqrt(), where I just plotted one of the closest pairs, and pdist2() where I plotted all 3 of the closest pairs.
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;
fprintf('Beginning to run %s.m ...\n', mfilename);
fontSize = 13;
% 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 color demo image.
folder = pwd; %'C:\Users\Lakshya\Documents\Temporary';
baseFileName = 'knee.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
grayImage = imread(fullFileName);
% Get the dimensions of the image. numberOfColorBands should be = 3.
[rows, columns, numberOfColorBands] = size(grayImage);
% Display the original image.
subplot(2, 2, 1);
imshow(grayImage);
axis on;
title('Original Gray Scale Image', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
% Binarize the image
binaryImage = imbinarize(grayImage);
% Display the image.
subplot(2, 2, 2);
imshow(binaryImage);
title('Binary Image', 'FontSize', fontSize);
% Fill the outline to make it solid so we don't get boundaries
% on both the inside of the shape and the outside of the shape.
binaryImage = imfill(binaryImage, 'holes');
% Display the image.
subplot(2, 2, 3);
imshow(binaryImage);
% bwboundaries() returns a cell array, where each cell contains the row/column coordinates for an object in the image.
% Plot the borders of all the coins on the original grayscale image using the coordinates returned by bwboundaries.
hold on;
boundaries = bwboundaries(binaryImage);
numberOfBoundaries = size(boundaries, 1);
for k = 1 : numberOfBoundaries
thisBoundary = boundaries{k};
plot(thisBoundary(:,2), thisBoundary(:,1), 'r', 'LineWidth', 3);
end
title('Filled Binary Image with Boundaries', 'FontSize', fontSize);
hold off;
% Define object boundaries
numberOfBoundaries = size(boundaries, 1)
boundary1 = boundaries{1};
boundary2 = boundaries{2};
boundary1x = boundary1(:, 2);
boundary1y = boundary1(:, 1);
x1=1;
y1=1;
x2=1;
y2=1;
overallMinDistance = inf; % Initialize.
index1 = 1;
index2 = 1;
for k = 1 : length(boundary2)
boundary2x = boundary2(k, 2);
boundary2y = boundary2(k, 1);
% For this blob, compute distances from boundaries to edge.
allDistances = sqrt((boundary1x - boundary2x).^2 + (boundary1y - boundary2y).^2);
% Find closest point, min distance.
[minDistance(k), indexOfMin] = min(allDistances);
if minDistance(k) < overallMinDistance
overallMinDistance = minDistance(k);
x1 = boundary1x(indexOfMin);
y1 = boundary1y(indexOfMin);
x2 = boundary2x;
y2 = boundary2y;
index2 = k;
index1 = indexOfMin;
end
end
% Report to command window.
fprintf('Min Distance from sqrt() method = %f at index %d of boundary 1 and index %d of boundary 2.\n', ...
overallMinDistance, index1, index2);
hFig = figure;
h1 = subplot(1, 2, 1);
imshow(binaryImage);
axis on;
title('Closest Distance from sqrt()', 'FontSize', fontSize);
h2 = subplot(1, 2, 2);
imshow(binaryImage);
axis on;
title('Closest Distances from pdist2()', 'FontSize', fontSize);
hFig.WindowState = 'maximized';
hold on;
% Draw a line between point 1 and 2
line(h1, [x1, x2], [y1, y2], 'Color', 'y', 'LineWidth', 3);
%======================================================================================
% For comparison, use pdist2()
allDistances2 = pdist2(boundary1, boundary2);
minDistance2 = min(allDistances2(:));
% Find all points that have that min distance - there may be several that have it.
[r, c] = find(allDistances2 == minDistance2)
boundary1x = boundary1(:, 2);
boundary1y = boundary1(:, 1);
boundary2x = boundary2(:, 2);
boundary2y = boundary2(:, 1);
for k = 1 : length(r)
% Report to command window.
index1 = r(k);
index2 = c(k);
fprintf('Min Distance from pdist2() method = %f at index %d of boundary 1 and index %d of boundary 2.\n', ...
minDistance2, index1, index2);
xLine = [boundary1x(index1), boundary2x(index2)];
yLine = [boundary1y(index1), boundary2y(index2)];
line(h2, xLine, yLine, 'Color', 'm', 'LineWidth', 1.5);
end
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/597295/image.png)
You can't see the magenta lines in the right image very well, but they're there.
Rhandrey Maestri
2022년 9월 17일
편집: Rhandrey Maestri
2022년 9월 17일
Following this code. Can you tell me how to detect the bubble border(black region) and the tube wall(this black line near the bubble). And then calculate the distance between the border of the bubble and the closest wall pixel near this bubble?
Also if I know that the inside diameter is 15 mm. Can I have this distances in mm?
Appreciate
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1127895/Screenshot_20220917-130734_Adobe%20Acrobat.jpg)
댓글 수: 1
Image Analyst
2022년 9월 17일
@Rhandrey Maestri, start a new question with this since it's sufficiently different. In there, show us how you processed the horizontal profile
horizontalProfile = mean(grayImage);
plot(horizontalProfile, 'b-');
to find the two valleys, and how you used thresholding to get the bubble, then how you use find() to find the edges, row-by-row, of the inner walls and the outer bubble walls. If you still need help, I'll help you there.
There is a variety of ways to do this, and it could be easier if you knew that the walls and bubbles were in approximately the same location every time, like just mask off the middle part and deal with only the central part, then threshold and take the 3 largest blobs with bwareafilt.
Fariya
2022년 12월 20일
댓글 수: 1
Image Analyst
2022년 12월 20일
Use imdistline to find the distance in pixels. Use the camera calibration functionality in the Computer Vision Toolbox if you need real world distances.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!