Chronic thromboembolic pulmonary hypertension and balloon pulmonary angioplasty - Where are we in 2024?

Pulmonary endarterectomy (PEA) is the first-line treatment for patients with chronic thromboembolic pulmonary hypertension (CTEPH). However, some patients with CTEPH are considered inoperable, and in the last decade, balloon pulmonary angioplasty (BPA) has emerged as a viable therapeutic option for these patients with prohibitive surgical risk or recurrent pulmonary hypertension following PEA. Numerous international centers have increased their procedural volume of BPA and have reported improvements in pulmonary hemodynamics, patient functional class and right ventricular function. Randomized controlled trials have also demonstrated similar findings. Recent refinements in procedural technique, increased operator experience and advancements in procedural technology have facilitated marked reduction in the risk of complications following BPA. Current guidelines recommend BPA for patients with inoperable CTEPH and persistent pulmonary hypertension following PEA. The pulmonary arterial endothelium plays a vital role in the pathophysiologic development and progression of CTEPH.

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.

The Hypoxic Tricuspid: A Flail Tricuspid Valve With Patent Foramen Ovale.

Etiologies of tricuspid regurgitation are often explored in patients with symptoms of right-sided heart failure. Blunt chest trauma is the major cause of traumatic tricuspid valve regurgitation (TTVR), a secondary type of tricuspid regurgitation. It is a rare condition; however, it may lead to severe consequences if not treated in a timely manner. TTVR should be considered in a patient presenting with chest trauma. In this case, we report a case of a young male who presented after a motor vehicle accident with secondary tricuspid valve regurgitation due to blunt chest trauma as well as a patent foramen ovale.

Dyskeratosis Congenita: A Case Report of a Patient With Coronary Artery Disease.

Clinical evidence demonstrates that patients with telomere biology disorders, such as dyskeratosis congenita, are more prone to coronary artery disease. We present the case of a 43-year-old female diagnosed with dyskeratosis congenita with critical cardiovascular disease. She underwent coronary artery bypass graft (CABG) with improvement of her cardiac function. Although this is a rare genetic disease, further studies are warranted to investigate the underlying pathophysiology of cardiovascular disease in patients with dyskeratosis congenita.

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.

Overview of the First NRG Oncology-National Cancer Institute Workshop on Dosimetry of Systemic Radiopharmaceutical Therapy.

In 2018, the National Cancer Institute and NRG Oncology partnered for the first time to host a joint workshop on systemic radiopharmaceutical therapy (RPT) to specifically address dosimetry issues and strategies for future clinical trials. The workshop focused on current dosimetric approaches for clinical trials, strategies under development that would optimize dose reporting, and future desired or optimized approaches for novel emerging radionuclides and carriers in development. In this article, we review the main approaches that are applied clinically to calculate the absorbed dose. These include absorbed doses calculated over a variety of spatial scales, including whole body, organ, suborgan, and voxel, the last 3 of which are achievable within the MIRD schema (S value) and can be calculated with analytic methods or Monte Carlo methods, the latter in most circumstances. This article will also contrast currently available methods and tools with those used in the past, to propose a pathway whereby dosimetry helps the field by optimizing the biologic effect of the treatment and trial design in the drug approval process to reduce financial and logistical costs. We also briefly discuss the dosimetric equivalent of biomarkers to help bring a precision medicine approach to RPT implementation when merited by evidence collected during early-phase trial investigations. Advances in the methodology and related tools have made dosimetry the optimum biomarker for RPT.

Specific absorbed fractions and radionuclide S-values for tumors of varying size and composition.

