Threshold segmentation for PET target volume delineation in radiation treatment planning: the role of target-to-background ratio and target size

A multivariable approach was adopted to study the dependence of the percentage threshold [TH(%)] used to define the boundaries of 18F-FDG positive tissue on emission scan duration (ESD) and activity at the start of acquisition (Aacq) for different target sizes and target-to-background (T/B) ratios. An anthropomorphic model, at least for counting rate characteristics, was used to study this dependence in conditions resembling the ones that can be encountered in the clinical studies. An annular ring of water bags of 3 cm thickness was fitted over an International Electro-technical Commission (IEC) phantom in order to obtain counting rates similar to those found in average patients. The scatter fraction of the modified IEC phantom was similar to the mean scatter fraction measured on patients, with a similar scanner. A supplemental set of microhollow spheres was positioned inside the phantom. The NEMA NU 2-2001 scatter phantom was positioned at the end of the IEC phantom to approximate the clinical situation of having activity that extends beyond the scanner field of view. The phantoms were filled with a solution of water and 18F (12 kBq/mL) and the spheres with various T/B ratios of 22.5, 10.3, and 3.6. Sequential imaging was performed to acquire PET images with varying background activity concentrations of about 12, 9, 6.4, 5.3, and 3.1 kBq/mL. The ESD was set to 60, 120, 180, and 240 s/bed. Data were fitted using two distinct multiple linear regression models for sphere ID < or = 10 mm and sphere ID > 10 mm. The fittings of both models were good with an R2 of 0.86 in both cases. Neither ESD nor Aacq resulted as significant predictors of the TH(%). For sphere ID < or =10 mm the target size was the most significant predictor of the TH(%), followed by the T/B ratio, while for sphere ID > 10 mm the explanatory power of the target size and T/B ratio were reversed, the T/B ratio being now the most important predictor of the TH(%). Both the target size and T/B ratio play a major role in explaining the variance of the TH(%), throughout the whole range of target sizes and T/B ratios examined. Thus, algorithms aimed at automatic threshold segmentation should incorporate both variables with a relative weight which critically depends on target size.     

Design and performance evaluation of a whole-body Ingenuity TF PET-MRI system

The Ingenuity TF PET-MRI is a newly released whole-body hybrid PET-MR imaging system with a Philips time-of-flight GEMINI TF PET and Achieva 3T X-series MRI system. Compared to PET-CT, modifications to the positron emission tomography (PET) gantry were made to avoid mutual system interference and deliver uncompromising performance which is equivalent to the standalone systems. The PET gantry was redesigned to introduce magnetic shielding for the photomultiplier tubes (PMTs). Stringent electromagnetic noise requirements of the MR system necessitated the removal of PET gantry electronics to be housed in the PET-MR equipment room. We report the standard NEMA measurements for the PET scanner. PET imaging and performance measurements were done at Geneva University Hospital as described in the NEMA Standards NU 2-2007 manual. The scatter fraction (SF) and noise equivalent count rate (NECR) measurements with the NEMA cylinder (20 cm diameter) were repeated for two larger cylinders (27 cm and 35 cm diameter), which better represent average and heavy patients. A NEMA/IEC torso phantom was used for overall assessment of image quality. The transverse and axial resolution near the center was 4.7 mm. Timing and energy resolution of the PET-MR system were measured to be 525 ps and 12%, respectively. The results were comparable to PET-CT systems demonstrating that the effect of design modifications required on the PET system to remove the harmful effect of the magnetic field on the PMTs was negligible. The absolute sensitivity of this scanner was 7.0 cps kBq(-1), whereas SF was 26%. NECR measurements performed with cylinders having three different diameters, and image quality measurements performed with IEC phantom yielded excellent results. The Ingenuity TF PET-MRI represents the first commercial whole-body hybrid PET-MRI system. The performance of the PET subsystem was comparable to the GEMINI TF PET-CT system using phantom and patient studies. It is conceived that advantages of hybrid PET-MRI will become more evident in the near future.     

Characterization of a single LSO crystal layer high resolution research tomograph

