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 901 records · Page 50

Generation of a suite of 3D computer-generated breast phantoms from a limited set of human subject data

Purpose: The authors previously reported on a three-dimensional computer-generated breast phantom, based on empirical human image data, including a realistic finite-element based compression model that was capable of simulating multimodality imaging data. The computerized breast phantoms are a hybrid of two phantom generation techniques, combining empirical breast CT (bCT) data with flexible computer graphics techniques. However, to date, these phantoms have been based on single human subjects. In this paper, the authors report on a new method to generate multiple phantoms, simulating additional subjects from the limited set of original dedicated breast CT data. The authors developed an image morphing technique to construct new phantoms by gradually transitioning between two human subject datasets, with the potential to generate hundreds of additional pseudoindependent phantoms from the limited bCT cases. The authors conducted a preliminary subjective assessment with a limited number of observers (n= 4) to illustrate how realistic the simulated images generated with the pseudoindependent phantoms appeared. Methods: Several mesh-based geometric transformations were developed to generate distorted breast datasets from the original human subject data. Segmented bCT data from two different human subjects were used as the 'base' and 'target' for morphing. Several combinations of transformations were applied to morph between the 'base' and 'target' datasets such as changing the breast shape, rotating the glandular data, and changing the distribution of the glandular tissue. Following the morphing, regions of skin and fat were assigned to the morphed dataset in order to appropriately assign mechanical properties during the compression simulation. The resulting morphed breast was compressed using a finite element algorithm and simulated mammograms were generated using techniques described previously. Sixty-two simulated mammograms, generated from morphing three human subject datasets, were used in a preliminary observer evaluation where four board certified breast radiologists with varying amounts of experience ranked the level of realism (from 1 ='fake' to 10 ='real') of the simulated images. Results: The morphing technique was able to successfully generate new and unique morphed datasets from the original human subject data. The radiologists evaluated the realism of simulated mammograms generated from the morphed and unmorphed human subject datasets and scored the realism with an average ranking of 5.87 {+-} 1.99, confirming that overall the phantom image datasets appeared more 'real' than 'fake.' Moreover, there was not a significant difference (p > 0.1) between the realism of the unmorphed datasets (6.0 {+-} 1.95) compared to the morphed datasets (5.86 {+-} 1.99). Three of the four observers had overall average rankings of 6.89 {+-} 0.89, 6.9 {+-} 1.24, 6.76 {+-} 1.22, whereas the fourth observer ranked them noticeably lower at 2.94 {+-} 0.7. Conclusions: This work presents a technique that can be used to generate a suite of realistic computerized breast phantoms from a limited number of human subjects. This suite of flexible breast phantoms can be used for multimodality imaging research to provide a known truth while concurrently producing realistic simulated imaging data.

60 APPLIED LIFE SCIENCES↗

High spatiotemporal resolution measurement of regional lung air volumes from 2D phase contrast x-ray images

Purpose: Described herein is a new technique for measuring regional lung air volumes from two-dimensional propagation-based phase contrast x-ray (PBI) images at very high spatial and temporal resolution. Phase contrast dramatically increases lung visibility and the outlined volumetric reconstruction technique quantifies dynamic changes in respiratory function. These methods can be used for assessing pulmonary disease and injury and for optimizing mechanical ventilation techniques for preterm infants using animal models. Methods: The volumetric reconstruction combines the algorithms of temporal subtraction and single image phase retrieval (SIPR) to isolate the image of the lungs from the thoracic cage in order to measure regional lung air volumes. The SIPR algorithm was used to recover the change in projected thickness of the lungs on a pixel-by-pixel basis (pixel dimensions {approx}16.2 {mu}m). The technique has been validated using numerical simulation and compared results of measuring regional lung air volumes with and without the use of temporal subtraction for removing the thoracic cage. To test this approach, a series of PBI images of newborn rabbit pups mechanically ventilated at different frequencies was employed. Results: Regional lung air volumes measured from PBI images of newborn rabbit pups showed on average an improvement of at least 20% in 16% of pixels within the lungs in comparison to that measured without the use of temporal subtraction. The majority of pixels that showed an improvement was found to be in regions occupied by bone. Applying the volumetric technique to sequences of PBI images of newborn rabbit pups, it is shown that lung aeration at birth can be highly heterogeneous. Conclusions: This paper presents an image segmentation technique based on temporal subtraction that has successfully been used to isolate the lungs from PBI chest images, allowing the change in lung air volume to be measured over regions as small as the pixel size. Using this technique, it is possible to measure changes in regional lung volume at high spatial and temporal resolution during breathing at much lower x-ray dose than would be required using computed tomography.

