A novel four-dimensional radiotherapy method for lung cancer: imaging, treatment planning and delivery

Computational framework for the generation of one-dimensional vascular models accounting for uncertainty in networks extracted from medical images

One-dimensional (1D) cardiovascular models offer a non-invasive method to answer medical questions, including predictions of wave-reflection, shear stress, functional flow reserve, vascular resistance and compliance. This model type can predict patient-specific outcomes by solving 1D fluid dynamics equations in geometric networks extracted from medical images. However, the inherent uncertainty in in vivo imaging introduces variability in network size and vessel dimensions, affecting haemodynamic predictions. Understanding the influence of variation in image-derived properties is essential to assess the fidelity of model predictions. Numerous programs exist to render three-dimensional surfaces and construct vessel centrelines. Still, there is no exact way to generate vascular trees from the centrelines while accounting for uncertainty in data. This study introduces an innovative framework employing statistical change point analysis to generate labelled trees that encode vessel dimensions and their associated uncertainty from medical images. To test this framework, we explore the impact of uncertainty in 1D haemodynamic predictions in a systemic and pulmonary arterial network. Simulations explore haemodynamic variations resulting from changes in vessel dimensions and segmentation; the latter is achieved by analysing multiple segmentations of the same images. Results demonstrate the importance of accurately defining vessel radii and lengths when generating high-fidelity patient-specific haemodynamics models. KEY POINTS: This study introduces novel algorithms for generating labelled directed trees from medical images, focusing on accurate junction node placement and radius extraction using change points to provide haemodynamic predictions with uncertainty within expected measurement error. Geometric features, such as vessel dimension (length and radius) and network size, significantly impact pressure and flow predictions in both pulmonary and aortic arterial networks. Standardizing networks to a consistent number of vessels is crucial for meaningful comparisons and decreases haemodynamic uncertainty. Change points are valuable to understanding structural transitions in vascular data, providing an automated and efficient way to detect shifts in vessel characteristics and ensure reliable extraction of representative vessel radii.

[Four-dimensional cone-beam CT reconstruction based on motion-compensated robust principal component analysis]

OBJECTIVE

To propose a motion compensation reconstruction method based on robust principal component analysis (RPCA) to reduce the influence of streak artifacts on accurate estimation of interphase motion deformation fields.

OBJECTIVE

We propose a RPCA motion compensation reconstruction algorithm to improve the estimation of motion deformation fields based on the traditional MC-FDK algorithm. RPCA was used to decompose the cone-beam computed tomography (CBCT) images into low-rank and sparse components, and the motion deformation fields between different phase images were then estimated using Horn and Schunck optical flow method from the low-rank images to reduce the influence of striping artifacts on the accuracy of estimation of interphase motion deformation fields. The performance of the algorithm was evaluated using simulation data and real data. The simulation phantom data was obtained by back-projection of 4D-CT images acquired from Philips 16-slice spiral CT using MATLAB software programming according to the scanning geometry of Varian Edge accelerator. The real patient data were obtained using the Elekta Synergy system of CBCT scanning system with half-fan mode CB projection data from lung cancer patients.

OBJECTIVE

Compared with images reconstructed using the traditional MC-FDK algorithm, the reconstructed image using the proposed method had clearer tissue boundaries with reduced motion artifact was reduced. The results of phantom data reconstruction showed that compared with the MC- FDK algorithm, the proposed algorithms resulted in improvements of PSNR by 25.4% and SSIM by 7.6%; compared with the FDK algorithm, PSNR was improved by 37.9% and SSIM by 17.6%.

OBJECTIVE

The proposed algorithm can achieve accurate estimation of inter-phase motion deformation fields and improve the quality of the reconstructed CBCT images.

Does Motion Assessment With 4-Dimensional Computed Tomographic Imaging for Non-Small Cell Lung Cancer Radiotherapy Improve Target Volume Coverage?

Introduction:

Modern radiotherapy with 4-dimensional computed tomographic (4D-CT) image acquisition for non–small cell lung cancer (NSCLC) captures respiratory-mediated tumor motion to provide more accurate target delineation. This study compares conventional 3-dimensional (3D) conformal radiotherapy (3DCRT) plans generated with standard helical free-breathing CT (FBCT) with plans generated on 4D-CT contoured volumes to determine whether target volume coverage is affected.

Materials and methods:

Fifteen patients with stage I to IV NSCLC were enrolled in the study. Free-breathing CT and 4D-CT data sets were acquired at the same simulation session and with the same immobilization. Gross tumor volume (GTV) for primary and/or nodal disease was contoured on FBCT (GTV_3D). The 3DCRT plans were obtained, and the patients were treated according to our institution’s standard protocol using FBCT imaging. Gross tumor volume was contoured on 4D-CT for primary and/or nodal disease on all 10 respiratory phases and merged to create internal gross tumor volume (IGTV)_4D. Clinical target volume margin was 5 mm in both plans, whereas planning tumor volume (PTV) expansion was 1 cm axially and 1.5 cm superior/inferior for FBCT-based plans to incorporate setup errors and an estimate of respiratory-mediated tumor motion vs 8 mm isotropic margin for setup error only in all 4D-CT plans. The 3DCRT plans generated from the FBCT scan were copied on the 4D-CT data set with the same beam parameters. GTV_3D, IGTV_4D, PTV, and dose volume histogram from both data sets were analyzed and compared. Dice coefficient evaluated PTV similarity between FBCT and 4D-CT data sets.

Results:

In total, 14 of the 15 patients were analyzed. One patient was excluded as there was no measurable GTV. Mean GTV_3D was 115.3 cm³ and mean IGTV_4D was 152.5 cm³ (P = .001). Mean PTV_3D was 530.0 cm³ and PTV_4D was 499.8 cm³ (P = .40). Both gross primary and nodal disease analyzed separately were larger on 4D compared with FBCT. D95 (95% isodose line) covered 98% of PTV_3D and 88% of PTV_4D (P = .003). Mean dice coefficient of PTV_3D and PTV_4D was 84%. Mean lung V20 was 24.0% for the 3D-based plans and 22.7% for the 4D-based plans (P = .057). Mean heart V40 was 12.1% for the 3D-based plans and 12.7% for the 4D-based plans (P = .53). Mean spinal cord D_max was 2517 and 2435 cGy for 3D-based and 4D-based plans, respectively (P = .019). Mean esophageal dose was 1580 and 1435 cGy for 3D and 4D plans, respectively (P = .13).

Conclusions:

IGTV_4D was significantly larger than GTV_3D for both primary and nodal disease combined or separately. Mean PTV_3D was larger than PTV_4D, but the difference was not statistically significant. The PTV_4D coverage with 95% isodose line was inferior, indicating the importance of incorporating the true size and shape of the target volume. Relatively less dose was delivered to spinal cord and esophagus with plans based on 4D data set. Dice coefficient analysis for degree of similarity revealed that 16% of PTVs from both data sets did not overlap, indicating different anatomical positions of the PTV due to tumor/nodal motion during a respiratory cycle. All patients with lung cancer planned for radical radiotherapy should have 4D-CT simulation to ensure accurate coverage of the target volumes.

Clinical and research priorities for combined modality therapy in stage III NSCLC

Stage III NSCLC comprises of a heterogeneous group of patients with regard to stage, extent of disease and prognosis. Disease presentation can often be complex in clinical practice, and the challenges are not well-defined in clinical trials or practice guidelines. In order to improve the therapeutic ratio of treatment, one needs to either increase the benefit of treatment or reduce toxicities, or both. Priorities in radiation therapy include dose escalation, imaging advances, improved target delineation, reduced planning target volume margin, intensity-modulated radiotherapy, image-guided radiotherapy, motion management and the use of cytoprotectants. Priorities in systemic therapy include the incorporation of molecularly targeted and immune-modulatory agents.

Visualisation of respiratory tumour motion and co-moving isodose lines in the context of respiratory gating, IMRT and flattening-filter-free beams

Respiratory motion during percutaneous radiotherapy can be considered based on respiration-correlated computed tomography (4DCT). However, most treatment planning systems perform the dose calculation based on a single primary CT data set, even though cine mode displays may allow for a visualisation of the complete breathing cycle. This might create the mistaken impression that the dose distribution were independent of tumour motion. We present a movie visualisation technique with the aim to direct attention to the fact that the dose distribution migrates to some degree with the tumour and discuss consequences for gated treatment, IMRT plans and flattening-filter-free beams. This is a feasibility test for a visualisation of tumour and isodose motion. Ten respiratory phases are distinguished on the CT, and the dose distribution from a stationary IMRT plan is calculated on each phase, to be integrated into a movie of tumour and dose motion during breathing. For one example patient out of the sample of five lesions, the plan is compared with a gated treatment plan with respect to tumour coverage and lung sparing. The interplay-effect for small segments in the IMRT plan is estimated. While the high dose rate, together with the cone-shaped beam profile, makes the use of flattening-filter-free beams more problematic for conformal and IMRT treatment, it can be the option of choice if gated treatment is preferred. The different effects of respiratory motion, dose build-up and beam properties (segments and flatness) for gated vs. un-gated treatment can best be considered if planning is performed on the full 4DCT data set, which may be an incentive for future developments of treatment planning systems.