The purpose of this study was to determine the performance of a single lutetium oxy-orthosilicate (LSO) crystal layer High Resolution Research Tomograph (HRRT) positron emission tomography (PET) scanner. The HRRT is a high resolution PET scanner designed for human brain and small animal imaging. The scanner consists of eight panel detectors, which have one layer of 2.1 x 2.1 x 7.5 mm thick LSO crystals. Several phantom studies were performed to determine scanner characteristics, such as resolution, scatter fraction, count rate and noise equivalent count rates (NECR). NECR curves were measured according to both NEMA NU2-1994 and NU2-2001 for three different energy windows, i.e. lower level discriminators (lld) of 350, 400 and 450 keV and an upper level discriminator (uld) of 650 keV. Accuracy of scatter and single photon attenuation corrections was evaluated according to NU2-1994. Data were acquired using a ring difference of 67 and a span of 9. Reconstructions were performed using FORE + 2D FBP or OSEM. Transaxial resolution varied from 2.7 to 2.9 mm FWHM between I and 10 cm off centre locations, and axial resolution varied from 3.2 to 4.4 mm FWHM. Scatter fractions (NU2-1994) equalled 0.31, 0.42 and 0.54 for lld of 450, 400 and 350 keV, respectively. NECR data were highest for an lid of 400 keV and showed a maximum of 46 kcps at 38 kBq cm(-3). Lower NECR values were observed according to NU2-2001, but were still optimal for an lld of 400 keV. After scatter and attenuation corrections, pixel values within water, air and teflon inserts of the NU2-1994 phantom were 14, 4 and 35% of the background activity, respectively. The single layer LSO HRRT scanner shows excellent spatial resolution, making it suitable for small animal studies. The low count rate performance, due to the small amount of LSO, prohibits studies of the human brain, but is sufficient for studies in small laboratory animals.     

Impact of Ge-68/Ga-68-based versus CT-based attenuation correction on PET

Transmission (Tx) scans are used in PET for attenuation correction (AC). For standalone PET this is typically done using Ge-68/Ga-68 sources, for PET-CT using CT. Therefore, standalone PET suffers from emission contamination during Tx scans, PET-CT does not. Here, we studied the effects of AC across the two systems. With a cylindrical phantom (Jaszczak Phantom, Data Spectrum Corp.) with hollow spheres (diameter 10-60 mm) two studies were performed. In the first study the hollow spheres were filled with 150 kBq/ml FDG and the background with 15 kBq/ml. In the second study we used 120 kBq/ml in the spheres and 50 kBq/ml in the background. Both a low and a high object-to-background ratio are studied this way. Multiple scans were acquired on a standalone PET and a PET-CT until 1% of the initial concentration remained. Activity concentration in the spheres and background was measured from the reconstructed images and compared to the actual concentration. For standalone PET, emission scans were reconstructed using hot Tx (emission contaminated) and cold Tx (not contaminated). Uniformity within the spheres was investigated by profile analysis. For PET-CT, the concentration in the big spheres (> 16 mm) was recovered. For the smaller spheres, recovery was insufficient due to partial volume effects. For standalone PET the recoveries of the spheres (> 16 mm) were 20% (first study) and 13% (second study) lower than the actual concentration. Using hot Tx, underestimation of activity concentration was up to > 50%. Nonuniformities within the biggest spheres were up to 35%, 12%, and 5% (first study), using standalone PET with hot Tx, cold Tx, and using PET-CT, respectively. Due to contamination of AC by emission photons, standalone PET results in a bias in the activity concentration and uniformity. Especially when patients get follow-up PET scans on both standalone PET and PET-CT, this may lead to misinterpretation.     

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

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

Effects of injected dose, BMI and scanner type on NECR and image noise in PET imaging

