Neural blind deconvolution for deblurring and supersampling PSMA PET.

Objective. To simultaneously deblur and supersample prostate specific membrane antigen (PSMA) positron emission tomography (PET) images using neural blind deconvolution.Approach. Blind deconvolution is a method of estimating the hypothetical 'deblurred' image along with the blur kernel (related to the point spread function) simultaneously. Traditionalmaximum a posterioriblind deconvolution methods require stringent assumptions and suffer from convergence to a trivial solution. A method of modelling the deblurred image and kernel with independent neural networks, called 'neural blind deconvolution' had demonstrated success for deblurring 2D natural images in 2020. In this work, we adapt neural blind deconvolution to deblur PSMA PET images while simultaneous supersampling to double the original resolution. We compare this methodology with several interpolation methods in terms of resultant blind image quality metrics and test the model's ability to predict accurate kernels by re-running the model after applying artificial 'pseudokernels' to deblurred images. The methodology was tested on a retrospective set of 30 prostate patients as well as phantom images containing spherical lesions of various volumes.Main results. Neural blind deconvolution led to improvements in image quality over other interpolation methods in terms of blind image quality metrics, recovery coefficients, and visual assessment. Predicted kernels were similar between patients, and the model accurately predicted several artificially-applied pseudokernels. Localization of activity in phantom spheres was improved after deblurring, allowing small lesions to be more accurately defined.Significance. The intrinsically low spatial resolution of PSMA PET leads to partial volume effects (PVEs) which negatively impact uptake quantification in small regions. The proposed method can be used to mitigate this issue, and can be straightforwardly adapted for other imaging modalities.

Development of the quantitative PET prostate phantom (Q3P) for improved quality assurance of 18F-PSMA PET imaging in metastatic prostate cancer.

BACKGROUND: Phantoms are commonly used to evaluate and compare the performance of imaging systems given the known ground truth. Positron emission tomography (PET) scanners are routinely validated using the NEMA image quality phantom, in which lesions are modeled using 10 to 37 mm fillable spheres. The NEMA phantom neglects, however, to model focal (3-10-mm), high-uptake lesions that are increasingly observed in prostate-specific membrane antigen (PSMA) PET images. PSMA-targeting radiopharmaceuticals allow for enhanced detection of metastatic prostate cancers. As such, there is significant need to develop an updated phantom which considers both the quantitative and lesion detectability of this new paradigm in oncological PET imaging. PURPOSE: In this work, we present the Quantitative PET Prostate Phantom (Q3P); a portable and modular phantom that can be used to improve and harmonize imaging protocols for 18F-PSMA PET scans. METHODS: A one-piece cylindrical phantom was designed effectively in two halves, which we call modules. Module 1 was designed to mimic lesions in the presence of background, and Module 2 mimicked very high contrast conditions (i.e., very low background) that can be observed in 18F-PSMA PET scans. Shell-less radioactive spheres (3-16-mm) were cast using epoxy resin mixed with sodium-22 (22Na), a long half-life positron emitter with positron range similar to 18F. To establish realistic lesion contrast, the 22Na spheres were mounted in a cylindrical chamber that can be filled with an 18F background (module 1). Thirteen exchangeable spherical cavity inserts (3-37-mm) were machined in two parts and solvent welded together, and filled with 18F (50 kBq/mL) to model lesions with very high contrast (module 2). Five 2.5-min PET scans were acquired on a 5-ring GE Discovery MI PET/CT scanner (General Electric, USA). Lesions were segmented using 41% of SUVmax fixed thresholding (41% FT) and recovery coefficients (RCs) were computed from 5 noise realizations. RESULTS: The manufactured phantom is portable (5.7 kg) and scan preparation takes less than 40 min. The total 22Na activity is 250 kBq, allowing it to be shipped as an exempt package under International Atomic Energy Agency (IAEA) regulations. Recovery coefficients, computed using PSF modeling and no post-reconstruction smoothing, were 130.3% (16 mm), 147.1% (10 mm), 87.2% (6 mm), and 7.0% (3 mm) for RCmax, which decreased to 91.1% (16 mm), 90.6% (10 mm), 53.2% (6 mm), and 3.6% (3 mm) for RCmean in the 22Na spheres. Comparatively, 18F sphere recovery was 110.7% (17 mm), 123.6% (10 mm), 106.5% (7 mm), and 23.3% (3 mm) for RCmax, which was reduced to 76.7% (17 mm), 77.7% (10 mm), 66.8% (7 mm), and 13.5% (3 mm), for RCmean. CONCLUSIONS: A standardized imaging phantom was developed for lesion quantification assessment in 18F-PSMA PET images. The phantom is configurable, providing users with the opportunity to modify background activity levels or sphere sizes according to clinical demands. Distributed to the community, the Q3P phantom has the potential to enable better assessment of lesion quantification and harmonization of 18F-PSMA PET imaging, which may lead to more robust predictive metrics and better outcome prediction in metastatic prostate cancer.