Accurate estimates of tumor absorbed dose are essential for the evaluation of treatment efficacy in radiopharmaceutical cancer therapy. Although tumor dosimetry via the MIRD schema has been previously investigated, prior studies have been limited to the consideration of soft-tissue tumors. In the present study, specific absorbed fractions (SAFs) for monoenergetic photons, electrons, and alpha particles in tumors of varying compositions were computed using Monte Carlo simulations in MCNPX after which self-irradiation S-values for 22 radionuclides (along with 14 additional alpha-emitter progeny) were generated for tumors of both varying size and tissue composition. The tumors were modeled as spheres with radii ranging from 0.10 cm to 6.0 cm and with compositions varying from 100% soft tissue (ST) to 100% mineral bone (MB). The energies of the photons and electrons were varied on a logarithm energy grid from 10 keV to 10 MeV. The energies of alpha particles were varied along a linear energy grid from 0.5 MeV to 12 MeV. In all cases, a homogenous activity distribution was assumed throughout the tumor volume. Furthermore, to assess the effect of tumor shape, several ellipsoidal tumors of different compositions were modeled and absorbed fractions were computed for monoenergetic electrons and photons. S-values were then generated using detailed decay data from the 2008 MIRD Monograph on Radionuclide Data and Decay Schemes. Our study results demonstrate that a soft-tissue model yields relative errors of 25% and 71% in the absorbed fraction assigned to uniform sources of 1.5 MeV electrons and 100 keV photons, respectively, localized within a 1 cm diameter tumor of MB. The data further show that absorbed fractions for moderate ellipsoids can be well approximated by a spherical shape of equal mass within a relative error of < 8%. S-values for 22 radionuclides (and their daughter progeny) were computed with results demonstrating how relative errors in SAFs could propagate to relative errors in tumor dose estimates as high as 86%. A comprehensive data set of radionuclide S-values by tumor size and tissue composition is provided for application of the MIRD schema for tumor dosimetry in radiopharmaceutical therapy.

Evaluation of different respiratory gating schemes for cardiac SPECT.

BACKGROUND: Respiratory gating reduces motion blurring in cardiac SPECT. Here we aim to evaluate the performance of three respiratory gating strategies using a population of digital phantoms with known truth and clinical data. METHODS: We analytically simulated 60 projections for 10 XCAT phantoms with 99mTc-sestamibi distributions using three gating schemes: equal amplitude gating (AG), equal count gating (CG), and equal time gating (TG). Clinical list-mode data for 10 patients who underwent 99mTc-sestamibi scans were also processed using the 3 gating schemes. Reconstructed images in each gate were registered to a reference gate, averaged and reoriented to generate the polar plots. For simulations, image noise, relative difference (RD) of averaged count for each of the 17 segment, and relative defect size difference (RSD) were analyzed. For clinical data, image intensity profile and FWHM were measured across the left ventricle wall. RESULTS: For simulations, AG and CG methods showed significantly lower RD and RSD compared to TG, while noise variation was more non-uniform through different gates for AG. In the clinical study, AG and CG had smaller FWHM than TG. CONCLUSIONS: AG and CG methods show better performance for motion reduction and are recommended for clinical respiratory gating SPECT implementation.

Interpolated CT for attenuation correction on respiratory gating cardiac SPECT/CT - A simulation study.

PURPOSE: Respiratory gated four-dimensional (4D) single photon emission computed tomography (SPECT) with phase-matched CT reduces respiratory blurring and attenuation correction (AC) artifacts in cardiac SPECT. This study aims to develop and investigate the effectiveness of an interpolated CT (ICT) method for improved cardiac SPECT AC using simulations. METHODS: We used the 4D XCAT phantom to simulate a population of ten patients varied in gender, anatomy, 99m Tc-sestamibi distribution, respiratory patterns, and disease states. Simulated 120 SPECT projection data were rebinned into six equal count gates. Activity and attenuation maps in each gate were averaged as gated SPECT and CT (GCT). Three helical CTs were simulated at end-inspiration (HCT-IN), end-expiration (HCT-EX), and mid-respiration (HCT-MID). The ICTs were obtained from HCT-EX and HCT-IN using the motion vector field generated between them from affine plus b-spline registration. Projections were reconstructed by OS-EM method, using GCT, ICT, and three HCTs for AC. Reconstructed images of each gate were registered to end-expiration and averaged to generate the polar plots. Relative difference for each segment and relative defect size were computed using images of GCT AC as reference. RESULTS: The average of maximum relative difference through ten phantoms was 7.93 ± 4.71%, 2.50 ± 0.98%, 3.58 ± 0.74%, and 2.14 ± 0.56% for noisy HCT-IN, HCT-MID, HCT-EX, and ICT AC data, respectively. The ICT showed closest defect size to GCT while the differences from HCTs can be over 40%. CONCLUSION: We conclude that the performance of ICT is similar to GCT. It improves the image quality and quantitative accuracy for respiratory-gated cardiac SPECT as compared to conventional HCT, while it can potentially further reduce the radiation dose of GCT.