Noise equivalent count rate (NECR) and image noise are two different but related metrics that have been used to predict and assess image quality, respectively. The aim of this study is to investigate, using patient studies, the relationships between injected dose (ID), body mass index (BMI) and scanner type on NECR and image noise measurements in PET imaging. Two groups of 90 patients each were imaged on a GE DSTE and a DRX PET/CT scanner, respectively. The patients in each group were divided into nine subgroups according to three BMI (20-24.9, 25-29.9, 30-45 kg m(-2)) and three ID (296-444, 444-555, 555-740 MBq) ranges, resulting in ten patients/subgroup. All PET data were acquired in 3D mode and reconstructed using the VuePoint HD? fully 3D OSEM algorithm (2 iterations, 21(DRX) or 20 (DSTE) subsets). NECR and image noise measurements for bed positions covering the liver were calculated for each patient. NECR was calculated from the trues, randoms and scatter events recorded in the DICOM header of each patient study, while image noise was determined as the standard deviation of 50 non-neighboring voxels in the liver of each patient. A t-test compared the NECR and image noise for different scanners but with the same BMI and ID. An ANOVA test on the other hand was used to compare the results of patients with different BMI but the same ID and scanner type as well as different ID but the same BMI and scanner type. As expected the t-test showed a significant difference in NECR between the two scanners for all BMI and ID subgroups. However, contrary to what is expected no such findings were observed for image noise measurement. The ANOVA results showed a statistically significant difference in both NECR and image noise among the different BMI for each ID and scanner subgroup. However, there was no statistically significant difference in NECR and image noise across different ID for each BMI and scanner subgroup. Although the GE DRX PET/CT scanner has better count rate performance than the GE DSTE PET/CT scanner, this improvement does not translate to a lower image noise when using OSEM reconstruction. Our results show that patients with larger BMI consistently generate poorer image quality. Dose reduction from >555 to 296-444 MBq has minimal impact on image quality independent of the scanner used. A reduction in ID decreases patient and technologist exposure and can potentially reduce the overall cost of the study.     

NEMA NU-04-based performance characteristics of the LabPET-8™ small animal PET scanner

The objective of this study is to characterize the performance of the preclinical avalanche photodiode (APD)-based LabPET-8? subsystem of the fully integrated trimodality PET/SPECT/CT Triumph? scanner using the National Electrical Manufacturers Association (NEMA) NU 04-2008 protocol. The characterized performance parameters include the spatial resolution, sensitivity, scatter fraction, counts rate performance and image-quality characteristics. The PET system is fully digital using APD-based detector modules with highly integrated electronics. The detector assembly consists of phoswich pairs of Lu(1.9)Y(0.1)SiO(5) (LYSO) and Lu(0.4)Gd(1.6)SiO(5) (LGSO) crystals with dimensions of 2 × 2 × 14 mm(3) having 7.5 cm axial and 10 cm transverse field of view (FOV). The spatial resolution and sensitivity were measured using a small (22)Na point source at different positions in the scanner's FOV. The scatter fraction and count rate characteristics were measured using mouse- and rat-sized phantoms fitted with an (18)F line source. The overall imaging capabilities of the scanner were assessed using the NEMA image-quality phantom and laboratory animal studies. The NEMA-based radial and tangential spatial resolution ranged from 1.7 mm at the center of the FOV to 2.59 mm at a radial offset of 2.5 cm and from 1.85 mm at the center of the FOV to 1.76 mm at a radial offset of 2.5 cm, respectively. Iterative reconstruction improved the spatial resolution to 0.84 mm at the center of the FOV. The total absolute system sensitivity is 12.74% for an energy window of 250-650 keV. For the mouse-sized phantom, the peak noise equivalent count rate (NECR) is 183 kcps at 2.07 MBq cc(-1), whereas the peak true count rate is 320 kcps at 2.5 MBq cc(-1) with a scatter fraction of 19%. The rat-sized phantom had a scatter fraction of 31%, with a peak NECR of 67 kcps at 0.23 MBq cc(-1) and a peak true count rate of 186 kcps at 0.27 MBq cc(-1). The average activity concentration and percentage standard deviation were 126.97 kBq ml(-1) and 7%, respectively. The performance of the LabPET-8? scanner was characterized based on the NEMA NU 04-2008 standards. The all in all performance demonstrates that the LabPET-8? system is able to produce high-quality and highly contrasted images in a reasonable time, and as such it is well suited for preclinical molecular imaging-based research.     

