Showing papers in "Medical Physics in 2005"

Respiratory correlated cone beam CT

Abstract: A cone beam computed tomography (CBCT) scanner integrated with a linear accelerator is a powerful tool for image guided radiotherapy. Respiratory motion, however, induces artifacts in CBCT, while the respiratory correlated procedures, developed to reduce motion artifacts in axial and helical CT are not suitable for such CBCT scanners. We have developed an alternative respiratory correlated procedure for CBCT and evaluated its performance. This respiratory correlated CBCT procedure consists of retrospective sorting in projection space, yielding subsets of projections that each corresponds to a certain breathing phase. Subsequently, these subsets are reconstructed into a four-dimensional (4D) CBCT dataset. The breathing signal, required for respiratory correlation, was directly extracted from the 2D projection data, removing the need for an additional respiratory monitor system. Due to the reduced number of projections per phase, the contrast-to-noise ratio in a 4D scan reduced by a factor 2.6-3.7 compared to a 3D scan based on all projections. Projection data of a spherical phantom moving with a 3 and 5 s period with and without simulated breathing irregularities were acquired and reconstructed into 3D and 4D CBCT datasets. The positional deviations of the phantoms center of gravity between 4D CBCT and fluoroscopy were small: 0.13 +/- 0.09 mm for the regular motion and 0.39 +/- 0.24 mm for the irregular motion. Motion artifacts, clearly present in the 3D CBCT datasets, were substantially reduced in the 4D datasets, even in the presence of breathing irregularities, such that the shape of the moving structures could be identified more accurately. Moreover, the 4D CBCT dataset provided information on the 3D trajectory of the moving structures, absent in the 3D data. Considerable breathing irregularities, however, substantially reduces the image quality. Data presented for three different lung cancer patients were in line with the results obtained from the phantom study. In conclusion, we have successfully implemented a respiratory correlated CBCT procedure yielding a 4D dataset. With respiratory correlated CBCT on a linear accelerator, the mean position, trajectory, and shape of a moving tumor can be verified just prior to treatment. Such verification reduces respiration induced geometrical uncertainties, enabling safe delivery of 4D radiotherapy such as gated radiotherapy with small margins.

Precise radiochromic film dosimetry using a flat-bed document scanner

Abstract: In this study, a measurement protocol is presented that improves the precision of dose measurements using a flat-bed document scanner in conjunction with two new GafChromic® film models, HS and Prototype A EBT exposed to 6MV photon beams. We established two sources of uncertainties in dose measurements, governed by measurement and calibration curve fit parameters contributions. We have quantitatively assessed the influence of different steps in the protocol on the overall dose measurement uncertainty. Applying the protocol described in this paper on the Agfa Arcus II flat-bed document scanner, the overall one-sigma dose measurement uncertainty for an uniform field amounts to 2% or less for doses above around 0.4Gy in the case of the EBT (Prototype A), and for doses above 5Gy in the case of the HS model GafChromic® film using a region of interest 2×2mm2 in size.

Four-dimensional computed tomography: Image formation and clinical protocol

Abstract: 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.

Assessment of display performance for medical imaging systems: executive summary of AAPM TG18 report.

Abstract: Digital imaging provides an effective means to electronically acquire, archive, distribute, and view medical images. Medical imaging display stations are an integral part of these operations. Therefore, it is vitally important to assure that electronic display devices do not compromise image quality and ultimately patient care. The AAPM Task Group 18 (TG18) recently published guidelines and acceptance criteria for acceptance testing and quality control of medical display devices. This paper is an executive summary of the TG18 report. TG18 guidelines include visual, quantitative, and advanced testing methodologies for primary and secondary class display devices. The characteristics, tested in conjunction with specially designed test patterns (i.e., TG18 patterns), include reflection, geometric distortion, luminance, the spatial and angular dependencies of luminance, resolution, noise, glare, chromaticity, and display artifacts. Geometric distortions are evaluated by linear measurements of the TG18-QC test pattern, which should render distortion coefficients less than 2%/5% for primary/secondary displays, respectively. Reflection measurements include specular and diffuse reflection coefficients from which the maximum allowable ambient lighting is determined such that contrast degradation due to display reflection remains below a 20% limit and the level of ambient luminance (Lamb) does not unduly compromise luminance ratio (LR) and contrast at low luminance levels. Luminance evaluation relies on visual assessment of low contrast features in the TG18-CT and TG18-MP test patterns, or quantitative measurements at 18 distinct luminance levels of the TG18-LN test patterns. The major acceptable criteria for primary/ secondary displays are maximum luminance of greater than 170/100 cd/m2, LR of greater than 250/100, and contrast conformance to that of the grayscale standard display function (GSDF) of better than 10%/20%, respectively. The angular response is tested to ascertain the viewing cone within which contrast conformance to the GSDF is better than 30%/60% and LR is greater than 175/70 for primary/secondary displays, or alternatively, within which the on-axis contrast thresholds of the TG18-CT test pattern remain discernible. The evaluation of luminance spatial uniformity at two distinct luminance levels across the display faceplate using TG18-UNL test patterns should yield nonuniformity coefficients smaller than 30%. The resolution evaluation includes the visual scoring of the CX test target in the TG18-QC or TG18-CX test patterns, which should yield scores greater than 4/6 for primary/secondary displays. Noise evaluation includes visual evaluation of the contrast threshold in the TG18-AFC test pattern, which should yield a minimum of 3/2 targets visible for primary/secondary displays. The guidelines also include methodologies for more quantitative resolution and noise measurements based on MTF and NPS analyses. The display glare test, based on the visibility of the low-contrast targets of the TG18-GV test pattern or the measurement of the glare ratio (GR), is expected to yield scores greater than 3/1 and GRs greater than 400/150 for primary/secondary displays. Chromaticity, measured across a display faceplate or between two display devices, is expected to render a u',v' color separation of less than 0.01 for primary displays. The report offers further descriptions of prior standardization efforts, current display technologies, testing prerequisites, streamlined procedures and timelines, and TG18 test patterns.

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

