Geometric calibration of a mobile C-arm for intraoperative cone-beam CT

A geometric calibration method that determines a complete description of source-detector geometry was adapted to a mobile C-arm for cone-beam computed tomography (CBCT). The non-iterative calibration algorithm calculates a unique solution for the positions of the source (X(s), Y(s), Z(s)), detector (X(d), Y(d), Z(d)), piercing point (U(o), V(o)), and detector rotation angles (phi, theta, eta) based on projections of a phantom consisting of two plane-parallel circles of ball bearings encased in a cylindrical acrylic tube. The prototype C-arm system was based on a Siemens PowerMobil modified to provide flat-panel CBCT for image-guided interventions. The magnitude of geometric nonidealities in the source-detector orbit was measured, and the short-term (approximately 4 h) and long-term (approximately 6 months) reproducibility of the calibration was evaluated. The C-arm exhibits large geometric nonidealities due to mechanical flex, with maximum departures from the average semicircular orbit of deltaU(o) = 15.8 mm and deltaV(o) = 9.8 mm (for the piercing point), deltaX and deltaY = 6-8 mm and deltaZ = 1 mm (for the source and detector), and deltaphi approximately 2.9 degrees, deltatheta approximately 1.9 degrees, and delta eta approximately 0.8 degrees (for the detector tilt/rotation). Despite such significant departures from a semicircular orbit, these system parameters were found to be reproducible, and therefore correctable by geometric calibration. Short-term reproducibility was < 0.16 mm (subpixel) for the piercing point coordinates, < 0.25 mm for the source-detector X and Y, < 0.035 mm for the source-detector Z, and < 0.02 degrees for the detector angles. Long-term reproducibility was similarly high, demonstrated by image quality and spatial resolution measurements over a period of 6 months. For example, the full-width at half-maximum (FWHM) in axial images of a thin steel wire increased slightly as a function of the time (delta) between calibration and image acquisition: FWHM=0.62, 0.63, 0.66, 0.71, and 0.72 mm at delta = 0 s, 1 h, 1 day, 1 month, and 6 months, respectively. For ongoing clinical trials in CBCT-guided surgery at our institution, geometric calibration is conducted monthly to provide sufficient three-dimensional (3D) image quality while managing time and workflow considerations of the calibration and quality assurance process. The sensitivity of 3D image quality to each of the system parameters was investigated, as was the tolerance to systematic and random errors in the geometric parameters, showing the most sensitive parameters to be the piercing point coordinates (U(o), V(o)) and in-plane positions of the source (X(s), Y(s)) and detector (X(d), Y(d)). Errors in the out-of-plane position of the source (Z(s)) and detector (Z(d)) and the detector angles (phi, theta, eta) were shown to have subtler effects on 3D image quality.