Optimally Stabilized PET Image Denoising Using Trilateral Filtering

Low-resolution and signal-dependent noise distribution in positron emission tomography (PET) images makes denoising process an inevitable step prior to qualitative and quantitative image analysis tasks. Conventional PET denoising methods either over-smooth small-sized structures due to resolution limitation or make incorrect assumptions about the noise characteristics. Therefore, clinically important quantitative information may be corrupted. To address these challenges, we introduced a novel approach to remove signal-dependent noise in the PET images where the noise distribution was considered as Poisson-Gaussian mixed. Meanwhile, the generalized Anscombe's transformation (GAT) was used to stabilize varying nature of the PET noise. Other than noise stabilization, it is also desirable for the noise removal filter to preserve the boundaries of the structures while smoothing the noisy regions. Indeed, it is important to avoid significant loss of quantitative information such as standard uptake value (SUV)-based metrics as well as metabolic lesion volume. To satisfy all these properties, we extended bilateral filtering method into trilateral filtering through multiscaling and optimal Gaussianization process. The proposed method was tested on more than 50 PET-CT images from various patients having different cancers and achieved the superior performance compared to the widely used denoising techniques in the literature.Comment: 8 pages, 3 figures; to appear in the Lecture Notes in Computer Science (MICCAI 2014

Evaluation of resistive-plate-chamber-based TOF-PET applied to in-beam particle therapy monitoring

Particle therapy is a highly conformal radiotherapy technique which reduces the dose deposited to the surrounding normal tissues. In order to fully exploit its advantages, treatment monitoring is necessary to minimize uncertainties related to the dose delivery. Up to now, the only clinically feasible technique for the monitoring of therapeutic irradiation with particle beams is Positron Emission Tomography (PET). In this work we have compared a Resistive Plate Chamber (RPC)-based PET scanner with a scintillation-crystal-based PET scanner for this application. In general, the main advantages of the RPC-PET system are its excellent timing resolution, low cost, and the possibility of building large area systems. We simulated a partial-ring scannerbeam monitoring, which has an intrinsically low positron yield compared to diagnostic PET. In addition, for in-beam PET there is a further data loss due to the partial ring configuration. In order to improve the performance of the RPC-based scanner, an improved version of the RPC detector (modifying the thickness of the gas and glass layers), providing a larger sensitivity, has been simulated and compared with an axially extended version of the crystal-based device. The improved version of the RPC shows better performance than the prototype, but the extended version of the crystal-based PET outperforms all other options. based on an RPC prototype under construction within the Fondazione per Adroterapia Oncologica (TERA). For comparison with the crystal-based PET scanner we have chosen the geometry of a commercially available PET scanner, the Philips Gemini TF. The coincidence time resolution used in the simulations takes into account the current achievable values as well as expected improvements of both technologies. Several scenarios (including patient data) have been simulated to evaluate the performance of different scanners. Initial results have shown that the low sensitivity of the RPC hampers its application to hadro

Optimization of PET protocol and interrater reliability of 18F-PSMA-11 imaging of prostate cancer

Background Several scan parameters for PET imaging with F-18-PSMA-11 such as dosage, acquisition time and scan duration were evaluated to determine the most appropriate scan protocol, as well as the effect of furosemide administration on lesion visualization. Forty-four patients were randomly assigned to a dosage group (2.0 +/- 0.2 or 4.0 +/- 0.4 MBq/kg F-18-PSMA-11). All patients received a full-body PET/CT 1 h and 3 h after radiotracer injection with a scan duration of 3 min/bed position. For comparison of the scan duration, images were reconstructed for 1.5 and 3 min/bed position. Patients were intravenously administered 0.5 mg/kg furosemide with a maximum dose of 40 mg. To evaluate the furosemide effect, 22 additional patients were recruited and received one full-body PET/CT 1 h after administration of 2.0 +/- 0.2 MBq/kg F-18-PSMA-11 with a scan duration of 3 min/bed position. To this group, no furosemide was administered. Images were scored on image quality using a 7-point scale and each suspicious lesion was described. To assess interrater reliability, two nuclear physicians scored all scans independently and described all observed suspicious lesions. Results The 4 MBq/kg group received for all reconstructed images (60 min p.i., 1.5 and 3 min/bed position and 180 min p.i., 1.5 and 3 min/bed position) the highest median image quality score compared to the 2 MBq/kg group (p values < 0.01). When comparing all reconstructed images, the highest image quality score was given to images at 60 min p.i., 3 min/bed position for both dosage groups (score 5 and 6 for 2 and 4 MBq/kg, respectively). The addition of furosemide administration decreased the interference score with one point (p = 0.01106) and facilitated the evaluation of lesions in proximity to the ureters. The interrater reliability for the comparison of each lesion separately after more than 40 F-18-PSMA-11 scan readings showed an increasing kappa value from 0.78 (95% CI, 0.65-0.92) to 0.94 (95% CI, 0.87-1). Conclusion Although the results indicate an administered activity of 4.0 +/- 0.4 MBq/kg, preference will be given to 2.0 +/- 0.2 MBq/kg due to the small difference in absolute score (max 1 point) and the ALARA principle. For evaluation of lesions in proximity to the ureters, the co-administration of a diuretic can be useful. The increase of the kappa value from 0.78 to 0.94 suggests a learning curve in the interpretation of F-18-PSMA-11 images

NEMA NU 2-2007 performance characteristics of GE Signa integrated PET/MR for different PET isotopes

