how to calculate the three sphere.

Hi all,

I have the image 20 images in subplot figure. then the image have three sphere.

Anyone can help me to calculate all the three sphere.

I used this command

T = regionprops('table', outt22,'Area','Centroid')

note: outt22 is the mask of my three volume. but it show only one area of the circle among the images.

T =

  1×2 table

    Area        Centroid    
    ____    ________________

     25     149.08    114.68

[KSSV]

Do you have each indivudal image? Or only the image which is shown?

[mohd akmal masud]

Yes I have individual images.

But if i view individual image, it will be have 20 figure.

So I have to calculate the volume every single image. With repeatition command.

What I supposed is just view 20 images in one subplot, then calculate for all three sphere volume

[KSSV]

You have to calculate that for each image and store into array. Calculating in the attched image is not feasible.

Answer 1

[Image Analyst]

You need to stack these images together into a 3-D volumetric image, and then use regionprops3() instead of regionprops().

[mohd akmal masud]

but my image is 2D sir. how could I stack it into 3D? can guide me?

[Image Analyst]

Like this

filePattern = fullfile(inputFolder, '*.png');
fileList = dir(filePattern)
numFiles = length(fileList);
image3d = zeros(rows, columns, numFiles, 'uint8');
for k = 1 : numFiles
    fullFileName = fullfile(fileList(k).folder, fileList(k).name)
    thisImage = imread(fullFileName)
    image3d(:, :, k) = thisImage;
end
% Make binary image
bw3 = image3d > someThreshold;
% Analyze it
props3 = regionprops3(bw3, 'Volume')

[mohd akmal masud]

ok, let me try first. tq sir. will get you back

[mohd akmal masud]

sir, if my name folder contain images is testimages, then my command like below?

filePattern = fullfile(testimages, '*.png');
fileList = dir(filePattern)
numFiles = length(fileList);
image3d = zeros(rows, columns, numFiles, 'uint8');
for k = 1 : numFiles
    fullFileName = fullfile(fileList(k).folder, fileList(k).name)
    thisImage = imread(fullFileName)
    image3d(:, :, k) = thisImage;
end
% Make binary image
bw3 = image3d > someThreshold;
% Analyze it
props3 = regionprops3(bw3, 'Volume')

[Image Analyst]

inputFolder = 'D:\My Pictures\testimages'; % Wherever it is...
if ~isfolder(inputFolder)
    message = sprintf('Input Folder does not exist:\n\n%s', inputFolder);
    uiwait(warndlg(message));
    return;
end
filePattern = fullfile(inputFolder, '*.png')

[mohd akmal masud]

I have try it. below is my command. but some error.

filePattern = fullfile('C:\Users\Akmal\Desktop\I-131 256 28.02.2020\I-131 SPECT NEMA VALIDATION 01112019 256X256\testimages', '*.png');
fileList = dir(filePattern)
numFiles = length(fileList);
image3d = zeros(rows, columns, numFiles, 'uint8');
for k = 1 : numFiles
    fullFileName = fullfile(fileList(k).folder, fileList(k).name)
    thisImage = imread(fullFileName)
    image3d(:, :, k) = thisImage;
end

ERROR

fileList = 

41×1 struct array with fields:

name
folder
date
bytes
isdir
datenum

Unrecognized function or variable 'rows'.

what is the meaning of Unrecognized function or variable 'rows'.

[mohd akmal masud]

I wrote this, then function

image3d = zeros(256, 256, numFiles, 'uint8');

[Image Analyst]

I thought that you knew in advance how many rows and columns your images have. If not, you can easily find out by just reading in one of them.

yourImage = imread(yourFileNameOfOneImage);
[rows, columns, numberOfColorChannels] = size(yourImage)