Development of scalable lymphatic system in the 4D XCAT phantom: Application to quantitative evaluation of lymphoma PET segmentations.

BACKGROUND: Digital anthropomorphic phantoms, such as the 4D extended cardiac-torso (XCAT) phantom, are actively used to develop, optimize, and evaluate a variety of imaging applications, allowing for realistic patient modeling and knowledge of ground truth. The XCAT phantom defines the activity and attenuation for a simulated patient, which includes a complete set of organs, muscle, bone, and soft tissue, while also accounting for cardiac and respiratory motion. However, the XCAT phantom does not currently include the lymphatic system, critical for evaluating medical imaging tasks such as sentinel node detection, node density measurement, and radiation dosimetry. PURPOSE: In this study, we aimed to develop a scalable lymphatic system in the XCAT phantom, to facilitate improved research of the lymphatic system in medical imaging. Using this scalable lymphatic system, we modeled the lymph node conglomerate pathology that is characteristically observed in primary mediastinal B-cell lymphoma (PMBCL). As an extended application, we evaluated positron emission tomography (PET) image quantification of metabolic tumor volume (MTV) and total lesion glycolysis (TLG) of these simulated lymphomas, though the phantoms may be applied to other imaging modalities and study design paradigms (e.g., image quality, detection). METHODS: A template model for the lymphatic system was developed based on anatomical data from the Visible Human Project of the National Library of Medicine. The segmented nodes and vessels were fit with non-uniform rational basis spline surfaces, and multichannel large deformation diffeomorphic metric mapping was used to propagate the template to different XCAT anatomies. To model conglomerates observed in PMBCL, lymph nodes were enlarged, converged within the mediastinum, and tracer concentration was increased. We used the phantoms as inputs to a PET simulation tool, which generated images using ordered subsets expectation maximization reconstruction with 2-8 mm Gaussian filters. Fixed thresholding (FT) and gradient segmentation were used to determine MTV and TLG. Percent bias (%Bias) and coefficient of variation (COV) were computed as measures of accuracy and precision, respectively, for each MTV and TLG measurement. RESULTS: Using the methodology described above, we introduced a scalable lymphatic system in the XCAT phantom, which allows for the radioactivity and attenuation ground truth to be generated in 116 ± 2.5 s using a 2.3 GHz processor. Within the Rhinoceros interface, lymph node anatomy and function were modified to create a cohort of 10 phantoms with lymph node conglomerates. Using the lymphoma phantoms to evaluate PET quantification of MTV, mean %Bias values were -9.3%, -41.3%, and 20.9%, while COV values were 4.08%, 7.6%, and 3.4% using 25% FT, 40% FT, and gradient segmentations, respectively. Comparatively for TLG, mean %Bias values were -27.4%, -45.8%, and -16.0%, while COV values were 1.9%, 5.7%, and 1.4%, for the 25% FT, 40% FT, and gradient segmentations, respectively. CONCLUSIONS: In this work, we upgraded the XCAT phantom to include a lymphatic system, comprised of a network of 276 scalable lymph nodes and corresponding vessels. As an application, we created a cohort of phantoms with lymph node conglomerates to evaluate lymphoma quantification in PET imaging, which highlights an important application of this work.

