is their any way to calculate automatically this parameter from the image ?

Hi all,

I am looking for the calculate the the parameter mentioned in the image, Is that possible to calculate from the image ?

(I have attached the original image and other one attached to understand you that which value i need to calculate, (Black width, highlighted with the lime colour))Thank you.

[Anish Khan]

@Image Analyst, can you please? Thank you in advance.

[Anish Khan]

Can anyone help ? Thanks

[Jenny Lee]

following this question

Answer 1

[Image Analyst]

Of course it's possible. Try this. Adapt as needed:

% Demo by Image Analyst, February, 2021.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clearvars;
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 16;
fprintf('Beginning to run %s.m ...\n', mfilename);

%-----------------------------------------------------------------------------------------------------------------------------------
% Read in image.
folder = [];
baseFileName = 'original image.JPG';
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;
	end
end
grayImage = imread(fullFileName);
[rows, columns, numberOfColorChannels] = size(grayImage)
if numberOfColorChannels == 3
	grayImage = rgb2gray(grayImage);
end
% Display the test image full size.
subplot(2, 3, 1);
imshow(grayImage, []);
axis('on', 'image');
caption = sprintf('Original Image : "%s"', baseFileName);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Set up figure properties:
% Enlarge figure to full screen.
hFig1 = gcf;
hFig1.Units = 'Normalized';
hFig1.WindowState = 'maximized';
% 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.
hFig1.Name = 'Demo by Image Analyst';

%--------------------------------------------------------------------------------------------------------
% Get histogram
subplot(2, 3, 2);
imhist(grayImage);
grid on;
xlabel('Gray Level', 'FontSize', fontSize);
ylabel('Count', 'FontSize', fontSize);

% Threshold image.
mask = grayImage < 75;
% Clean up the mask to get only the stripes we want.
mask(:, 1 : 762) = false;
mask(:, 2615 : end) = false;
% Get rid any blob not 1500 wide
labeledImage = bwlabel(mask);
props = regionprops(labeledImage, 'BoundingBox');
allBB = vertcat(props.BoundingBox);
widths = allBB(:, 3);
wideBlobs = find(widths > 1500);
% Extract only those that are wide enough.
mask = ismember(labeledImage, wideBlobs);

% Display mask image.
subplot(2, 3, 3);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Label each blob with 8-connectivity, so we can make measurements of it
[labeledImage, numberOfBlobs] = bwlabel(mask, 8);
% Apply a variety of pseudo-colors to the regions.
coloredLabelsImage = label2rgb (labeledImage, 'hsv', 'k', 'shuffle');
% Display the pseudo-colored image.
subplot(2, 3, 4);
imshow(coloredLabelsImage);
title('Labeled Blobs', 'FontSize', fontSize, 'Interpreter', 'None');
% Get all the blob properties.  Can only pass in originalImage in version R2008a and later.
blobMeasurements = regionprops(labeledImage, grayImage, 'all');
numberOfBlobs = size(blobMeasurements, 1);

% For each blob, for each line, get the left edge and right edge then widths.
leftEdges = zeros(1, columns);
rightEdges = zeros(1, columns);
widths = zeros(columns, numberOfBlobs);
for blob = 1 : numberOfBlobs
	thisBlob = ismember(labeledImage, blob);
	% Display mask image.
	subplot(2, 3, 5);
	imshow(thisBlob, []);
	axis('on', 'image');
	caption = sprintf('Mask of Blob #%d of %d', blob, numberOfBlobs);
	title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
	drawnow;
	hp = impixelinfo(); % Set up status line to see values when you mouse over the image.
	for col = 1 : columns
		thisCol = thisBlob(:, col);
		t = find(thisCol, 1, 'first');
		if ~isempty(t)
			leftEdges(col) = t;
			rightEdges(col) = find(thisCol, 1, 'last');
		end
	end
	widths(:, blob) = (rightEdges - leftEdges)';
	subplot(2, 3, 6);
	plot(widths, '-', 'LineWidth', 1);
	hold on;
	grid on;
	drawnow;
	legendStrings{blob} = sprintf('Blob #%d', blob);
end
grid on;
title('Widths', 'FontSize', fontSize, 'Interpreter', 'None');
xlabel('Row', 'FontSize', fontSize);
ylabel('Width in Pixels', 'FontSize', fontSize);
legend(legendStrings);

% Display mask image again.
subplot(2, 3, 5);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

fprintf('Done running %s.m\n', mfilename);

[Anish Khan]

Thanks @Image Analyst

[Image Analyst]

So did it solve it?

[Anish Khan]

Yes it solved my problem, thanks a lot @Image Analyst. Your code explained everything very well. It require a little further modification but that purpose i don't want to disturb great expert like you.

Thanks

[Image Analyst]

