Introducing a novel sub‐millimeter lung CT image registration error quantitation tool
Authors: Boyle Peter , Naumann Louise , Lauria Michael , Miller Claudia , Andosca Ryan , Savjani Ricky , O'Connell Dylan , Moghanaki Drew , Barjaktarevic Igor , Goldin Jonathan , Low Daniel A. 
Published: December 2024

Background

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.

Purpose

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.

Methods

λ 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 mm3 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 95th 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.

Results

We found θ to be coincident with the image gradient. We found that ’s spatial component tracked the vessel shifts. We determined the spatial sensitivity limit of to be < 0.2 mm. The noise level and smoothing kernel influenced sensitivity, worsening with increasing noise, and improving with increasing smoothing. For the clinical images, we observed to qualitatively match the absolute difference of intensity in the image pairs and to restrict itself to high gradient regions or regions of visually apparent errors.

Conclusion

detected sub‐millimeter positioning errors between simulated vessels in the presence of typical CT noise. The noise magnitude and choice of noise smoothing kernel were inversely related to sensitivity, implying that study‐specific tuning of the pre‐smoothing kernel may be required. The demonstrated ability in geometric tests of to detect subvoxel DIR errors warrants further evaluation and testing.