Navigator channel adaptation to reconstruct three dimensional heart volumes from two dimensional radiotherapy planning data
© Ng et al; licensee BioMed Central Ltd. 2012
Received: 26 May 2011
Accepted: 18 January 2012
Published: 18 January 2012
Biologically-based models that utilize 3D radiation dosimetry data to estimate the risk of late cardiac effects could have significant utility for planning radiotherapy in young patients. A major challenge arises from having only 2D treatment planning data for patients with long-term follow-up. In this study, we evaluate the accuracy of an advanced deformable image registration (DIR) and navigator channels (NC) adaptation technique to reconstruct 3D heart volumes from 2D radiotherapy planning images for Hodgkin's Lymphoma (HL) patients.
Planning CT images were obtained for 50 HL patients who underwent mediastinal radiotherapy. Twelve image sets (6 male, 6 female) were used to construct a male and a female population heart model, which was registered to 23 HL "Reference" patients' CT images using a DIR algorithm, MORFEUS. This generated a series of population-to-Reference patient specific 3D deformation maps. The technique was independently tested on 15 additional "Test" patients by reconstructing their 3D heart volumes using 2D digitally reconstructed radiographs (DRR). The technique involved: 1) identifying a matching Reference patient for each Test patient using thorax measurements, 2) placement of six NCs on matching Reference and Test patients' DRRs to capture differences in significant heart curvatures, 3) adapting the population-to-Reference patient-specific deformation maps to generate population-to-Test patient-specific deformation maps using linear and bilinear interpolation methods, 4) applying population-to-Test patient specific deformation to the population model to reconstruct Test-patient specific 3D heart models. The percentage volume overlap between the NC-adapted reconstruction and actual Test patient's true heart volume was calculated using the Dice coefficient.
The average Dice coefficient expressed as a percentage between the NC-adapted and actual Test model was 89.4 ± 2.8%. The modified NC adaptation technique made significant improvements to the population deformation heart models (p = 0.01). As standard evaluation, the residual Dice error after adaptation was comparable to the volumetric differences observed in free-breathing heart volumes (p = 0.62).
The reconstruction technique described generates accurate 3D heart models from limited 2D planning data. This development could potentially be used to retrospectively calculate delivered dose to the heart for historically treated patients and thereby provide a better understanding of late radiation-related cardiac effects.
Radiation-induced cardiac toxicity is a significant cause of morbidity, following treatment of Hodgkin's Lymphoma (HL) . Laboratory and autopsy studies have demonstrated that radiation therapy (RT) produces a range of adverse cardiac effects, including atherosclerosis, cardiomyopathy, and valvular damage [2, 3]. Several studies have investigated the relationship between late cardiac effects and radiation dose and/or radiated volume [4–7]. Emami et al. attributed the tolerance doses of cardiac tissues to therapeutic irradiation, but limited to only one type of radiation injury, pericarditis . Similarly, Stewart and Fajardo reported a steep dose response (pericarditis) in HL patients with an RT field that included > 50% of the external heart contour . In Hancock et al. an increased risk of coronary artery and cardiac diseases was reported when mediastinal radiation reached 40-45 Gy . Furthermore, existing normal-tissue complication probability models have been reported limiting due to the assumption that all regions of the organ at risk are of equal radiosensitivity . Therefore, a refinement of existing RT-related 3D dose-risk model of cardiac toxicity could potentially allow a better understanding of cardiac radiation tolerances, risk estimation of late cardiac effects, and optimization of RT treatment to minimize cardiac toxicity. However, due to the delayed onset of cardiac effects, the major barrier to the application of biophysical models remains the absence of 3D dosimetry on historically treated patients, for whom late toxicity outcomes are available, but were treated with a 2D treatment planning system. In order to obtain corresponding dose-volume and reliable late toxicity data, a 3D model reconstruction system is needed to extract 3D volume information from historical 2D planning datasets.
Various techniques have been described to reconstruct 3D organ volumes from 2D imaging data for the planning of RT and surgery [8–12]. The finite-element based deformable image registration system (FEM-DIR) has been utilized for 3D estimation of cardiac deformation , changes in brain interstitial pressure , and temporal changes in liver and lung shape due to breathing motion [10–12]. In particular, a recent development of the Navigator Channels (NC) has been used in conjunction with a biomechanical-based FEM-DIR platform, MORFEUS, to predict patient specific liver movement by adapting liver motion detected from coronal 2D images to a population liver motion model [13, 14]. This method can quickly associate 1D information obtained from 2D images with 3D organ information, potentially reproducing patient-specific 3D reconstructions of a given organ. In a previous investigation, single NCs were applied on 2D digitally reconstructed radiographs (DRRs) to capture differences in patients' lung structures . A linear adaptation approach was then used to accurately reconstruct patient-specific 3D lung volumes based on the structural differences . The current investigation builds on this technique, using an advanced NC-adaptation technique based on a series of NCs, and a combination of linear and bilinear adaptation approaches to reconstruct cardiac volumes directly from 2D images. This technique offers the potential to retrospectively correlate 3D dosimetry results with delayed cardiac toxicities for historically treated patients. Having a better understanding of dose-volume associations could in turn provide valuable insights into estimating late cardiac toxicities, and improve future HL radiation treatments.
The advanced NC adaptation technique combines the use of a biomechanical model-based deformable image registration system with multiple NCs to generate 3D heart volumes from 2D radiotherapy planning images. This technique was tested on 15 HL patients (Test patients) whose mediastinal RT was planned with a 3D CT-based treatment planning system. This technique expands on the previous lung reconstruction technique by utilizing multiple NCs on the right and left side of the heart (vs. single NCs for lung reconstruction) to delineate the deforming curvatures of the organ. While the heart is confined within the thoracic cage, its motion within its immediate boundaries (mediastinal pleura, adjacent to the lungs) is more variable than the lungs, which is more rigidly confined within the parietal pleura adjacent to the thoracic musculoskeletal wall. Therefore, in order to capture variability in heart shape, multiple NCs are used in this study to account for curvature differences at the circumflexes of the heart. Two-dimensional digitally reconstructed radiographs (DRRs), generated from 3D CT images, were used to simulate 2D fluoroscopic planning images available from historical 2D treatment plans.
B. Patient Data
A total of 50 sets of thoracic CT images from recently treated HL patients were obtained. Twelve image sets (6 male and 6 female Base patients) served as the basis for constructing the population heart models, 23 sets (Reference patients) were used to generate population-to-Reference patient deformation templates, and 15 additional sets (Test patients) were used to test the accuracy of the adaptation technique. Test patients that were treated with 3D treatment planning system between 1999 and 2008 were selected consecutively to avoid possible bias. This was then followed by consecutive selection of Reference patients that were also treated during this time period. The CT images were acquired at 120 kVp, 400 mA, with a voxel size of 0.1 × 0.1 × 0.3 cm. The patient was imaged with arms on hips, supine in an immobilization bag, under free breathing conditions. The acquired images were imported into a 3D CT-based treatment planning system (PINNACLE3 v7.6, Philips Radiation Oncology Systems, Madison, WI). Two-dimensional anterior DRRs were reconstructed from the CT images at an optimized resolution of 512 × 512 pixels, and source-to-axis distance (SAD) of 100 cm. The heart was highlighted on the DRR based on manual heart contours, to facilitate the placement of NCs during the adaptation process. For retrospective studies on historically treated patients, the heart would also be highlighted on the digitized 2D fluoroscopic planning images. Two-dimensional fluoroscopic images has higher soft tissue constrast than DRRs as they were generally taken at a lower kVp setting (80-100 kVp) at the time of conventional simulation. A graticule was included for image scaling of the DRRs.
C. Population-to-Patient Specific Deformation Models
Using an in-house FEM-DIR platform, MORFEUS [10, 15], a male and a female biomechanical-based finite-element population heart model were constructed from combining manually delineated heart contours on the CT images of 6 male or 6 female HL patients. This process involved 1) exporting heart contours as binary masks from the treatment planning system, 2) summation of the binary mask files using a logical OR function, and 3) conversion of the data into a 3D tetrahedral finite element volume mesh, using the Interactive Data Language software (IDL v6.3, ITT, Boulder, CO). The population models, which encompass all original Base patients' hearts, served as an initial 3D representation of the general shape and volume of the heart.
The male or female population heart model was then deformed into each of the 35 patients (12 Base patients that made up the population models and 23 Reference patients to generate 3D population-to-Reference patient deformation templates). Each deformation was carried out by first deforming the surface nodes of the FEM population heart model onto a surface defined by the secondary Reference heart model. This is followed by deformation of the internal structure of the heart according to the biomechanical properties assigned to the heart model (Young's modulus: 7.8 kPa and Poisson's ratio: 0.45) . Accuracy of each deformation was visually verified using the HyperMesh software (HYPERMESH v 9.0, Altair Engineering, Troy, MI). The 23 resultant population-to-Reference patient deformation template maps were then adapted to rebuild 3D heart models of 15 additional Test patients using only 2D DRRs, as described below.
D. Identification of Matching Patients
Matching Reference-Test pairs (RTPs) were identified based on lowest LSD scores (highest similarity) within their own gender group. A total of 15 RTPs were selected for this study. It is important to note that thorax measurements were used to identify patients with similar heart sizes, because radiographic film studies have expressed a strong relationship between thoracic measurements and cardiac sizes [16, 17]. In particular, the cardiothoracic ratio has traditionally been used as a basis for distinguishing age groups, ethnic groups, and patients with cardiac abnormality [18–20].
E. Quantification of Structural Differences
To refine selected Reference patient's deformation template maps for rebuilding Test patients' 3D heart models, the structural and positional differences between each RTP was quantified. The overall quantification process involved the following steps.
First, provided that the Reference and Test DRRs were initially reconstructed at an SAD of 100 cm, the user selected two points that corresponded to 10 scales on the graticule first on the Reference patient image then on the Test patient image. The software then automatically scaled the Test image to correspond to the Reference image (ie: 1 cm on the Test image corresponds to 1 cm on the Reference image).
Second, the anterior DRRs of the Reference and Test patients were rigidly aligned using the mid sternal notch (T4/T5), and the mid spinous processes of T6 and T8. This places the two images in the same comparable spatial location. The user then sets the x and y coordinates of the crosshair on the Reference image to the x and y coordinates of the isocenter on the Reference CT image as defined in the treatment plan. This provides a positional linkage between the 2D image and the 3D population-to-Reference patient deformation map.
Third, NCs were applied onto the two images. A NC was first placed on the Reference DRR by the user, then a spatially corresponding NC was automatically placed on the Test DRR by the NC software. A NC is a rectangular region of interest placed on two images that captures the intensity value differences within the region and generates an image intensity plot . The locations for placing the NCs were selected based on clearly defined boundaries of the heart as shown on the Reference and Test patient's radiograph (Figure 2b). To account for differences in the superior and inferior edge of the heart, a NC was placed superiorly and inferiorly. For Test cases where the inferior edge of the heart is difficult to outline on the 2D image as a result of inherent limitations of using a 2D radiographic film, the dome of the diaphragm was used to define that inferior edge. For Test cases where the superior edge is difficult to define due to blurring of adjacent structures with similar tissue densities (aorta, superior vena cava, pulmonary veins and arteries), the inferior edge of the carina, which is at the anatomical level of the pulmonary artery (located bilaterally and carries deoxygenated blood from the heart to the lungs) was used to define the superior edge of the heart. To account for lateral edge differences between Reference and individual Test patient's heart, locations where there is a significant change in curvatures between the Test and Reference patient heart structure were selected. Four NCs were generally placed at the major circumflexes of the heart, one at the edge of each of the four heart chambers (Right/Left atrium, Right/Left ventricles). In total, six NCs were placed at the edges of the heart: 2 vertical NCs (0.5 × 4.0 cm) for quantifying edge differences in the superior and inferior directions, 4 horizontal NCs (2.0 × 0.5 cm) for quantifying edge differences in the left and right directions at the 4 major circumflexes of the heart (Figure 2b).
Fourth, the image intensity within each NC was converted into a 1D intensity plot function. Finally, using a second derivative calculation, the difference or shift in intensity plots of each pair of NCs were calculated based on the differences of their points of inflection (POI) [13, 14]. POI occurs where the rate of change in the intensity values is the greatest, thus indicating the edge of the organ. Figure 2c illustrates this process for one of the NCs in the superior medial edge position.
F. Adaptation of Population Heart Model
The longitudinal adjustment of a node displacement (y 0 ) from population-to-Test patient specific, ∂(y 0 ), is equal to the predicted population-to-Reference patient deformation at the same node, Δ(y 0 ). This adjustment is refined by a linear weighting of the difference in the heart edges between the Reference and Test patient DRRs determined at each NC, ∂NC 5 (superior edge), and ∂NC 6 , (inferior edge), against the predicted displacement of the population-to-Reference patient deformation map at the NC node positions, ΔNC 5 and ΔNC 6 , respectively. The adjustment is further weighted by the relative distance of the adapting node to the two NC positions, k y and (1-k y ), where k y varies between 0 and 1 (0 is most superior and 1 is most inferior edge of the heart).
Equations (3), (4), and (5) account for the lateral adjustments of the nodes within the upper, middle, and lower regions of the heart, respectively. Four NCs were used because they help quantify major lateral edge differences of the four heart chambers (atrium & ventricles). At the same time, these 4 NCs provided a division of the heart that allowed for more accurate nodal adjustments, by deforming each node based on the NCs that were closest to it. Equations 3 and 5 are linearly weighted equations that were derived from the linear interpolation of the shifts calculated at its respective navigator channel positions. The position of the nodes in the upper region (Eq. 3) was calculated based on its relative position to NC3 and NC4. The position of the nodes in the lower region (Eq. 5) was calculated based on its relative position to NC1 and NC2. Based on the lateral shifts detected by its respective NC positions, equation 3 and 5 updates the position of each node in the population-Reference heart model Δ(y0) to become Test patient specific ∂(y0). Equation 4 is a bilinearly weighted equation that is derived from bilinear interpolating the shifts calculated at the four lateral NC positions, NC1, NC2, NC3, and NC4. Unlike linear interpolation, bilinear interpolation considers the closest 2 × 2 neighboring NC positions surrounding the unknown nodal position, Δ(x0). It then takes a weighted average of these four positions to arrive at its final, interpolated value, ∂(x0). The weight on each of the 4 NC positions is based on the computed nodal position (in 2D space) from each of the known NC positions. A bilinear equation is necessary, because it provides a more accurate deformation of the population-Reference model by accounting all possible shifts as a result of the differences detected by the four navigator channels.
In Equation (3), the lateral adjustment of a node displacement from population-to-Test patient specific, ∂(x 0 ), is equal to the predicted population-to-Reference patient displacement, Δ(x 0 ). This is adjusted by a linear weighting of the shifts calculated at NC3, ∂NC 3 , and NC4, ∂NC 4 , against the predicted population-to-Reference patient deformation at the NC3 and NC4 positions, ΔNC 3 and ΔNC 4 . The difference weighting was further adjusted by the relative distance of the adapting node to the two NC positions, k x and (1- k x ), where k x varies between 0 and 1 (0 is the most right and 1 is the most left edge of the heart).
In Equation (4), the deformation at a node, ∂(x 0 ), is equal to the population-to-Reference patient predicted deformation, Δ(x 0 ). This calculation is adjusted using bilinear interpolation of the edge differences between the Reference and Test DRRs, determined at the four NC positions, ∂NC 1 , ∂NC 2 , ∂NC 3 , and ∂NC 4 , against the predicted displacement of the population-to-Reference patient deformation map at the NC node positions, ΔNC 1 , ΔNC 2 , ΔNC 3 , and ΔNC 4 , respectively. The shift weighting is determined by the relative distance of the node to the four NC positions, k x , k y , (1- k x ), and (1- k x ) (where 0 is the most right and superior edge of heart).
Based on these equations, the dense population-to-Reference deformation field, which represents all the nodal coordinates of the Reference heart model was modified according to the position of each node relative to each NC location in order to calculate the deformation from population-to-Test patient. To ensure continuity, the interpolated displacements were then applied to each of the nodal points that make up the finite-element based population-Reference deformation heart model to reconstruct Test-patient-specific 3D heart models (Figure 2d).
G. Accuracy of Adapted Heart Models
The Dice coefficient was chosen in favor of simple matching similarity measure, because it does not "over-estimate" the percentage overlap of false-positive volumes . For comparison, the Dice coefficient expressed as a percentage between the original Reference (without NC adaptation) and its corresponding Test volume was calculated. Results with and without NC adaptation were compared using a Student's paired T-test.
To provide context for the accuracy of the reconstruction method, the Dice results were compared to the heart volume variation due to free-breathing lungs and diaphragm motion. 4D-CT images of 9 randomly selected lung cancer patients were analyzed. These CT images were acquired during synchronized respiration, in cine-mode, using a four-slice fan-beam CT scanner (Discovery LS, GE, Waukesha, WI). Sets of end-inhale and exhale images were transferred to the treatment planning system to assess the maximum variation in heart volume and position. Heart contours generated from both image sets were exported from the treatment plan as binary mask files and used to calculate their percentage volume overlap using the Dice coefficient (V1 = Inhale volume, V2 = Exhale volume). A paired T-test was performed between the NC-adapted and respiratory heart Dice results. This comparison was made, because the largest displacement occurs in the superior-inferior direction, as a result of free-breathing respiratory motion during radiotherapy . It is important to note that a rigid registration of the inhale and exhale volumes was not performed here because it is the intention of this study to determine how the accuracy of the NC reconstruction technique compare to the inherent positional and volumetric heart variation due to respiratory motion.
A. Population Deformations
B. Identification of Matching Reference-Test Pairs
Summary of LSD scores for the 15 RTPs
C. Accuracy of Adapted Heart Models
Summary of percentage volume overlap or Dice results
Volume Overlap Tests *
Ref. (Not NC Adapted) & Test
NC-adapted & Test
FB inhale & exhale subset
The reconstructed 3D heart volumes using the advanced NC adaptation process match well with the volumes obtained from CT. The Dice results demonstrate that the heart volume can be reconstructed with an accuracy of about 90%, which is comparable to the volumetric and/or positional variation due to respiration (p = 0.62). Statistically, it has been shown that the NCs made a significant improvement to the population deformation models (p = 0.01).
In this study, 3D heart volumes were reconstructed without lateral depth information. Although a "similar" Reference patient was identified (by thorax measurements) to approximate the lateral depth of the Test patient's heart, there remains an observable discrepancy in the anterior-posterior direction of the reconstructed models. This difference is largely due to the limitation that only anterior-posterior simulator films are available from historical 2D treatment plans. This limitation could introduce error in the computation of heart doses, as variations in heart-tissue density may cause the primary photons to become more or less attenuated than expected. However, current study results also indicated that this reconstruction error is no greater than the uncertainty imposed by a patient's free-breathing motion. Additionally, variation in day-to-day patient positioning will affect dose actually delivered to the heart during treatment [25, 26]. A previous study on the reproducibility of treatment setup for mantle irradiation in HL patients, using sequential port films, found that 29% of the films had minor deviations from original simulation film and 5% were considered unacceptable . Errors were noted in the superior and inferior mediastinum and could only be eliminated by resimulation . Taken together, we believe that the inaccuracies of the reconstruction technique described here are no greater than that caused by organ motion due to free-breathing and other setup uncertainties that are part of clinical practice.
The current NC adaptation technique is entirely based on soft-tissue alignment to allow for complete heart deformation. The manual contours of the heart were included on the DRR to facilitate NC placement. For retrospective studies, the same contours can be used for defining the heart on the Reference DRR. On digitized 2D fluoroscopic films of retrospective patients, the heart shape could also be highlighted as the soft tissue heart contour would in fact be better enhanced by the use of lower kVp at the time of conventional simulation. To further enhance tissue contrast on a 2D simulation film, a high-resolution diagnostic film scanner, combined with film processing software with an adjustable grayscale feature could be used. While it is recognized that overlapping boundaries (heart & ribs) may sometimes be difficult to define on radiographic films, the NC technique overcome this limitation by utilizing only sections of the organ edge in the adaption process. This is one of the major advantages of the NC technique over common segmentation techniques, which require complete organ boundary definition [27, 28]. Specifically, six NCs were applied only to sections of the edges that would most significantly change the shape of the heart (superiorly, inferiorly, & circumflexes of the four chambers) and the entire heart volume was quickly regenerated using the adaptation process.
Another key advantage of the advanced NC technique is its ability to delineate organ shape curvatures using a combination of linear and bilinear adaptation. This technique incorporates a bilinear method to allow multiple NCs to detect shape differences along the same organ edge. Potentially, cubic or bicubic interpolation methods could be used with multiple NCs to further improve reconstruction accuracy [21, 29]. However, the number of NCs employed should be organ-dependent, and calculation times should be balanced with accuracy outcomes. For example, heart volume overlap was computed using two additional NCs in the superior and inferior direction. The additional NCs did not result in a significant increase in volumetric accuracy.
The capacity to recreate 3D organ volumes from 2D planning data potentially allows for better correlation between organ dose and late onset toxicity. Investigation is underway to 1) reconstruct the 3D CT volumes (lungs & heart) and treatment plans of current patients for dosimetric validation of the current adaptation technique, and 2) reproduce the 3D dosimetry of retrospective 2D treatment plans in order to assess and refine existing dose, volume, and outcome relationships [4–7]. The latter process involves combining the use of the traditional film digitizer (ADAC Numonics-Accugrid Digitizer) with the 3D CT-based treatment planning system in order to incorporate shielding (lungs, heart, larynx, C-spine, humeral heads) from a 2D radiograph onto the adapted CT image.
The current development is a preliminary study, limited to the reconstruction of the overall heart volume. While this technique can offer a general dose-volume relationship of the heart, there remains several constraints when combined with 3D dosimetry data. First, localized radiation effects of specific cardiac structures remain unknown. Existing models tend to discard organ-specific spatial information, assuming all regions of the heart are of equal radiosensitivity [4–7]. However, this assumption is inaccurate, as the heart can be defined both anatomically (cardiac chambers, pericardium) and physiologically (cardiac cycles) [7, 30]. Second, the long term cardiac effects of radiation may be confounded with other treatment factors that may require further investigation. For instance, adjuvant chemotherapy is increasingly common for treatment of HL. Patients treated with chemotherapy (anthracyclines) alone were found to have increased cardiac toxicity [31, 32]. In contrast, non-chemo lipid-lowering treatments (statin) have demonstrated risk reduction in cardiac death [33, 34]. Therefore, in order to provide more accurate representation of the organ, future studies could aim at localization of the sub-regions of the heart, modeling the interaction of these sub-regions, and investigation of the combined late effects of radiation and other treatment regimens.
In summary, an advanced NC adaptation technique to reconstruct 3D heart volumes from 2D planning images has been developed. The reconstruction accuracy is comparable to the uncertainties observed in respiratory motion and thus would not be the limiting factor in estimating heart tissue exposure. A better understanding of the relationships between dose, volume and cardiac toxicity from conventional treatments can provide means for evaluation and refinement of existing dose-tissue constraints and NTCP models, in order to make improvements to modern radiotherapy techniques.
Deformable image registration
Digitally reconstructed radiograph
Finite-element based model
Least squares difference
Normal-tissue complication probability
Point of inflection
Inter-vertebral space between second and third thoracic vertebrae
Inter-vertebral space between ninth and tenth thoracic vertebrae.
Project funding provided by the Canadian Institutes for Health Research (Grant no. 103405). D. C. Hodgson and K. K. Brock are supported by Research Chairs from Cancer Care Ontario.
- Gustavsson A, Osterman B, Cavallin-Stahl E: A systematic overview of radiation therapy effects in Hodgkin's lymphoma. Acta Oncol. 2003, 42: 589-604. 10.1080/02841860310013346.View ArticlePubMed
- Gaya AM, Ashford RF: Cardiac complications of radiation therapy. Clin Oncol (R Coll Radiol). 2005, 17: 153-9. 10.1016/j.clon.2004.09.016.View Article
- Yeh ET, Bickford CL: Cardiovascular complications of cancer therapy: incidence, pathogenesis, diagnosis, and management. J Am Coll Cardiol. 2009, 53: 2231-47. 10.1016/j.jacc.2009.02.050.View ArticlePubMed
- Emami B, Lyman J, Brown A, Coia L, Goitein M, Munzenrider JE, Shank B, Solin LJ, Wesson M: Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys. 1991, 21: 109-22.View ArticlePubMed
- Stewart JR, Fajardo LF: Dose response in human and experimental radiation-induced heart disease. Application of the nominal standard dose (NSD) concept. Radiology. 1971, 99: 403-8.View ArticlePubMed
- Hancock SL, Donaldson SS, Hoppe RT: Cardiac disease following treatment of Hodgkin's disease in children and adolescents. J Clin Oncol. 1993, 11: 1208-15.PubMed
- Marks LB, Yorke ED, Jackson A, Ten Haken RK, Constine LS, Eisbruch A, Bentzen SM, Nam J, Deasy JO: Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys. 76: S10-9.
- Papademetris X, Sinusas AJ, Dione DP, Constable RT, Duncan JS: Estimation of 3-D left ventricular deformation from medical images using biomechanical models. IEEE Trans Med Imaging. 2002, 21: 786-800. 10.1109/TMI.2002.801163.View ArticlePubMed
- Miga MI, Paulsen KD, Hoopes PJ, Kennedy FE, Hartov A, Roberts DW: In vivo modeling of interstitial pressure in the brain under surgical load using finite elements. J Biomech Eng. 2000, 122: 354-63. 10.1115/1.1288207.View ArticlePubMed
- Brock KK, Hollister SJ, Dawson LA, Balter JM: Technical note: creating a four-dimensional model of the liver using finite element analysis. Med Phys. 2002, 29: 1403-5. 10.1118/1.1485055.View ArticlePubMed
- Brock KK, McShan DL, Ten Haken RK, Hollister SJ, Dawson LA, Balter JM: Inclusion of organ deformation in dose calculations. Med Phys. 2003, 30: 290-5. 10.1118/1.1539039.View ArticlePubMed
- Ng A, Nguyen T-N, Moseley LJ, Hodgson CD, Sharpe BM, Brock KK: Reconstruction of 3D Lung Models from 2D Planning Datasets for Hodgkin's Lymphoma Patients Using Combined Deformable Image Registration and Navigator Channels. Med Phys. 2010, 37: 12-10.1118/1.3263615.View Article
- Nguyen TN, Moseley JL, Dawson LA, Jaffray DA, Brock KK: Adapting liver motion models using a navigator channel technique. Med Phys. 2009, 36: 1061-73. 10.1118/1.3077923.PubMed CentralView ArticlePubMed
- Nguyen TN, Moseley JL, Dawson LA, Jaffray DA, Brock KK: Adapting population liver motion models for individualized online image-guided therapy. Proceedings of the Annual Conference of IEEE Engineering in Medicine and Biology Society: Vancouver. 2008, 3945-8.
- Brock KK, Sharpe MB, Dawson LA, Kim SM, Jaffray DA: Accuracy of finite element model-based multi-organ deformable image registration. Med Phys. 2005, 32: 1647-59. 10.1118/1.1915012.View ArticlePubMed
- Danzer CS: The cardio-thoracic ratio. An index of cardiac enlargement. Am J Med Sci. 1919, 157: 513-521. 10.1097/00000441-191904000-00007.View Article
- Bedford DE, Treadgold HA: The size of the healthy heart and its measurement. Lancet. 1931, 2: 836-842.
- Comeau WJ, White PD: A critical analysis of standard methods of estimating heart size from Roentgen measurements. Am J Roentgenol. 1942, 47: 665-677.
- Nickol K, Wade AJ: Radiographic heart size and cardiothoracic ratio in three ethnic groups: a basis for a simple screening test for cardiac enlargement in men. Br J Radiol. 1982, 55: 399-403. 10.1259/0007-1285-55-654-399.View ArticlePubMed
- Ungerleider HE, Gubner R: Evaluation of heart size measurements. Am Heart J. 1942, 24: 494-510. 10.1016/S0002-8703(42)90966-9.View Article
- Meijering E: A chronology of interpolation: from ancient astronomy to modern signal and image processing. Proceedings of the IEEE. 2002, 90: 319-342. 10.1109/5.993400.View Article
- Dice LR: Measures of the amount of ecologic association between species. Ecology. 1945, 26: 297-302. 10.2307/1932409.View Article
- Babalola KO, Patenaude B, Aljabar P, Schnabel J, Kennedy D, Crum W, Smith S, Cootes T, Jenkinson M, Rueckert D: An evaluation of four automatic methods of segmenting the subcortical structures in the brain. Neuroimage. 2009, 47: 1435-1447. 10.1016/j.neuroimage.2009.05.029.View ArticlePubMed
- Jagsi R, Moran JM, Kessler ML, Marsh RB, Balter JM, Pierce LJ: Respiratory motion of the heart and positional reproducibility under active breathing control. Int J Radiat Oncol Biol Phys. 2007, 68: 253-8. 10.1016/j.ijrobp.2006.12.058.PubMed CentralView ArticlePubMed
- McCord DL, Million RR, Northrop MF, Kavanaugh HV: Daily reproducibility of lung blocks in the mantle technique. Radiology. 1973, 109: 735-6.View ArticlePubMed
- Taylor BW, Mendenhall NP, Million RR: Reproducibility of mantle irradiation with daily imaging films. Int J Radiat Oncol Biol Phys. 1990, 19: 149-51.View ArticlePubMed
- Zhuang X, Rhode K, Arridge S, Razavi R, Hill D, Hawkes D, Ourselin S: An atlas-based segmentation propagation framework locally affine registration--application to automatic whole heart segmentation. Med Image Comput Comput Assist Interv. 2008, 11: 425-33.PubMed
- Ragan D, Starkschall G, McNutt T, Kaus M, Guerrero T, Stevens CW: Semiautomated four-dimensional computed tomography segmentation using deformable models. Med Phys. 2005, 32: 2254-61. 10.1118/1.1929207.View ArticlePubMed
- Pratt WK: Digital Image Processing. 1991, New York: Wiley
- Gagliardi G, Constine LS, Moiseenko V, Correa C, Pierce LJ, Allen AM, Marks LB: Radiation dose-volume effects in the heart. Int J Radiat Oncol Biol Phys. 76 (Suppl 3): 77-85.
- Elbl L, Vasova I, Navratil M, Tomaskova I, Jedlicka F, Chaloupka V, Vorlicek J: Late cardiotoxicity in patients with malignant lymphoma treated with doxorubicin chemotherapy. Vnitr Lek. 2006, 52: 328-38.PubMed
- Smith LA, Cornelius VR, Plummer CJ, Levitt G, Verrill M, Canney P, Jones A: Cardiotoxicity of anthracycline agents for the treatment of cancer: systematic review and meta-analysis of randomised controlled trials. BMC Cancer. 2010, 10: 337-10.1186/1471-2407-10-337.PubMed CentralView ArticlePubMed
- Chen AB, Punglia RS, Kuntz KM, Mauch PM, Ng AK: Cost effectiveness and screening interval of lipid screening in Hodgkin's lymphoma survivors. J Clin Oncol. 2009, 27: 5383-9. 10.1200/JCO.2009.22.8460.PubMed CentralView ArticlePubMed
- Prosser LA, Stinnett AA, Goldman PA, Williams LW, Hunink MG, Goldman L, Weinstein MC: Cost-effectiveness of cholesterol-lowering therapies according to selected patient characteristics. Ann Intern Med. 2000, 132: 769-79.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1756-6649/12/1/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.