I deleted the attachment. Original posters can "Accept" and "vote" for Answers to award "reputation points" to the answerer. Others can "Vote" for answers to award reputation points (but can't Accept the answer if they're not the original poster). Thanks @Anish Khan for Accepting the Answer. 🙂

[Image Analyst]

I created a generic image of blurry stripes of slightly varying widths. No one will be able to tell what it is, especially with the aliasing going on when it's subsampled for display. It's just a bunch of stripes.

% Demo by Image Analyst, March, 2021.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clearvars;
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 16;
fprintf('Beginning to run %s.m ...\n', mfilename);

%-----------------------------------------------------------------------------------------------------------------------------------
% Read in image.
folder = [];
baseFileName = 'original image3.JPG';
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;
	end
end
grayImage = imread(fullFileName);
[rows, columns, numberOfColorChannels] = size(grayImage)
if numberOfColorChannels == 3
	grayImage = rgb2gray(grayImage);
end
% Display the test image full size.
subplot(2, 3, 1);
imshow(grayImage, []);
axis('on', 'image');
caption = sprintf('Original Image : "%s"', baseFileName);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Set up figure properties:
% Enlarge figure to full screen.
hFig1 = gcf;
hFig1.Units = 'Normalized';
hFig1.WindowState = 'maximized';
% 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.
hFig1.Name = 'Demo by Image Analyst';

%--------------------------------------------------------------------------------------------------------
% Get histogram
subplot(2, 3, 2);
imhist(grayImage);
grid on;
xlabel('Gray Level', 'FontSize', fontSize);
ylabel('Count', 'FontSize', fontSize);
title('Histogram of Gray Scale Image', 'FontSize', fontSize, 'Interpreter', 'None');

% Threshold image.
mask = grayImage < 75;
% Clean up the mask to get only the stripes we want.
mask(:, 1 : 762) = false;
mask(:, 2615 : end) = false;
% Get rid any blob not 1500 wide
labeledImage = bwlabel(mask);
props = regionprops(labeledImage, 'BoundingBox');
allBB = vertcat(props.BoundingBox);
widths = allBB(:, 3);
wideBlobs = find(widths > 1500);
% Extract only those that are wide enough.
mask = ismember(labeledImage, wideBlobs);

% Display mask image.
subplot(2, 3, 3);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Label each blob with 8-connectivity, so we can make measurements of it
[labeledImage, numberOfBlobs] = bwlabel(mask, 8);
% Apply a variety of pseudo-colors to the regions.
coloredLabelsImage = label2rgb (labeledImage, 'hsv', 'k', 'shuffle');
% Display the pseudo-colored image.
subplot(2, 3, 4);
imshow(coloredLabelsImage);
caption = sprintf('%d Labeled Blobs', numberOfBlobs);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
% Get all the blob properties.  Can only pass in originalImage in version R2008a and later.
blobMeasurements = regionprops(labeledImage, grayImage, 'all');
numberOfBlobs = size(blobMeasurements, 1);

% For each blob, for each line, get the left edge and right edge then widths.
leftEdges = zeros(1, columns);
rightEdges = zeros(1, columns);
widths = zeros(columns, numberOfBlobs);
for blob = 1 : numberOfBlobs
	thisBlob = ismember(labeledImage, blob);
	% Display mask image.
	subplot(2, 3, 5);
	imshow(thisBlob, []);
	axis('on', 'image');
	caption = sprintf('Mask of Blob #%d of %d', blob, numberOfBlobs);
	title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
	drawnow;
	hp = impixelinfo(); % Set up status line to see values when you mouse over the image.
	for col = 1 : columns
		thisCol = thisBlob(:, col);
		t = find(thisCol, 1, 'first');
		if ~isempty(t)
			leftEdges(col) = t;
			rightEdges(col) = find(thisCol, 1, 'last');
		end
	end
	widths(:, blob) = (rightEdges - leftEdges)';
	subplot(2, 3, 6);
	plot(widths, '-', 'LineWidth', 1);
	hold on;
	grid on;
	drawnow;
	legendStrings{blob} = sprintf('Blob #%d', blob);
end
grid on;
title('Widths', 'FontSize', fontSize, 'Interpreter', 'None');
xlabel('Row', 'FontSize', fontSize);
ylabel('Width in Pixels', 'FontSize', fontSize);
legend(legendStrings, 'Location', 'south');

% Display mask image again.
subplot(2, 3, 5);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

fprintf('Done running %s.m\n', mfilename);

[Image Analyst]

@Mehbuba Shaikh, I made a script, test5.m, to completely synthesize an image of horizontal stripes that no one on earth could possibly say that it has anything to do with your original images because I made it myself out of nothing but just the cos() and linspace() of MATLAB. It's just a completely generic arbitrary image of wavy stripes. No more changes, okay?

