Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Image Processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Anatomical noise in contrast-enhanced digital mammography. Part II. Dual-energy imaging

Purpose: Dual-energy (DE) contrast-enhanced digital mammography (CEDM) uses an iodinated contrast agent in combination with digital mammography (DM) to evaluate lesions on the basis of tumor angiogenesis. In DE imaging, low-energy (LE) and high-energy (HE) images are acquired after contrast administration and their logarithms are subtracted to cancel the appearance of normal breast tissue. Often there is incomplete signal cancellation in the subtracted images, creating a background “clutter” that can impair lesion detection. This is the second component of a two-part report on anatomical noise in CEDM. In Part I the authors characterized the anatomical noise for single-energy (SE) temporal subtraction CEDM by a power law, with model parameters α and β. In this work the authors quantify the anatomical noise in DE CEDM clinical images and compare this with the noise in SE CEDM. The influence on the anatomical noise of the presence of iodine in the breast, the timing of imaging postcontrast administration, and the x-ray energy used for acquisition are each evaluated.Methods: The power law parameters, α and β, were measured from unprocessed LE and HE images and from DE subtracted images to quantify the anatomical noise. A total of 98 DE CEDM cases acquired in a previous clinical pilot study were assessed. Conventional DM images from 75 of the women were evaluated for comparison with DE CEDM. The influence of the imaging technique on anatomical noise was determined from an analysis of differences between the power law parameters as measured in DM, LE, HE, and DE subtracted images for each subject.Results: In DE CEDM, weighted image subtraction lowers β to about 1.1 from 3.2 and 3.1 in LE and HE unprocessed images, respectively. The presence of iodine has a small but significant effect in LE images, reducing β by about 0.07 compared to DM, with α unchanged. Increasing the x-ray energy, from that typical in DM to a HE beam, significantly decreases α by about 2 × 10{sup −5} mm{sup 2}, and lowers β by about 0.14 compared to LE images. A comparison of SE and DE CEDM at 4 min postcontrast shows equivalent power law parameters in unprocessed images, and lower α and β by about 3 × 10{sup −5} mm{sup 2} and 0.50, respectively, in DE versus SE subtracted images.Conclusions: Image subtraction in both SE and DE CEDM reduces β by over a factor of 2, while maintaining α below that in DM. Given the equivalent α between SE and DE unprocessed CEDM images, and the smaller anatomical noise in the DE subtracted images, the DE approach may have an advantage over SE CEDM. It will be necessary to test this potential advantage in future lesion detectability experiments, which account for realistic lesion signals. The authors' results suggest that LE images could be used in place of DM images in CEDM exam interpretation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Calibration of time-of-flight cameras for accurate intraoperative surface reconstruction

Purpose: In image-guided surgery (IGS) intraoperative image acquisition of tissue shape, motion, and morphology is one of the main challenges. Recently, time-of-flight (ToF) cameras have emerged as a new means for fast range image acquisition that can be used for multimodal registration of the patient anatomy during surgery. The major drawbacks of ToF cameras are systematic errors in the image acquisition technique that compromise the quality of the measured range images. In this paper, we propose a calibration concept that, for the first time, accounts for all known systematic errors affecting the quality of ToF range images. Laboratory and in vitro experiments assess its performance in the context of IGS.Methods: For calibration the camera-related error sources depending on the sensor, the sensor temperature and the set integration time are corrected first, followed by the scene-specific errors, which are modeled as function of the measured distance, the amplitude and the radial distance to the principal point of the camera. Accounting for the high accuracy demands in IGS, we use a custom-made calibration device to provide reference distance data, the cameras are calibrated too. To evaluate the mitigation of the error, the remaining residual error after ToF depth calibration was compared with that arising from using the manufacturer routines for several state-of-the-art ToF cameras. The accuracy of reconstructed ToF surfaces was investigated after multimodal registration with computed tomography (CT) data of liver models by assessment of the target registration error (TRE) of markers introduced in the livers.Results: For the inspected distance range of up to 2 m, our calibration approach yielded a mean residual error to reference data ranging from 1.5 ± 4.3 mm for the best camera to 7.2 ± 11.0 mm. When compared to the data obtained from the manufacturer routines, the residual error was reduced by at least 78% from worst calibration result to most accurate manufacturer data. After registration of the CT data with the ToF data, the mean TRE ranged from 3.7 ± 2.1 mm for point-based and 5.7 ± 1.9 mm for surface-based registration for the best camera to 6.2 ± 3.4 and 11.1 ± 2.8 mm, respectively. Compared to data provided by the manufacturer, the mean TRE decreased by 8%–60% for point-based and by 18%–74% for surface-based registration.Conclusions: Using the proposed calibration approach improved the measurement accuracy of all investigated ToF cameras. Although evaluated in the context of intraoperative image acquisition, the proposed calibration procedure can easily be applied to other medical applications using ToF cameras, such as patient positioning or respiratory motion tracking in radiotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

