On Wed, Aug 1, 2012 at 11:10 PM, 龍龍 <anjingdelonglong-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:

HI ALL

i'm doing my project using ITK. but i got a problem about the extract of 3D image.

i'm using vtkRegionofInterestImageFilter. follow the example:

ImageType::IndexType start;
start[0] = …
start[1] = …
start[2] = …

ImageType::SizeType size;
size[0] = …
size[1] = …
size[2] = …

ImageType::RegionType desiredRegion;
desiredRegion.SetSize(size);
desiredRegion.SetIndex(start);

BUT, i hope that the start point that i gave can be the center point of the cube. not the cornor.

what should i do . thank u.

--
Ryuu

You just have to subtract half of the 'size' from the 'center' you want so you can pass 'start' as the corner:

itk::Index<2> corner = center - size/2; (I think you'll have to do this element-wise)
ImageType::RegionType desiredRegion(corner, size);