I’ve been using itk::ImageSeriesReader to read various DICOM studies, and code similar to that described in the “Process a 2D Slice of a 3D Image” example to extract any of the slices in a series. My code works for all CT and MRI series I’ve thrown at it, with one exception. When I try to extract any images from a multiplanar scout from an MRI examination, itk::ExtractImageFilter throws an exception:
ITK ERROR: ExtractImageFilter(00000186D95E84E0): Invalid submatrix extracted for collapsed direction.
This makes sense, in that the scout series has 19 mixed multiplanar (axial/sagittal/coronal) images, rather than a volume of axial (or coronal or sagittal) images.
Is there a way to tell from the DICOM header, or preferably from the itk::Image<> object, that this is a multiplanar scout? Or is there some other way I should be extracting the images in a series, besides assuming they’re from a 3D volume?
Thanks for any advice,