Abstract: As more pretreatment imaging becomes integrated into the treatment planning process and full three-dimensional image-guidance becomes part of the treatment delivery the need for a deformable image registration technique becomes more apparent. A novel finite element model-based multiorgan deformable image registration method, MORFEUS, has been developed. The basis of this method is twofold: first, individual organ deformation can be accurately modeled by deforming the surface of the organ at one instance into the surface of the organ at another instance and assigning the material properties that allow the internal structures to be accurately deformed into the secondary position and second, multi-organ deformable alignment can be achieved by explicitly defining the deformation of a subset of organs and assigning surface interfaces between organs. The feasibility and accuracy of the method was tested on MR thoracic and abdominal images of healthy volunteers at inhale and exhale. For the thoracic cases, the lungs and external surface were explicitly deformed and the breasts were implicitly deformed based on its relation to the lung and external surface. For the abdominal cases, the liver, spleen, and external surface were explicitly deformed and the stomach and kidneys were implicitly deformed. The average accuracy (average absolute error) of the lung and liver deformation, determined by tracking visible bifurcations, was 0.19 (s.d.: 0.09), 0.28 (s.d.: 0.12) and 0.17 (s.d.: 0.07) cm, in the LR, AP, and IS directions, respectively. The average accuracy of implicitly deformed organs was 0.11 (s.d.: 0.11), 0.13 (s.d.: 0.12), and 0.08 (s.d.: 0.09) cm, in the LR, AP, and IS directions, respectively. The average vector magnitude of the accuracy was 0.44 (s.d.: 0.20) cm for the lung and liver deformation and 0.24 (s.d.: 0.18) cm for the implicitly deformed organs. The two main processes, explicit deformation of the selected organs and finite element analysis calculations, require less than 120 and 495 s, respectively. This platform can facilitate the integration of deformable image registration into online image guidance procedures, dose calculations, and tissue response monitoring as well as performing multi-modality image registration for purposes of treatment planning.

Electrical Impedance Tomography: Methods, History and Applications

Abstract: This article reviews Electrical Impedance Tomography: Methods, History and Applications by David S. Holder . , Bristol, UK, 2005. ISBN 0-7503-0952-0 (hardcover), 456 pp. Price: $120.00.

Four-dimensional radiotherapy planning for DMLC-based respiratory motion tracking

Abstract: Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy. Temporal anatomic changes can occur for many reasons, though the focus of the current investigation is respiration motion for lung tumors. The aim of this study was to develop 4D radiotherapy treatment-planning methodology for DMLC-based respiratory motion tracking. A 4D computed tomography (CT) scan consisting of a series of eight 3D CT image sets acquired at different respiratory phases was used for treatment planning. Deformable image registration was performed to map each CT set from the peak-inhale respiration phase to the CT image sets corresponding to subsequent respiration phases. Deformable registration allows the contours defined on the peak-inhale CT to be automatically transferred to the other respiratory phase CT image sets. Treatment planning was simultaneously performed on each of the eight 3D image sets via automated scripts in which the MLC-defined beam aperture conforms to the PTV (which in this case equaled the GTV due to CT scan length limitations) plus a penumbral margin at each respiratory phase. The dose distribution from each respiratory phase CT image set was mapped back to the peak-inhale CT image set for analysis. The treatment intent of 4D planning is that the radiation beam defined by the DMLC tracks the respiration-induced target motion based on a feedback loop including the respiration signal to a real-time MLC controller. Deformation with respiration was observed for the lung tumor and normal tissues. This deformation was verified by examining the mapping of high contrast objects, such as the lungs and cord, between image sets. For the test case, dosimetric reductions for the cord, heart, and lungs were found for 4D planning compared with 3D planning. 4D radiotherapy planning for DMLC-based respiratory motion tracking is feasible and may offer tumor dose escalation and/or a reduction in treatment-related complications. However, 4D planning requires new planning tools, such as deformable registration and automated treatment planning on multiple CT image sets.

Volume CT with a flat-panel detector on a mobile, isocentric C-arm: Pre-clinical investigation in guidance of minimally invasive surgery

Abstract: A mobile isocentric C-arm (Siemens PowerMobil) has been modified in our laboratory to include a large area flat-panel detector (in place of the x-ray image intensifier), providing multi-mode fluoroscopy and cone-beam computed tomography(CT)imaging capability. This platform represents a promising technology for minimally invasive, image-guided surgical procedures where precision in the placement of interventional tools with respect to bony and soft-tissue structures is critical. The image quality and performance in surgical guidance was investigated in pre-clinical evaluation in image-guided spinal surgery. The control, acquisition, and reconstruction system are described. The reproducibility of geometric calibration, essential to achieving high three-dimensional (3D) image quality, is tested over extended time scales ( 7 months ) and across a broad range in C-arm angulation (up to 45°), quantifying the effect of improper calibration on spatial resolution, soft-tissue visibility, and image artifacts. Phantom studies were performed to investigate the precision of 3D localization (viz., fiber optic probes within a vertebral body) and effect of lateral projection truncation (limited field of view) on soft-tissue detectability in image reconstructions. Pre-clinical investigation was undertaken in a specific spinal procedure (photodynamic therapy of spinal metastases) in five animal subjects (pigs). In each procedure, placement of fiber optic catheters in two vertebrae (L1 and L2) was guided by fluoroscopy and cone-beam CT. Experience across five procedures is reported, focusing on 3D image quality, the effects of respiratory motion, limited field of view, reconstruction filter, and imaging dose. Overall, the intraoperative cone-beam CTimages were sufficient for guidance of needles and catheters with respect to bony anatomy and improved surgical performance and confidence through 3D visualization and verification of transpedicular trajectories and tool placement. Future investigation includes improvement in image quality, particularly regarding x-ray scatter, motion artifacts and field of view, and integration with optical tracking and navigation systems.