The radiobiological P(+) index for pretreatment plan assessment with emphasis on four-dimensional radiotherapy modalities

PURPOSE

Radiation treatment modalities will continue to emerge that promise better clinical outcomes albeit technologically challenging to implement. An important question facing the radiotherapy community then is the need to justify the added technological effort for the clinical return. Mobile tumor radiotherapy is a typical example, where 4D tumor tracking radiotherapy (4DTRT) has been proposed over the simpler conventional modality for better results. The modality choice per patient can depend on a wide variety of factors. In this work, we studied the complication-free tumor control probability (P(+)) index, which combines the physical complexity of the treatment plan with the radiobiological characteristics of the clinical case at hand and therefore found to be useful in evaluating different treatment techniques and estimating the expected clinical effectiveness of different radiation modalities.

METHODS

4DCT volumes of 18 previously treated lung cancer patients with tumor motion and size ranging from 2 mm to 15 mm and from 4 cc to 462 cc, respectively, were used. For each patient, 4D treatment plans were generated to extract the 4D dose distributions, which were subsequently used with clinically derived radiobiological parameters to compute the P(+) index per modality.

RESULTS

The authors observed, on average, a statistically significant increase in P(+) of 3.4% ± 3.8% (p < 0.003) in favor of 4DTRT. There was high variability among the patients with a <0.5% up to 13.4% improvement in P(+).

CONCLUSIONS

The observed variability in the improvement of the clinical effectiveness suggests that the relative benefit of tracking should be evaluated on a per patient basis. Most importantly, this variability could be effectively captured in the computed P(+). The index can thus be useful to discriminate and hence point out the need for a complex modality like 4DTRT over another. Besides tumor mobility, a wide range of other factors, e.g., size, location, fractionation, etc., can affect the relative benefits. Application of the P(+) objective is a simple and effective way to combine these factors in the evaluation of a treatment plan.

Thoracic target volume delineation using various maximum-intensity projection computed tomography image sets for radiotherapy treatment planning

PURPOSE

Four-dimensional computed tomography (4D-CT) is commonly used to account for respiratory motion of target volumes in radiotherapy to the thorax. From the 4D-CT acquisition, a maximum-intensity projection (MIP) image set can be created and used to help define the tumor motion envelope or the internal gross tumor volume (iGTV). The purpose of this study was to quantify the differences in automatically contoured target volumes for usage in the delivery of stereotactic body radiation therapy using MIP data sets generated from one of the four methods: (1) 4D-CT phase-binned (PB) based on retrospective phase calculations, (2) 4D-CT phase-corrected phase-binned (PC-PB) based on motion extrema, (3) 4D-CT amplitude-binned (AB), and (4) cine CT built from all available images.

METHODS

MIP image data sets using each of the four methods were generated for a cohort of 28 patients who had prior thoracic 4D-CT scans that exhibited lung tumor motion of at least 1 cm. Each MIP image set was automatically contoured on commercial radiation treatment planning system. Margins were added to the iGTV to observe differences in the final simulated planning target volumes (PTVs).

RESULTS

For all patients, the iGTV measured on the MIP generated from the entire cine CT data set (iGTVcine) was the largest. Expressed as a percentage of iGTVcine, 4D-CT iGTV (all sorting methods) ranged from 83.8% to 99.1%, representing differences in the absolute volume ranging from 0.02 to 4.20 cm3; the largest average and range of 4D-CT iGTV measurements was from the PC-PB data set. Expressed as a percentage of PTVcine (expansions applied to iGTVeine), the 4D-CT PTV ranged from 87.6% to 99.6%, representing differences in the absolute volume ranging from 0.08 to 7.42 cm3. Regions of the measured respiratory waveform corresponding to a rapid change of phase or amplitude showed an increased susceptibility to the selection of identical images for adjacent bins. Duplicate image selection was most common in the AB implementation, followed by the PC-PB method. The authors also found that the image associated with the minimum amplitude measurement did not always correlate with the image that showed maximum tumor motion extent.

