MATLAB: How to draw a best fit line based on pixels, on the image itself

digital image processingregression

I have picked certain pixel points on an image and I wanted to draw a best fit line based on the distribution of those points right on the image itself.

Edit:

My input is the white pixels you see on the image. They all deviate to a certain degree horizontally from the center of the image, being zero. So for my data set I have values that are column numbers in reference to the center. This is done for each row.

For example, first row the white pixel may be at +3, that is 3 columns to the right of the center Second row the pixel may be at -27, 27 columns to the left of the center.

This continues all the way down to the image going row by row. I want to now draw a best fit line of these points on the image itself.

Best Answer

Please post your image and annotate it to indicate which pixels the line should be drawn through and indicate if you already know their (x,y) or (row,column) locations. In the meantime, see my attached polyfit demo.

And explain what "the distribution" means to you. Is it the histogram (probability density function), or is it the spatial arrangement/layout of where the pixels are?

Related Solutions

MATLAB: How to display the input image, after when the threshold values exceeds certain limit

It's just 3 lines. Just mask the original image:

binaryImage = grayImage > 150;
maskedImage = zeros(size(grayImage)); % Initialize.

maskedImage(binaryImage) = grayImage(binaryImage);

See the full demo below:

% Demo by ImageAnalyst
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables.
workspace;  % Make sure the workspace panel is showing.
fontSize = 14;
% Read in a standard MATLAB gray scale demo image.
folder = fullfile(matlabroot, '\toolbox\images\imdemos');
baseFileName = 'coins.png';
% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);
% Check if file exists.
if ~exist(fullFileName, 'file')
  % File doesn't exist -- 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 in the search path folders.', 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 Grayscale Image', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
% Give a name to the title bar.
set(gcf,'name','Demo by ImageAnalyst','numbertitle','off') 
% Let's compute and display the histogram.
[pixelCount grayLevels] = imhist(grayImage);
subplot(2, 2, 2); 
bar(pixelCount);
grid on;
title('Histogram of original image', 'FontSize', fontSize);
xlim([0 grayLevels(end)]); % Scale x axis manually.
yl = ylim();
line([150 150], [0 yl(2)], 'Color', 'r');
% Threshold at 150
binaryImage = grayImage > 150;
% Display the original gray scale image.
subplot(2, 2, 3);
imshow(binaryImage, []);
title('Binary Image', 'FontSize', fontSize);
% Mask the image
maskedImage = zeros(size(grayImage)); % Initialize.
% Assign values from grayImage that are > 150.
maskedImage(binaryImage) = grayImage(binaryImage);
% Display the original gray scale image.
subplot(2, 2, 4);
imshow(maskedImage, []);
title('Masked Image', 'FontSize', fontSize);

MATLAB: Crop dark area and display the face only

Convert to grayscale and then threshold and find the bounding box

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. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
% Read in the color image.
folder = pwd;
baseFileName = 'm2.bmp';
% 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
rgbImage = imread(fullFileName);
% Get the dimensions of the image.  numberOfColorBands should be = 3.
[rows, columns, numberOfColorBands] = size(rgbImage);
% Display the original color image.
subplot(2, 3, 1);
imshow(rgbImage, []);
axis on;
title('Original Color Image', 'FontSize', fontSize);
drawnow;
% Enlarge figure to full screen.
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
% Histogram the red channel, but there is a huge spike at 255 so let's not count those.
subplot(2, 3, 2);
grayImage = rgb2gray(rgbImage);
imshow(grayImage, []);
title('Gray Scale version', 'FontSize', fontSize);
subplot(2, 3, 3);
histogram(grayImage(grayImage<255))
grid on;
title('Histogram of Gray Scale Image', 'FontSize', fontSize);
someThresholdValue = 150;
binaryImage = grayImage > someThresholdValue;
% Take largest blob
binaryImage = bwareafilt(binaryImage, 1);
subplot(2, 3, 4);
imshow(binaryImage);
grid on;
title('Binary Image', 'FontSize', fontSize);
% Label the binary image.
labeledImage = bwlabel(binaryImage);
props = regionprops(labeledImage, 'BoundingBox');
bbox = props.BoundingBox
croppedImageRGB = imcrop(rgbImage, bbox);
subplot(2, 3, 5);
imshow(croppedImageRGB);
grid on;
title('Cropped RGB Image', 'FontSize', fontSize);
% Gray scale version of cropped image
croppedImageGray = imcrop(grayImage, bbox);
subplot(2, 3, 6);
imshow(croppedImageGray);
grid on;
title('Cropped Gray Scale Image', 'FontSize', fontSize);

Best Answer

Related Solutions

MATLAB: How to display the input image, after when the threshold values exceeds certain limit

MATLAB: Crop dark area and display the face only

Related Question