Intraprocedural C-arm dual-phase cone-beam enhancement patterns correlate with tumor absorbed dose after radioembolization.

BACKGROUND: Recent studies have shown a clear relationship between absorbed dose and tumor response to treatment after hepatic radioembolization. These findings help to create more personalized treatment planning and dosimetry. However, crucial to this goal is the ability to predict the dose distribution prior to treatment. The microsphere distribution is ultimately determined by (i) the hepatic vasculature and the resulting blood flow dynamics and (ii) the catheter position. PURPOSE: To show that pretreatment, intra-procedural imaging of blood flow patterns, as quantified by catheter-directed intra-arterial contrast enhancement, correlate with posttreatment microsphere accumulation and, consequently, absorbed dose. MATERIALS AND METHODS: Patients who participated in a clinical trial (NCT01177007) and for whom both a pretreatment dual-phase contrast-enhanced cone-beam CT (CBCT) and a posttreatment 90Y PET/CT scan were available were included in this retrospective study. Tumors and perfused volumes were manually delineated on the CBCT by an experienced radiologist. The mean, sum, and standard deviation of the voxels in each volume were recorded. The delineations were transferred to the PET-based absorbed dose maps by coregistration of the corresponding CTs. Linear multiple regression was used to correlate pretreatment CBCT enhancement to posttreatment 90Y PET/CT-based absorbed dose in each region. Leave-one-out cross-validation and Bland-Altman analyses were performed on the predicted versus measured absorbed doses. RESULTS: Nine patients, with a total of 23 tumors were included. All presented with hepatocellular carcinoma (HCC). Visually, all patients had a clear correspondence between CBCT enhancement and absorbed dose. The correlation between CBCT enhancement and posttherapy absorbed tumor dose based was strong (R2 = 0.91), and moderate for the non-tumor liver tissue (R2 = 0.61). Limits of agreement were approximately ±55 Gray for tumor tissue. CONCLUSION: There is a linear relationship between pretreatment blood dynamics in HCC tumors and posttreatment absorbed dose, which, if shown to be generalizable, allows for pretreatment tumor absorbed dose prediction.

Girth-based administered activity for pediatric 99m Tc-DMSA SPECT.

BACKGROUND: Pediatric molecular imaging requires a balance between administering an activity that will yield sufficient diagnostic image quality while maintaining patient radiation exposure at acceptable levels. In current clinical practice, this balance is arrived at by the current North American Consensus Guidelines in which patient weight is used to recommend the administered activity (AA). PURPOSE: We have previously demonstrated that girth (waist circumference at the level of the kidneys) is better at equalizing image quality than patient weight for pediatric Tc-99m DMSA renal function imaging. However, the correlation between image quality (IQ), AA, and patient girth has not been rigorously and systematically developed. In this work, we generate a series of curves showing the tradeoff between AA and IQ as a function of patient girth, providing the data for standards bodies to develop the next generation of dosing guideline for pediatric DMSA SPECT. METHODS: An anthropomorphic phantom series that included variations in age (5, 10, and 15 years), gender (M, F), local body morphometry (5, 10, 50, 90, and 95th girth percentiles), and kidney size (±15% standard size), was used to generate realistic SPECT projections. A fixed and clinically challenging defect-to-organ volume percentage (0.49% of renal cortex value) was used to model a focal defect with zero uptake (i.e., full local loss of renal function). Task-based IQ assessment methods were used to rigorously measure IQ in terms of renal perfusion defect detectability. This assessment was performed at multiple count levels (corresponding to various AAs) for groups of patients that had similar girths and defect sizes. Receiver-operating characteristics (ROC) analysis was applied; the area under the ROC curve (AUC) was used as a figure-of-merit for task performance. Curves showing the tradeoff between AUC and AA were generated for these groups of phantoms. RESULTS: Overall, the girth-based dosing method suggested different amounts of AA compared to weight-based dosing for the phantoms that had a relatively large body weight but a small girth or phantoms with relatively small bodyweight but large girth. Reductions of AA to 62.9% compared to weight-based dosing guidelines can potentially be realized while maintaining a baseline (AUC = 0.80) IQ for certain 15-year-olds who have a relatively small girth and large defect size. Note that the task-based IQ results are heavily dependent on the simulated defect size for the defect detection task and the appropriate AUC value must be decided by the physicians for this diagnostic task. These results are based purely on simulation and are subject to future clinical validation. CONCLUSIONS: The study provides simulation-based IQ-AA data for a girth-based dosing method for pediatric renal SPECT, suggesting that patient waist circumference at the level of kidneys should be considered in selecting the AA needed to achieve an acceptable IQ. This data may be useful for standards bodies to develop girth-based dosing guidelines.

Lack of repeatability of radiomic features derived from PET scans: Results from a 18 F-DCFPyL test-retest cohort.