60 APPLIED LIFE SCIENCES↗

The feasibility of a regional CTDI{sub vol} to estimate organ dose from tube current modulated CT exams

Purpose: In AAPM Task Group 204, the size-specific dose estimate (SSDE) was developed by providing size adjustment factors which are applied to the Computed Tomography (CT) standardized dose metric, CTDI{sub vol}. However, that work focused on fixed tube current scans and did not specifically address tube current modulation (TCM) scans, which are currently the majority of clinical scans performed. The purpose of this study was to extend the SSDE concept to account for TCM by investigating the feasibility of using anatomic and organ specific regions of scanner output to improve accuracy of dose estimates. Methods: Thirty-nine adult abdomen/pelvis and 32 chest scans from clinically indicated CT exams acquired on a multidetector CT using TCM were obtained with Institutional Review Board approval for generating voxelized models. Along with image data, raw projection data were obtained to extract TCM functions for use in Monte Carlo simulations. Patient size was calculated using the effective diameter described in TG 204. In addition, the scanner-reported CTDI{sub vol} (CTDI{sub vol,global}) was obtained for each patient, which is based on the average tube current across the entire scan. For the abdomen/pelvis scans, liver, spleen, and kidneys were manually segmented from the patient datasets; for the chest scans, lungs and for female models only, glandular breast tissue were segmented. For each patient organ doses were estimated using Monte Carlo Methods. To investigate the utility of regional measures of scanner output, regional and organ anatomic boundaries were identified from image data and used to calculate regional and organ-specific average tube current values. From these regional and organ-specific averages, CTDI{sub vol} values, referred to as regional and organ-specific CTDI{sub vol}, were calculated for each patient. Using an approach similar to TG 204, all CTDI{sub vol} values were used to normalize simulated organ doses; and the ability of each normalized dose to correlate with patient size was investigated. Results: For all five organs, the correlations with patient size increased when organ doses were normalized by regional and organ-specific CTDI{sub vol} values. For example, when estimating dose to the liver, CTDI{sub vol,global} yielded a R{sup 2} value of 0.26, which improved to 0.77 and 0.86, when using the regional and organ-specific CTDI{sub vol} for abdomen and liver, respectively. For breast dose, the global CTDI{sub vol} yielded a R{sup 2} value of 0.08, which improved to 0.58 and 0.83, when using the regional and organ-specific CTDI{sub vol} for chest and breasts, respectively. The R{sup 2} values also increased once the thoracic models were separated for the analysis into females and males, indicating differences between genders in this region not explained by a simple measure of effective diameter. Conclusions: This work demonstrated the utility of regional and organ-specific CTDI{sub vol} as normalization factors when using TCM. It was demonstrated that CTDI{sub vol,global} is not an effective normalization factor in TCM exams where attenuation (and therefore tube current) varies considerably throughout the scan, such as abdomen/pelvis and even thorax. These exams can be more accurately assessed for dose using regional CTDI{sub vol} descriptors that account for local variations in scanner output present when TCM is employed.

60 APPLIED LIFE SCIENCES↗

A virtual trial framework for quantifying the detectability of masses in breast tomosynthesis projection data

