about world2vox method

[fzpr]

Hello Dr.mirsky.
I have a problem in world2vox method and this line

CT-GAN/utils/dicom_utils.py

Line 118 in 5601ca1

world_coord = np.dot(np.linalg.inv(np.dot(orientation, np.diag(spacing))), world_coord - origin)

orientation value is [1. 0. 0. 0. 1. 0.] and spacing value is in format [z , y , x] which after np.diag that value is

[[z 0 0]
 [0 y 0]
 [0 0 x]]

so in np.dot(orientation, np.diag(spacing)) we cannot calculate np.dot because orientation and spacing cannot be multiplied and then throw exception with this message:
ValueError: shapes (6,) and (3,3) not aligned: 6 (dim 0) != 3 (dim 0)
Can you help me about this?

[ymirsky]

Here are some resources to help you convert the dicom 6 element orientation array to the 3x3 orientation array. https://nipy.org/nibabel/dicom/dicom_orientation.html http://dicomiseasy.blogspot.com/2013/06/getting-oriented-using-image-plane.html?m=1 You may just be able to use the identity matrix if it is all aligned for all scans.

…

On Sat, Jul 10, 2021, 12:55 fzpr ***@***.***> wrote: Hello Dr.mirsky. I have a problem in world2vox method and this line https://github.com/ymirsky/CT-GAN/blob/5601ca1648d1b2c1ca7bffaac3e6dce803ded961/utils/dicom_utils.py#L118 orientation value is [1. 0. 0. 0. 1. 0.] and spacing value is in format [z , y , x] which after np.diag that value is [[z 0 0] [0 y 0] [0 0 x]] so in np.dot(orientation, np.diag(spacing)) we cannot calculate np.dot because orientation and spacing cannot be multiplied and then throw exception with this message: ValueError: shapes (6,) and (3,3) not aligned: 6 (dim 0) != 3 (dim 0) Can you help me about this? — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#3>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ACYEV2ZJC63AWO25LKFGTSTTXB3RPANCNFSM5AEMSGSA> .