Accuracy of finite element model-based multi-organ deformable image registration

Estimation of the delivered patient dose in lung IMRT treatment based on deformable registration of 4D-CT data and Monte Carlo simulations

The purpose of this study is to accurately estimate the difference between the planned and the delivered dose due to respiratory motion and free breathing helical CT artefacts for lung IMRT treatments, and to estimate the impact of this difference on clinical outcome. Six patients with representative tumour motion, size and position were selected for this retrospective study. For each patient, we had acquired both a free breathing helical CT and a ten-phase 4D-CT scan. A commercial treatment planning system was used to create four IMRT plans for each patient. The first two plans were based on the GTV as contoured on the free breathing helical CT set, with a GTV to PTV expansion of 1.5 cm and 2.0 cm, respectively. The third plan was based on the ITV, a composite volume formed by the union of the CTV volumes contoured on free breathing helical CT, end-of-inhale (EOI) and end-of-exhale (EOE) 4D-CT. The fourth plan was based on GTV contoured on the EOE 4D-CT. The prescribed dose was 60 Gy for all four plans. Fluence maps and beam setup parameters of the IMRT plans were used by the Monte Carlo dose calculation engine MCSIM for absolute dose calculation on both the free breathing CT and 4D-CT data. CT deformable registration between the breathing phases was performed to estimate the motion trajectory for both the tumour and healthy tissue. Then, a composite dose distribution over the whole breathing cycle was calculated as a final estimate of the delivered dose. EUD values were computed on the basis of the composite dose for all four plans. For the patient with the largest motion effect, the difference in the EUD of CTV between the planed and the delivered doses was 33, 11, 1 and 0 Gy for the first, second, third and fourth plan, respectively. The number of breathing phases required for accurate dose prediction was also investigated. With the advent of 4D-CT, deformable registration and Monte Carlo simulations, it is feasible to perform an accurate calculation of the delivered dose, and compare our delivered dose with doses estimated using prior techniques.

Feasibility of a novel deformable image registration technique to facilitate classification, targeting, and monitoring of tumor and normal tissue

PURPOSE

To investigate the feasibility of a biomechanical-based deformable image registration technique for the integration of multimodality imaging, image guided treatment, and response monitoring.

METHODS AND MATERIALS

A multiorgan deformable image registration technique based on finite element modeling (FEM) and surface projection alignment of selected regions of interest with biomechanical material and interface models has been developed. FEM also provides an inherent method for direct tracking specified regions through treatment and follow-up.

RESULTS

The technique was demonstrated on 5 liver cancer patients. Differences of up to 1 cm of motion were seen between the diaphragm and the tumor center of mass after deformable image registration of exhale and inhale CT scans. Spatial differences of 5 mm or more were observed for up to 86% of the surface of the defined tumor after deformable image registration of the computed tomography (CT) and magnetic resonance images. Up to 6.8 mm of motion was observed for the tumor after deformable image registration of the CT and cone-beam CT scan after rigid registration of the liver. Deformable registration of the CT to the follow-up CT allowed a more accurate assessment of tumor response.

CONCLUSIONS

This biomechanical-based deformable image registration technique incorporates classification, targeting, and monitoring of tumor and normal tissue using one methodology.

A material sensitivity study on the accuracy of deformable organ registration using linear biomechanical modelsa)

Model-based deformable organ registration techniques using the finite element method (FEM) have recently been investigated intensively and applied to image-guided adaptive radiotherapy (IGART). These techniques assume that human organs are linearly elastic material, and their mechanical properties are predetermined. Unfortunately, the accurate measurement of the tissue material properties is challenging and the properties usually vary between patients. A common issue is therefore the achievable accuracy of the calculation due to the limited access to tissue elastic material constants. In this study, we performed a systematic investigation on this subject based on tissue biomechanics and computer simulations to establish the relationships between achievable registration accuracy and tissue mechanical and organ geometrical properties. Primarily we focused on image registration for three organs: rectal wall, bladder wall, and prostate. The tissue anisotropy due to orientation preference in tissue fiber alignment is captured by using an orthotropic or a transversely isotropic elastic model. First we developed biomechanical models for the rectal wall, bladder wall, and prostate using simplified geometries and investigated the effect of varying material parameters on the resulting organ deformation. Then computer models based on patient image data were constructed, and image registrations were performed. The sensitivity of registration errors was studied by perturbating the tissue material properties from their mean values while fixing the boundary conditions. The simulation results demonstrated that registration error for a subvolume increases as its distance from the boundary increases. Also, a variable associated with material stability was found to be a dominant factor in registration accuracy in the context of material uncertainty. For hollow thin organs such as rectal walls and bladder walls, the registration errors are limited. Given 30% in material uncertainty, the registration error is limited to within 1.3 mm. For a solid organ such as the prostate, the registration errors are much larger. Given 30% in material uncertainty, the registration error can reach 4.5 mm. However, the registration error distribution for prostates shows that most of the subvolumes have a much smaller registration error. A deformable organ registration technique that uses FEM is a good candidate in IGART if the mean material parameters are available.

Deformable registration for image-guided radiation therapy

Ths paper examines several applications of deformable registration algorithms in the field of image-guided radiotherapy. The first part focuses on the description of input and output of deformable registration algorithms, with a brief review of conventional and most current methods. The typical applications of deformable registration are then reviewed on the basis of four practical examples. The first two sets of examples deal with the fusion of images obtained from the same patient (inter-fraction registration), with time intervals of several days between each image. The other two examples deal with the fusion of images obtained in immediate sequence (intra-fraction registration); in this case, the focus is the displacement during image acquisition or patient treatment (mainly due to respiratory movement), with time intervals in the order of magnitude of tenths of seconds. Finally, the registration of images of different patients (inter-patient registration) is also discussed. In conclusion, deformable registration has become a fundamental tool for image analysis in radiotherapy. Although extensive validation of the numerous existing methods is required before extending its clinical use, deformable registration is expected to become a standard methodology in the treatment planning systems in the near future.