M Szegedi

4D CT image reconstruction with diffeomorphic motion model.

Four-dimensional (4D) respiratory correlated computed tomography (RCCT) has been widely used for studying organ motion. Most current RCCT imaging algorithms use binning techniques that are susceptible to artifacts and challenge the quantitative analysis of organ motion. In this paper, we develop an algorithm for analyzing organ motion which uses the raw, time-stamped imaging data to reconstruct images while simultaneously estimating deformation in the subjects anatomy. This results in reduction of artifacts and facilitates a reduction in dose to the patient during scanning while providing equivalent or better image quality as compared to RCCT. The framework also incorporates fundamental physical properties of organ motion, such as the conservation of local tissue volume. We demonstrate that this approach is accurate and robust against noise and irregular breathing patterns. We present results for a simulated cone beam CT phantom, as well as a detailed real porcine liver phantom study demonstrating accuracy and robustness of the algorithm. An example of applying this algorithm to real patient image data is also presented.

Evaluation of alignment error due to a speed artifact in stereotactic ultrasound image guidance

Ultrasound (US) image guidance systems used in radiotherapy are typically calibrated for soft tissue applications, thus introducing errors in depth-from-transducer representation when used in media with a different speed of sound propagation (e.g. fat). This error is commonly referred to as the speed artifact. In this study we utilized a standard US phantom to demonstrate the existence of the speed artifact when using a commercial US image guidance system to image through layers of simulated body fat, and we compared the results with calculated/predicted values. A general purpose US phantom (speed of sound (SOS) = 1540 m s(-1)) was imaged on a multi-slice CT scanner at a 0.625 mm slice thickness and 0.5 mm x 0.5 mm axial pixel size. Target-simulating wires inside the phantom were contoured and later transferred to the US guidance system. Layers of various thickness (1-8 cm) of commercially manufactured fat-simulating material (SOS = 1435 m s(-1)) were placed on top of the phantom to study the depth-related alignment error. In order to demonstrate that the speed artifact is not caused by adding additional layers on top of the phantom, we repeated these measurements in an identical setup using commercially manufactured tissue-simulating material (SOS = 1540 m s(-1)) for the top layers. For the fat-simulating material used in this study, we observed the magnitude of the depth-related alignment errors resulting from the speed artifact to be 0.7 mm cm(-1) of fat imaged through. The measured alignment errors caused by the speed artifact agreed with the calculated values within one standard deviation for all of the different thicknesses of fat-simulating material studied here. We demonstrated the depth-related alignment error due to the speed artifact when using US image guidance for radiation treatment alignment and note that the presence of fat causes the target to be aliased to a depth greater than it actually is. For typical US guidance systems in use today, this will lead to delivery of the high dose region at a position slightly posterior to the intended region for a supine patient. When possible, care should be taken to avoid imaging through a thick layer of fat for larger patients in US alignments or, if unavoidable, the spatial inaccuracies introduced by the artifact should be considered by the physician during the formulation of the treatment plan.

A proto-type design of a real-tissue phantom for the validation of deformation algorithms and 4D dose calculations

The purpose of this study is to design a real-tissue phantom for use in the validation of deformation algorithms. A phantom motion controller that runs sinusoidal and non-regular patient-based breathing pattern, via a piston, was applied to porcine liver tissue. It was regulated to simulate movement ranges similar to recorded implanted liver markers from patients. 4D CT was applied to analyze deformation. The suitability of various markers in the liver and the position reproducibility of markers and of reference points were studied. The similarity of marker motion pattern in the liver phantom and in real patients was evaluated. The viability of the phantom over time and its use with electro-magnetic tracking devices were also assessed. High contrast markers, such as carbon markers, implanted in the porcine liver produced less image artifacts on CT and were well visualized compared to metallic ones. The repositionability of markers was within a measurement accuracy of +/-2 mm. Similar anatomical patient motions were reproducible up to elongations of 3 cm for a time period of at least 90 min. The phantom is compatible with electro-magnetic tracking devices and 4D CT. The phantom motion is reproducible and simulates realistic patient motion and deformation. The ability to carry out voxel-based tracking allows for the evaluation of deformation algorithms in a controlled environment with recorded patient traces. The phantom is compatible with all therapy devices clinically encountered in our department.

