Guillaume Janssens

Evaluation of nonrigid registration models for interfraction dose accumulation in radiotherapy.

PURPOSE Interfraction dose accumulation is necessary to evaluate the dose distribution of an entire course of treatment by adding up multiple dose distributions of different treatment fractions. This accumulation of dose distributions is not straightforward as changes in the patient anatomy may occur during treatment. For this purpose, the accuracy of nonrigid registration methods is assessed for dose accumulation based on the calculated deformations fields. METHODS A phantom study using a deformable cubic silicon phantom with implanted markers and a cylindrical silicon phantom with MOSFET detectors has been performed. The phantoms were deformed and images were acquired using a cone-beam CT imager. Dose calculations were performed on these CT scans using the treatment planning system. Nonrigid CT-based registration was performed using two different methods, the Morphons and Demons. The resulting deformation field was applied on the dose distribution. For both phantoms, accuracy of the registered dose distribution was assessed. For the cylindrical phantom, also measured dose values in the deformed conditions were compared with the dose values of the registered dose distributions. Finally, interfraction dose accumulation for two treatment fractions of a patient with primary rectal cancer has been performed and evaluated using isodose lines and the dose volume histograms of the target volume and normal tissue. RESULTS A significant decrease in the difference in marker or MOSFET position was observed after nonrigid registration methods (p < 0.001) for both phantoms and with both methods, as well as a significant decrease in the dose estimation error (p < 0.01 for the cubic phantom and p < 0.001 for the cylindrical) with both methods. Considering the whole data set at once, the difference between estimated and measured doses was also significantly decreased using registration (p < 0.001 for both methods). The patient case showed a slightly underdosed planning target volume and an overdosed bladder volume due to anatomical deformations. CONCLUSIONS Dose accumulation using nonrigid registration methods is possible using repeated CT imaging. This opens possibilities for interfraction dose accumulation and adaptive radiotherapy to incorporate possible differences in dose delivered to the target volume and organs at risk due to anatomical deformations.

Diffeomorphic registration of images with variable contrast enhancement

Nonrigid image registration is widely used to estimate tissue deformations in highly deformable anatomies. Among the existing methods, nonparametric registration algorithms such as optical flow, or Demons, usually have the advantage of being fast and easy to use. Recently, a diffeomorphic version of the Demons algorithm was proposed. This provides the advantage of producing invertible displacement fields, which is a necessary condition for these to be physical. However, such methods are based on the matching of intensities and are not suitable for registering images with different contrast enhancement. In such cases, a registration method based on the local phase like the Morphons has to be used. In this paper, a diffeomorphic version of the Morphons registration method is proposed and compared to conventional Morphons, Demons, and diffeomorphic Demons. The method is validated in the context of radiotherapy for lung cancer patients on several 4D respiratory-correlated CT scans of the thorax with and without variable contrast enhancement.

First Clinical Investigation of Cone Beam Computed Tomography and Deformable Registration for Adaptive Proton Therapy for Lung Cancer

PURPOSE An adaptive proton therapy workflow using cone beam computed tomography (CBCT) is proposed. It consists of an online evaluation of a fast range-corrected dose distribution based on a virtual CT (vCT) scan. This can be followed by more accurate offline dose recalculation on the vCT scan, which can trigger a rescan CT (rCT) for replanning. METHODS AND MATERIALS The workflow was tested retrospectively for 20 consecutive lung cancer patients. A diffeomorphic Morphon algorithm was used to generate the lung vCT by deforming the average planning CT onto the CBCT scan. An additional correction step was applied to account for anatomic modifications that cannot be modeled by deformation alone. A set of clinical indicators for replanning were generated according to the water equivalent thickness (WET) and dose statistics and compared with those obtained on the rCT scan. The fast dose approximation consisted of warping the initial planned dose onto the vCT scan according to the changes in WET. The potential under- and over-ranges were assessed as a variation in WET at the targets distal surface. RESULTS The range-corrected dose from the vCT scan reproduced clinical indicators similar to those of the rCT scan. The workflow performed well under different clinical scenarios, including atelectasis, lung reinflation, and different types of tumor response. Between the vCT and rCT scans, we found a difference in the measured 95% percentile of the over-range distribution of 3.4 ± 2.7 mm. The limitations of the technique consisted of inherent uncertainties in deformable registration and the drawbacks of CBCT imaging. The correction step was adequate when gross errors occurred but could not recover subtle anatomic or density changes in tumors with complex topology. CONCLUSIONS A proton therapy workflow based on CBCT provided clinical indicators similar to those using rCT for patients with lung cancer with considerable anatomic changes.

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation. METHODS Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT. RESULTS The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy]. CONCLUSIONS In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

