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 847 records · Page 47

Image quality evaluation of breast tomosynthesis with synchrotron radiation

Purpose: This study investigates the image quality of tomosynthesis slices obtained from several acquisition sets with synchrotron radiation using a breast phantom incorporating details that mimic various breast lesions, in a heterogeneous background. Methods: A complex Breast phantom (MAMMAX) with a heterogeneous background and thickness that corresponds to 4.5 cm compressed breast with an average composition of 50% adipose and 50% glandular tissue was assembled using two commercial phantoms. Projection images using acquisition arcs of 24 Degree-Sign , 32 Degree-Sign , 40 Degree-Sign , 48 Degree-Sign , and 56 Degree-Sign at incident energy of 17 keV were obtained from the phantom with the synchrotron radiation for medical physics beamline at ELETTRA Synchrotron Light Laboratory. The total mean glandular dose was set equal to 2.5 mGy. Tomograms were reconstructed with simple multiple projection algorithm (MPA) and filtered MPA. In the latter case, a median filter, a sinc filter, and a combination of those two filters were applied on the experimental data prior to MPA reconstruction. Visual inspection, contrast to noise ratio, contrast, and artifact spread function were the figures of merit used in the evaluation of the visualisation and detection of low- and high-contrast breast features, as a function of the reconstruction algorithm and acquisition arc. To study the benefits of using monochromatic beams, single projection images at incident energies ranging from 14 to 27 keV were acquired with the same phantom and weighted to synthesize polychromatic images at a typical incident x-ray spectrum with W target. Results: Filters were optimised to reconstruct features with different attenuation characteristics and dimensions. In the case of 6 mm low-contrast details, improved visual appearance as well as higher contrast to noise ratio and contrast values were observed for the two filtered MPA algorithms that exploit the sinc filter. These features are better visualized at extended arc length, as the acquisition arc of 56 Degree-Sign with 15 projection images demonstrates the highest image reconstruction quality. For microcalcifications, filtered MPA implemented with a combination of median and sinc filters indicates better feature appearance due to efficient suppression of background tissue. The image quality of these features is less sensitive to the acquisition arc. Calcifications with size ranging from 170 to 500 {mu}m, like the ones presently studied, are well identified and visualized for all arcs used. The comparison of single projection images obtained under different beam conditions showed that the use of monochromatic beam can produce an image with higher contrast and contrast to noise ratio compared to an image corresponding to a polychromatic beam even when the latter is acquired with double incident exposure. Conclusions: Filter optimization in respect to the type of feature characteristics is important before the reconstruction. The MPA combined with median and sinc filters results in improved reconstruction of microcalcifications and low-contrast features. The latter are better visualized at extended arc length, while microcalcifications are less sensitive to this acquisition parameter. Use of monochromatic beams may result in tomographic images with higher contrast acquired at lower incident exposures.

60 APPLIED LIFE SCIENCES↗

Target-specific optimization of four-dimensional cone beam computed tomography

Purpose: Under-sampling artifacts are a major problem in four-dimensional cone-beam computed tomography (4D-CBCT) and may compromise evaluation of target motion. Several studies have addressed scan parameter selection for minimizing under-sampling artifacts; however, the role of the target characteristics in scan parameter selection has not been investigated. In this work, the authors evaluated 4D-CBCT performance by assessing the accuracy of target motion measurements for various target sizes and motions. The results may serve as patient-specific guidelines for scan parameters selection in 4D-CBCT. Methods: The authors acquired 4D-CBCT scans of a moving phantom consisting of six water-filled sphere targets of sizes 10-37 mm, using various scan times ranging from 30 s to 3 min., setting the motion to 3-s and 6-s periods. The authors used automatic image registration to extract the target motion trajectories and evaluated these measurements for various target sizes and motions over various combinations of scan parameters including scan time, detector configuration, number of respiration phases, and reconstruction filters. Results: The most important object parameter to 4D-CBCT performance was the period of motion. Measurements for the 6-s motion were always systematically less accurate than measurements for the 3-s motion for 34 of 36 objects of various sizes and periods of motion. The 6-s motion required a greater scan time than the 3-s motion for equivalent measurement accuracy. The second most influential parameter was the target size. For the 3-s period of motion, objects larger than 13 mm were tracked with sub-millimeter accuracy with a 1-min scan time. For the 6-s period of motion, objects larger than 22-mm were tracked with sub-millimeter accuracy with a 1.5-min scan time. For all sizes and motions, temporal blurring was observed when the number of phases was fewer than 8. Offset detector configuration provided the same performance as centered detector except for small targets (20 mm.) at short scan times ({<=}1 min). Finally, reconstruction filtering and number of respiratory phases did not affect performance. Conclusions: Scan time should be set according to target size and motion. The authors have provided figures that provide the minimum scan time needed to achieve the particular motion measurement accuracy for a particular size and motion period.

60 APPLIED LIFE SCIENCES↗

A study of respiration-correlated cone-beam CT scans to correct target positioning errors in radiotherapy of thoracic cancer