In-line phase shift tomosynthesis

Purpose: The purpose of this work is to (1) demonstrate laboratory measurements of phase shift images derived from in-line phase-contrast radiographs using the attenuation-partition based algorithm (APBA) of Yan et al.[Opt. Express 18(15), 16074–16089 (2010)], (2) verify that the APBA reconstructed images obey the linearity principle, and (3) reconstruct tomosynthesis phase shift images from a collection of angularly sampled planar phase shift images.Methods: An unmodified, commercially available cabinet x-ray system (Faxitron LX-60) was used in this experiment. This system contains a tungsten anode x-ray tube with a nominal focal spot size of 10 μm. The digital detector uses CsI/CMOS with a pixel size of 50 × 50 μm. The phantoms used consisted of one acrylic plate, two polystyrene plates, and a habanero pepper. Tomosynthesis images were reconstructed from 51 images acquired over a ±25° arc. All phase shift images were reconstructed using the APBA.Results: Image contrast derived from the planar phase shift image of an acrylic plate of uniform thickness exceeded the contrast of the traditional attenuation image by an approximate factor of two. Comparison of the planar phase shift images from a single, uniform thickness polystyrene plate with two polystyrene plates demonstrated an approximate linearity of the estimated phase shift with plate thickness (−1600 rad vs −2970 rad). Tomographic phase shift images of the habanero pepper exhibited acceptable spatial resolution and contrast comparable to the corresponding attenuation image.Conclusions: This work demonstrated the feasibility of laboratory-based phase shift tomosynthesis and suggests that phase shift imaging could potentially provide a new imaging biomarker. Further investigation will be needed to determine if phase shift contrast will be able to provide new tissue contrast information or improved clinical performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Influence of cold walls on PET image quantification and volume segmentation: A phantom study

Purpose: Commercially available fillable plastic inserts used in positron emission tomography phantoms usually have thick plastic walls, separating their content from the background activity. These “cold” walls can modify the intensity values of neighboring active regions due to the partial volume effect, resulting in errors in the estimation of standardized uptake values. Numerous papers suggest that this is an issue for phantom work simulating tumor tissue, quality control, and calibration work. This study aims to investigate the influence of the cold plastic wall thickness on the quantification of 18F-fluorodeoxyglucose on the image activity recovery and on the performance of advanced automatic segmentation algorithms for the delineation of active regions delimited by plastic walls.Methods: A commercial set of six spheres of different diameters was replicated using a manufacturing technique which achieves a reduction in plastic walls thickness of up to 90%, while keeping the same internal volume. Both sets of thin- and thick-wall inserts were imaged simultaneously in a custom phantom for six different tumor-to-background ratios. Intensity values were compared in terms of mean and maximum standardized uptake values (SUVs) in the spheres and mean SUV of the hottest 1 ml region (SUV{sub max}, SUV{sub mean}, and SUV{sub peak}). The recovery coefficient (RC) was also derived for each sphere. The results were compared against the values predicted by a theoretical model of the PET-intensity profiles for the same tumor-to-background ratios (TBRs), sphere sizes, and wall thicknesses. In addition, ten automatic segmentation methods, written in house, were applied to both thin- and thick-wall inserts. The contours obtained were compared to computed tomography derived gold standard (“ground truth”), using five different accuracy metrics.Results: The authors' results showed that thin-wall inserts achieved significantly higher SUV{sub mean}, SUV{sub max}, and RC values (up to 25%, 16%, and 25% higher, respectively) compared to thick-wall inserts, which was in agreement with the theory. This effect decreased with increasing sphere size and TBR, and resulted in substantial (>5%) differences between thin- and thick-wall inserts for spheres up to 30 mm diameter and TBR up to 4. Thinner plastic walls were also shown to significantly improve the delineation accuracy for the majority of the segmentation methods tested, by increasing the proportion of lesion voxels detected, although the errors in image quantification remained non-negligible.Conclusions: This study quantified the significant effect of a 90% reduction in the thickness of insert walls on SUV quantification and PET-based boundary detection. Mean SUVs inside the inserts and recovery coefficients were particularly affected by the presence of thick cold walls, as predicted by a theoretical approach. The accuracy of some delineation algorithms was also significantly improved by the introduction of thin wall inserts instead of thick wall inserts. This study demonstrates the risk of errors deriving from the use of cold wall inserts to assess and compare the performance of PET segmentation methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Heel effect adaptive flat field correction of digital x-ray detectors

