Search

The aim of this chapter is twofold. On the one hand, it describes a complete simulation method to evaluate the performance of a Doppler LIDAR for aircraft wake vortices detection near the ground. The principal interest of the simulation presented here is to combine a LASER beam propagation methodwith fluid dynamics simulations. On the other hand, it proposes effective wind velocity map reconstruction algorithms.

Purpose Ventilation-induced tumour motion remains a challenge for the accuracy of proton therapy treatments in lung patients. We investigated the feasibility of using a 4D virtual CT (4D-vCT) approach based on deformable image registration (DIR) and motion-aware 4D CBCT reconstruction (MA-ROOSTER) to enable accurate daily proton dose calculation using a gantry-mounted CBCT scanner tailored to proton therapy. Methods Ventilation correlated data of 10 breathing phases were acquired from a porcine ex-vivo functional lung phantom using CT and CBCT. 4D-vCTs were generated by (1) DIR of the mid-position 4D-CT to the mid-position 4D-CBCT (reconstructed with the MA-ROOSTER) using a diffeomorphic Morphons algorithm and (2) subsequent propagation of the obtained mid-position vCT to the individual 4D-CBCT phases. Proton therapy treatment planning was performed to evaluate dose calculation accuracy of the 4D-vCTs. A robust treatment plan delivering a nominal dose of 60 Gy was generated on the average intensity image of the 4D-CT for an approximated internal target volume (ITV). Dose distributions were then recalculated on individual phases of the 4D-CT and the 4D-vCT based on the optimized plan. Dose accumulation was performed for 4D-vCT and 4D-CT using DIR of each phase to the mid position, which was chosen as reference. Dose based on the 4D-vCT was then evaluated against the dose calculated on 4D-CT both, phase-by-phase as well as accumulated, by comparing dose volume histogram (DVH) values (Dmean, D2%, D98%, D95%) for the ITV, and by a 3D-gamma index analysis (global, 3%/3 mm, 5 Gy, 20 Gy and 30 Gy dose thresholds). Results Good agreement was found between the 4D-CT and 4D-vCT-based ITV-DVH curves. The relative differences ((CT-vCT)/CT) between accumulated values of ITV Dmean, D2%, D95% and D98% for the 4D-CT and 4D-vCT-based dose distributions were −0.2%, 0.0%, −0.1% and −0.1%, respectively. Phase specific values varied between −0.5% and 0.2%, −0.2% and 0.5%, −3.5% and 1.5%, and −5.7% and 2.3%. The relative difference of accumulated Dmean over the lungs was 2.3% and Dmean for the phases varied between −5.4% and 5.8%. The gamma pass-rates with 5 Gy, 20 Gy and 30 Gy thresholds for the accumulated doses were 96.7%, 99.6% and 99.9%, respectively. Phase-by-phase comparison yielded pass-rates between 86% and 97%, 88% and 98%, and 94% and 100%. Conclusions Feasibility of the suggested 4D-vCT workflow using proton therapy specific imaging equipment was shown. Results indicate the potential of the method to be applied for daily 4D proton dose estimation.

Purpose: The aim of this work is to propose a general and simple procedure for the calibration and validation of kilo-voltage cone-beam CT (kV CBCT) models against experimental data. Methods: The calibration and validation of the CT model is a two-step procedure: the source model then the detector model. The source is described by the direction dependent photon energy spectrum at each voltage while the detector is described by the pixel intensity value as a function of the direction and the energy of incident photons. The measurements for the source consist of a series of dose measurements in air performed at each voltage with varying filter thicknesses and materials in front of the x-ray tube. The measurements for the detector are acquisitions of projection images using the same filters and several tube voltages. The proposed procedure has been applied to calibrate and assess the accuracy of simple models of the source and the detector of three commercial kV CBCT units. If the CBCT system models had been calibrated differently, the current procedure would have been exclusively used to validate the models. Several high-purity attenuation filters of aluminum, copper, and silver combined with a dosimeter which is sensitive to the range of voltages of interest were used. A sensitivity analysis of the model has also been conducted for each parameter of the source and the detector models. Results: Average deviations between experimental and theoretical dose values are below 1.5% after calibration for the three x-ray sources. The predicted energy deposited in the detector agrees with experimental data within 4% for all imaging systems. Conclusions: The authors developed and applied an experimental procedure to calibrate and validate any model of the source and the detector of a CBCT unit. The present protocol has been successfully applied to three x-ray imaging systems. The minimum requirements in terms of material and equipment would make its implementation suitable in most clinical environments.