Purpose: There is increasingly widespread usage of cone-beam CT (CBCT) for guiding radiation treatment in advanced-stage lung tumors, but difficulties associated with daily CBCT in conventionally fractionated treatments include imaging dose to the patient, increased workload and longer treatment times. Respiration-correlated cone-beam CT (RC-CBCT) can improve localization accuracy in mobile lung tumors, but further increases the time and workload for conventionally fractionated treatments. This study investigates whether RC-CBCT-guided correction of systematic tumor deviations in standard fractionated lung tumor radiation treatments is more effective than 2D image-based correction of skeletal deviations alone. A second study goal compares respiration-correlated vs respiration-averaged images for determining tumor deviations. Methods: Eleven stage II-IV nonsmall cell lung cancer patients are enrolled in an IRB-approved prospective off-line protocol using RC-CBCT guidance to correct for systematic errors in GTV position. Patients receive a respiration-correlated planning CT (RCCT) at simulation, daily kilovoltage RC-CBCT scans during the first week of treatment and weekly scans thereafter. Four types of correction methods are compared: (1) systematic error in gross tumor volume (GTV) position, (2) systematic error in skeletal anatomy, (3) daily skeletal corrections, and (4) weekly skeletal corrections. The comparison is in terms of weighted average of the residual GTV deviations measured from the RC-CBCT scans and representing the estimated residual deviation over the treatment course. In the second study goal, GTV deviations computed from matching RCCT and RC-CBCT are compared to deviations computed from matching respiration-averaged images consisting of a CBCT reconstructed using all projections and an average-intensity-projection CT computed from the RCCT. Results: Of the eleven patients in the GTV-based systematic correction protocol, two required no correction, seven required a single correction, one required two corrections, and one required three corrections. Mean residual GTV deviation (3D distance) following GTV-based systematic correction (mean {+-} 1 standard deviation 4.8 {+-} 1.5 mm) is significantly lower than for systematic skeletal-based (6.5 {+-} 2.9 mm, p= 0.015), and weekly skeletal-based correction (7.2 {+-} 3.0 mm, p= 0.001), but is not significantly lower than daily skeletal-based correction (5.4 {+-} 2.6 mm, p= 0.34). In two cases, first-day CBCT images reveal tumor changes-one showing tumor growth, the other showing large tumor displacement-that are not readily observed in radiographs. Differences in computed GTV deviations between respiration-correlated and respiration-averaged images are 0.2 {+-} 1.8 mm in the superior-inferior direction and are of similar magnitude in the other directions. Conclusions: An off-line protocol to correct GTV-based systematic error in locally advanced lung tumor cases can be effective at reducing tumor deviations, although the findings need confirmation with larger patient statistics. In some cases, a single cone-beam CT can be useful for assessing tumor changes early in treatment, if more than a few days elapse between simulation and the start of treatment. Tumor deviations measured with respiration-averaged CT and CBCT images are consistent with those measured with respiration-correlated images; the respiration-averaged method is more easily implemented in the clinic.

60 APPLIED LIFE SCIENCES↗

Reduction of prostate intrafraction motion using gas-release rectal balloons

Purpose: To analyze prostate intrafraction motion using both non-gas-release (NGR) and gas-release (GR) rectal balloons and to evaluate the ability of GR rectal balloons to reduce prostate intrafraction motion. Methods: Twenty-nine patients with NGR rectal balloons and 29 patients with GR balloons were randomly selected from prostate patients treated with proton therapy at University of Florida Proton Therapy Institute (Jacksonville, FL). Their pretreatment and post-treatment orthogonal radiographs were analyzed, and both pretreatment setup residual error and intrafraction-motion data were obtained. Population histograms of intrafraction motion were plotted for both types of balloons. Population planning target-volume (PTV) margins were calculated with the van Herk formula of 2.5{Sigma}+ 0.7{sigma} to account for setup residual errors and intrafraction motion errors. Results: Pretreatment and post-treatment radiographs indicated that the use of gas-release rectal balloons reduced prostate intrafraction motion along superior-inferior (SI) and anterior-posterior (AP) directions. Similar patient setup residual errors were exhibited for both types of balloons. Gas-release rectal balloons resulted in PTV margin reductions from 3.9 to 2.8 mm in the SI direction, 3.1 to 1.8 mm in the AP direction, and an increase from 1.9 to 2.1 mm in the left-right direction. Conclusions: Prostate intrafraction motion is an important uncertainty source in radiotherapy after image-guided patient setup with online corrections. Compared to non-gas-release rectal balloons, gas-release balloons can reduce prostate intrafraction motion in the SI and AP directions caused by gas buildup.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A method of dose reconstruction for moving targets compatible with dynamic treatments

