How to convert dicom image world coordinates into Cartesian coordinates

33 次查看（过去 30 天）

Saman Tariq 2019-12-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495272-how-to-convert-dicom-image-world-coordinates-into-cartesian-coordinates

评论： Rik 2019-12-8

I have coordinates of arteries seed points in world coordinate system how can I convert them to the Cartesian ones, by using the DICOM header information of CTA image.

E.g I have 512x512x299volume of dicom images and Ihave a seedpoint represented in world coordinates say (54.3057 56.3043 78.9999), how could I get the corresponding cartesian coordinate of this world coordinate using header information in Matlab??

5 个评论
显示 3更早的评论隐藏 3更早的评论

Saman Tariq 2019-12-8

Rik Thank you for your help but after subtracting from the ImagePositionPatient tag it still not gives integer values of coordinates.

Rik 2019-12-8

Are you certain they should be integers?

Contact the person who created those coordinates.

请先登录，再进行评论。

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

回答（1 个）

Walter Roberson 2019-12-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/495272-how-to-convert-dicom-image-world-coordinates-into-cartesian-coordinates#answer_405375

See https://www.mathworks.com/matlabcentral/answers/215612-how-can-i-convert-pixels-to-mm-for-a-specific-dicom-image#answer_178301 including the comment that links to extending to 3D.

I suspect that when you talk about cartesian coordinates, you are referring to row, column, layer coordinates.

Note that typically for images, what you care about is coordinates relative to corner of the volume, rather than absolute coordinates. In such cases, you can typically translate coordinate systems just by knowing the x, y, and z resolution. However, you specifically asked for world coordinates, so you will have to pay attention to the Image Position (0020,0032) and Image Orientation (0020,0037). I have to say, though, that it is very rare to see an imaging scanner with an attached GPS, as would be necessary to work with world coordinates.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Rik 2019-12-8

CT is generally the easiest to work with. You only need to worry about the physical dimensions of the voxels. Then you can simply use one corner voxel as the origin, and use the voxel dimensions to translate coordinates in voxel offsets.

Note that the slice thickness sounds like the distance between two slice, but that information is actually encoded in the slice interval (which is an optional tag, so you might have to calculate it from the z position of two slices).

Saman Tariq 2019-12-8

@Rik I only wanted to change coordinates my seedpoint (54.3057 56.3043 78.9999) so I don"t know which tag from dicom hearder is used and how it will be done.

请先登录，再进行评论。

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