Purpose: Digital breast tomosynthesis (DBT) is a promising breast cancer screening tool that has already begun making inroads into clinical practice. However, there is ongoing debate over how to quantitatively evaluate and optimize these systems, because different definitions of image quality can lead to different optimal design strategies. Powerful and accurate tools are desired to extend our understanding of DBT system optimization and validate published design principles. Methods: The authors developed a virtual trial framework for task-specific DBT assessment that uses digital phantoms, open-source x-ray transport codes, and a projection-space, spatial-domain observer model for quantitative system evaluation. The authors considered evaluation of reconstruction algorithms as a separate problem and focused on the information content in the raw, unfiltered projection images. Specifically, the authors investigated the effects of scan angle and number of angular projections on detectability of a small (3 mm diameter) signal embedded in randomly-varying anatomical backgrounds. Detectability was measured by the area under the receiver-operating characteristic curve (AUC). Experiments were repeated for three test cases where the detectability-limiting factor was anatomical variability, quantum noise, or electronic noise. The authors also juxtaposed the virtual trial framework with other published studies to illustrate its advantages and disadvantages. Results: The large number of variables in a virtual DBT study make it difficult to directly compare different authors' results, so each result must be interpreted within the context of the specific virtual trial framework. The following results apply to 25% density phantoms with 5.15 cm compressed thickness and 500 {mu}m{sup 3} voxels (larger 500 {mu}m{sup 2} detector pixels were used to avoid voxel-edge artifacts): 1. For raw, unfiltered projection images in the anatomical-variability-limited regime, AUC appeared to remain constant or increase slightly with scan angle. 2. In the same regime, when the authors fixed the scan angle, AUC increased asymptotically with the number of projections. The threshold number of projections for asymptotic AUC performance depended on the scan angle. In the quantum- and electronic-noise dominant regimes, AUC behaviors as a function of scan angle and number of projections sometimes differed from the anatomy-limited regime. For example, with a fixed scan angle, AUC generally decreased with the number of projections in the electronic-noise dominant regime. These results are intended to demonstrate the capabilities of the virtual trial framework, not to be used as optimization rules for DBT. Conclusions: The authors have demonstrated a novel simulation framework and tools for evaluating DBT systems in an objective, task-specific manner. This framework facilitates further investigation of image quality tradeoffs in DBT.

60 APPLIED LIFE SCIENCES↗

Upgrade and benchmarking of a 4D treatment planning system for scanned ion beam therapy

Purpose: Upgrade and benchmarking of a research 4D treatment planning system (4DTPS) suitable for realistic patient treatment planning and treatment simulations taking into account specific requirements for scanned ion beam therapy, i.e., modeling of dose heterogeneities due to interplay effects and range changes caused by patient motion and dynamic beam delivery. Methods: The 4DTPS integrates data interfaces to 4D computed tomography (4DCT), deformable image registration and clinically used motion monitoring devices. The authors implemented a novel data model for 4D image segmentation using Boolean mask volume datasets and developed an algorithm propagating a manually contoured reference contour dataset to all 4DCT phases. They further included detailed treatment simulation and dose reconstruction functionality, based on the irregular patient motion and the temporal structure of the beam delivery. The treatment simulation functionality was validated against experimental data from irradiation of moving radiographic films in air, 3D moving ionization chambers in a water phantom, and moving cells in a biological phantom with a scanned carbon ion beam. The performance of the program was compared to results obtained with predecessor programs. Results: The measured optical density distributions of the radiographic films were reproduced by the simulations to (-2 {+-} 12)%. Compared to earlier versions of the 4DTPS, the mean agreement improved by 2%, standard deviations were reduced by 7%. The simulated dose to the moving ionization chambers in water showed an agreement with the measured dose of (-1 {+-} 4)% for the typical beam configuration. The mean deviation of the simulated from the measured biologically effective dose determined via cell survival was (617 {+-} 538) mGy relative biological effectiveness corresponding to (10 {+-} 9)%. Conclusions: The authors developed a research 4DTPS suitable for realistic treatment planning on patient data and capable of simulating dose delivery to a moving patient geometry for scanned ion beams. The accuracy and reliability of treatment simulations improved considerably with respect to earlier versions of the 4DTPS.

60 APPLIED LIFE SCIENCES↗

A real-time regional adaptive exposure method for saving dose-area product in x-ray fluoroscopy

Purpose: Reduction of radiation dose in x-ray imaging has been recognized as a high priority in the medical community. Here the authors show that a regional adaptive exposure method can reduce dose-area product (DAP) in x-ray fluoroscopy. The authors' method is particularly geared toward providing dose savings for the pediatric population. Methods: The scanning beam digital x-ray system uses a large-area x-ray source with 8000 focal spots in combination with a small photon-counting detector. An imaging frame is obtained by acquiring and reconstructing up to 8000 detector images, each viewing only a small portion of the patient. Regional adaptive exposure was implemented by varying the exposure of the detector images depending on the local opacity of the object. A family of phantoms ranging in size from infant to obese adult was imaged in anteroposterior view with and without adaptive exposure. The DAP delivered to each phantom was measured in each case, and noise performance was compared by generating noise arrays to represent regional noise in the images. These noise arrays were generated by dividing the image into regions of about 6 mm{sup 2}, calculating the relative noise in each region, and placing the relative noise value of each region in a one-dimensional array (noise array) sorted from highest to lowest. Dose-area product savings were calculated as the difference between the ratio of DAP with adaptive exposure to DAP without adaptive exposure. The authors modified this value by a correction factor that matches the noise arrays where relative noise is the highest to report a final dose-area product savings. Results: The average dose-area product saving across the phantom family was (42 {+-} 8)% with the highest dose-area product saving in the child-sized phantom (50%) and the lowest in the phantom mimicking an obese adult (23%). Conclusions: Phantom measurements indicate that a regional adaptive exposure method can produce large DAP savings without compromising the noise performance in the image regions with highest noise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dosimetric comparison of patient setup strategies in stereotactic body radiation therapy for lung cancer

