Rotating DICOM file around z-axis

5 views (last 30 days)
Anouk Baeten
Anouk Baeten on 8 Sep 2022
Commented: Cris LaPierre on 8 Sep 2022
clear variables
az = 3.1415926536;
Rz = [cos(az) -sin(az) 0; sin(az) cos(az) 0; 0 0 1];
Y = dicomread('C:\Users\20182002\OneDrive - TU Eindhoven\Documents\AFSTUDEREN\Matlab\horizontal_flip\EPI02\IMG-0001-00015.dcm');
Y_rotation = Y*(Rz);
dicomwrite(Y_rotation, ("Image_flip.dcm"));
Hello,
I have a DICOM image (see attachment) which i love to rotate about the Z axis. For this I wrote a code (see above) but when I run it, I get the error: 'MTIMES (*) is not fully supported for integer classes. At least one argument must be scalar.'
So I was looking for a solution and found that you can use double() or single() to fix it. For me it would be double(Rz).
But this don't seem to work for me, can somebody help me?
  4 Comments
Cris LaPierre
Cris LaPierre on 8 Sep 2022
As I said in my answer, you could use double, but that leads to a new error about a dimension mismatch error when multiplying your two arrays.

Sign in to comment.

Answers (1)

Cris LaPierre
Cris LaPierre on 8 Sep 2022
Y is uint16 while Rz is double. However, converting Y to double will just lead to a dimension error with your matrix multiplication, since Y is 384x384 while Rz is 3x3.
Perhaps you want to use imrotate or imrotate3?
az = 3.1415926536;
Rz = [cos(az) -sin(az) 0; sin(az) cos(az) 0; 0 0 1];
unzip('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1120005/IMG-0001-00015.dcm.zip');
Y = dicomread('IMG-0001-00015.dcm');
imshow(Y,[])
Y_rotation = imrotate(Y,rad2deg(az));
figure
imshow(Y_rotation,[])
  3 Comments
Cris LaPierre
Cris LaPierre on 8 Sep 2022
Where are you putting the Z axis? Whether it's in the center or a corner, I feel like visually the result is the same.

Sign in to comment.

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!

Translated by