Purpose: To develop a method that allows a commercial treatment planning system (TPS) to perform accurate dose reconstruction for rigidly moving targets and to validate the method in phantom measurements for a range of treatments including intensity modulated radiation therapy (IMRT), volumetric arc therapy (VMAT), and dynamic multileaf collimator (DMLC) tracking. Methods: An in-house computer program was developed to manipulate Dicom treatment plans exported from a TPS (Eclipse, Varian Medical Systems) such that target motion during treatment delivery was incorporated into the plans. For each treatment, a motion including plan was generated by dividing the intratreatment target motion into 1 mm position bins and construct sub-beams that represented the parts of the treatment that were delivered, while the target was located within each position bin. For each sub-beam, the target shift was modeled by a corresponding isocenter shift. The motion incorporating Dicom plans were reimported into the TPS, where dose calculation resulted in motion including target dose distributions. For experimental validation of the dose reconstruction a thorax phantom with a moveable lung equivalent rod with a tumor insert of solid water was first CT scanned. The tumor insert was delineated as a gross tumor volume (GTV), and a planning target volume (PTV) was formed by adding margins. A conformal plan, two IMRT plans (step-and-shoot and sliding windows), and a VMAT plan were generated giving minimum target doses of 95% (GTV) and 67% (PTV) of the prescription dose (3 Gy). Two conformal fields with MLC leaves perpendicular and parallel to the tumor motion, respectively, were generated for DMLC tracking. All treatment plans were delivered to the thorax phantom without tumor motion and with a sinusoidal tumor motion. The two conformal fields were delivered with and without portal image guided DMLC tracking based on an embedded gold marker. The target dose distribution was measured with a radiochromic film in the moving rod and compared with the reconstructed doses using gamma tests. Results: Considerable interplay effects between machine motion and target motion were observed for the treatments without tracking. For nontracking experiments, the mean 2 mm/2% gamma pass rate over all investigated scenarios was 99.6% between calculated and measured doses. For tracking experiments, the mean gamma pass rate was 99.4%. Conclusions: A method for accurate dose reconstruction for moving targets with dynamic treatments was developed and experimentally validated in a variety of delivery scenarios. The method is suitable for integration into TPSs, e.g., for reconstruction of the dose delivered to moving tumors or calculation of target doses delivered with DMLC tracking.

60 APPLIED LIFE SCIENCES↗

Robust methods for automatic image-to-world registration in cone-beam CT interventional guidance

Purpose: Real-time surgical navigation relies on accurate image-to-world registration to align the coordinate systems of the image and patient. Conventional manual registration can present a workflow bottleneck and is prone to manual error and intraoperator variability. This work reports alternative means of automatic image-to-world registration, each method involving an automatic registration marker (ARM) used in conjunction with C-arm cone-beam CT (CBCT). The first involves a Known-Model registration method in which the ARM is a predefined tool, and the second is a Free-Form method in which the ARM is freely configurable. Methods: Studies were performed using a prototype C-arm for CBCT and a surgical tracking system. A simple ARM was designed with markers comprising a tungsten sphere within infrared reflectors to permit detection of markers in both x-ray projections and by an infrared tracker. The Known-Model method exercised a predefined specification of the ARM in combination with 3D-2D registration to estimate the transformation that yields the optimal match between forward projection of the ARM and the measured projection images. The Free-Form method localizes markers individually in projection data by a robust Hough transform approach extended from previous work, backprojected to 3D image coordinates based on C-arm geometric calibration. Image-domain point sets were transformed to world coordinates by rigid-body point-based registration. The robustness and registration accuracy of each method was tested in comparison to manual registration across a range of body sites (head, thorax, and abdomen) of interest in CBCT-guided surgery, including cases with interventional tools in the radiographic scene. Results: The automatic methods exhibited similar target registration error (TRE) and were comparable or superior to manual registration for placement of the ARM within {approx}200 mm of C-arm isocenter. Marker localization in projection data was robust across all anatomical sites, including challenging scenarios involving the presence of interventional tools. The reprojection error of marker localization was independent of the distance of the ARM from isocenter, and the overall TRE was dominated by the configuration of individual fiducials and distance from the target as predicted by theory. The median TRE increased with greater ARM-to-isocenter distance (e.g., for the Free-Form method, TRE increasing from 0.78 mm to 2.04 mm at distances of {approx}75 mm and 370 mm, respectively). The median TRE within {approx}200 mm distance was consistently lower than that of the manual method (TRE = 0.82 mm). Registration performance was independent of anatomical site (head, thorax, and abdomen). The Free-Form method demonstrated a statistically significant improvement (p= 0.0044) in reproducibility compared to manual registration (0.22 mm versus 0.30 mm, respectively). Conclusions: Automatic image-to-world registration methods demonstrate the potential for improved accuracy, reproducibility, and workflow in CBCT-guided procedures. A Free-Form method was shown to exhibit robustness against anatomical site, with comparable or improved TRE compared to manual registration. It was also comparable or superior in performance to a Known-Model method in which the ARM configuration is specified as a predefined tool, thereby allowing configuration of fiducials on the fly or attachment to the patient.

60 APPLIED LIFE SCIENCES↗

Positron emission mammography (PEM): Effect of activity concentration, object size, and object contrast on phantom lesion detection