Purpose: Anode heel effect renders large-scale background nonuniformities in digital radiographs. Conventional offset/gain calibration is performed at mono source-to-image distance (SID), and disregards the SID-dependent characteristic of heel effect. It results in a residual nonuniform background in the corrected radiographs when the SID settings for calibration and correction differ. In this work, the authors develop a robust and efficient computational method for digital x-ray detector gain correction adapted to SID-variant heel effect, without resorting to physical filters, phantoms, complicated heel effect models, or multiple-SID calibration and interpolation.Methods: The authors present the Duo-SID projection correction method. In our approach, conventional offset/gain calibrations are performed only twice, at the minimum and maximum SIDs of the system in typical clinical use. A fast iterative separation algorithm is devised to extract the detector gain and basis heel patterns from the min/max SID calibrations. The resultant detector gain is independent of SID, while the basis heel patterns are parameterized by the min- and max-SID. The heel pattern at any SID is obtained from the min-SID basis heel pattern via projection imaging principles. The system gain desired at a specific acquisition SID is then constructed using the projected heel pattern and detector gain map.Results: The method was evaluated for flat field and anatomical phantom image corrections. It demonstrated promising improvements over interpolation and conventional gain calibration/correction methods, lowering their correction errors by approximately 70% and 80%, respectively. The separation algorithm was able to extract the detector gain and heel patterns with less than 2% error, and the Duo-SID corrected images showed perceptually appealing uniform background across the detector.Conclusions: The Duo-SID correction method has substantially improved on conventional offset/gain corrections for digital x-ray imaging in an SID-variant environment. The technique is relatively simple, and can be easily incorporated into multiple-point gain calibration/correction techniques. It offers a potentially valuable tool for preprocessing digital x-ray images to boost image quality of mammography, chest and cardiac radiography, as well as automated computer aided diagnostic radiology.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Investigation of energy weighting using an energy discriminating photon counting detector for breast CT

Purpose: Breast CT is an emerging imaging technique that can portray the breast in 3D and improve visualization of important diagnostic features. Early clinical studies have suggested that breast CT has sufficient spatial and contrast resolution for accurate detection of masses and microcalcifications in the breast, reducing structural overlap that is often a limiting factor in reading mammographic images. For a number of reasons, image quality in breast CT may be improved by use of an energy resolving photon counting detector. In this study, the authors investigate the improvements in image quality obtained when using energy weighting with an energy resolving photon counting detector as compared to that with a conventional energy integrating detector.Methods: Using computer simulation, realistic CT images of multiple breast phantoms were generated. The simulation modeled a prototype breast CT system using an amorphous silicon (a-Si), CsI based energy integrating detector with different x-ray spectra, and a hypothetical, ideal CZT based photon counting detector with capability of energy discrimination. Three biological signals of interest were modeled as spherical lesions and inserted into breast phantoms; hydroxyapatite (HA) to represent microcalcification, infiltrating ductal carcinoma (IDC), and iodine enhanced infiltrating ductal carcinoma (IIDC). Signal-to-noise ratio (SNR) of these three lesions was measured from the CT reconstructions. In addition, a psychophysical study was conducted to evaluate observer performance in detecting microcalcifications embedded into a realistic anthropomorphic breast phantom.Results: In the energy range tested, improvements in SNR with a photon counting detector using energy weighting was higher (than the energy integrating detector method) by 30%–63% and 4%–34%, for HA and IDC lesions and 12%–30% (with Al filtration) and 32%–38% (with Ce filtration) for the IIDC lesion, respectively. The average area under the receiver operating characteristic curve (AUC) for detection of microcalcifications was higher by greater than 19% (for the different energy weighting methods tested) as compared to the AUC obtained with an energy integrating detector.Conclusions: This study showed that breast CT with a CZT photon counting detector using energy weighting can provide improvements in pixel SNR, and detectability of microcalcifications as compared to that with a conventional energy integrating detector. Since a number of degrading physical factors were not modeled into the photon counting detector, this improvement should be considered as an upper bound on achievable performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A measurement technique to determine the calibration accuracy of an electromagnetic tracking system to radiation isocenter

Purpose: To present and characterize a measurement technique to quantify the calibration accuracy of an electromagnetic tracking system to radiation isocenter.Methods: This technique was developed as a quality assurance method for electromagnetic tracking systems used in a multi-institutional clinical hypofractionated prostate study. In this technique, the electromagnetic tracking system is calibrated to isocenter with the manufacturers recommended technique, using laser-based alignment. A test patient is created with a transponder at isocenter whose position is measured electromagnetically. Four portal images of the transponder are taken with collimator rotations of 45° 135°, 225°, and 315°, at each of four gantry angles (0°, 90°, 180°, 270°) using a 3 × 6 cm{sup 2} radiation field. In each image, the center of the copper-wrapped iron core of the transponder is determined. All measurements are made relative to this transponder position to remove gantry and imager sag effects. For each of the 16 images, the 50% collimation edges are identified and used to find a ray representing the rotational axis of each collimation edge. The 16 collimator rotation rays from four gantry angles pass through and bound the radiation isocenter volume. The center of the bounded region, relative to the transponder, is calculated and then transformed to tracking system coordinates using the transponder position, allowing the tracking system's calibration offset from radiation isocenter to be found. All image analysis and calculations are automated with inhouse software for user-independent accuracy. Three different tracking systems at two different sites were evaluated for this study.Results: The magnitude of the calibration offset was always less than the manufacturer's stated accuracy of 0.2 cm using their standard clinical calibration procedure, and ranged from 0.014 to 0.175 cm. On three systems in clinical use, the magnitude of the offset was found to be 0.053 ± 0.036, 0.121 ± 0.023, and 0.093 ± 0.013 cm.Conclusions: The method presented here provides an independent technique to verify the calibration of an electromagnetic tracking system to radiation isocenter. The calibration accuracy of the system was better than the 0.2 cm accuracy stated by the manufacturer. However, it should not be assumed to be zero, especially for stereotactic radiation therapy treatments where planning target volume margins are very small.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Spatial and dose–response analysis of fibrotic lung changes after stereotactic body radiation therapy