A comparison of resampling schemes for estimating model observer performance with small ensembles.

In objective assessment of image quality, an ensemble of images is used to compute the 1st and 2nd order statistics of the data. Often, only a finite number of images is available, leading to the issue of statistical variability in numerical observer performance. Resampling-based strategies can help overcome this issue. In this paper, we compared different combinations of resampling schemes (the leave-one-out (LOO) and the half-train/half-test (HT/HT)) and model observers (the conventional channelized Hotelling observer (CHO), channelized linear discriminant (CLD) and channelized quadratic discriminant). Observer performance was quantified by the area under the ROC curve (AUC). For a binary classification task and for each observer, the AUC value for an ensemble size of 2000 samples per class served as a gold standard for that observer. Results indicated that each observer yielded a different performance depending on the ensemble size and the resampling scheme. For a small ensemble size, the combination [CHO, HT/HT] had more accurate rankings than the combination [CHO, LOO]. Using the LOO scheme, the CLD and CHO had similar performance for large ensembles. However, the CLD outperformed the CHO and gave more accurate rankings for smaller ensembles. As the ensemble size decreased, the performance of the [CHO, LOO] combination seriously deteriorated as opposed to the [CLD, LOO] combination. Thus, it might be desirable to use the CLD with the LOO scheme when smaller ensemble size is available.

Use of Sub-Ensembles and Multi-Template Observers to Evaluate Detection Task Performance for Data That are Not Multivariate Normal.

The Hotelling Observer (HO) is widely used to evaluate image quality in medical imaging. However, applying it to data that are not multivariate-normally (MVN) distributed is not optimal. In this paper, we apply two multi-template linear observer strategies to handle such data. First, the entire data ensemble is divided into sub-ensembles that are exactly or approximately MVN and homoscedastic. Next, a different linear observer template is estimated for and applied to each sub-ensemble. The first multi-template strategy, adapted from previous work, applies the HO to each sub-ensemble, calculates the area under the receiver operating characteristics curve (AUC) for each sub-ensemble, and averages the AUCs from all the sub-ensembles. The second strategy applies the Linear Discriminant (LD) to estimate test statistics for each sub-ensemble and calculates a single global AUC using the pooled test statistics from all the sub-ensembles. We show that this second strategy produces the maximum AUC when only shifting of the HO test statistics is allowed. We compared these strategies to the use of a single HO template for the entire data ensemble by applying them to the non-MVN data obtained from reconstructed images of a realistic simulated population of myocardial perfusion SPECT studies with the goal of optimizing the reconstruction parameters. Of the strategies investigated, the multi-template LD strategy yielded the highest AUC for any given set of reconstruction parameters. The optimal reconstruction parameters obtained by the two multi-template strategies were comparable and produced higher AUCs for each sub-ensemble than the single-template HO strategy.

Factors affecting the normality of channel outputs of channelized model observers: an investigation using realistic myocardial perfusion SPECT images.

The channelized Hotelling observer (CHO) uses the first- and second-order statistics of channel outputs under both hypotheses to compute test statistics used in binary classification tasks. If these input data deviate from a multivariate normal (MVN) distribution, the classification performance will be suboptimal compared to an ideal observer operating on the same channel outputs. We conducted a comprehensive investigation to rigorously study the validity of the MVN assumption under various kinds of background and signal variability in a realistic population of phantoms. The study was performed in the context of myocardial perfusion SPECT imaging; anatomical, uptake (intensity), and signal variability were simulated. Quantitative measures and graphical approaches applied to the outputs of each channel were used to investigate the amount and type of deviation from normality. For some types of background and signal variations, the channel outputs, under both hypotheses, were non-normal (i.e., skewed or multimodal). This indicates that, for realistic medical images in cases where there is signal or background variability, the normality of the channel outputs should be evaluated before applying a CHO. Finally, the different degrees of departure from normality of the various channels are explained in terms of violations of the central limit theorem.

Collimator optimization in myocardial perfusion SPECT using the ideal observer and realistic background variability for lesion detection and joint detection and localization tasks.