CONCLUSIONS

The authors identified cases in which the MIP generated from a 4D-CT sorting process under-represented the iGTV by more than 10% or up to 4.2 cm3 when compared to the iGTVcine. They suggest utilization of a MIP generated from the full cine CT data set to ensure maximum inclusive tumor extent.

A radiobiological analysis of the effect of 3D versus 4D image-based planning in lung cancer radiotherapy

Dose distributions generated on a static anatomy may differ significantly from those delivered to temporally varying anatomy such as for abdominal and thoracic tumors, due largely in part to the unavoidable organ motion and deformation effects stemming from respiration. In this work, the degree of such variation for three treatment techniques, namely static conventional, gating and target tracking radiotherapy, was investigated. The actual delivered dose was approximated by planning all the phases of a 4DCT image set. Data from six (n = 6) previously treated lung cancer patients were used for this study with tumor motion ranging from 2 to 10 mm. Complete radiobiological analyses were performed to assess the clinical significance of the observed discrepancies between the 3D and 4DCT image-based dose distributions. Using the complication-free tumor control probability (P+) objective, we observed small differences in P+ between the 3D and 4DCT image-based plans (<2.0% difference on average) for the gating and static conventional regimens and higher differences in P+ (4.0% on average) for the tracking regimen. Furthermore, we observed, as a general trend, that the 3D plan underestimated the P+ values. While it is not possible to draw any general conclusions from a small patient cohort, our results suggest that there exists a patient population in which 4D planning does not provide any additional benefits beyond that afforded by 3D planning for static conventional or gated radiotherapy. This statement is consistent with previous studies based on physical dosimetric evaluations only. The higher differences observed with the tracking technique suggest that individual patient plans should be evaluated on a case-by-case basis to assess if 3D or 4D imaging is appropriate for the tracking technique.

Delivery of four-dimensional radiotherapy with TrackBeam for moving target using an AccuKnife dual-layer MLC: dynamic phantoms study

Respiratory motion has been considered a clinical challenge for lung tumor treatments due to target motion. In this study, we aimed to perform an experimental evaluation based on dynamic phantoms using MLC‐based beam tracking. TrackBeam, a prototype real‐time beam tracking system, has been assembled and evaluated in our clinic. TrackBeam includes an orthogonal dual‐layer micro multileaf collimator (DmMLC), an on‐board mega‐voltage (MV) portal imaging device, and an image processing workstation. With a fiducial marker implanted in a moving target, the onboard imaging device can capture the motion. The TrackBeam workstation processes the online MV fluence and detects and predicts tumor motion. The DmMLC system then dynamically repositions each leaf to form new beam apertures based on the movement of the fiducial marker. In this study, a dynamic phantom was used for the measurements. Three delivery patterns were evaluated for dosimetric verification based on radiographic films: no‐motion lung‐tumor (NMLT), three‐dimensional conformal radiotherapy (3DCRT), and four‐dimensional tracking radiotherapy (4DTRT). The displacement between the DmMLC dynamic beam isocenter and the fiducial marker was in the range of 0.5 mm to 1.5 mm. With radiographic film analysis, the planar dose histogram difference between 3DCRT and NLMT was 48.6% and 38.0% with dose difference tolerances of 10% and 20%, respectively. The planar dose histogram difference between 4DTRT and NLMT was 15.2% and 4.0%, respectively. Based on dose volume histogram analysis, 4DTRT reduces the mean dose for the surrounding tissue from 35.4 Gy to 19.5 Gy, reduces the relative volume of the total lung from 28% to 18% at V20, and reduces the amount of dose from 35.2 Gy to 15.0 Gy at D20. The experimental results show that MLC‐based real‐time beam tracking delivery provides a potential solution to respiratory motion control. Beam tracking delivers a highly conformal dose to a moving target, while sparing surrounding normal tissue.

PACS number: 87.55.de, 87.55.ne, 87.56.nk

Influence of increased target dose inhomogeneity on margins for breathing motion compensation in conformal stereotactic body radiotherapy

Background