Purpose: Stereotactic body radiation therapy (SBRT) is becoming the standard of care for early stage nonoperable lung cancers. Accurate dose–response modeling is challenging for SBRT because of the decreased number of clinical toxicity events. As a surrogate for a clinical toxicity endpoint, studies have proposed to use radiographic changes in follow up computed tomography (CT) scans to evaluate lung SBRT normal tissue effects. The purpose of the current study was to use local fibrotic lung regions to spatially and dosimetrically evaluate lung changes in patients that underwent SBRT.Methods: Forty seven SBRT patients treated at our institution from 2003 to 2009 were used for the current study. Our patient cohort had a total of 148 follow up CT scans ranging from 3 to 48 months post-therapy. Post-treatment scans were binned into intervals of 3, 6, 12, 18, 24, 30, and 36 months after the completion of treatment. Deformable image registration was used to align the follow up CT scans with the pretreatment CT and dose distribution. Areas of visible fibrotic changes were contoured. The centroid of each gross tumor volume (GTV) and contoured fibrosis volume was calculated and the fibrosis volume location and movement (magnitude and direction) relative to the GTV and 30 Gy isodose centroid were analyzed. To perform a dose–response analysis, each voxel in the fibrosis volume was sorted into 10 Gy dose bins and the average CT number value for each dose bin was calculated. Dose–response curves were generated by plotting the CT number as a function of dose bin and time posttherapy.Results: Both fibrosis and GTV centroids were concentrated in the upper third of the lung. The average radial movement of fibrosis centroids relative to the GTV centroids was 2.6 cm with movement greater than 5 cm occurring in 11% of patients. Evaluating dose–response curves revealed an overall trend of increasing CT number as a function of dose. The authors observed a CT number plateau at doses ranging from 30 to 50 Gy for the 3, 6, and 12 months posttherapy time points. There was no evident plateau for the dose–response curves generated using data from the 18, 24, 30, and 36 months posttherapy time points.Conclusions: Regions of local fibrotic lung changes in patients that underwent SBRT were evaluated spatially and dosimetrically. The authors found that the average fibrosis movement was 2.6 cm with movement greater than 5 cm possible. Evaluating dose–response curves revealed an overall trend of increasing CT number as a function of dose. Furthermore, our dose–response data also suggest that one of the possible explanations of the CT number plateau effect may be the time posttherapy of the acquired data. Understanding normal tissue dose–response is important for reducing toxicity after SBRT, especially in cases where larger tumors are treated. The methods presented in the current work build on prior quantitative studies and further enhance the understanding of normal lung dose–response after SBRT.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

CT head-scan dosimetry in an anthropomorphic phantom and associated measurement of ACR accreditation-phantom imaging metrics under clinically representative scan conditions