In SPECT imaging, collimators are a major factor limiting image quality and largely determine the noise and resolution of SPECT images. In this paper, we seek the collimator with the optimal tradeoff between image noise and resolution with respect to performance on two tasks related to myocardial perfusion SPECT: perfusion defect detection and joint detection and localization. We used the Ideal Observer (IO) operating on realistic background-known-statistically (BKS) and signal-known-exactly (SKE) data. The areas under the receiver operating characteristic (ROC) and localization ROC (LROC) curves (AUCd, AUCd+l), respectively, were used as the figures of merit for both tasks. We used a previously developed population of 54 phantoms based on the eXtended Cardiac Torso Phantom (XCAT) that included variations in gender, body size, heart size and subcutaneous adipose tissue level. For each phantom, organ uptakes were varied randomly based on distributions observed in patient data. We simulated perfusion defects at six different locations with extents and severities of 10% and 25%, respectively, which represented challenging but clinically relevant defects. The extent and severity are, respectively, the perfusion defect's fraction of the myocardial volume and reduction of uptake relative to the normal myocardium. Projection data were generated using an analytical projector that modeled attenuation, scatter, and collimator-detector response effects, a 9% energy resolution at 140 keV, and a 4 mm full-width at half maximum (FWHM) intrinsic spatial resolution. We investigated a family of eight parallel-hole collimators that spanned a large range of sensitivity-resolution tradeoffs. For each collimator and defect location, the IO test statistics were computed using a Markov Chain Monte Carlo (MCMC) method for an ensemble of 540 pairs of defect-present and -absent images that included the aforementioned anatomical and uptake variability. Sets of test statistics were computed for both tasks and analyzed using ROC and LROC analysis methodologies. The results of this study suggest that collimators with somewhat poorer resolution and higher sensitivity than those of a typical low-energy high-resolution (LEHR) collimator were optimal for both defect detection and joint detection and localization tasks in myocardial perfusion SPECT for the range of defect sizes investigated. This study also indicates that optimizing instrumentation for a detection task may provide near-optimal performance on the more challenging detection-localization task.

Collimator optimization and collimator-detector response compensation in myocardial perfusion SPECT using the ideal observer with and without model mismatch and an anthropomorphic model observer.

The collimator is the primary factor that determines the spatial resolution and noise tradeoff in myocardial perfusion SPECT images. In this paper, the goal was to find the collimator that optimizes the image quality in terms of a perfusion defect detection task. Since the optimal collimator could depend on the level of approximation of the collimator-detector response (CDR) compensation modeled in reconstruction, we performed this optimization for the cases of modeling the full CDR (including geometric, septal penetration and septal scatter responses), the geometric CDR, or no model of the CDR. We evaluated the performance on the detection task using three model observers. Two observers operated on data in the projection domain: the Ideal Observer (IO) and IO with Model-Mismatch (IO-MM). The third observer was an anthropomorphic Channelized Hotelling Observer (CHO), which operated on reconstructed images. The projection-domain observers have the advantage that they are computationally less intensive. The IO has perfect knowledge of the image formation process, i.e. it has a perfect model of the CDR. The IO-MM takes into account the mismatch between the true (complete and accurate) model and an approximate model, e.g. one that might be used in reconstruction. We evaluated the utility of these projection domain observers in optimizing instrumentation parameters. We investigated a family of 8 parallel-hole collimators, spanning a wide range of resolution and sensitivity tradeoffs, using a population of simulated projection (for the IO and IO-MM) and reconstructed (for the CHO) images that included background variability. We simulated anterolateral and inferior perfusion defects with variable extents and severities. The area under the ROC curve was estimated from the IO, IO-MM, and CHO test statistics and served as the figure-of-merit. The optimal collimator for the IO had a resolution of 9-11 mm FWHM at 10 cm, which is poorer resolution than typical collimators used for MPS. When the IO-MM and CHO used a geometric or no model of the CDR, the optimal collimator shifted toward higher resolution than that obtained using the IO and the CHO with full CDR modeling. With the optimal collimator, the IO-MM and CHO using geometric modeling gave similar performance to full CDR modeling. Collimators with poorer resolution were optimal when CDR modeling was used. The agreement of rankings between the IO-MM and CHO confirmed that the IO-MM is useful for optimization tasks when model mismatch is present due to its substantially reduced computational burden compared to the CHO.

Optimization of energy window and evaluation of scatter compensation methods in myocardial perfusion SPECT using the ideal observer with and without model mismatch and an anthropomorphic model observer.

