Management of Respiratory Motion in Radiation Therapy



Management of Respiratory Motion in Radiation Therapy


Gig S. Mageras

Paul J. Keall



Increased interest in managing respiratory motion in radiation treatments has been prompted by several factors. First, there has been limited ability to control tumors in the thorax and abdomen with standard radiotherapy techniques. Second, anatomic movement with respiration, at least in some circumstances, limits the accuracy with which radiation can be delivered to tumor-bearing tissue. The resultant larger treatment volumes required to accommodate target mobility may limit the tumoricidal dose due to the larger amounts of surrounding normal tissue exposure, particularly for larger treatment volumes. Conversely, underestimation of the required treatment margins may result in marginal misses. Therefore, it seems desirable to limit respiratory motion for tumor sites exhibiting large excursions. Third, technologic advances have spawned new capabilities for measuring and reducing respiratory motion. This chapter discusses the problems associated with respiratory motion and surveys different strategies for managing respiratory motion during radiation treatment. A more detailed discussion can be found in the report of Task Group 76 of the American Association of Physicists in Medicine (AAPM).1

Respiratory motion affects image acquisition, treatment planning, and radiation delivery. Computed tomography (CT) images appear distorted because CT study times are long with respect to respiratory motion timescales, leading to a wrong appreciation of target and organ-at-risk (OAR) shapes and extents. The left panel of Figure 5.1 illustrates the effect of respiratory motion on standard CT acquisition.2 The object moves inferiorly-superiorly through the CT scan plane; hence, consecutive CT images occur at disordered positions in the object along the superior-inferior direction, and the object appears distorted. The right panel of Figure 5.1 illustrates this effect in a clinical image. The motion artifacts give rise to a systematic error in tumor position and extent (i.e., the tumor observed in the image is not the respiration-averaged position).3, 4, 5

Respiratory motion also limits target localization accuracy while imaging at treatment. A limitation of cone beam CT (CBCT) systems is image degradation caused by patient motion artifacts, which is a consequence of the limited gantry speed and the resultant long (approximately 1 minute) acquisition times. Respiratory motion, in particular, degrades images in the thorax and abdomen. Distortion in CBCT is completely different than for standard CT: Projection images occur at different respiration phases; thus, sagittal, coronal, and axial reconstructions are blurred, and axial images show streaking artifacts. The resultant difference between CT and CBCT leads to a systematic error in observed tumor position and shape.6

Treatment planning to account for respiratory motion poses additional problems. There is no motion information available from standard CT, so one often resorts to using generic margins rather than patient-specific ones. However, patients exhibit widely differing amounts of tumor mobility with respiration; thus, the choice of generic margins is likely to be too large or too small.

Respiratory motion during radiation delivery blurs the planned dose distribution, which in standard clinical practice is calculated without explicit inclusion of motion. For treatment fields with uniform radiation intensity within the field, respiratory motion will broaden the dose distribution in the anatomy moving near the beam edges. For fields delivered with physical compensators, the dose gradients in moving tissue will broaden and become less steep.

Intensity-modulated radiation treatment (IMRT) delivered with a multileaf collimator (MLC) poses additional considerations because of the interplay between moving MLC and respiratory motion.7 An intensity-modulated field is composed of many small fields that are delivered temporally; thus, the dose actually received by moving tissue may be less than, or greater than, the planned dose. Initial studies reported large dose variations, exceeding 20%, for a single field.8,9 Subsequent studies have found that for multiple field treatments with 30 fractions and assuming that the start of each treatment field occurs at a random point in the breathing cycle, dose variations average out to produce distributions that are the same as for treatments delivered with a physical compensator.10, 11, 12 Consideration of interplay effects should be given for limited-fraction IMRT and with scanning particle beams. Respiratory motion in liver can have a clinically important effect on proton fields (Fig. 5.2), causing underdose in the distal portion of the target volume or overdose to normal tissue beyond the distal portion.13

Interfraction variations in tumor position result in a displacement of the dose distribution. There are few published measurements on interfraction variation of tumor position in
lung because the intrafraction respiration-induced motion is more dominant and obscures the interfraction component.14,15 Interfraction variations can be decomposed into systematic and random components, similar to analyses of other geometric errors. A systematic error is one in which a quantity observed in the planning image set, such as target position, differs from its mean position over the course of treatment, whereas a random error is the deviation of the quantity from its mean position on any given treatment fraction. For a similar magnitude of error, a systematic displacement of the dose distribution has a larger effect on the dose to the target, whereas the effects of breathing motion and random displacements are smaller.6,15,16