Purpose: To measure radiation absorbed dose and its distribution in an anthropomorphic head phantom under clinically representative scan conditions in three widely used computed tomography (CT) scanners, and to relate those dose values to metrics such as high-contrast resolution, noise, and contrast-to-noise ratio (CNR) in the American College of Radiology CT accreditation phantom.Methods: By inserting optically stimulated luminescence dosimeters (OSLDs) in the head of an anthropomorphic phantom specially developed for CT dosimetry (University of Florida, Gainesville), we measured dose with three commonly used scanners (GE Discovery CT750 HD, Siemens Definition, Philips Brilliance 64) at two different clinical sites (Walter Reed National Military Medical Center, National Institutes of Health). The scanners were set to operate with the same data-acquisition and image-reconstruction protocols as used clinically for typical head scans, respective of the practices of each facility for each scanner. We also analyzed images of the ACR CT accreditation phantom with the corresponding protocols. While the Siemens Definition and the Philips Brilliance protocols utilized only conventional, filtered back-projection (FBP) image-reconstruction methods, the GE Discovery also employed its particular version of an adaptive statistical iterative reconstruction (ASIR) algorithm that can be blended in desired proportions with the FBP algorithm. We did an objective image-metrics analysis evaluating the modulation transfer function (MTF), noise power spectrum (NPS), and CNR for images reconstructed with FBP. For images reconstructed with ASIR, we only analyzed the CNR, since MTF and NPS results are expected to depend on the object for iterative reconstruction algorithms.Results: The OSLD measurements showed that the Siemens Definition and the Philips Brilliance scanners (located at two different clinical facilities) yield average absorbed doses in tissue of 42.6 and 43.1 mGy, respectively. The GE Discovery delivers about the same amount of dose (43.7 mGy) when run under similar operating and image-reconstruction conditions, i.e., without tube current modulation and ASIR. The image-metrics analysis likewise showed that the MTF, NPS, and CNR associated with the reconstructed images are mutually comparable when the three scanners are run with similar settings, and differences can be attributed to different edge-enhancement properties of the applied reconstruction filters. Moreover, when the GE scanner was operated with the facility's scanner settings for routine head exams, which apply 50% ASIR and use only approximately half of the 100%-FBP dose, the CNR of the images showed no significant change. Even though the CNR alone is not sufficient to characterize the image quality and justify any dose reduction claims, it can be useful as a constancy test metric.Conclusions: This work presents a straightforward method to connect direct measurements of CT dose with objective image metrics such as high-contrast resolution, noise, and CNR. It demonstrates that OSLD measurements in an anthropomorphic head phantom allow a realistic and locally precise estimation of magnitude and spatial distribution of dose in tissue delivered during a typical CT head scan. Additional objective analysis of the images of the ACR accreditation phantom can be used to relate the measured doses to high contrast resolution, noise, and CNR.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Gantry-angle resolved VMAT pretreatment verification using EPID image prediction

Purpose: Pretreatment verification of volumetric modulated arc therapy (VMAT) dose delivery with electronic portal imaging device (EPID) uses images integrated over the entire delivery or over large subarcs. This work aims to develop a new method for gantry-angle-resolved verification of VMAT dose delivery using EPID.Methods: An EPID dose prediction model was used to calculate EPID images as a function of gantry angle for eight prostate patient deliveries. EPID image frames at 7.5 frames per second were acquired during delivery via a frame-grabber system. The gantry angle for each image was encoded in kV frames which were synchronized to the MV frames. Gamma analysis results as a function of gantry angle were assessed by integrating the frames over 2° subarcs with an angle-to-agreement tolerance of 0.5° about the measured image angle.Results: The model agreed with EPID images integrated over the entire delivery with average Gamma pass-rates at 2%, 2 mm of 99.7% (10% threshold). The accuracy of the kV derived gantry angle for each image was found to be 0.1° (1 SD) using a phantom test. For the gantry-resolved analysis all Gamma pass-rates were greater than 90% at 3%, 3 mm criteria (with only two exceptions), and more than 90% had a 95% pass-rate, with an average of 97.3%. The measured gantry angle lagged behind the predicted angle by a mean of 0.3°± 0.3°, with a maximum lag of 1.3°.Conclusions: The method provides a comprehensive and highly efficient pretreatment verification of VMAT delivery using EPID. Dose delivery accuracy is assessed as a function of gantry angle to ensure accurate treatment.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Characterization of a novel EPID designed for simultaneous imaging and dose verification in radiotherapy

Purpose: Standard amorphous silicon electronic portal imaging devices (a-Si EPIDs) are x-ray imagers used frequently in radiotherapy that indirectly detect incident x-rays using a metal plate and phosphor screen. These detectors may also be used as two-dimensional dosimeters; however, they have a well-characterized nonwater-equivalent dosimetric response. Plastic scintillating (PS) fibers, on the other hand, have been shown to respond in a water-equivalent manner to x-rays in the energy range typically encountered during radiotherapy. In this study, the authors report on the first experimental measurements taken with a novel prototype PS a-Si EPID developed for the purpose of performing simultaneous imaging and dosimetry in radiotherapy. This prototype employs an array of PS fibers in place of the standard metal plate and phosphor screen. The imaging performance and dosimetric response of the prototype EPID were evaluated experimentally and compared to that of the standard EPID.Methods: Clinical 6 MV photon beams were used to first measure the detector sensitivity, linearity of dose response, and pixel noise characteristics of the prototype and standard EPIDs. Second, the dosimetric response of each EPID was evaluated relative to a reference water-equivalent dosimeter by measuring the off-axis and field size response in a nontransit configuration, along with the off-axis, field size, and transmission response in a transit configuration using solid water blocks. Finally, the imaging performance of the prototype and standard EPIDs was evaluated quantitatively by using an image quality phantom to measure the contrast to noise ratio (CNR) and spatial resolution of images acquired with each detector, and qualitatively by using an anthropomorphic phantom to acquire images representative of human anatomy.Results: The prototype EPID's sensitivity was 0.37 times that of the standard EPID. Both EPIDs exhibited responses that were linear with delivered dose over a range of 1–100 monitor units. Over this range, the prototype and standard EPID central axis responses agreed to within 1.6%. Images taken with the prototype EPID were noisier than those taken with the standard EPID, with fractional uncertainties of 0.2% and 0.05% within the central 1 cm{sup 2}, respectively. For all dosimetry measurements, the prototype EPID exhibited a near water-equivalent response whereas the standard EPID did not. The CNR and spatial resolution of images taken with the standard EPID were greater than those taken with the prototype EPID.Conclusions: A prototype EPID employing an array of PS fibers has been developed and the first experimental measurements are reported. The prototype EPID demonstrated a much morewater-equivalent dose response than the standard EPID. While the imaging performance of the standard EPID was superior to that of the prototype, the prototype EPID has many design characteristics that may be optimized to improve imaging performance. This investigation demonstrates the feasibility of a new detector design for simultaneous imaging and dosimetry treatment verification in radiotherapy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Joint surface reconstruction and 4D deformation estimation from sparse data and prior knowledge for marker-less Respiratory motion tracking