OBJECTIVES: PET-based radiomic metrics are increasingly utilized as predictive image biomarkers. However, the repeatability of radiomic features on PET has not been assessed in a test-retest setting. The prostate-specific membrane antigen-targeted compound 18 F-DCFPyL is a high-affinity, high-contrast PET agent that we utilized in a test-retest cohort of men with metastatic prostate cancer (PC). METHODS: Data of 21 patients enrolled in a prospective clinical trial with histologically proven PC underwent two 18 F-DCFPyL PET scans within 7 days, using identical acquisition and reconstruction parameters. Sites of disease were segmented and a set of 29 different radiomic parameters were assessed on both scans. We determined repeatability of quantification by using Pearson's correlations, within-subject coefficient of variation (wCOV), and Bland-Altman analysis. RESULTS: In total, 230 lesions (177 bone, 38 lymph nodes, 15 others) were assessed on both scans. For all investigated radiomic features, a broad range of inter-scan correlation was found (r, 0.07-0.95), with acceptable reproducibility for entropy and homogeneity (wCOV, 16.0% and 12.7%, respectively). On Bland-Altman analysis, no systematic increase or decrease between the scans was observed for either parameter (±1.96 SD: 1.07/-1.30, 0.23/-0.18, respectively). The remaining 27 tested radiomic metrics, however, achieved unacceptable high wCOV (≥21.7%). CONCLUSION: Many common radiomic features derived from a test-retest PET study had poor repeatability. Only Entropy and homogeneity achieved good repeatability, supporting the notion that those image biomarkers may be incorporated in future clinical trials. Those radiomic features based on high frequency aspects of images appear to lack the repeatability on PET to justify further study.

TransMorph: Transformer for unsupervised medical image registration.

In the last decade, convolutional neural networks (ConvNets) have been a major focus of research in medical image analysis. However, the performances of ConvNets may be limited by a lack of explicit consideration of the long-range spatial relationships in an image. Recently, Vision Transformer architectures have been proposed to address the shortcomings of ConvNets and have produced state-of-the-art performances in many medical imaging applications. Transformers may be a strong candidate for image registration because their substantially larger receptive field enables a more precise comprehension of the spatial correspondence between moving and fixed images. Here, we present TransMorph, a hybrid Transformer-ConvNet model for volumetric medical image registration. This paper also presents diffeomorphic and Bayesian variants of TransMorph: the diffeomorphic variants ensure the topology-preserving deformations, and the Bayesian variant produces a well-calibrated registration uncertainty estimate. We extensively validated the proposed models using 3D medical images from three applications: inter-patient and atlas-to-patient brain MRI registration and phantom-to-CT registration. The proposed models are evaluated in comparison to a variety of existing registration methods and Transformer architectures. Qualitative and quantitative results demonstrate that the proposed Transformer-based model leads to a substantial performance improvement over the baseline methods, confirming the effectiveness of Transformers for medical image registration.

Automated segmentation of lung, liver, and liver tumors from Tc-99m MAA SPECT/CT images for Y-90 radioembolization using convolutional neural networks.

PURPOSE: 90 Y selective internal radiation therapy (SIRT) has become a safe and effective treatment option for liver cancer. However, segmentation of target and organ-at-risks is labor-intensive and time-consuming in 90 Y SIRT planning. In this study, we developed a convolutional neural network (CNN)-based method for automated lungs, liver, and tumor segmentation on 99m Tc-MAA SPECT/CT images for 90 Y SIRT planning. METHODS: 99m Tc-MAA SPECT/CT images and corresponding clinical segmentations were retrospectively collected from 56 patients who underwent 90 Y SIRT. The collected data were used to train three CNN-based segmentation algorithms for lungs, liver, and tumor segmentation. Segmentation performance was evaluated using the Dice similarity coefficient (DSC), surface DSC, and average symmetric surface distance (ASSD). Dosimetric parameters (volume, counts, and lung shunt fraction) were measured from the segmentation results and were compared with clinical reference segmentations. RESULTS: The evaluation results show that the method can accurately segment lungs, liver, and tumor with median [interquartile range] DSCs of 0.98 [0.97-0.98], 0.91 [0.83-0.93], and 0.85 [0.71-0.88]; surface DSCs of 0.99 [0.97-0.99], 0.86 [0.77-0.93], and 0.85 [0.62-0.93], and ASSDs of 0.91 [0.69-1.5], 4.8 [2.6-8.4], and 4.7 [3.5-9.2] mm, respectively. Dosimetric parameters from the three segmentation networks show relationship with those from the reference segmentations. The overall segmentation took about 1 min per patient on an NVIDIA RTX-2080Ti GPU. CONCLUSION: This work presents CNN-based algorithms to segment lungs, liver, and tumor from 99m Tc-MAA SPECT/CT images. The results demonstrated the potential of the proposed CNN-based segmentation method for assisting 90 Y SIRT planning while drastically reducing operator time.

Dosimetric impact of Ac-227 in accelerator-produced Ac-225 for alpha-emitter radiopharmaceutical therapy of patients with hematological malignancies: a pharmacokinetic modeling analysis.