Purpose: To characterize the relationship between lesion detection sensitivity and injected activity as a function of lesion size and contrast on the PEM (positron emission mammography) Flex Solo II scanner using phantom experiments. Methods: Phantom lesions (spheres 4, 8, 12, 16, and 20 mm diameter) were randomly located in uniform background. Sphere activity concentrations were 3 to 21 times the background activity concentration (BGc). BGc was a surrogate for injected activity; BGc ranged from 0.44-4.1 kBq/mL, corresponding to 46-400 MBq injections. Seven radiologists read 108 images containing zero, one, or two spheres. Readers used a 5-point confidence scale to score the presence of spheres. Results: Sensitivity was 100% for lesions {>=}12 mm under all conditions except for one 12 mm sphere with the lowest contrast and lowest BGc (60% sensitivity). Sensitivity was 100% for 8 mm spheres when either contrast or BGc was high, and 100% for 4 mm spheres only when both contrast and BGc were highest. Sphere contrast recovery coefficients (CRC) were 49%, 34%, 26%, 14%, and 2.8% for the largest to smallest spheres. Cumulative specificity was 98%. Conclusions: Phantom lesion detection sensitivity depends more on sphere size and contrast than on BGc. Detection sensitivity remained {>=}90% for injected activities as low as 100 MBq, for lesions {>=}8 mm. Low CRC in 4 mm objects results in moderate detection sensitivity even for 400 MBq injected activity, making it impractical to optimize injected activity for such lesions. Low CRC indicates that when lesions <8 mm are observed on PEM images they are highly tracer avid with greater potential of clinical significance. High specificity (98%) suggests that image statistical noise does not lead to false positive findings. These results apply to the 85 mm thick object used to obtain them; lesion detectability should be better (worse) for thinner (thicker) objects based on the reduced (increased) influence of photon attenuation.

60 APPLIED LIFE SCIENCES↗

Liver motion during cone beam computed tomography guided stereotactic body radiation therapy

Purpose: Understanding motion characteristics of liver such as, interfractional and intrafractional motion variability, difference in motion within different locations in the organ, and their complex relationship with the breathing cycles are particularly important for image-guided liver SBRT. The purpose of this study was to investigate such motion characteristics based on fiducial markers tracked with the x-ray projections of the CBCT scans, taken immediately prior to the treatments. Methods: Twenty liver SBRT patients were analyzed. Each patient had three fiducial markers (2 Multiplication-Sign 5-mm gold) percutaneously implanted around the gross tumor. The prescription ranged from 2 to 8 fractions per patient. The CBCT projections data for each fraction ({approx}650 projections/scan), for each patient, were analyzed and the 2D positions of the markers were extracted using an in-house algorithm. In total, >55 000 x-ray projections were analyzed from 85 CBCT scans. From the 2D extracted positions, a 3D motion trajectory of the markers was constructed, from each CBCT scans, resulting in left-right (LR), anterior-posterior (AP), and cranio-caudal (CC) location information of the markers with >55 000 data points. The authors then analyzed the interfraction and intrafraction liver motion variability, within different locations in the organ, and as a function of the breathing cycle. The authors also compared the motion characteristics against the planning 4DCT and the RPM Trade-Mark-Sign (Varian Medical Systems, Palo Alto, CA) breathing traces. Variations in the appropriate gating window (defined as the percent of the maximum range at which 50% of the marker positions are contained), between fractions were calculated as well. Results: The range of motion for the 20 patients were 3.0 {+-} 2.0 mm, 5.1 {+-} 3.1 mm, and 17.9 {+-} 5.1 mm in the planning 4DCT, and 2.8 {+-} 1.6 mm, 5.3 {+-} 3.1 mm, and 16.5 {+-} 5.7 mm in the treatment CBCT, for LR, AP, and CC directions, respectively. The range of respiratory period was 3.9 {+-} 0.7 and 4.2 {+-} 0.8 s during the 4DCT simulation and the CBCT scans, respectively. The authors found that breathing-induced AP and CC motions are highly correlated. That is, all markers moved cranially also moved posteriorly and vice versa, irrespective of the location. The LR motion had a more variable relationship with the AP/CC motions, and appeared random with respect to the location. That is, when the markers moved toward cranial-posterior direction, 58% of the markers moved to the patient-right, 22% of the markers moved to the patient-left, and 20% of the markers had minimal/none motion. The absolute difference in the motion magnitude between the markers, in different locations within the liver, had a positive correlation with the absolute distance between the markers (R{sup 2}= 0.69, linear-fit). The interfractional gating window varied significantly for some patients, with the largest having 29.4%-56.4% range between fractions. Conclusions: This study analyzed the liver motion characteristics of 20 patients undergoing SBRT. A large variation in motion was observed, interfractionally and intrafractionally, and that as the distance between the markers increased, the difference in the absolute range of motion also increased. This suggests that marker(s) in closest proximity to the target be used.

60 APPLIED LIFE SCIENCES↗

2D lag and signal nonlinearity correction in an amorphous silicon EPID and their impact on pretreatment dosimetric verification