Purpose: The intraprocedural tracking of respiratory motion has the potential to substantially improve image-guided diagnosis and interventions. The authors have developed a sparse-to-dense registration approach that is capable of recovering the patient's external 3D body surface and estimating a 4D (3D + time) surface motion field from sparse sampling data and patient-specific prior shape knowledge.Methods: The system utilizes an emerging marker-less and laser-based active triangulation (AT) sensor that delivers sparse but highly accurate 3D measurements in real-time. These sparse position measurements are registered with a dense reference surface extracted from planning data. Thereby a dense displacement field is recovered, which describes the spatio-temporal 4D deformation of the complete patient body surface, depending on the type and state of respiration. It yields both a reconstruction of the instantaneous patient shape and a high-dimensional respiratory surrogate for respiratory motion tracking. The method is validated on a 4D CT respiration phantom and evaluated on both real data from an AT prototype and synthetic data sampled from dense surface scans acquired with a structured-light scanner.Results: In the experiments, the authors estimated surface motion fields with the proposed algorithm on 256 datasets from 16 subjects and in different respiration states, achieving a mean surface reconstruction accuracy of ±0.23 mm with respect to ground truth data—down from a mean initial surface mismatch of 5.66 mm. The 95th percentile of the local residual mesh-to-mesh distance after registration did not exceed 1.17 mm for any subject. On average, the total runtime of our proof of concept CPU implementation is 2.3 s per frame, outperforming related work substantially.Conclusions: In external beam radiation therapy, the approach holds potential for patient monitoring during treatment using the reconstructed surface, and for motion-compensated dose delivery using the estimated 4D surface motion field in combination with external-internal correlation models.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Digital tomosynthesis (DTS) for verification of target position in early stage lung cancer patients

Purpose: The ability to verify intrafraction tumor position is clinically useful for hypofractionated treatments. Short arc kV digital tomosynthesis (DTS) could facilitate more frequent target verification. The authors used DTS combined with triangulation to determine the mean temporal position of small-volume lung tumor targets treated with stereotactic radiotherapy. DTS registration results were benchmarked against online clinical localization using registration between free-breathing cone-beam computed tomography (CBCT) and the average intensity projection (AvIP) of the planning 4DCT.Methods: In this retrospective study, 76 sets of kV-projection images from online CBCT scans of 13 patients were used to generate DTS image slices (CB-DTS) with nonclinical research software (DTS Toolkit, Varian Medical Systems). Three-dimensional tumor motion was 1.3–4 mm in six patients and 6.1–25.4 mm in seven patients on 4DCT (significant difference in the mean of the groups, P < 0.01). The 4DCT AvIP was used to digitally reconstruct the Reference-DTS. DTS registration and DTS registration combined with triangulation were investigated. Progressive shortening of total DTS arc lengths from 95° to 35° around 0° gantry position was evaluated for different scenarios: DTS registration using the entire arc; DTS registration plus triangulation using two nonoverlapping arcs; and for 55° and 45° total gantry rotation, DTS registration plus triangulation using two overlapping arcs. Finally, DTS registration plus triangulation performed at eight gantry angles, each separated by 45° was evaluated using full fan kV projection data for one patient with an immobile tumor and five patients with mobile tumors.Results: For DTS registration alone, shortening arc length did not influence accuracy in X- and Y-directions, but in Z-direction, mean deviations from online CBCT localization systematically increased for shorter arc length (P < 0.05). For example, using a 95° arc mean DTS-CBCT difference was 0.8 mm (1 SD = 0.6 mm) and for a 35° arc the mean was 2.4 mm (1 SD = 1.7 mm). DTS plus triangulation using nonoverlapping-arcs increased accuracy in Z-direction for tested arc lengths ≤55° (P < 0.01). Overlapping arcs increased accuracy in Y-direction for tumors with motion >4 mm (P < 0.02) but increased Z-direction accuracy was only observed with 55° total gantry rotation. The 95th percentile deviations with this overlapping technique in X-, Y-, and Z-directions were 1.3, 2.0, and 2.5 mm, respectively. For the five patients with mobile tumors where DTS + triangulation was performed with 45° intervals, the pooled deviation from online CBCT correction showed, for X-, Y-, and Z-directions, mean of 1.1 mm, standard deviations (SD) of 0.9, 1.0, and 0.9 mm, respectively. The mean + 2 SD was <3 mm for each direction.Conclusions: Short-arc DTS verification of time averaged lung tumor position is feasible using free-breathing kV projection data and the AvIP of the 4DCT as a reference. Observed differences between DTS and online CBCT registration with AvIP were ≤3 mm (mean + 2 SD), however, the increased temporal resolution of DTS + triangulation also identified short period deviations from the average target position on the CBCT. Short-arc DTS appears promising for intrafraction tumor position monitoring during stereotactic lung radiotherapy delivered with a rotational technique.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A system for EPID-based real-time treatment delivery verification during dynamic IMRT treatment

