Impact of beam‐hardening corrections on proton relative stopping power estimates from single‐ and dual‐energy CT

A major contributing factor to proton range uncertainty is the conversion of computed tomography (CT) Hounsfield units (HU) to proton relative stopping power (RSP). This uncertainty is heightened in the presence of X‐ray beam‐hardening artifact (BHA), which has two manifestations: cupping and streaking, especially in and near bone tissue. This uncertainty can affect the accuracy of proton RSP calculation for treatment planning in proton radiotherapy. Dual‐energy CT (DECT) and iterative beam‐hardening correction (iBHC) both show promise in mitigating CT BHA. This present work attempts to analyze the relative robustness of iBHC and DECT techniques on both manifestations of BHA. The stoichiometric method for HU to RSP conversion was used for single‐energy CT (SECT) and DECT‐based monochromatic techniques using a tissue substitute phantom. Cupping BHA was simulated by measuring the HU of a bone substitute plug in wax/3D‐printed phantoms of increasing size. Streaking BHA was simulated by placing a solid water plug between two bone plugs in a wax phantom. Finally, the effect of varying calibration phantom size on RSP was calculated in an anthropomorphic head phantom. The RSP decreased −0.002 cm(–1) as phantom size increased for SECT but remained largely constant when iBHC applied or with DECT techniques. The RSP varied a maximum of 2.60% in the presence of streaking BHA in SECT but was reduced to 1.40% with iBHC. For DECT techniques, the maximum difference was 2.40%, reduced to 0.6% with iBHC. Comparing calibration phantoms of 20‐ and 33‐cm diameter, maximum voxel differences of 5 mm in the water‐equivalent thickness were observed in the skull but reduced to 1.3 mm with iBHC. The DECT techniques excelled in mitigating cupping BHA, but streaking BHA still could be observed. The use of iBHC reduced RSP variation with BHA in both SECT and DECT techniques.