I am attaching the image . Please help me how to remove the outermost semicicrlces.Please anyone can say the method , how to remove it?
MATLAB: How to remove Outermost semicircle on the grayscale image
image processing
Related Solutions
Justin - are you still there? Perhaps you're still struggling over this so I made a complete and fancy demo for you. Here it is:
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 = 25;%===============================================================================
% Get the name of the image the user wants to use.
baseFileName = 'concrete_inverted.png';% Get the full filename, with path prepended.
folder = []; % Determine where demo folder is (works with all versions).
fullFileName = fullfile(folder, baseFileName);%===============================================================================% Read in a demo image.
grayImage = 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, 2, 1);imshow(grayImage, []);axis on;axis image;caption = sprintf('Original Gray Scale Image');title(caption, 'FontSize', fontSize, 'Interpreter', 'None');drawnow;hp = impixelinfo();% 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')drawnow;% Binarize the image by thresholding.
mask = grayImage > 128;% For some reason, the top two lines of his image are all white. Blacken those 2 lines:
mask(1:2, :) = false;% Display the mask image.
subplot(2, 2, 2);imshow(mask);axis on;axis image; % Make sure image is not artificially stretched because of screen's aspect ratio.
title('Binary Image Mask', 'fontSize', fontSize);drawnow;% Get rid of small blobs.
mask = bwareaopen(mask, 500);% Label the image.
labeledImage = bwlabel(mask);% Find the areas and perimeters
props = regionprops(labeledImage, 'Area', 'Perimeter');allAreas = [props.Area];sortedAreas = sort(allAreas, 'descend')allPerimeters = [props.Perimeter];% Compute circularities
circularities = allPerimeters .^ 2 ./ (4 * pi * allAreas)sortedC = sort(circularities, 'descend')% Keep conly blobs that are nowhere close to circular or compact.
minAllowableCircularity = 10;keeperIndexes = find(circularities >= minAllowableCircularity);mask = ismember(labeledImage, keeperIndexes);% Display the mask image.subplot(2, 2, 3);imshow(mask);axis on;axis image; % Make sure image is not artificially stretched because of screen's aspect ratio.title('Intermediate Cleaned Mask', 'FontSize', fontSize);drawnow;% Get rid of black islands (holes) in struts without filling large black areas.
subplot(2, 2, 4);mask = ~bwareaopen(~mask, 1000);imshow(mask);axis on;axis image; % Make sure image is not artificially stretched because of screen's aspect ratio.title('Final Cleaned Mask', 'FontSize', fontSize);drawnow;
I think the results look pretty good, like what I think you're expecting. Don't you? If not, let me know what I could do to improve it. Basically the same original crack lines are in the image (no change in exterior shape), except that the small black islands in them are filled in, all all the little isolated white island are removed. Is that what you want? There are a few parameters/variables in there, with obvious names, that you can adjust if you want to. Please get back to me.
No need to complicate things with MSER when a simple cropping will do. Simply find the bright panel, and go in a known percentage to crop away the stuff you don't want (everything except the 7 digit number). Trivial. Here's a full demo.
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 = 15;%--------------------------------------------------------------------------------------------------------
% READ IN IMAGE
folder = pwd;baseFileName = '64.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; endend% Read in the image from disk. If storedColorMap is not empty, it's an indexed image with a stored colormap.
[grayImage, storedColorMap] = imread(fullFileName);if ~isempty(storedColorMap) grayImage = ind2rgb(grayImage, storedColorMap);end% 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.
hFig = figure;subplot(2, 3, 1);imshow(grayImage, []);title('Original Grayscale Image', 'FontSize', fontSize, 'Interpreter', 'None');hFig.WindowState = 'maximized'; % May not work in earlier versions of MATLAB.
drawnow;% Display the histogram
subplot(2, 3, 2);imhist(grayImage);grid on;%--------------------------------------------------------------------------------------------------------% SEGMENTATION
% Binarize the image
binaryImage = imbinarize(grayImage);% Display the image.subplot(2, 3, 3);imshow(binaryImage, []);title('Initial Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');axis('on', 'image');% Fill holes.
binaryImage = imfill(binaryImage, 'holes');% Extract the largest blob only.
binaryImage = bwareafilt(binaryImage, 1);% Display the image.subplot(2, 3, 4);imshow(binaryImage, []);title('Filled Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');axis('on', 'image');hp = impixelinfo();% Get the bounding box
props = regionprops(binaryImage, 'BoundingBox');% Crop the image
croppedImage = imcrop(grayImage, props.BoundingBox);% Plot them over the inverted gray scale image.
subplot(2, 3, 5);imshow(croppedImage, []);impixelinfocaption = sprintf('Cropped Image');title(caption, 'FontSize', fontSize, 'Interpreter', 'None');%--------------------------------------------------------------------------------------------------------% EXTRACTION OF ROI
% Now crop away the stuff that is not the 7 digit number we want.
pctH1 = 18;pctH2 = 86;pctV1 = 10;pctV2 = 62;% 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(croppedImage)row1 = round(pctV1 * rows / 100)row2 = round(pctV2 * rows / 100)col1 = round(pctH1 * columns / 100)col2 = round(pctH2 * columns / 100)finalImage = croppedImage(row1:row2, col1:col2);% Plot them over the inverted gray scale image.subplot(2, 3, 6);imshow(finalImage, []);impixelinfoaxis('on', 'image');caption = sprintf('Final Image');title(caption, 'FontSize', fontSize, 'Interpreter', 'None');hold on;
If the panel may be rotated, you can scan the image to find the top of the panel in the binary image then use ransac, fitPolynomialRANSAC(), to find the top edge and its angle, then use imrotate() to rotate it level.
Best Answer