Creating approximation line between objects.
显示 更早的评论
Hello, I have image with a few objects I measured centroids of each one, now I want to create approximation line using centroids. Do you have any ideas how could I do that? Thank you in advance.
采纳的回答
Image Analyst
2016-11-19
0 个投票
Compute the centroids with regionprops(). Then pass those into polyfit(). Then use polyval() to get coordinates along the line.
4 个评论
Marcin Heba
2016-11-27
Thank you for yor answer. I wrote a code but results are not so good. Could you look at this and give me some advice?
Here's my code:
clear all
close all
obraz_org = imread('image1.png');
obraz_bw = im2bw(obraz_org);
obraz_label = bwlabel(obraz_bw, 8);
bb = regionprops(obraz_label,'Centroid');
centroids = cat(1, bb.Centroid);
imshow(obraz_bw);
axis on;
hold on
[Y,I]=sort(centroids(:,2));
B=centroids(I,:);
x = B(:,1)
y = B(:,2)
plot(x, y, 'r+')
[p,~,mu] = polyfit(x, y, 1);
f = polyval(p,x,[],mu);
plot(x,f)
hold off
Final image:
Orginal image:
John D'Errico
2016-11-27
编辑：John D'Errico
2016-11-27
Note that polyfit is NOT the right tool for this fit. You need a code that can fit a line with errors in BOTH variables. If you choose to use the wrong tool (i.e., polyfit) then you WILL get a biased line. The slope will be wrong. Always.
My comment above should be a hint as to why polyfit gave the wrong slope. Again, polyfit is the WRONG tool here.
Look on the file exchange. There are lots of such tools. I would start with the keyword "orthogonal regression".
Be careful though, because not all orthogonal regressions were done at all well. (I just looked, and I found a fair amount of crap.) You need to use a tool that has SVD at its core. This one should suffice:
Image Analyst
2016-11-27
Marcin: you got confused when you did the fit as to which coordinate was the independent variable and which was the dependent variable. Fixed code is below:
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 a gray scale demo image.
folder = pwd;
baseFileName = 'image1.png';
% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);
grayImage = imread(fullFileName);
% Get the dimensions of the image.
% numberOfColorBands 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.
% Convert it to gray scale by taking only the green channel,
% which in a typical snapshot will be the least noisy channel.
% grayImage = grayImage(:, :, 2); % Take green channel.
end
% Display the image.
% subplot(2, 2, 1);
imshow(grayImage, []);
axis on;
axis image;
title('Original Grayscale Image', 'FontSize', fontSize, 'Interpreter', 'None');
% Set up figure properties:
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
% Get rid of tool bar and pulldown menus that are along top of figure.
set(gcf, 'Toolbar', 'none', 'Menu', 'none');
% Give a name to the title bar.
set(gcf, 'Name', 'Demo by ImageAnalyst', 'NumberTitle', 'Off')
binaryImage = im2bw(grayImage);
labeledImage = bwlabel(binaryImage, 8);
props = regionprops(labeledImage, 'Centroid');
centroids = [props.Centroid];
xCenters = centroids(1:2:end);
yCenters = centroids(2:2:end);
hold on;
plot(xCenters, yCenters, 'r+', 'LineWidth', 2, 'MarkerSize', 20);
% Fit to a polynomial
coefficients = polyfit(yCenters, xCenters, 1);
% Get fit over all y (row) values.
rowFit = 1 : rows;
% Get an x value (columns value) for each row.
colFit = polyval(coefficients, rowFit);
% Plot it over the image.
plot(colFit, rowFit, 'm-', 'LineWidth', 2);
hold off
Marcin Heba
2016-11-28
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Explore and Edit Images with Image Viewer App 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
