I want to know the code. i'm not good at coding sir.Thankyou.
MATLAB: How to draw a circle in this image from the edges that I marked and calculate the angle from the circle like this
anklefeetimage processingImage Processing Toolbox
Related Solutions
I see several possible places that could be a trapezoid shape, not only the ones you outlined in red. What I would try to do it this, assuming it's fairly vertical
- Find the orientation and centroid of the whole shape with regionprops()
- Divide the shape in two and find the orientation of only the lower part below the centroid.
- Rotate the image by that orientation so that now the bottom part is pretty vertical.
- Scan down from top to bottom using find() to get the left and right boundary of each line.
- Find the place where the width goes from decreasing to fairly stable. This is the bottom of the trapezoid.
- Define some width that will define how wide the "top" of the trapezoid should be.
- Go up from the bottom until the width is the top width.
- Now you have all 4 points of the trapezoid.
Give it a try yourself, because we can't because you didn't attach your two initial binary images. Even if you did, it's good to try it yourself first. Post your code if you can't figure it out, and I'll try to help.
Use regionprops():
% Code to detect angle of a blob.
% By Image Analyst
clc; % Clear the command window.
fprintf('Beginning to run %s.m ...\n', mfilename);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 = 15;%--------------------------------------------------------------------------------------------------------
% READ IN IMAGE
folder = pwd;baseFileName = 'ZoomIn_100f.jpg';% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);% Check if file exists.
if ~exist(fullFileName, 'file') % The file doesn't exist -- didn't find it there in that folder.
% Check the entire search path (other folders) for the file by stripping off the folder.
fullFileNameOnSearchPath = baseFileName; % No path this time.
if ~exist(fullFileNameOnSearchPath, '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; endendgrayImage = imread(fullFileName);% Get the dimensions of the image.
% numberOfColorChannels 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.
% Use weighted sum of ALL channels to create a gray scale image.
% grayImage = rgb2gray(grayImage);
% ALTERNATE METHOD: Convert it to gray scale by taking only the green channel,
% which in a typical snapshot will be the least noisy channel.
grayImage = grayImage(:, :, 1); % Take red channel.
end% Crop it to the right lines.
% grayImage = grayImage(1750:2250, 800:end);
subplot(2, 2, 1);imshow(grayImage, []);impixelinfo;title('Red Channel Image', 'FontSize', fontSize, 'Interpreter', 'None');hFig = gcf;hFig.WindowState = 'maximized'; % May not work in earlier versions of MATLAB.
drawnow;subplot(2, 2, 2);imhist(grayImage);grid on;title('Histogram of Original Image', 'FontSize', fontSize, 'Interpreter', 'None');%--------------------------------------------------------------------------------------------------------% SEGMENTATION OF IMAGE
% Get a binary image by interactively thresholding using the function at
% https://www.mathworks.com/matlabcentral/fileexchange/29372-thresholding-an-image
% 91 seems good for this image.
lowThreshold = 91;highThreshold = 255;% [lowThreshold, highThreshold] = threshold(lowThreshold, highThreshold, grayImage)
mask = grayImage > lowThreshold & grayImage < highThreshold;mask = imclearborder(mask);mask = imfill(mask, 'holes');% Take largest blobs.
mask = bwareafilt(mask, 1);subplot(2, 2, 3);imshow(mask, []);impixelinfo;caption = sprintf('Mask of White SemiCircles');title(caption, 'FontSize', fontSize, 'Interpreter', 'None');% Measure centroid and orientation.
props = regionprops(mask, 'Centroid', 'Orientation')xCenter = props.Centroid(1);yCenter = props.Centroid(2);hold on;plot(xCenter, yCenter, 'r+', 'MarkerSize', 40, 'LineWidth', 2);% Display the masked image.
subplot(2, 2, 4);overlayImage = imoverlay(grayImage, mask, 'r');imshow(overlayImage, []);impixelinfo;caption = sprintf('Masked Image. Angle = %.2f degrees.', props.Orientation);title(caption, 'FontSize', fontSize, 'Interpreter', 'None');hFig = gcf;hFig.WindowState = 'maximized'; % May not work in earlier versions of MATLAB.drawnow;% Make a line at that angle at the centroid.
slope = tan(props.Orientation)x = 1 : columns;y = slope * (x - xCenter) + yCenter;hold on;plot(x, y, 'r-', 'LineWidth', 2);fprintf('Done running %s.m ...\n', mfilename);
Best Answer