MATLAB: I have attached two photos of component. In one of the piece there is a portion cut. Is it possible to produce output with or without comparison with the default using image processing. can someone provide the code
image processingImage Processing Toolboximage segmentation
Related Solutions
Do you mean like this:
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 = 22;%--------------------------------------------------------------------------------------------------------
% READ IN IMAGE
folder = pwd;baseFileName = 'Q4_4.png';% 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(:, :, 2); % Take green channel.
end% Display the image.
subplot(2, 3, 1);imshow(grayImage);axis('on', 'image');title('Original Grayscale Image', 'FontSize', fontSize, 'Interpreter', 'None');impixelinfo;hFig = gcf;hFig.WindowState = 'maximized'; % May not work in earlier versions of MATLAB.
drawnow;%--------------------------------------------------------------------------------------------------------% Binarize the image.
binaryImage = ~imbinarize(grayImage);% Get rid of particles less than 200 in area.
binaryImage = bwareaopen(binaryImage, 200);% Eliminate partial blobs - those touching the edge of the image.
binaryImage = imclearborder(binaryImage);subplot(2, 3, 2);imshow(binaryImage, []);axis('on', 'image');title('Mask Image', 'FontSize', fontSize, 'Interpreter', 'None');impixelinfo;drawnow;% Find blobs
props = regionprops(binaryImage, 'Centroid', 'Area');xy = vertcat(props.Centroid)allAreas = [props.Area];subplot(2, 3, 3);histogram(allAreas, 10);grid on;title('Histogram of Blob Areas', 'FontSize', fontSize, 'Interpreter', 'None');xlabel('Area', 'FontSize', fontSize, 'Interpreter', 'None');ylabel('Count', 'FontSize', fontSize, 'Interpreter', 'None');% Get blobs smaller than 2500
[labeledImage, numBlobs] = bwlabel(binaryImage);smallIndexes = find(allAreas < 2500);largeIndexes = find(allAreas >= 2500);% Extract the small blobs with ismember()
smallBlobs = ismember(labeledImage, smallIndexes);subplot(2, 3, 4);imshow(smallBlobs);title('Small Blobs', 'FontSize', fontSize, 'Interpreter', 'None');% Extract the large blobs with ismember()
largeBlobs = ismember(labeledImage, largeIndexes);subplot(2, 3, 5);imshow(largeBlobs);title('Large Blobs', 'FontSize', fontSize, 'Interpreter', 'None');
See my Image Processing Tutorial in my File Exchange: https://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862&sort=downloads_desc
Then run this code:
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 gray scale demo image.
folder = pwd; % Determine where demo folder is (works with all versions).
baseFileName = 'abc.png';% 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; endendrgbImage = 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(rgbImage)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(rgbImage);
% 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 = rgbImage(:, :, 2); % Take green channel.
else grayImage = rgbImage; % It's already gray scale.
end% Now it's gray scale with range of 0 to 255.
% Display the image.
subplot(2, 2, 1);imshow(grayImage, []);title('Original Image', 'FontSize', fontSize, 'Interpreter', 'None');axis('on', 'image');hp = impixelinfo();%------------------------------------------------------------------------------
% Set up figure properties:
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 1, 0.96]);% 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')drawnow;% Binarize the image
binaryImage = ~imbinarize(uint8(grayImage));% Display the image.subplot(2, 2, 2);imshow(binaryImage, []);title('Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');axis('on', 'image');hp = impixelinfo();drawnow;% Label the image
labeledImage = bwlabel(binaryImage);% Label each blob with 8-connectivity, so we can make measurements of it
[labeledImage, numberOfBlobs] = bwlabel(binaryImage, 8);% Apply a variety of pseudo-colors to the regions.
coloredLabelsImage = label2rgb (labeledImage, 'hsv', 'k', 'shuffle');% Display the pseudo-colored image.
subplot(2, 2, 3);imshow(coloredLabelsImage, []);title('Labeled Image', 'FontSize', fontSize, 'Interpreter', 'None');axis('on', 'image');hp = impixelinfo();% Get all the blob properties. Can only pass in originalImage in version R2008a and later.
props = regionprops(labeledImage, 'all');numberOfBlobs = size(props, 1);% Get all the areas
allAreas = [props.Area]% Get all the Perimeters
allPerimeters = [props.Perimeter]% Get all the Solidities
allSolidities = [props.Solidity]% Show histogram of areas:
subplot(2, 2, 4);histogram(allAreas);grid on;title('Area Distribution', 'FontSize', fontSize, 'Interpreter', 'None');xlabel('Area', 'FontSize', fontSize);ylabel('Count', 'FontSize', fontSize);
Best Answer