Purpose: In this work, the authors retrospectively compared the accumulated dose over the treatment course for stereotactic body radiation therapy (SBRT) of lung cancer for three patient setup strategies. Methods: Ten patients who underwent lung SBRT were selected for this study. At each fraction, patients were immobilized using a vacuum cushion and were CT scanned. Treatment plans were performed on the simulation CT. The planning target volume (PTV) was created by adding a 5-mm uniform margin to the internal target volume derived from the 4DCT. All plans were normalized such that 99% of the PTV received 60 Gy. The plan parameters were copied onto the daily CT images for dose recalculation under three setup scenarios: skin marker, bony structure, and soft tissue based alignments. The accumulated dose was calculated by summing the dose at each fraction along the trajectory of a voxel over the treatment course through deformable image registration of each CT with the planning CT. The accumulated doses were analyzed for the comparison of setup accuracy. Results: The tumor volume receiving 60 Gy was 91.7 {+-} 17.9%, 74.1 {+-} 39.1%, and 99.6 {+-} 1.3% for setup using skin marks, bony structures, and soft tissue, respectively. The isodose line covering 100% of the GTV was 55.5 {+-} 7.1, 42.1 {+-} 16.0, and 64.3 {+-} 7.1 Gy, respectively. The corresponding average biologically effective dose of the tumor was 237.3 {+-} 29.4, 207.4 {+-} 61.2, and 258.3 {+-} 17.7 Gy, respectively. The differences in lung biologically effective dose, mean dose, and V20 between the setup scenarios were insignificant. Conclusions: The authors' results suggest that skin marks and bony structure are insufficient for aligning patients in lung SBRT. Soft tissue based alignment is needed to match the prescribed dose delivered to the tumors.

60 APPLIED LIFE SCIENCES↗

Toward correcting drift in target position during radiotherapy via computer-controlled couch adjustments on a programmable Linac

Purpose: Real-time tracking of respiratory target motion during radiation therapy is technically challenging, owing to rapid and possibly irregular breathing variations. The authors report on a method to predict and correct respiration-averaged drift in target position by means of couch adjustments on an accelerator equipped with such capability. Methods: Dose delivery is broken up into a sequence of 10 s field segments, each followed by a couch adjustment based on analysis of breathing motion from an external monitor as a surrogate of internal target motion. Signal averaging over three respiratory cycles yields a baseline representing target drift. A Kalman filter predicts the baseline position 5 s in advance, for determination of the couch correction. The method's feasibility is tested with a motion phantom programmed according to previously recorded patient signals. Computed couch corrections are preprogrammed into a research mode of an accelerator capable of computer-controlled couch translations synchronized with the motion phantom. The method's performance is evaluated with five cases recorded during hypofractionated treatment and five from respiration-correlated CT simulation, using a root-mean-squared deviation (RMSD) of the baseline from the treatment planned position. Results: RMSD is reduced in all 10 cases, from a mean of 4.9 mm (range 2.7-9.4 mm) before correction to 1.7 mm (range 0.7-2.3 mm) after correction. Treatment time is increased {approx}5% relative to that for no corrections. Conclusions: This work illustrates the potential for reduction in baseline respiratory drift with periodic adjustments in couch position during treatment. Future treatment machine capabilities will enable the use of 'on-the-fly' couch adjustments during treatment.

60 APPLIED LIFE SCIENCES↗

Resolution enhancement of lung 4D-CT data using multiscale interphase iterative nonlocal means

