Multiframe Image arrays 4-D DICOM

Margarida Moura

2018 6 19

0 个回答

31 次查看（30 天）

显示更早的评论

0 个投票

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
imag_bw_label=bwlabel(imag_bw);
imshow(label2rgb(imag_bw_label));
propriedades= regionprops(imag_bw_label);

6 个评论
显示 4更早的评论隐藏 4更早的评论

Bader 2025-1-18

编辑：Bader 2025-1-18

Thank you, Crise, for your answer. By "fuse," I mean in the context of multi-modal imaging, such as PET-MR.

My DICOM image in multi frame images (4D uint16 DICOM) with 48 frame and the CT in anotther folder in 3D (512x512x130 uint16).

I have tried this the the next code and i got this erorr [[Unable to perform assignment because the size of the left side is 128-by-128 and the size of the right side is 128-by-128-by-1-by-128]]

......................................................................................

imgSize = [info.Rows, info.Columns, 1, numFiles];

dicomImages = zeros(imgSize, 'like', dicomread(fullfile(dicomFolder, dicomFiles(1).name)));

% Read each DICOM file

for i = 1:numFiles

dicomImages(:,:,1,i) = dicomread(fullfile(dicomFolder, dicomFiles(i).name));

end

................................................................................................

Is there any code available to combine my data?

Many Thanks