Purpose: This investigation provides measurements of signal lag and nonlinearity separately for the Varian aS500 electronic portal imaging device (EPID), and an algorithm to correct for these effects in 2D; their potential impact on intensity modulated radiation therapy (IMRT) verification is also investigated. The authors quantify lag, as a function of both delivered monitor units (MU) and time, by using a range of MUs delivered at a clinically used rate of 400 MU/min. Explicit cumulative lag curves are thus determined for a range of MUs and times between the end of irradiation and the end of image acquisition. Signal nonlinearity is also investigated as a function of total MUs delivered. The family of cumulative lag curves and signal nonlinearity are then used to determine their effects on dynamic multileaf collimator (MLC) (IMRT) deliveries, and to correct for theses effects in 2D. Methods: Images acquired with an aS500 EPID and Varis Portal-Vision software were used to quantify detector lag and signal-nonlinearity. For the signal lag investigation, Portal-Vision's service monitor was used to acquire EPID images at a rate of 8 frames/s. The images were acquired during irradiation and 66 s thereafter, by inhibiting the M-holdoff-In signal of the Linac for a range of 4.5-198.5 MUs. Relative cumulative lag was calculated by integrating the EPID signal for a time after beam-off, and normalizing this to the integrated EPID signal accumulated during radiation. Signal nonlinearity was studied by acquiring 10 Multiplication-Sign 10 cm{sup 2} open-field EPID images in 'integrated image' mode for a range of 2-500 MUs, and normalized to the 100 MU case. All data were incorporated into in-house written software to create a 2D correction map for these effects, using the field's MLC file and a field-specific calculated 2D 'time-map,' which keeps track of the time elapsed from the last fluence delivered at each given point in the image to the end of the beam delivery. Results: Relative cumulative lag curves reveal that the lag alone can deviate the EPID's perceived dose by as large as 6% (1 MU delivery, 60 s postirradiation). For signal nonlinearity relative to 100 MU, EPID signals per MU of 0.84 and 1.01 were observed for 2 and 500 MUs, respectively. Correction maps were applied to a 1 cm sweeping-window 14 Multiplication-Sign 14 cm{sup 2} field and clinical head-and-neck IMRT field. A mean correction of 1.028 was implemented in the head-and-neck field, which significantly reduced lag-related asymmetries in the EPID images, and restored linearity to the EPID imager's dose response. Corrections made to the sweeping-field showed good agreement with the treatment planning system-predicted field, yielding an average percent difference of 0.05%{+-} 0.91%, compared to the -1.32%{+-} 1.02% before corrections, or 1.75%{+-} 1.04% when only a signal nonlinearity correction is made. Conclusions: Lag and signal-nonlinearity have been quantified for an aS500 EPID imager, and an effective 2D correction method has been developed which effectively removes nonlinearity and lag effects. Both of these effects were shown to negatively impact IMRT verifications. Especially fields that involve prolonged irradiation and small overall MUs should be corrected for in 2D.

60 APPLIED LIFE SCIENCES↗

Observation of super-resolution in digital breast tomosynthesis

Purpose: Digital breast tomosynthesis (DBT) is a 3D x-ray imaging modality in which tomographic sections of the breast are generated from a limited range of tube angles. Because oblique x-ray incidence shifts the image of an object in subpixel detector element increments with each increasing projection angle, it is demonstrated that DBT is capable of super-resolution (i.e., subpixel resolution). Methods: By convention, DBT reconstructions are performed on planes parallel to the breast support at various depths of the breast volume. In order for resolution in each reconstructed slice to be comparable to the detector, the pixel size should match that of the detector elements; hence, the highest frequency that can be resolved in the plane of reconstruction is the alias frequency of the detector. This study considers reconstruction grids with much smaller pixelation to visualize higher frequencies. For analytical proof of super-resolution, a theoretical framework is developed in which the reconstruction of a high frequency sinusoidal input is calculated using both simple backprojection (SBP) and filtered backprojection. To study the frequency spectrum of the reconstruction, its Fourier transform is also determined. The experimental feasibility of super-resolution was investigated by acquiring images of a bar pattern phantom with frequencies higher than the detector alias frequency. Results: Using analytical modeling, it is shown that the central projection cannot resolve frequencies exceeding the detector alias frequency. The Fourier transform of the central projection is maximized at a lower frequency than the input as evidence of aliasing. By contrast, SBP reconstruction can resolve the input, and its Fourier transform is correctly maximized at the input frequency. Incorporating filters into the reconstruction smoothens pixelation artifacts in the spatial domain and reduces spectral leakage in the Fourier domain. It is also demonstrated that the existence of super-resolution is dependent on position in the reconstruction and on the directionality of the input frequency. Consistent with the analytical results, experimental reconstructions of bar patterns showed visibility of frequencies greater than the detector alias frequency. Super-resolution was present at positions predicted from analytical modeling. Conclusions: This work demonstrates the existence of super-resolution in DBT. Super-resolution has the potential to impact the visualization of fine structural details in the breast, such as microcalcifications and other subtle signs of cancer.

60 APPLIED LIFE SCIENCES↗

Automated size-specific CT dose monitoring program: Assessing variability in CT dose