Figure 5.1. A: Schematic of the effect of respiratory motion on standard computed tomography (CT) acquisition. (From Caldwell CB, Mah K, Skinner M, et al. Can PET provide the 3D extent of tumor motion for individualized internal target volumes? A phantom study of the limitations of CT and the promise of PET. Int J Radiat Oncol Biol Phys. 2003;55:1381-1393, with permission.) B: Clinical example of respiration-induced distortions.






Figure 5.2. Axial computed tomography (CT) images at the same table position but at two different points in the respiratory cycle of a patient receiving proton treatment of a liver lesion. Left image is at end inspiration (0% phase); right image is at end expiration (50% phase). Arrows indicate the beam path for right lateral (RL) and posterior (PA) treatment fields. From 0% to 50% phase, the water-equivalent depth increases by 3 mm and 12 mm for the RL and PA fields, respectively. (From Lu HM, Brett R, Sharp G, et al. A respiratory-gated treatment system for proton therapy. Med Phys. 2007;34:3273-3278, with permission.)


MEASUREMENT OF RESPIRATORY MOTION

Published observations of the extent of respiration-induced organ motion include lung, esophagus, liver, pancreas, breast,
kidney, spleen, and prostate.1,17,18 Studies have tracked the movement of tumor, host organ, radiographic fiducials implanted at the tumor site, and surrogate organs such as the diaphragm, which are assumed to correlate with the tumor. Patient breathing patterns can vary in magnitude, period, and regularity during and between treatment sessions, although the breathing pattern within a session is usually more stable than between sessions (Fig. 5.3). A study of three-dimensional (3D) trajectories of fiducials implanted in or near pulmonary tumors of 20 patients revealed hysteresis in half of the cases.19 The hysteresis amounted to a 1- to 5-mm difference in the trajectories between inhalation and exhalation (Fig. 5.4).






Figure 5.3. Motion amplitudes in the left-right (x), superior-inferior (y) and anterior-posterior (z) direction of a 2 mm gold fiducial implanted in a lung tumor. The fiducial was tracked in each treatment port of a patient treated with 4 ports per day for 4 days. Total recorded time was approximately 80 min. (From Shirato H, Seppenwoolde Y, Kitamura K, et al., Intrafractional tumor motion: lung and liver, Semin Radiat Oncol 2004;14:10-18 with permission.)

Recent studies have used respiration-correlated CT (RCCT; described further in the next section), which provides 3D information on anatomic position at several points in the respiratory cycle, to measure organ motion in the thorax and abdomen.20, 21, 22, 23, 24, 25 In the largest such statistical analysis, Liu et al.25 examined tumor motion in 152 patients with stage III or IV non-small-cell lung carcinoma (NSCLC). The analysis used rigid soft tissue-based image registration to measure tumor motion, which was found to be primarily in the superior-inferior (SI) direction. Ninety-five percent of tumors moved less than 1.3 cm SI, 0.4 cm left-right (LR), and 0.6 cm anterior-posterior (AP). Tumor motion was highly correlated with diaphragm excursion and tumor location in the SI direction. Brandner et al.21 reported on the motion of abdominal organs in 13 patients using RCCT. Motion was largest in the SI direction, with mean SI displacement during tidal breathing of 1.3 cm for the liver, 1.3 cm for the spleen, and 1.2 cm for the kidneys. Motion of up to 2.5 cm inferiorly was observed in all organs, and motion of up to 1.2 cm anteriorly was observed in liver and kidneys. The numerous studies indicate that patients exhibit a broad range of respiratory patterns and amplitudes of motion.


ACCOUNTING FOR RESPIRATORY MOTION AT SIMULATION

Different CT methods have been proposed for measuring respiration-induced tumor motion as a means of designing the treatment plan to encompass the motion. It is important to keep in mind that in all of these methods, respiratory patterns and, hence, tumor motion extent can change between simulation and treatment sessions.


SLOW COMPUTED TOMOGRAPHY SCANNING AND BREATH-HOLD COMPUTED TOMOGRAPHY

