clc;
close all;
imtool close all;
clear;
workspace;
format longg;
format compact;
fontSize = 20;
folder = 'C:\Users\mark\Documents\Temporary';
baseFileName = '36543017.png';
fullFileName = fullfile(folder, baseFileName);
if ~exist(fullFileName, 'file')
errorMessage = sprintf('Error: %s does not exist in the %s folder.', baseFileName, folder);
uiwait(warndlg(errorMessage));
fullFileName = baseFileName;
if ~exist(fullFileName, 'file')
errorMessage = sprintf('Error: %s does not exist in the search path folders either.', fullFileName);
uiwait(warndlg(errorMessage));
return;
end
end
grayImage = imread(fullFileName);
[rows columns numberOfColorBands] = size(grayImage);
if numberOfColorBands > 1
grayImage = grayImage(:,:,1);
end
subplot(2, 2, 1);
imshow(grayImage, []);
title('Original Grayscale Image', 'FontSize', fontSize);
set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
set(gcf,'name','Demo by ImageAnalyst','numbertitle','off')
binaryImage = grayImage < 65000;
subplot(2, 2, 2);
imshow(binaryImage, []);
title('Binary Image', 'FontSize', fontSize);
binaryImage = imclose(binaryImage, true(3));
subplot(2, 2, 3);
imshow(binaryImage, []);
title('After morphological closing', 'FontSize', fontSize);
interiors = imclearborder(~binaryImage);
interiors = bwareaopen(interiors, 60);
subplot(2, 2, 4);
imshow(interiors, []);
title('Interiors', 'FontSize', fontSize);
porosity = sum(interiors(:))/numel(interiors);
message = sprintf('The porosity (area fraction) of the interiors is %.4f', porosity);
uiwait(helpdlg(message));
Best Answer