Diffuse optical tomography of breast cancer during neoadjuvant chemotherapy: a case study with comparison to MRI.

Abstract: We employ diffuse optical tomography (DOT) to track treatment progress in a female subject presenting with locally advanced invasive carcinoma of the breast during neoadjuvant chemotherapy. Three-dimensional images of total hemoglobin concentration and scattering identified the tumor. Our measurements reveal tumor shrinkage during the course of chemotherapy, in reasonable agreement with magnetic resonance images of the same subject. A decrease in total hemoglobin concentration contrast between tumor and normal tissue was also observed over time. The results demonstrate the potential of DOT for measuring physiological parameters of breast lesions during chemotherapy.

A simple, direct method for x-ray scatter estimation and correction in digital radiography and cone-beam CT.

Abstract: X-ray scatter poses a significant limitation to image quality in cone-beam CT (CBCT), resulting in contrast reduction, image artifacts, and lack of CT number accuracy. We report the performance of a simple scatter correction method in which scatter fluence is estimated directly in each projection from pixel values near the edge of the detector behind the collimator leaves. The algorithm operates on the simple assumption that signal in the collimator shadow is attributable to x-ray scatter, and the 2D scatter fluence is estimated by interpolating between pixel values measured along the top and bottom edges of the detector behind the collimator leaves. The resulting scatter fluence estimate is subtracted from each projection to yield an estimate of the primary-only images for CBCT reconstruction. Performance was investigated in phantom experiments on an experimental CBCT bench-top, and the effect on image quality was demonstrated in patient images (head, abdomen, and pelvis sites) obtained on a preclinical system for CBCT-guided radiation therapy. The algorithm provides significant reduction in scatter artifacts without compromise in contrast-to-noise ratio (CNR). For example, in a head phantom, cupping artifact was essentially eliminated, CT number accuracy was restored to within 3%, and CNR (breast-to-water) was improved by up to 50%. Similarly in a body phantom, cupping artifact was reduced by at least a factor of 2 without loss in CNR. Patient images demonstrate significantly increased uniformity, accuracy, and contrast, with an overall improvement in image quality in all sites investigated. Qualitative evaluation illustrates that soft-tissue structures that are otherwise undetectable are clearly delineated in scatter-corrected reconstructions. Since scatter is estimated directly in each projection, the algorithm is robust with respect to system geometry, patient size and heterogeneity, patient motion, etc. Operating without prior information, analytical modeling, or Monte Carlo, the technique is easily incorporated as a preprocessing step in CBCT reconstruction to provide significant scatter reduction.

Accurate technique for complete geometric calibration of cone-beam computed tomography systems.

Abstract: Cone-beam computed tomography systems have been developed to provide in situ imaging for the purpose of guiding radiation therapy. Clinical systems have been constructed using this approach, a clinical linear accelerator (Elekta Synergy RP) and an iso-centric C-arm. Geometric calibration involves the estimation of a set of parameters that describes the geometry of such systems, and is essential for accurate image reconstruction. We have developed a general analytic algorithm and corresponding calibration phantom for estimating these geometric parameters in cone-beam computed tomography (CT) systems. The performance of the calibration algorithm is evaluated and its application is discussed. The algorithm makes use of a calibration phantom to estimate the geometric parameters of the system. The phantom consists of 24 steel ball bearings (BBs) in a known geometry. Twelve BBs are spaced evenly at 30 deg in two plane-parallel circles separated by a given distance along the tube axis. The detector (e.g., a flat panel detector) is assumed to have no spatial distortion. The method estimates geometric parameters including the position of the x-ray source, position, and rotation of the detector, and gantry angle, and can describe complex source-detector trajectories. The accuracy and sensitivity of the calibration algorithm was analyzed. The calibration algorithm estimates geometric parameters in a high level of accuracy such that the quality of CT reconstruction is not degraded by the error of estimation. Sensitivity analysis shows uncertainty of 0.01 degrees (around beam direction) to 0.3 degrees (normal to the beam direction) in rotation, and 0.2 mm (orthogonal to the beam direction) to 4.9 mm (beam direction) in position for the medical linear accelerator geometry. Experimental measurements using a laboratory bench Cone-beam CT system of known geometry demonstrate the sensitivity of the method in detecting small changes in the imaging geometry with an uncertainty of 0.1 mm in transverse and vertical (perpendicular to the beam direction) and 1.0 mm in the longitudinal (beam axis) directions. The calibration algorithm was compared to a previously reported method, which uses one ball bearing at the isocenter of the system, to investigate the impact of more precise calibration on the image quality of cone-beam CT reconstruction. A thin steel wire located inside the calibration phantom was imaged on the conebeam CT lab bench with and without perturbations in source and detector position during the scan. The described calibration method improved the quality of the image and the geometric accuracy of the object reconstructed, improving the full width at half maximum of the wire by 27.5% and increasing contrast of the wire by 52.8%. The proposed method is not limited to the geometric calibration of cone-beam CT systems but can be used for many other systems, which consist of one or more point sources and area detectors such as calibration of megavoltage (MV) treatment system (focal spot movement during the beam delivery, MV source trajectory versus gantry angle, the axis of collimator rotation, and couch motion), cross calibration between Kilovolt imaging and MV treatment system, and cross calibration between multiple imaging systems. Using the complete information of the system geometry, it was demonstrated that high image quality in CT reconstructions is possible even in systems with large geometric nonidealities.