Slow CT scanning has been proposed for peripheral lung tumors,26 such that one or more respiratory cycles are recorded per CT slice; hence, the image of the tumor, at least when surrounded by lower density lung, shows the full extent of respiratory motion. An additional advantage is that the patient geometry is a respiration-averaged one; thus, the dose calculation is more representative of treatment than is fast CT. A disadvantage is motion blurring, which renders organ delineation in the images more difficult. In addition, the method can underestimate the motion extent of lung tumors attached to the mediastinum, chest wall, or diaphragm.







Figure 5.4. Orthogonal projections (left: coronal; right: sagittal) of trajectories of 21 fiducials implanted in or near pulmonary tumors. Circles denote tumors attached to skeletal structures. (From Seppenwoolde Y, Shirato H, Kitamura K, et al. Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. Int J Radiat Oncol Biol Phys. 2002;53:822-834, with permission.)

Two sets of breath-hold CT scans, one each at tidal end expiration and end inspiration, are often used to estimate tumor respiratory excursion.27 The method relies on the patient’s ability to hold his or her breath reproducibly and at the correct inspiration level; thus, some means of verifying inspiration level is desirable, such as spirometry or optical monitoring of abdominal displacement. Because the breath-hold scans will overestimate or underestimate lung volumes at treatment, a free-breathing or slow CT scan may also be acquired for treatment planning purposes.


RESPIRATION-CORRELATED COMPUTED TOMOGRAPHY/FOUR-DIMENSIONAL COMPUTED TOMOGRAPHY

RCCT, which is often referred to as four-dimensional (4D) CT, addresses both problems associated with free-breathing CT—image distortion caused by respiratory motion and lack of motion information.28, 29, 30, 31 The approach common to these methods is to acquire sufficient data for generating CT images at all phases of the respiratory cycle while simultaneously recording respiration, then retrospectively correlating the CT images with phase.

The concept of using cine CT acquisition for correlating with respiration was reported independently by Low et al.,29 who used spirometry for generating a respiratory signal, and Pan et al.,30 who used abdominal displacement. In both approaches, repeat CT images are acquired over slightly more than one respiratory cycle with the couch stationary while recording patient respiration; the couch is then incremented, and the process is repeated (Fig. 5.5). Following acquisition, the images are sorted with respect to the respiratory signal, leading to a set of volume images at different respiration points in the cycle.

Helical acquisition is performed by substantially reducing the pitch routinely used to approximately 0.1 to 0.5 and adjusting the acquisition parameters such that the CT beam will be on for at least one respiratory cycle at each couch position. This means that the couch velocity is much lower than for a conventional scan. With the CT beam on for at least one respiratory cycle at each couch position, there is sufficient oversampling of the respiratory cycle to create images, either reconstructed during acquisition or retrospectively using sinogram interpolation, at a number (typically eight to ten) of different states of the respiratory cycle. RCCT was originally demonstrated on single-slice scanners in helical mode.22,28,31 Subsequently, helical mode multislice RCCT acquisition has been demonstrated32,33 and is now commercially available.

Figure 5.6 shows an early clinical example of RCCT using cine acquisition. Artifacts are clearly reduced compared to standard CT, and anatomic motion is evident; however, some artifacts are still present, as indicated by the arrows, particularly at mid-expiration and mid-inspiration where motion is most rapid. Artifacts can be partly caused by too few samples per cycle. At a given desired point of reconstruction, at each table position, one chooses the image occurring at the nearest neighboring phase. Thus at each table position, the nearest neighbor image occurs at a slightly different phase. Early RCCT protocols with ten images per table position often caused artifacts. Current cine CT protocols acquire 20 to 25 images per table position to minimize this effect but at a price of typically 10 times the imaging dose of a standard CT, even with a reduced tube current.

In some commercial systems, CT images are sorted based on a respiration phase angle assigned by an algorithm that determines the period behavior of the respiratory signal. Phase-based sorting assumes repeatable breathing cycles (i.e., that the anatomy is at the same position for a given phase in every cycle). It has been recognized early on that this assumption is known to be violated for normal breathing, with its varying cycle-to-cycle amplitude variations, and has led to investigations into sorting based on the displacement of the
respiratory signal, in what is termed displacement-based sorting.34, 35, 36, 37, 38 Figure 5.7 shows the respiratory signal (tidal volume) of a patient that exhibits considerable variation in amplitude among breath cycles. Data points selected at end inspiration based on displacement (blue points in upper row) occur at similar values of tidal volume, whereas tidal volumes of data points based on phase (red points in lower row) show large variations. The resulting images reconstructed with displacement sorting (Fig. 5.8, upper row) show less artifacts than those with phase sorting (lower row).






