Rotating DICOM file around z-axis

3 次查看（过去 30 天）

Anouk Baeten 2022-9-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1800550-rotating-dicom-file-around-z-axis

评论： Cris LaPierre 2022-9-8

IMG-0001-00015.dcm.zip

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 个评论
显示 2更早的评论隐藏 2更早的评论

Stephen23 2022-9-8

"Or how would you do this rotation?"

IMROTATE

Cris LaPierre 2022-9-8

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.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Cris LaPierre 2022-9-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1800550-rotating-dicom-file-around-z-axis#answer_1048685

在 MATLAB Online 中打开

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,[])