Purpose: Four-dimensional computer tomography (4D-CT) has been widely used in lung cancer radiotherapy due to its capability in providing important tumor motion information. However, the prolonged scanning duration required by 4D-CT causes considerable increase in radiation dose. To minimize the radiation-related health risk, radiation dose is often reduced at the expense of interslice spatial resolution. However, inadequate resolution in 4D-CT causes artifacts and increases uncertainty in tumor localization, which eventually results in extra damages of healthy tissues during radiotherapy. In this paper, the authors propose a novel postprocessing algorithm to enhance the resolution of lung 4D-CT data. Methods: The authors' premise is that anatomical information missing in one phase can be recovered from the complementary information embedded in other phases. The authors employ a patch-based mechanism to propagate information across phases for the reconstruction of intermediate slices in the longitudinal direction, where resolution is normally the lowest. Specifically, the structurally matching and spatially nearby patches are combined for reconstruction of each patch. For greater sensitivity to anatomical details, the authors employ a quad-tree technique to adaptively partition the image for more fine-grained refinement. The authors further devise an iterative strategy for significant enhancement of anatomical details. Results: The authors evaluated their algorithm using a publicly available lung data that consist of 10 4D-CT cases. The authors' algorithm gives very promising results with significantly enhanced image structures and much less artifacts. Quantitative analysis shows that the authors' algorithm increases peak signal-to-noise ratio by 3-4 dB and the structural similarity index by 3%-5% when compared with the standard interpolation-based algorithms. Conclusions: The authors have developed a new algorithm to improve the resolution of 4D-CT. It outperforms the conventional interpolation-based approaches by producing images with the markedly improved structural clarity and greatly reduced artifacts.

60 APPLIED LIFE SCIENCES↗

A multimodality vascular imaging phantom of an abdominal aortic aneurysm with a visible thrombus

Purpose: With the continuous development of new stent grafts and implantation techniques, it has now become technically feasible to treat abdominal aortic aneurysms (AAA) with challenging anatomy using endovascular repair with standard, fenestrated, or branched stent-grafts. In vitro experimentations are very useful to improve stent-graft design and conformability or imaging guidance for stent-graft delivery or follow-up. Vascular replicas also help to better understand the limitation of endovascular approaches in challenging anatomy and possibly improve surgical planning or training by practicing high risk clinical procedures in the laboratory to improve outcomes in the operating room. Most AAA phantoms available have a very basic anatomy, which is not representative of the clinical reality. This paper presents a method of fabrication of a realistic AAA phantom with a visible thrombus, as well as some mechanical properties characterizing such phantom. Methods: A realistic AAA geometry replica of a real patient anatomy taken from a multidetector computed tomography (CT) scan was manufactured. To demonstrate the multimodality imaging capability of this new phantom with a thrombus visible in magnetic resonance (MR) angiography, CT angiography (CTA), digital subtraction angiography (DSA), and ultrasound, image acquisitions with all these modalities were performed by using standard clinical protocols. Potential use of this phantom for stent deployment was also tested. A rheometer allowed defining hyperelastic and viscoelastic properties of phantom materials. Results: MR imaging measurements of SNR and CNR values on T1 and T2-weighted sequences and MR angiography indicated reasonable agreement with published values of AAA thrombus and abdominal components in vivo. X-ray absorption also lay within normal ranges of AAA patients and was representative of findings observed on CTA, fluoroscopy, and DSA. Ultrasound propagation speeds for developed materials were also in concordance with the literature for vascular and abdominal tissues. Conclusions: The mimicked abdominal tissues, AAA wall, and surrounding thrombus were developed to match imaging features of in vivo MR, CT, and ultrasound examinations. This phantom should be of value for image calibration, segmentation, and testing of endovascular devices for AAA endovascular repair.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Feasibility of an image planning system for kilovoltage image-guided radiation therapy