Figure 5.5. Cine computed tomography (CT) scanning and image reconstruction. The duration of the x-ray on-interval at each couch position is greater than the respiration cycle plus the time to acquire data for one image. Dots on the respiratory signal (abdominal displacement) indicate the mid-scan time of a CT image; four such images are simultaneously acquired on a four-slice CT. There is an x-ray off-interval for couch increment to the next position. (From Pan T, Lee TY, Rietzel E, et al. 4D-CT imaging of a volume influenced by respiratory motion on multi-slice CT. Med Phys. 2004;31:333-340, with permission.)






Figure 5.6. Clinical example of respiration-correlated computed tomography (CT) using cine acquisition (LightSpeed GX/I; GE Medical Systems, Waukesha, Wis), with 10 images acquired per respiratory cycle to obtain 10 volume image sets. Numbers indicate percent phase calculated with a periodic filter algorithm (Real-time Position Management System; Varian Medical Systems, Palo Alto, Calif), where 0% corresponds to end inspiration.







Figure 5.7. Comparison of displacement-based versus phase-based sorting. Respiration signal is tidal volume measured with spirometry. Blue data points in the upper plot are classified as end inspiration (EI) by displacement sorting; red data points in the lower plot are classified by phase-based sorting. (From Lu W, Parikh PJ, Hubenschmidt JP, et al. A comparison between amplitude sorting and phase-angle sorting using external respiratory measurement for 4D CT. Med Phys. 2006;33:2964-2974, with permission.)

However, displacement-based correlation also has its limitations with irregular breathing. There may be no data available for a given displacement at some table positions, for example, during the shallow peaks in the upper row of Figure 5.7. A recent area of investigation has been the development of methods to further reduce irregular breathing artifacts in RCCT image sets by using deformable image registration to interpolate or extrapolate from (relatively) artifact-free images in the RCCT set.39






Figure 5.8. Images reconstructed using displacement-based sorting (upper row) and phase-based sorting (lower row) for the patient corresponding to Figure 5.7. Columns from left to right: end expiration, mid-inspiration, end inspiration, and mid-expiration. Arrows indicate locations of artifacts in the phase-sorted images. (From Lu W, Parikh PJ, Hubenschmidt JP, et al. A comparison between amplitude sorting and phase-angle sorting using external respiratory measurement for 4D CT. Med Phys. 2006;33:2964-2974, with permission.)


RESPIRATORY MOTION MEASUREMENT WITH POSITRON EMISSION TOMOGRAPHY AND MAGNETIC RESONANCE IMAGING

A number of different approaches have been examined to reduce the effects of respiratory motion on image quality in positron emission tomography (PET). Respiratory motion causes image blurring, reducing the signal-to-background contrast, decreasing the activity concentration per voxel within the lesion, and causing overestimation of lesion size. An initial approach to reducing motion is to use respiration-triggered or gated PET, in which PET data are acquired into discrete bins within each respiratory cycle.40 Clinical examples of gated PET have demonstrated its ability to reduce apparent lesion size and increase activity concentration (Fig. 5.9), albeit at the cost of reduced counts in each bin. Respiration-correlated PET has also been demonstrated, yielding similar results to gated PET but with the advantage of increased flexibility through retrospective data sorting.41 A disadvantage is that it requires considerably more computer memory and processing time than does gated PET. More recent approaches have studied deblurring of PET images by deconvolution methods,42 using deformable registration to fuse respiration-binned PET data into a single image43,44 and incorporating patient-specific motion model into PET image reconstruction.45







Figure 5.9. Axial fluorodeoxyglucose positron emission tomography (FDG-PET) image through a patient lesion in (A) nongated and (B) gated acquisition. Blue line in B indicates the active area of nongated lesion observed in A. (From Nehmeh SA, Erdi YE, Ling CC, et al. Effect of respiratory gating on quantifying PET images of lung cancer. J Nucl Med. 2002;43:876-881, with permission.)