BACKGROUND: Actinium-225 is an alpha-particle emitter under investigation for use in radiopharmaceutical therapy. To address limited supply, accelerator-produced 225Ac has been recently made available. Accelerator-produced 225Ac via 232Th irradiation (denoted 225/7Ac) contains a low percentage (0.1-0.3%) of 227Ac (21.77-year half-life) activity at end of bombardment. Using pharmacokinetic modeling, we have examined the dosimetric impact of 227Ac on the use of accelerator-produced 225Ac for radiopharmaceutical therapy. We examine the contribution of 227Ac and its daughters to tissue absorbed doses. The dosimetric analysis was performed for antibody-conjugated 225/7Ac administered intravenously to treat patients with hematological cancers. Published pharmacokinetic models are used to obtain the distribution of 225/7Ac-labeled antibody and also the distribution of either free or antibody-conjugated 227Th. RESULTS: Based on our modeling, the tissue specific absorbed dose from 227Ac would be negligible in the context of therapy, less than 0.02 mGy/MBq for the top 6 highest absorbed tissues and less than 0.007 mGy/MBq for all other tissues. Compared to that from 225Ac, the absorbed dose from 227Ac makes up a very small component (less than 0.04%) of the total absorbed dose delivered to the 6 highest dose tissues: red marrow, spleen, endosteal cells, liver, lungs and kidneys when accelerator produced 225/7Ac-conjugated anti-CD33 antibody is used to treat leukemia patients. For all tissues, the dominant contributor to the absorbed dose arising from the 227Ac is 227Th, the first daughter of 227Ac which has the potential to deliver absorbed dose both while it is antibody-bound and while it is free. CONCLUSIONS: These results suggest that the absorbed dose arising from 227Ac to normal organs would be negligible for an 225/7Ac-labeled antibody that targets hematological cancer.

Prospective SPECT-CT Organ Dosimetry-Driven Radiation-Absorbed Dose Escalation Using the In-111 (111In)/Yttrium 90 (90Y) Ibritumomab Tiuxetan (Zevalin®) Theranostic Pair in Patients with Lymphoma at Myeloablative Dose Levels.

PURPOSE: We prospectively evaluated the feasibility of SPECT-CT/planar organ dosimetry-based radiation dose escalation radioimmunotherapy in patients with recurrent non-Hodgkin's lymphoma using the theranostic pair of 111In and 90Y anti-CD20 ibritumomab tiuxetan (Zevalin®) at myeloablative radiation-absorbed doses with autologous stem cell support. We also assessed acute non-hematopoietic toxicity and early tumor response in this two-center outpatient study. METHODS: 24 patients with CD20-positive relapsed or refractory rituximab-sensitive, low-grade, mantle cell, or diffuse large-cell NHL, with normal organ function, platelet counts > 75,000/mm3, and <35% tumor involvement in the marrow were treated with Rituximab (375 mg/m2) weekly for 4 consecutive weeks, then one dose of cyclophosphamide 2.5 g/m2 with filgrastim 10 mcg/kg/day until stem cell collection. Of these, 18 patients with successful stem cell collection (at least 2 × 106 CD34 cells/kg) proceeded to RIT. A dosimetric administration of 111In ibritumomab tiuxetan (185 MBq) followed by five sequential quantitative planar and one SPECT/CT scan was used to determine predicted organ radiation-absorbed dose. Two weeks later, 90Y ibritumomab tiuxetan was administered in an outpatient setting at a cohort- and patient-specific predicted organ radiation-absorbed dose guided by a Continuous Response Assessment (CRM) methodology with the following cohorts for dose escalation: 14.8 MBq/kg, and targeted 18, 24, 28, and 30.5 Gy to the liver. Autologous stem cell infusion occurred when the estimated marrow radiation-absorbed dose rate was predicted to be <1 cGy/h. Feasibility, short-term toxicities, and tumor response were assessed. RESULTS: Patient-specific hybrid SPECT/CT + planar organ dosimetry was feasible in all 18 cases and used to determine the patient-specific therapeutic dose and guide dose escalation (26.8 ± 7.3 MBq/kg (mean), 26.3 MBq/kg (median) of 90Y (range: 12.1-41.4 MBq/kg)) of ibritumomab tiuxetan that was required to deliver 10 Gy to the liver. Infused stem cells engrafted rapidly. The most common treatment-related toxicities were hematological and were reversible following stem cell infusion. No significant hepatotoxicity was seen. One patient died from probable treatment-related causes-pneumonia at day 27 post-transplant. One patient at dose level 18 Gy developed myelodysplastic syndrome (MDS), 4 patients required admission post-90Y RIT for febrile neutropenia, 16/18 patients receiving 90Y ibritumomab tiuxetan (89%) responded to the therapy, with 13 CR (72%) and 3/18 PR (17%), at 60 days post-treatment. Two patients had progressive disease at sixty days. One patient was lost to follow-up. Median time to progression was estimated to be at least 13 months. MTD to the liver is greater than 28 Gy, but the MTD was not reached as the study was terminated due to unexpected discontinuation of availability of the therapeutic agent. CONCLUSIONS: Patient-specific outpatient 90Y ibritumomab tiuxetan RIT with myeloablative doses of RIT up to a targeted 30.5 Gy to the liver is feasible, guided by prospective SPECT/CT + planar imaging with the theranostic pair of 111In and 90Y anti-CD20, with outpatient autologous stem cell transplant support. Administered activity over 5 times the standard FDA-approved activity was well-tolerated. The non-hematopoietic MTD in this study exceeds 28 Gy to the liver. Initial tumor responses were common at all dose levels. This study supports the feasibility of organ dosimetry-driven patient-specific dose escalation in the treatment of NHL with stem cell transplant and provides additional information on the radiation tolerance of the normal liver to radiopharmaceutical therapy.

Renal 99mTc-DMSA pharmacokinetics in pediatric patients.