Individualized margins for prostate patients using a wireless localization and tracking system.

This study investigates the dosimetric benefits of designing patient‐specific margins for prostate cancer patients based on 4D localization and tracking. Ten prostate patients, each implanted with three radiofrequency transponders, were localized and tracked for 40 fractions. “Conventional margin” (CM) planning target volumes (PTV) and PTVs resulting from uniform margins of 5 mm (5M) and 7 mm (7M) were explored. Through retrospective review of each patients tracking data, an individualized margin (IM) design for each patient was determined. IMRT treatment plans with identical constraints were generated for all four margin strategies and compared. The IM plans generally created the smallest PTV volumes. For similar PTV coverage, the IM plans had a lower mean bladder (rectal) dose by an average of 3.9% (2.5%), 8.5% (5.7%) and 16.2 % (9.8%) compared to 5M, 7M and CM plans, respectively. The IM plan had the lowest gEUD value of 23.8 Gy for bladder, compared to 35.1, 28.4 and 25.7, for CM, 7M and 5M, respectively. Likewise, the IM plan had the lowest NTCP value for rectum of 0.04, compared to 0.07, 0.06 and 0.05 for CM, 7M and 5M, respectively. Individualized margins can lead to significantly reduced PTV volumes and critical structure doses, while still ensuring a minimum delivered CTV dose equal to 95% of the prescribed dose. PACS numbers: 87.53.Kn, 87.55.D

Dosimetric evaluation of a Monte Carlo IMRT treatment planning system incorporating the MIMiC

The high dose per fraction delivered to lung lesions in stereotactic body radiation therapy (SBRT) demands high dose calculation and delivery accuracy. The inhomogeneous density in the thoracic region along with the small fields used typically in intensity-modulated radiation therapy (IMRT) treatments poses a challenge in the accuracy of dose calculation. In this study we dosimetrically evaluated a pre-release version of a Monte Carlo planning system (PEREGRINE 1.6b, NOMOS Corp., Cranberry Township, PA), which incorporates the modeling of serial tomotherapy IMRT treatments with the binary multileaf intensity modulating collimator (MIMiC). The aim of this study is to show the validation process of PEREGRINE 1.6b since it was used as a benchmark to investigate the accuracy of doses calculated by a finite size pencil beam (FSPB) algorithm for lung lesions treated on the SBRT dose regime via serial tomotherapy in our previous study. Doses calculated by PEREGRINE were compared against measurements in homogeneous and inhomogeneous materials carried out on a Varian 600C with a 6 MV photon beam. Phantom studies simulating various sized lesions were also carried out to explain some of the large dose discrepancies seen in the dose calculations with small lesions. Doses calculated by PEREGRINE agreed to within 2% in water and up to 3% for measurements in an inhomogeneous phantom containing lung, bone and unit density tissue.

Planning for mARC treatments with the Eclipse treatment planning system

While modulated arc (mARC) capabilities have been available on Siemens linear accelerators for almost two years now, there was, until recently, only one treatment planning system capable of planning these treatments. The Eclipse treatment planning system now offers a module that can plan for mARC treatments. The purpose of this work was to test the module to determine whether it is capable of creating clinically acceptable plans. A total of 23 plans were created for various clinical sites and all plans delivered without anomaly. The average 3%/3 mm gamma pass rate for the plans was 98.0%, with a standard deviation of 1.7%. For a total of 14 plans, an equivalent static gantry IMRT plan was also created to compare delivery time. In all but two cases, the mARC plans delivered significantly faster than the static gantry plan. We have confirmed the successful creation of mARC plans that are deliverable with high fidelity on an ARTISTE linear accelerator, thus demonstrating the successful implementation of the Eclipse mARC module. PACS numbers: 87.55.D‐, 87.55.ne, 87.57.uq,