We used the ideal observer (IO) and IO with model mismatch (IO-MM) applied in the projection domain and an anthropomorphic channelized Hotelling observer (CHO) applied to reconstructed images to optimize the acquisition energy window width and to evaluate various scatter compensation methods in the context of a myocardial perfusion single-photon emission computed tomography (SPECT) defect detection task. The IO has perfect knowledge of the image formation process and thus reflects the performance with perfect compensation for image-degrading factors. Thus, using the IO to optimize imaging systems could lead to suboptimal parameters compared with those optimized for humans interpreting SPECT images reconstructed with imperfect or no compensation. The IO-MM allows incorporating imperfect system models into the IO optimization process. We found that with near-perfect scatter compensation, the optimal energy window for the IO and CHO was similar; in its absence, the IO-MM gave a better prediction of the optimal energy window for the CHO using different scatter compensation methods. These data suggest that the IO-MM may be useful for projectiondomain optimization when MM is significant and that the IO is useful when followed by reconstruction with good models of the image formation process.

Optimization and comparison of simultaneous and separate acquisition protocols for dual isotope myocardial perfusion SPECT.

Dual-isotope simultaneous-acquisition (DISA) rest-stress myocardial perfusion SPECT (MPS) protocols offer a number of advantages over separate acquisition. However, crosstalk contamination due to scatter in the patient and interactions in the collimator degrade image quality. Compensation can reduce the effects of crosstalk, but does not entirely eliminate image degradations. Optimizing acquisition parameters could further reduce the impact of crosstalk. In this paper we investigate the optimization of the rest Tl-201 energy window width and relative injected activities using the ideal observer (IO), a realistic digital phantom population and Monte Carlo (MC) simulated Tc-99m and Tl-201 projections as a means to improve image quality. We compared performance on a perfusion defect detection task for Tl-201 acquisition energy window widths varying from 4 to 40 keV centered at 72 keV for a camera with a 9% energy resolution. We also investigated 7 different relative injected activities, defined as the ratio of Tc-99m and Tl-201 activities, while keeping the total effective dose constant at 13.5 mSv. For each energy window and relative injected activity, we computed the IO test statistics using a Markov chain Monte Carlo (MCMC) method for an ensemble of 1,620 triplets of fixed and reversible defect-present, and defect-absent noisy images modeling realistic background variations. The volume under the 3-class receiver operating characteristic (ROC) surface (VUS) was estimated and served as the figure of merit. For simultaneous acquisition, the IO suggested that relative Tc-to-Tl injected activity ratios of 2.6-5 and acquisition energy window widths of 16-22% were optimal. For separate acquisition, we observed a broad range of optimal relative injected activities from 2.6 to 12.1 and acquisition energy window of widths 16-22%. A negative correlation between Tl-201 injected activity and the width of the Tl-201 energy window was observed in these ranges. The results also suggested that DISA methods could potentially provide image quality as good as that obtained with separate acquisition protocols. We compared observer performance for the optimized protocols and the current clinical protocol using separate acquisition. The current clinical protocols provided better performance at a cost of injecting the patient with approximately double the injected activity of Tc-99m and Tl-201, resulting in substantially increased radiation dose.

Optimization and evaluation of reconstruction-based compensation methods and reconstruction parameters for Tc-99m MIBI parathyroid SPECT.

The value of Tc-99m MIBI parathyroid SPECT for localizing parathyroid hyperplasia in chronic renal failure patients remains inconclusive due to limited image quality. Advanced reconstruction methods to improve image quality have been developed but require optimization and evaluation. The goal of this study was to optimize and evaluate compensation methods and reconstruction parameters for Tc-99m MIBI parathyroid SPECT. A phantom population and projection data that modelled clinically realistic variations found in patients were simulated. The 3D OS-EM reconstruction with compensation for attenuation, detector response and scatter in various combinations were studied. For each compensation, the number of updates for OS-EM and the cutoff frequency of a 3D Butterworth filter were optimized and evaluated using anthropomorphic model observer. With optimal parameters, the method with compensation for attenuation and detector response, with or without the addition of scatter compensation, provided the highest lesion detectability for Tc-99m MIBI parathyroid SPECT.

Design of a digital phantom population for myocardial perfusion SPECT imaging research.