BackgroundFully integrated PET/MR systems are being used frequently in clinical research and routine. National Electrical Manufacturers Association (NEMA) characterization of these systems is generally done with F-18 which is clinically the most relevant PET isotope. However, other PET isotopes, such as Ga-68 and Y-90, are gaining clinical importance as they are of specific interest for oncological applications and for follow-up of Y-90-based radionuclide therapy. These isotopes have a complex decay scheme with a variety of prompt gammas in coincidence. Ga-68 and Y-90 have higher positron energy and, because of the larger positron range, there may be interference with the magnetic field of the MR compared to F-18. Therefore, it is relevant to determine the performance of PET/MR for these clinically relevant and commercially available isotopes.MethodsNEMA NU 2-2007 performance measurements were performed for characterizing the spatial resolution, sensitivity, image quality, and the accuracy of attenuation and scatter corrections for F-18, Ga-68, and Y-90. Scatter fraction and noise equivalent count rate (NECR) tests were performed using F-18 and Ga-68. All phantom data were acquired on the GE Signa integrated PET/MR system, installed in UZ Leuven, Belgium.Results(18)F, Ga-68, and Y-90 NEMA performance tests resulted in substantially different system characteristics. In comparison with F-18, the spatial resolution is about 1mm larger in the axial direction for Ga-68 and no significative effect was found for Y-90. The impact of this lower resolution is also visible in the recovery coefficients of the smallest spheres of Ga-68 in image quality measurements, where clearly lower values are obtained. For Y-90, the low number of counts leads to a large variability in the image quality measurements. The primary factor for the sensitivity change is the scale factor related to the positron emission fraction. There is also an impact on the peak NECR, which is lower for Ga-68 than for F-18 and appears at higher activities.ConclusionsThe system performance of GE Signa integrated PET/MR was substantially different, in terms of NEMA spatial resolution, image quality, and NECR for Ga-68 and Y-90 compared to F-18. But these differences are compensated by the PET/MR scanner technologies and reconstructions methods

Performance of three-photon PET imaging: Monte Carlo simulations

We have recently introduced the idea of making use of three-photon positron annihilations in positron emission tomography. In this paper the basic characteristics of the three-gamma imaging in PET are studied by means of Monte Carlo simulations and analytical computations. Two typical configurations of human and small animal scanners are considered. Three-photon imaging requires high energy resolution detectors. Parameters currently attainable by CdZnTe semiconductor detectors, the technology of choice for the future development of radiation imaging, are assumed. Spatial resolution is calculated as a function of detector energy resolution and size, position in the field of view, scanner size, and the energies of the three gamma annihilation photons. Possible ways to improve the spatial resolution obtained for nominal parameters: 1.5 cm and 3.2 mm FWHM for human and small animal scanners, respectively, are indicated. Counting rates of true and random three-photon events for typical human and small animal scanning configurations are assessed. A simple formula for minimum size of lesions detectable in the three-gamma based images is derived. Depending on the contrast and total number of registered counts, lesions of a few mm size for human and sub mm for small animal scanners can be detected

Texture analysis of aggressive and nonaggressive lung tumor CE CT images

This paper presents the potential for fractal analysis of time sequence contrast-enhanced (CE) computed tomography (CT) images to differentiate between aggressive and nonaggressive malignant lung tumors (i.e., high and low metabolic tumors). The aim is to enhance CT tumor staging prediction accuracy through identifying malignant aggressiveness of lung tumors. As branching of blood vessels can be considered a fractal process, the research examines vascularized tumor regions that exhibit strong fractal characteristics. The analysis is performed after injecting 15 patients with a contrast agent and transforming at least 11 time sequence CE CT images from each patient to the fractal dimension and determining corresponding lacunarity. The fractal texture features were averaged over the tumor region and quantitative classification showed up to 83.3% accuracy in distinction between advanced (aggressive) and early-stage (nonaggressive) malignant tumors. Also, it showed strong correlation with corresponding lung tumor stage and standardized tumor uptake value of fluoro deoxyglucose as determined by positron emission tomography. These results indicate that fractal analysis of time sequence CE CT images of malignant lung tumors could provide additional information about likely tumor aggression that could potentially impact on clinical management decisions in choosing the appropriate treatment procedure

MITK-ModelFit: A generic open-source framework for model fits and their exploration in medical imaging -- design, implementation and application on the example of DCE-MRI

Many medical imaging techniques utilize fitting approaches for quantitative parameter estimation and analysis. Common examples are pharmacokinetic modeling in DCE MRI/CT, ADC calculations and IVIM modeling in diffusion-weighted MRI and Z-spectra analysis in chemical exchange saturation transfer MRI. Most available software tools are limited to a special purpose and do not allow for own developments and extensions. Furthermore, they are mostly designed as stand-alone solutions using external frameworks and thus cannot be easily incorporated natively in the analysis workflow. We present a framework for medical image fitting tasks that is included in MITK, following a rigorous open-source, well-integrated and operating system independent policy. Software engineering-wise, the local models, the fitting infrastructure and the results representation are abstracted and thus can be easily adapted to any model fitting task on image data, independent of image modality or model. Several ready-to-use libraries for model fitting and use-cases, including fit evaluation and visualization, were implemented. Their embedding into MITK allows for easy data loading, pre- and post-processing and thus a natural inclusion of model fitting into an overarching workflow. As an example, we present a comprehensive set of plug-ins for the analysis of DCE MRI data, which we validated on existing and novel digital phantoms, yielding competitive deviations between fit and ground truth. Providing a very flexible environment, our software mainly addresses developers of medical imaging software that includes model fitting algorithms and tools. Additionally, the framework is of high interest to users in the domain of perfusion MRI, as it offers feature-rich, freely available, validated tools to perform pharmacokinetic analysis on DCE MRI data, with both interactive and automatized batch processing workflows.Comment: 31 pages, 11 figures URL: http://mitk.org/wiki/MITK-ModelFi