[Insight-users] OrientedImageAdaptor or ExtractOrientedImageFilter?

[Frank Miller]

Constantine,

What I do is to represent "slices" with an

  itk::OrientedImage< TPixel, 3 >

with size[2] == 1. You can resample the volume with a "slice" just as
you describe in your first way and then there is no need to use the
ExtractImageFilter.

Note, it is not possible to accurately represent slices with an

  itk::OrientedImage< TPixel, 2 >

because then the origin is in R^2 and the direction matrix is 2x2. So
when you use the ExtractImageFilter you throw away information.

Frank

On Mon, Mar 16, 2009 at 07:43:55PM +0000, Constantine Zakkaroff wrote:
> Dear ALL,
>
> I'm looking for a way of extracting slices of arbitrary orientation 
> specified by any given orientation matrix. At the moment I can do it in two 
> different ways, but I suspect there should be a better way.
>
> My first way is to resample the volume with the given orientation matrix 
> and then use ExtractImageFilter to get a slice.
>
> The other way I tried is a bit more involved but it is closer to what I 
> need. First I construct an OrientedImage with the given orientation matrix 
> and then use ImageRegionIteratorWithIndex to traverse the pixels and for 
> each index obtained from the iterator I get the value of a pixel in the 
> corresponding physical location in the volume.
>
> I'm wondering if it is possible (and sensible) to use something like 
> SpatialOrientationAdapter, but with an ability to set arbitrary orientation 
> matrix so that the image appears already resampled if I want to extract a 
> slice form it?
>
> Or is there a way to cook up something like ExtractImageFilter that would 
> let me extract slices under any arbitrary orientation specified through the 
> given orientation matrix?
>
> Thank you very much,
> Constantine
>
>
>
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> 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://www.itk.org/mailman/listinfo/insight-users