99mTc-DMSA is one of the most commonly used pediatric nuclear medicine imaging agents. Nevertheless, there are no pharmacokinetic (PK) models for 99mTc-DMSA in children, and currently available pediatric dose estimates for 99mTc-DMSA use pediatric S values with PK data derived from adults. Furthermore, the adult PK data were collected in the mid-70's using quantification techniques and instrumentation available at the time. Using pediatric imaging data for DMSA, we have obtained kinetic parameters for DMSA that differ from those applicable to adults. METHODS: We obtained patient data from a retrospective re-evaluation of clinically collected pediatric SPECT images of 99mTc-DMSA in 54 pediatric patients from Boston's Children Hospital (BCH), ranging in age from 1 to 16 years old. These were supplemented by prospective data from twenty-three pediatric patients (age range: 4 months to 6 years old). RESULTS: In pediatric patients, the plateau phase in fractional kidney uptake occurs at a fractional uptake value closer to 0.3 than the value of 0.5 reported by the International Commission on Radiological Protection (ICRP) for adult patients. This leads to a 27% lower time-integrated activity coefficient in pediatric patients than in adults. Over the age range examined, no age dependency in uptake fraction at the clinical imaging time was observed. Female pediatric patients had a 17% higher fractional kidney uptake at the clinical imaging time than males (P < 0.001). CONCLUSIONS: Pediatric 99mTc-DMSA kinetics differ from those reported for adults and should be considered in pediatric patient dosimetry. Alternatively, the differences obtained in this study could reflect improved quantification methods and the need to re-examine DMSA kinetics in adults.

Learning fuzzy clustering for SPECT/CT segmentation via convolutional neural networks.

PURPOSE: Quantitative bone single-photon emission computed tomography (QBSPECT) has the potential to provide a better quantitative assessment of bone metastasis than planar bone scintigraphy due to its ability to better quantify activity in overlapping structures. An important element of assessing the response of bone metastasis is accurate image segmentation. However, limited by the properties of QBSPECT images, the segmentation of anatomical regions-of-interests (ROIs) still relies heavily on the manual delineation by experts. This work proposes a fast and robust automated segmentation method for partitioning a QBSPECT image into lesion, bone, and background. METHODS: We present a new unsupervised segmentation loss function and its semi- and supervised variants for training a convolutional neural network (ConvNet). The loss functions were developed based on the objective function of the classical Fuzzy C-means (FCM) algorithm. The first proposed loss function can be computed within the input image itself without any ground truth labels, and is thus unsupervised; the proposed supervised loss function follows the traditional paradigm of the deep learning-based segmentation methods and leverages ground truth labels during training. The last loss function is a combination of the first and the second and includes a weighting parameter, which enables semi-supervised segmentation using deep learning neural network. EXPERIMENTS AND RESULTS: We conducted a comprehensive study to compare our proposed methods with ConvNets trained using supervised, cross-entropy and Dice loss functions, and conventional clustering methods. The Dice similarity coefficient (DSC) and several other metrics were used as figures of merit as applied to the task of delineating lesion and bone in both simulated and clinical SPECT/CT images. We experimentally demonstrated that the proposed methods yielded good segmentation results on a clinical dataset even though the training was done using realistic simulated images. On simulated SPECT/CT, the proposed unsupervised model's accuracy was greater than the conventional clustering methods while reducing computation time by 200-fold. For the clinical QBSPECT/CT, the proposed semi-supervised ConvNet model, trained using simulated images, produced DSCs of 0.75 and 0.74 for lesion and bone segmentation in SPECT, and a DSC of 0.79 bone segmentation of CT images. These DSCs were larger than that for standard segmentation loss functions by > 0.4 for SPECT segmentation, and > 0.07 for CT segmentation with P-values < 0.001 from a paired t-test. CONCLUSIONS: A ConvNet-based image segmentation method that uses novel loss functions was developed and evaluated. The method can operate in unsupervised, semi-supervised, or fully-supervised modes depending on the availability of annotated training data. The results demonstrated that the proposed method provides fast and robust lesion and bone segmentation for QBSPECT/CT. The method can potentially be applied to other medical image segmentation applications.

DeepAMO: a multi-slice, multi-view anthropomorphic model observer for visual detection tasks performed on volume images.

Purpose: We propose a deep learning-based anthropomorphic model observer (DeepAMO) for image quality evaluation of multi-orientation, multi-slice image sets with respect to a clinically realistic 3D defect detection task. Approach: The DeepAMO is developed based on a hypothetical model of the decision process of a human reader performing a detection task using a 3D volume. The DeepAMO is comprised of three sequential stages: defect segmentation, defect confirmation (DC), and rating value inference. The input to the DeepAMO is a composite image, typical of that used to view 3D volumes in clinical practice. The output is a rating value designed to reproduce a human observer's defect detection performance. In stages 2 and 3, we propose: (1) a projection-based DC block that confirms defect presence in two 2D orthogonal orientations and (2) a calibration method that "learns" the mapping from the features of stage 2 to the distribution of observer ratings from the human observer rating data (thus modeling inter- or intraobserver variability) using a mixture density network. We implemented and evaluated the DeepAMO in the context of Tc 99 m -DMSA SPECT imaging. A human observer study was conducted, with two medical imaging physics graduate students serving as observers. A 5 × 2 -fold cross-validation experiment was conducted to test the statistical equivalence in defect detection performance between the DeepAMO and the human observer. We also compared the performance of the DeepAMO to an unoptimized implementation of a scanning linear discriminant observer (SLDO). Results: The results show that the DeepAMO's and human observer's performances on unseen images were statistically equivalent with a margin of difference ( Δ AUC ) of 0.0426 at p < 0.05 , using 288 training images. A limited implementation of an SLDO had a substantially higher AUC (0.99) compared to the DeepAMO and human observer. Conclusion: The results show that the DeepAMO has the potential to reproduce the absolute performance, and not just the relative ranking of human observers on a clinically realistic defect detection task, and that building conceptual components of the human reading process into deep learning-based models can allow training of these models in settings where limited training images are available.