Image reconstruction and image quality evaluation for a 64-slice CT scanner with z-flying focal spot

Abstract: We present a theoretical overview and a performance evaluation of a novel z-sampling technique for multidetector row CT (MDCT), relying on a periodic motion of the focal spot in the longitudinal direction (z-flying focal spot) to double the number of simultaneously acquired slices. The z-flying focal spot technique has been implemented in a recently introduced MDCT scanner. Using 32 x 0.6 mm collimation, this scanner acquires 64 overlapping 0.6 mm slices per rotation in its spiral (helical) mode of operation, with the goal of improved longitudinal resolution and reduction of spiral artifacts. The longitudinal sampling distance at isocenter is 0.3 mm. We discuss in detail the impact of the z-flying focal spot technique on image reconstruction. We present measurements of spiral slice sensitivity profiles (SSPs) and of longitudinal resolution, both in the isocenter and off-center. We evaluate the pitch dependence of the image noise measured in a centered 20 cm water phantom. To investigate spiral image quality we present images of an anthropomorphic thorax phantom and patient scans. The full width at half maximum (FWHM) of the spiral SSPs shows only minor variations as a function of the pitch, measured values differ by less than 0.15 mm from the nominal values 0.6, 0.75, 1, 1.5, and 2 mm. The measured FWHM of the smallest slice ranges between 0.66 and 0.68 mm at isocenter, except for pitch 0.55 (0.72 mm). In a centered z-resolution phantom, bar patterns up to 15 lp/cm can be visualized independent of the pitch, corresponding to 0.33 mm longitudinal resolution. 100 mm off-center, bar patterns up to 14 lp/cm are visible, corresponding to an object size of 0.36 mm that can be resolved in the z direction. Image noise for constant effective mAs is almost independent of the pitch. Measured values show a variation of less than 7% as a function of the pitch, which demonstrates correct utilization of the applied radiation dose at any pitch. The product of image noise and square root of the slice width (FWHM of the respective SSP) is the same constant for all slices except 0.6 mm. For the thinnest slice, relative image noise is increased by 17%. Spiral windmill-type artifacts are effectively suppressed with the z-flying focal spot technique, which has the potential to maintain a low artifact level up to pitch 1.5, in this way increasing the maximum volume coverage speed that can be clinically used.

Performance characterization of megavoltage computed tomography imaging on a helical tomotherapy unit.

Abstract: Helical tomotherapy is an innovative means of delivering IGRT and IMRT using a device that combines features of a linear accelerator and a helical computed tomography (CT) scanner. The HI-ART II can generate CT images from the same megavoltage x-ray beam it uses for treatment. These megavoltage CT (MVCT) images offer verification of the patient position prior to and potentially during radiation therapy. Since the unit uses the actual treatment beam as the x-ray source for image acquisition, no surrogate telemetry systems are required to register image space to treatment space. The disadvantage to using the treatment beam for imaging, however, is that the physics of radiation interactions in the megavoltage energy range may force compromises between the dose delivered and the image quality in comparison to diagnostic CT scanners. The performance of the system is therefore characterized in terms of objective measures of noise, uniformity, contrast, and spatial resolution as a function of the dose delivered by the MVCT beam. The uniformity and spatial resolutions of MVCT images generated by the HI-ART II are comparable to that of diagnostic CT images. Furthermore, the MVCT scan contrast is linear with respect to the electron density of material imaged. MVCT images do not have the same performance characteristics as state-of-the art diagnostic CT scanners when one objectively examines noise and low-contrast resolution. These inferior results may be explained, at least partially, by the low doses delivered by our unit; the dose is 1.1 cGy in a 20 cm diameter cylindrical phantom. In spite of the poorer low-contrast resolution, these relatively low-dose MVCT scans provide sufficient contrast to delineate many soft-tissue structures. Hence, these images are useful not only for verifying the patient's position at the time of therapy, but they are also sufficient for delineating many anatomic structures. In conjunction with the ability to recalculate radiotherapy doses on these images, this enables dose guidance as well as image guidance of radiotherapy treatments.

A phantom evaluation of a stereo-vision surface imaging system for radiotherapy patient setup.

Abstract: 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.

Feasibility study of helical tomotherapy for total body or total marrow irradiation.