The relative stopping power (RSP) uncertainty is the largest contributor to the range uncertainty in proton therapy. The purpose of this work was to develop a systematic method that yields accurate and patient-specific RSPs by combining (1) pre-treatment x-ray CT and (2) daily proton radiography of the patient. The method was formulated as a penalized least squares optimization problem (argmin([Formula: see text])). The parameter A represents the cumulative path-length crossed by the proton in each material, separated by thresholding on the HU. The material RSPs (water equivalent thickness/physical thickness) are denoted by x. The parameter b is the list-mode proton radiography produced using Geant4 simulations. The problem was solved using a non-negative linear-solver with [Formula: see text]. A was computed by superposing proton trajectories calculated with a cubic or linear spline approach to the CT. The material's RSP assigned in Geant4 were used for reference while the clinical HU-RSP calibration curve was used for comparison. The Gammex RMI-467 phantom was first investigated. The standard deviation between the estimated material RSP and the calculated RSP is 0.45%. The robustness of the techniques was then assessed as a function of the number of projections and initial proton energy. Optimization with two initial projections yields precise RSP (⩽1.0%) for 330 MeV protons. 250 MeV protons have shown higher uncertainty (⩽2.0%) due to the loss of precision in the path estimate. Anthropomorphic phantoms of the head, pelvis, and lung were subsequently evaluated. Accurate RSP has been obtained for the head ([Formula: see text]), the lung ([Formula: see text]) and the pelvis ([Formula: see text]). The range precision has been optimized using the calibration curves obtained with the algorithm, yielding a mean [Formula: see text] difference to the reference of 0.11 ±0.09%, 0.28 ± 0.34% and [Formula: see text] in the same order. The solution's accuracy is limited by the assumed HU/RSP bijection, neglecting inherent degeneracy. The proposed formulation of the problem with prior knowledge x-ray CT demonstrates potential to increase the accuracy of present RSP estimates.

Multiple Coulomb scattering (MCS) is the largest contributor to blurring in proton imaging. In this work, we developed a maximum likelihood least squares estimator that improves proton radiography's spatial resolution. The water equivalent thickness (WET) through projections defined from the source to the detector pixels were estimated such that they maximizes the likelihood of the energy loss of every proton crossing the volume. The length spent in each projection was calculated through the optimized cubic spline path estimate. The proton radiographies were produced using Geant4 simulations. Three phantoms were studied here: a slanted cube in a tank of water to measure 2D spatial resolution, a voxelized head phantom for clinical performance evaluation as well as a parametric Catphan phantom (CTP528) for 3D spatial resolution. Two proton beam configurations were used: a parallel and a conical beam. Proton beams of 200 and 330 MeV were simulated to acquire the radiography. Spatial resolution is increased from 2.44 lp cm-1 to 4.53 lp cm-1 in the 200 MeV beam and from 3.49 lp cm-1 to 5.76 lp cm-1 in the 330 MeV beam. Beam configurations do not affect the reconstructed spatial resolution as investigated between a radiography acquired with the parallel (3.49 lp cm-1 to 5.76 lp cm-1) or conical beam (from 3.49 lp cm-1 to 5.56 lp cm-1). The improved images were then used as input in a photon tomography algorithm. The proton CT reconstruction of the Catphan phantom shows high spatial resolution (from 2.79 to 5.55 lp cm-1 for the parallel beam and from 3.03 to 5.15 lp cm-1 for the conical beam) and the reconstruction of the head phantom, although qualitative, shows high contrast in the gradient region. The proposed formulation of the optimization demonstrates serious potential to increase the spatial resolution (up by 65[Formula: see text]) in proton radiography and greatly accelerate proton computed tomography reconstruction.

