[ITK Community] OrientImageFilter ASL >> RAI wrong output?

Matt McCormick matt.mccormick at kitware.com
Sun Dec 1 21:10:30 EST 2013


Hi Pascal,

Welcome to the list and good first post!

The OrientImageFilter class is deprecated and not supported.

For more information, see [1] [2]

Thanks,
Matt

[1] http://www.slicer.org/slicerWiki/index.php/Slicer:VTK_ITK_Coordinates
[2] http://www.slicer.org/slicerWiki/index.php/Coordinate_systems


On Thu, Nov 28, 2013 at 10:32 AM, <p.wettmann at hispeed.ch> wrote:

> hi,
>
> I intend to read an arbitrary MRI file (e.g. RAI, ASL, ...) and convert it
> to RAI for further use. If i enter RAI, use orientImageFilter, output is
> correct RAI. However, if i enter ASL, output is not RAI. L/R axis is
> inverse. If I use LAI or a axis flip afterwards, the image seems to be
> correct, despite the fact that it is shifted.
>
> Here a code snipets:
>
>          typedef short PixelType;
>          const unsigned int VolumeDimension = 3;
>          typedef itk::Image< PixelType, VolumeDimension > ImageType;
>          typedef itk::ImageSeriesReader< ImageType >    ReaderType;
>
>          [...]
>
>          ReaderType::Pointer reader = ReaderType::New();
>
>          itk::OrientImageFilter<ImageType,ImageType>::Pointer orienter =
>          itk::OrientImageFilter<ImageType,ImageType>::New();
>
>          FlipImageFilterType::Pointer flipFilter
>             = FlipImageFilterType::New ();
>
>          [...]
>
>          reader->SetFileNames(fileNamesNew);
>          reader->Update();
>
>          rval = reader->GetOutput();
>
>          orienter->UseImageDirectionOn();
>
>  orienter->SetDesiredCoordinateOrientation(itk::SpatialOrientation::ITK_COORDINATE_ORIENTATION_RAI);
>          //orienter->SetDesiredCoordinateOrientationToAxial();
>          orienter->SetInput(rval);
>          orienter->Update();
>
>          rval = orienter->GetOutput();
>
> //         itk::FixedArray<bool, 3> flipAxes;
> //         flipAxes[0] = true;
> //         flipAxes[1] = false;
> //         flipAxes[2] = false;
>
> //         flipFilter->SetInput(rval);
> //         flipFilter->SetFlipAxes(flipAxes);
>
>
> Am I missing something? This is my first e-mail to a mailing list, so
> please let me know any improvements of posting or style. :)
>
> Thanks for any help!
>
> Pascal
>
> _______________________________________________
> Community mailing list
> Community at itk.org
> http://public.kitware.com/cgi-bin/mailman/listinfo/community
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/community/attachments/20131201/e048b905/attachment.html>


More information about the Community mailing list