Development and first use of a novel cylindrical ball bearing phantom for 9-DOF geometric calibrations of flat panel imaging devices used in image-guided ion beam therapy

Reduction of cone-beam CT artifacts in a robotic CBCT device using saddle trajectories with integrated infrared tracking

BACKGROUND

Cone beam computed tomography (CBCT) is widely used in many medical fields. However, conventional CBCT circular scans suffer from cone beam (CB) artifacts that limit the quality and reliability of the reconstructed images due to incomplete data.

PURPOSE

Saddle trajectories in theory might be able to improve the CBCT image quality by providing a larger region with complete data. Therefore, we investigated the feasibility and performance of saddle trajectory CBCT scans and compared them to circular trajectory scans.

METHODS

We performed circular and saddle trajectory scans using a novel robotic CBCT scanner (Mobile ImagingRing (IRm); medPhoton, Salzburg, Austria). For the saddle trajectory, the gantry executed yaw motion up to± 10 ∘ $\pm 10^{\circ }$ using motorized wheels driving on the floor. An infrared (IR) tracking device with reflective markers was used for online geometric calibration correction (mainly floor unevenness). All images were reconstructed using penalized least-squares minimization with the conjugate gradient algorithm from RTK with0.5 × 0.5 × 0.5 mm 3 $0.5 \times 0.5\times 0.5 \text{ mm}^3$ voxel size. A disk phantom and an Alderson phantom were scanned to assess the image quality. Results were correlated with the local incompleteness value represented bytan ( ψ ) $\tan (\psi)$ , which was calculated at each voxel as a function of the source trajectory and the voxel's 3D coordinates. We assessed the magnitude of CB artifacts using the full width half maximum (FWHM) of each disk profile in the axial center of the reconstructed images. Spatial resolution was also quantified by the modulation transfer function at 10% (MTF10).

RESULTS

When using the saddle trajectory, the region without CB artifacts was increased from 43 to 190 mm in the SI direction compared to the circular trajectory. This region coincided with low values fortan ( ψ ) $\tan (\psi)$ . Whentan ( ψ ) $\tan (\psi)$ was larger than 0.02, we found there was a linear relationship between the FWHM andtan ( ψ ) $\tan (\psi)$ . For the saddle, IR tracking allowed the increase of MTF10 from 0.37 to 0.98 lp/mm.

CONCLUSIONS

We achieved saddle trajectory CBCT scans with a novel CBCT system combined with IR tracking. The results show that the saddle trajectory provides a larger region with reliable reconstruction compared to the circular trajectory. The proposed method can be used to evaluate other non-circular trajectories.

Reference-free calibration method for asynchronous rotation in robotic CT

BACKGROUND

Geometry calibration for robotic CT system is necessary for obtaining acceptable images under the asynchrony of two manipulators.

OBJECTIVE

We aim to evaluate the impact of different types of asynchrony on images and propose a reference-free calibration method based on a simplified geometry model.

METHODS

We evaluate the impact of different types of asynchrony on images and propose a novel calibration method focused on asynchronous rotation of robotic CT. The proposed method is initialized with reconstructions under default uncalibrated geometry and uses grid sampling of estimated geometry to determine the direction of optimization. Difference between the re-projections of sampling points and the original projection is used to guide the optimization direction. Images and estimated geometry are optimized alternatively in an iteration, and it stops when the difference of residual projections is close enough, or when the maximum iteration number is reached.

RESULTS

In our simulation experiments, proposed method shows better performance, with the PSNR increasing by 2%, and the SSIM increasing by 13.6% after calibration. The experiments reveal fewer artifacts and higher image quality.

CONCLUSION

We find that asynchronous rotation has a more significant impact on reconstruction, and the proposed method offers a feasible solution for correcting asynchronous rotation.

Online Calibration of a Linear Micro Tomosynthesis Scanner

In a linear tomosynthesis scanner designed for imaging histologic samples of several centimeters size at 10 µm resolution, the mechanical instability of the scanning stage (±10 µm) exceeded the resolution of the image system, making it necessary to determine the trajectory of the stage for each scan to avoid blurring and artifacts in the images that would arise from the errors in the geometric information used in 3D reconstruction. We present a method for online calibration by attaching a layer of randomly dispersed micro glass beads or calcium particles to the bottom of the sample stage. The method was based on a parametric representation of the rigid body motion of the sample stage-marker layer assembly. The marker layer was easy to produce and proven effective in the calibration procedure.

Online Geometric Calibration of a Hybrid CT System for Ultrahigh-Resolution Imaging

A hybrid imaging system consisting of a standard computed tomography (CT) scanner and a low-profile photon-counting detector insert in contact with the patient's body has been used to produce ultrahigh-resolution images in a limited volume in chest scans of patients. The detector insert is placed on the patient bed as needed and not attached. Thus, its position and orientation in the scanner is dependent on the patient's position and scan settings. To allow accurate image reconstruction, we devised a method of determining the relative geometry of the detector insert and the CT scanner for each scan using fiducial markers. This method uses an iterative registration algorithm to align the markers in the reconstructed volume from the detector insert to that of the concurrent CT scan. After obtaining precise geometric information of the detector insert relative to the CT scanner, the two complementary sets of images are summed together to create a detailed image with reduced artifacts.

Auto-calibration of cone beam geometries from arbitrary rotating markers using a vector geometry formulation of projection matrices