Characterization, prediction, and correction of geometric distortion in 3 T MR images

The work presented herein describes our methods and results for predicting, measuring and correcting geometric distortions in a 3 T clinical magnetic resonance (MR) scanner for the purpose of image guidance in radiation treatment planning. Geometric inaccuracies due to both inhomogeneities in the background field and nonlinearities in the applied gradients were easily visualized on the MR images of a regularly structured three-dimensional (3D) grid phantom. From a computed tomography scan, the locations of just under 10 000 control points within the phantom were accurately determined in three dimensions using a MATLAB-based computer program. MR distortion was then determined by measuring the corresponding locations of the control points when the phantom was imaged using the MR scanner. Using a reversed gradient method, distortions due to gradient nonlinearities were separated from distortions due to inhomogeneities in the background B0 field. Because the various sources of machine-related distortions can be individually characterized, distortions present in other imaging sequences (for which 3D distortion cannot accurately be measured using phantom methods) can be predicted negating the need for individual distortion calculation for a variety of other imaging sequences. Distortions were found to be primarily caused by gradient nonlinearities and maximum image distortions were reported to be less than those previously found by other researchers at 1.5 T. Finally, the image slices were corrected for distortion in order to provide geometrically accurate phantom images.     

Effect of motion on tracer activity determination in CT attenuation corrected PET images: a lung phantom study

Respiratory motion is known to affect the quantitation of 18FDG uptake in lung lesions. The aim of the study was to investigate the magnitude of errors in tracer activity determination due to motion, and its dependence upon CT attenuation at different phases of the motion cycle. To estimate these errors we have compared maximum activity concentrations determined from PET/CT images of a lung phantom at rest and under simulated respiratory motion. The NEMA 2001 IEC body phantom, containing six hollow spheres with diameters 37, 28, 22, 17, 13, and 10 mm, was used in this study. To mimic lung tissue density, the phantom (excluding spheres) was filled with low density polystyrene beads and water. The phantom spheres were filled with 18FDG solution setting the target-to-background activity concentration ratio at 8:1. PET/CT data were acquired with the phantom at rest, and while it was undergoing periodic motion along the longitudinal axis of the scanner with a range of displacement being 2 cm, and a period of 5 s. The phantom at rest and in motion was scanned using manufacturer provided standard helical/clinical protocol, a helical CT scan followed by a PET emission scan. The moving phantom was also scanned using a 4D-CT protocol that provides volume image sets at different phases of the motion cycle. To estimate the effect of motion on quantitation of activities in six spheres, we have examined the activity concentration data for (a) the stationary phantom, (b) the phantom undergoing simulated respiratory motion, and (c) a moving phantom acquired with PET/4D-CT protocol in which attenuation correction was performed with CT images acquired at different phases of motion cycle. The data for the phantom at rest and in motion acquired with the standard helical/clinical protocol showed that the activity concentration in the spheres can be underestimated by as much as 75%, depending on the sphere diameter. We have also demonstrated that fluctuations in sphere's activity concentration from one PET/CT scan to another acquired with standard helical/clinical protocol can arise as a consequence of spatial mismatch between the sphere's location in PET emission and the CT data.     

ROC analysis for assessment of lesion detection performance in 3D PET: influence of reconstruction algorithms