Four-dimensional tissue deformation reconstruction (4D TDR) validation using a real tissue phantom

Calculation of four‐dimensional (4D) dose distributions requires the remapping of dose calculated on each available binned phase of the 4D CT onto a reference phase for summation. Deformable image registration (DIR) is usually used for this task, but unfortunately almost always considers only endpoints rather than the whole motion path. A new algorithm, 4D tissue deformation reconstruction (4D TDR), that uses either CT projection data or all available 4D CT images to reconstruct 4D motion data, was developed. The purpose of this work is to verify the accuracy of the fit of this new algorithm using a realistic tissue phantom. A previously described fresh tissue phantom with implanted electromagnetic tracking (EMT) fiducials was used for this experiment. The phantom was animated using a sinusoidal and a real patient‐breathing signal. Four‐dimensional computer tomography (4D CT) and EMT tracking were performed. Deformation reconstruction was conducted using the 4D TDR and a modified 4D TDR which takes real tissue hysteresis (4D TDRHysteresis) into account. Deformation estimation results were compared to the EMT and 4D CT coordinate measurements. To eliminate the possibility of the high contrast markers driving the 4D TDR, a comparison was made using the original 4D CT data and data in which the fiducials were electronically masked. For the sinusoidal animation, the average deviation of the 4D TDR compared to the manually determined coordinates from 4D CT data was 1.9 mm, albeit with as large as 4.5 mm deviation. The 4D TDR calculation traces matched 95% of the EMT trace within 2.8 mm. The motion hysteresis generated by real tissue is not properly projected other than at endpoints of motion. Sinusoidal animation resulted in 95% of EMT measured locations to be within less than 1.2 mm of the measured 4D CT motion path, enabling accurate motion characterization of the tissue hysteresis. The 4D TDRHysteresis calculation traces accounted well for the hysteresis and matched 95% of the EMT trace within 1.6 mm. An irregular (in amplitude and frequency) recorded patient trace applied to the same tissue resulted in 95% of the EMT trace points within less than 4.5 mm when compared to both the 4D CT and 4D TDRHysteresis motion paths. The average deviation of 4D TDRHysteresis compared to 4D CT datasets was 0.9 mm under regular sinusoidal and 1.0 mm under irregular patient trace animation. The EMT trace data fit to the 4D TDRHysteresis was within 1.6 mm for sinusoidal and 4.5 mm for patient trace animation. While various algorithms have been validated for end‐to‐end accuracy, one can only be fully confident in the performance of a predictive algorithm if one looks at data along the full motion path. The 4D TDR, calculating the whole motion path rather than only phase‐ or endpoints, allows us to fully characterize the accuracy of a predictive algorithm, minimizing assumptions. This algorithm went one step further by allowing for the inclusion of tissue hysteresis effects, a real‐world effect that is neglected when endpoint‐only validation is performed. Our results show that the 4D TDRHysteresis correctly models the deformation at the endpoints and any intermediate points along the motion path. PACS numbers: 87.55.km, 87.55.Qr, 87.57.nf, 87.85.Tu

SU-FF-J-128: Characterization of Liver Motion Based On Implanted Markers

Purpose: The purpose of this work is to evaluate the use of implanted markers to characterize motion and deformation of livertumors and tissue.Method and Materials: We retrospectively analyzed 4DCT data sets of livercancer patients treated with stereotactic body radiation therapy at our center. All of these patients had either gold‐ or carbon‐markers (Civco medical solutions, IA) implanted in the surrounding tissue of the lesion. Using a Varian RPM system, ten phase‐CT‐image sets were acquired for treatment planning and prior to every treatment. Two references points, one in bony anatomy and one external were used for analysis. The movements of markers with reference to each other and the reference‐locations for each phase as well as the shift of the geometric center of the implanted markers were studied. The intra session and inter session motion range and geometric deformation was investigated. Results: The average motion‐range of the three markers for these patients differed in location, direction and magnitude for each day as well as from day to day. All motion vectors of these markers show a predominance in SI motion, but deviated for some degrees from each other. The spread of motion vectors additionally indicates a shift of the centroid of the virtual marker triangle for each phase, showing the deformation of tissue and the change in position in LR and AP direction. The greatest difference in motion‐range between the treatment days for one marker was 10.4 mm for one marker, which was located centrally the liver.Conclusion: Motion of liver and livertumors can be visualized with implanted markers and is similar to previous publications, comparing organ outlines/edges for motion deduction. Markers additionally provide information on different rates of motion at point locations and deformation within the region of interest in the organ, helping deformation algorithm validation.

