MATLAB Answers

0

how to convert 2d dicom images to polar maps

Asked by andhavarapu lokesh on 30 Nov 2016
Latest activity Commented on by andhavarapu lokesh on 3 Dec 2016
how to convert the 2d dicom images each containg 72 frames of each how do i convert these each dicom images containg 72 frames to polar maps where each arc has different regions. I need like this image.
projectdir = 'E:\SHIVA BACKUP\THYROID\P1\newcodes\data1\13002';
% y = length(projectdir);
y=72;
X = zeros(128, 128, 1, 72, y, 'uint8');
% Read the series of images.
for p=1:1:y
thisfile = sprintf('IM_%d.dcm', p);
filename = fullfile( projectdir, thisfile );
imdata = dicomread(filename);
imsize = size(imdata);
if ~isequal( imsize, [128 128 1 72] )
fprintf('file is unexpected size %s instead of [128 128 1 72], skipping "%s"\n', mat2str(imsize), filename);
else
X(:, :, :, :, p) = imdata;
end
end
figure();
% Display the first image stack.
z=montage(X(:,:,:,:,71),[]);
axis image
title('First file, all layers')
imcontrast(z)

  3 Comments

sir how to use this code since my dicom images are already gray scale which i need to convert in to polar form and display
where to put this code exactly in my code

Sign in to comment.