Overcoming Barriers to Radiopharmaceutical Therapy (RPT): An Overview From the NRG-NCI Working Group on Dosimetry of Radiopharmaceutical Therapy.

Radiopharmaceutical therapy (RPT) continues to demonstrate tremendous potential in improving the therapeutic gains in radiation therapy by specifically delivering radiation to tumors that can be well assessed in terms of dosimetry and imaging. Dosimetry in external beam radiation therapy is standard practice. This is not the case, however, in RPT. This NRG (acronym formed from the first letter of the 3 original groups: National Surgical Adjuvant Breast and Bowel Project, the Radiation Therapy Oncology Group, and the Gynecologic Oncology Group)-National Cancer Institute Working Group review describes some of the challenges to improving RPT. The main priorities for advancing the field include (1) developing and adopting best practice guidelines for incorporating patient-specific dosimetry for RPT that can be used at both large clinics with substantial resources and more modest clinics that have limited resources, (2) establishing and improving strategies for introducing new radiopharmaceuticals for clinical investigation, (3) developing approaches to address the radiophobia that is associated with the administration of radioactivity for cancer therapy, and (4) solving the financial and logistical issues of expertise and training in the developing field of RPT.

Dosimetric considerations of 99mTc-MDP uptake within the epiphyseal plates of the long bones of pediatric patients.

Skeletal scintigraphy is most performed in pediatric patients using the radiopharmaceutical 99mTc labelled methylene diphosphonate (99mTc-MDP). Reference biokinetic models for 99mTc-MDP indicate 50% of the administered activity is uniformly localized to the interior bone surfaces (trabecular and cortical regions), yet imaging data clearly show some preferential uptake to the epiphyseal growth plates of the long bones. To explore the dosimetric consequences of these regional activity concentrations, we have modified mesh-type computational phantoms of the International Commission on Radiological Protection (ICRP) reference pediatric series to explicitly include geometric models of the epiphyseal growth plates (2 mm in thickness) within the left/right, distal/proximal ends of the humeri, radii, ulnae, femora, tibia, and fibulae. Bone mineral activity from the ICRP Publication 128 biokinetic model for 99mTc-MDP (ICRP 2015) was then partitioned to the growth plates at values of 0.5%, 4.4%, 8.3%, 12.2%, 16.1%, and 20%. Radiation transport simulations were performed to compute 99mTc S-values and organ dose coefficients to the soft tissues and to bone site-specific regions of spongiosa. As the percentage of bone activity assigned to the growth plates was increased (from 0.5% to 20%), absorbed doses to the soft tissue organs, active bone marrow, bone endosteum (BE), as well as the detriment-weighted dose, were shown to decrease from their nominal values (no substantial growth plate activity), while epiphyseal plate self-doses increased. In the 15 year old male phantom, moving from 0.5% to 20% relative bone activity within the epiphyseal plates resulted in a 15% reduction in active marrow (AM) and BE dose, a 10% reduction in mean soft tissue and detriment-weighted dose, and a 6.3-fold increase in epiphyseal plate self-dose. In the newborn female phantom, we observed a 18% decrease in AM and BE dose, a 10% decrease in mean soft tissue dose, a 15% decrease in detriment-weighted dose, and 12.8-fold increase in epiphyseal plate self-dose. Increases (to 3 mm) and decreases (to 1 mm) in the assumed growth plate thickness of our models were shown to impact only the growth plate self-dose. Future work in differential quantification of 99mTc-MDP activity-growth plates versus other bone surfaces-is required to provide clinically realistic data on activity partitioning as a function of patient age, and perhaps skeletal site. The phantom series presented here may be used to develop more optimized age-related guidance on 99mTc-MDP administered activities to children.

Body morphometry appropriate computational phantoms for dose and risk optimization in pediatric renal imaging with Tc-99m DMSA and Tc-99m MAG3.

Current guidelines for administered activity (AA) in pediatric nuclear medicine imaging studies are based on a 2016 harmonization of the 2010 North American Consensus guidelines and the 2007 European Association of Nuclear Medicine pediatric dosage card. These guidelines assign AA scaled to patient body mass, with further constraints on maximum and minimum values of radiopharmaceutical activity. These guidelines, however, are not formulated based upon a rigor-ous evaluation of diagnostic image quality. In a recent study of the renal cortex imaging agent 99mTc-DMSA (Li Y et al 2019), body mass-based dosing guidelines were shown to not give the same level of image quality for patients of differing body mass. Their data suggest that patient girth at the level of the kidneys may be a better morphometric parameter to consider when selecting AA for renal nuclear medicine imaging. The objective of the present work was thus to develop a dedicated series of computational phantoms to support image quality and organ dose studies in pediatric renal imaging using 99mTc-DMSA or 99mTc-MAG3. The final library consists of 50 male and female phantoms of ages 0 to 15 years, with percentile variations (5th to 95th) in waist circumference (WC) at each age. For each phantom, nominal values of kidney volume, length, and depth were incorporated into the phantom design. Organ absorbed doses, detriment-weighted doses, and stochastic risks were assessed using ICRP reference biokinetic models for both agents. In Monte Carlo radiation transport simulations, organ doses for these agents yielded detriment-weighted dose coefficients (mSv/MBq) that were in general larger than current ICRP values of the effective dose coefficients (age and WC-averaged ratios of eDW/e were 1.40 for the male phantoms and 1.49 for the female phantoms). Values of risk index (ratio of radiation-induced to natural background cancer incidence risk x 100) varied between 0.062 (newborns) to 0.108 (15-year-olds) for 99mTc-DMSA and between 0.026 (newborns) to 0.122 (15-year-olds) for 99mTc-MAG3. Using tallies of photon exit fluence as a rough surrogate for uniform image quality, our study demonstrated that through body region-of-interest optimization of AA, there is the potential for further dose and risk reductions of between factors of 1.5 to 3.0 beyond simple weight-based dosing guidance.