Image quality in positron emission tomography (PET) can be assessed with physical parameters, as spatial resolution and signal-to-noise ratio, or using psychophysical approaches, which include the observer performance and the considered task (ROC analysis). For PET in oncology, such a task is the detection of hot lesions. The aim of the present study was to assess the lesion detection performance due to adequate modeling of the scanner and the measurement process in the image reconstruction process. We compared the standard OSEM software of the manufacturer with a sophisticated fully 3D iterative reconstruction technique (USC MAP). A rectangular phantom with 6 oblique line sources in a homogeneous background (2.6 kBq/ml 18F) was imaged dynamically with an ECAT EXACT HR+ scanner in 3D mode. Reconstructed activity contrasts varied between 15 and 0, as the line sources were filled with 11C (3.2 MBq/ml). Measured attenuation and standard randoms, dead time, and scatter corrections of the manufacturer were employed. For the ROC analysis, a software tool presented a cut-out of the phantom (15 x 15 pixels) to two observers. These cut-outs were rated (5 classes) and the area Az under the ROC curve was determined as a measure of detection performance. The improvement for Az with USC MAP compared to the OSEM reconstructions ranged between 0.02 and 0.23 for signal-to-noise ratios of the background between 2.8 and 3.1 and lesion contrast between 2.1 and 4.2. This study demonstrates that adequate modeling of the measurement process in the reconstruction algorithm improves the detection of small hot lesions markedly.     

Use of an in-field-of-view shield to improve count rate performance of the single crystal layer high-resolution research tomograph PET scanner for small animal brain scans

The count rate performance of the single LSO crystal layer high-resolution research tomograph (HRRT-S) PET scanner is limited by the processing speed of its electronics. Therefore, the feasibility of using an in-field-of-view (in-FOV) shield to improve the noise equivalent count rates (NECR) for small animal brain studies was investigated. The in-FOV shield consists of a lead tube of 12 cm length, 6 cm inner diameter and 9 mm wall thickness. It is large enough to shield the activity in the body of a rat or mouse. First, the effect of this shield on NECR was studied. Secondly, a number of experiments were performed to assess the effects of the shield on the accuracy of transmission scan data and, next, on reconstructed activity distribution in the brain. For activities below 150 MBq NECR improved only by 5-10%. For higher activities NECR maxima of 1.2E4 cps at 200 MBq and 2.2E4 cps at 370 MBq were found without and with shield, respectively. Listmode data taken without shield, however, were corrupted for activities above 75 MBq due to data overrun problems (time tag losses) of the electronics. When the shield was used data overrun was avoided up to activities of 150 MBq. For the unshielded part of the phantom, transmission scan data were the same with and without shield. The estimated scatter contribution was approximately 8.5% without and 5.5% with shield. Reconstructed emission data showed a difference up to 5% in the unshielded part of the phantom at 5 mm or more from the edge of the shielding. Of this 5% about 3% results from the difference in the uncorrected scatter contribution. In conclusion, an in-FOV shield can be used successfully in an HRRT PET scanner to improve NECR and accuracy of small animal brain studies. The latter is especially important when high activities are required for tracers with low brain uptake or when multiple animals are scanned simultaneously.     

Determining the minimum number of detectable cardiac-transplanted 111In-tropolone-labelled bone-marrow-derived mesenchymal stem cells by SPECT

In this work, we determined the minimum number of detectable 111In-tropolone-labelled bone-marrow-derived stem cells from the maximum activity per cell which did not affect viability, proliferation and differentiation, and the minimum detectable activity (MDA) of 111In by SPECT. Canine bone marrow mesenchymal cells were isolated, cultured and expanded. A number of samples, each containing 5x10(6) cells, were labelled with 111In-tropolone from 0.1 to 18 MBq, and cell viability was measured afterwards for each sample for 2 weeks. To determine the MDA, the anthropomorphic torso phantom (DataSpectrum Corporation, Hillsborough, NC) was used. A point source of 202 kBq 111In was placed on the surface of the heart compartment, and the phantom and all compartments were then filled with water. Three 111In SPECT scans (duration: 16, 32 and 64 min; parameters: 128x128 matrix with 128 projections over 360 degrees) were acquired every three days until the 111In radioactivity decayed to undetectable quantities. 111In SPECT images were reconstructed using OSEM with and without background, scatter or attenuation corrections. Contrast-to-noise ratio (CNR) in the reconstructed image was calculated, and MDA was set equal to the 111In activity corresponding to a CNR of 4. The cells had 100% viability when incubated with no more than 0.9 MBq of 111In (80% labelling efficiency), which corresponded to 0.14 Bq per cell. Background correction improved the detection limits for 111In-tropolone-labelled cells. The MDAs for 16, 32 and 64 min scans with background correction were observed to be 1.4 kBq, 700 Bq and 400 Bq, which implies that, in the case where the location of the transplantation is known and fixed, as few as 10,000, 5000 and 2900 cells respectively can be detected.     

