Hello!
I need to work with collections of images related by time or view, such as magnetic resonance imaging (MRI) slices or movie frames. (attached)

I need to get the ROI of the image but in different times.
So, I have a multi frame image (4D uint16 DICOM) with 36 frames and I also have my code but it only shows the sum of the 36 frames…

P.S. It is an image of a 3D SPECT scan I need to check the values of ROI in different times to get the intensities/ countings, so that I can get the curve Activity vs Time.

Could please somebody help me?

%%OPEN MULTI FRAME
image_data =dicomread(‘FLOW001_DS.dcm’);
sum2d = sum(image_data,4);
sum2d=256-sum2d;
imshow(sum2d,[])

%Otsu to the generated image
A=otsu(sum2d,2); % testar com outros valores (3 por exemplo);

%Cany countour to select the ROI
BW1 = edge(A,’Canny’);
imshow(BW1)

%Sum original image+ cany ;
sumf=imadd(sum2d,im2double(BW1).*250);
imshow(sumf,[])
imagesc(sumf)

%binary image
umbral=graythresh(sumf);
imag_bw=im2bw(sumf,umbral);
imshow(imag_bw)

%object identification in image

if true
% code
end
imag_bw_label=bwlabel(imag_bw);
imshow(label2rgb(imag_bw_label));
propriedades= regionprops(imag_bw_label);Hello!
I need to work with collections of images related by time or view, such as magnetic resonance imaging (MRI) slices or movie frames. (attached)

I need to get the ROI of the image but in different times.
So, I have a multi frame image (4D uint16 DICOM) with 36 frames and I also have my code but it only shows the sum of the 36 frames…

P.S. It is an image of a 3D SPECT scan I need to check the values of ROI in different times to get the intensities/ countings, so that I can get the curve Activity vs Time.

Could please somebody help me?

%%OPEN MULTI FRAME
image_data =dicomread(‘FLOW001_DS.dcm’);
sum2d = sum(image_data,4);
sum2d=256-sum2d;
imshow(sum2d,[])

%Otsu to the generated image
A=otsu(sum2d,2); % testar com outros valores (3 por exemplo);

%Cany countour to select the ROI
BW1 = edge(A,’Canny’);
imshow(BW1)

%Sum original image+ cany ;
sumf=imadd(sum2d,im2double(BW1).*250);
imshow(sumf,[])
imagesc(sumf)

%binary image
umbral=graythresh(sumf);
imag_bw=im2bw(sumf,umbral);
imshow(imag_bw)

%object identification in image

if true
% code
end
imag_bw_label=bwlabel(imag_bw);
imshow(label2rgb(imag_bw_label));
propriedades= regionprops(imag_bw_label); Hello!
I need to work with collections of images related by time or view, such as magnetic resonance imaging (MRI) slices or movie frames. (attached)

I need to get the ROI of the image but in different times.
So, I have a multi frame image (4D uint16 DICOM) with 36 frames and I also have my code but it only shows the sum of the 36 frames…

P.S. It is an image of a 3D SPECT scan I need to check the values of ROI in different times to get the intensities/ countings, so that I can get the curve Activity vs Time.

Could please somebody help me?

%%OPEN MULTI FRAME
image_data =dicomread(‘FLOW001_DS.dcm’);
sum2d = sum(image_data,4);
sum2d=256-sum2d;
imshow(sum2d,[])

%Otsu to the generated image
A=otsu(sum2d,2); % testar com outros valores (3 por exemplo);

%Cany countour to select the ROI
BW1 = edge(A,’Canny’);
imshow(BW1)

%Sum original image+ cany ;
sumf=imadd(sum2d,im2double(BW1).*250);
imshow(sumf,[])
imagesc(sumf)

%binary image
umbral=graythresh(sumf);
imag_bw=im2bw(sumf,umbral);
imshow(imag_bw)

%object identification in image

if true
% code
end
imag_bw_label=bwlabel(imag_bw);
imshow(label2rgb(imag_bw_label));
propriedades= regionprops(imag_bw_label); multi frame, dynamic, time MATLAB Answers — New Questions

​