Search

NIR fluorescence imaging using bisphosphonate-Indocyanine green has been indicated for early interproximal caries detection. This study assessed diagnostic accuracy of caries detection by NIR fluorescence imaging with OsteoSense 750® (OS750) in vitro and ex vivo, and to analyze the therapeutic efficacy of a bisphosphonate (Etidronate) in inhibiting enamel caries progression in vitro. Methods: Four experiments were conducted using extracted human teeth; 1) to calculate the infiltration rate of OS750 into interproximal white spot lesions using fluorescence microscope, 2) to assess diagnostic accuracy of interproximal natural white spot lesions using desktop NIR fluorescence imaging device in vitro setting, 3) to assess diagnostic accuracy of artificially created deeper enamel carious lesion (0.5 mm~1.0 mm) using NIR fluorescence image through the head-mount display in ex vivo setting, 4) to compare the progression on the enamel caries lesions treated by Etidronate, NaF and distilled-water. Diagnostic accuracy was analyzed using sensitivity, specificity and receiver operating curves (ROC). The caries progression was calculated with micro-CT and was statistically analyzed using a two-way ANOVA and the Tukey HDS post-hoc test. Results: 1) The infiltration rate of OS750 was 101.83% ± 8.66 (Min: 90.10%, Max: 133.94%). 2) The average of sensitivity and specificity in vitro setting experiments were 86.7% ± 4.4% and 70% ± 11%, respectively. The average of area under the ROC curves (AUC) was 0.883 ± 0.059 indicating excellent performance. 3) The mean sensitivity and specificity in ex vivo setting was 82.97% ± 15% and 76.78% ± 13.27% respectively. 4) The carious lesion volume treated by Etidronate was significantly smaller at post treatment-1 (p

Weekly serial 4DCT scans were acquired under free breathing conditions to assess water-equivalent path length (WEL) variations due to both intrafractional and interfractional changes in tissue thickness and density and to calculate proton dose distributions resulting from anatomical variations observed in serial 4DCT. A template of region of interests (ROIs) was defined on the anterior–posterior (AP) beam's eye view, and WEL measurements were made over these ROIs to quantify chest wall thickness variations. Interfractional proton dose distributions were calculated to assess changes in the expected dose distributions caused by range variations. Mean intrafractional chest wall WEL changes during respiration varied by: −4.1 mm (

The superior dose conformation from protons is attributed to the Bragg peak near the end of the proton range. One challenge in proton cancer treatment is to assess the proton range fluctuations due to organ motion such as respiration. A time-resolved proton range telescope that measures coordinates, direction cosines, and the residual range of each proton can be useful in detecting and quantifying variations in radiological path length during the course of proton radiotherapy. In this paper, the Monte Carlo N-Particle eXtended (MCNPX) code was used to simulate the range telescope and study the image quality. To validate the MCNPX simulations, a simulated proton radiograph was compared with an experimentally acquired film for the same phantom. In addition, four quality assurance phantoms were simulated to investigate the quality of simulated proton radiography. Finally, the methods were applied to one phase of a patient four-dimensional computed tomography (4DCT) data set for proton radiography sim...

Purpose: Monte Carlo methods are used to simulate and optimize a time-resolved proton range telescope (TRRT) in localization of intrafractional and interfractional motions of lung tumor and in quantification of proton range variations. Methods: The Monte Carlo N-Particle eXtended (MCNPX) code with a particle tracking feature was employed to evaluate the TRRT performance, especially in visualizing and quantifying proton range variations during respiration. Protons of 230 MeV were tracked one by one as they pass through position detectors, patient 4DCT phantom, and finally scintillator detectors that measured residual ranges. The energy response of the scintillator telescope was investigated. Mass density and elemental composition of tissues were defined for 4DCT data. Results: Proton water equivalent length (WEL) was deduced by a reconstruction algorithm that incorporates linear proton track and lateral spatial discrimination to improve the image quality. 4DCT data for three patients were used to visualize and measure tumor motion and WEL variations. The tumor trajectories extracted from the WEL map were found to be within {approx}1 mm agreement with direct 4DCT measurement. Quantitative WEL variation studies showed that the proton radiograph is a good representation of WEL changes from entrance to distal of the target. Conclusions: MCNPX simulation results showedmore » that TRRT can accurately track the motion of the tumor and detect the WEL variations. Image quality was optimized by choosing proton energy, testing parameters of image reconstruction algorithm, and comparing to ground truth 4DCT. The future study will demonstrate the feasibility of using the time resolved proton radiography as an imaging tool for proton treatments of lung tumors.« less

The purpose of this study is to assess the temporal and reconstruction accuracy of a surface imaging system, the GateCT under ideal conditions, and compare the device with a commonly used respiratory surrogate: the Varian RPM. A clinical CT scanner, run in cine mode, was used with two optical devices, GateCT and RPM, to detect respiratory motion. A radiation detector, GM‐10, triggers the X‐ray on/off to GateCT system, while the RPM is directly synchronized with the CT scanner through an electronic connection. Two phantoms were imaged: the first phantom translated on a rigid plate along the anterior–posterior (AP) direction, and was used to assess the temporal synchronization of each optical system with the CT scanner. The second phantom, consisting of five spheres translating 3 cm peak‐to‐peak in the superior–inferior direction, was used to assess the quality of rebinned images created by GateCT and RPM. Calibration assessment showed a nearly perfect synchronization with the scanner for both the RPM and GateCT systems, thus demonstrating the good performance of the radiation detector. Results for the volume rebinning test showed discrepancies in volumes for the 3D reconstruction (compared to ground truth) of up to 36% for GateCT and up to 40% for RPM. No statistical difference was proven between the two systems in volume sorting. Errors are mainly due to phase detection inaccuracies and to the large motion of the phantom. This feasibility study assessed the consistency of two optical systems in synchronizing the respiratory signal with the image acquisition. A new patient protocol based on both RPM and GateCT will be soon started. PACS number: 87

Purpose: To estimate the position and volume errors in 4D-CT caused by irregular breathing. Methods: A virtual 4D-CT scanner was designed to reproduce axial mode scans with retrospective resorting. This virtual scanner creates an artificial spherical tumor based on the specifications of the user, and recreates images that might be produced by a 4D-CT scanner using a patient breathing waveform. 155 respiratory waveforms of patients were used to test the variability of 4D-CT scans. Each breathing waveform was normalized and scaled to 1, 2, and 3 cm peak-to-peak motion, and artificial tumors with 2 and 4 cm radius were simulated for each scaled waveform. The center of mass and volume of resorted 4D-CT images were calculated and compared to the expected values of center of mass and volume for the artificial tumor. Intrasubject variability was investigated by running the virtual scanner over different subintervals of each patient’s breathing waveform. Results: The average error in the center of mass location of an artificial tumor was less than 2 mm standard deviation for 2 cm motion. The corresponding average error in volume was less than 4%. In the worst-case scenarios, a center of mass error of 1.0 cm standard deviation and volume errors of 30%–60% at inhale were found. Systematic errors were observed in a subset of patients due to irregular breathing, and these errors were more pronounced when the tumor volume is smaller. Conclusions : Irregular breathing during 4D-CT simulation causes systematic errors in volume and center of mass measurements. These errors are small but depend on the tumor size, motion amplitude, and degree of breathing irregularity.

In this study, the authors investigated the feasibility of using a 3D liquid scintillator (LS) detector system for the verification and characterization of proton beams in real time for intensity and energy-modulated proton therapy. A plastic tank filled with liquid scintillator was irradiated with pristine proton Bragg peaks. Scintillation light produced during the irradiation was measured with a CCD camera. Acquisition rates of 20 and 10 frames per second (fps) were used to image consecutive frame sequences. These measurements were then compared to ion chamber measurements and Monte Carlo simulations. The light distribution measured from the images acquired at rates of 20 and 10 fps have standard deviations of 1.1% and 0.7%, respectively, in the plateau region of the Bragg curve. Differences were seen between the raw LS signal and the ion chamber due to the quenching effects of the LS and due to the optical properties of the imaging system. The authors showed that this effect can be accounted for and corrected by Monte Carlo simulations. The liquid scintillator detector system has a good potential for performing fast proton beam verification and characterization.

Surfaceimaging is in use in radiotherapy clinical practice for patient setup optimization and monitoring. Breast alignment is accomplished by searching for a tentative spatial correspondence between the reference and daily surface shape models. In this study, the authors quantify whole breast shape alignment by relying on texture features digitized on 3D surfacemodels. Texture feature localization was validated through repeated measurements in a silicone breast phantom, mounted on a high precision mechanical stage. Clinical investigations on breast shape alignment included 133 fractions in 18 patients treated with accelerated partial breast irradiation. The breast shape was detected with a 3D video based surfaceimagingsystem so that breathing was compensated. An in-house algorithm for breast alignment, based on surface fitting constrained by nipple matching (constrained surface fitting), was applied. Results were compared with a commercial software where no constraints are utilized (unconstrained surface fitting). Texture feature localization was validated within 2 mm in each anatomical direction. Clinical data show that unconstrained surface fitting achieves adequate accuracy in most cases, though nipple mismatch is considerably higher than residual surface distances (3.9 mm vs 0.6 mm on average). Outliers beyond 1 cm can be experienced as the result of a degenerate surface fit, where unconstrained surface fitting is not sufficient to establish spatial correspondence. In the constrained surface fitting algorithm, average surface mismatch within 1 mm was obtained when nipple position was forced to match in the [1.5; 5] mm range. In conclusion, optimal results can be obtained by trading off the desired overall surface congruence vs matching of selected landmarks (constraint). Constrained surface fitting is put forward to represent an improvement in setup accuracy for those applications where whole breast positional reproducibility is an issue.

Intrafractional motion and interfractional changes affect the accuracy of the delivered dose in radiotherapy, particularly in charged-particle radiotherapy. Most recent studies are focused on intrafractional motion (respiratory motion). Here, we report a quantitative simulation analysis of the effects of interfractional changes on water-equivalent pathlength (WEL) in charged-particle lung therapy. Serial four-dimensional (4D) CT scans were performed under free breathing conditions; the time span between the first and second 4DCT scans was five weeks. We quantified WEL changes between the first and second CT scans due to interfractional changes (tumor shrinkage and tissue density changes) and compared the particle-beam-stopping point between the serial 4DCT scans with use of the same initial bolus. Both tumor-shrinkage and lung-density changes were observed in a single patient over the course of therapy. The lung density decreased by approximately 0.1 g/cm(3) between the first and second-CT scans, resulting in a 1.5 cm WEL changes. Tumor shrinkage resulted in approximately 3 cm WEL changes. If the same initial bolus and plan were used through the treatment course, an unexpected significant beam overshoot would occur by interfractional changes due to tumor shrinkage and lung density variation.

Four-dimensional radiotherapy requires deformable registration to track delivered dose across varying anatomical states. Deformable registration based on B-splines was implemented to register 4D computed tomography data to a reference respiratory phase. To assess registration performance, anatomical landmarks were selected across ten respiratory phases in five patients. These point landmarks were transformed according to global registration parameters between different respiratory phases. Registration uncertainties were computed by subtraction of transformed and reference landmark positions. The selection of appropriate registration masks to separate independently moving anatomical subunits is crucial to registration performance. The average registration error for five landmarks for each of five patients was 2.1 mm. This level of accuracy is acceptable for most radiotherapy applications.

The susceptibility of IMRT dose distributions to intrafraction organ motion: An investigation into smoothing filters derived from four dimensional computed tomography data This study investigated the sensitivity of static planning of intensity-modulated beams (IMBs) to intrafraction deformable organ motion and assessed whether smoothing of the IMBs at the treatment-planning stage can reduce this sensitivity. The study was performed with a 4D computed tomography (CT) data set for an IMRT treatment of a patient with liver cancer. Fluence profiles obtained from inverse-planning calculations on a standard reference CT scan were redelivered on a CT scan from the 4D data set at a different part of the breathing cycle. The use of a nonrigid registration model on the 4D data set additionally enabled detailed analysis of the overall intrafraction motion effects on the IMRT delivery during free breathing. Smoothing filters were then applied to the beam profiles within the optimization process to investigate whether this could reduce the sensitivity of IMBs to intrafraction organ motion. In addition, optimal fluence profiles from calculations on each individual phase of the breathing cycle were averaged to mimic the convolution of a static dose distribution with a motion probability kernel and assess its usefulness. Results from nonrigid registrations of the CT scan data showed a maximum liver motion of 7 mm in superior-inferior direction for this patient. Dose-volume histogram (DVH) comparison indicated a systematic shift when planning treatment on a motion-frozen, standard CT scan but delivering over a full breathing cycle. The ratio of the dose to 50% of the normal liver to 50% of the planning target volume (PTV) changed up to 28% between different phases. Smoothing beam profiles with a median-window filter did not overcome the substantial shift in dose due to a difference in breathing phase between planning and delivery of treatment. Averaging of optimal beam profiles at different phases of the breathing cycle mainly resulted in an increase in dose to the organs at risk (OAR) and did not seem beneficial to compensate for organ motion compared with using a large margin. Additionally, the results emphasized the need for 4D CT scans when aiming to reduce the internal margin (IM). Using only a single planning scan introduces a systematic shift in the dose distribution during delivery. Smoothing beam profiles either based on a single scan or over the different breathing phases was not beneficial for reducing this shift. (C) 2006 American Association of Physicists in Medicine.

Purpose: To assess the utility of surface imaging on patient setup for accelerated partial-breast irradiation (APBI). Methods and Material: A photogrammetry system was used in parallel to APBI setup by laser and portal imaging. Surface data were acquired after laser and port-film setup for 9 patients. Surfaces were analyzed in comparison to a reference surface from the first treatment session by use of rigid transformations. The surface model after laser setup was used in a simulated photogrammetry setup procedure. In addition, breathing data were acquired by surface acquisition at a frame rate of 7 Hz. Results: Mean 3D displacement was 7.3 mm (SD, 4.4 mm) and 7.6 mm (SD, 4.2 mm) for laser and port film, respectively. Simulated setup with the photogrammetry system yielded mean displacement of 1 mm (SD, 1.2 mm). Distance analysis resulted in mean distances of 3.7 mm (SD, 4.9 mm), 4.3 mm (SD, 5.6 mm), and 1.6 mm (SD, 2.4 mm) for laser, port film, and photogrammetry, respectively. Breathing motion at isocenter was smaller than 3.7 mm, with a mean of 1.9 mm (SD, 1.1 mm). Conclusions: Surface imaging for PBI setup appears promising. Alignment of the 3D breast surface achieved by stereo-photogrammetry shows greater breast topology congruence than when patients are set up by laser or portal imaging. A correlation of breast surface and CTV must be quantitatively established.

External beam irradiation requires precise positioning of the target relative to the treatment planning coordinate system. A three-dimensional (3D) surface imaging system for patient positioning has recently been installed in one of our linear accelerator (linac) rooms. The device utilizes close-range photogrammetry to generate a 3D model of the patient's surface. This geometric model can be made to look like a digital camera image if wrapped with a gray-level image (texture mapping) that shows surface coloration. The system is calibrated to the linac coordinate system and has been designed as a patient setup device. To reproduce patient position in fractionated radiotherapy, the daily patient surface model is registered to a previously recorded reference surface. Using surface registration, the system calculates the rigid-body transformation that minimizes the distance between the treatment and the reference surface models in a region-of-interest (ROI). This transformation is expressed as a set of new couch coordinates at which the patient position best matches with the reference data. If respiratory motion is a concern, the surface can be obtained with a gated acquisition at a specified phase of the respiratory cycle. To analyze the accuracy of the system, we performed several experiments with phantoms to assess stability, alignment accuracy, precision of the gating function, and surface topology. The reproducibility of surface measurements was tested for periods up to 57 h. Each recorded frame was registered to the reference surface to calculate the required couch adjustment. The system stability over this time period was better than 0.5 mm. To measure the accuracy of the system to detect and quantify patient shift relative to a reference image, we compared the shift detected by the surface imaging system with known couch transitions in a phantom study. The maximum standard deviation was 0.75 mm for the three translational degrees of freedom, and less than 0.1 degrees for each rotation. Surface model precision was tested against computed tomography (CT)-derived surface topology. The root-mean-square rms of the distance between the surfaces was 0.65 mm, excluding regions where beam hardening caused artifacts in the CT data. Measurements were made to test the gated acquisition mode. The time-dependent amplitude was measured with the surface imaging system and an established respiratory gating system based on infrared (IR)-marker detection. The measured motion trajectories from both systems were compared to the known trajectory of the stage. The standard deviations of the amplitude differences to the motor trajectory were 0.04 and 0.15 mm for the IR-marker system and the 3D surface imaging system, respectively. A limitation of the surface-imaging device is the frame rate of 6.5 Hz, because rapid changes of the motion trajectory cannot be detected. In conclusion, the system is accurate and sufficiently stable to be used in the clinic. The errors computed when comparing the surface model with CT geometry were submillimeter, and deviations in the alignment and gating-signal tests were of the same magnitude.

Respiratory gating based on external surrogates is performed in many clinics. We have developed a new technique for treatment verification using an electronic portal imaging device (EPID) in cine mode for gated 3D conformal therapy. Implanted radiopaque fiducial markers inside or near the target are required for this technique. The markers are contoured on the planning CT set, enabling us to create digitally reconstructed radiographs (DRRs) for each treatment beam. During the treatment, a sequence of EPID images can be acquired without disrupting the treatment. Implanted markers are visualized in the images and their positions in the beam's eye view are calculated off-line and compared to the reference position by matching the field apertures in corresponding EPID and DRR images. The precision of the patient set-up, the placement of the beam-gating window, as well as the residual tumour motion can be assessed for each treatment fraction. This technique has been demonstrated with a case study patient, who had three markers implanted in his liver. For this patient, the intra-fractional variation of all marker positions in the gating window had a 95% range of 4.8 mm in the SI direction (the primary axis of motion). This was about the same (5 mm) as the residual motion considered in the planning process. The inter-fractional variation of the daily mean positions of the markers, which indicates the uncertainty in the set-up procedure, was within +8.3 mm/-4.5 mm (95% range) in the SI direction for this case.

Purpose: The correlation of the respiratory motion of external patient markers and abdominal tumors was examined. Data of this type are important for image-guided therapy techniques, such as respiratory gating, that monitor the movement of external fiducials. Methods and Materials: Fluoroscopy sessions for 4 patients with internal, radiopaque tumor fiducial clips were analyzed by computer vision techniques. The motion of the internal clips and the external markers placed on the patient’s abdominal skin surface were quantified and correlated. Results: In general, the motion of the tumor and external markers were well correlated. The maximum amount of peak-to-peak craniocaudal tumor motion was 2.5 cm. The ratio of tumor motion to external-marker motion ranged from 0.85 to 7.1. The variation in tumor position for a given external-marker position ranged from 2 to 9 mm. The period of the breathing cycle ranged from 2.7 to 4.5 seconds, and the frequency patterns for both the tumor and the external markers were similar. Conclusions: Although tumor motion generally correlated well with external fiducial marker motion, relatively large underlying tumor motion can occur compared with external-marker motion and variations in the tumor position for a given marker position. Treatment margins should be determined on the basis of a detailed understanding of tumor motion, as opposed to relying only on external-marker information. © 2005 Elsevier Inc. Image-guided therapy, Gating, Fluoroscope, Abdomen, Tracking.

Respiratory motion can introduce significant errors in radiotherapy. Conventional CT scans as commonly used for treatment planning can include severe motion artifacts that result from interplay effects between the advancing scan plane and object motion. To explicitly include organ/target motion in treatment planning and delivery, time-resolved CT data acquisition (4D Computed Tomography) is needed. 4DCT can be accomplished by oversampled CT data acquisition at each slice. During several CT tube rotations projection data are collected in axial cine mode for the duration of the patient's respiratory cycle (plus the time needed for a full CT gantry rotation). Multiple images are then reconstructed per slice that are evenly distributed over the acquisition time. Each of these images represents a different anatomical state during a respiratory cycle. After data acquisition at one couch position is completed, x rays are turned off and the couch advances to begin data acquisition again until full coverage of the scan length has been obtained. Concurrent to CT data acquisition the patient's abdominal surface motion is recorded in precise temporal correlation. To obtain CT volumes at different respiratory states, reconstructed images are sorted into different spatio-temporally coherent volumes based on respiratory phase as obtained from the patient's surface motion. During binning, phase tolerances are chosen to obtain complete volumetric information since images at different couch positions are reconstructed at different respiratory phases. We describe 4DCT image formation and associated experiments that characterize the properties of 4DCT. Residual motion artifacts remain due to partial projection effects. Temporal coherence within resorted 4DCT volumes is dominated by the number of reconstructed images per slice. The more images are reconstructed, the smaller phase tolerances can be for retrospective sorting. From phantom studies a precision of about 2.5 mm for quasiregular motion and typical respiratory periods could be concluded. A protocol for 4DCT scanning was evaluated and clinically implemented at the MGH. Patient data are presented to elucidate how additional patient specific parameters can impact 4DCT imaging.

Summary Background and Purpose Clinical target volumes of the thorax and abdomen are typically expanded to account for inter- and intrafractional organ motion. Usually, such expansions are based on clinical experience and planar observations of target motion during simulation. More precise, 4- dimensional motion margins for a specific patient may improve dose coverage of mobile targets and yet limit unnecessarily large field expansions. We are studying approaches to targeting moving tumors throughout the entire treatment process, from treatment planning to beam delivery. Material and Methods Radio-opaque markers were implanted under CT guidance in the liver at the gross tumor periphery. Organ motion during light respiration was volumetrically imaged by 4D Computed Tomography. Marker motion was also acquired by fluoroscopy and compared with 4DCT data. During treatment, daily diagnostic x-ray images were captured at end-exhale and -inhale for patient set-up and target localization. Results Based on the time-resolved CT data, target volumes can be designed to account for respiratory motion during treatment. Motion of the tumor as derived from 4DCT was consistent with fluoroscopic motion analysis. Radiographs acquired in the treatment room enabled millimeter-level patient set-up and assessment of target position relative to bony anatomy. Daily positional variations between bony anatomy and implanted markers were observed. Conclusions Image guided therapy, based on 4DCT imaging, fluoroscopic imaging studies, and daily gated diagonstic energy set-up radiographs is being developed to improve beam delivery precision. Monitoring internal target motion throughout the entire treatment process will ensure adequate dose coverage of the target while sparing the maximum healthy tissue.

To develop a dose-response relationship for the occurrence of permanent alopecia after cranial irradiation and to analyze potential confounding variables that may contribute to this unwanted and often unavoidable complication of treatment.Twenty-six patients were enrolled in this study. Three reviewers independently assessed 61 scalp regions and assigned a score for the degree of alopecia in each region using a 4-point scale. Patient and treatment data were collected using a patient questionnaire and outpatient medical chart review. The hair follicle dose was calculated for each scalp region and correlated with the alopecia score for that region. A dose-response relationship was established using the data from these correlations.Permanent alopecia correlated significantly with the follicle dose only (p0.001). A personal history of alopecia and the use of chemotherapy correlated with permanent alopecia with borderline statistical significance (p = 0.059 and p = 0.068, respectively). Patient age, family history of baldness, gender, tobacco use, diabetes, and beam energy did not correlate with alopecia.We report the first human dose-response relationship describing the effect of the follicle dose on the subsequent development of permanent scalp alopecia after cranial irradiation. This information will assist the radiation oncologist, physicist, and dosimetrist in designing a treatment plan that might minimize the risk of this untoward side effect of cranial irradiation.

We propose a new scanning protocol for generating 4D-CT image data sets influenced by respiratory motion. A cine scanning protocol is used during data acquisition, and two registration methods are used to sort images into temporal phases. A volume is imaged in multiple acquisitions of 1 or 2 cm length along the cranial-caudal direction. In each acquisition, the scans are continuously acquired for a time interval greater than or equal to the average respiratory cycle plus the duration of the data for an image reconstruction. The x ray is turned off during CT table translation and the acquisition is repeated until the prescribed volume is completely scanned. The scanning for 20 cm coverage takes about 1 min with an eight-slice CT or 2 mins with a four-slice CT. After data acquisition, the CT data are registered into respiratory phases based on either an internal anatomical match or an external respiratory signal. The internal approach registers the data according to correlation of anatomy in the CT images between two adjacent locations in consecutive respiratory cycles. We have demonstrated the technique with ROIs placed in the region of diaphragm. The external approach registers the image data according to an externally recorded respiratory signal generated by the Real-Time Position Management (RPM) Respiratory Gating System (Varian Medical Systems, Palo Alto, CA). Compared with previously reported prospective or retrospective imaging of the respiratory motion with a single-slice or multi-slice CT, the 4D-CT method proposed here provides (1) a shorter scan time of three to six times faster than the single-slice CT with prospective gating; (2) a shorter scan time of two to four times improvement over a previously reported multi-slice CT implementation, and (3) images over all phases of a breathing cycle. We have applied the scanning and registration methods on phantom, animal and patients, and initial results suggest the applicability of both the scanning and the registration methods.

The design of an integrated radiotherapy imaging system (IRIS), consisting of gantry mounted diagnostic (kV) x-ray tubes and fast read-out flat-panel amorphous-silicon detectors, has been studied. The system is meant to be capable of three main functions: radiographs for three-dimensional (3D) patient set-up, cone-beam CT and real-time tumour/marker tracking. The goal of the current study is to determine whether one source/panel pair is sufficient for real-time tumour/marker tracking and, if two are needed, the optimal position of each relative to other components and the isocentre. A single gantry-mounted source/imager pair is certainly capable of the first two of the three functions listed above and may also be useful for the third, if combined with prior knowledge of the target's trajectory. This would be necessary because only motion in two dimensions is visible with a single imager/source system. However, with previously collected information about the trajectory, the third coordinate may be derived from the other two with sufficient accuracy to facilitate tracking. This deduction of the third coordinate can only be made if the 3D tumour/marker trajectory is consistent from fraction to fraction. The feasibility of tumour tracking with one source/imager pair has been theoretically examined here using measured lung marker trajectory data for seven patients from multiple treatment fractions. The patients' selection criteria include minimum mean amplitudes of the tumour motions greater than 1 cm peak-to-peak. The marker trajectory for each patient was modelled using the first fraction data. Then for the rest of the data, marker positions were derived from the imager projections at various gantry angles and compared with the measured tumour positions. Our results show that, due to the three dimensionality and irregular trajectory characteristics of tumour motion, on a fraction-to-fraction basis, a 'monoscopic' system (single source/imager) is inadequate for consistent real-time tumour tracking, even with prior knowledge. We found that, among the seven patients studied with peak-to-peak marker motion greater than 1 cm, five cases have mean localization errors greater than 2 mm and two have mean errors greater than 3 mm. Because of this uncertainty associated with a monoscopic system, two source/imager pairs are necessary for robust 3D target localization. Dual orthogonal x-ray source/imager pairs mounted on the linac gantry are chosen for the IRIS. We further studied the placement of the x-ray sources/panel based on the geometric specifications of the Varian 21EX Clinac. The best configuration minimizes the localization error while maintaining a large field of view and avoiding collisions with the floor/ceiling or couch.

Target volumes in the thorax and abdomen are commonly computed tomography (CT) scanned during light respiration. In this article, we analyze the distortions introduced in helical scanning of moving objects. Objects of known geometry are placed on a moving sled and scanned in a multirow helical CT scanner. The motion of the sled approximates the magnitude and velocity of organ movement in patients during light respiration (amplitude 1 cm, period 4 seconds). Scans of the phantom are obtained in high speed mode at incremental phases of respiration, and the resulting images are compared with scans obtained when the phantom is static. Computer simulations of the scan process are also performed to interpret the results and extend the analysis to a greater range of parameters in scanning, motion, and object size. Resulting scans show that spherical test objects can be shortened by as much as 2 cm or twice the periodic motion amplitude. Object shape was significantly distorted, and the geometric center of the object was displaced by as much as +/-0.8 cm. Computer simulation results qualitatively agree with the experimentally observed phantom images. These simulations predict that the effect is clearly observable even if the amplitude is decreased to 0.5 cm. Implications of scanning moving objects on treatment planning are discussed.

Purpose After vertebral or paravertebral tumor resection, tumor cells may remain on the dura. Because a tumoricidal dose is difficult to achieve using external beam radiotherapy without exceeding the spinal cord tolerance, we developed intraoperative applicators to deliver additional dose to the dura. Methods and materials Eight patients with vertebral or paravertebral tumor underwent conformal external beam radiotherapy, tumor resection, and intraoperative radiotherapy to the dura involved by tumor. At surgery, vertebra, soft tissue, and epidural tumor were resected. A radioactive applicator plaque was placed on the dura to deliver 7.5–15 Gy, and then removed. Vertebral reconstruction and stabilization was completed. Chemotherapy was administered for large, high-grade sarcomas. Results We progressed through three plaque designs, initially 192 Ir, subsequently liquid 90 Y, and finally 90 Y foil in a semicylindrical polycarbonate plaque, in the treatment of 8 patients. The low-energy 90 Y β-emissions provided a more attractive depth dose profile than that achievable with iridium and gave negligible staff radiation exposure. The 90 Y depth dose measured 29% at 2 mm and 9% at 4 mm from the surface of the foil plaque, with acceptable surface dose homogeneity. The average surface dose rate ranged from 18.7 to 47.6 cGy/min for the iridium plaques and 45.2 to 187.5 cGy/min for the 90 Y plaques. The treatments have been without acute or late neurologic complications. The disease of 6 of 8 patients was locally controlled at median potential follow-up of 24 months. Conclusions The 90 Y foil applicator is technically elegant, easy to use, and superior to the earlier models. It has been incorporated into a protocol for spinal tumor treatment.

Respiration-induced tumour motion can potentially compromise the use of intensity-modulated radiotherapy (IMRT) as a dose escalation tool for lung tumour treatment. We have experimentally investigated the intra-fractional orga nm otion effects in lung IMRT treatments delivered by multi-leaf collimator (MLC). An in-house made motor-driven platform, which moves sinusoidally with an amplitude of 1 cm and a period of 4 s, was used to mimic tumour motion. Tumour motion was simulated along cranial-caudal direction while MLC leaves moved across the patient from left to right, as in most clinical cases. The dose to a point near the centre of the tumour mass was measured according to geometric and dosimetri cp arameters from two five-field lung IMRT plans. For each field, measurement was done for two dose rates (300 and 500 MU min −1 ), three MLC delivery modes (sliding window, step-and-shoot with 10 and 20 intensity levels) and eight equally spaced starting phases of tumour motion. The dose to the measurement point delivered from all five field sw as derived for both a single fraction and 30 fractions by randomly sampling from measured dose values of each field at different initial phases. It was found that the mean dose to a moving tumour differs slightly (

Radiation treatment planning (RTP), historically an image-intensive discipline and one of the first areas in which 3D information from imaging was clinically applied, has become even more critically dependent on accurate 3D definition of target and non-target structures in recent years with the advent of conformal radiation therapy. In addition to the interactive display of wireframe or shaded surface models of anatomic objects, proposed radiation beams, beam modifying devices, and calculated dose distributions, recently significant use has been made of direct visualization of relevant anatomy from image data. Dedicated systems are commercially available for the purpose of geometrically optimizing beam placement, implementing in virtual reality the functionality of standard radiation therapy simulators. Such “CT simulation” systems rely heavily on 3D visualization and on reprojection of image data to produce simulated radiographs for comparison with either diagnostic-quality radiographs made on a simulato...

Respiratory correlated CT is commercially available, and we have implemented its routine clinical use in planning lungtumor patients. Its value is determined by the fidelity of the spatiotemporal data set after processing the acquired reconstructed slices. Retrospective sorting of reconstructed slices is based on respiratory phase. However, the existing commercial software inadequately models respiratory phase for about 30% of the patients, mainly due to irregularities in the respiratory cycle. We have developed software that improves phase determination and consequently leads to an improvement of retrospective data sorting to make 4DCT data acquisition feasible for routine clinical use. Peak inhalation and exhalation respiratory states are selected manually; intermediate phases are interpolated. Residual motion artifacts in the resulting 4DCT volumes are reduced and allow use of the 4D imaging studies for treatment planning.

Purpose: We have developed and tested an interactive video system that utilizes image subtraction techniques to enablehigh precision patient repositioning using surface features. We report quantitative measurements of system performance characteristics. Methods and Materials: Video images can provide a high precision, low cost measure of patient position. Image subtraction techniques enable one to incorporate detailed information contained in the image of a carefully verified reference position into real-time images. We have developed a system using video cameras providing orthogonal images of the treatment setup. The images are acquired, processed and viewed using an inexpensive frame grabber and a PC. The subtraction images provide the interactive guidance needed to quickly and accurately place a patient in the same position for each treatment session. We describe the design and implementation of our system, and its quantitative performance, using images both to measure changes in position, and to achieve accurate setup reproducibility. Results: Under clinical conditions (60 cm field of view, 3.6 m object distance), the position of static, high contrast objects could be measured with a resolution of 0.04 mm (rms) in each of two dimensions. The two-dimensional position could be reproduced using the real-time image display with a resolution of 0.15 mm (rms). Two-dimensional measurement resolution of the head of a patient undergoing treatment for head and neck cancer was 0.1 mm (rms), using a lateral view, measuring the variation in position of the nose and the ear over the course of a single radiation treatment. Three-dimensional repositioning accuracy of the head of a healthy volunteer using orthogonal camera views was less than 0.7 mm (systematic error) with an rms variation of 1.2 mm. Setup adjustments based on the video images were typically performed within a few minutes. The higher precision achieved using the system to measure objects than to reposition them suggests that the variability in repositioning is dominated by the ability of the therapist to make small, controlled changes in the position of the patient. Conclusion: Using affordable, off-the-shelf technology, we have developed a patient positioning system that achieves repositioning accuracy normally associated with fractionated stereotactic systems. The technique provides real-time guidance and can be used to easily and quickly correct patient setup before every treatment, thus significantly reducing overall random positioning error. This improved positioning capability provides the precision required to realize the potential gains of conformal radiotherapy. 0 1997 Elsevier Science Inc.

An efficient technique to optimize beam weights and wedge angles in radiotherapytreatment planning has been developed. Based on the fact that a wedged field can be regarded as a superposition of an open field and a nominal wedged field, this approach reduces the problem of finding J beam weights and the corresponding wedge angles to optimizing a linear system with 2J unknowns (weights of J open beams and J nominal wedged beams), where J is the total number of incident beam directions. Two iterative algorithms similar to the iterative-least-square technique in image reconstruction are used to optimize the system. Application of the algorithms to two specific examples shows that this technique can reduce treatment planning time and effort, and promises to create a better solution for an arbitrarily complex treatment plan.

Functional dose - volume histograms are proposed as an extension of the conventional dose - volume histograms, for quantitative assessment of three-dimensional radiation dose coverage of functionally heterogeneous normal organs. Examples are given to illustrate possible applications of this approach to the treatment of a brain tumour or a lung tumour, in which cases the distribution of the normal organ function can be obtained from functional imaging modalities. It is shown that a significant difference exists between the functional dose - volume histograms and the conventional dose - volume histograms when the normal organ function is non-uniformly distributed within the organ. Utilization of functional dose - volume histograms as the input for the calculation of normal tissue complication probabilities is discussed for different normal tissue structures.

Purpose: Explore the use of functional imaging data in radiation treatment planning of brain lesions. Methods and Materials: Compare the treatment-planning process with and without the use of functional brain imaging for clinical cases where functional studies using either single photon emission computed tomography or magnetic resonance imaging are available. Results: A method to register functional image data with planning image studies is needed for functional treatment planning. Functional volumes are not simply connected regions. One activation study may produce many isolated functional areas. After finding the functional volumes and registering the functional information with the planning imaging data, the tools used for conventional three-dimensional treatment planning are sufficient for functional treatment planning. However, the planning system must provide dose-volume histograms for volumes of interest that consist of isolated pieces. Treatment plans that spare functional brain while providing identical target coverage can be constructed for lesions situated near the functional volume. However, the dose to other areas of the brain may be increased. Conclusions: Functional imaging will make determination of dose response of eloquent areas of the brain possible when combined with volumetric dose information and neuropsychological prior to and after radiation therapy. Realizing the full potential of functional imaging studies will require improved delineation of activated volumes and determination of the uncertainties in functional volume delineation. Optimization of treatment plans by minimizing dose to volumes activated during functional imaging studies should be used cautiously, because the dose to "silent," but possibly eloquent, brain may be increased.

Conventional source localization algorithms require a one-to-one source correspondence between films. This requirement makes source localization cumbersome and error prone because multiple sources must be carefully digitized and some sources can be obscured or missed. A new source localization algorithm is described in this paper. The algorithm fits a ribbon or needle image on film to a linear-quadratic equation, then analytically determines the 3-D ribbon locus by its image on the other projection, and finally localizes the sources in the ribbon by tracing along the ribbon image. Only three points per ribbon per film are required, and corresponding points need not be identified on the other film. Phantom experiments and tests on clinical cases demonstrate that the source localization algorithm can increase the efficiency by a factor of up to 5, improve accuracy to about 1 mm, and reconstruct obscured or shifted sources without decreased accuracy and efficiency. The simplicity and minimal entry of data make this technique desirable for clinical use.

To determine the variability of patient positioning during three-dimensional conformal radiotherapy (3D-CRT) for prostate cancer treated with no immobilization or one of four immunobilization devices, and to determine the effects of patient body habitus and pelvic circumference on patient movement with each individual inmobilization technique.To see whether our immobilization techniques have improved day-to-day patient movement, a retrospective analysis was carried out. A total of 62 patients treated at one facility on a single machine with 3D-CRT via a four-field box technique (anterior-posterior and opposed laterals) in the supine position with either no immobilization or one of four immobilization devices. Five groups of patients were compared: (a) group 1-no immobilization; (b) group 2-alpha cradle from the waist to upper thigh; (c) group 3-alpha cradle from waist to below the knees; (d) group 4-styrofoam leg immobilizer (below knees); and (e) group 5-aquaplast cast encompassing the entire abdomen and pelvis to midthigh with alpha cradle immobilization to their lower legs and feet. Prior to starting radiotherapy, portal films of all four treatment fields were obtained 1 day before treatment. Subsequently, portal films were then obtained at least once a week. Portal films were compared with the simulation films and appropriate changes were made and verified on the next day prior to treatment. A deviation of greater than 0.5 cm or greater was considered to be clincally significant in our analysis. We studied the difference among the types of immobilization and no immobilization by looking at the frequency of movements (overall, and on each of the three axes) that a patient had during the course of his treatment. Using a logistic regression model, the probability of overall and individual directional movement for each group was obtained. In addition, the effects of patient body habitus and pelvic circumference on movement were analyzed.The maximum deviation was 2 cm and the median deviation was 1.2 cm. For each patient, the probability of movement ranged from 0 to 76%, with a mean of 39%. There was no significant difference seen in overall movement with any of the immobilzation devices compared to no immobilization, but there was less vertical (9 vs. 18%; p = 0.03) and AP (6 vs. 15%; p = 0.14) movement with the aquaplast than any other group. However, when examining the lateral direction, the aquaplast had significantly more movement (32 vs. 9%; p0.001). When accounting for body habitus and pelvic circumference, no immobilization device was effective in reducing movement in obese patients or in patients with pelvic circumference greater than 105 cm. The aquaplast group had a significantly increased amount of lateral movement with obesity (42 vs. 23%; p0.05), and with pelvic circumference105 cm (33 vs. 29%; p0.05).There was no significant reduction in overall patient movement noted with any of the immobilization devices compared to no immobilization. The aquaplast group had reduced vertical and AP movement of greater than 0.5 cm. There was significantly more lateral movement with aquaplast appreciated in obese patients or patients with pelvic circumferences greater than 105 cm. The aquaplast immobilization appears to be useful in reducing movement in two very clinicaly important dimensions (AP and vertical). Despite our findings, other immobilization may still be useful especially in the treatment of nonobese patients. It is clear that the optimal immobilization technique and patient positioning are yet to be determined.

Purpose : To document the size and location of the prostate, seminal vesicles, bladder, and rectum throughout the course of external beam radiotherapy. The frequency and range of motion of these organs are quantified. Methods and Materials : Ten patients with localized carcinoma of the prostate had conventional simulation followed immediately by a treatment planning computed tomography scan (TPCT 0 ). Once treatment was initiated, each patient had a weekly CT (TPCT 1-N ) before or after his daily treatment. Anatomical structures from CT were delineated on a computer workstation for analysis. The serial CT sets were spatially registered to the initial scan using image correlation software that brings into congruence the bony pelvis of the different scans. The location of the prostate, seminal vesicles, bladder, and rectum on subsequent scans were compared to TPCT 0 , as well as to each other. Results : Prostate volumes were observed to vary by an average of ±10% during the course of radiation therapy, while the seminal vesicle volumes varied by as much as 100%. Bladder and rectal volumes varied by ±30%. Compared to TPCT 0 , movement of the prostate was demonstrated in all patients. Quantitation of the center-of-mass (CM) showed motion of less than 1 mm in the left-right direction, while motion ranging from 0 to ± 1 cm was observed in the anterior-posterior and superior-inferior directions. The individual standard deviations of these motions varied from approximately 1-5 mm. These variations were correlated to changes in the dimensions of the bladder and rectum. Conclusions : Changes in the location of the prostate, seminal vesicles, and normal tissue volumes during the course of radiation therapy occur and have dosimetric consequences that may impact tumor control and normal tissue complication probabilities. Conformal therapy for prostate cancer will require the incorporation of knowledge of the anatomic relationships of these structures as a function of time. Therefore, these uncertainties must be taken into account when designing treatment plans and in considering dose escalation trials.

Purpose: To develop a method of analyzing rectal surface area irradiated and rectal complications in prostate conformal radiotherapy. Methods and Materials: Dose-surface histograms of the rectum, which state the rectal surface area irradiated to any given dose, were calculated for a group of 27 patients treated with a four-field box technique to a total (tumor minimum) dose ranging from 68 to 70 Gy. Occurrences of rectal toxicities as defined by the Radiation Therapy Oncology Group (RTOG) were recorded and examined in terms of dose and rectal surface area irradiated. For a specified end point of rectal complication, the complication probability was analyzed as a function of dose irradiated to a fixed rectal area, and as a function of area receiving a fixed dose. Lyman’s model of normal tissue complication probability (NTCP) was used to fit the data. Results: The observed occurrences of rectal complications appear to depend on the rectal surface area irradiated to a given dose level. The patient distribution of each toxicity grade exhibits a maximum as a function of percentage surface area irradiated, and the maximum moves to higher values of percentage surface area as the toxicity grade increases. The dependence of the NTCP for the specified end point on dose and percentage surface area irradiated was fitted to Lyman’s NTCP model with a set of parameters. The curvature of the NTCP as a function of the surface area suggests that the rectum is a parallel structured organ. Conclusions: The described method of analyzing rectal surface area irradiated yields interesting insight into understanding rectal complications in prostate conformal radiotherapy. Application of the method to a larger patient data set has the potential to facilitate the construction of a full dose-surface-complication relationship, which would be most useful in guiding clinical practice. Normal tissue complication probability (NTCP), Dose-surface histogram, Conformal radiotherapy, Prostate cancer.

A method of unfolding patient movement between biplane radiographs is presented. This method requires biplane projections of two points with a known separation to determine a translation motion or biplane projections of two or three point pairs to determine a general patient motion. Such a point pair can be provided by two ends of a line source or two radiopaque marks, with a fixed distance between them. The point pairs can be taped on skin, implanted in the target volume, or embedded in the bony structures. After obtaining the point locations on each film, the three-dimensional coordinates of points of interest at the instant of each view can be reconstructed, with an error of less than 0.5 mm. The motion parameters, characterized by a pseudotranslation vector and three axial rotational angles for the region of interest relative to the machine coordinate system, are also provided. In comparing with some other reconstruction algorithms, the point pair method can significantly improve source localization or point reconstruction in the presence of patient motion. Computer simulations and phantom experiments demonstrates the feasibility and accuracy of the method.

A computer-aided method for reconstruction of the source positions of a Fletcher-Suit applicator has been developed. The tandem source positions are determined by digitizing of the tip and two arbitrary points on the shaft from each of two orthogonal simulator films. A colpostat source position is reconstructed by digitizing of a single point on the endcap and three arbitrary points on the cylindrical sidewalls of the colpostat on the films. This computer-aided method considers the true projection geometry and applicator shape and permits localization of the source positions to within a mean error of less than 1 mm. Compared with the conventional method, the new approach (1) is more time efficient because only a few easily identified points are digitized, (2) allows localization when the tandem sources are shifted by a spacer or when colpostat sources are difficult to visualize on the lateral film, and (3) is more accurate than the conventional technique because no manual drawing of source positions on films is involved.

Purpose: With on-line portal imaging devices and image registration tools, the verification of radiation field position prior to each treatment becomes technically feasible. In this paper, we analyze the impact of pre-treatment verification and field position adjustment on target coverage and normal tissue sparing. Methods and Materials: Port films were compared with corresponding simulation films to determine the magnitude of setup variations in patients treated for prostate cancer. From these data, an analytic function was determined between geometric coverage of the target and field margin size. A paradigm for on-line patient repositioning was employed to generate a new relationship between margin and target coverage. Margins were selected for the situations of normal treatment and on-line repositioning to ensure target coverage. Dose-volume histograms were generated for a typical prostate treatment using these margins. Results: On-line repositioning, when setup errors exceed 1 cm, results in a 6 mm reduction in margin, suggesting that 10% of the volume of bladder and rectum may be spared of high dose. Conclusion: The use of on-line imaging and image registration to guide adjustment of patient setup may lead to a reduction in the volume of normal tissues irradiated, and possibly improve the probability of complication-free survival in future treatments.

Intraperitoneal and intracavitary radioimmunotherapy differ from other approaches of radioimmunotherapy in that high activity and dose gradients exist near the solution/tumor interface. Dose to tumor and normal tissue at the interface is a function of depth and is due to three major components: (1) the activity concentration of the administered radiolabeled antibody solution as a function of time within the compartment; (2) the spatial distribution of antibody/radionuclide complex as a function of depth and time as the biomolecules bind to and permeate tumor/normal tissues; and (3) the physical characteristics of the radionuclide in relation to depth of antibody penetration. In this review, the biological and physical aspects of intraperitoneally administered radiolabeled antibodies are discussed, and the state of experimental and calculational studies for this site is described. Areas requiring future investigation are examined, and recommendations are made regarding the type of measurements and calculations which are required for accurate dosimetry.

No definite objective guidelines for selecting patients for pelvic nodal or seminal-vesicle irradiation exist in radical radiotherapy for prostate cancer. Developing such guidelines will aid in decrea

The chronic and acute toxicities associated with conventional radiotherapy of localized prostate cancer are well documented. However, the degree and incidence of toxicities with conformal techniques are not known. Studying side effects associated with modern radiotherapeutic techniques is more important now since there has been a general trend to use computerized tomography-based techniques in recent years; beam's eye view-based conformal techniques are also becoming more commonplace. It is possible that the local disease control can be improved with the delivery of higher doses than currently used. Conformation of the treatment volume to the target volume may facilitate such dose-escalation. However, prior to such dose-escalation, it is important to know the toxicities associated with such techniques with conventional doses.We have compared week-by-week acute toxicities associated with conventional (Group A, 16 patients), computerized tomography-based, manual (Group B, 57 patients) and beam's eye view-based (Group C, 43 patients) techniques during 7 weeks of radiotherapy. Group B and C patients were treated contemporaneously (1988-1990).Acute side effects gradually increased from week 1 through weeks 4-5 and generally declined or plateaued after that. The incidence of acute toxicities was significantly less with the beam's eye view/based technique than with the other two methods. For instance, the percentages of Grade 2 acute genitourinary toxicities for Groups A, B, and C were as follows: Week 1-0, 0, 0; Week 2-6, 0, 0; Week 3-6, 9, 2; Week 4-12, 14, 9; Week 5-35, 14, 9; Week 6-31, 16, 7; Week 7-33, 8, 8, respectively. The p values associated with differences in acute genitourinary toxicities for Weeks 1-7 using chi-square test were 0.072, 0.627, 0.389, 0.538, 0.123, 0.06, and 0.012; the p values for acute gastrointestinal toxicities were 0.512, 0.09, 0.031, 0.031, 0.003,0.0001, and 0.004, respectively. Pairwise comparison (Wilcoxon rank-sum test) showed statistically significant lower acute toxicity in Group C than Group B (e.g., p values, Weeks 1-7 for gastrointestinal toxicity: 0.633, 0.056, 0.010, 0.014,0.0001,0.0001, and0.0001, respectively) in the latter part of the treatment course. No correlation was found between the extent of toxicity and the patient age or the overall treatment time. Also, no correlation was found between the degree of toxicity and the radiation dose and fraction size, within the narrow ranges used (65-70 Gy and 180-200 cGy, respectively). A trend suggesting increased severity of toxicity with increase in the volume of treatment was seen.The findings in this retrospective study need to be confirmed by other prospective studies.

In this work we tabulate and histogram the image segmentation time for beam’s eye view (BEV) treatment planning in our center. The average time needed to generate contours on CT images delineating normal structures and treatment target volumes is calculated using a data base containing over 500 patients’ BEV plans. The average number of contours and total image segmentation time needed for BEV plans in three common treatment sites, namely, head/neck, lung/chest, and prostate, were estimated.

Purpose To characterize and quantify multiorgan respiration-induced motion in the abdomen in liver and pancreatic cancer patients. Methods and Materials Four-dimensional computed tomography scans were acquired for 18 patients treated for abdominal tumors. Contours of multiple abdominal organs were drawn by the radiation oncologist at one respiratory phase; these contours were propagated to other respiratory phases by deformable registration. Three-dimensional organ models were generated from the resulting contours at each phase. Motions of the bounding box and center of mass were extracted and analyzed for the clinical target volume and organs at risk. Results On average, the center of mass motion for liver clinical target volumes was 9.7 mm (SD 5 mm) in the superior–inferior direction, with a range of 3 to 18 mm; for pancreatic tumors, the average was 5 mm (SD 1 mm) m with a range of 3 to 7 mm. Abdominal organs move in unison, but with varying amplitudes. Gating near exhale (T40–T60) reduces the range of motion by a factor of ∼10. Conclusion We have used deformable registration to calculate the trajectories of abdominal organs in four dimensions, based on center of mass and bounding box motion metrics. Our results are compared with previously reported studies. Possible reasons for differences are discussed.

This paper examines the uncertainty in estimating lung motion from external surrogates for lung cancer patients with regular and irregular breathing. 4DCT data sets were analyzed using a template matching algorithm to track the spatial movement of vessel bifurcations in 12 patients. The detected internal movement of features in 3D was retrospectively synchronized with the RPM surrogate signal, and the correlation index R2 and the prediction error were computed. Patients were classified into two groups depending on the presence or not of irregularities in their breathing pattern. Peak-to-peak values of feature motion in the SI direction ranged from 0.8 mm (upper lung) to 25.3 mm (lower lung). Some patients exhibited large motion also in the latero-lateral (10.6 mm) and anterior-posterior (12.2 mm) directions. The median ± quartile of R2 in SI direction was 0.89 ± 0.09. Prediction error values were up to 4.2 mm (95th percentile) with a maximum value of 4.9 mm. Statistical differences between regular and irregular breathers were found for R2, while prediction error depended only on the range of motion. This study is relevant for image guided radiotherapy methods that rely on external surrogates to monitor motion.

Prior to epilepsy surgery, subdural electrodes are often implanted and monitored for a few days to identify the focus of abnormal electrical activity. During the implantation and subsequent brain resection, there may be uncertainty about the exact location of the electrodes with respect to features of brain anatomy such as specific gyral convolutions or lesions. In experiments with a phantom and patients, implanted electrodes were imaged with multiplanar skull radiographs (or CT scans). After retrospective registration with pre-implantation MR data, the electrodes were mapped from these studies onto an MR-derived three-dimensional brain model. The resulting multimodality displays showed the relationship of the electrodes to brain anatomy. In one patient the position of each electrode with respect to a metabolic lesion was also displayed by mapping preimplantation PET data onto the same brain model. This new display of electrode positions may strengthen the interpretation of subdural electrical recordings and thereby reduce uncertainty in planning the resection of epileptic tissue.