The effects of respiratory motion can degrade the accuracy of free-breathing CT-derived attenuation corrections on combined PET/CT units, resulting in biased activity estimates in the PET images.46 The use of RCCT for applying attenuation corrections to gated PET images, in which CT-derived attenuation corrections correspond to the respiratory bin of the PET image, has been shown to improve co-registration between PET and CT images.47,48 In situations where respiration-averaged PET is standard clinical practice, use of respiration-averaged CT for calculating attenuation corrections improves tumor identification in PET images.49 For diagnostic imaging practice, Nehmeh et al.50 have reported on a deep inspiration breath-hold (DIBH) technique for acquiring PET/CT acquisition at a single breathing phase (Fig. 5.10).

Fast magnetic resonance imaging (MRI) acquisition techniques have been used for dynamic visualization of respiratory motion in the thorax and abdomen (Fig. 5.11).51, 52, 53 An advantage over RCCT is the continuous imaging over multiple breath cycles without ionizing radiation, thereby providing an assessment of variability in tumor mobility. von Siebenthal et al.54 describe a method of respiration-correlated volumetric MRI, using anatomic structures detected in a two-dimensional (2D) image (e.g., a sagittal slice at one location) for retrospectively sorting dynamic 2D images over the patient volume (Fig. 5.12).


FOUR-DIMENSIONAL PATIENT MODELING AND TREATMENT PLANNING

RCCT has several applications, including tumor motion measurement and motion management decision making, planning for respiratory-gated radiotherapy, planning for motion-inclusive radiotherapy, and planning for tumor tracking radiotherapy. In all such applications, it is important to note that even though useful motion information is included in the RCCT, the scan typically only includes information from a single breathing cycle per anatomic position. The breathing pattern will change with time, the tumor will change position with respect to the skeletal anatomy and other organs during a course of treatment, and tumor may change size and shape, all of which need to be accounted for with an appropriate treatment strategy.

The application of RCCT for tumor motion measurement is described earlier. However the magnitude of the observed motion can affect clinical decision making. If there is little motion observed or little normal tissue sparing advantage (e.g., constant cardiac-target distance with respiratory phase), explicit respiratory motion management procedures, such as gated or 4D planning techniques, may not be necessary. The AAPM Task Group 76 report1 gives a guideline that, based on other errors in radiotherapy, for motion less than 5 mm, explicit respiratory motion management may not be necessary. Figure 5.6 in the AAPM report gives a clinical decision process for a patient for whom respiratory motion is a concern during the radiotherapy process.

Planning for motion-inclusive, respiratory-gated, and tumor tracking radiotherapy can be put in a common theoretical framework, shown schematically in Figure 5.13. The mathematical description of the problem is to find the motion sequence M as a function of phase θ, M(θ), and gating beam on/off status H(θ), where the dose distribution to be delivered, D, is


where λ(θ) is the time spent in each respiratory phase, d is the dose distribution computed on each phase, I(θ) is a CT set for a given RCCT phase, Iref is the image set used as a reference, and u(θ) is the deformable image registration vector field linking anatomy in the different image sets. The vector field is particularly useful for 4D planning because it allows information related to the anatomy to be transferred between image sets of the 4D CT scan. For example, contours drawn on one phase can be transferred to the appropriate anatomy on other phases, or dose distributions from phases can be deformably summed to estimate the composite dose. The use of vector fields to automatically compute contours on subsequent phases or sum dose distributions is necessary given that 4D CT has an order of magnitude more data and an order of magnitude more manual effort in the absence of vector fields. However, deformable registration algorithms are an area of intense study, and variations with human observers invariably exist. The algorithms should be used with caution. At the time of press, deformable registration algorithms are not widely available. However, it is likely that in 2 to 3 years, all major treatment planning system vendors will offer some form of deformable registration to facilitate 4D planning.







Figure 5.10. Axial, coronal, and sagittal sections of (A) clinical free breathing and (B) deep inspiration breath-hold (DIBH) fused positron emission tomography (PET) and computed tomography (CT) images. Arrows indicate a lesion in rib, identified from the CT images. Lesion in the clinical PET images appears partially in lung and not fully aligned with the CT lesion because of respiratory motion. DIBH technique results in improved coregistration of the lesion between PET and CT. (From Nehmeh SA, Erdi YE, Meirelles GS, et al. Deep-inspiration breath-hold PET/CT of the thorax. J Nucl Med. 2007;48:22-26, with permission.)