% Demo by Image Analyst, February, 2021.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clearvars;
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 16;
fprintf('Beginning to run %s.m ...\n', mfilename);

%-----------------------------------------------------------------------------------------------------------------------------------
% Read in image.
folder = [];
baseFileName = 'original image3.JPG';
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;
	end
end
grayImage = imread(fullFileName);
[rows, columns, numberOfColorChannels] = size(grayImage)
if numberOfColorChannels == 3
	grayImage = rgb2gray(grayImage);
end
% Display the test image full size.
subplot(2, 3, 1);
imshow(grayImage, []);
axis('on', 'image');
caption = sprintf('Original Image : "%s"', baseFileName);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Set up figure properties:
% Enlarge figure to full screen.
hFig1 = gcf;
hFig1.Units = 'Normalized';
hFig1.WindowState = 'maximized';
% 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.
hFig1.Name = 'Demo by Image Analyst';

%--------------------------------------------------------------------------------------------------------
% Get histogram
subplot(2, 3, 2);
imhist(grayImage);
grid on;
xlabel('Gray Level', 'FontSize', fontSize);
ylabel('Count', 'FontSize', fontSize);
title('Histogram of Gray Scale Image', 'FontSize', fontSize, 'Interpreter', 'None');

% Threshold image.
mask = grayImage < 75;
% Clean up the mask to get only the stripes we want.
mask(:, 1 : 762) = false;
mask(:, 2615 : end) = false;
% Get rid any blob not 1500 wide
labeledImage = bwlabel(mask);
props = regionprops(labeledImage, 'BoundingBox');
allBB = vertcat(props.BoundingBox);
widths = allBB(:, 3);
wideBlobs = find(widths > 1500);
% Extract only those that are wide enough.
mask = ismember(labeledImage, wideBlobs);

% Display mask image.
subplot(2, 3, 3);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

% Label each blob with 8-connectivity, so we can make measurements of it
[labeledImage, numberOfBlobs] = bwlabel(mask, 8);
% Apply a variety of pseudo-colors to the regions.
coloredLabelsImage = label2rgb (labeledImage, 'hsv', 'k', 'shuffle');
% Display the pseudo-colored image.
subplot(2, 3, 4);
imshow(coloredLabelsImage);
caption = sprintf('%d Labeled Blobs', numberOfBlobs);
title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
% Get all the blob properties.  Can only pass in originalImage in version R2008a and later.
blobMeasurements = regionprops(labeledImage, grayImage, 'all');
numberOfBlobs = size(blobMeasurements, 1);

% For each blob, for each line, get the left edge and right edge then widths.
leftEdges = zeros(1, columns);
rightEdges = zeros(1, columns);
widths = zeros(columns, numberOfBlobs);
for blob = 1 : numberOfBlobs
	thisBlob = ismember(labeledImage, blob);
	% Display mask image.
	subplot(2, 3, 5);
	imshow(thisBlob, []);
	axis('on', 'image');
	caption = sprintf('Mask of Blob #%d of %d', blob, numberOfBlobs);
	title(caption, 'FontSize', fontSize, 'Interpreter', 'None');
	drawnow;
	hp = impixelinfo(); % Set up status line to see values when you mouse over the image.
	for col = 1 : columns
		thisCol = thisBlob(:, col);
		t = find(thisCol, 1, 'first');
		if ~isempty(t)
			leftEdges(col) = t;
			rightEdges(col) = find(thisCol, 1, 'last');
		end
	end
	widths(:, blob) = (rightEdges - leftEdges)';
	subplot(2, 3, 6);
	plot(widths, '-', 'LineWidth', 1);
	hold on;
	grid on;
	drawnow;
	legendStrings{blob} = sprintf('Blob #%d', blob);
end
grid on;
title('Widths', 'FontSize', fontSize, 'Interpreter', 'None');
xlabel('Row', 'FontSize', fontSize);
ylabel('Width in Pixels', 'FontSize', fontSize);
legend(legendStrings, 'Location', 'south');

% Display mask image again.
subplot(2, 3, 5);
imshow(mask, []);
axis('on', 'image');
title('Mask', 'FontSize', fontSize, 'Interpreter', 'None');
drawnow;
hold on;
hp = impixelinfo(); % Set up status line to see values when you mouse over the image.

fprintf('Done running %s.m\n', mfilename);

[Mehbuba Shaikh]

Thanks a lot @Image Analyst

Mehbuba

[Mehbuba Shaikh]

Thanks a lot @Image Analyst

Mehbuba

[Anish Khan]

Thanks @Image Analyst

Anish

Answer 2

[Hannah Paul]

I tried with segmentation but not able to overcome with the problem. Seems like its not possible.