Purpose: To design and develop a real-time electronic portal imaging device (EPID)-based delivery verification system for dynamic intensity modulated radiation therapy (IMRT) which enables detection of gross treatment delivery errors before delivery of substantial radiation to the patient.Methods: The system utilizes a comprehensive physics-based model to generate a series of predicted transit EPID image frames as a reference dataset and compares these to measured EPID frames acquired during treatment. The two datasets are using MLC aperture comparison and cumulative signal checking techniques. The system operation in real-time was simulated offline using previously acquired images for 19 IMRT patient deliveries with both frame-by-frame comparison and cumulative frame comparison. Simulated error case studies were used to demonstrate the system sensitivity and performance.Results: The accuracy of the synchronization method was shown to agree within two control points which corresponds to approximately ∼1% of the total MU to be delivered for dynamic IMRT. The system achieved mean real-time gamma results for frame-by-frame analysis of 86.6% and 89.0% for 3%, 3 mm and 4%, 4 mm criteria, respectively, and 97.9% and 98.6% for cumulative gamma analysis. The system can detect a 10% MU error using 3%, 3 mm criteria within approximately 10 s. The EPID-based real-time delivery verification system successfully detected simulated gross errors introduced into patient plan deliveries in near real-time (within 0.1 s).Conclusions: A real-time radiation delivery verification system for dynamic IMRT has been demonstrated that is designed to prevent major mistreatments in modern radiation therapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2013

Determination of lung segments in computed tomography images using the Euclidean distance to the pulmonary artery

Purpose: Computed tomography (CT) imaging is the modality of choice for lung cancer diagnostics. With the increasing number of lung interventions on sublobar level in recent years, determining and visualizing pulmonary segments in CT images and, in oncological cases, reliable segment-related information about the location of tumors has become increasingly desirable. Computer-assisted identification of lung segments in CT images is subject of this work.Methods: The authors present a new interactive approach for the segmentation of lung segments that uses the Euclidean distance of each point in the lung to the segmental branches of the pulmonary artery. The aim is to analyze the potential of the method. Detailed manual pulmonary artery segmentations are used to achieve the best possible segment approximation results. A detailed description of the method and its evaluation on 11 CT scans from clinical routine are given.Results: An accuracy of 2–3 mm is measured for the segment boundaries computed by the pulmonary artery-based method. On average, maximum deviations of 8 mm are observed. 135 intersegmental pulmonary veins detected in the 11 test CT scans serve as reference data. Furthermore, a comparison of the presented pulmonary artery-based approach to a similar approach that uses the Euclidean distance to the segmental branches of the bronchial tree is presented. It shows a significantly higher accuracy for the pulmonary artery-based approach in lung regions at least 30 mm distal to the lung hilum.Conclusions: A pulmonary artery-based determination of lung segments in CT images is promising. In the tests, the pulmonary artery-based determination has been shown to be superior to the bronchial tree-based determination. The suitability of the segment approximation method for application in the planning of segment resections in clinical practice has already been verified in experimental cases. However, automation of the method accompanied by an evaluation on a larger number of test cases is required before application in the daily clinical routine.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Real-time soft tissue motion estimation for lung tumors during radiotherapy delivery