Three-dimensional regions-of-interest-based intra-operative four-dimensional soft tissue perfusion imaging using a standard x-ray system with no gantry rotation: A simulation study for a proof of concept.

PURPOSE: Many interventional procedures aim at changing soft tissue perfusion or blood flow. One problem at present is that soft tissue perfusion and its changes cannot be assessed in an interventional suite because cone-beam computed tomography is too slow (it takes 4-10 s per volume scan). In order to address the problem, we propose a novel method called IPEN for Intra-operative four-dimensional soft tissue PErfusion using a standard x-ray system with No gantry rotation. METHODS: IPEN uses two input datasets: (a) the contours and locations of three-dimensional regions-of-interest (ROIs) such as arteries and sub-sections of cancerous lesions, and (b) a series of x-ray projection data obtained from an intra-arterial contrast injection to contrast enhancement to wash-out. IPEN then estimates a time-enhancement curve (TEC) for each ROI directly from projections without reconstructing cross-sectional images by maximizing the agreement between synthesized and measured projections with a temporal roughness penalty. When path lengths through ROIs are known for each x-ray beam, the ROI-specific enhancement can be accurately estimated from projections. Computer simulations are performed to assess the performance of the IPEN algorithm. Intra-arterial contrast-enhanced liver scans over 25 s were simulated using XCAT phantom version 2.0 with heterogeneous tissue textures and cancerous lesions. The following four sub-studies were performed: (a) The accuracy of the estimated TECs with overlapped lesions was evaluated at various noise (dose) levels with either homogeneous or heterogeneous lesion enhancement patterns; (b) the accuracy of IPEN with inaccurate ROI contours was assessed; (c) we investigated how overlapping ROIs and noise in projections affected the accuracy of the IPEN algorithm; and (d) the accuracy of the perfusion indices was assessed. RESULTS: The TECs estimated by IPEN were sufficiently accurate at a reference dose level with the root-mean-square deviation (RMSD) of 0.0027 ± 0.0001 cm-1 or 13 ± 1 Hounsfield unit (mean ± standard deviation) for the homogeneous lesion enhancement and 0.0032 ± 0.0005 cm-1 for the heterogeneous enhancement (N = 20 each). The accuracy was degraded with decreasing doses: The RMSD with homogeneous enhancement was 0.0220 ± 0.0003 cm-1 for 20% of the reference dose level. Performing 3 × 3 pixel averaging on projection data improved the RMSDs to 0.0051 ± 0.0002 cm-1 for 20% dose. When the ROI contours were inaccurate, smaller ROI contours resulted in positive biases in TECs, whereas larger ROI contours produced negative biases. The bias remained small, within ± 0.0070 cm-1 , when the Sorenson-Dice coefficients (SDCs) were larger than 0.81. The RMSD of the TEC estimation was strongly associated with the condition of the problem, which can be empirically quantified using the condition number of a matrix A z that maps a vector of ROI enhancement values z to projection data and a weighted variance of projection data: a linear correlation coefficient (R) was 0.794 (P < 0.001). The perfusion index values computed from the estimated TECs agreed well with the true values (R ≥ 0.985, P < 0.0001). CONCLUSION: The IPEN algorithm can estimate ROI-specific TECs with high accuracy especially when 3 × 3 pixel averaging is applied, even when lesion enhancement is heterogeneous, or ROI contours are inaccurate but the SDC is at least 0.81.

Generating anthropomorphic phantoms using fully unsupervised deformable image registration with convolutional neural networks.

PURPOSE: Computerized phantoms have been widely used in nuclear medicine imaging for imaging system optimization and validation. Although the existing computerized phantoms can model anatomical variations through organ and phantom scaling, they do not provide a way to fully reproduce the anatomical variations and details seen in humans. In this work, we present a novel registration-based method for creating highly anatomically detailed computerized phantoms. We experimentally show substantially improved image similarity of the generated phantom to a patient image. METHODS: We propose a deep-learning-based unsupervised registration method to generate a highly anatomically detailed computerized phantom by warping an XCAT phantom to a patient computed tomography (CT) scan. We implemented and evaluated the proposed method using the NURBS-based XCAT phantom and a publicly available low-dose CT dataset from TCIA. A rigorous tradeoff analysis between image similarity and deformation regularization was conducted to select the loss function and regularization term for the proposed method. A novel SSIM-based unsupervised objective function was proposed. Finally, ablation studies were conducted to evaluate the performance of the proposed method (using the optimal regularization and loss function) and the current state-of-the-art unsupervised registration methods. RESULTS: The proposed method outperformed the state-of-the-art registration methods, such as SyN and VoxelMorph, by more than 8%, measured by the SSIM and less than 30%, by the MSE. The phantom generated by the proposed method was highly detailed and was almost identical in appearance to a patient image. CONCLUSIONS: A deep-learning-based unsupervised registration method was developed to create anthropomorphic phantoms with anatomies labels that can be used as the basis for modeling organ properties. Experimental results demonstrate the effectiveness of the proposed method. The resulting anthropomorphic phantom is highly realistic. Combined with realistic simulations of the image formation process, the generated phantoms could serve in many applications of medical imaging research.

