ITK Release 4/Image Class Hierarchy Refactoring

Image Class Hierarchy Refactoring

Motivation

The goal of this proposal is to add support for additional image types in ITK.

ITK currently supports several types of image:

• Image
• VectorImage
• SpecialCoordinatesImage
  • PhasedArray3DSpecialCoordinatesImage
• BloxImage
• SparseImage
• LabelMap
• ImageAdaptors and subclasses

With the exception of the PhasedArray3DSpecialCoordinatesImage, all the image types assume that the image topology is an n-dimensional lattice and that voxels have regular spacing in each dimension. While it makes sense to preserve the topology assumption in this refactoring, the assumption of regular spacing should be relaxed wherever possible. The regular spacing assumption covers many cases in biomedical imaging, but it hinders ITK's use in certain applications such as microscopy and remote sensing images. For example, certain motorized stages used in confocal microscopes fail to achieve regular spacing in z when acquiring a 3D image (often called a stack of 2D images at different focal depths). Helpfully, these stages report the positions of z-planes acquired during image stack collection. Assuming a regular z-spacing in place of the actual z-plane positions may result in significant errors in image analysis algorithms run on the image.

Goals

In particular, this proposal aims to provide several new kinds of image classes:

• Digital Image - no physical embedding; raw lattice of voxel values
• Regular Image - analog of ITK's current Image with regular spacing
• Rectilinear Image - a physical image with support for irregular spacing in x, y, and z

Secondary Goals

Merge VectorImages into Images. An Image is merely a VectorImage with a single component per voxel.

Requirements

• Make changes fully backwards compatible (all examples and tests should compile and run without modification and all tests should pass)

Challenges and Potential Pitfalls

• ImageAdaptors should present the correct interface for the image type they adapt. This can be accomplished with partial template specialization of the ImageAdaptor class to avoid modifying any of the existing adaptors.