A method for assessing voxel correspondence in longitudinal tumor imaging

Tumor characterization employing a voxel-wise analysis of image signal facilitates the determination of the spatial distribution of tumor attributes, and when employed for therapy response assessment offers the promise of greater sensitivity to change than conventional approaches. However, the accuracy of a voxel-wise analysis of change is limited by local registration uncertainties that can disrupt the spatiotemporal correspondence between assessed voxels. We present a method for assessing voxel correspondence strength using a multiresolution local histogram-based measure of image structure similarity. When employed in a longitudinal tumor imaging context, a voxel similarity measure must be robust to intensity variations that can arise from the image acquisition, treatment effects, or changes in underlying disease processes. Consequently, the local histogram-based similarity measure is evaluated for sensitivity to structural change and robustness to intensity variation and is compared against normalized mutual information and normalized cross-correlation.

T1-weighted (T1W) magnetic resonance (MR) images of glioblastoma acquired as part of a longitudinal response assessment study are first rigidly registered, and then similarity between spatially corresponding voxels is evaluated using multiresolution local histograms. Region-based and nonuniform intensity changes of varying magnitude as well as deformations to image structure are applied individually and in combination to the test images. Statistical analysis is used to test for interaction effects between the applied perturbations and the value of the local histogram similarity function. Pair-wise voxel similarity maps are computed for pairs of longitudinal clinical MR image volumes and compared with observed patterns of change on conventional imaging.

The simulations demonstrated that the local histogram measure was robust to intensity modulations applied to increasing region sizes and exhibited a strong negative correlation with the magnitude of local deformation. No statistically significant interaction effects were observed upon the value of the local histogram similarity function when deformation was applied in conjunction with a nonuniform intensity change. Pair-wise voxel similarity maps were consistent with image change observed on T1W MR imaging and revealed patterns of change not apparent in conventional image sequences.

A measure of local histogram image structure similarity can be used to assess the strength of voxel to voxel correspondences independently of intensity nonuniformities. The metric can provide a local estimate of the limits of achievable correspondence underlying the registration and voxel-wise comparison of signal in longitudinal imaging used for assessing tumor response to treatment.