Motion correction of 3D dynamic contrast-enhanced ultrasound imaging without anatomical B-Mode images: Pilot evaluation in eight patients.

Abstract

Dynamic contrast-enhanced ultrasound (DCE-US) is highly susceptible to motion artifacts arising from patient movement, respiration, and operator handling and experience. Motion artifacts can be especially problematic in the context of perfusion quantification. In conventional 2D DCE-US, motion correction (MC) algorithms take advantage of accompanying side-by-side anatomical B-Mode images that contain time-stable features. However, current commercial models of 3D DCE-US do not provide side-by-side B-Mode images, which makes MC challenging. This work introduces a novel MC algorithm for 3D DCE-US and assesses its efficacy when handling clinical data sets. In brief, the algorithm uses a pyramidal approach whereby short temporal windows consisting of three consecutive frames are created to perform local registrations, which are then registered to a master reference derived from a weighted average of all frames. We applied the algorithm to imaging studies from eight patients with metastatic lesions in the liver and assessed improvements in original versus motion corrected 3D DCE-US cine using: (i) frame-to-frame volumetric overlap of segmented lesions, (ii) normalized correlation coefficient (NCC) between frames (similarity analysis), and (iii) sum of squared errors (SSE), root-mean-squared error (RMSE), and r-squared (R2) quality-of-fit from fitted time-intensity curves (TIC) extracted from a segmented lesion. We noted improvements in frame-to-frame lesion overlap across all patients, from 68%±13% without correction to 83%±3% with MC (p=0.023). Frame-to-frame similarity as assessed by NCC also improved on two different sets of time points from 0.694±0.057 (original cine) to 0.862±0.049 (corresponding MC cine) and 0.723±0.066 to 0.886±0.036 (p ≤ 0.001 for both). TIC analysis displayed a significant decrease in RMSE (p=0.018) and a significant increase in R2 goodness-of-fit (p=0.029) for the patient cohort. Overall, results suggest decreases in 3D DCE-US motion after applying the proposed algorithm.