how to convert 2d dicom images to polar maps

1 view (last 30 days)
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
andhavarapu lokesh
andhavarapu lokesh on 3 Dec 2016
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
andhavarapu lokesh
andhavarapu lokesh on 3 Dec 2016
where to put this code exactly in my code

Sign in to comment.

Answers (1)

Walter Roberson
Walter Roberson on 3 Dec 2016

Categories

Find more on DICOM Format in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!