CT spine saggital image

Hi all,

i have CT spine images and i need to work on both saggital and axial planes , but when i used dicominfo followed by dicomread commands i got the actual ct image which is in axial plane only. is there anyway to extract the images in saggital plane as i need to apply segmentation on it?

waiting any reply

thanks in advance

Answer 1

[Amith Kamath]

Assuming your data has at least 3 dimensions, you only need to index it correctly so as to get a 2D matrix for displaying. For the want of a dataset to demonstrate this, I have used mri. Try running this:

load mri
imshow(D(:,:,1,1)) %transverse plane, first slice.
imshow(imrotate(squeeze(D(:,55,1,:)),90),[]) %sagittal plane, 55th slice.

This is done using SQUEEZE which removes all the singleton dimensions, and IMROTATE which rotates the image appropriately.