Breathing motion should be considered for stereotactic body radiotherapy (SBRT) of lung tumors. Four-dimensional computer tomography (4D-CT) offers detailed information of tumor motion. The aim of this work is to evaluate the influence of inhomogeneous dose distributions in the presence of breathing induced target motion and to calculate margins for motion compensation.

Methods

Based on 4D-CT examinations, the probability density function of pulmonary tumors was generated for ten patients. The time-accumulated dose to the tumor was calculated using one-dimensional (1D) convolution simulations of a 'static' dose distribution and target probability density function (PDF). In analogy to stereotactic body radiotherapy (SBRT), different degrees of dose inhomogeneity were allowed in the target volume: minimum doses of 100% were prescribed to the edge of the target and maximum doses varied between 102% (P102) and 150% (P150). The dose loss due to breathing motion was quantified and margins were added until this loss was completely compensated.

Results

With the time-weighted mean tumor position as the isocentre, a close correlation with a quadratic relationship between the standard deviation of the PDF and the margin size was observed. Increased dose inhomogeneity in the target volume required smaller margins for motion compensation: margins of 2.5 mm, 2.4 mm and 1.3 mm were sufficient for compensation of 11.5 mm motion range and standard deviation of 3.9 mm in P105, P125 and P150, respectively. This effect of smaller margins for increased dose inhomogeneity was observed for all patients. Optimal sparing of the organ-at-risk surrounding the target was achieved for dose prescriptions P105 to P118. The internal target volume concept over-compensated breathing motion with higher than planned doses to the target and increased doses to the surrounding normal tissue.

Conclusion

Treatment planning with inhomogeneous dose distributions in the target volume required smaller margins for compensation of breathing induced target motion with the consequence of lower doses to the surrounding organs-at-risk.

Development of a novel post-processing treatment planning platform for 4D radiotherapy

The aim of this study is to develop an Automatic Post-processing Tool for four-dimensional (4D) treatment planning (APT4D) that enables the user to perform some necessary procedures related to 4D treatment planning, such as automated image registration, automatic propagation of regions of interest, and dose distribution transformation. Demons-based deformable registrations were performed to map the moving phase images (such as the end-inhalation phase or 0% phase) to the reference phase (typically the end-exhalation fixed phase or 50% phase). Contours were automatically propagated into the moving phase using the image registration results. The dose distribution of each moving phase was transformed to the fixed phase and subsequently was summed as an average with equal weighting factor. To validate the application of APT4D utility, the 4D computed tomography (CT) images of a lung cancer patient and an abdominal cancer patient were acquired and resorted into ten respiratory phases. 4D plans based on the 4D CT images were developed. The correlation coefficient ranged from 0.992 to 0.999 for the re-sampled deformed moving phase image against the fixed phase image for the lung patient plan and from 0.977 to 0.999 for the abdominal patient plan. For all the organs, the match indices between the manual contours and automatic contour propagation results were around 0.92 to 0.95. The 4D composite dose-volume histogram showed dosimetric reductions for liver and kidneys in the high dose region. The APT4D adds automation, efficiency, and functionality, while integrating the whole process of 4D treatment planning.

Cine computed tomography without respiratory surrogate in planning stereotactic radiotherapy for non-small-cell lung cancer

PURPOSE

To determine whether cine computed tomography (CT) can serve as an alternative to four-dimensional (4D)-CT by providing tumor motion information and producing equivalent target volumes when used to contour in radiotherapy planning without a respiratory surrogate.

METHODS AND MATERIALS

Cine CT images from a commercial CT scanner were used to form maximum intensity projection and respiratory-averaged CT image sets. These image sets then were used together to define the targets for radiotherapy. Phantoms oscillating under irregular motion were used to assess the differences between contouring using cine CT and 4D-CT. We also retrospectively reviewed the image sets for 26 patients (27 lesions) at our institution who had undergone stereotactic radiotherapy for Stage I non-small-cell lung cancer. The patients were included if the tumor motion was >1 cm. The lesions were first contoured using maximum intensity projection and respiratory-averaged CT image sets processed from cine CT and then with 4D-CT maximum intensity projection and 10-phase image sets. The mean ratios of the volume magnitude were compared with intraobserver variation, the mean centroid shifts were calculated, and the volume overlap was assessed with the normalized Dice similarity coefficient index.

RESULTS