Purpose: The potential health risks associated with low levels of ionizing radiation have created a movement in the radiology community to optimize computed tomography (CT) imaging protocols to use the lowest radiation dose possible without compromising the diagnostic usefulness of the images. Despite efforts to use appropriate and consistent radiation doses, studies suggest that a great deal of variability in radiation dose exists both within and between institutions for CT imaging. In this context, the authors have developed an automated size-specific radiation dose monitoring program for CT and used this program to assess variability in size-adjusted effective dose from CT imaging. Methods: The authors radiation dose monitoring program operates on an independent health insurance portability and accountability act compliant dosimetry server. Digital imaging and communication in medicine routing software is used to isolate dose report screen captures and scout images for all incoming CT studies. Effective dose conversion factors (k-factors) are determined based on the protocol and optical character recognition is used to extract the CT dose index and dose-length product. The patient's thickness is obtained by applying an adaptive thresholding algorithm to the scout images and is used to calculate the size-adjusted effective dose (ED{sub adj}). The radiation dose monitoring program was used to collect data on 6351 CT studies from three scanner models (GE Lightspeed Pro 16, GE Lightspeed VCT, and GE Definition CT750 HD) and two institutions over a one-month period and to analyze the variability in ED{sub adj} between scanner models and across institutions. Results: No significant difference was found between computer measurements of patient thickness and observer measurements (p= 0.17), and the average difference between the two methods was less than 4%. Applying the size correction resulted in ED{sub adj} that differed by up to 44% from effective dose estimates that were not adjusted by patient size. Additionally, considerable differences were noted in ED{sub adj} distributions between scanners, with scanners employing iterative reconstruction exhibiting significantly lower ED{sub adj} (range: 9%-64%). Finally, a significant difference (up to 59%) in ED{sub adj} distributions was observed between institutions, indicating the potential for dose reduction. Conclusions: The authors developed a robust automated size-specific radiation dose monitoring program for CT. Using this program, significant differences in ED{sub adj} were observed between scanner models and across institutions. This new dose monitoring program offers a unique tool for improving quality assurance and standardization both within and across institutions.

60 APPLIED LIFE SCIENCES↗

A fully 3D approach for metal artifact reduction in computed tomography

Purpose: In computed tomography imaging metal objects in the region of interest introduce inconsistencies during data acquisition. Reconstructing these data leads to an image in spatial domain including star-shaped or stripe-like artifacts. In order to enhance the quality of the resulting image the influence of the metal objects can be reduced. Here, a metal artifact reduction (MAR) approach is proposed that is based on a recomputation of the inconsistent projection data using a fully three-dimensional Fourier-based interpolation. The success of the projection space restoration depends sensitively on a sensible continuation of neighboring structures into the recomputed area. Fortunately, structural information of the entire data is inherently included in the Fourier space of the data. This can be used for a reasonable recomputation of the inconsistent projection data. Methods: The key step of the proposed MAR strategy is the recomputation of the inconsistent projection data based on an interpolation using nonequispaced fast Fourier transforms (NFFT). The NFFT interpolation can be applied in arbitrary dimension. The approach overcomes the problem of adequate neighborhood definitions on irregular grids, since this is inherently given through the usage of higher dimensional Fourier transforms. Here, applications up to the third interpolation dimension are presented and validated. Furthermore, prior knowledge may be included by an appropriate damping of the transform during the interpolation step. This MAR method is applicable on each angular view of a detector row, on two-dimensional projection data as well as on three-dimensional projection data, e.g., a set of sequential acquisitions at different spatial positions, projection data of a spiral acquisition, or cone-beam projection data. Results: Results of the novel MAR scheme based on one-, two-, and three-dimensional NFFT interpolations are presented. All results are compared in projection data space and spatial domain with the well-known one-dimensional linear interpolation strategy. Conclusions: In conclusion, it is recommended to include as much spatial information into the recomputation step as possible. This is realized by increasing the dimension of the NFFT. The resulting image quality can be enhanced considerably.

60 APPLIED LIFE SCIENCES↗

Fast transit portal dosimetry using density-scaled layer modeling of aSi-based electronic portal imaging device and Monte Carlo method

Purpose: Fast and accurate transit portal dosimetry was investigated by developing a density-scaled layer model of electronic portal imaging device (EPID) and applying it to a clinical environment. Methods: The model was developed for fast Monte Carlo dose calculation. The model was validated through comparison with measurements of dose on EPID using first open beams of varying field sizes under a 20-cm-thick flat phantom. After this basic validation, the model was further tested by applying it to transit dosimetry and dose reconstruction that employed our predetermined dose-response-based algorithm developed earlier. The application employed clinical intensity-modulated beams irradiated on a Rando phantom. The clinical beams were obtained through planning on pelvic regions of the Rando phantom simulating prostate and large pelvis intensity modulated radiation therapy. To enhance agreement between calculations and measurements of dose near penumbral regions, convolution conversion of acquired EPID images was alternatively used. In addition, thickness-dependent image-to-dose calibration factors were generated through measurements of image and calculations of dose in EPID through flat phantoms of various thicknesses. The factors were used to convert acquired images in EPID into dose. Results: For open beam measurements, the model showed agreement with measurements in dose difference better than 2% across open fields. For tests with a Rando phantom, the transit dosimetry measurements were compared with forwardly calculated doses in EPID showing gamma pass rates between 90.8% and 98.8% given 4.5 mm distance-to-agreement (DTA) and 3% dose difference (DD) for all individual beams tried in this study. The reconstructed dose in the phantom was compared with forwardly calculated doses showing pass rates between 93.3% and 100% in isocentric perpendicular planes to the beam direction given 3 mm DTA and 3% DD for all beams. On isocentric axial planes, the pass rates varied between 95.8% and 99.9% for all individual beams and they were 98.2% and 99.9% for the composite beams of the small and large pelvis cases, respectively. Three-dimensional gamma pass rates were 99.0% and 96.4% for the small and large pelvis cases, respectively. Conclusions: The layer model of EPID built for Monte Carlo calculations offered fast (less than 1 min) and accurate calculation for transit dosimety and dose reconstruction.