First test of the prompt gamma ray timing method with heterogeneous targets at a clinical proton therapy facility.

Ion beam therapy promises enhanced tumour coverage compared to conventional radiotherapy, but particle range uncertainties significantly blunt the achievable precision. Experimental tools for range verification in real-time are not yet available in clinical routine. The prompt gamma ray timing method has been recently proposed as an alternative to collimated imaging systems. The detection times of prompt gamma rays encode essential information about the depth-dose profile thanks to the measurable transit time of ions through matter. In a collaboration between OncoRay, Helmholtz-Zentrum Dresden-Rossendorf and IBA, the first test at a clinical proton accelerator (Westdeutsches Protonentherapiezentrum Essen, Germany) with several detectors and phantoms is performed. The robustness of the method against background and stability of the beam bunch time profile is explored, and the bunch time spread is characterized for different proton energies. For a beam spot with a hundred million protons and a single detector, range differences of 5 mm in defined heterogeneous targets are identified by numerical comparison of the spectrum shape. For higher statistics, range shifts down to 2 mm are detectable. A proton bunch monitor, higher detector throughput and quantitative range retrieval are the upcoming steps towards a clinically applicable prototype. In conclusion, the experimental results highlight the prospects of this straightforward verification method at a clinical pencil beam and settle this novel approach as a promising alternative in the field of in vivo dosimetry.

Real-time proton beam range monitoring by means of prompt-gamma detection with a collimated camera

Prompt-gamma profile was measured at WPE-Essen using 160 MeV protons impinging a movable PMMA target. A single collimated detector was used with time-of-flight (TOF) to reduce the background due to neutrons. The target entrance rise and the Bragg peak falloff retrieval precision was determined as a function of incident proton number by a fitting procedure using independent data sets. Assuming improved sensitivity of this camera design by using a greater number of detectors, retrieval precisions of 1 to 2 mm (rms) are expected for a clinical pencil beam. TOF improves the contrast-to-noise ratio and the performance of the method significantly.

Measurement of prompt gamma profiles in inhomogeneous targets with a knife-edge slit camera during proton irradiation.

Proton and ion beam therapies become increasingly relevant in radiation therapy. To fully exploit the potential of this irradiation technique and to achieve maximum target volume conformality, the verification of particle ranges is highly desirable. Many research activities focus on the measurement of the spatial distributions of prompt gamma rays emitted during irradiation. A passively collimating knife-edge slit camera is a promising option to perform such measurements. In former publications, the feasibility of accurate detection of proton range shifts in homogeneous targets could be shown with such a camera. We present slit camera measurements of prompt gamma depth profiles in inhomogeneous targets. From real treatment plans and their underlying CTs, representative beam paths are selected and assembled as one-dimensional inhomogeneous targets built from tissue equivalent materials. These phantoms have been irradiated with monoenergetic proton pencil beams. The accuracy of range deviation estimation as well as the detectability of range shifts is investigated in different scenarios. In most cases, range deviations can be detected within less than 2 mm. In close vicinity to low-density regions, range detection is challenging. In particular, a minimum beam penetration depth of 7 mm beyond a cavity is required for reliable detection of a cavity filling with the present setup. Dedicated data post-processing methods may be capable of overcoming this limitation.

Analytical computation of prompt gamma ray emission and detection for proton range verification