Physical and clinical performance of the mCT time-of-flight PET/CT scanner

Time-of-flight (TOF) measurement capability promises to improve PET image quality. We characterized the physical and clinical PET performance of the first Biograph mCT TOF PET/CT scanner (Siemens Medical Solutions USA, Inc.) in comparison with its predecessor, the Biograph TruePoint TrueV. In particular, we defined the improvements with TOF. The physical performance was evaluated according to the National Electrical Manufacturers Association (NEMA) NU 2-2007 standard with additional measurements to specifically address the TOF capability. Patient data were analyzed to obtain the clinical performance of the scanner. As expected for the same size crystal detectors, a similar spatial resolution was measured on the mCT as on the TruePoint TrueV. The mCT demonstrated modestly higher sensitivity (increase by 19.7 ± 2.8%) and peak noise equivalent count rate (NECR) (increase by 15.5 ± 5.7%) with similar scatter fractions. The energy, time and spatial resolutions for a varying single count rate of up to 55 Mcps resulted in 11.5 ± 0.2% (FWHM), 527.5 ± 4.9 ps (FWHM) and 4.1 ± 0.0 mm (FWHM), respectively. With the addition of TOF, the mCT also produced substantially higher image contrast recovery and signal-to-noise ratios in a clinically-relevant phantom geometry. The benefits of TOF were clearly demonstrated in representative patient images.     

Optimal whole-body PET scanner configurations for different volumes of LSO scintillator: a simulation study

The axial field of view (AFOV) of the current generation of clinical whole-body PET scanners range from 15-22?cm, which limits sensitivity and renders applications such as whole-body dynamic imaging or imaging of very low activities in whole-body cellular tracking studies, almost impossible. Generally, extending the AFOV significantly increases the sensitivity and count-rate performance. However, extending the AFOV while maintaining detector thickness has significant cost implications. In addition, random coincidences, detector dead time, and object attenuation may reduce scanner performance as the AFOV increases. In this paper, we use Monte Carlo simulations to find the optimal scanner geometry (i.e. AFOV, detector thickness and acceptance angle) based on count-rate performance for a range of scintillator volumes ranging from 10 to 93 l with detector thickness varying from 5 to 20?mm. We compare the results to the performance of a scanner based on the current Siemens Biograph mCT geometry and electronics. Our simulation models were developed based on individual components of the Siemens Biograph mCT and were validated against experimental data using the NEMA NU-2 2007 count-rate protocol. In the study, noise-equivalent count rate (NECR) was computed as a function of maximum ring difference (i.e. acceptance angle) and activity concentration using a 27?cm diameter, 200?cm uniformly filled cylindrical phantom for each scanner configuration. To reduce the effect of random coincidences, we implemented a variable coincidence time window based on the length of the lines of response, which increased NECR performance up to 10% compared to using a static coincidence time window for scanners with a large maximum ring difference values. For a given scintillator volume, the optimal configuration results in modest count-rate performance gains of up to 16% compared to the shortest AFOV scanner with the thickest detectors. However, the longest AFOV of approximately 2?m with 20?mm thick detectors resulted in performance gains of 25-31?times higher NECR relative to the current Siemens Biograph mCT scanner configuration.     

Effect of area x-ray beam equalization on image quality and dose in digital mammography

In mammography, thick or dense breast regions persistently suffer from reduced contrast-to-noise ratio (CNR) because of degraded contrast from large scatter intensities and relatively high noise. Area x-ray beam equalization can improve image quality by increasing the x-ray exposure to under-penetrated regions without increasing the exposure to other breast regions. Optimal equalization parameters with respect to image quality and patient dose were determined through computer simulations and validated with experimental observations on a step phantom and an anthropomorphic breast phantom. Three parameters important in equalization digital mammography were considered: attenuator material (Z = 13-92), beam energy (22-34 kVp) and equalization level. A Mo/Mo digital mammography system was used for image acquisition. A prototype 16 x 16 piston driven equalization system was used for preparing patient-specific equalization masks. Simulation studies showed that a molybdenum attenuator and an equalization level of 20 were optimal for improving contrast, CNR and figure of merit (FOM = CNR2/dose). Experimental measurements using these parameters showed significant improvements in contrast, CNR and FOM. Moreover, equalized images of a breast phantom showed improved image quality. These results indicate that area beam equalization can improve image quality in digital mammography.     

