[ITK] [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/community/attachments/20160626/bb543e87/attachment.html>
-------------- next part --------------
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/insight-users