Abstract: Total body radiation (TBI) has been used for many years as a preconditioning agent before bone marrow transplantation. Many side effects still plague its use. We investigated the planning and delivery of total body irradiation (TBI) and selective total marrow irradiation (TMI) and a reduced radiation dose to sensitive structures using image-guided helical tomotherapy. To assess the feasibility of using helical tomotherapy, (A) we studied variations in pitch, field width, and modulation factor on total body and total marrow helical tomotherapy treatments. We varied these parameters to provide a uniform dose along with a treatment times similar to conventional TBI (15-30 min). (B) We also investigated limited (head, chest, and pelvis) megavoltage CT (MVCT) scanning for the dimensional pretreatment setup verification rather than total body MVCT scanning to shorten the overall treatment time per treatment fraction. (C) We placed thermoluminescent detectors (TLDs) inside a Rando phantom to measure the dose at seven anatomical sites, including the lungs. A simulated TBI treatment showed homogeneous dose coverage (+/-10%) to the whole body. Doses to the sensitive organs were reduced by 35%-70% of the target dose. TLD measurements on Rando showed an accurate dose delivery (+/-7%) to the target and critical organs. In the TMI study, the dose was delivered conformally to the bone marrow only. The TBI and TMI treatment delivery time was reduced (by 50%) by increasing the field width from 2.5 to 5.0 cm in the inferior-superior direction. A limited MVCT reduced the target localization time 60% compared to whole body MVCT. MVCT image-guided helical tomotherapy offers a novel method to deliver a precise, homogeneous radiation dose to the whole body target while reducing the dose significantly to all critical organs. A judicious selection of pitch, modulation factor, and field size is required to produce a homogeneous dose distribution along with an acceptable treatment time. In addition, conformal radiation to the bone marrow appears feasible in an external radiation treatment using image-guided helical tomotherapy.

Spectral discrimination of Čerenkov radiation in scintillating dosimeters

Abstract: Radiation therapy accelerators require highly accurate dose deposition and the output must be monitored frequently and regularly. Ionization chambers are the primary tool for this control, but their size, their high voltage needed, and the correction needed for electrons make them unsuitable for use during patient treatment. We have developed a small (1-mm-diam and 1-mm-long active part), flexible, and water-equivalent dosimeter. It is suitable for photon and electron beams without corrections, and performs on line dose measurements. This detector is based on only one scintillating fiber and a CCD camera. A new signal processing is used to remove the effect of Cerenkov radiation background, which only requires a preliminary calibration. Central-axis depth-dose distribution comparisons have been achieved with standard ionization chambers, over a range from 8 to 25 MV photons and from 6 to 21 MeV electrons in order to validate this calibration. Results show a very good agreement, with less than 1% difference between the two detectors.

Technical and dosimetric aspects of respiratory gating using a pressure-sensor motion monitoring system.

Abstract: This work introduces a gating technique that uses 4DCT to determine gating parameters and to plan gated treatment, and employs a Siemens linear accelerator to deliver the gated treatment. Because of technology incompatibility, the 4DCT scanner (LightSpeed, GE) and the Siemens accelerator require two different motion-monitoring systems. The motion monitoring system (AZ-773V, Anzai Med.) used for the gated delivery utilizes a pressure sensor to detect the external respiratory motion (pressure change) in real time. Another system (RPM, Varian) used for the 4DCT scanner (LightSpeed, GE) is based on an infrared camera to detect motion of external markers. These two motion monitoring systems (RPM and Anzai systems) were found to correlate well with each other. The depth doses and profile measured for gated delivery (with a duty cycle of 25% or 50%) were found to agree within 1.0% with those measured for ungated delivery, indicating that gating did not significantly alter beam characteristics. The measurement verified also that the MU linearity and beam output remained unchanged (within 0.3%). A practical method of using 4DCT to plan a gated treatment was developed. The duty cycle for either phase or amplitude gating can be determined based on 4DCT with consideration of set-up error and delivery efficiency. The close-loop measurement involving the entire gating process (imaging, planning, and delivery) showed that the measured isodose distributions agreed with those intended, validating the accuracy and reliability of the gating technique. Based these observations, we conclude that the gating technique introduced in this work, integrating Siemens linear accelerator and Anzai pressure sensor device with GE/Varian RPM 4DCT, is reliable and effective, and it can be used clinically to account for respiratory motion during radiation therapy.

Energy dependence of response of new high sensitivity radiochromic films for megavoltage and kilovoltage radiation energies

Abstract: The purpose of this paper is to evaluate the energy dependence of the response of two new high sensitivity models of radiochromic films EBT and XR-QA. We determined the dose response curves of these films for four different radiation sources, namely, 6 MV photon beams (6 MVX), Ir-192, I-125, and Pd-103. The first type (EBT) is designed for intensity modulated radiation therapy (IMRT) dosimetry, and the second type (XR-QA) is designed for kilovoltage dosimetry. All films were scanned using red (665 nm) and green (520 nm) light sources in a charge-coupled device-based densitometer. The dose response curves [net optical density (NOD) versus dose] were plotted and compared for different radiation energies and light sources. Contrary to the early GAFCHROMIC film types (such as models XR, HS, MD55-2, and HD810), the net optical densities of both EBT and XR-QA were higher with a green (520 nm) than those with a red (665 nm) light source due to the different absorption spectrum of the new radiochromic emulsion. Both film types yield measurable optical densities for doses below 2 Gy. EBT film response is nearly independent of radiation energy, within the uncertainty of measurement. The NOD values of EBT film at 1 and 2 Gy are 0.13 and 0.25 for green, and 0.1 and 0.17 for red, respectively. In contrast, the XR-QA film sensitivity varies with radiation energy. The doses required to produce NOD of 0.5 are 6.9, 5.4, 0.7, and 0.9 Gy with green light and 19, 13, 1.7, and 1.5 Gy with red light, for 6 MVX, Ir-192, I -125, and Pd-103, respectively. EBT film was found to have minimal photon energy dependence of response for the energies tested and is suitable for dosimetry of radiation with a wide energy spectrum, including primary and scattered radiation. XR-QA film is promising for kilovoltage sources with a narrow energy spectra. The new high sensitivity radiochromic films are promising tools in radiation dosimetry.

Measurement accuracy and cerenkov removal for high performance, high spatial resolution scintillation dosimetry.