The phantom studies demonstrated that cine CT captured a greater extent of irregular tumor motion than did 4D-CT, producing a larger tumor volume. The patient studies demonstrated that the gross tumor defined using cine CT imaging was similar to, or slightly larger than, that defined using 4D-CT.

CONCLUSION

The results of our study have shown that cine CT is a promising alternative to 4D-CT for stereotactic radiotherapy planning.

Assessment of intrafraction mediastinal and hilar lymph node movement and comparison to lung tumor motion using four-dimensional CT

PURPOSE

To quantify the amount of free-breathing motion measured using Four-dimensional (4D) CT scans of mediastinal and hilar lymph nodes and to compare this motion to the primary lung tumor motion.

METHODS AND MATERIALS

Twenty patients with primary lung cancer, radiographically positive lymph nodes, and prior 4D CT scans were retrospectively analyzed. The 4D CT data sets were divided into four respiratory phases, and the primary tumor and radiographically positive nodes were contoured. Geometric and volumetric analysis was performed to analyze the motion of the primary tumors and the lymph nodes.

RESULTS

The mean lymph node motion was 2.6 mm in the mediolateral direction, 2.5 mm in the anterior-posterior direction, and 5.2 mm in the cranial-caudal direction with a maximum of 14.4 mm. All lymph nodes were found to move inferiorly during inspiration, with 12.5% of nodes moving more than 1 cm. Lymph nodes located below the carina showed significantly more motion than those above the carina (p = 0.01). In comparing the primary tumor motion to the lymph node motion, no correlation was identified.

CONCLUSIONS

Four-dimensional CT scans can be used to measure the motion of the primary lung tumor and pathologic lymph nodes encountered during the respiratory cycle. Both the primary lung tumor and the lymph node must to be examined to assess their individual degree of motion. This study demonstrates the need for individualized plans to assess the heterogeneous motion encountered in both primary lung tumors and among lymph node stations.

Defining internal target volume (ITV) for hepatocellular carcinoma using four-dimensional CT

BACKGROUND AND PURPOSE

To define individualized internal target volume (ITV) for hepatocellular carcinoma using four-dimensional computed tomography (4DCT).

MATERIALS AND METHODS

Gross tumor volumes (GTVs) and clinical target volumes (CTVs) were contoured on all 10 respiratory phases of 4DCT scans in 10 patients with hepatocellular carcinoma. The 3D and 4D treatment plans were performed for each patient using two different planning target volumes (PTVs): (1) PTV(3D) was derived from a single CTV plus conventional margins; (2) PTV(4D) was derived from ITV(4D), which encompassed all 10 CTVs plus setup margins (SMs). The volumes of PTVs and dose distribution were compared between the two plans.

RESULTS

The average PTV volume of the 4D plans (328.4+/-152.2cm(3)) was less than 3D plans (407.0+/-165.6cm(3)). The 4D plans spared more surrounding normal tissues than 3D plans, especially normal liver. Compared with 3D plans, the mean dose to normal liver (MDTNL) decreased from 22.7 to 20.3Gy. Without increasing the normal tissue complication probability (NTCP), the 4D plans allowed for increasing the calculated dose from 50.4+/-1.3 to 54.2+/-2.6Gy, an average increase of 7.5% (range 4.0-16.0%).

CONCLUSIONS

The conventional 3D plans can result in geometric miss and include excess normal tissues. The 4DCT-based plans can reduce the target volumes to spare more normal tissues and allow dose escalation compared with 3D plans.

Automatic segmentation of phase-correlated CT scans through nonrigid image registration using geometrically regularized free-form deformation