Figure 5.11. Magnetic resonance imaging of tumor mobility of breathing cycle in coronal plane. After tidal respiration, the patient performed a deep inspiration followed by deep expiration. Small artifact from a metal clip is visible in tumor region. (From Plathow C, Ley S, Fink C, et al. Analysis of intrathoracic tumor mobility during whole breathing cycle by dynamic MRI. Int J Radiat Oncol Biol Phys. 2004;59:952-959, with permission.)


MOTION-INCLUSIVE TREATMENT PLANNING

Most treatment rooms do not currently have respiratory motion management devices. Therefore, treatment planning solutions are required for these situations. Based on the RCCT scan, a common approach is to determine the gross tumor volume (GTV) in each phase or a subset of the scan phases. The final GTV used for further treatment planning is the combination of each of the individual phases.55,56 Some researchers have noted that it may not be necessary to cover the full extent of tumor motion to derive an appropriate plan using inverse techniques.57 As with all respiratory motion management methods, the variation of breathing from cycle to cycle and the potential change in tumor position from day to day with respect to bony anatomy are challenges that warrant additional image guidance procedures.







Figure 5.12. Schematic of respiration-correlated volumetric magnetic resonance imaging. A: Sagittal slices covering the volume of interest. One dedicated slice N is used as a navigator for image sorting. B: Interleaved acquisition of data and navigator slices. (From von Siebenthal M, Szekely G, Gamper U, et al. 4D MR imaging of respiratory organ motion and its variability. Phys Med Biol. 2007;52:1547-1564, with permission.)






Figure 5.13. A Venn diagram showing the solution space for planning for the variables of gating (allowing the beam to turn on and off during treatment) and motion tracking (allowing the beam-target alignment to vary as a function of respiratory phase).


RESPIRATORY-GATED TREATMENT PLANNING

For respiratory-gated delivery,58,59 often the respiratory phases to be used for treatment are chosen based on a review of the RCCT, and phases in which the target position are in similar positions are used to compute both the treatment plan and duty cycle (delivery efficiency with gating). The duty cycle is typically planned to be 30% to 50%. During delivery, depending on breathing variability, the actual duty cycle, which is determined by the ratio of gate duration to the respiratory period, can vary from that planned. Often for thoracic and abdominal tumors, the exhalation phases are used for the treatment plan and treatment because the typical patient will spend more time in exhalation than inhalation and the exhalation position is more reproducible than inhalation. The advantage to gating at inhalation is that critical structures are further from the treatment field. During inhalation, the heart moves inferiorly and posteriorly, and thus, gated breast treatments at inhalation can reduce cardiac dose. Also, the lung expansion at inhalation (and especially for deep inhalation) means that the fraction of lung irradiated for a given beam aperture is smaller at inhalation than for other phases (discussed further later in this chapter).


MOTION-ENCOMPASSING TREATMENT METHODS

Depending on clinical treatment goals, it may be sufficient not to take any measures to control respiratory motion during treatment, but instead to use motion-encompassing treatment methods. In these circumstances, it is important that respiratory motion is properly accounted for during imaging for treatment planning, as discussed in the previous sections, as well as during imaging at treatment. Chapter 3 discusses margin determination.


A common approach in commercial systems is to define a target volume that encompasses the range of tumor excursion in all images of the RCCT. Software tools allow the user to loop through reconstructed slices in the three principal directions to evaluate the encompassing target volume. For lung tumors, a maximum intensity projection (MIP) tool is available in some systems, which produces an image in which each voxel is set to the maximum CT number for that voxel in the RCCT image set (Fig. 5.14).60 The MIP image provides a tumor motion-encompassing volume in situations where the tumor is surrounded by low-density lung tissue but will underestimate the encompassing volume in locations where the tumor is attached to mediastinum, chest wall, or diaphragm (Fig. 5.15).

Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

May 23, 2016 | Posted by in GENERAL RADIOLOGY | Comments Off on Management of Respiratory Motion in Radiation Therapy

Full access? Get Clinical Tree

Get Clinical Tree app for offline access