Dosimetric impact of the 160 MLC on head and neck IMRT treatments.

The purpose of this work is to investigate if the change in plan quality with the finer leaf resolution and lower leakage of the 160 MLC would be dosimetrically significant for head and neck intensity‐modulated radiation therapy (IMRT) treatment plans. The 160 MLC consisting of 80 leaves of 0.5 cm on each bank, a leaf span of 20 cm, and leakage of less than 0.37% without additional backup jaws was compared against the 120 Millennium MLC with 60 leaves of 0.5 and 1.0 cm, a leaf span of 14.5 cm, and leakage of 2.0%. CT image sets of 16 patients previously treated for stage III and IV head and neck carcinomas were replanned on Prowess 5.0 and Eclipse 11.0 using the 160 MLC and the 120 MLC. IMRT constraints for both sets of 6 MV plans were identical and based on RTOG 0522. Dose‐volume histograms (DVHs), minimum dose, mean dose, maximum dose, and dose to 1 cc to the organ at risks (OAR) and the planning target volume, as recommended by QUANTEC 2010, were compared. Both collimators were able to achieve the target dose to the PTVs. The dose to the organs at risk (brainstem, spinal cord, parotids, and larynx) were 1%–12% (i.e., 0.5–8 Gy for a 70 Gy prescription) lower with the 160 MLC compared to the 120 MLC, depending on the proximity of the organ to the target. The large field HN plans generated with the 160 MLC were dosimetrically advantageous for critical structures, especially those located further away from the central axis, without compromising the target volume. PACS number: 87.55 D‐

SU‐E‐U‐08: Presentation of a New Intrafractional Prostate Monitoring Method with Ultrasound Image Guidance During Radiotherapy Treatment

Purpose: Clinical presentation of a new intrafractional ultrasound image guidance (USIG) prostate monitoring/tracking method. Methods: Clarity ultrasound system has recently released a new feature to monitor/track intrafractional prostate motion using a trans‐perineal image acquisition position. After initial localization and image guidance correction, the ultrasound probe remains fixed in the sagittal‐plane, trans‐perineal imaging position via a couch‐mounted arm‐support system. The software then enters a live monitoring mode, where the scanning ultrasound probe continuously acquires fanned sagittal images for tracking of during‐treatment prostate position. Compared to other tracking technologies (e.g. RF implanted beacon tracking, or real‐time fluoro monitoring), the ultrasound system has the unique advantages of live 3D image display, without invasive procedure or imaging radiation dose. We present our initial experience using this monitoring feature for 8 patients and 218 treatment sessions. Results: The monitoring software functioned as expected during treatment with consistent reporting of prostate deviation from isocenter location, and notification of instances where positional error exceeded our tolerance of 3 mm for 5 second. Comparisons of prostate tracking data streams obtained from Clarity USIG with data streams obtained in our clinic from Calypso RF tracking of a similar patient population show similarity of recorded motion information for both methods. Average during‐treatment prostate motion for both tracking methods was on the order of 1–2 mm in 3 cardinal directions, with roughly 5% of instances where motion exceeded 3 mm, 5 second tolerance level. When Monitoring reported that tolerance was exceeded we confirmed this by performing redundant, static imaging Alignment, with confirmation that target position had, indeed, changed as reported. Conclusion: The newly released trans‐perineal ultrasound Monitoring approach was confirmed to function well clinically, and to report accurate intrafractional monitoring data. Acquired motion streams were consistent with our Calypso monitoring experience obtained on a similar patient population in our own clinic.