Abstract: With highly conformal radiation therapy techniques such as intensity-modulated radiation therapy, radiosurgery, and tomotherapy becoming more common in clinical practice, the use of these narrow beams requires a higher level of precision in quality assurance and dosimetry. Plastic scintillators with their water equivalence, energy independence, and dose rate linearity have been shown to possess excellent qualities that suit the most complex and demanding radiation therapy treatment plans. The primary disadvantage of plastic scintillators is the presence of Cerenkov radiation generated in the light guide, which results in an undesired stem effect. Several techniques have been proposed to minimize this effect. In this study, we compared three such techniques-background subtraction, simple filtering, and chromatic removal-in terms of reproducibility and dose accuracy as gauges of their ability to remove the Cerenkov stem effect from the dose signal. The dosimeter used in this study comprised a 6-mm(3) plastic scintillating fiber probe, an optical fiber, and a color charge-coupled device camera. The whole system was shown to be linear and the total light collected by the camera was reproducible to within 0.31% for 5-s integration time. Background subtraction and chromatic removal were both found to be suitable for precise dose evaluation, with average absolute dose discrepancies of 0.52% and 0.67%, respectively, from ion chamber values. Background subtraction required two optical fibers, but chromatic removal used only one, thereby preventing possible measurement artifacts when a strong dose gradient was perpendicular to the optical fiber. Our findings showed that a plastic scintillation dosimeter could be made free of the effect of Cerenkov radiation.

Density resolution of proton computed tomography.

Abstract: Conformal proton radiation therapy requires accurate prediction of the Bragg peak position. Protons may be more suitable than conventional x rays for this task since the relative electron density distribution can be measured directly with proton computed tomography (CT). However, proton CT has its own limitations, which need to be carefully studied before this technique can be introduced into routine clinical practice. In this work, we have used analytical relationships as well as the Monte Carlo simulation tool GEANT4 to study the principal resolution limits of proton CT. The noise level observed in proton CT images of a cylindrical water phantom with embedded tissue-equivalent density inhomogeneities, which were generated based on GEANT4 simulations, compared well with predictions based on Tschalar's theory of energy loss straggling. The relationship between phantom thickness, initial energy, and the relative electron density resolution was systematically investigated to estimate the proton dose needed to obtain a given density resolution. We show that a reasonable density resolution can be achieved with a relatively small dose, which is comparable to or even lower than that of x-ray CT.

Development of ultrasound tomography for breast imaging: Technical assessment

Abstract: Ultrasound imaging is widely used in medicine because of its benign characteristics and real-time capabilities. Physics theory suggests that the application of tomographic techniques may allow ultrasound imaging to reach its full potential as a diagnostic tool allowing it to compete with other tomographic modalities such as x-ray computer tomography, and MRI. This paper describes the construction and use of a prototype tomographic scanner and reports on the feasibility of implementing tomographic theory in practice and the potential of ultrasound (US) tomography in diagnostic imaging. Data were collected with the prototype by scanning two types of phantoms and a cadaveric breast. A specialized suite of algorithms was developed and utilized to construct images of reflectivity and sound speed from the phantom data. The basic results can be summarized as follows. (i) A fast, clinically relevant US tomography scanner can be built using existing technology. (ii) The spatial resolution, deduced from images of reflectivity, is 0.4 mm. The demonstrated 10 cm depth-of-field is superior to that of conventional ultrasound and the image contrast is improved through the reduction of speckle noise and overall lowering of the noise floor. (iii) Images of acoustic properties such as sound speed suggest that it is possible to measure variations in the sound speed of 5 m/s. An apparent correlation with x-ray attenuation suggests that the sound speed can be used to discriminate between various types of soft tissue. (iv) Ultrasound tomography has the potential to improve diagnostic imaging in relation to breast cancer detection.

The stability of mechanical calibration for a kV cone beam computed tomography system integrated with linear accelerator.

Abstract: The geometric accuracy and precision of an image-guided treatment system were assessed. Image guidance is performed using an x-ray volume imaging (XVI) system integrated with a linear accelerator and treatment planning system. Using an amorphous silicon detector and x-ray tube, volumetric computed tomography images are reconstructed from kilovoltage radiographs by filtered backprojection. Image fusion and assessment of geometric targeting are supported by the treatment planning system. To assess the limiting accuracy and precision of image-guided treatment delivery, a rigid spherical target embedded in an opaque phantom was subjected to 21 treatment sessions over a three-month period. For each session, a volumetric data set was acquired and loaded directly into an active treatment planning session. Image fusion was used to ascertain the couch correction required to position the target at the prescribed iso-center. Corrections were validated independently using megavoltage electronic portal imaging to record the target position with respect to symmetric treatment beam apertures. An initial calibration cycle followed by repeated image-guidance sessions demonstrated the XVI system could be used to relocate an unambiguous object to within less than 1 mm of the prescribed location. Treatment could then proceed within the mechanical accuracy and precision of the delivery system. The calibrationmore » procedure maintained excellent spatial resolution and delivery precision over the duration of this study, while the linear accelerator was in routine clinical use. Based on these results, the mechanical accuracy and precision of the system are ideal for supporting high-precision localization and treatment of soft-tissue targets.« less

The application of the sinusoidal model to lung cancer patient respiratory motion.

