Simulation of x‐ray‐induced acoustic imaging for absolute dosimetry: Accuracy of image reconstruction methods

Purpose Three‐dimensional in‐vivo dose verification is one of the standing challenges in radiation therapy. X‐ray‐induced acoustic tomography has recently been proposed as an imaging method for use in in‐vivo dosimetry. The aim of this study was to investigate the accuracy of reconstructing three‐dimensional (3D) absolute dose using x‐ray‐induced acoustic tomography. We performed this investigation using two different tomographic dose reconstruction techniques. Methods Two examples of 3D dose reconstruction techniques for x‐ray acoustic imaging are investigated. Dose distributions are calculated for varying field sizes using a clinical treatment planning system. The induced acoustic pressure waves which are generated by the increase in temperature due to the absorption of pulsed MV x‐rays are simulated using an advanced numerical modeling package for acoustic wave propagation in the time domain. Two imaging techniques, back projection and iterative time reversal, are used to reconstruct the 3D dose distribution in a water phantom with open fields. Image analysis is performed and reconstructed depth dose curves from x‐ray acoustic imaging are compared to the depth dose curves calculated from the treatment planning system. Calculated field sizes from the reconstructed dose profiles by back projection and time reversal are compared to the planned field size to determine their accuracy. The iterative time reversal imaging technique is also used to reconstruct dose in an example clinical dose distribution. Image analysis of this clinical test case is performed using the gamma passing rate. In addition, gamma passing rates are used to validate the stopping criteria in the iterative time reversal method. Results Water phantom simulations showed that back projection does not adequately reconstruct the shape and intensity of the depth dose. When compared to the depth of maximum dose calculated by a treatment planning system, the maximum dose depth by back projection is shifted deeper by 55 and 75 mm for 4 × 4 cm and 10 × 10 cm field sizes, respectively. The reconstructed depth dose by iterative time reversal accurately agrees with the planned depth dose for a 4 × 4 cm field size and is shifted deeper by 12 mm for the 10 × 10 cm field size. When reconstructing field sizes, the back projection method leads to 18% and 35% larger sizes for the 4 × 4 cm and 10 × 10 cm fields, respectively, whereas the iterative time reversal method reconstructs both field sizes with