Fisher information analysis of list-mode SPECT emission data for joint estimation of activity and attenuation distribution.

The potential to perform attenuation and scatter compensation (ASC) in single-photon emission computed tomography (SPECT) imaging without a separate transmission scan is highly significant. In this context, attenuation in SPECT is primarily due to Compton scattering, where the probability of Compton scatter is proportional to the attenuation coefficient of the tissue and the energy of the scattered photon and the scattering angle are related. Based on this premise, we investigated whether the SPECT scattered-photon data acquired in list-mode (LM) format and including the energy information can be used to estimate the attenuation map. For this purpose, we propose a Fisher-information-based method that yields the Cramer-Rao bound (CRB) for the task of jointly estimating the activity and attenuation distribution using only the SPECT emission data. In the process, a path-based formalism to process the LM SPECT emission data, including the scattered-photon data, is proposed. The Fisher information method was implemented on NVIDIA graphics processing units (GPU) for acceleration. The method was applied to analyze the information content of SPECT LM emission data, which contains up to first-order scattered events, in a simulated SPECT system with parameters modeling a clinical system using realistic computational studies with 2-D digital synthetic and anthropomorphic phantoms. The method was also applied to LM data containing up to second-order scatter for a synthetic phantom. Experiments with anthropomorphic phantoms simulated myocardial perfusion and dopamine transporter (DaT)-Scan SPECT studies. The results show that the CRB obtained for the attenuation and activity coefficients was typically much lower than the true value of these coefficients. An increase in the number of detected photons yielded lower CRB for both the attenuation and activity coefficients. Further, we observed that systems with better energy resolution yielded a lower CRB for the attenuation coefficient. Overall, the results provide evidence that LM SPECT emission data, including the scattered photons, contains information to jointly estimate the activity and attenuation coefficients.

Verification of a method to detect glass microspheres via micro-CT.

INTRODUCTION: The ability to determine the microscopic distribution of glass microspheres in 90 Y radioembolization has important applications in post-treatment microdosimetry and cluster analysis. Current methods are time-intensive and labor-intensive and thus are typically only applied to small samples. MATERIALS AND METHODS: A high-resolution micro-CT image with a voxel size of 8.74 µm was acquired of phantoms containing ~25 µm-diameter glass microspheres embedded in tissue-equivalent materials that were optically transparent, which allowed true microsphere locations to be determined using transmission light microscopy. A 3-stage algorithm was developed to estimate the number and locations of microspheres in tissue regions. The stages are thresholding the CT image and discarding regions with insufficient voxels, estimating the number of microspheres in each region using the values of the detected and neighboring region voxels and estimating locations for each microsphere using the outputs of the previous two stages. Two different methods for estimating the number of microspheres in each region were derived, as were five methods for localizing microspheres. Metrics for each stage were computed, and the mean absolute error (MAE) between the dose to 72 µm voxels of the true and estimated dose maps created from the microsphere locations was used as the figure of merit for overall algorithm performance. Microsphere locations identified in the optical micrograph were used as the gold standard for the metrics of all stages. The method's utility was then demonstrated using a specimen from a human neuroendocrine tumor (NET) treated with glass 90 Y microspheres. RESULTS: The stage detecting regions containing microspheres found 100% of microspheres inside regions. The number of incorrectly detected regions without microspheres was 1.5% of the total number of regions. In stage 2, with these parameters, nearly 94% of the actual number of spheres in each region was correctly counted, and only 5% of the estimated sphere quantities in each region were false positives. The MAE between the true dose maps and dose maps estimated using the full algorithm with optimal parameter and method choices was 4.2%. A total of 5,713 glass microspheres were identified as being distributed heterogeneously in the NET specimen with a maximum tumor dose of >2500 Gy and 46% of the specimen receiving <20 Gy. CONCLUSIONS: This work developed and evaluated a method to detect and estimate the three-dimensional locations of glass microspheres in whole tissue samples that require less manual effort than traditional methods. This method could be used to gain important insights into the heterogeneity of microsphere distributions that would be useful for improving radioembolization treatment planning.

Current pediatric administered activity guidelines for 99m Tc-DMSA SPECT based on patient weight do not provide the same task-based image quality.