Abstract: Accurate modeling of the respiratory cycle is important to account for the effect of organ motion on dose calculation for lung cancer patients. The aim of this study is to evaluate the accuracy of a respiratory model for lung cancer patients. Lujan et al. [Med. Phys. 26(5), 715-720 (1999)] proposed a model, which became widely used, to describe organ motion due to respiration. This model assumes that the parameters do not vary between and within breathing cycles. In this study, first, the correlation of respiratory motion traces with the model f(t) as a function of the parameter n (n = 1, 2, 3) was undertaken for each breathing cycle from 331 four-minute respiratory traces acquired from 24 lung cancer patients using three breathing types: free breathing, audio instruction, and audio-visual biofeedback. Because cos2 and cos4 had similar correlation coefficients, and cos2 and cos1 have a trigonometric relationship, for simplicity, the cos1 value was consequently used for further analysis in which the variations in mean position (z0), amplitude of motion (b) and period (tau) with and without biofeedback or instructions were investigated. For all breathing types, the parameter values, mean position (z0), amplitude of motion (b), and period (tau) exhibited significant cycle-to-cycle variations. Audio-visual biofeedback showed the least variations for all three parameters (z0, b, and tau). It was found that mean position (z0) could be approximated with a normal distribution, and the amplitude of motion (b) and period (tau) could be approximated with log normal distributions. The overall probability density function (pdf) of f(t) for each of the three breathing types was fitted with three models: normal, bimodal, and the pdf of a simple harmonic oscillator. It was found that the normal and the bimodal models represented the overall respiratory motion pdfs with correlation values from 0.95 to 0.99, whereas the range of the simple harmonic oscillator pdf correlation values was 0.71 to 0.81. This study demonstrates that the pdfs of mean position (z0), amplitude of motion (b), and period (tau) can be used for sampling to obtain more realistic respiratory traces. The overall standard deviations of respiratory motion were 0.48, 0.57, and 0.55 cm for free breathing, audio instruction, and audio-visual biofeedback, respectively.

Penalized-likelihood sinogram smoothing for low-dose CT

Abstract: We have developed a sinogram smoothing approach for low-dose computed tomography (CT) that seeks to estimate the line integrals needed for reconstruction from the noisy measurements by maximizing a penalized-likelihood objective function. The maximization is performed by an algorithm derived by use of the separable paraboloidal surrogates framework. The approach overcomes some of the computational limitations of a previously proposed spline-based penalized-likelihood sinogram smoothing approach, and it is found to yield better resolution-variance tradeoffs than this spline-based approach as well an existing adaptive filtering approach. Such sinogram smoothing approaches could be valuable when applied to the low-dose data acquired in CT screening exams, such as those being considered for lung-nodule detection.

Validation of GATE Monte Carlo simulations of the GE Advance/Discovery LS PET scanners

Abstract: The recently developed GATE (GEANT4 application for tomographic emission) Monte Carlo package, designed to simulate positron emission tomography (PET) and single photon emission computed tomography (SPECT) scanners, provides the ability to model and account for the effects of photon noncollinearity, off-axis detector penetration, detector size and response, positron range, photon scatter, and patient motion on the resolution and quality of PET images. The objective of this study is to validate a model within GATE of the General Electric (GE) Advance/Discovery Light Speed (LS) PET scanner. Our three-dimensional PET simulation model of the scanner consists of 12 096 detectors grouped into blocks, which are grouped into modules as per the vendor's specifications. The GATE results are compared to experimental data obtained in accordance with the National Electrical Manufactures Association/Society of Nuclear Medicine (NEMA/SNM), NEMA NU 2-1994, and NEMA NU 2-2001 protocols. The respective phantoms are also accurately modeled thus allowing us to simulate the sensitivity, scatter fraction, count rate performance, and spatial resolution. In-house software was developed to produce and analyze sinograms from the simulated data. With our model of the GE Advance/Discovery LS PET scanner, the ratio of the sensitivities with sources radially offset 0 and 10 cm from the scanner's main axis are reproduced to within 1% of measurements. Similarly, the simulated scatter fraction for the NEMA NU 2-2001 phantom agrees to within less than 3% of measured values (the measured scatter fractions are 44.8% and 40.9 +/- 1.4% and the simulated scatter fraction is 43.5 +/- 0.3%). The simulated count rate curves were made to match the experimental curves by using deadtimes as fit parameters. This resulted in deadtime values of 625 and 332 ns at the Block and Coincidence levels, respectively. The experimental peak true count rate of 139.0 kcps and the peak activity concentration of 21.5 kBq/cc were matched by the simulated results to within 0.5% and 0.1% respectively. The simulated count rate curves also resulted in a peak NECR of 35.2 kcps at 10.8 kBq/cc compared to 37.6 kcps at 10.0 kBq/cc from averaged experimental values. The spatial resolution of the simulated scanner matched the experimental results to within 0.2 mm.

Brachytherapy needle deflection evaluation and correction

Abstract: In prostate brachytherapy, an 18-gauge needle is used to implant radioactive seeds. This thin needle can be deflected from the preplanned trajectory in the prostate, potentially resulting in a suboptimum dose pattern and at times requiring repeated needle insertion to achieve optimal dosimetry. In this paper, we report on the evaluation of brachytherapy needle deflection and bending in test phantoms and two approaches to overcome the problem. First we tested the relationship between needle deflection and insertion depth as well as whether needle bending occurred. Targeting accuracy was tested by inserting a brachytherapy needle to target 16 points in chicken tissue phantoms. By implanting dummy seeds into chicken tissue phantoms under 3D ultrasound guidance, the overall accuracy of seed implantation was determined. We evaluated methods to overcome brachytherapy needle deflection with three different insertion methods: constant orientation, constant rotation, and orientation reversal at half of the insertion depth. Our results showed that needle deflection is linear with needle insertion depth, and that no noticeable bending occurs with needle insertion into the tissue and agar phantoms. A 3D principal component analysis was performed to obtain the population distribution of needle tip and seed position relative to the target positions. Our resultsmore » showed that with the constant orientation insertion method, the mean needle targeting error was 2.8 mm and the mean seed implantation error was 2.9 mm. Using the constant rotation and orientation reversal at half insertion depth methods, the deflection error was reduced. The mean needle targeting errors were 0.8 and 1.2 mm for the constant rotation and orientation reversal methods, respectively, and the seed implantation errors were 0.9 and 1.5 mm for constant rotation insertion and orientation reversal methods, respectively.« less