Quantitative evaluation of PSMA PET imaging using a realistic anthropomorphic phantom and shell-less radioactive epoxy lesions.

BACKGROUND: Positron emission tomography (PET) with prostate specific membrane antigen (PSMA) have shown superior performance in detecting metastatic prostate cancers. Relative to [18F]fluorodeoxyglucose ([18F]FDG) PET images, PSMA PET images tend to visualize significantly higher-contrast focal lesions. We aim to evaluate segmentation and reconstruction algorithms in this emerging context. Specifically, Bayesian or maximum a posteriori (MAP) image reconstruction, compared to standard ordered subsets expectation maximization (OSEM) reconstruction, has received significant interest for its potential to reach convergence with minimal noise amplifications. However, few phantom studies have evaluated the quantitative accuracy of such reconstructions for high contrast, small lesions (sub-10 mm) that are typically observed in PSMA images. In this study, we cast 3 mm-16-mm spheres using epoxy resin infused with a long half-life positron emitter (sodium-22; 22Na) to simulate prostate cancer metastasis. The anthropomorphic Probe-IQ phantom, which features a liver, bladder, lungs, and ureters, was used to model relevant anatomy. Dynamic PET acquisitions were acquired and images were reconstructed with OSEM (varying subsets and iterations) and BSREM (varying ß parameters), and the effects on lesion quantitation were evaluated. RESULTS: The 22Na lesions were scanned against an aqueous solution containing fluorine-18 (18F) as the background. Regions-of-interest were drawn with MIM Software using 40% fixed threshold (40% FT) and a gradient segmentation algorithm (MIM's PET Edge+). Recovery coefficients (RCs) (max, mean, peak, and newly defined "apex"), metabolic tumour volume (MTV), and total tumour uptake (TTU) were calculated for each sphere. SUVpeak and SUVapex had the most consistent RCs for different lesion-to-background ratios and reconstruction parameters. The gradient-based segmentation algorithm was more accurate than 40% FT for determining MTV and TTU, particularly for lesions [Formula: see text] 6 mm in diameter (R2 = 0.979-0.996 vs. R2 = 0.115-0.527, respectively). CONCLUSION: An anthropomorphic phantom was used to evaluate quantitation for PSMA PET imaging of metastatic prostate cancer lesions. BSREM with ß = 200-400 and OSEM with 2-5 iterations resulted in the most accurate and robust measurements of SUVmean, MTV, and TTU for imaging conditions in 18F-PSMA PET/CT images. SUVapex, a hybrid metric of SUVmax and SUVpeak, was proposed for robust, accurate, and segmentation-free quantitation of lesions for PSMA PET.

An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 1: Design and Early Results of the SNMMI Dosimetry Challenge.

In this work, we present details and initial results from a 177Lu dosimetry challenge that has been designed to collect data from the global nuclear medicine community aiming at identifying, understanding, and quantitatively characterizing the consequences of the various sources of variability in dosimetry. Methods: The challenge covers different approaches to performing dosimetry: planar, hybrid, and pure SPECT. It consists of 5 different and independent tasks to measure the variability of each step in the dosimetry workflow. Each task involves the calculation of absorbed doses to organs and tumors and was meant to be performed in sequential order. The order of the tasks is such that results from a previous one would not affect subsequent ones. Different sources of variability are removed as the participants advance through the challenge by giving them the data required to begin the calculations at different steps of the dosimetry workflow. Data from 2 patients after a therapeutic administration of 177Lu-DOTATATE were used for this study. The data are hosted in Deep Blue Data, a data repository service run by the University of Michigan. Participants submit results in standardized spreadsheets and with a short description summarizing their methods. Results: In total, 178 participants have signed up for the challenge, and 119 submissions have been received. Sixty percent of submissions have used voxelized dose methods, with 47% of those using commercial software. In initial analysis, the volume of organs showed a variability of up to 49.8% whereas for lesions this was up to 176%. Variability in time-integrated activity was up to 192%. Mean absorbed doses varied up to 57.7%. Segmentation is the step that required the longest time to complete, with a median of 43 min. The median total time to perform the full calculation was 89 min. Conclusion: To advance dosimetry and encourage its routine use in radiopharmaceutical therapy applications, it is critical that dosimetry results be reproducible across centers. Our initial results provide insights into the variability associated with performing dose calculations. It is expected that this dataset, including results from future stages, will result in efforts to standardize and harmonize methods and procedures.

Design of an anthropomorphic PET phantom with elastic lungs and respiration modeling.

PURPOSE: Respiratory motion during positron emission tomography (PET) scans can be a major detriment to image quality in oncological imaging. The impact of motion on lesion quantification and detectability can be assessed using phantoms with realistic anatomy representation and motion modeling. In this work, we develop an anthropomorphic phantom for PET imaging that combines anatomic fidelity and a realistic breathing mechanism with deformable lungs. METHODS: We start from a previously developed anatomically accurate but static phantom of a human torso, and add elastic lungs with a highly controllable actuation mechanism which replicates the physics of breathing. The space outside the lungs is filled with a radioactive water solution. To maintain anatomical accuracy and realistic gamma ray attenuation in the torso, all motion mechanisms and actuators are positioned outside of the phantom compartment. The actuation mechanism can produce custom respiratory waveforms with breathing rates up to 25 breaths per minute and tidal volumes up to 1200 mL. RESULTS: Several tests were performed to validate the performance of the phantom assembly, in which the phantom was filled with water and given respiratory waveforms to execute. All parts demonstrated expected performance. Force requirements were not exceeded and no leaks were detected, although continued use of the phantom is required to evaluate wear. The motion of the lungs was determined to be within a reasonable realistic range. CONCLUSIONS: The full mechanical design is described in this paper, as well as a software application with graphical user interface which was developed to plan and visualize respiratory patterns. Both are available online as open source files. The developed phantom will facilitate future work in evaluating the impact of respiratory motion on lesion quantification and detectability in clinical practice.

Monte Carlo simulation of 6-MV dynamic wave VMAT deliveries by Vero4DRT linear accelerator using EGSnrc moving sources.

The commissioning and benchmark of a Monte Carlo (MC) model of the 6-MV Brainlab-Mitsubishi Vero4DRT linear accelerator for the purpose of quality assurance of clinical dynamic wave arc (DWA) treatment plans is reported. Open-source MC applications based on EGSnrc particle transport codes are used to simulate the medical linear accelerator head components. Complex radiotherapy irradiations can be simulated in a single MC run using a shared library format combined with BEAMnrc "source20." Electron energy tuning is achieved by comparing measured vs simulated percentage depth doses (PDDs) for MLC-defined field sizes in a water phantom. Electron spot size tuning is achieved by comparing measured and simulated inplane and crossplane beam profiles. DWA treatment plans generated from RayStation (RaySearch) treatment planning system (TPS) are simulated on voxelized (2.5 mm3 ) patient CT datasets. Planning target volume (PTV) and organs at risk (OAR) dose-volume histograms (DVHs) are compared to TPS-calculated doses for clinically deliverable dynamic volumetric modulated arc therapy (VMAT) trajectories. MC simulations with an electron beam energy of 5.9 MeV and spot size FWHM of 1.9 mm had the closest agreement with measurement. DWA beam deliveries simulated on patient CT datasets results in DVH agreement with TPS-calculated doses. PTV coverage agreed within 0.1% and OAR max doses (to 0.035 cc volume) agreed within 1 Gy. This MC model can be used as an independent dose calculation from the TPS and as a quality assurance tool for complex, dynamic radiotherapy treatment deliveries. Full patient CT treatment simulations are performed in a single Monte Carlo run in 23 min. Simulations are run in parallel using the Condor High-Throughput Computing software1 on a cluster of eight servers. Each server has two physical processors (Intel Xeon CPU E5-2650 0 @2.00 GHz), with 8 cores per CPU and two threads per core for 256 calculation nodes.