Purpose: To provide real-time lung tumor motion estimation during radiotherapy treatment delivery without the need for implanted fiducial markers or additional imaging dose to the patient.Methods: 2D radiographs from the therapy beam's-eye-view (BEV) perspective are captured at a frame rate of 12.8 Hz with a frame grabber allowing direct RAM access to the image buffer. An in-house developed real-time soft tissue localization algorithm is utilized to calculate soft tissue displacement from these images in real-time. The system is tested with a Varian TX linear accelerator and an AS-1000 amorphous silicon electronic portal imaging device operating at a resolution of 512 × 384 pixels. The accuracy of the motion estimation is verified with a dynamic motion phantom. Clinical accuracy was tested on lung SBRT images acquired at 2 fps.Results: Real-time lung tumor motion estimation from BEV images without fiducial markers is successfully demonstrated. For the phantom study, a mean tracking error <1.0 mm [root mean square (rms) error of 0.3 mm] was observed. The tracking rms accuracy on BEV images from a lung SBRT patient (≈20 mm tumor motion range) is 1.0 mm.Conclusions: The authors demonstrate for the first time real-time markerless lung tumor motion estimation from BEV images alone. The described system can operate at a frame rate of 12.8 Hz and does not require prior knowledge to establish traceable landmarks for tracking on the fly. The authors show that the geometric accuracy is similar to (or better than) previously published markerless algorithms not operating in real-time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

4D tumor centroid tracking using orthogonal 2D dynamic MRI: Implications for radiotherapy planning

Purpose: Current pretreatment, 4D imaging techniques are suboptimal in that they sample breathing motion over a very limited “snapshot” in time. Heretofore, long-duration, 4D motion characterization for radiotherapy planning, margin optimization, and validation have been impractical for safety reasons, requiring invasive markers imaged under x-ray fluoroscopy. To characterize 3D tumor motion and associated variability over durations more consistent with treatments, the authors have developed a practical dynamic MRI (dMRI) technique employing two orthogonal planes acquired in a continuous, interleaved fashion.Methods: 2D balanced steady-state free precession MRI was acquired continuously over 9–14 min at approximately 4 Hz in three healthy volunteers using a commercial 1.5 T system; alternating orthogonal imaging planes (sagittal, coronal, sagittal, etc.) were employed. The 2D in-plane pixel resolution was 2 × 2 mm{sup 2} with a 5 mm slice profile. Simultaneous with image acquisition, the authors monitored a 1D surrogate respiratory signal using a device available with the MRI system. 2D template matching-based anatomic feature registration, or tracking, was performed independently in each orientation. 4D feature tracking at the raw frame rate was derived using spline interpolation.Results: Tracking vascular features in the lung for two volunteers and pancreatic features in one volunteer, the authors have successfully demonstrated this method. Registration error, defined here as the difference between the sagittal and coronal tracking result in the SI direction, ranged from 0.7 to 1.6 mm (1σ) which was less than the acquired image resolution. Although the healthy volunteers were instructed to relax and breathe normally, significantly variable respiration was observed. To demonstrate potential applications of this technique, the authors subsequently explored the intrafraction stability of hypothetical tumoral internal target volumes and 3D spatial probability distribution functions. The surrogate respiratory information allowed the authors to show how this technique can be used to study correlations between internal and external (surrogate) information over these prolonged durations. However, compared against the gold standard of the time stamps in the dMRI frames, the temporal synchronization of the surrogate 1D respiratory information was shown to be likely unreliable.Conclusions: The authors have established viability of a novel and practical pretreatment, 4D tumor centroid tracking method employing a commercially available dynamic MRI sequence. Further developments from the vendor are likely needed to provide a reliably synchronized surrogate 1D respiratory signal, which will likely broaden the utility of this method in the pretreatment radiotherapy planning context.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Low dose tomographic fluoroscopy: 4D intervention guidance with running prior

Purpose: Today's standard imaging technique in interventional radiology is the single- or biplane x-ray fluoroscopy which delivers 2D projection images as a function of time (2D+T). This state-of-the-art technology, however, suffers from its projective nature and is limited by the superposition of the patient's anatomy. Temporally resolved tomographic volumes (3D+T) would significantly improve the visualization of complex structures. A continuous tomographic data acquisition, if carried out with today's technology, would yield an excessive patient dose. Recently the authors proposed a method that enables tomographic fluoroscopy at the same dose level as projective fluoroscopy which means that if scanning time of an intervention guided by projective fluoroscopy is the same as that of an intervention guided by tomographic fluoroscopy, almost the same dose is administered to the patient. The purpose of this work is to extend authors' previous work and allow for patient motion during the intervention.Methods: The authors propose the running prior technique for adaptation of a prior image. This adaptation is realized by a combination of registration and projection replacement. In a first step the prior is deformed to the current position via affine and deformable registration. Then the information from outdated projections is replaced by newly acquired projections using forward and backprojection steps. The thus adapted volume is the running prior. The proposed method is validated by simulated as well as measured data. To investigate motion during intervention a moving head phantom was simulated. Real in vivo data of a pig are acquired by a prototype CT system consisting of a flat detector and a continuously rotating clinical gantry.Results: With the running prior technique it is possible to correct for motion without additional dose. For an application in intervention guidance both steps of the running prior technique, registration and replacement, are necessary. Reconstructed volumes based on the running prior show high image quality without introducing new artifacts and the interventional materials are displayed at the correct position.Conclusions: The running prior improves the robustness of low dose 3D+T intervention guidance toward intended or unintended patient motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