A prompt gamma (PG) slit camera prototype recently demonstrated that Bragg Peak position in a clinical proton scanned beam could be measured with 1-2 mm accuracy by comparing an expected PG detection profile to a measured one. The computation of the expected PG detection profile in the context of a clinical framework is challenging but must be solved before clinical implementation. Obviously, Monte Carlo methods (MC) can simulate the expected PG profile but at prohibitively long calculation times. We implemented a much faster method that is based on analytical processing of precomputed MC data that would allow practical evaluation of this range monitoring approach in clinical conditions. Reference PG emission profiles were generated with MC simulations (PENH) in targets consisting of either (12)C, (14)N, (16)O, (31)P or (40)Ca, with 10% of (1)H. In a given geometry, the local PG emission can then be derived by adding the contribution of each element, according to the local energy of the proton obtained by continuous slowing down approximation and the local composition. The actual incident spot size is taken into account using an optical model fitted to measurements and by super sampling the spot with several rays (up to 113). PG transport in the patient/camera geometries and the detector response are modelled by convolving the PG production profile with a transfer function. The latter is interpolated from a database of transfer functions fitted to MC data (PENELOPE) generated for a photon source in a cylindrical phantom with various radiuses and a camera placed at various positions. As a benchmark, the analytical model was compared to MC and experiments in homogeneous and heterogeneous phantoms. Comparisons with MC were also performed in a thoracic CT. For all cases, the analytical model reproduced the prediction of the position of the Bragg peak computed with MC within 1 mm for the camera in nominal configuration. When compared to measurements, the shape of the profiles was well reproduced and agreement for the estimation of the position of the Bragg peak was within 2.7 mm on average (1.4 mm standard deviation). On a non-optimized MATLAB code, computation time with the analytical model is between 0.3 to 10 s depending on the number of rays simulated per spot. The analytical model can be further used to determine which spots are the best candidates to evaluate the range in clinical conditions and eventually correct for over- and under-shoots depending on the acquired PG profiles.

Residual metabolic tumor activity after chemo-radiotherapy is mainly located in initially high FDG uptake areas in rectal cancer.

PURPOSE Recent literature suggests that tumor cells and areas within tumors with a high initial FDG uptake might be more resistant to (chemo)radiotherapy ((C)RT). This study was undertaken to test this hypothesis in rectal cancer using rigid and non-rigid image registration. PATIENTS AND METHODS Twenty-eight patients, diagnosed with locally advanced rectal cancer and referred for pre-operative treatment with CRT were included in this study. All patients underwent FDG-PET-CT imaging prior to and after CRT. Rigid and non-rigid image registration was performed to compensate organ deformations between the pre- and post-treatment PET-CT scans. The tumor was contoured on both PET-scans using SUV iso-contouring based on the SBR-method. The voxels with residual increased FDG uptake were studied and correlated to their pre-treatment FDG uptake level. Two SUV-volume-histograms were made based on the pre-treatment PET-data, one for the voxels within the pre-treatment tumor PET-based iso-contour and one for the voxels within the PET-based iso-contour of the residual tumor non-rigidly registered onto the pre-treatment scan. RESULTS For the voxels with a pre-treatment FDG uptake of >50% of SUV(max), 70.6±5.6% of the voxels were still metabolic active in the residual tumor, whereas for voxels with an FDG uptake of <50% of SUV(max) only 51.1±6.7% were present in the metabolic active residual tumor. CONCLUSION This study presents areas in rectal tumors with an initially high FDG uptake to be most likely to show residual disease after CRT. This could indicate a higher (C)RT-resistance for tumor regions with a high FDG uptake prior to treatment.

Assessment of tumor motion reproducibility with audio-visual coaching through successive 4D CT sessions

This study aimed to compare combined audio‐visual coaching with audio coaching alone and assess their respective impact on the reproducibility of external breathing motion and, one step further, on the internal lung tumor motion itself, through successive sessions. Thirteen patients with NSCLC were enrolled in this study. The tumor motion was assessed by three to four successive 4D CT sessions, while the breathing signal was measured from magnetic sensors positioned on the epigastric region. For all sessions, the breathing was regularized with either audio coaching alone (AC, n=5) or combined with a real‐time visual feedback (A/VC, n=8) when tolerated by the patients. Peak‐to‐peak amplitude, period and signal shape of both breathing and tumor motions were first measured. Then, the correlation between the respiratory signal and internal tumor motion over time was evaluated, as well as the residual tumor motion for a gated strategy. Although breathing and tumor motions were comparable between AC and AV/C groups, A/VC approach achieved better reproducibility through sessions than AC alone (mean tumor motion of 7.2 mm±1 vs. 8.6 mm±1.8 mm, and mean breathing motion of 14.9 mm±1.2 mm vs. 13.3 mm±3.7 mm, respectively). High internal/external correlation reproducibility was achieved in the superior‐inferior tumor motion direction for all patients. For the anterior‐posterior tumor motion direction, better correlation reproducibility has been observed when visual feedback has been used. For a displacement‐based gating approach, A/VC might also be recommended, since it led to smaller residual tumor motion within clinically relevant duty cycles. This study suggests that combining real‐time visual feedback with audio coaching might improve the reproducibility of key characteristics of the breathing pattern, and might thus be considered in the implementation of lung tumor radiotherapy. PACS number: 87