[ITK-users] SimpleITK - Resampling Coronal Image in the Axial View

[Sara Gh]

Hello experts,

I have two MRI scans: Axial and Coronal
# *imageA* : Axial MRI scan, size : 256x256x56 (256mmx256mmx168mm) ,
 spacing: 1,1,3  ( slice gap : 1  =>  slice thickness : 2)
# *imageC* : Coronal MRI scan, size : 256x256x68 (256mmx256mmx204mm) ,
 spacing: 1,1,3  ( slice gap : 1  =>  slice thickness : 2)

I am trying to view Coronal image from the axial view i.e. resampling
Coronal image using the direction of Axial image. Here is my code:

resample = sitk.ResampleImageFilter()
resample.SetReferenceImage(imageC)
resample.SetOutputDirection(imageA.GetDirection())
resample.SetOutputOrigin(imageA.GetOrigin())
imageCA = resample.Execute(imageC)
sitk.Show(imageCA)

# Result: imageCA: size : 256x256x68 (256mmx256mmx204mm), spacing: 1,1,3

I believe I have missed something because, in my idea, the size of
resampled coronal image (imageCA) should be 68x256x256 (204mmx256mmx256mm) with
the spacing [3,1,1].

I also tried to set the OutputSpacing to [3,1,1], but the size of the
imageCA changes to 256x256x68 (768mmx256mmx68mm).

I would be thankful if someone could help me in this regard.

Best regards,
Sara Gh
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/insight-users/attachments/20160626/bb543e87/attachment.html>