# How to create ROI then know the total pixel number in ROI

15 views (last 30 days)
mohd akmal masud on 26 Feb 2022
Commented: yanqi liu on 1 Mar 2022
Hai Everyone, Anyone know how to create ROI (as red circle) in my image dicom (3D) (as attached), then I can know the total pixel counts in ROI it self
This is my coding to view the images. also the function of imshow3D
clc
clear all
close all
info = dicominfo('I-13125610N1.dcm');
gp=info.SliceThickness;
ps=info.PixelSpacing;
spect=(squeeze(spect));%smooth3
aa=size(spect);aa=aa(3);
figure, imshow3D(spect)
yanqi liu on 1 Mar 2022
may be upload I-13125610N1.dcm to analysis

AndresVar on 27 Feb 2022
Edited: AndresVar on 27 Feb 2022
You can use drawcircle to create a circular roi and then make a mask form it. Then you can sum or find the number of non-zero elements in the mask: Circular region of interest - MATLAB (mathworks.com)
This is an example using matlab's peppers image. When you run it and the image is displayed just click to draw a circle. Modify the circle to show the number of pixels.
clear;
close all;
imshow3D(I);
roi = drawcircle; % interactive
roi.HandleVisibility='off';
function allevents(src,evt)
% when the roi was modified update label with the number of pixel inside
src.Label = sprintf('%g',count);
end
BTW by total pixel counts you mean the number of pixel? or the total intensity counts? if you want intensity counts then you can apply the mask to the image and just sum. sum(im(mask),'all')
AndresVar on 27 Feb 2022
Here is a verision of allevents that gets the intensity counts also
function allevents(src,evt)
% when the roi was modified update label with the number of pixel inside
% find the image data in the imshow3d axis
imdata = findall(src.Parent,'Type','Image');
imdata = imdata(1).CData; % there should be only 1 anyway
src.Label = sprintf('%g (%g)',count,intcount);
end

R2021b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!