The aim of this work is to propose a general and simple procedure for the calibration and validation of kilo-voltage cone-beam CT (kV CBCT) models against experimental data.The calibration and validation of the CT model is a two-step procedure: the source model then the detector model. The source is described by the direction dependent photon energy spectrum at each voltage while the detector is described by the pixel intensity value as a function of the direction and the energy of incident photons. The measurements for the source consist of a series of dose measurements in air performed at each voltage with varying filter thicknesses and materials in front of the x-ray tube. The measurements for the detector are acquisitions of projection images using the same filters and several tube voltages. The proposed procedure has been applied to calibrate and assess the accuracy of simple models of the source and the detector of three commercial kV CBCT units. If the CBCT system models had been calibrated differently, the current procedure would have been exclusively used to validate the models. Several high-purity attenuation filters of aluminum, copper, and silver combined with a dosimeter which is sensitive to the range of voltages of interest were used. A sensitivity analysis of the model has also been conducted for each parameter of the source and the detector models.Average deviations between experimental and theoretical dose values are below 1.5% after calibration for the three x-ray sources. The predicted energy deposited in the detector agrees with experimental data within 4% for all imaging systems.The authors developed and applied an experimental procedure to calibrate and validate any model of the source and the detector of a CBCT unit. The present protocol has been successfully applied to three x-ray imaging systems. The minimum requirements in terms of material and equipment would make its implementation suitable in most clinical environments.

In this paper, we present the development of an axial aircraft wake vortex light detection and ranging (LIDAR) sensor, working in Mie scattering regime, based on pulsed 1.5-mu m high-brightness large-core fiber amplifier. An end-to-end Doppler heterodyne LIDAR simulator is used for the LIDAR design. The simulation includes the observation geometry, the wake vortex velocity image, the scanning pattern, the LIDAR instrument, the wind turbulence outside the vortex, and the signal processing. An innovative high-brightness pulsed 1.5-mum laser source is described, based on a master oscillator power fiber amplifier (MOPFA) architecture with a large-core fiber. The obtained beam quality is excellent (M 2 = 1.3), and achieved pulsed energy is 120 muJ with a pulse repetition frequency of 12 kHz and a pulse duration of 800 ns. A Doppler heterodyne LIDAR is developed based on this laser source with a high-isolation free-space circulator. The LIDAR includes a real-time display of the wind field. Wind dispersion is postprocessed. Field tests carried out at Orly airport in April 2008 are reported. Axial aircraft wake vortex signatures have been successfully observed and acquired at a range of 1.2 km with axial resolution of 75 m for the first time with fiber laser source.

