Introducing a novel sub-millimeter lung CT image registration error quantitation tool

Lung computed tomography (CT) scan image registration is being used for lung function analysis such as ventilation. Given the high sensitivity of functional analyses to image registration errors, an image registration error scoring tool that can measure submillimeter image registration errors is needed. To propose an image registration error scoring tool, termed λ, whose spatial sensitivity can be used to quantify image registration errors in steep image gradient regions under realistic noise conditions. λ compares two images, termed reference and evaluated. The HU and distance scales of both images are normalized by user-selected scaling criteria. For each voxel in the reference image, the 4D Euclidian distances between the reference voxel and the nearby evaluated voxels are calculated, and the minimum of these distances is . We tested in simulated individual blood vessels comprised of 1, 3, and 5 mm diameter cylinders in 1 × 1 × 1 mm voxel images, which were blurred to simulate CT scanner intrinsic resolution and volume averaging. We placed the simulated vessels in a homogeneous background simulating parenchymal tissue density and injected 20, 40, and 60 HU standard deviation Gaussian noise. We used isotropic Gaussian filters with 0.5, 1.0, and 1.5 mm standard deviation kernels to smooth the simulated images. We assessed using reference-evaluated vessel shifts of -1.0 to 1.0 mm in 0.05 mm steps via rigid translational and rotational deformations. We examined whether tracked the translation vector via its internal spatial component. We restricted to voxels using the angle, termed , between the vector and the normalized spatial-distance axes, terming the results the restricted- , , where was hypothesized to be a proxy for image gradient. We determined whether was coincident with the image gradient by examining if the voxels with tracked the evaluated vessels. We used the 95 percentile of , , to determine spatial sensitivity, which we took as a conservative estimate of registration error, by fitting to a modified absolute-value function for each tested rigid translation, noise level, smoothing kernel, and vessel radius combination. We demonstrated the use of on a clinical example consisting of a set of 25 deformably registered free-breathing thoracic CT scans. We visually compared the and results against the HU differences between each clinical image pair. We found θ to be coincident with the image gradient. We found that 's spatial component tracked the ve