Quantification of lesion size, depth, and uptake using a dual-head molecular breast imaging system

A method to perform quantitative lesion analysis in molecular breast imaging (MBI) was developed using the opposing views from a novel dual-head dedicated gamma camera. Monte Carlo simulations and phantom models were used to simulate MBI images with known lesion parameters. A relationship between the full widths at 25%, 35%, and 50% of the maximum of intensity profiles through lesions and the true lesion diameter as a function of compressed breast thickness was developed in order to measure lesion diameter. Using knowledge of compressed breast thickness and the attenuation of gamma rays in soft tissue, a method was developed to measure the depth of the lesion to the collimator face. Using the measured lesion diameter and measurements of counts in the lesion and background breast region, relative radiotracer uptake or tumor to background ratio (T/B ratio) was calculated. Validation of the methods showed that the size, depth, and T/B ratio can be accurately measured for a range of small breast lesions with T/B ratios between 10:1 and 40:1 in breasts with compressed thicknesses between 4 and 10 cm. Future applications of this work include providing information about lesion location in patients for performing a biopsy of site and the development of a threshold for the T/B ratio that can distinguish benign from malignant disease.     

Performance of a novel collimator for high-sensitivity brain SPECT

We assessed improvements in performance in detection and estimation tasks due to a novel brain single photon computed tomography collimator. Data were acquired on the CeraSPECT scanner using both new and standard collimators. The new variable focusing collimator SensOgrade samples the projections unequally, with central regions more heavily represented, to compensate for attenuation of counts from central brain structures. Furthermore, it utilizes more of the cylindrical crystal surface. Two phantom studies were performed. The first phantom was a 21-cm-diameter cylindrical background containing nine spheres ranging from 0.5 to 5 cm3 in volume. 99mTc sphere to background activity ratio was 10:1. Twenty-nine 10-min datasets were acquired with each collimator. The second phantom was the Radiology Support Devices (Long Beach, CA) striatal phantom with striatal-background ratios of 10:1 on the left and 5:1 on the right. Twenty-nine 4-min datasets were acquired with each collimator. Perfusion imaging using 99mTc-HMPAO was also performed in three healthy volunteers using both collimators under identical simulations. Projections were reconstructed by filtered backprojection with an unwindowed ramp filter. The nonprewhitening matched filter signal-to-noise ratio (NPW-SNR) was computed as a surrogate for human performance in detecting spherical lesions. Sphere activity concentration, radius, and location coordinates were simultaneously estimated by fitting images to an assumed model using an iterative nonlinear algorithm. Resolution recovery was implicit in the estimation procedure, as the point spread function was incorporated into the model. NPW-SNR for sphere detection was 1.5 to 2 times greater with the new collimator; for the striatal phantom the improvement in SNR was 54%. The SNR for estimating sphere activity concentration improved by 46 to 89% for spheres located more than 5 cm from the phantom center. Images acquired with the standard collimator were too noisy in the central regions to allow estimation of sphere activity. In 99mTc-HMPAO human studies, SNR was improved by 21 to 41% in the cortex, 66% in the basal ganglia, and 74% in the thalamus. The new collimator leads to substantially improved detection and estimation performance throughout the brain. The higher sensitivity will be particularly important for dynamic imaging.     

The potential role of positron emission mammography for detection of breast cancer. A phantom study