International audience; Purpose: To develop an open-source toolkit for fast cone-beam CT reconstruction based on the Insight Toolkit. Methods: We have started the Reconstruction Toolkit (RTK, http://www.openrtk.org), an open-source toolkit for cone-beam CT reconstruction, based on the Insight Toolkit (ITK, http://www.itk.org/) and using GPU code extracted from Plastimatch (http://www.plastimatch.org/). RTK is developed by an open consortium (see affiliations) under the non- contaminating Apache 2.0 license. The quality of the platform is daily checked with regression tests in partnership with Kitware which already supports ITK. Results: Several features are already available: Elekta, Varian and IBA inputs, multi-threaded Feldkamp-David-Kress reconstruction on CPU and GPU, Parker short scan weighting, multi-threaded CPU and GPU forward projectors, etc. Each feature is either accessible through command line tools or C++ classes that can be included in independent software. A MIDAS community (http://midas3.kitware.com) has been opened to provide CatPhan datasets of several vendors (Elekta, Varian and IBA). RTK will be used in the upcoming cone-beam CT scanner developed by IBA for proton therapy rooms. Many features are under development: new input format support, iterative reconstruction, hybrid Monte Carlo / deterministic CBCT simulation, etc. Conclusions: RTK has been built to freely share tomographic reconstruction development between researchers and is open for new contributions.

The uncertainties in water equivalent thickness (WET) and accuracy of dose estimation using a virtual CT (vCT), generated from deforming the planning CT (pCT) onto the daily cone-beam CT (CBCT), were comprehensively evaluated in the context of lung malignancies and passive scattering proton therapy. The validation methodology utilized multiple CBCT datasets to generate the vCTs of twenty lung cancer patients. A correction step was applied to the vCTs to account for anatomical modifications that could not be modeled by deformation alone. The CBCT datasets included a regular CBCT (rCBCT) and synthetic CBCTs created from the rCBCT and rescan CT (rCT), which minimized the variation in setup between the vCT and the gold-standard image (i.e., rCT). The uncertainty in WET was defined as the voxelwise difference in WET between vCT and rCT, and calculated in 3D (planning target volume, PTV) and 2D (distal and proximal surfaces). The uncertainty in WET based dose warping was defined as the difference between the warped dose and a forward dose recalculation on the rCT. The overall root mean square (RMS) uncertainty in WET was 3.6 ± 1.8, 2.2 ± 1.4 and 3.3 ± 1.8 mm for the distal surface, proximal surface and PTV, respectively. For the warped dose, the RMS uncertainty of the voxelwise dose difference was 6% ± 2% of the maximum dose (%mD), using a 20% cut-off. The rCBCT resulted in higher uncertainties due to setup variability with the rCT; the uncertainties reported with the two synthetic CBCTs were similar. The vCT followed by a correction step was found to be an accurate alternative to rCT.

This study concerns the modeling and the design of a monostatic heterodyne pulsed LIDAR. The heart of the system is constituted of a 1.55 mum able to produce high pulse repetition frequency. The aim of this work is to assess its efficiency to perform accurate wind speed measurements in the low atmospheric boundary layer, from a 2-D scanning pattern, in the presence of refractive turbulence. A complete LIDAR numerical simulation technique has been developed. Its main originality is the integration of both optical and fluid dynamics numerical methods to take into account the signal coherence loss due to refractive turbulence and speckle effect as well as the fine structures of the wind field. The wind speed profiles along each line-of-sight are retrieved from the return signal using a low-order autoregressive model. An adequate averaging model is then used estimate horizontal components of the wind speed for altitudes up to 150 m.

Purpose: To develop a novel technique to enhance the image contrast of clinical cone beam CT projections and extract respiratory signals based on anatomical motion using the modified Amsterdam Shroud (AS) method to benefit image guided radiation therapy. Methods: Thoracic cone beam CT projections acquired prior to treatment were preprocessed to increase their contrast for better respiratory signal extraction. Air intensity on raw images was firstly estimated and then applied to correct the projections to generate new attenuation images that were subsequently improved with deeper anatomy feature enhancement through taking logarithm operation, derivative along superior-inferior direction, respectively. All pixels on individual post-processed two dimensional images were horizontally summed to one column and all projections were combined side by side to create an AS image from which patient’s respiratory signal was extracted. The impact of gantry rotation on the breathing signal rendering was also investigated. Ten projection image sets from five lung cancer patients acquired with the Varian Onboard Imager on 21iX Clinac (Varian Medical Systems, Palo Alto, CA) were employed to assess the proposed technique. Results: Application of the air correction on raw projections showed that more than an order of magnitude of contrast enhancement was achievable. The typical contrastmore » on the raw projections is around 0.02 while that on attenuation images could greater than 0.5. Clear and stable breathing signal can be reliably extracted from the new images while the uncorrected projection sets failed to yield clear signals most of the time. Conclusion: Anatomy feature plays a key role in yielding breathing signal from the projection images using the AS technique. The air correction process facilitated the contrast enhancement significantly and attenuation images thus obtained provides a practical solution to obtaining markerless breathing motion tracking.« less

Purpose: To demonstrate that the discontinuous nature of the CT number to stopping power ratio (SPR) calibration curve, combined with the presence of uncorrelated zero-mean Gaussian CT noise, leads to non-negligible and tissue-dependent systematic errors in SPRs and proton range, typically not taken into account in usual safety margins for proton therapy. Methods: Increased systematic errors with noise standard deviation have first been observed in proton range Monte-Carlo simulations with stoichiometric calibrations, whereas only zero-mean random errors were expected. Their existence has then been proved analytically for arbitrary calibration curves and material distributions along the proton path and validated through continuous slowing down approximation (CSDA) simulations. Their importance relative to the other sources of uncertainty has then been estimated in head-and-neck, lung, and pelvis patient data for multiple beam orientations. CT noise has first been reduced using a double-pass median filtering approach and a Gaussian noise has then been added to obtain total standard deviations between 10 to 40 HU. Results: This study provides close form equations for the systematic error and uncertainty on SPR and proton range due to uncorrelated noise. They have shown to accurately match CSDA simulation results with realistic calibration curves and material distributions. Depending on the tissue distribution and the position of the discontinuities along the curve the resulting effect on range varies but has shown never to cancel out completely as opposed to common beliefs. The analysis performed on patient data with clinical calibration curves has confirmed that fact with estimated systematic range errors of 0.2–0.5% and uncertainties (4 σ) between 0.5 and 1% with typical CT noise levels. Conclusion: A new source of SPR and range systematic errors has been highlighted and proved not to be negligible compared to the 3.5% uncertainty reference value used for safety margin design This study is linked to a public partnership between UCL and IBA funded by the Walloon region under convention number 1017266 and 1217662

Purpose: The material relative stopping power (RSP) uncertainty is the highest contributor to the range uncertainty in proton therapy. The purpose of this work is to develop a robust and systematic method that yields accurate, patient specific, RSP by combining 1) pre-treatment x-ray CT and 2) daily proton radiograph of the patient. Methods: The method is formulated as a linear least-square optimization problem (min||Ax-B||2). The parameter A represents the pathlength crossed by the proton in each material. The RSPs for the materials (water equivalent thickness (WET)/physical thickness) are denoted by x. B is the proton radiograph expressed as WET crossed. The problem is minimized using a convex-conic optimization algorithm with xi