Purpose: Image guidance has become a standard of care for many treatment scenarios in radiation therapy. This is most typically accomplished by use of kV x-ray devices mounted onto the linear accelerator (Linac) gantry that yield planar, fluoroscopic, and cone-beam computed tomography (CBCT) images. Image acquisition parameters are chosen via preset techniques that rely on broad categorizations in patient anatomy and imaging goal. However, the optimal imaging technique results in detectability of the features of interest while exposing the patient to minimum dose. Herein, the authors present an investigation into the feasibility of developing an image planning system (IPS) for radiotherapy.Methods: In this first phase, the authors focused on developing an algorithm to predict tissue contrast produced by a common radiotherapy planar imaging chain. Input parameters include a CT dataset and simulated planar imaging technique settings that include kV and mAs. Energy-specific attenuation through each voxel of the CT dataset was calculated in the algorithm to derive a net transmitted intensity. The response of the flat panel detector was integrated into the image simulation algorithm. Verification was conducted by comparing simulated and measured images using four phantoms. Comparisons were made in both high and low contrast settings, as well as changes in the geometric appearance due to image saturation. Results: The authors studied a lung nodule test object to assess the planning system's ability to predict object contrast and detectability. Verification demonstrated that the slope of the pixel intensities is similar, the presence of the nodule is evident, and image saturation at high mAs values is evident in both images. The appearance of the lung nodule is a function of the image detector saturation. The authors assessed the dimensions of the lung nodule in measured and simulated images. Good quantitative agreement affirmed the algorithm's predictive capabilities. The invariance of contrast with kVp and mAs prior to saturation was predicted, as well as the gradual loss of object detectability as saturation was approached. Small changes in soft tissue density were studied using a mammography step wedge phantom. Data were acquired at beam qualities of 80 and 120 kVp and over exposure values ranging from 0.04 to 500 mAs. The data showed good agreement in terms of the absolute value of pixel intensities predicted, as well as small variations across the step wedge pattern. The saturation pixel intensity was consistent between the two beam qualities studied. Boney tissue contrast was assessed using two abdominal phantoms. Measured and calculated values agree in terms of predicting the mAs value at which detector saturation, and subsequent loss of contrast occurs. The lack of variation in contrast over mAs values lower than 10 suggests that there is wide latitude for minimizing patient dose. Conclusions: The authors developed and tested an algorithm that can be used to assist in kV imaging technique selection during localization for radiotherapy. Phantom testing demonstrated the algorithm's predictive accuracy for both low and high contrast imaging scenarios. Detector saturation with subsequent loss of imaging detail, both in terms of object size and contrast were accurately predicted by the algorithm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Transconvolution and the virtual positron emission tomograph-A new method for cross calibration in quantitative PET/CT imaging

Purpose: Positron emission tomography (PET)/computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET/CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET/CT in the context of multicenter trials. Methods: To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET/CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET/CT systems, a dedicated solid-state phantom incorporating {sup 68}Ge/{sup 68}Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET/CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. Results: The proposed Transconvolution method matched different PET/CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. Conclusions: By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Grid artifact reduction for direct digital radiography detectors based on rotated stationary grids with homomorphic filtering

Purpose: Grid artifacts are caused when using the antiscatter grid in obtaining digital x-ray images. In this paper, research on grid artifact reduction techniques is conducted especially for the direct detectors, which are based on amorphous selenium. Methods: In order to analyze and reduce the grid artifacts, the authors consider a multiplicative grid image model and propose a homomorphic filtering technique. For minimal damage due to filters, which are used to suppress the grid artifacts, rotated grids with respect to the sampling direction are employed, and min-max optimization problems for searching optimal grid frequencies and angles for given sampling frequencies are established. The authors then propose algorithms for the grid artifact reduction based on the band-stop filters as well as low-pass filters. Results: The proposed algorithms are experimentally tested for digital x-ray images, which are obtained from direct detectors with the rotated grids, and are compared with other algorithms. It is shown that the proposed algorithms can successfully reduce the grid artifacts for direct detectors. Conclusions: By employing the homomorphic filtering technique, the authors can considerably suppress the strong grid artifacts with relatively narrow-bandwidth filters compared to the normal filtering case. Using rotated grids also significantly reduces the ringing artifact. Furthermore, for specific grid frequencies and angles, the authors can use simple homomorphic low-pass filters in the spatial domain, and thus alleviate the grid artifacts with very low implementation complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MR-CT registration using a Ni-Ti prostate stent in image-guided radiotherapy of prostate cancer

Purpose: In image-guided radiotherapy of prostate cancer defining the clinical target volume often relies on magnetic resonance (MR). The task of transferring the clinical target volume from MR to standard planning computed tomography (CT) is not trivial due to prostate mobility. In this paper, an automatic local registration approach is proposed based on a newly developed removable Ni-Ti prostate stent.Methods: The registration uses the voxel similarity measure mutual information in a two-step approach where the pelvic bones are used to establish an initial registration for the local registration.Results: In a phantom study, the accuracy was measured to 0.97 mm and visual inspection showed accurate registration of all 30 data sets. The consistency of the registration was examined where translation and rotation displacements yield a rotation error of 0.41 Degree-Sign {+-} 0.45 Degree-Sign and a translation error of 1.67 {+-} 2.24 mm.Conclusions: This study demonstrated the feasibility for an automatic local MR-CT registration using the prostate stent.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A method for removing arm backscatter from EPID images