An efficient method for the determination of the projection geometry of cone beam micro computed tomography systems based on two or more fiducial markers of unknown position within the field of view is derived. By employing the projection matrix formalism commonly used in computer graphics, a very clear presentation of the resulting self consistent calibration problem can be given relating the sought-for matrix to observable parameters of the markers' projections. Both an easy to implement solution procedure for both the unknown projection matrix and the marker assembly as well as the mapping from projection matrices to real space positions and orientations of source and detector relative to the rotational axis are provided. The separate treatment of the calibration problem in terms of projection matrices on the one hand and the independent transformation to a more intuitive geometry representation on the other hand proves to be very helpful with respect to the discussion of the ambiguities occurring in reference-free calibration. In particular, a link between methods based on knowledge on the sample and those based on knowledge solely on the detector geometry can be drawn. This further provides another intuitive view on the often reported difficulty in the estimation of the detector tilt towards the rotational axis. A simulation study considering 10⁶randomly generated cone beam imaging configurations and fiducial marker distributions within a range of typical scenarios is performed in order to assess the stability of the proposed technique. An experimental example supports the simulation results.

Evaluation of the inter- and intrafraction displacement for head patients treated at the particle therapy centre MedAustron based on the comparison of different commercial immobilisation devices

In December 2016 the clinical operation has started at the particle therapy centre MedAustron, Wiener Neustadt, Austria. Different commercial immobilisation devices are used for head patients. These immobilisation devices are a combination of table tops (Qfix BoS™ Headframe, Elekta HeadStep™), pillows (BoS™ Standard pillow, Moldcare®, HeadStep™ pillow) and thermoplastic masks (Klarity Green™, Qfix Fibreplast™, HeadStep™ iCAST double). For each patient image-guided radiotherapy (IGRT) is performed by acquiring orthogonal X-ray imaging and 2D3D registration and the application of the resulting 6-degree of freedom (DOF) position correction on the robotic couch. The inter- and intrafraction displacement of 101 adult head patients and 27 paediatric sedated head patients were evaluated and compared among each other regarding reproducibility during the entire treatment and stability during each fraction. For the comparison, statistical methods (Shapiro-Wilk test, Mann-Whitney U-test) were applied on the position corrections as well as on the position verifications. The actual planning target volume margins of 3mm (adults) and 2mm (children) were evaluated by applying the van Herk formula on the intrafraction displacement results and performing treatment plan robustness simulations of twelve different translational offset scenarios including a HU uncertainty of 3.5%. Statistically significant differences between the immobilisation devices were found, but they turned out to be clinically irrelevant. The margin calculation for adult head patients resulted in 0.8mm (lateral), 1.2mm (cranio-caudal) and 0.6mm (anterior-posterior), and for paediatric head patients under anaesthesia in 0.8mm (lateral), 0.5mm (cranio-caudal) and 0.9mm (anterior-posterior). Based on these values, robustness evaluations of selected adult head patients and sedated children showed the validity of the currently used PTV margins.

An EPID-based method to determine mechanical deformations in a linear accelerator

PURPOSE

Medical linear accelerators (linac) are delivering increasingly complex treatments using modern techniques in radiation therapy. Complete and precise mechanical QA of the linac is therefore necessary to ensure that there is no unexpected deviation from the gantry's planned course. However, state-of-the-art EPID-based mechanical QA procedures often neglect some degrees of freedom (DOF) like the in-plane rotations of the gantry and imager or the source movements inside the gantry head. Therefore, the purpose of this work is to characterize a 14 DOF method for the mechanical QA of linacs. This method seeks to measure every mechanical deformation in a linac, including source movements, in addition to relevant clinical parameters like mechanical and radiation isocenters.

METHODS

A widely available commercial phantom and a custom-made accessory inserted in the linac's interface mount are imaged using the electronic portal imaging device (EPID) at multiple gantry angles. Then, simulated images are generated using the nominal geometry of the linac and digitized models of the phantoms. The nominal geometry used to generate these images can be modified using 14 DOF (3 rigid rotations and 3 translations for the imager and the gantry, and 2 in-plane translations of the source) and any change will modify the simulated image. The set of mechanical deformations that minimizes the differences between the simulated and measured image is found using a genetic algorithm coupled with a gradient-descent optimizer. Phantom mispositioning and gantry angular offset were subsequently calculated and extracted from the results. Simulations of the performances of the method for different levels of noise in the phantom models were performed to calculate the absolute uncertainty of the measured mechanical deformations. The measured source positions and the center of collimation were used to define the beam central axis and calculate the radiation isocenter position and radius.

RESULTS

After the simultaneous optimization of the 14 DOF, the average distance between the center of the measured and simulated ball bearings on the imager was 0.086 mm. Over the course of a full counter-clockwise gantry rotation, all mechanical deformations were measured, showing sub-millimeter translations and rotations smaller than 1° along every axis. The average absolute uncertainty of the 14 DOF (1 SD) was 0.15 mm or degree. Phantom positioning errors were determined with more than 0.1 mm precision. Errors introduced in the experimental setup like phantom positioning errors, source movements or gantry angular offsets were all successfully detected by our QA method. The mechanical deformations measured are shown to be reproducible over the course of a few weeks and are not sensitive to the experimental setup.

CONCLUSION

This work presents of new method for an accurate mechanical QA of the linacs. It features a 14 DOF model of the mechanical deformations that is both more complete and precise than other available methods. It has demonstrated sub-millimeter accuracy through simulation and experimentation. Introduced errors were successfully detected with high precision.