Purpose: This study assesses the potential use of CBCT images in adaptive protontherapy by estimating the contribution of the main sources of noise and calibration errors to the proton range uncertainty. Methods: Measurements intended to highlight each particular source have been achieved by adapting either the testbench configuration, e.g. use of filtration, fan-beam collimation, beam stop arrays, phantoms and detector reset light, or the sequence of correction algorithms including water precorrection. Additional Monte-Carlo simulations have been performed to complement these measurements, especially for the beam hardening and the scatter cases. Simulations of proton beams penetration through the resulting images have then been carried out to quantify the range change due to these effects. The particular case of a brain irradiation is considered mainly because of the multiple effects that the skull bones have on the internal soft tissues. Results: On top of the range error sources is the undercorrection of scatter. Its influence has been analyzed from a comparison of fan-beam and full axial FOV acquisitions. In this case, large range errors of about 12 mm can be reached if the assumption is made that the scatter has only a constant contribution over the projection images. Even the detector lag,more » which a priori induces a much smaller effect, has been shown to contribute for up to 2 mm to the overall error if its correction only aims at reducing the skin artefact. This last result can partially be understood by the larger interface between tissues and bones inside the skull. Conclusion: This study has set the basis of a more systematical analysis of the effect CBCT noise on range uncertainties based on a combination of measurements, simulations and theoretical results. With our method, even more subtle effects such as the cone-beam artifact or the detector lag can be assessed. SBR and JOR are financed by iMagX, a public-private partnership between the region Wallone of Belgium and IBA under convention #1217662.« less

Purpose/Objective: To correct for scatter in cone-beam CT images and to calibrate the image intensities in Hounsfield units (HU) using fast and accurate simulation of x-ray images of the planning CT.