PURPOSE: In the current clinical practice, administered activity (AA) for pediatric molecular imaging is often based on the North American expert consensus guidelines or the European Association of Nuclear Medicine dosage card, both of which were developed based on the best clinical practice. These guidelines were not formulated using a rigorous evaluation of diagnostic image quality (IQ) relative to AA. In the guidelines, AA is determined by a weight-based scaling of the adult AA, along with minimum and maximum AA constraints. In this study, we use task-based IQ assessment methods to rigorously evaluate the efficacy of weight-based scaling in equalizing IQ using a population of pediatric patients of different ages and body weights. METHODS: A previously developed projection image database was used. We measured task-based IQ, with respect to the detection of a renal functional defect at six different AA levels (AA relative to the AA obtained from the guidelines). IQ was assessed using an anthropomorphic model observer. Receiver-operating characteristics (ROC) analysis was applied; the area under the ROC curve (AUC) served as a figure-of-merit for task performance. In addition, we investigated patient girth (circumference) as a potential improved predictor of the IQ. RESULTS: The data demonstrate a monotonic and modestly saturating increase in AUC with increasing AA, indicating that defect detectability was limited by quantum noise and the effects of object variability were modest over the range of AA levels studied. The AA for a given value of the AUC increased with increasing age. The AUC vs AA plots for all the patient ages indicate that, for the current guidelines, the newborn and 10- and 15-yr phantoms had similar IQ for the same AA suggested by the North American expert consensus guidelines, but the 5- and 1-yr phantoms had lower IQ. The results also showed that girth has a stronger correlation with the needed AA to provide a constant AUC for 99m Tc-DMSA renal SPECT. CONCLUSIONS: The results suggest that (a) weight-based scaling is not sufficient to equalize task-based IQ for patients of different weights in pediatric 99m Tc-DMSA renal SPECT; and (b) patient girth should be considered instead of weight in developing new administration guidelines for pediatric patients.

Development of a Customizable Hepatic Arterial Tree and Particle Transport Model for Use in Treatment Planning.

Optimal treatment planning for radioembolization of hepatic cancers produces sufficient dose to tumors for control and dose to normal liver parenchyma that is below the threshold for toxicity. The non-uniform distribution of particles in liver microanatomy complicates the planning process as different functional regions receive different doses. Having realistic and patient-specific models of the arterial tree and microsphere trapping would be useful for developing more optimal treatment plans. We propose a macrocell-based growth method to generate models of the hepatic arterial tree from the proper hepatic artery to the terminal arterioles supplying the capillaries in the parenchyma. We show how these trees can be adapted to match patient values of pressure, flow, and vessel diameters while still conforming to laws controlling vessel bifurcation, changes in pressure, and blood flow. We also introduce a method to model particle transport within the tree that accounts for vessel and particle diameter distributions and show the non-uniform microsphere deposition pattern that results. Potential applications include investigating dose heterogeneity and microsphere deposition patterns.

A projection image database to investigate factors affecting image quality in weight-based dosing: application to pediatric renal SPECT.

Balancing the tradeoff between radiation dose, acquisition duration and diagnostic image quality is essential for medical imaging modalities involving ionizing radiation. Lower administered activities to the patient can reduce absorbed dose, but can result in reduced diagnostic image quality or require longer acquisition durations. In pediatric nuclear medicine, it is desirable to use the lowest amount of administered radiopharmaceutical activity and the shortest acquisition duration that gives sufficient image quality for clinical diagnosis. However, diagnostic image quality is a complex function of patient factors including body morphometry. In this study, we present a digital population of 90 computational anatomic phantoms that model realistic variations in body morphometry and internal anatomy. These phantoms were used to generate a large database of projection images modeling pediatric SPECT imaging using a 99mTc-DMSA tracer. We used an analytic projection code that models attenuation, spatially varying collimator-detector response, and object-dependent scatter to generate the projections. The projections for each organ were generated separately and can be subsequently scaled by parameters extracted from a pharmacokinetics model to simulate realistic tracer biodistribution, including variations in uptake, inside each relevant organ or tissue structure for a given tracer. Noise-free projection images can be obtained by summing these individual organ projections and scaling by the system sensitivity and acquisition duration. We applied this database in the context of 99mTc-DMSA renal SPECT, the most common nuclear medicine imaging procedure in pediatric patients. Organ uptake fractions based on literature values and patient studies were used. Patient SPECT images were used to verify that the sum of counts in the simulated projection images was clinically realistic. For each phantom, 384 uptake realizations, modeling random variations in the uptakes of organs of interest, were generated, producing 34 560 noise-free projection datasets (384 uptake realizations times 90 phantoms). Noisy images modeling various count levels (corresponding to different products of acquisition duration and administered activity) were generated by appropriately scaling these images and simulating Poisson noise. Acquisition duration was fixed; six count levels were simulated corresponding to projection images acquired using 25%, 50%, 75%, 100%, 125%, and 150% of the original weight-based administrated activity as computed using the North American Guidelines (Gelfand et al 2011 J. Nucl. Med. 52 318-22). Combined, a total number of 207 360 noisy projection images were generated, creating a realistic projection database for use in renal pediatric SPECT imaging research. The phantoms and projection datasets were used to calculate three surrogate indices for factors affecting image quality: renal count density, average radius of rotation, and scatter-to-primary ratio. Differences in these indices were seen across the phantoms for dosing based on current guidelines, and especially for the phantom modeling the newborn. We also performed an image quality study using an anthropomorphic model observer that demonstrates that the weight-based dose scaling does not equalize image quality as measured by the area under the receiver-operating characteristics curve. These studies suggest that a dosing procedure beyond weight-based scaling of administered activities is required to equalize image quality in pediatric renal SPECT.