how to convert 2d dicom images to polar maps

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)

[Walter Roberson]

Merged with later duplicate http://www.mathworks.com/matlabcentral/answers/315204-how-to-convert-the-frames-in-dicom-image-into-polar-maps

[andhavarapu lokesh]

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]

where to put this code exactly in my code

Answer 1

[Walter Roberson]

https://www.mathworks.com/matlabcentral/fileexchange/17933-polar-to-from-rectangular-transform-of-images