Conventional radiotherapy is planned using free-breathing computed tomography (CT), ignoring the motion and deformation of the anatomy from respiration. New breath-hold-synchronized, gated, and four-dimensional (4D) CT acquisition strategies are enabling radiotherapy planning utilizing a set of CT scans belonging to different phases of the breathing cycle. Such 4D treatment planning relies on the availability of tumor and organ contours in all phases. The current practice of manual segmentation is impractical for 4D CT, because it is time consuming and tedious. A viable solution is registration-based segmentation, through which contours provided by an expert for a particular phase are propagated to all other phases while accounting for phase-to-phase motion and anatomical deformation. Deformable image registration is central to this task, and a free-form deformation-based nonrigid image registration algorithm will be presented. Compared with the original algorithm, this version uses novel, computationally simpler geometric constraints to preserve the topology of the dense control-point grid used to represent free-form deformation and prevent tissue fold-over. Using mean squared difference as an image similarity criterion, the inhale phase is registered to the exhale phase of lung CT scans of five patients and of characteristically low-contrast abdominal CT scans of four patients. In addition, using expert contours for the inhale phase, the corresponding contours were automatically generated for the exhale phase. The accuracy of the segmentation (and hence deformable image registration) was judged by comparing automatically segmented contours with expert contours traced directly in the exhale phase scan using three metrics: volume overlap index, root mean square distance, and Hausdorff distance. The accuracy of the segmentation (in terms of radial distance mismatch) was approximately 2 mm in the thorax and 3 mm in the abdomen, which compares favorably to the accuracies reported elsewhere. Unlike most prior work, segmentation of the tumor is also presented. The clinical implementation of 4D treatment planning is critically dependent on automatic segmentation, for which is offered one of the most accurate algorithms yet presented.

Investigation of the dosimetric effect of respiratory motion using four-dimensional weighted radiotherapy

We have developed a four-dimensional weighted radiotherapy (4DW-RT) technique. This method involves designing the motion of the linear accelerator beam to coincide with the tumour motion determined from 4D-CT imaging while including a weighting factor to account for irregular motion and limitations of the delivery system. Experiments were conducted with a moving phantom to assess limitations of the delivery system when applying this method. Although the multi-leaf collimator motion remains within the tolerance of the linear accelerator, the extent of motion was less than 1 mm larger than the designed one, and there was a net system latency of approximately 0.2 s. The dose distributions were measured and simulated using different weighting factors and motion scenarios. The breathing characteristics (period, extent of motion, drift and standard deviations) of 32 patients were evaluated using the Varian RPM system. Breathing variability was assessed by plotting the average breathing motion as a function of the breathing phase. Simulations were carried out to determine the optimal weighting factor based on typical patient breathing characteristics. These results establish that the 4DW-RT method demonstrates potential for dose escalation without increasing exposure to healthy tissue.

Analysis of organ motion effects on the effective fluences for liver IMRT

An effective fluence concept was employed to make forward dose calculations to investigate the effects of a distorted fluence map on dose plans. Fluence changes caused by organ motion were calculated using Chui's algorithm (2003 Med. Phys. 30 1736). In two test cases with various fluence maps, the effects of motion were simulated using a maximal displacement from 5 mm to 25 mm; 108 fluence maps that were calculated from 16 IMRT plans for eight liver cancer patients were analyzed and compared with and without gating. Fluoroscopic measurements were made of a moving diaphragm in this study. Fluence changes associated with superior-inferior organ motion, perpendicular to the moving MLC, were also examined. The effects of motion on the fluence maps were evaluated from both the fluence differences between static and motion and the chi function. The maximum displacements of the organs in all of these cases were analyzed and correlated with the change in fluence generated from the liver IMRT plans. The dosimetric effects on the target coverage were evaluated for each plan. The results indicate that, for the same fluence map, the mean fluence intensity error or the percentage of the fluence points that have an unacceptable error is linearly related to the extent of motion. For different fluence maps, the degree to which the fluence is distorted by motion is strongly related to the product of the motion extent and the fluence gradient in the direction of diaphragm motion. For eight liver patients and 16 IMRT plans in this work (with gated technique, motion extent from 0.5 cm to 1.0 cm; without gated technique, motion extent from 0.9 cm to 1.8 cm), the fluence modulations are mild, such that the respiratory motion of each patient did not strongly affect the CTV coverage. The mean dose error is 1.5% for free motion (0.9-1.8 cm) and is around 1% for gated motion (0.5-1 cm).

Imaging tumour motion for radiotherapy planning using MRI

Novel technology has made dynamic magnetic resonance imaging (MRI) of lung motion and lung tumour mobility during continuous respiration feasible. This might be beneficial for planning of radiotherapy of lung tumours, especially when using high precision techniques. This paper describes the recent developments to analyze and visualize pulmonary nodules during continuous respiration using MRI. Besides recent dynamic two-dimensional approaches to quantify motion of pulmonary nodules during respiration novel three-dimensional techniques are presented. Beyond good correlation to pulmonary function tests MRI also provides regional information about differences between tumour-bearing and non-tumour bearing lung and the restrictive effects of radiotherapy as well as the compensation by the contralateral lung.