Purpose: To develop a method for removing the support arm backscatter from images acquired using current Varian electronic portal imaging devices (EPIDs).Methods: The effect of arm backscatter on EPID images was modeled using a kernel convolution method. The parameters of the model were optimized by comparing on-arm images to off-arm images. The model was used to develop a method to remove the effect of backscatter from measured EPID images. The performance of the backscatter removal method was tested by comparing backscatter corrected on-arm images to measured off-arm images for 17 rectangular fields of different sizes and locations on the imager. The method was also tested using on- and off-arm images from 42 intensity modulated radiotherapy (IMRT) fields.Results: Images generated by the backscatter removal method gave consistently better agreement with off-arm images than images without backscatter correction. For the 17 rectangular fields studied, the root mean square difference of in-plane profiles compared to off-arm profiles was reduced from 1.19% (standard deviation 0.59%) on average without backscatter removal to 0.38% (standard deviation 0.18%) when using the backscatter removal method. When comparing to the off-arm images from the 42 IMRT fields, the mean {gamma} and percentage of pixels with {gamma} < 1 were improved by the backscatter removal method in all but one of the images studied. The mean {gamma} value (1%, 1 mm) for the IMRT fields studied was reduced from 0.80 to 0.57 by using the backscatter removal method, while the mean {gamma} pass rate was increased from 72.2% to 84.6%.Conclusions: A backscatter removal method has been developed to estimate the image acquired by the EPID without any arm backscatter from an image acquired in the presence of arm backscatter. The method has been shown to produce consistently reliable results for a wide range of field sizes and jaw configurations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Real-time x-ray fluoroscopy-based catheter detection and tracking for cardiac electrophysiology interventions

Purpose: X-ray fluoroscopically guided cardiac electrophysiology (EP) procedures are commonly carried out to treat patients with arrhythmias. X-ray images have poor soft tissue contrast and, for this reason, overlay of a three-dimensional (3D) roadmap derived from preprocedural volumetric images can be used to add anatomical information. It is useful to know the position of the catheter electrodes relative to the cardiac anatomy, for example, to record ablation therapy locations during atrial fibrillation therapy. Also, the electrode positions of the coronary sinus (CS) catheter or lasso catheter can be used for road map motion correction.Methods: In this paper, the authors present a novel unified computational framework for image-based catheter detection and tracking without any user interaction. The proposed framework includes fast blob detection, shape-constrained searching and model-based detection. In addition, catheter tracking methods were designed based on the customized catheter models input from the detection method. Three real-time detection and tracking methods are derived from the computational framework to detect or track the three most common types of catheters in EP procedures: the ablation catheter, the CS catheter, and the lasso catheter. Since the proposed methods use the same blob detection method to extract key information from x-ray images, the ablation, CS, and lasso catheters can be detected and tracked simultaneously in real-time.Results: The catheter detection methods were tested on 105 different clinical fluoroscopy sequences taken from 31 clinical procedures. Two-dimensional (2D) detection errors of 0.50 {+-} 0.29, 0.92 {+-} 0.61, and 0.63 {+-} 0.45 mm as well as success rates of 99.4%, 97.2%, and 88.9% were achieved for the CS catheter, ablation catheter, and lasso catheter, respectively. With the tracking method, accuracies were increased to 0.45 {+-} 0.28, 0.64 {+-} 0.37, and 0.53 {+-} 0.38 mm and success rates increased to 100%, 99.2%, and 96.5% for the CS, ablation, and lasso catheters, respectively. Subjective clinical evaluation by three experienced electrophysiologists showed that the detection and tracking results were clinically acceptable.Conclusions: The proposed detection and tracking methods are automatic and can detect and track CS, ablation, and lasso catheters simultaneously and in real-time. The accuracy of the proposed methods is sub-mm and the methods are robust toward low-dose x-ray fluoroscopic images, which are mainly used during EP procedures to maintain low radiation dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Investigation of statistical iterative reconstruction for dedicated breast CT