60 APPLIED LIFE SCIENCES↗

Dedicated breast CT: Fibroglandular volume measurements in a diagnostic population

Purpose: To determine the mean and range of volumetric glandular fraction (VGF) of the breast in a diagnostic population using a high-resolution flat-panel cone-beam dedicated breast CT system. This information is important for Monte Carlo-based estimation of normalized glandular dose coefficients and for investigating the dependence of VGF on breast dimensions, race, and pathology. Methods: Image data from a clinical trial investigating the role of dedicated breast CT that enrolled 150 women were retrospectively analyzed to determine the VGF. The study was conducted in adherence to a protocol approved by the institutional human subjects review boards and written informed consent was obtained from all study participants. All participants in the study were assigned BI-RADS{sup Registered-Sign} 4 or 5 as per the American College of Radiology assessment categories after standard diagnostic work-up and underwent dedicated breast CT exam prior to biopsy. A Gaussian-kernel based fuzzy c-means algorithm was used to partition the breast CT images into adipose and fibroglandular tissue after segmenting the skin. Upon determination of the accuracy of the algorithm with a phantom, it was applied to 137 breast CT volumes from 136 women. VGF was determined for each breast and the mean and range were determined. Pathology results with classification as benign, malignant, and hyperplasia were available for 132 women, and were used to investigate if the distributions of VGF varied with pathology. Results: The algorithm was accurate to within {+-}1.9% in determining the volume of an irregular shaped phantom. The study mean ({+-} inter-breast SD) for the VGF was 0.172 {+-} 0.142 (range: 0.012-0.719). VGF was found to be negatively correlated with age, breast dimensions (chest-wall to nipple length, pectoralis to nipple length, and effective diameter at chest-wall), and total breast volume, and positively correlated with fibroglandular volume. Based on pathology, pairwise statistical analysis (Mann-Whitney test) indicated that at the 0.05 significance level, there was no significant difference in distributions of VGF without adjustment for age between malignant and nonmalignant breasts (p= 0.41). Pairwise comparisons of the distributions of VGF in increasing order of mammographic breast density indicated all comparisons were statistically significant (p < 0.002). Conclusions: This study used a different clinical prototype breast CT system than that in previous studies to image subjects from a different geographical region, and used a different algorithm for analysis of image data. The mean VGF estimated from this study is within the range reported in previous studies, indicating that the choice of 50% glandular weight fraction to represent an average breast for Monte Carlo-based estimation of normalized glandular dose coefficients in mammography needs revising. In the study, the distributions of VGF did not differ significantly with pathology.

60 APPLIED LIFE SCIENCES↗

A patient alignment solution for lung SBRT setups based on a deformable registration technique

Purpose: In this work, the authors propose a novel registration strategy for translation-only correction scenarios of lung stereotactic body radiation therapy setups, which can achieve optimal dose coverage for tumors as well as preserve the consistency of registrations with minimal human interference. Methods: The proposed solution (centroid-to-centroidor CTC solution) uses the average four-dimensional CT (A4DCT) as the reference CT. The cone-beam CT (CBCT) is deformed to acquire a new centroid for the internal target volume (ITV) on the CBCT. The registration is then accomplished by simply aligning the centroids of the ITVs between the A4DCT and the CBCT. Sixty-seven cases using 64 patients (each case is associated with separate isocenters) have been investigated with the CTC method and compared with the conventional gray-value (G) mode and bone (B) mode registration methods. Dosimetric effects among the tree methods were demonstrated by 18 selected cases. The uncertainty of the CTC method has also been studied. Results: The registration results demonstrate the superiority of the CTC method over the other two methods. The differences in the D99 and D95 ITV dose coverage between the CTC method and the original plan is small (within 5%) for all of the selected cases except for one for which the tumor presented significant growth during the period between the CT scan and the treatment. Meanwhile, the dose coverage differences between the original plan and the registration results using either the B or G method are significant, as tumor positions varied dramatically, relative to the rib cage, from their positions on the original CT. The largest differences between the D99 and D95 dose coverage of the ITV using the B or G method versus the original plan are as high as 50%. The D20 differences between any of the methods versus the original plan are all less than 2%. Conclusions: The CTC method can generate optimal dose coverage to tumors with much better consistency compared with either the G or B method, and it is especially useful when the tumor position varies greatly from its position on the original CT, relative to the rib cage.

