how to convert 2d dicom images to polar maps

Question

andhavarapu lokesh 2016-11-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/314747-how-to-convert-2d-dicom-images-to-polar-maps

评论： andhavarapu lokesh 2016-12-3

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)