Positron emission mammography (PEM) is a new, specialized imaging modality utilizing PET radiopharmaceuticals to detect breast cancer. The capabilities and limitations of PEM in detecting breast tumors were investigated with a series of phantom experiments. The PEM imager was mounted on a standard Lorad biopsy table (separated by 18 cm). In the initial phase of the investigation, basic scanner parameters (resolution, sensitivity, and scatter fraction) were measured. The effects of a number of breast imaging parameters (length of acquisition, breast thickness, and breast density) on detection of breast lesions were then explored utilizing special phantoms. Moderately compressed breasts were simulated with a block of gelatin containing amounts of FDG consistent with 370 MBq injections. Lesions were simulated with four hollow spheres (inner diameters=5 mm, 8 mm, 12 mm, and 15 mm) filled with amounts of FDG representative of uptake in malignant breast tumors (target-to-background concentration ratio=8.5:1). Resolution at the center of the imager was 3.9 mm, sensitivity was 0.059 kcps/kBq/ml and the Compton scatter fraction was approximately 12%. Objects as small as 8 mm in diameter could be detected after 30 s of data acquisition; 5 mm spheres were detectable after 300 s. Object detection capabilities were reduced with increasing breast thickness. In thin compressed breasts (2 cm) even the smallest sphere (5 mm in diameter) could be detected; increasing breast thickness increased the minimum detectable sphere diameter to 8 mm. Increased background activity caused by FDG uptake in metabolically active normal tissue more prevalent in radiodense breasts compared to "fatty" breasts was simulated and shown to reduce the minimum detectable lesion size to 12 mm for the densest breasts. These results demonstrate the potential of PEM for the detection of breast lesions. The addition of the system to a standard biopsy apparatus indicates its potential for use to guide some core biopsies of breast cancers.     

Effect of 176Lu background on singles transmission for LSO-based PET cameras

We explore how the radioactive background from naturally occurring 176Lu affects single photon transmission imaging for lutetium orthosilicate (LSO) scintillator-based PET cameras by estimating the transmission noise equivalent count rate (NECR) including this background. Assuming a typical PET camera geometry (80 cm detector ring diameter), we use a combination of measurement and analytic computation to estimate the counting rates due to transmission, scatter and background events as a function of singles transmission source strength. We then compute a NECR for singles transmission. We find that the presence of radiation from the naturally occurring 176Lu reduces the NECR by 60% or higher for source strengths less than 10 mCi, and that a 25% reduction of the NECR can occur even with a source strength of 40 mCi.     

An investigation of entrance surface dose calculations for diagnostic radiology using Monte Carlo simulations and radiotherapy dosimetry formalisms

Our aim in this work was to investigate the methodology used in the determination of the entrance surface dose (ESD) in diagnostic radiology. In kV x-rays for low-energy photons (tube potential up to 160 kV, HVL: 1-8 mm Al), the ESD is based on the use of the ratio of mass-energy absorption coefficients and backscatter factors. A full simulation of the photon and electron transport in a kilovoltage x-ray unit, using the Monte Carlo code BEAM/EGS4, was performed to obtain an accurate beam phase space for use in dose calculation. The modelled phase space was experimentally validated for the beam qualities (measured HVL: 3.3 mm Al-2.2 mm Cu) and showed good agreement between calculated and measured HVLs, air kerma and relative dose distributions. We have computed the conversion factors from air kerma to water or soft tissue absorbed dose at the surface of a phantom for beam qualities (HVL: 3.3-8.35 mm Al). The same model was also used to calculate the ESD in water and in soft tissue for the low-energy photon range considered. The results show that the numerical differences between the air kerma and the water kerma based backscatter factors are insignificant. The same conclusion was reached for the (mu(en)/rho) ratios, for soft tissue to air, evaluated using either the primary photon spectra or the spectra at the surface of a phantom. Furthermore, the good agreement obtained for the computation of the conversion factors with a full BEAM/EGS4 model confirms the previous studies which are based on different sources for the spectral distribution and different beam geometries (pencil beam or point source assumptions). On the other hand, the ESD in water or soft tissue is well described either with the B(air) or the B(w) formalism. Conversion factors from air kerma to ESD in these media are proposed in this work for several beam qualities in diagnostic radiology.