60 APPLIED LIFE SCIENCES↗

Accuracy of volume measurement using 3D ultrasound and development of CT-3D US image fusion algorithm for prostate cancer radiotherapy

Purpose: To evaluate the accuracy of measuring volumes using three-dimensional ultrasound (3D US), and to verify the feasibility of the replacement of CT-MR fusion images with CT-3D US in radiotherapy treatment planning. Methods: Phantoms, consisting of water, contrast agent, and agarose, were manufactured. The volume was measured using 3D US, CT, and MR devices. A CT-3D US and MR-3D US image fusion software was developed using the Insight Toolkit library in order to acquire three-dimensional fusion images. The quality of the image fusion was evaluated using metric value and fusion images. Results: Volume measurement, using 3D US, shows a 2.8 {+-} 1.5% error, 4.4 {+-} 3.0% error for CT, and 3.1 {+-} 2.0% error for MR. The results imply that volume measurement using the 3D US devices has a similar accuracy level to that of CT and MR. Three-dimensional image fusion of CT-3D US and MR-3D US was successfully performed using phantom images. Moreover, MR-3D US image fusion was performed using human bladder images. Conclusions: 3D US could be used in the volume measurement of human bladders and prostates. CT-3D US image fusion could be used in monitoring the target position in each fraction of external beam radiation therapy. Moreover, the feasibility of replacing the CT-MR image fusion to the CT-3D US in radiotherapy treatment planning was verified.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A digitally reconstructed radiograph algorithm calculated from first principles

Purpose: To develop an algorithm for computing realistic digitally reconstructed radiographs (DRRs) that match real cone-beam CT (CBCT) projections with no artificial adjustments. Methods: The authors used measured attenuation data from cone-beam CT projection radiographs of different materials to obtain a function to convert CT number to linear attenuation coefficient (LAC). The effects of scatter, beam hardening, and veiling glare were first removed from the attenuation data. Using this conversion function the authors calculated the line integral of LAC through a CT along rays connecting the radiation source and detector pixels with a ray-tracing algorithm, producing raw DRRs. The effects of scatter, beam hardening, and veiling glare were then included in the DRRs through postprocessing. Results: The authors compared actual CBCT projections to DRRs produced with all corrections (scatter, beam hardening, and veiling glare) and to uncorrected DRRs. Algorithm accuracy was assessed through visual comparison of projections and DRRs, pixel intensity comparisons, intensity histogram comparisons, and correlation plots of DRR-to-projection pixel intensities. In general, the fully corrected algorithm provided a small but nontrivial improvement in accuracy over the uncorrected algorithm. The authors also investigated both measurement- and computation-based methods for determining the beam hardening correction, and found the computation-based method to be superior, as it accounted for nonuniform bowtie filter thickness. The authors benchmarked the algorithm for speed and found that it produced DRRs in about 0.35 s for full detector and CT resolution at a ray step-size of 0.5 mm. Conclusions: The authors have demonstrated a DRR algorithm calculated from first principles that accounts for scatter, beam hardening, and veiling glare in order to produce accurate DRRs. The algorithm is computationally efficient, making it a good candidate for iterative CT reconstruction techniques that require a data fidelity term based on the matching of DRRs and projections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GPU-based relative fuzzy connectedness image segmentation

Purpose:Recently, clinical radiological research and practice are becoming increasingly quantitative. Further, images continue to increase in size and volume. For quantitative radiology to become practical, it is crucial that image segmentation algorithms and their implementations are rapid and yield practical run time on very large data sets. The purpose of this paper is to present a parallel version of an algorithm that belongs to the family of fuzzy connectedness (FC) algorithms, to achieve an interactive speed for segmenting large medical image data sets. Methods: The most common FC segmentations, optimizing an Script-Small-L {sub {infinity}}-based energy, are known as relative fuzzy connectedness (RFC) and iterative relative fuzzy connectedness (IRFC). Both RFC and IRFC objects (of which IRFC contains RFC) can be found via linear time algorithms, linear with respect to the image size. The new algorithm, P-ORFC (for parallel optimal RFC), which is implemented by using NVIDIA's Compute Unified Device Architecture (CUDA) platform, considerably improves the computational speed of the above mentioned CPU based IRFC algorithm. Results: Experiments based on four data sets of small, medium, large, and super data size, achieved speedup factors of 32.8 Multiplication-Sign , 22.9 Multiplication-Sign , 20.9 Multiplication-Sign , and 17.5 Multiplication-Sign , correspondingly, on the NVIDIA Tesla C1060 platform. Although the output of P-ORFC need not precisely match that of IRFC output, it is very close to it and, as the authors prove, always lies between the RFC and IRFC objects. Conclusions: A parallel version of a top-of-the-line algorithm in the family of FC has been developed on the NVIDIA GPUs. An interactive speed of segmentation has been achieved, even for the largest medical image data set. Such GPU implementations may play a crucial role in automatic anatomy recognition in clinical radiology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