Purpose: Dedicated breast CT has great potential for improving the detection and diagnosis of breast cancer. Statistical iterative reconstruction (SIR) in dedicated breast CT is a promising alternative to traditional filtered backprojection (FBP). One of the difficulties in using SIR is the presence of free parameters in the algorithm that control the appearance of the resulting image. These parameters require tuning in order to achieve high quality reconstructions. In this study, the authors investigated the penalized maximum likelihood (PML) method with two commonly used types of roughness penalty functions: hyperbolic potential and anisotropic total variation (TV) norm. Reconstructed images were compared with images obtained using standard FBP. Optimal parameters for PML with the hyperbolic prior are reported for the task of detecting microcalcifications embedded in breast tissue.Methods: Computer simulations were used to acquire projections in a half-cone beam geometry. The modeled setup describes a realistic breast CT benchtop system, with an x-ray spectra produced by a point source and an a-Si, CsI:Tl flat-panel detector. A voxelized anthropomorphic breast phantom with 280 μm microcalcification spheres embedded in it was used to model attenuation properties of the uncompressed woman's breast in a pendant position. The reconstruction of 3D images was performed using the separable paraboloidal surrogates algorithm with ordered subsets. Task performance was assessed with the ideal observer detectability index to determine optimal PML parameters.Results: The authors' findings suggest that there is a preferred range of values of the roughness penalty weight and the edge preservation threshold in the penalized objective function with the hyperbolic potential, which resulted in low noise images with high contrast microcalcifications preserved. In terms of numerical observer detectability index, the PML method with optimal parameters yielded substantially improved performance (by a factor of greater than 10) compared to FBP. The hyperbolic prior was also observed to be superior to the TV norm. A few of the best-performing parameter pairs for the PML method also demonstrated superior performance for various radiation doses. In fact, using PML with certain parameter values results in better images, acquired using 2 mGy dose, than FBP-reconstructed images acquired using 6 mGy dose.Conclusions: A range of optimal free parameters for the PML algorithm with hyperbolic and TV norm-based potentials is presented for the microcalcification detection task, in dedicated breast CT. The reported values can be used as starting values of the free parameters, when SIR techniques are used for image reconstruction. Significant improvement in image quality can be achieved by using PML with optimal combination of parameters, as compared to FBP. Importantly, these results suggest improved detection of microcalcifications can be obtained by using PML with lower radiation dose to the patient, than using FBP with higher dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Implementation and commissioning of an integrated micro-CT/RT system with computerized independent jaw collimation

Purpose: To design, construct, and commission a set of computer-controlled motorized jaws for a micro-CT/RT system to perform conformal image-guided small animal radiotherapy.Methods: The authors designed and evaluated a system of custom-built motorized orthogonal jaws, which allows the delivery of off-axis rectangular fields on a GE eXplore CT 120 preclinical imaging system. The jaws in the x direction are independently driven, while the y-direction jaws are symmetric. All motors have backup encoders, verifying jaw positions. Mechanical performance of the jaws was characterized. Square beam profiles ranging from 2 × 2 to 60 × 60 mm{sup 2} were measured using EBT2 film in the center of a 70 × 70 × 22 mm{sup 3} solid water block. Similarly, absolute depth dose was measured in a solid water and EBT2 film stack 50 × 50 × 50 mm{sup 3}. A calibrated Farmer ion chamber in a 70 × 70 × 20 mm{sup 3} solid water block was used to measure the output of three field sizes: 50 × 50, 40 × 40, and 30 × 30 mm{sup 2}. Elliptical target plans were delivered to films to assess overall system performance. Respiratory-gated treatment was implemented on the system and initially proved using a simple sinusoidal motion phantom. All films were scanned on a flatbed scanner (Epson 1000XL) and converted to dose using a fitted calibration curve. A Monte Carlo beam model of the micro-CT with the jaws has been created using BEAMnrc for comparison with the measurements. An example image-guided partial lung irradiation in a rat is demonstrated.Results: The averaged random error of positioning each jaw is less than 0.1 mm. Relative output factors measured with the ion chamber agree with Monte Carlo simulations within 2%. Beam profiles and absolute depth dose curves measured from the films agree with simulations within measurement uncertainty. Respiratory-gated treatments applied to a phantom moving with a peak-to-peak amplitude of 5 mm showed improved beam penumbra (80%–20%) from 3.9 to 0.8 mm.Conclusions: A set of computer-controlled motorized jaws for a micro-CT/RT system were constructed with position reliably better than a tenth of a millimeter. The hardware system is ready for image-guided conformal radiotherapy for small animals with capability of respiratory-gated delivery.

62 RADIOLOGY AND NUCLEAR MEDICINE↗