Design and application of an assessment protocol for electromagnetic tracking systems.

Abstract: This paper defines a simple protocol for competitive and quantified evaluation of electromagnetic tracking systems such as the NDI Aurora (A) and Ascension microBIRD with dipole transmitter (B). It establishes new methods and a new phantom design which assesses the reproducibility and allows comparability with different tracking systems in a consistent environment. A machined base plate was designed and manufactured in which a 50 mm grid of holes was precisely drilled for position measurements. In the center a circle of 32 equispaced holes enables the accurate measurement of rotation. The sensors can be clamped in a small mount which fits into pairs of grid holes on the base plate. Relative positional/orientational errors are found by subtracting the known distances/rotations between the machined locations from the differences of the mean observed positions/rotation. To measure the influence of metallic objects we inserted rods made of steel (SST 303, SST 416), aluminum, and bronze into the sensitive volume between sensor and emitter. We calculated the fiducial registration error and fiducial location error with a standard stylus calibration for both tracking systems and assessed two different methods of stylus calibration.The positional jitter amounted to 0.14 mm(A) and 0.08 mm(B). A relative positional error of 0.96 mm ± 0.68 mm , range −0.06 mm; 2.23 mm(A) and 1.14 mm ± 0.78 mm , range −3.72 mm; 1.57 mm(B) for a given distance of 50 mm was found. The relative rotation error was found to be 0.51° (A)/0.04° (B). The most relevant distortion caused by metallic objects results from SST 416. The maximum error 4.2 mm ( A ) ∕ ⩾ 100 mm ( B ) occurs when the rod is close to the sensor(20 mm). While (B) is more sensitive with respect to metallic objects, (A) is less accurate concerning orientation measurements. (B) showed a systematic error when distances are calculated.

Comparison of helical and cine acquisitions for 4D-CT imaging with multislice CT

Abstract: We proposed a data sufficiency condition (DSC) for four-dimensional-CT (4D-CT) imaging on a multislice CT scanner, designed a pitch factor for a helical 4D-CT, and compared the acquisition time, slice sensitivity profile (SSP), effective dose, ability to cope with an irregular breathing cycle, and gating technique (retrospective or prospective) of the helical 4D-CT and the cine 4D-CT on the General Electric (GE) LightSpeed RT (4-slice), Plus (4-slice), Ultra (8-slice) and 16 (16-slice) multislice CT scanners. To satisfy the DSC, a helical or cine 4D-CT acquisition has to collect data at each location for the duration of a breathing cycle plus the duration of data acquisition for an image reconstruction. The conditions for the comparison were 20 cm coverage in the cranial-caudal direction, a 4 s breathing cycle, and half-scan reconstruction. We found that the helical 4D-CT has the advantage of a shorter scan time that is 10% shorter than that of the cine 4D-CT, and the disadvantages of 1.8 times broadening of SSP and requires an additional breathing cycle of scanning to ensure an adequate sampling at the start and end locations. The cine 4D-CT has the advantages of maintaining the same SSP as slice collimation (e.g., 8 x 2.5 mm slice collimation generates 2.5 mm SSP in the cine 4D-CT as opposed to 4.5 mm in the helical 4D-CT) and a lower dose by 4% on the 8- and 16-slice systems, and 8% on the 4-slice system. The advantage of faster scanning in the helical 4D-CT will diminish if a repeat scan at the location of a breathing irregularity becomes necessary. The cine 4D-CT performs better than the helical 4D-CT in the repeat scan because it can scan faster and is more dose efficient.

Radiation Physics for Medical Physicists

Abstract: This article reviews Radiation Physics for Medical Physicists by Ervin B. Podgorsak 437+xxi pp. , New York, 2005. $139. ISBN 3-540-25041-7.

On using an adaptive neural network to predict lung tumor motion during respiration for radiotherapy applications.

Abstract: In this study we address the problem of predicting the position of a moving lung tumor during respiration on the basis of external breathing signals--a technique used for beam gating, tracking, and other dynamic motion management techniques in radiation therapy. We demonstrate the use of neural network filters to correlate tumor position with external surrogate markers while simultaneously predicting the motion ahead in time, for situations in which neither the breathing pattern nor the correlation between moving anatomical elements is constant in time. One pancreatic cancer patient and two lung cancer patients with mid/upper lobe tumors were fluoroscopically imaged to observe tumor motion synchronously with the movement of external chest markers during free breathing. The external marker position was provided as input to a feed-forward neural network that correlated the marker and tumor movement to predict the tumor position up to 800 ms in advance. The predicted tumor position was compared to its observed position to establish the accuracy with which the filter could dynamically track tumor motion under nonstationary conditions. These results were compared to simplified linear versions of the filter. The two lung cancer patients exhibited complex respiratory behavior in which the correlation between surrogate marker and tumor position changed with each cycle of breathing. By automatically and continuously adjusting its parameters to the observations, the neural network achieved better tracking accuracy than the fixed and adaptive linear filters. Variability and instability in human respiration complicate the task of predicting tumor position from surrogate breathing signals. Our results show that adaptive signal-processing filters can provide more accurate tumor position estimates than simpler stationary filters when presented with nonstationary breathing motion.