Digital phantoms and Monte Carlo (MC) simulations have become important tools for optimizing and evaluating instrumentation, acquisition and processing methods for myocardial perfusion SPECT (MPS). In this work, we designed a new adult digital phantom population and generated corresponding Tc-99m and Tl-201 projections for use in MPS research. The population is based on the three-dimensional XCAT phantom with organ parameters sampled from the Emory PET Torso Model Database. Phantoms included three variations each in body size, heart size, and subcutaneous adipose tissue level, for a total of 27 phantoms of each gender. The SimSET MC code and angular response functions were used to model interactions in the body and the collimator-detector system, respectively. We divided each phantom into seven organs, each simulated separately, allowing use of post-simulation summing to efficiently model uptake variations. Also, we adapted and used a criterion based on the relative Poisson effective count level to determine the required number of simulated photons for each simulated organ. This technique provided a quantitative estimate of the true noise in the simulated projection data, including residual MC simulation noise. Projections were generated in 1 keV wide energy windows from 48-184 keV assuming perfect energy resolution to permit study of the effects of window width, energy resolution, and crosstalk in the context of dual isotope MPS. We have developed a comprehensive method for efficiently simulating realistic projections for a realistic population of phantoms in the context of MPS imaging. The new phantom population and realistic database of simulated projections will be useful in performing mathematical and human observer studies to evaluate various acquisition and processing methods such as optimizing the energy window width, investigating the effect of energy resolution on image quality and evaluating compensation methods for degrading factors such as crosstalk in the context of single and dual isotope MPS.

Optimization of energy window for 90Y bremsstrahlung SPECT imaging for detection tasks using the ideal observer with model-mismatch.

PURPOSE: In yttrium-90 ((90)Y) microsphere brachytherapy (radioembolization) of unresectable liver cancer, posttherapy (90)Y bremsstrahlung single photon emission computed tomography (SPECT) has been used to document the distribution of microspheres in the patient and to help predict potential side effects. The energy window used during projection acquisition can have a significant effect on image quality. Thus, using an optimal energy window is desirable. However, there has been great variability in the choice of energy window due to the continuous and broad energy distribution of (90)Y bremsstrahlung photons. The area under the receiver operating characteristic curve (AUC) for the ideal observer (IO) is a widely used figure of merit (FOM) for optimizing the imaging system for detection tasks. The IO implicitly assumes a perfect model of the image formation process. However, for (90)Y bremsstrahlung SPECT there can be substantial model-mismatch (i.e., difference between the actual image formation process and the model of it assumed in reconstruction), and the amount of the model-mismatch depends on the energy window. It is thus important to account for the degradation of the observer performance due to model-mismatch in the optimization of the energy window. The purpose of this paper is to optimize the energy window for (90)Y bremsstrahlung SPECT for a detection task while taking into account the effects of the model-mismatch. METHODS: An observer, termed the ideal observer with model-mismatch (IO-MM), has been proposed previously to account for the effects of the model-mismatch on IO performance. In this work, the AUC for the IO-MM was used as the FOM for the optimization. To provide a clinically realistic object model and imaging simulation, the authors used a background-known-statistically and signal-known-statistically task. The background was modeled as multiple compartments in the liver with activity parameters independently following a Gaussian distribution; the signal was modeled as a tumor with a Gaussian-distributed activity parameter located randomly with equal probability at one of three positions. The IO test statistics (i.e., likelihood ratios) were estimated using Markov-chain Monte Carlo methods. The authors realistically modeled human anatomy using a digital phantom code, and realistically simulated (90)Y bremsstrahlung SPECT imaging with a clinical SPECT system and typical imaging parameters using a previously validated Monte Carlo bremsstrahlung simulation method. Model-mismatch was included by modeling image formation process in the calculation of IO test statistics using an analytic modeling method previously developed for quantitative (90)Y bremsstrahlung SPECT. To demonstrate the effects of the model-mismatch on the detection task, the authors optimized the energy window both with and without model-mismatch included in the IO. RESULTS: For all the energy windows, the AUC values for the IO-MM were smaller than that for the IO. The optimal windows for the IO-MM and the IO were 80-180 and 60-400 keV, respectively. CONCLUSIONS: The authors have demonstrated the degradation of the ideal performance due to model-mismatch and optimized the energy window for (90)Y bremsstrahlung SPECT for detection tasks by accounting for the effects of the model-mismatch. The obtained optimal window was much narrower when taking into account the model-mismatch and similar to that obtained previously for estimation tasks.